www.pudn.com > antennas-and-radio-wave-propagation.rar > zxz.m, change:2002-01-01,size:308b


function y=zxz(sita); 
global n posai sita d; 
f1=abs(sin(n*posai/2)); 
eps=2.2204e-016; 
f2=abs(sin(posai/2)); 
for j=1:length(posai); 
    if f1(j)<eps&f2(j)<eps; 
        f1(j)=abs(n/2*cos(n/2*posai(j))); 
        f2(j)=abs(1/2*cos(posai(j)/2)); 
    end; 
end; 
y=f1./f2; 
y=y/max(y); 
y=y.*y.*sin(sita);