19路光干涉.rar_19路激光介绍_干涉_激光_激光分布_激光干涉

clc,clear a=10;lambda=1.080;nc=1.47;ng=1.45; %设定光束腰直径a，波长lambda，纤芯和包层折射率ng和nc V=(2*pi*a/lambda)*sqrt(nc^2-ng^2); %设定光纤V值 m = input('输入0到4的m值：'); %输入m和n值 n = input('输入1到3的n值：'); U=[2.4,5.5,8.6;3.8,7,10;5.1,8.4,11;6.3,9,13;7.5,11,14]; %对于相应的mn组合值给出对应的U值 W=V^2-U(m+1,n)^2; %利用UVW三者关系计算W值 Npoint=400; %取样点数 r1=1;r2=5; %给出纤芯直径和包层直径 R1=linspace(0,r1,Npoint); %设定极坐标点 R2=linspace(r1,r2,Npoint); Theta1=linspace(0,2*pi,Npoint); Theta2=linspace(0,2*pi,Npoint); E1=zeros(Npoint,Npoint); %生成两个零矩阵 E2=zeros(Npoint,Npoint); I1=E1;I2=E2; %% 特征方程 for i=1:Npoint %计算纤芯内的光场振幅和光强分布 for j=1:Npoint E1(i,j)=besselj(m,U(m+1,n)*R1(i))./besselj(m,U(m+1,n))*cos(m*Theta1(j)); I1(i,j)=E1(i,j).^2; end end for i=1:Npoint %计算纤芯到包层的光场振幅和光强分布 for j=1:Npoint E2(i,j)=besselk(m,W.*R2(i))./besselk(m,W)*cos(m*Theta1(j)); I2(i,j)=E2(i,j).^2; end end phase1=angle(E1)-pi/2*ones(Npoint,Npoint); %计算相位 phase2=zeros(Npoint,Npoint); E2=E2/max(max(E1)); %振幅和光强归一化 E1=E1/max(max(E1)); I2=I2/max(max(I1)); I1=I1/max(max(I1)); [Theta1, R1]=meshgrid(Theta1,R1); %生成二维网格矩阵 [Theta2, R2]=meshgrid(Theta2,R2); [X1, Y1]=pol2cart(Theta1,R1); %极坐标换直角坐标 [X2, Y2]=pol2cart(Theta2,R2); %% 作图 figure %绘制振幅分布图 pcolor(X1,Y1,E1) shading interp hold on pcolor(X2,Y2,E2) shading interp xlabel('x'),ylabel('y'),zlabel('E') axis([-2 2 -2 2]); colorbar; title(['Amplitude distribution of LP',num2str(m),num2str(n),'Mode']) figure %绘制相位分布图 pcolor(X1,Y1,phase1) shading interp hold on pcolor(X2,Y2,phase2) shading interp xlabel('x'),ylabel('y'),zlabel('phase') axis([-2 2 -2 2]); colorbar; title(['Phase distribution of LP',num2str(m),num2str(n),'Mode']) figure %绘制光强分布图 pcolor(X1,Y1,I1) shading interp hold on pcolor(X2,Y2,I2) shading interp xlabel('x'),ylabel('y'),zlabel('I') axis([-2 2 -2 2]); colorbar; title(['Light intensity distribution of LP',num2str(m),num2str(n),' Mode']) % figure, surf(X1,Y1,E1) % shading interp % hold on % surf(X2,Y2,E2) % shading interp % xlabel('x'),ylabel('y'),zlabel('E') % axis([-2 2 -2 2]); colorbar; % title(['Amplitude distribution of LP',num2str(m),num2str(n),'Mode, V = ' num2str(V)]) % figure,surf(X1,Y1,I1) % shading interp % hold on % surf(X2,Y2,I2) % shading interp % xlabel('x'),ylabel('y'),zlabel('I') % axis([-2 2 -2 2]); colorbar; % title(['Light intensity distribution of LP',num2str(m),num2str(n),' Mode, V = ' num2str(V)])