资源描述：
《信号与系统_实验6new》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、实验六利用MATLAB实现周期信号的分解与合成一、实验目的1．复习巩固信号与系统课程中关于周期信号展开成傅里叶级数等有关理论知识；2．验证、设计并掌握利用MATLAB实现周期信号的分解与合成的方法。二、实验内容1．问题一周期方波信号如图所示，求出该信号的三角函数形式的傅里叶级数，并用MATLAB编程实现其各次谐波叠加情况的观察与分析。三角函数形式的傅里叶级数展开：该信号的周期T=2π，基频为因为该方波信号为奇函数，所以而所以程序代码：%观察周期方波信号的分解与合成%m：傅里叶级数展开的项数display('Ple
2、aseinputthevalueofm（傅里叶级数展开的项数）');%在命令窗口中显示提示m=input('m=');%键盘输入傅里叶级数展开的项数t=-2*pi:0.01:2*pi%时域波形的时间范围-2π～2π，采样间隔0.01n=round(length(t)/4);%根据周期方波信号的周期，计算1/2周期的数据点数f=[ones(n,1);-1*ones(n,1);ones(n,1);-1*ones(n+1,1)];%构造周期方波信号y=zeros(m+1,max(size(t)));y(m+1,:)=f
3、';figure(1);plot(t/pi,y(m+1,:));%绘制方波信号gridon;axis([-22-1.51.5]);title('周期方波');xlabel('单位：pi','Fontsize',8);x=zeros(size(t));kk='1';fork=1:2:2*m-1%循环显示谐波叠加图形pause;x=x+sin(k*t)/k;y((k+1)/2,:)=4/pi*x;%计算各次谐波叠加和plot(t/pi,y(m+1,:));holdon;plot(t/pi,y((k+1)/2,:));
4、%绘制谐波叠加信号holdoff;grid;axis([-22-1.51.5]);title(strcat('第',kk,'次谐波叠加'));xlabel('单位：pi','Fontsize',8);kk=strcat(kk,'、',num2str(k+2));endpause;plot(t/pi,y(1:m+1,:))grid;axis([-22-1.51.5]);title('各次谐波叠加波形');xlabel('单位：pi','Fontsize',8);%end执行结果：观察实验结果，理解：“当所包含的谐波
5、分量越多时，合成波形越接近于原来的周期信号”。结合实验过程，理解上述程序中的命令和相应显示结果的对应关系。2．问题二求下图所示矩形脉冲信号f(t)的傅里叶级数表示式，并用MATLAB方法求出N=7和N=21时的合成图。傅里叶级数展开：信号f(t)的离散频谱函数（复傅里叶系数）（参考教材97页式4－8）得傅里叶级数表达式即f(t)的前N项合成表达式为t=-2:0.001:2;N=input('N=');F0=0.5;fN=F0*ones(1,length(t));forn=1:2:NfN=fN+cos(pi*n*t
6、)*sinc(n/2);endplot(t,fN);title(['N='num2str(N)])axis([-22-0.21.2]);程序代码：执行结果：3．设计信号f(t)如图所示，设T=2s、A=1，试求其复指数形式的傅里叶级数，设计一个MATLAB程序，通过运行程序观察此信号的分解与合成。三、作业1．进一步理解周期信号分解为傅里叶级数的理论；2．完成问题二中f(t)的前N项合成表达式的详细推导。附加说明◇关于dsolve函数调用格式yp=dsolve(’方程1，方程2，…’，’条件1，条件2，…’，’变量
7、’)简要说明（1）函数返回微分方程1（或方程1、方程2、…构成的微分方程组）的解；（2）各个微分方程中的一阶导数符号用‘D’表示，D2、D3、…Dn分别表示二阶、三阶、…n阶导数；（3）各个微分方程中的自变量默认为‘t’，否则，可用格式中的’变量’重新指明；（4）各个’条件’项可用于指明初始条件。若不指明或指明不全，则函数的返回值中可能出现C1、C2、…等待定常数。◇关于tf函数tf函数根据描述系统的微分方程的系数生成系统函数对象（TF对象）。调用格式sys=tf(b,a)简要说明输入参量a－微分方程左边多项式系
8、数构成的行向量；输入参量b－微分方程右边多项式系数构成的行向量；输出参量sys－系统函数对象对一般形式的系统微分方程:或a和b应分别定义成：、。a=[132];b=[-12];sys=tf(b,a)%调用tf函数生成系统函数对象sys例如，对由微分方程描述的系统，由tf函数生成其系统对象sys的命令为Transferfunction:-s+2----------s^2+3
