www.pudn.com > antennas-and-radio-wave-propagation.rar > 矩形同相平面口径的辐射.m, change:2002-01-01,size:1232b


%矩形同相平面口径的辐射,图6-2-3 
a=2; 
b=3; 
%矩形口径均匀分布 
%平面极坐标方向图 
delta=eps:2*pi/90:2*pi; 
fai1=0.5*2*pi*b*sin(delta); 
fai2=0.5*2*pi*a*sin(delta); 
fe1=abs(0.5*(1+cos(delta)).*sin(fai1)./fai1);femax1=max(max(fe1)); 
fh1=abs(0.5*(1+cos(delta)).*sin(fai2)./fai2);fhmax1=max(max(fh1)); 
subplot(1,2,1);polar(delta,fe1/femax1);title('fe1');  
subplot(1,2,2);polar(delta,fh1/fhmax1);title('fh1'); 
%矩形口径余弦分布 
%平面极坐标方向图 
figure(2); 
delta=eps:2*pi/90:2*pi; 
fai1=0.5*2*pi*b*sin(delta); 
fai2=0.5*2*pi*a*sin(delta); 
fe2=abs(0.5*(1+cos(delta)).*sin(fai1)./fai1);femax2=max(max(fe2)); 
fh2=abs(0.5*(1+cos(delta)).*cos(fai2)./(1-(2/pi*fai2).^2));fhmax2=max(max(fh2)); 
subplot(1,2,1);polar(delta,fe2/femax2);title('fe2');  
subplot(1,2,2);polar(delta,fh2/fhmax2);title('fh2'); 
%主平面直角坐标方向图 
figure(3); 
angle=-pi:pi/180:pi; 
fai1=0.5*2*pi*b*sin(angle); 
fai2=0.5*2*pi*a*sin(angle); 
fe1=abs(0.5*(1+cos(angle)).*sin(fai1)./fai1); 
fh1=abs(0.5*(1+cos(angle)).*sin(fai2)./fai2); 
fe2=abs(0.5*(1+cos(angle)).*sin(fai1)./fai1); 
fh2=abs(0.5*(1+cos(angle)).*cos(fai2)./(1-(2/pi*fai2).^2)); 
plot(angle,fe1,'-');hold on; 
plot(angle,fh1,'--');hold on; 
plot(angle,fe2,'.');hold on; 
plot(angle,fh2,'-.');