方向图clear lamda=10;%自由空间的波长 n=1.125; l=n*lamda; k=2*pi/lamda;%自由空间的相移常数 theta0=0.0001:0.1:360; theta=theta0*pi/180; for i=1:length(theta0) fe(i)=abs(cos(k*l*cos(theta(i)-cos(k*l)/sin(theta(i); % fh(i)=1-cos(k*l); end % figure polar(theta,fe/max(fe); xlabel(l=1.125); % figure % polar(theta,fh/max(fh); 辐射阻抗clear lamda=10;%自由空间的波长 n=0.1:0.01:0.9; l=n*lamda; k=2*pi/lamda;%自由空间的相移常数 gama=0.5772156; for i=1:length(n) Rr1(i)=sin(2*k*l(i)*(Si(4*k*l(i)
