数字信号处理C5上机实验_28511

1、数字信号处理上机实验数字信号处理上机实验主要内容主要内容实验一： 信号、 系统及系统响应实验二： 用FFT 作谱分析实验三： 用窗函 数法设计FIR数字滤波器实验一： 信号、 系统及系统响应实验目的实验目的(1) 熟悉连续信号经理想采样前后的频谱变化关系，熟悉连续信号经理想采样前后的频谱变化关系， 加深对时域采样定理的理解。加深对时域采样定理的理解。(2) 熟悉时域离散系统的时域特性。熟悉时域离散系统的时域特性。(3) 利用卷积方法观察分析系统的时域特性。利用卷积方法观察分析系统的时域特性。(4) 掌握序列傅里叶变换的计算机实现方法，掌握序列傅里叶变换的计算机实现方法， 利用利用序列的傅里叶变

2、换对连续信号、序列的傅里叶变换对连续信号、 离散信号及系离散信号及系统响应进行频域分析。统响应进行频域分析。 实验一： 信号、 系统及系统响应实验原理与方法实验原理与方法(1) 时域采样。时域采样。(2) LTI系统的输入输出关系。系统的输入输出关系。 实验内容及步骤实验内容及步骤(1) 认真复习采样理论、认真复习采样理论、 离散信号与系统、离散信号与系统、 线性线性卷积、卷积、 序列的傅里叶变换及性质等有关内容，序列的傅里叶变换及性质等有关内容， 阅读本实验原理与方法。阅读本实验原理与方法。(2) 编制实验用主程序及相应子程序。编制实验用主程序及相应子程序。 信号产生子程序，信号产生子程序，

3、 用于产生实验中要用用于产生实验中要用到的下列信号序列：到的下列信号序列： xa(t)=Ae-at sin(0t)u(t) 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) 有限长序列线性卷积子程序有限长序列线性卷积子程序 用于完成两个给定长度的序列的卷积。用于完成两个给定长度的序列的卷积。 可以直可以

4、直接调用接调用MATLAB语言中的卷积函数语言中的卷积函数conv。 conv用于两个有限长度序列的卷积，用于两个有限长度序列的卷积， 它假定它假定两个序列都从两个序列都从n=0 开始。开始。 调用格式如下：调用格式如下： y=conv (x, h)l调通并运行实验程序，调通并运行实验程序， 完成下述实验内容：完成下述实验内容： 分析采样序列的特性。分析采样序列的特性。 a. 取采样频率取采样频率fs=1 kHz, 即即T=1 ms。 b. 改变采样频率，改变采样频率， fs=300 Hz， 观察观察|X(ej)|的变化，的变化， 并做记录并做记录(打印曲线打印曲线)； 进一进一步降低采样频率

5、，步降低采样频率， fs=200 Hz， 观察频谱混叠观察频谱混叠是否明显存在，是否明显存在， 说明原因，说明原因， 并记录并记录(打印打印)这时这时的的|X(ej)|曲线。曲线。 时域离散信号、时域离散信号、 系统和系统响应分析。系统和系统响应分析。 a. 观察信号观察信号xb(n)和系统和系统hb(n)的时域和频域的时域和频域特性；特性； 利用线性卷积求信号利用线性卷积求信号xb(n)通过系统通过系统hb(n)的响应的响应y(n)， 比较所求响应比较所求响应y(n)和和hb(n)的时域及的时域及频域特性，频域特性， 注意它们之间有无差别，注意它们之间有无差别， 绘图说明，绘图说明， 并用所

6、学理论解释所得结果。并用所学理论解释所得结果。 b. 观察系统观察系统ha(n)对信号对信号xc(n)的响应特性。的响应特性。 卷积定理的验证。卷积定理的验证。 l思考题思考题 (1) 在分析理想采样序列特性的实验中，在分析理想采样序列特性的实验中， 采样频率不同时，采样频率不同时， 相应理想采样序列的傅里相应理想采样序列的傅里叶变换频谱的数字频率度量是否都相同叶变换频谱的数字频率度量是否都相同? 它们它们所对应的模拟频率是否相同所对应的模拟频率是否相同? 为什么为什么? (2) 在卷积定理验证的实验中，在卷积定理验证的实验中， 如果选如果选用不同的频域采样点数用不同的频域采样点数M值，值，

7、例如，例如， 选选M=10和和M=20， 分别做序列的傅里叶变换，分别做序列的傅里叶变换， 求得的求得的结果有无差异？为什么？结果有无差异？为什么？l实验报告要求实验报告要求 (1) 简述实验目的及实验原理。简述实验目的及实验原理。 (2) 按实验步骤附上实验过程中的信号序按实验步骤附上实验过程中的信号序列、列、 系统单位脉冲响应及系统响应序列的时系统单位脉冲响应及系统响应序列的时域和幅频特性曲线，域和幅频特性曲线， 并对所得结果进行分析并对所得结果进行分析和解释。和解释。 (3) 总结实验中的主要结论。总结实验中的主要结论。 (4) 简要回答思考题。简要回答思考题。实验二： 用FFT作谱分析

8、实验目的实验目的(1) 进一步加深进一步加深DFT算法原理和基本性质的理解算法原理和基本性质的理解(因因为为FFT只是只是DFT的一种快速算法，的一种快速算法， 所以所以FFT的运的运算结果必然满足算结果必然满足DFT的基本性质的基本性质)。(2) 熟悉熟悉FFT算法原理和算法原理和FFT子程序的应用。子程序的应用。(3) 学习用学习用FFT对连续信号和时域离散信号进行谱分对连续信号和时域离散信号进行谱分析的方法，析的方法， 了解可能出现的分析误差及其原因，了解可能出现的分析误差及其原因， 以便在实际中正确应用以便在实际中正确应用FFT。l实验步骤实验步骤 (1) 复习复习DFT的定义、的定义

9、、 性质和用性质和用DFT作作谱分析的有关内容。谱分析的有关内容。 (2) 复习复习FFT算法原理与编程思想，算法原理与编程思想， 并并对照对照DIT-FFT运算流图和程序框图，运算流图和程序框图， 读懂本读懂本实验提供的实验提供的FFT子程序。子程序。 (3) 编制信号产生子程序，编制信号产生子程序， 产生以下典产生以下典型信号供谱分析用：型信号供谱分析用： 1423( )( )1,03( )8470403( )3470 x nR nnnx nnnnnx nnn 其它n 其它n 456( )cos4( )sin8( )cos8cos16cos20 x nnx nnx nttt(4) 编写主程

10、序。编写主程序。下图给出了主程序框图，下图给出了主程序框图， 供参考。供参考。 本实验提本实验提供供FFT子程序和通用绘图子程序。子程序和通用绘图子程序。(5) 按实验内容要求，按实验内容要求， 上机实验，上机实验， 并写出实验并写出实验报告。报告。 读入长度N开始调用信号产生子程序产生实验信号调用绘图子程序(函数)绘制时间序列波形图调用FFT子程序(函数)计算信号的DFT调用绘图子程序(函数)绘制 X(k) 曲线结束l实验内容实验内容(1) 对对 2 中所给出的信号逐个进行谱分析。中所给出的信号逐个进行谱分析。 (2) 令令x(n)=x4(n)+x5(n)， 用用FFT计算计算 8 点和点和

11、 16 点离点离散傅里叶变换，散傅里叶变换， X(k)=DFTx(n)(3) 令令x(n)=x4(n)+jx5(n)， 重复重复(2)。l思考题思考题 (1) 在在N=8时，时， x2(n)和和x3(n)的幅频特性会相同吗的幅频特性会相同吗? 为什么为什么? N=16呢呢? (2) 如果周期信号的周期预先不知道，如果周期信号的周期预先不知道， 如何用如何用FFT进行谱分析进行谱分析?l实验报告要求实验报告要求 (1) 简述实验原理及目的。简述实验原理及目的。 (2) 结合实验中所得给定典型序列幅频特性曲线，结合实验中所得给定典型序列幅频特性曲线， 与理论结果比较，与理论结果比较， 并分析说明误

12、差产生的原因以及用并分析说明误差产生的原因以及用FFT作谱分析时有关参数的选择方作谱分析时有关参数的选择方法。法。 (3) 总结实验所得主要结论。总结实验所得主要结论。 (4) 简要回答思考题。简要回答思考题。实验三： 用窗函数法设计FIR数字滤波器实验目的实验目的(1) 掌握用窗函数法设计掌握用窗函数法设计FIR数字滤波器的原理和方数字滤波器的原理和方法。法。(2) 熟悉线性相位熟悉线性相位FIR数字滤波器特性。数字滤波器特性。(3) 了解各种窗函数对滤波特性的影响。了解各种窗函数对滤波特性的影响。l实验内容及步骤实验内容及步骤 (1) 复习用窗函数法设计复习用窗函数法设计FIR数字滤波器一

13、节数字滤波器一节内容，内容， 阅读本实验原理，阅读本实验原理， 掌握设计步骤。掌握设计步骤。 (2) 编写程序。编写程序。 编写能产生矩型窗、编写能产生矩型窗、 升余弦窗、升余弦窗、 改进升改进升余弦窗和二阶升余弦窗的窗函数子程序。余弦窗和二阶升余弦窗的窗函数子程序。 编写主程序。编写主程序。 其中幅度特性要求用其中幅度特性要求用dB表表示。示。 用窗函数法设计滤波器主程序框图 读入窗口长度N开始调用窗函数子程序求w(n)调用子程序(函数)计算H(k) DFTh(n)调用绘图子程序(函数)绘制H(k)幅度相位曲线结束计算hd(n)计算h(n) hd(n)w(n)上机实验内容上机实验内容 ,()01211( )()22sin()()ccj acjdcjj nj aj nddceHeNah nHeeeednana其中 l思考题思考题 (1) 如果给定通带截止频率和阻带截止频率如果给定通带截止频率和阻带截止频率以及阻带最小衰减，以及阻带最小衰减， 如何用窗函数法设计线性相如何用窗函数法设计线性相位低通滤波器位低通滤波器? 写出设计步骤。写出设计步骤。 (2) 如果要求用窗函数法设计带通滤波器，如果要求用窗函数法设计带通滤波器， 且给定上、且给定上、 下边带截止频率为下边带截止频率为1和和2，试求理，试求理想带通的单位脉