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

昨天试着编了一下泰勒线源综合的程序，结果如下，输入副瓣电平SLL，阵元个数n，等副瓣数n1，就可以得到电流分布I（i） O3@DU#N&s  
SLL=input('请输入SLL=');n=input('请输入n=');n1=input('请输入n1='); vSR&>Q%X  
R0=10^(SLL/-20);A=(1/pi)*acosh(R0);n0=2*A^2+4; $KbZ4bB[Bo  
x=ones(1,n);lamda_g=10.73; d=lamda_g/2;L=(n)*d; a~7D4G  
for i=1:n R>O_2`c  
    if i<=n1-1 |-x-CSN  
        x(i)=(n1)^2*(A^2+(i-1/2)^2)/(A^2+(n1-1/2)^2); n7fhc*}:`  
    else s{@R|5  
         x(i)=i^2; Z?{\34lPj  
    end _#sy  
     g8N"-j&@  
end Q4m> 3I  
for j=1:n &N^~=y^`C'  
    if j<n1 BQ0?B*yqd  
       for k=1:n1-1 >8_y-74  
           B(j,k)=1-j^2/x(k); rS_G;}Zr  
       end; zuLW'a6F-
 
       B1=B(j,:); bI &<L O  
       s(j)=(factorial(n1-1))^2/(factorial(n1-1+j)*factorial(n1-1-j))*prod(B1); (U.Go/A#wE  
    else Mwnr4$]  
        s(j)=0; ^$:w  
    end iSCkV2  
end p qz~9y~
 
if mod(n,2)>0;   =$J(]KPv!?  
    for o=1:(n-1)/2+1 !^%3  
        z(o)=(o-1)*d; LZn'+{\`  
    end !G-+O#W`  
       for l=1:(n-1)/2+1 kM506U<g  
           for m=1:n1-1 \:4*h  
                z0(l,m)=s(m)*cos(2*m*pi*z(l)/L); #gQn3.PX+y  
           end oRCc8&  
            I1(l)=1+2*sum(z0(l,:)); btuG%D{a^  
       end h.X4x2(.  
    max1=I1(1); oh
8:1E,I  
    for i=1:(n-1)/2 L*Z.T^h  
        I(i)=I1((n-1)/2+2-i)/max1;I(n+1-i)=I(i); $XhMI;h  
    end eLSzGbKf  
    I((n-1)/2+1)=1; Ma|4nLC}  
else   @!P2f   
     for o=1:n/2 f\hMTebma$  
        z(o)=(2*o-1)/2*d; 5`.CzQVb  
    end 7<xnE]jdq  
       for l=1:n/2 A0Mjk  
           for m=1:n1-1 2'?'dfj  
                z0(l,m)=s(m)*cos(2*m*pi*z(l)/L); &YXJ{<s  
           end XLn9NBT4K  
        &n .. ^xa, r#N:V  
,,{;G'R|  

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

结果怎么出不来。。