资源描述：
《信号和线性系统分析基于MATLAB方法和实现.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、已知周期半波余弦信号和周期全波余弦信号的波形如图所示，用MATLAB编程求出它的傅立叶系数，绘出其直流、一次、二次、三次、四次及五次谐波叠加后的波形图，并将其与原周期信号的时域波形进行比较，观察周期信号的分解与合成过程。%dm09101%观察周期方波信号的分解与合成%m：傅里叶级数展开的项数display('Pleaseinputthevalueofm(傅里叶级数展开的项数)');%在命令窗口显示提示信息m=input('m=');%键盘输入傅里叶级数展开的项数t=-2*pi:0.01:2*pi;%时域波形的时间范围-2π~2π，采样间隔0.01n=round(leng
2、th(t)/4);%根据周期方波信号的周期，计算1/2周期的数据点数f=cos(t).*(t<(-3*pi/2))+cos(t).*((t>(-pi/2))&(t<(pi/2)))+cos(t).*(t>3*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([-22-0.51.5]);%指定图形显示的横坐标范围和纵坐标范围title('周期信号');%给显示的图形加上标题xla
3、bel('单位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);holdon;plot(t/pi,y(1,:),'r','LineWidth',3);%绘制谐波叠加信号holdoff;grid;axis([-22-0.51.5]);title(strcat(kk));xlabel('单位pi','Fontsize',8);kk=strcat(kk,'和第',num2str(
4、k));pause;k=k+1;x=x+0.5.*cos(t);y(2,:)=x;%计算各次谐波叠加和plot(t/pi,y(m+1,:),'LineWidth',2);holdon;plot(t/pi,y(2,:),'r','LineWidth',3);%绘制谐波叠加信号holdoff;grid;axis([-22-0.51.5]);title(strcat(kk,'次谐波叠加'));xlabel('单位pi','Fontsize',8);kk=strcat(kk,'、',num2str(k));pause;x=x+2/(3.*pi).*cos(2.*t);y(3,:
5、)=x;%计算各次谐波叠加和plot(t/pi,y(m+1,:),'LineWidth',2);holdon;plot(t/pi,y(3,:),'r','LineWidth',3);%绘制谐波叠加信号holdoff;grid;axis([-22-0.51.5]);title(strcat(kk,'次谐波叠加'));xlabel('单位pi','Fontsize',8);kk=strcat(kk,'、',num2str(k+2));pause;x=x-2/(15.*pi).*cos(4.*t);y(4,:)=x;%计算各次谐波叠加和plot(t/pi,y(m+1,:),'
6、LineWidth',2);holdon;plot(t/pi,y(4,:),'r','LineWidth',3);%绘制谐波叠加信号holdoff;grid;axis([-22-0.51.5]);title(strcat(kk,'次谐波叠加'));xlabel('单位pi','Fontsize',8);pause;plot(t/pi,y(1:m+1,:),'LineWidth',3);grid;axis([-22-0.51.5]);title('各次谐波叠加波形');xlabel('单位pi','Fontsize',8);%End%dm09201%绘制周期信号的频谱特性
7、functionCTFS_RP%以周期矩形脉冲信号为例，计算其频谱特性%Nf：级数分解的谐波次数，由键盘输入%Nn：输出数据的准确位数%a0：直流项系数%an：第1,2,3,...次谐波余弦项展开系数%bn：第1,2,3,...次谐波正弦项展开系数%tao：周期矩形脉冲信号脉宽，由键盘输入%T：周期矩形脉冲信号周期，由键盘输入display('PleaseinputthevalueofT,taoandNf');%命令窗口提示用户输入参数T=input('T=');tao=input('tao=');Nf=input('Nf=');syms
