数字信号处理实验报告_7.docVIP

数字信号处理实验报告 专业班级 电信1101 姓 名 曾文 学 号 20111186020 指导老师 吴莉 华中科技大学武昌分校 2014年 5 月 21 日 实验一 信号、系统及系统响应 实验目的 (1)加深对离散线性移不变(LSI)系统基本理论的理解，明确差分方程与系统函数之 间的关系。 (2)初步了解用MATLAB语言进行离散时间系统研究的基本方法。 (3)掌握求解离散时间系统冲激响应和阶跃响应程序的编写方法，了解常用子函数。 (4)通过实验进一步理解卷积定理，了解卷积的过程。 (5)了解MATLAB中有关卷积的子函数及其应用方法。 实验原理 (1)离散LSI系统的响应与激励 由离散时间系统的时域和频域分析方法可知，一个线性移不变离散系统可以用线性常系数差分方程表示： 系统函数H(z)反映了系统响应与激励的关系。一旦上式中的bm 和ak 的数据确定了，则系统的性质也就确定了。其中特别注意：a0 必须进行归一化处理， 即a0 ＝1。 对于复杂信号激励下的线性系统，可以将激励信号在时域中分解为单位脉冲序列或单位阶跃序列，把这些单元激励信号分别加于系统求其响应，然后把这些响应叠加，即可得到复杂信号加于系统的零状态响应。因此，求解系统的冲激响应和阶跃响应尤为重要。由图1－1可以看出一个离散LSI系统响应与激励的关系。 同时，图1－1显示了系统时域分析方法和z变换域分析法的关系。如果已知系统的冲激响应h(n)，则对它进行z变换即可求得系统函数H(z)；反之，知道了系统函数H(z)，对其进行z逆变换，即可求得系统的冲激响应h(n)。 y(n)=x(n)*h(n)x(n)h(n) y(n)=x(n)*h(n) x(n) h(n) H(z) H(z) Y(z)=X(z)H(z) Y(z)=X(z)H(z) X(z) X(z) 图 1－1 离散LSI系统响应与激励的关系 图 1－1 离散LSI系统响应与激励的关系 (2) 离散LSI系统的线性卷积 由理论学习我们已知，对于线性移不变离散系统，任意的输入信号x(n)可以用 及其位移的线性组合来表示，即 当输入为 时，系统的输出y(n)=h(n)，由系统的线性移不变性质可以得到系统对x(n)的响应y(n)为 称为离散系统的线性卷积，简记为 也就是说，如果已知系统的冲激响应，将输入信号与系统的冲激响应进行卷积运算，即可求得系统的响应。 实验内容 (1)已知一个因果系统的差分方程为 满足初始条件y(-1)=0, x(-1)=0, 求系统的冲激响应和阶跃响应。 编写仿真程序，并调试得到结果，进行分析。 (2)已知两个信号序列： f1 = 0.8n (0n20) f2 = u(n) (0n10) 求两个序列的卷积和。 编写仿真程序，并调试得到结果，进行分析。 程序：①a=[1,0,1/3,0];%分母多项式系数 b=[1/6,1/2,1/2,1/6];%分子多项式系数 N=32; n=0:N-1;%一维数组，含N个分量，间隔为1 hn=impz(b,a,n);%单位冲激响应 gn=dstep(b,a,n);%单位阶跃响应 subplot(1,2,1),stem(n,hn,k); title(系统的单位冲激响应); ylabel(h(n));xlabel(n); subplot(1,2,2),stem(n,gn,k); title(系统的单位阶跃响应); ylabel(g(n));xlabel(n); ②nf1=0:20; f1=0.8.^nf1; subplot(2,2,1);stem(nf1,f1,filled); title(f1(n)); nf2=0:10; lf2=length(nf2); f2=ones(1,lf2); subplot(2,2,2);stem(nf2,f2,filled); title(f2(n)); y=conv(f1,f2); subplot(2,1,2);stem(y,filled);title(y(n)); 。 实验二 频域采样 1. 实验目的： (1) 掌握频率域采样会引起时域周期化的概念，以及频率域采样定理及其对频域采样点数选择的指导作用。 (2) 会用MATLAB语言进行频域抽样与恢复时程序的编写方法。 实验原理： 了解频域采样定理的要点，掌握采样理论的结论：“频域采样时域信号周期延拓”。 实验内容： （1）频域采样理论的验证。 给定信号如下： （2）编写程序分别对频谱函数在区间上等间隔采样32 和16点，得到，再分别对进行32点和16点IFFT，得到。