改一下圆口径泰勒的程序

新手刚编程，希望大侠们指点一番
function I=Ip(rho,a)
%Ip为圆口径泰勒分布电流幅度
%a为圆口径半径
R0dB=30;
a=0.25;
R0=10^(R0dB/20);
A=log(R0+sqrt(R0^2-1))/pi;
n1=6;
r1=[1.2197 ;2.2331; 3.2383; 4.2411; 5.2428 ;6.2439;];%一阶bessel的前6个根/pi
phi=6.244/sqrt(A^2+(n1-0.5)^2);%波瓣展宽因子
for n=1:n1
un=phi*sqrt(A^2+(n-0.5)^2);
end
rho=1:0.01:a;
g=1;
if m==0
Bm=1;
else m=1:n1-1
R11=1;
R22=1;
for n=1:n1-1
R11=1-(r1(m,:)/un)^2;
if n=m
R22=1-(r1(m,:)/r1(1,n))^2;
end
R2=R22*R;
end
R1=R11*R;
Bm(m)=-R1/R2/besselj(0,r1(1,m)*pi);
g=g+Bm(m)*besselj(0,r1(1,m)*rho*pi/a);
end

 

已经解决了

 

能发那程序给我吗

申明：网友回复良莠不齐，仅供参考。如需专业帮助，请学习业界专家讲授的天线设计视频培训教程。