信号与线性系统分析（基于MATLAB的方法与实现)第九章.docVIP

专业资料 word完美格式 已知周期半波余弦信号和周期全波余弦信号的波形如图所示，用MATLAB编程求出它的傅立叶系数，绘出其直流、一次、二次、三次、四次及五次谐波叠加后的波形图，并将其与原周期信号的时域波形进行比较，观察周期信号的分解与合成过程。 % dm09101 % 观察周期方波信号的分解与合成 % m：傅里叶级数展开的项数 display(Please input the value of m (傅里叶级数展开的项数)); % 在命令窗口显示提示信息 m = input(m = ); % 键盘输入傅里叶级数展开的项数 t = -2*pi:0.01:2*pi; % 时域波形的时间范围-2π~2π，采样间隔0.01 n = round(length(t)/4); % 根据周期方波信号的周期，计算1/2周期的数据点数 f = cos(t).*(t(-3*pi/2))+cos(t).*((t(-pi/2))(t(pi/2)))+cos(t).*(t3*pi/2); %构造周期方波信号 y = zeros(m+1,max(size(t))); y(m+1,:) = f; figure(1); plot(t/pi,y(m+1,:),LineWidth,2); %绘制方波信号 grid; %在图形中加入栅格 axis([-2 2 -0.5 1.5]); %指定图形显示的横坐标范围和纵坐标范围 title(周期信号); %给显示的图形加上标题 xlabel(单位pi,Fontsize, 8); %显示横坐标单位 x = zeros(size(t)); kk = 直流分量; pause; k=1; x = x+1/pi; y(1,:) = x; %计算各次谐波叠加和 plot(t/pi,y(m+1,:),LineWidth,2); hold on; plot(t/pi,y(1,:),r,LineWidth,3); %绘制谐波叠加信号 hold off; grid; axis([-2 2 -0.5 1.5]); title(strcat(kk)); xlabel(单位pi,Fontsize, 8); kk = strcat(kk,和第,num2str(k)); pause; k=k+1; x = x+0.5.*cos(t); y(2,:) = x; %计算各次谐波叠加和 plot(t/pi,y(m+1,:),LineWidth,2); hold on; plot(t/pi,y(2,:),r,LineWidth,3); %绘制谐波叠加信号 hold off; grid; axis([-2 2 -0.5 1.5]); title(strcat(kk,次谐波叠加)); xlabel(单位pi,Fontsize, 8); kk = strcat(kk,、,num2str(k)); pause; x = x+2/(3.*pi).*cos(2.*t); y(3,:) = x; %计算各次谐波叠加和 plot(t/pi,y(m+1,:),LineWidth,2); hold on; plot(t/pi,y(3,: