lp模matlab,有会用MATLAB做LP模场分布的吗

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

clc,clear,close all

V=5.500;

U=3.468;

W=4.269;

Npoint=201;

R1=linspace(0,1,Npoint);

R2=linspace(1,5,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(1,U*R1(i))*cos(Theta1(j));

I1(i,j)=E1(i,j).^2;

end

end

for i=1:Npoint

for j=1:Npoint

E2(i,j)=besselj(1,U).*besselk(1,W.*R2(i))./besselk(1,W)*cos(Theta1(j));

I2(i,j)=E2(i,j).^2;

end

end

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,surf(X1,Y1,E1)

shading interp

hold on

surf(X2,Y2,E2)

shading interp

xlabel('x'),ylabel('y'),zlabel('E')

title(['LP_{11} 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')

title(['LP_{11} Mode, V = ' num2str(V)])