lp模matlab,LP02模和LP22模的V-w曲线怎么修改

如何根据书上修改现在想仿真出这个曲线分析

LP01

clear

close all

tic

Vmax = 10;

N =100;

for j = 1:N

V(j) = j/N*Vmax;

Vtemp = V(j);

Utemp = NaN;

i = 0;

while (isnan(Utemp)  && i

init = 2.5*(i+1)/N-eps;

try

Utemp = fzero(@(Utemp) ...

besselj(0,Utemp)/(Utemp*besselj(1,Utemp)) - ...

besselk(0,sqrt(Vtemp^2-Utemp^2))/ ...

(sqrt(Vtemp^2-Utemp^2)*besselk(1,sqrt(Vtemp^2-Utemp^2))),init);

catch

end

i = i+1;

end

U(j) = Utemp;

end

W = sqrt(V.^2-U.^2);

Ymax = ceil(max([U,W]));

figure

subplot(1,2,1)

plot(V,U,'r');

axis equal

axis([0 Vmax 0 Ymax])

xlabel('V')

ylabel('U')

title('LP_{01}  V-U')

subplot(1,2,2)

plot(V,W,'r');

axis equal

axis([0 Vmax 0 Ymax])

xlabel('V')

ylabel('W')

title('LP_{01}  V-W')

toc

LP11模的

clc

clear

close all

tic

Vmax = 10;

N =100;

for j = 1:N

V(j) = j/N*Vmax;

Vtemp = V(j);

Utemp = NaN;

i = 0;

while (isnan(Utemp)  && i

init = 3.5*(N-i)/N;

try

Utemp = fzero(@(Utemp) ...

besselj(1,Utemp)/(Utemp*besselj(0,Utemp)) + ...

besselk(1,sqrt(Vtemp^2-Utemp^2))/(sqrt(Vtemp^2-Utemp^2)* ...

besselk(0,sqrt(Vtemp^2-Utemp^2))),init);

catch

end

i = i+1;

end

U(j) = Utemp;

end

W = sqrt(V.^2-U.^2);

Ymax = ceil(max([U,W]));

figure

subplot(1,2,1)

plot(V,U);

axis equal

axis([0 Vmax 0 Ymax])

xlabel('V')

ylabel('U')

title('LP_{1m} V-U')

subplot(1,2,2)

plot(V,W);

axis equal

axis([0 Vmax 0 Ymax])

xlabel('V')

ylabel('W')

title('LP_{1m} V-W')

toc