problem_3_2_c.m

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

clear all;
iif=0;
for f=0:.001:2;
    iif=iif+1;
    if f==0;X3(iif)=2;else X3(iif)=(-1/(4*pi^2*f^2))*(1-j*4*pi*f*exp(-j*4*pi*f)-exp(-j*4*pi*f));end;
end,

figure (1);f=0:.001:2;plot(f,abs(X3),'r');grid;
xlabel('f');ylabel('|X_3(f)|');grid on;
figure (2);f=0:.001:2;plot(f,angle(X3),'r');grid;
xlabel('f');ylabel('Phase of X_3(f)');grid on;
%