数字信号处理实验报告1.docxVIP

实验1、离散时间系统的时域特性分析姓名：倪烽 学号级周一678节实验目的线性时不变（LTI）离散时间系统在时域中可以通过常系数线性差分方程来描述，冲激响应列可以刻画时域特性。本次实验通过使用MATLAB函数研究离散时间系统的时域特性，以加深对离散时间系统的差分方程、冲激响应和系统的线性和时不变性的理解。基本原理一个离散时间系统是将输入序列变换成输出序列的一种运算。离散时间系统中最重要、最常用的是“线性时不变系统”。线性系统满足叠加原理的系统称为线性系统，即若某一输入是由N个信号的加权和组成的，则输出就是系统对这几个信号中每一个输入的响应的加权和。即那么当且仅当系统同时满足和时，系统是线性的。在证明一个系统是线性系统时，必须证明此系统同时满足可加性和比例性，而且信号以及任何比例系数都可以是复数。时不变系统系统的运算关系在整个运算过程中不随时间（也即序列的先后）而变化，这种系统称为时不变系统（或称移不变系统）。若输入的输出为，则将输入序列移动任意位后，其输出序列除了跟着位移外，数值应该保持不变，即则满足以上关系的系统称为时不变系统。常系数线性差分方程线性时不变离散系统的输入、输出关系可用以下常系数线性差分方程描述：当输入为单位冲激序列时，输出即为系统的单位冲激响应。当时，是有限长度的，称系统为有限长单位冲激响应（FIR）系统；反之，则称系统为无限长单位冲激响应（IIR）系统。4．计算系统单位冲激相应的两种函数impz和filter实验内容和要求实验内容考虑如下差分方程描述的两个离散时间系统：系统1：系统2：输入：编程求上述两个系统的输出，并画出系统的输入与输出波形。编程求上述两个系统的冲激响应序列，并画出波形。若系统的初始状态为零，判断系统2是否为时不变的？是否为线性的？实验结果及分析解释（1）编程求上述两个系统的输出，并画出系统的输入、输出波形。源代码：%系统1和2输入输出波形n=0:299;N=300;a1=[1];b1=[0.5 0.27 0.77];a2=[1 -0.53 0.46];b2=[0.45 0.5 0.45];x=cos(20.*n.*pi./256)+cos(200.*n.*pi./256);y1=impz(b1,a1,N);y2=impz(b2,a2,N);subplot(3,1,1);stem(n,x);title(系统输入波形);subplot(3,1,2);stem(n,y1);title(系统1输出波形);subplot(3,1,3);stem(n,y2);title(系统2输出波形);实验结果如右图所示，（2）编程求上述两个系统冲激响应序列，并画出波形。源代码：%系统1和2冲激响应序列波形n=0:49;N=50;a1=[1];b1=[0.5 0.27 0.77];y1=impz(b1,a1,N);subplot(2,1,1);stem(n,y1);title(系统1冲激响应序列波形);a2=[1 -0.53 0.46];b2=[0.45 0.5 0.45];y2=impz(b2,a2,N);subplot(2,1,2);stem(n,y2);title(系统2冲激响应序列波形);（3）若系统的初始状态为零，判断系统2是否为时不变的？是否为线性的？1.应用时延差值来判断系统2是否为时不变系统，即证明源代码：%系统2时不变验证clear allclcn=0:299;D=10;a=3;b=-2;x1=cos(2*pi*0.1*n);x2=cos(2*pi*0.4*n);x=a*x1+b*x2;xd=[zeros(1,D) x];%生成新序列，延时D个单位。xd长度比x长D=10num=[0.45 0.5 0.45];den=[1 -0.53 0.46];ic=[0 0]; %设置0初始条件y=filter(num,den,x,ic);yd=filter(num,den,xd,ic);N=length(y);d=y-yd(1+D:N+D); %求y和yd的差值subplot(3,1,1);stem(n,y);ylabel(信号幅度);title(输出y[n]);grid;subplot(3,1,2);stem(n,yd(1:length(yd-D)));%xd长度比x长10ylabel(信号幅度);title(由于输入延时而产生的输出y[n]);subplot(3,1,3);stem(n,d);xlabel(时间序号n);ylabel(信号幅度);title(差值信号);grid;实验结果如右图所示，从实验结果的第3幅图“差值信号”的幅度为0可以得出，系统2为时不变系统。2、在证明一个系统是线性系统时，必须证明此系统同时满足可加性和比例性：即当且仅