云南大学信息学院 信号与系统实验三.docVIP

9.1(1) display(Please input the value of m (傅里叶级数展开的项数)); % 在命令窗口显示提示信息 m = input(m = ); % 键盘输入傅里叶级数展开的项数 t = -4*pi:0.01:4*pi; % 时域波形的时间范围-2π~2π，采样间隔0.01 f=(cos(t)0).*cos(t); y=zeros(m+3,max(size(t))); y(m+3,:) = f; figure(1); plot(t/pi,y(m+3,:)); %绘制方波信号 grid; %在图形中加入栅格 axis([-5 5 -0.5 1.1]); %指定图形显示的横坐标范围和纵坐标范围 title(半波整流); %给显示的图形加上标题 xlabel(单位pi,Fontsize, 8); pause; x=1/pi*ones(1,length(t)); y(m+2,:)=x; plot(t/pi,y(m+2,:)); grid; title(直流分量); xlabel(单位pi,Fontsize, 8); pause; x=x+1/2*cos(t); y(m+1,:)=x; plot(t/pi,y(m+1,:)); title(strcat(第1次谐波叠加)); grid; kk = 1; for k=2:2:2*m %循环显示谐波叠加图形 pause; x = x-2/pi*(1/(k^2-1))*cos(k*pi/2+k*t); y(k/2,:) = x; plot(t/pi,y(m+3,:)); hold on; %计算各次谐波叠加和 plot(t/pi,y(k/2,:)); %绘制谐波叠加信号 hold off; grid; axis([-5 5 -0.5 1.1]); title(strcat(第,kk,次谐波叠加)); xlabel(单位pi,Fontsize, 8); kk = strcat(kk,、,num2str(k+2)); end pause; plot(t/pi,y(1:m+3,:)); grid; axis([-5 5 -0.5 1.1]); title(各次谐波叠加波形); xlabel(单位pi,Fontsize, 8); % End Please input the value of m (傅里叶级数展开的项数) m = %9.1(2) display(Please input the value of m (傅里叶级数展开的项数)); % 在命令窗口显示提示信息 m = input(m = ); % 键盘输入傅里叶级数展开的项数 t = -4*pi:0.01:4*pi; % 时域波形的时间范围-2π~2π，采样间隔0.01 f=abs(cos(t)); y=zeros(m+1,max(size(t))); y(m+1,:) = f; figure(1); plot(t/pi,y(m+1,:)); %绘制方波信号 grid; %在图形中加入栅格 axis([-5 5 -0.1 1.1]); %指定图形显示的横坐标范围和纵坐标范围 title(全波整流); %给显示的图形加上标题 xlabel(单位pi,Fontsize, 8); pause; x=2/pi*ones(1,length(t)); y(m+2,:)=x; plot(t/pi,y(m+2,:)); grid; title(直流分量); xlabel(单位pi,Fontsize, 8); kk = 1; for k=2:2:2*m %循环显示谐波叠加图形 pau