《信号与系统实验3.docVIP

9.1已知周期半波余弦信号和周期全波余弦信号的波形分别如图9-15（a）和（b）所示，用MATLAB编程求出他们的傅里叶系数，绘出其，直流、一次、二次、三次、四次及五次谐波叠加后的波形图，并将其与、原周期信号的时域波形进行比较，观察周期信号的分解与合成过程。 解：由图可以看出（a）的表达式f(t)=0.5*(cos(t)+|cos(t)|)，周期为2*pi；（b）的表达式f(t)=|cos(t)|，周期为pi。 按傅里叶级数的定义，f（t）可以由三角函数的线性组合来表示 其中 display(Please input the value of m (傅里叶级数展开的项数)); m = input(m = ); t = -3*pi:0.01:3*pi; n = round(length(t)/4); f =cos(t).*(heaviside(t+2.5*pi)-heaviside(t+1.5*pi)+heaviside(t+0.5*pi)-heaviside(t-0.5*pi)+heaviside(t-1.5*pi)-heaviside(t-2.5*pi)); y = zeros(m+1,max(size(t))); y(m+1,:) = f; figure(1); plot(t/pi,y(m+1,:)); grid; axis([-3 3 -1 1.5]); title(半波余弦信号); xlabel(单位pi,Fontsize, 8); x = zeros(size(t)); kk = 1; y(1,:)=1/pi; y(2,:)=1/pi+1/2*cos(t); x=1/pi+1/2*cos(t); for k=2:6 pause; x = x+(-2*cos(pi*k/2)/pi/(k^2-1))*cos(k*t); y((k+1),:) = x; plot(t/pi,y(m+1,:)); hold on; plot(t/pi,y((k+1),:)); hold off; grid; axis([-3 3 -1 1.5]); title(strcat(第,kk,次谐波叠加)); xlabel(单位pi,Fontsize, 8); kk = strcat(kk,、,num2str(k)); end pause; plot(t/pi,y(1:m+1,:)); grid; axis([-3 3 -1 1.5]); title(谐波叠加); xlabel(单位pi,Fontsize, 8); (2) display(Please input the value of Nf); Nf = input(Nf = ); an=zeros(Nf+1,1); cn=zeros(Nf+1,1); phase=zeros(Nf+1,1); an(1)=2/pi; for i=1:Nf an(i+1)=((-2)/(pi*(i*i-1)))*cos(i*pi/2); cn(i+1)=abs(an(i+1)); end for i=1:Nf if an(i)=0 phase(i)=0; else phase(i)=pi; end end T=2*pi; t=-5*T:0.01:5*T; x=cos(t); d=