fig_2_11_ab.m

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

clear all
tau=5e-10;f_1=4e9;
E=1;
K=sqrt(E/tau/sqrt(pi/2));
it=0;
for t=-2e-9:1e-12:2e-9;
    it=it+1;
    s_j(it)=sqrt(2)*K*exp(-t^2/tau^2)*cos(2*pi*f_1*t);
end,
figure (1);t=-2e-9:1e-12:2e-9;plot(t,s_j,'b');
xlabel('Time, sec');ylabel('Amplitude');grid;

iif=0;
for f=0:1e8:2e10;
    iif=iif+1;
    S_j(iif)=sqrt(2)*.5*K*tau*sqrt(pi)*exp(-(pi*tau*(f-f_1))^2);
end,
figure (3);f=0:1e8:2e10;plot(f,S_j,'r')
xlabel('Frequency, Hz');ylabel('Magnitude');grid;