《信号与系统实验6.docVIP

实验六 利用MATLAB实现周期信号的分解与合成 一、实验目的 1．复习巩固信号与系统课程中关于周期信号展开成傅里叶级数等有关理论知识； 2．验证、设计并掌握利用MATLAB实现周期信号的分解与合成的方法。 二、实验内容 1．问题一 周期方波信号如图所示，求出该信号的三角函数形式的傅里叶级数，并用MATLAB编程实现其各次谐波叠加情况的观察与分析。 三角函数形式的傅里叶级数展开： 该信号的周期T=2π，基频为 因为该方波信号为奇函数，所以 而 所以 程序代码： 执行结果： 观察实验结果，理解：“当所包含的谐波分量越多时，合成波形越接近于原来的周期信号”。结合实验过程，理解上述程序中的命令和相应显示结果的对应关系。 2．问题二 求下图所示矩形脉冲信号f(t)的傅里叶级数表示式，并用MATLAB方法求出N=7和N=21时的合成图。 傅里叶级数展开： 信号f(t)的离散频谱函数（复傅里叶系数） （参考教材97页式4－8） 得傅里叶级数表达式 即f(t)的前N项合成表达式为 程序代码： 执行结果： 3．设计 信号f(t)如图所示，设T=2s、A=1，试求其 复指数形式的傅里叶级数，设计一个MATLAB程序， 通过运行程序观察此信号的分解与合成。 三、作业 1．进一步理解周期信号分解为傅里叶级数的理论； 2．完成问题二中f(t)的前N项合成表达式的详细推导。 附加说明 ◇ 关于dsolve函数 调用格式 yp=dsolve(’方程1，方程2，…’，’条件1，条件2，…’，’变量’) 简要说明 （1）函数返回微分方程1（或方程1、方程2、…构成的微分方程组）的解； （2）各个微分方程中的一阶导数符号用‘D’表示，D2、D3、…Dn分别表示二阶、三阶、…n阶导数； （3）各个微分方程中的自变量默认为‘t’，否则，可用格式中的 ’变量’ 重新指明； （4）各个’条件’项可用于指明初始条件。若不指明或指明不全，则函数的返回值中可能出现C1、C2、…等待定常数。 ◇ 关于tf函数 tf函数根据描述系统的微分方程的系数生成系统函数对象（TF对象）。 调用格式 sys=tf(b,a) 简要说明 输入参量a － 微分方程左边多项式系数构成的行向量； 输入参量b － 微分方程右边多项式系数构成的行向量； 输出参量sys － 系统函数对象 对一般形式的系统微分方程: 或 a和b应分别定义成：、。 例如，对由微分方程描述的系统，由tf函数生成其系统对象sys的命令为 上述命令运行结果为 调用tf函数生成系统函数对象sys,并用向量f和t定义了系统激励信号后，即可调用lsim函数对连续系统的响应进行仿真。 ◇ 关于lsim函数 调用格式 lsim函数的一种调用格式为： lsim(sys,f,t) 简要说明 这种调用格式对向量t定义的时间范围内的系统响应进行仿真，即绘制LTI连续系统的响应的时域波形，同时还绘出系统的激励信号对应的时域波形。 在这种调用格式中，输入参量f和t是两个表示输入信号的行向量，其中t为表示输入信号时间范围的向量，f则是输入信号在向量t定义的时间点上的采样值。 例如，命令 t=0:0.01:5; f=cos(2*t); 就定义了0－5秒时间范围内的余弦输入信号cos(2*t)（采样时间间隔为0.01秒）。 输入参量sys是由tf函数根据描述系统的微分方程的系数生成的系统函数对象。 ◇ 关于impulse函数 Impulse函数主要用于求连续系统冲激响应并绘制其时域波形。 调用格式 （1）impulse(b,a) （2）impulse(b,a,t) （3）impulse(b,a,t1:ts:t2) （4）y=impulse(b,a,t1:ts:t2) 简要说明 （1）该调用格式以默认的方式绘制由向量a和b定义的LTI连续系统的冲激响应时域波形，并提供交互式功能查看任意时刻冲激响应的信号样值。其中输入参量a和b分别是描述系统的微分方程左边和右边系数构成的行向量。 例如，对由微分方程描述的系统，则定义该系统的向量a和b应使用命令： a=[1 0 2 4]; b=[-1 6]; 。 （2）该调用格式绘制由向量a和b定义的连续系统在时间参数t所指定的0～t时间范围内系统冲激响应的时域波形。 （3