l《数字信号处理》实验.docVIP

姓名： 班级： 学号： 指导老师： 目录 目录 1 实验一 信号、 系统及系统响应 1 (一) 实验目的 1 (二) 实验原理与方法 2 (三) 实验内容及步骤 3 (四) 实验程序及运行结果 6 实验二： 用FFT作谱分析 14 (一) 实验目的 14 (二) 实验步骤 14 (三) 上机实验内容 15 (四) 程序实现及运行结果 15 实验三：设计IIR数字滤波器 25 (一) 实验目的 25 (二) 实验内容 25 (三) 实验步骤 26 (四) 程序实现及运行结果 27 实验四： 用窗函数法设计FIR数字滤波器 36 (一) 实验目的 36 (二) 实验原理与方法 36 (三) 实验内容及步骤 36 (四) 程序实现及运行结果 38 实验心得： 45 实验一 信号、 系统及系统响应 实验目的 1. 掌握常用序列。 2. 熟悉连续信号经理想采样前后的频谱变化关系， 加深对时域采样定理的理解。 3. 熟悉时域离散系统的时域特性。 4. 利用卷积方法观察分析系统的时域特性。 5. 掌握序列傅里叶变换的计算机实现方法， 利用序列的傅里叶变换对连续信号、 离散信号及系统响应进行频域分析。 实验原理与方法 采样是连续信号数字处理的第一个关键环节。 对采样过程的研究不仅可以了解采样前后信号时域特性和频域特性发生的变化以及信号信息不丢失的条件，而且可以加深对傅立业变换、Z变换和序列傅立业变换之间关系式的理解。 我们知道，对一个连续信号xa(t)进行理想采样的过程可用(1)式表示。 其中 为xa(t)的理想采样，p(t)为周期冲激脉冲， 即得到（２）式： 的傅里叶变换 为（３）式： 此式表明： 为X(jΩ)的周期延拓，其延拓周期为采样角频（）。只有满足采样定理时，才不会发生频率混叠失真。 现将（２）式带入（１）式并进行傅立叶变换的（４）式： 式中的xa(nT)就是采样后得到的序列x(n)，即 x(n)的傅里叶变换为（５）式： 比较(5)和(4)可知：得到（６）式 这说明两者之间只在频率量度上差一个常数因子Ｔ。 为了在数字计算机上观察分析各种序列的频域特性，通常对X(ejω)在［0,2π］上进行M点采样来观察分析。对长度为N的有限长序列x(n), 有（７）式： 其中 通常Ｍ要大于Ｎ． 一个时域离散线性非移变系统的输入/输出关系为（８）式 上述卷积运算也可以在频域实现 实验内容及步骤 1、编写常用信号的程序： （１）：单位采样序列δ（ｎ）。 （２）：单位阶跃序列ｕ（ｎ）。 （３）：单位矩形序列。 （4）：正弦序列。 （5）：。 （6）：。 2、认真复习采样理论、 离散信号与系统、 线性卷积、 序列的傅里叶变换及性质等有关内容， 阅读本实验原理与方法。 观察并分析采用不同频率时，对函数的频谱影响。 　　　　（ａ）：以，对其进行采样得到。 　　　　（ｂ）：以，对其进行采样得到。 （ｃ）：对、采用理想内插函数重建原始信号。 注：ＺＯＨ表示零阶保持内插，ｓｐｌｉｎｅ表示三次样条内插（）． 实现两序列的卷积。利用ｙ＝ｃｏｎｖ（ｘ，ｈ）．例如：求解 试求：． 在完成了（１）、（２）的基础上，编写（３）： 　　（３）　编制实验用主程序及相应子程序。 ① 信号产生子程序， 用于产生实验中要用到的下列信号序列： xa(t)=Ae-at sin(Ω0t)u(t) 进行采样， 可得到采样序列 xa(n)=xa(nT)=Ae-anT sin(Ω0nT)u(n), 0≤n50 其中A为幅度因子， a为衰减因子， Ω0是模拟角频率， T为采样间隔。 这些参数都要在实验过程中由键盘输入， 产生不同的xa(t)和xa(n)。 b. 单位脉冲序列： xb(n)=δ(n) c. 矩形序列： xc(n)=RN(n), N=10 ② 系统单位脉冲响应序列产生子程序。 本实验要用到两种FIR系统。 a. ha(n)=R10(n); b. hb(n)=δ(n)+2.5δ(n-1)+2.5δ(n-2)+δ(n-3) ③ 有限长序列线性卷积子程序， 用于完成两个给定长度的序列的卷积。 可以直接调用MATLAB语言中的卷积函数conv。 conv用于两个有限长度序列的卷积， 它假定两个序列都从n=0 开始。 调用格式如下： y=conv (x, h) 图10.3.1 实验一的主程序框图 图10.3.2 xa(t)的幅频特性