zmn.m

来自「用矩量法和正交模法分析四元八木天线的特性。两种方法的对比。用E面、H面和方向性图」· M 代码 · 共 40 行

function z=zmn(b,c,d,e,f,g)
n=9;     
a=5*10^(-3);        
k=2*pi;    
w=3*10^8*2*pi;   
E=8.854*10^(-12);  
u=4*pi*10^(-7);      
for p=1:n
    for q=1:n
        if p==q
            s(p,q)=log(d/a)/(2*pi*d)-i*k/(4*pi);
        else
            s(p,q)=exp(-i*k*e((p+b*n),(q+c*n)))/(4*pi*e((p+b*n),(q+c*n)));
        end
    end
end
for p=1:n
    for q=1:n
        if p==q+1
            s1(p,q)=log(d/a)/(2*pi*d)-i*k/(4*pi);
        else
             s1(p,q)=exp(-i*k*f((p+b*n),(q+c*n)))/(4*pi*f((p+b*n),(q+c*n)));
         end
     end
 end
 for p=1:n
     for q=1:n
         if p==q-1
              s2(p,q)=log(d/a)/(2*pi*d)-i*k/(4*pi);
          else
               s2(p,q)=exp(-i*k*g((p+b*n),(q+c*n)))/(4*pi*g((p+b*n),(q+c*n)));
           end
       end
   end
   for p=1:n
       for q=1:n
           z(p,q)=i*w*u*d*d*s(p,q)+(2*s(p,q)-s1(p,q)-s2(p,q))/(i*w*E);
       end
   end