sytx1.m

我做了很久的

r=0.02;
t=-5:r:5;
N=200;
W=2*pi*1;
k=-N:N
w=k*W/N
f1=1/2*exp(-2*(t-0.3)).*Heaviside(t-0.3);%定义法f(t-0.3)
F=r*f1*exp(-j*t'*w);%求f(t-1)的傅立叶变换
F1=abs(F)
P1=angle(F)
subplot(311)
plot(t,f1)
grid;
xlabel('t')
ylabel('f(t)')
title('f1(t-0.3)')
subplot(312);
plot(w,F1)
xlabel('w')
grid;%添加网格
ylabel('F(jw)')
subplot(313);
plot(w,P1*180/pi);
grid;
xlabel('w')
ylabel('相位P(度)');