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');```