pinyi.m

我做了很久的

R=0.02;
t=-2:R:2
f=Heaviside(t+1)-Heaviside(t-1);
f1=f.*exp(-j*20*t)
f2=f.*exp(j*20*t)
W1=2*pi*5
N=500;k=-N:N;W=k*W1/N;
F=f*exp(-j*t'*W)*R;%求f(t)的傅立叶变换
F1=f1*exp(-j*t'*W)*R%求f1(t)的傅立叶变换
F2=f2*exp(-j*t'*W)*R%求f2(t)的傅立叶变换
F=real(F)
F1=real(F1)
F2=real(F2)
subplot(2,2,1)%图形显示坐标
plot(W,F)
xlabel('w')
ylabel('F(jw)');
title('F(w)的频谱F(jw)');
subplot(2,2,2)
plot(W,F1)
xlabel('w')
ylabel('F1(jw)');
title('F(w)左移w=20处的频谱F1(jw)');
subplot(2,2,3)
plot(W,F2)
xlabel('w')
ylabel('F2(jw)');
title('F(w)右移w=20处的频谱F2(jw)');