信号及系统的谱分析.pdfVIP

数字信号处理实验一：信号及系统的谱分析 学号 姓名 注：1）此次实验作为《数字信号处理》课程实验成绩的重要依据，请同学们认真、独立完成，不得抄袭。 2 ）请在授课教师规定的时间内完成； 3 ）完成作业后，请以word 格式保存，文件名为：学号+姓名 4 ）请通读全文，依据第2 及第 3 两部分内容，认真填写第 4 部分所需的实验数据，并完成实验分析。 5 ）需将这次实验的内容给出一个纸质报告（31-40 号）。全体将报告的电子版交给班长以便实验结束 后刻成光盘 1. 实验目的 (1) 熟练利用 DFT 计算公式对信号进行谱分析， 加深 DFT 算法原理和基本性质的理解。 (2) 利用卷积方法计算信号经过离散系统输出响应，并观察输出信号的频谱变化。 (3) 熟悉 FFT 算法原理和 FFT 子程序的应用，掌握利用函数 fft.m 对离散信号及系统响应进行 频域分析。 (4) 理解并掌握利用 FFT 实现线性卷积的方法。了解可能出现的分析误差及其原因， 以便 在实际中正确应用 FFT 。 2. 实验原理与方法 1 ）离散傅里叶变换（DFT ）的基本原理 离散傅里叶变换（DFT ）是分析有限长序列频谱成分的重要工具，在信号处理的理 论上有重要意义。由于其可以在计算机上实现谱分析、 卷积、相关等主要的信号频谱分 析过程，因此 DFT 的快速算法得到了广泛的应用。 实现 DFT 的基本计算公式如下： 2 ）系统响应信号的时域分析（卷积运算） 离散信号输入离散系统后，若系统起始状态为 0 ，则系统的响应输出是 其方框图表示如下： 离散系统 h （n ） 图 1 在 matlab 中 计算卷积的函数为 y=conv(x,h) 。 3 ）FFT 实现线性卷积的快速计算 设一离散线性移不变系统的冲激响应为 ，长度为L 点；其输入信号为 ， 长度为 M 点；其输出为 ，长度为M+L-1 点。 当满足一定条件 时，有限长序列的线性卷积可用圆周卷积和来代替， 而圆周卷积可用 FFT 来计算，从而可以大大提高运算速度。 用 FFT 实现线性卷积计算的具体步骤： (1)有限长序列 和 补零值点,至长度为大于或等于 M+L-1 点，且 为 ， r 为整数。 (2)求 ，N 点 DFT ，用FFT 快速算法实现; (3)求 ，N 点 DFT ，用FFT 快速算法实现; (4)计算 ; (5)求 N 点 IDFT ，用IFFT 快速算法完成。 3. 实验内容及步骤 sin[0.15 (n  500) h(n)  R1001(n)  (n  500) 某系统的单位样值响应为： ， 信号 x （n）= （SIN （ω1n）+COS （ω2n））R1023 （n）， 输入该系统后，输出的响应信号为 y(n)。 请认真复习离散信号与系统、 线性卷积、 序列的傅里叶变换及性质等有关内容， 阅读上述实验原理与方法，编制 2 个程序文件完成如下 2 部分实验内容。 一） 利用函数 y=conv(x,h)求解响 开始 应信