复件 antemom.asv

本程序采用矩量法计算全波振子的电流分布情况。

a=0.005;
sum=0;
lanbta=1;
m=21;
v=[0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0];
l=1;
j=sqrt(-1);
for n=1:21;
    for k=1:21;
        h1=exp(-j*2*pi*(l*(abs((k-n)+0.005))/(lanbta*m)));
        h4=exp(-j*2*pi*(l*(abs((k-n)+0.005))/(lanbta*m)));
        
        if (n>k)
            h2=exp(-j*2*pi*(l*(abs(k-n)-1)/(lanbta*m)));
        end
        if (n==k)
            h2=exp(-j*2*pi*(l/(lanbta*m)));
        end
        if (n<k)
            h2=exp(-j*2*pi*(l*(abs((k-n)+1))/(lanbta*m)));
        end
        
       if (n>k)
            h3=exp(-j*2*pi*(l*(abs((k-n)+1))/(lanbta*m)));
        end
        if (n==k)
            h3=exp(-j*2*pi*(l/(lanbta*m)));
        end
        if (n<k)
            h3=exp(-j*2*pi*(l*(abs(k-n)-1)/(lanbta*m)));
        end
        
        r1=abs(n-k)+0.005;
        r4=abs(n-k)+0.005;
        
        if (n>k) r2=(abs(k-n)-1)+0.005;      end
        if (n==k) r2=1;                      end
        if (n<k) r2=abs(n-k)+1;              end
        
        if (n>k) r3=abs(n-k)+1;              end
        if (n==k) r3=1;                      end
        if (n<k) r3=(abs(n-k)-1)+0.005;      end
        
        if (n==k)
            z(n,k)=j.*377*l/(2*lanbta*m)*(h1/r1)-j*4.775*m*(lanbta/l)*(h1/r1-h2/r2-h3/r3+h4/r4);
        else
            z(n,k)=j*377*l/(2*lanbta*m)*(log((abs(n-k)+1/2)/(abs(n-k)-1/2))-j*((2*pi*l)/(lanbta*m)))-j*4.775*m*(lanbta/l)*(h1/r1-h2/r2-h3/r3+h4/r4);
        end
    end    
end
z(11,11);

p=11;
y=inv(z);
for q=1:21;
    c(q)=abs((y(q,11)))
end   
q=1:21; 

plot((q-9)/32,1.15*10^4*c(q),'k--p');
hold on

i=1:21;
plot((i-9)/32,abs(sin((i-1)/3.2)),'b')
ylabel('归一化电流')
title('全波振子电流分布')
text(0.2,0.3,'矩量法')
text(0,0.9,'sin(x)')