数字信号处理实验指导.doc

数字信号处理实验指导 实验一、离散时间系统的时域分析 （一）实验要求 掌握离散卷积计算方法； 学会差分方程的迭代解法； 了解全响应、零输入响应、零状态响应和初始状态的物理意义和具体求法； （二）实验内容 信号的加 数学描述 MATLAB实现 设[ x10=[1 0.7 0.4 0.1 0]; x20=[0.1 0.3 0.5 0.7 0.9 1];] 参考程序： n1=1:5; x10=[1 0.7 0.4 0.1 0]; n2=3:8; x20=[0.1 0.3 0.5 0.7 0.9 1]; n=1:8; x1=[x10 zeros(1,8-length(n1))]; x2=[zeros(1,8-length(n2)) x20]; x=x1+ x2; Subplot(3,1,1); Stem(n,x1); Subplot(3,1,2); Stem(n,x2); Subplot(3,1,3); Stem(n,x); 信号的乘 数学描述 MATLAB实现 设[ x10=[1 0.7 0.4 0.1 0]; x20=[0.1 0.3 0.5 0.7 0.9 1];] 参考程序： n1=1:5; x10=[1 0.7 0.4 0.1 0]; n2=3:8; x20=[0.1 0.3 0.5 0.7 0.9 1]; n=1:8; x1=[x10 zeros(1,8-length(n1))]; x2=[zeros(1,8-length(n2)) x20]; x=x1.*x2; Subplot(3,1,1); Stem(n,x1); Subplot(3,1,2); Stem(n,x2); Subplot(3,1,3); Stem(n,x); ==================================================================== n1=1:5; x10=[1 0.7 0.4 0.1 0]; n2=3:8; x20=[0.1 0.3 0.5 0.7 0.9 1]; n=1:8; x1=[x10 zeros(1,8-length(n1))]; x2=[zeros(1,8-length(n2)) x20]; x=x1.*x2 ===================================================================== 计算卷积 用MATLAB计算序列{-2 0 1 –1 3}和序列{1 2 0 -1}的离散卷积。 首先用手工计算，然后用MATLAB编程验证。 MATLAB参考程序如下： a=[-2 0 1 -1 3]; b=[1 2 0 -1]; c=conv(a,b); M=length(c)-1; n=0:1:M; stem(n,c); xlabel(n); ylabel(幅度); （三）实验报告要求 1、简述实验目的和实验原理，先人工计算求出你选定的序列x(n)、h(n)的卷积结果并与计算机计算结果相比较，给出笔算和计算机计算的差分方程结果对照表。 2、编写教材P26(7t)的运算的MATLAB程序。 实验二、 离散系统的Z域分析 （一）实验要求 学习和掌握离散系统的频率特性及其幅度特性、相位特性的物理意义； 深入理解离散系统频率特性的对称性和周期性； 认识离散系统频率特性与系统参数之间的关系； 通过阅读、修改并调试本实验系统所给源程序，加强计算机编程能力； （二）实验内容 1、计算差分方程 （1）用MATLAB计算差分方程 当输入序列为 时的输出结果 。 MATLAB程序如下： N=41; a=[0.8 -0.44 0.36 0.22]; b=[1 0.7 -0.45 -0.6]; x=[1 zeros(1,N-1)]; k=0:1:N-1; y=filter(a,b,x); stem(k,y) xlabel(n);ylabel(幅度) ????????请给出了该差分方程的前41个样点的输出，即该系统的单位脉冲响应。 2、 用MATLAB计算差分方程 所对应的系统函数的DTFT。 差分方程所对应的系统函数为： ?????????????????????????????? ??????? 其DTFT为 ?????????????????????????????? 用MATLAB计算的程序如下： k=256; num=[0.8 -0.44 0.36 0.02]; den=[1 0.7 -0.45 -0.6]; w=0:pi/k:pi; h=freqz(num,den,w); subplot(2,2,1); plot(w/pi,re