fig_6_13_15.m

来自「Ultra.Wideband.Signals.and.Systems.in.Co」· M 代码 · 共 28 行

clear all
t0=10;
K=10;
fl=4000000000;fh=8000000000;
c=300000000;
d=c/2/fh;
for f=fl:(fh-fl)/4:fh;
    for t=-90:90;
        tet=t*pi/180;
        tt=t+91;
        for k=1:K,
            al=-(2*pi*f/c)*d*sin(t0*pi/180);
            V(k)=1/1*exp(j*k*al);
            F1(k)=V(k)*exp(j*(2*pi*f/c)*k*d*sin(tet));
        end;
        F(tt)=sum(F1);
    end;
    AA = [-90:1:90];   
    plot(AA,20*log10(abs(F)),'b');hold on;
end;
grid on;
axis([-90 90 -60 30]);
xlabel('Angle, degrees');ylabel('Magnitude, dB');
legend('f=4 GHz','f=5 GHz','f=6 GHz','f=7 GHz','f=8 GHz');