资源描述：
《利用matlab,仿真偏振控制器(利用邦加球算法)》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、利用MATLAB，仿真偏振控制器（利用邦加球算法），验证固定偏振态到任意偏振态的转换，其结果布满整个邦加球-MATLABsimulationusingpolarimetriccontroller(with),fixedballalgorithmtoanypolarizationtransformationofpolarization,andtheresultisfullofwholenationsplustheballclearclc%QHQsubplot(1,2,1)[x,y,z]=sphere;mesh(x,y,z);x
2、label('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),'*')%画出入射光
3、的偏振态位置%各波片的方位角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);%求出射光的偏振态holdonfori=1:5000a=theta1(i);b=theta2(
4、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
5、);%[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);%正
6、态分布均值为0方差为0.005Delta_alpha=normrnd(0,0.005,1,1000);Delta_beta=normrnd(0,0.005,1,1000);fori=1:1000%当波片的方位角有误差存在时，求出其所对应的偏振态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(Delt
7、a_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)*c
8、os(theta_3)+sin(theta_3)*sin(theta_3)))-xi_out);subplot(2,3,1)xlabel('Deltaalpha/rad'),ylabel('Deltatheta/rad'),plot(Delta_alpha(i),Delta_the
