离散时间信号通过线性时不变系统实验报告..doc

离散时间信号通过线性时不变系统 1.1 实验目的 1、通过本实验，进一步加深对离散线性时不变系统的理解。 2、掌握利用线性卷积求解离散线性时不变系统输出的方法。 3、掌握利用差分方程求解离散线性时不变系统输出的方法。 1.2 实验原理与方法 离散时间系统的输出信号与输入信号之间的关系为： 离散时间线性时不变系统的特点是系统具有线性性质和线性时不变性质。即满足线性叠加定理： 系统对输入信号的运算关系在整个过程中不随时间变化，即系统是时不变系统： 对于LTI系统，设系统的单位脉冲响应为h（n），则系统输入输出满足线性卷积关系： 还可以用一个N阶线性常系数差分方程来表示： 当a=0时，系统输入输出满足线性卷积的关系。 1.3 实验内容及步骤 1、编制信号产生、卷积、求解差分方程和显示子程序 设模拟正弦信号：，现以T为周期，对其进行采样，得到正弦序列为： 其中f为数字频率，w为数字角频率。 产生两个信号: N取大于500. 设两个离散时间线性时不变系统分别为： 系统1： 系统2： 其中： 2、实验内容 （1）令分别通过1、2（注意：通过系统1可以用线性卷积求解，通过系统2必须用差分方程求解），注意系统1、2输出信号的区别和共性。 Matlab程序： n=0:500; xa=sin(0.0364*pi*n); figure(1); stem(n,xa,.); xlabel(n);ylabel(xa(n)); hn(1:50)=1; y1=conv(xa,hn); figure(2); stem(y1,.); xlabel(n);ylabel(y1(n)); a0=1;a1=-1.7007;a2=0.7613; b0=0.02;b1=0.01;b2=0.0223; ys=1; B=[b0,b1,b2];A=[a0,a1,a2]; xi=filtic(B,A,ys); y2=filter(B,A,x1,xi); figure(3); stem(y2,.); xlabel(n);ylabel(y2(n)); 运行结果： 结论：单频信号通过系统1和系统2输出结果的频率相同。 统１存在过冲，系统2效果更好。? （2）令分别通过1、2，注意系统1、2输出信号的区别和共性与（1）有何不同之处。 Matlab程序： n=0:499; xb=sin(0.0364*pi*n)+sin(0.364*pi*n); figure(1); stem(n,xb,.); xlabel(n);ylabel(xb(n))； hn(1:50)=1; y1=conv(xb,hn); figure(2); stem(y1,.); a0=1;a1=-1.7007;a2=0.7613; b0=0.0223;b1=0.01;b2=0.0223; ys=1; B=[b0,b1,b2];A=[a0,a1,a2]; xi=filtic(B,A,ys); y2=filter(B,A,xb,xi); figure(3); stem(y2,.); 运行结果: 结论：1和f2，f1与频率相同，过系统1和系统2后，输出结果与过系统1和系统2类似，但过系统后有毛刺，曲线不光滑，且过系统2曲线比过系统1曲线光滑。说明，这两个系统相当于一个低通滤波器。 1.4 实验总结 通过本次实验离散时间信号通过线性系统的原理有了的了解。使用ＭＡＴＬＡＢ求解线性卷积的方法。单频率双频信号相同系统输出的的不同也有了更好的认识。