www.pudn.com > antennas-and-radio-wave-propagation.rar > 行波单导线方向图.m, change:2002-01-01,size:1075b


%行波单导线方向图 
%立体图 
sita=meshgrid(eps:pi/180:pi); 
fai=meshgrid(eps:2*pi/180:2*pi)'; 
t1=1.0; 
r=abs(sin(sita).*sin(pi*t1*(1-cos(sita)))./(1-cos(sita)+eps)); 
rmax=max(max(r)); 
[x,y,z]=sph2cart(fai,pi/2-sita,r/rmax); 
subplot(2,3,1);mesh(x,y,z);title('t1=1.0'); 
t1=1.5; 
r=abs(sin(sita).*sin(pi*t1*(1-cos(sita)))./(1-cos(sita)+eps)); 
rmax=max(max(r)); 
[x,y,z]=sph2cart(fai,pi/2-sita,r/rmax); 
subplot(2,3,2);mesh(x,y,z);title('t1=1.5'); 
t1=3; 
r=abs(sin(sita).*sin(pi*t1*(1-cos(sita)))./(1-cos(sita)+eps)); 
rmax=max(max(r)); 
[x,y,z]=sph2cart(fai,pi/2-sita,r/rmax); 
subplot(2,3,3);mesh(x,y,z);title('t1=3'); 
%平面图 
delta=0:2*pi/90:2*pi; 
t1=1; 
f=abs(sin(delta).*sin(pi*t1*(1-cos(delta)))./(1-cos(delta)+eps)); 
F=f./max(max(f)); 
subplot(2,3,4);polar(delta,F);title('t1=1'); 
t1=1.5; 
f=abs(sin(delta).*sin(pi*t1*(1-cos(delta)))./(1-cos(delta)+eps)); 
F=f./max(max(f)); 
subplot(2,3,5);polar(delta,F);title('t1=1.5'); 
t1=3; 
f=abs(sin(delta).*sin(pi*t1*(1-cos(delta)))./(1-cos(delta)+eps)); 
F=f./max(max(f)); 
subplot(2,3,6);polar(delta,F);title('t1=3');