信号与系统实验3 连续信号的频域分析

问题1 求实验图所示的周期信号(T=2 =1)的傅里叶级数，做出其前3、9、21、45项谐波的合成波形并与原信号做比较，作出其单边幅度谱和相位谱

t = -6:0.01:6;
d = -6:2:6;
fxx = pulstran(t,d,'tripuls'); %三角脉冲
f1 = fourierseries(3,t);    %fourierseries 自定义函数
f2 = fourierseries(9,t);
f3 = fourierseries(21,t);
f4 = fourierseries(45,t);
N = length(fxx);
f5 = fft(fxx,N);
f = (0:N-1)'*100/N;
f6 = angle(f5);
subplot(3,2,1);plot(t,fxx,'r',t,f1,'b');xlabel('t');ylabel('f(t)');title('前3项谐波');
grid on
axis([-6 6 -0.1 1.1]);
subplot(3,2,2);plot(t,fxx,'r',t,f2,'b');xlabel('t');ylabel('f(t)');title('前9项谐波');
grid on
axis([-6 6 -0.1 1.1]);
subplot(3,2,3);plot(t,fxx,'r',t,f3,'b');xlabel('t');ylabel('f(t)');title('前21项谐波');
grid on
axis([-6 6 -0.1 1.1]);
subplot(3,2,4);plot(t,fxx,'r',t,f4,'b');xlabel('t');ylabel('f(t)');title('前45项谐波');
grid on
axis([-6 6 -0.1 1.1]);
subplot(3,2,5);stem(f(1:600),abs(f5(1:600)),'.');xlabel('w');ylabel('F(jw)');title('单边频谱图');
grid on
subplot(3,2,6);stem(f6,'.');xlabel('w');ylabel('f6');title('相位谱');
grid on
axis([0 1200 -200 200]);
print(gcf,'-dpng','3-1.png');