数字信号处理实验实验一.doc

数字信号处理实验报告 实验名称：离散时间系统的时域特性分析 学生姓名： z 学生学号： 学生班级： 上课时间： 周二上午 指导老师： 实验目的 线性时不变离散时间系统在时域中可以通过常系数线性差分方程来描述，冲激响应序列可以刻画其时域特性。本实验通过使用MATLAB函数研究离散时间系统的时域特性，以加深对离散时间系统的差分方程、冲激响应和系统的线性和时不变特性的理解。 实验原理 1.线性系统?满足叠加原理的系统称为线性系统，即若某一输入是由N个信号的加权和组成的，则输出就是系统对这几个信号中每一个输入的响应的加权和。?即：如果系统在x1（n）和x2（n）输入时对应的输出分别为y1（n）和y2（n），当对任意常数a1和a2，式T[a1x1(n)+a2x2(n)]=a1T[x1(n)]+a2[x2(n)]=a1y1(n)+a2y2(n)成立，则该系统是线性系统。?2.时不变系统?若输入x(n)的输出为y(n)，则将输入序列移动任意位后，其输出序列除了跟着位移外，数值应保持不变。即：当T[x(n)]=y(n)，满足T[x(n-m)]=y(n-m)?(m为任意整数)时，则该系统就称为时不变系统。?3.常系数线性差分方程?线性时不变离散系统的输入、输出关系可以用以下常系数线性差分描述：??????y(n)=-?∑aky(n-k)+?∑brx(n-r)?当输入x(n)为单位冲激序列时，输出y(n)即为系统的单位冲击响应h(n)。n=0:299; clf; n=0:299; x1=cos(20/256*pi*n); x2=cos(200/256*pi*n); x=x1+x2; num1=[0.5 0.27 0.77]; num2=[0.45 0.5 0.45]; den1=[1 0 0]; den2=[1 -0.53 0.46]; y1=filter(num1,den1,x1); y2=filter(num1,den1,x2); y=filter(num1,den1,x); yt=y1+y2; y3=filter(num2,den2,x1); y4=filter(num2,den2,x2); yn=filter(num2,den2,x); ytt=y3+y4; subplot(4,1,1); stem(n,x); ylabel(信号幅度); title(输入信号); subplot(4,1,2); stem(n,yt); ylabel(信号幅度); title(输出信号); subplot(4,1,3); stem(n,ytt); ylabel(信号幅度); title(输出信号); 系统的输入与输出波形如下图形所示： （2）解：MATLAB程序如下： N=300； num1=[0.5 0.27 0.77]; num2=[0.45 0.5 0.45]; den1=[1 0 0]; den2=[1 -0.53 0.46]; y1=impz(num1,den1,N); subplot(2,1,1); stem(y1); xlabel(时间序列); ylabel(信号幅度); title(输入信号);y1=impz(num1,den1,N); y2=impz(num2,den2,N); subplot(2,1,2); stem(y2); xlabel(时间序列); ylabel(信号幅度); title(输入信号); 系统1与系统2的冲激响应序列的波形如下图形所示： (3)ⅰ 解：MATLAB程序如下： n=0:40;D=5;x=cos(20*pi*n/256)+cos(200*pi*n/256);xd=[zeros(1,D)?x];?%生成新序列xd=x(n-D),延时D个单位num=[0.45?0.5?0.45];den=[1?-0.53?0.46];ic=[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(yd(1:length(yd)))?;ylabel(信号幅度)?;subplot(3,1,3)?;stem(n,d)?;xlabel(