problem_6_5_iv.m

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

clear all
t0=10;
K=10;
fl=5000000000;fh=7000000000;
c=300000000;
d=(c/2/fh)*8;
for f=fl:(fh-fl)/20:fh;
    it=0;
    for t=-90:.1:90;
        tet=t*pi/180;
        it=it+1;
        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(it)=sum(F1);
    end;
    t=-90:.1:90;plot(t,20*log10(abs(F)),'b');hold on;
end;
    axis([-90 90 -60 30]);
    xlabel('angle [deg]');ylabel('dB');grid;