利用MATLAB-仿真偏振控制器(利用邦加球算法)

______________________________________________________________________________________________________________精品資料利用MATLAB，仿真偏振控制器（利用邦加球算法），驗證固定偏振態(tài)到任意偏振態(tài)的轉(zhuǎn)換，其結(jié)果布滿整個邦加球-MATLABsimulationusingpolarimetriccontroller(with),fixedballalgorithmtoanypolarizationtransformationofpolarization,andtheresultisfullofwholenationsplustheballclearclc%QHQsubplot(1,2,1)[x,y,z]=sphere;mesh(x,y,z);xlabel('S1')ylabel('S2')zlabel('S3')boxoncolormap(bone);holdon%入射光%I=[1;0;0;1];%右旋圓偏振光(可見)%I=[1;0;0;-1];%左旋圓偏振光%I=[1;1;0;0];%水平線偏振光%I=[1;-1;0;0];%垂直線偏振光(可見)%I=[1;0;1;0];%45度線偏振光I=[1;0;-1;0];%-45度線偏振光(可見)%I=[1;-0.2;-0.5;0.8994];%某一橢圓偏振光plot3(I(2),I(3),I(4),'*')%畫出入射光的偏振態(tài)位置%各波片的方位角theta1=rand(1,5000)*pi-pi/2;theta2=rand(1,5000)*pi-pi/2;theta3=rand(1,5000)*pi-pi/2;subplot(1,2,2)[x,y,z]=sphere;mesh(x,y,z);xlabel('S1')ylabel('S2')zlabel('S3')axis([-1,1,-1,1,-1,1]);boxoncolormap(bone);%求出射光的偏振態(tài)holdonfori=1:5000a=theta1(i);b=theta2(i);c=theta3(i);M3=Muller(a,pi/2);M2=Muller(b,pi);M1=Muller(c,pi/2);M=M3*M2*M1;O=M*I;O_1=O(2);O_2=O(3);O_3=O(4);plot3(O_1,O_2,O_3,'r*')holdonend1;偏振控制器，利用綁架球算法對QHQ型進行方位角的誤差分析%已知輸入輸出，分析方位角誤差clearclc%I=[1;0.2;0.4;0.8944];%采用斯托克斯方法輸入的情況。%O=[1;0;0;1];%[I1,I2,I3]=sop_3(I);%[O1,O2,O3]=sop_3(O);%theta_in=atan(I2/I1);%theta_out=atan(O2/O1);%xi_in=asin(I3/sqrt(I1^2+I2^2+I3^2));%xi_out=asin(O3/sqrt(O1^2+O2^2+O3^2));theta_in=30*pi/180;%直接輸入的情況。xi_in=45*pi/180;theta_out=135*pi/180;xi_out=22.5*pi/180;Delta_gamma=normrnd(0,0.005,1,1000);%正態(tài)分布均值為0方差為0.005Delta_alpha=normrnd(0,0.005,1,1000);Delta_beta=normrnd(0,0.005,1,1000);fori=1:1000%當(dāng)波片的方位角有誤差存在時，求出其所對應(yīng)的偏振態(tài)theta_3=theta_out-xi_out-Delta_alpha(i)-atan(tan(xi_in)/cos(Delta_alpha(i)))+xi_in+2*Delta_beta(i);xi_3=-acos(sqrt(cos(Delta_alpha(i))*cos(Delta_alpha(i))*cos(xi_in)*cos(xi_in)+sin(xi_in)*sin(xi_in)));Delta_theta(i)=1/2*(Delta_gamma(i)+atan(tan(xi_3)/cos(theta_out+Delta_gamma(i)-theta_3)));Delta_xi(i)=1/2*(acos(sqrt(cos(theta_out+Delta_gamma(i)-theta_3)*cos(theta_out+Delta_gamma(i)-theta_3)*cos(theta_3)*cos(theta_3)+sin(theta_3)*sin(theta_3)))-xi_out);subplot(2,3,1)xlabel('\Delta\alpha/rad'),ylabel('\Delta\theta/rad'),plot(Delta_alpha(i),Delta_theta(i),'r')holdonsubplot(2,3,2)xlabel('\Delta\beta/rad'),ylabel('\Delta\theta/rad'),plot(Delta_beta(i),Delta_theta(i),'g')holdonsubplot(2,3,3)xlabel('\Delta\gamma/rad'),ylabel('\Delta\theta/rad'),plot(Delta_gamma(i),Delta_theta(i))holdonsubplot(2,3,4)xlabel('\Delta\alpha/rad'),ylabel('\Delta\epsilon/rad'),plot(Delta_alpha(i),Delta_xi(i),'r')holdonsubplot(2,3,5)xlabel('\Delta\beta/rad'),ylabel('\Delta\epsilon/rad'),plot(Delta_beta(i),Delta_xi(i),'g')holdonsubplot(2,3,6)xlabel('\Delta\gamma/rad'),ylabel('\Delta\epsilon/rad'),plot(Delta_gamma(i),Delta_xi(i))holdonend...偏振控制器，利用綁架球?qū)HQ型進行波長誤差分析-Polarization...clearclc%I=[1;0.2;0.4;0.8944];%采用斯托克斯方法輸入的情況。%O=[1;0;0;1];%[I1,I2,I3]=sop_3(I);%[O1,O2,O3]=sop_3(O);%theta_in=atan(I2/I1);%theta_out=atan(O2/O1);%xi_in=asin(I3/sqrt(I1^2+I2^2+I3^2));%xi_out=asin(O3/sqrt(O1^2+O2^2+O3^2));theta_in=30*pi/180;%直接輸入的情況。xi_in=45*pi/180;theta_out=135*pi/180;xi_out=22.5*pi/180;Delta_gamma=normrnd(0,0.005,1,1000);%正態(tài)分布均值為0方差為0.005Delta_alpha=normrnd(0,0.005,1,1000);Delta_beta=normrnd(0,0.005,1,1000);fori=1:1000theta_3=theta_out-xi_out-Delta_alpha(i)-atan(tan(xi_in)/cos(Delta_alpha(i)))+xi_in+2*Delta_beta(i);xi_3=-acos(sqrt(cos(Delta_alpha(i))*cos(Delta_alpha(i))*cos(xi_in)*cos(xi_in)+sin(xi_in)*sin(xi_in)));Delta_theta=1/2*(Delta_gamma(i)+atan(tan(xi_3)/cos(theta_out+Delta_gamma(i)-theta_3)));Delta_xi=1/2*(acos(sqrt(cos(theta_out+Delta_gamma(i)-theta_3)*cos(theta_out+Delta_gamma(i)-theta_3)*cos(theta_3)*cos(theta_3)+sin(theta_3)*sin(theta_3)))-xi_out);subplot(3,2,1)xlabel('\Delta\alpha/rad'),ylabel('\Delta\theta/rad'),plot(Delta_alpha(i),Delta_theta,'r')holdonsubplot(3,2,3)xlabel('\Delta\beta/rad'),ylabel('\Delta\theta/rad'),plot(Delta_beta(i),Delta_theta,'g')holdonsubplot(3,2,5)xlabel('\Delta\gamma/rad'),ylabel('\Delta\theta/rad'),plot(Delta_gamma(i),Delta_theta)holdonsubplot(3,2,2)xlabel('\Delta\alpha/rad'),ylabel('\Delta\epsilon/rad'),plot(Delta_alpha(i),Delta_xi,'r')holdonsubplot(3,2,4)xlabel('\Delta\beta/rad'),ylabel('\Delta\epsilon/rad'),plot(Delta_beta(i),Delta_xi,'g')holdonsubplot(3,2,6)xlabel('\Delta\gamma/rad'),ylabel('\Delta\epsilon/rad'),plot(Delta_gamma(i),Delta_xi)holdonend...figure('color','k','Position',[1311280696]);axes1=axes('position',[0.20.60.30.3]);axis(axes1,'equal');axis(axes1,'off');holdonrl=1;rs=0.5;N=30;t1=linspace(0,2*pi,N);xl=rl*cos(t1);yl=rl*sin(t1);t2=linspace(0,2*pi,100);xs=rs*cos(t2);ys=rs*sin(t2);fori=1:Nplot(xl(i)+xs,yl(i)+ys);endaxes2=axes('position',[0.50.60.30.3]);axis(axes2,'equal');axis(axes2,'off');holdonrl=1;rs=1;N=30;t1=linspace(0,2*pi,N);xl=rl*cos(t1);yl=rl*sin(t1);t2=linspace(0,2*pi,100);xs=rs*cos(t2);ys=rs*sin(t2);fori=1:Nplot(xl(i)+xs,yl(i)+ys);endaxes3=axes('position',[0.350.200.30.3]);axis(axes3,'equal');axis(axes3,'off');holdonrl=1;rs=1.5;N=30;t1=linspace(0,2*pi,N);xl=rl*cos(t1);yl=rl*sin(t1);t2=linspace(0,2*pi,100);xs=rs*cos(t2);ys=rs*sin(t2);fori=1:Nplot(xl(i)+xs,yl(i)+ys);end%存腳本文件heart.m%%一個半徑為R的圓周等分成N份，然后一等分點為圓心