[RFEDA原创]分享一个泰勒线源综合的matlab程序

昨天试着编了一下泰勒线源综合的程序，结果如下，输入副瓣电平SLL，阵元个数n，等副瓣数n1，就可以得到电流分布I（i） S?z j&XY3  
SLL=input('请输入SLL=');n=input('请输入n=');n1=input('请输入n1='); _G8y9!J  
R0=10^(SLL/-20);A=(1/pi)*acosh(R0);n0=2*A^2+4;
_<~Vxz9  
x=ones(1,n);lamda_g=10.73; d=lamda_g/2;L=(n)*d; T+{'W  
for i=1:n #?d>S;)+  
    if i<=n1-1 &b]KMAo3  
        x(i)=(n1)^2*(A^2+(i-1/2)^2)/(A^2+(n1-1/2)^2); 3z. >b  
    else h|>n3-k|p  
         x(i)=i^2; e$+? v2.  
    end n\)f.}YD8d  
     =Q0)t_z_  
end m?CjYqvf  
for j=1:n McI4oD~"  
    if j<n1 ] j?Fk$C  
       for k=1:n1-1 `DY4d$!4  
           B(j,k)=1-j^2/x(k); NU\t3JaR  
       end; e5\1k#@  
       B1=B(j,:); =&G|} M  
       s(j)=(factorial(n1-1))^2/(factorial(n1-1+j)*factorial(n1-1-j))*prod(B1); RKZk/ly  
    else x=r6vOj  
        s(j)=0; }Mlz\'{  
    end e&7}N Za
 
end v__Go kj-  
if mod(n,2)>0;   &aht K}u  
    for o=1:(n-1)/2+1 #^[N4uV  
        z(o)=(o-1)*d; +fR`@HI  
    end H.S|njn:r  
       for l=1:(n-1)/2+1 fOfz^W  
           for m=1:n1-1 zW&O>H  
                z0(l,m)=s(m)*cos(2*m*pi*z(l)/L);
b,V=B{(~  
           end 'GO..m"G  
            I1(l)=1+2*sum(z0(l,:)); O6 J<Lqgh  
       end M@86u^80  
    max1=I1(1); |,~A9  
    for i=1:(n-1)/2 !8M'ms>s=  
        I(i)=I1((n-1)/2+2-i)/max1;I(n+1-i)=I(i); WD1$"}R
 
    end O g~"+IGp  
    I((n-1)/2+1)=1; /9#jv]C:  
else   <iqyDPj  
     for o=1:n/2 TZ]Gl4@  
        z(o)=(2*o-1)/2*d; <NsT[r~C  
    end qk"=nAJX  
       for l=1:n/2 gV>\lMc[-%  
           for m=1:n1-1 i-W2!;G  
                z0(l,m)=s(m)*cos(2*m*pi*z(l)/L); IY6S\Gn  
           end /[T8/7;_l  
        &n .. 9r*T3=u.S  
@Oay$gP{T  

未注册仅能浏览部分内容，查看全部内容及附件请先 登录 或 注册

结果怎么出不来。。