参考资料-dsp实验报告一.docVIP

实验一 信号系统及系统响应 一、实验目的 1、 熟悉理想采样的性质，了解信号采样前后的频谱变化，加深对采样定理的理解。 2、 熟悉离散信号和系统的时域特性。 3、 熟悉线性卷积的计算编程方法：利用卷积的方法，观察、分析系统响应的时域特性。 4、 掌握序列傅氏变换的计算机实现方法，利用序列的傅氏变换对离散信号、系统及系统响应进行频域分析。 二、实验原理 （一）连续时间信号的采样 采样是指按一定的频率从模拟信号抽样获得数字信号。采样是从连续时间信号到离散时间信号的过渡桥梁。 对一个连续时间信号进行理想采样的过程可以表示为该信号的一个周期冲激脉冲的乘积，即 （1） 其中连续信号的理想采样，是周期冲激脉冲 （2） 它也可以用傅立叶级数表示为： （3) 其中T为采样周期，Ω是采样角频率。设是连续时间信号的双边拉氏变换，即有： （4） 此时理想采样信号的拉氏变换为 （5） 作为拉氏变换的一种特例，信号理想采样的傅立叶变换 （6） 由式(5)和式(6)可知，信号理想采样后的频谱是原信号频谱的周期延拓，其延拓周期等于采样频率。根据Shannon取样定理，如果原信号是带限信号，且采样频率高于原信号最高频率分量的2倍，则采样以后不会发生频谱混淆现象。 计算机处理是利用序列的傅立叶变换计算信号的频谱，定义序列，根据Z变换的定义，可以得到序列x(n)的Z变换为： （7） 以代替上式中的z，就可以得到序列x(n)的傅立叶变换  （8） 式（6）和式（8）具有如下关系： （9） 由式（9）可知，在分析一个连续时间信号的频谱时，可以通过取样将有关的计算转化为序列傅立叶变换的计算。 （二）有限长序列分析 一般来说，在计算机上不可能，也不必要处理连续的曲线，通常，我们只要观察、分析在某些频率点上的值。对于长度为N的有限长序列 (10) 一般只需要在0?2π之间均匀地取M个频率点，计算这些点上的序列傅立叶变换 （11） 其中，=2πk/M,k=0,1,..,M-1。是一个复函数，它的模就是幅频特性曲线。（本实验一律采取M=50）。 （三）信号卷积 一个线性时不变离散系统的响应y(n)可以用它的单位冲激响应h(n)和输入信号x(n)的卷积来表示： (12) 根据傅立叶变换和Z变换的性质，与式（12）对应应该有 （13） （14） 式（12）告诉我们可以通过对两个序列的移位、相乘、累加计算信号响应；而式（14）告诉我们卷积运算也可以在频域上用乘积实现。 三、实验内容及步骤、结果分析 （一）编制实验用主程序及相应子程序 1、信号产生子程序： （1）理想采样信号序列：对信号进行理想采样，可以得到一个理想的采样信号序列：，其中A为幅度因子，α是衰减因子，Ω是频率。T为采样周期。 （2）单位脉冲序列 ??? （3）矩形序列 ，其中N=10 2、系统单位脉冲响应序列产生子程序，本实验中用到两种FIR系统： （1） （2） 3、有限长序列线性卷积子程序，用于计算两个给定长度（分别是M和N）的序列的卷积，输出序列长度为L=N+M-1。 程序见源代码。 （二）上机实验内容与结果分析 1、分析理想采样信号序列的特性。 产生理想采样信号序列，使A=128.444,α=50=Ω0。（1）首先选用采样频率为1000Hz，T=1/1000，观察所得理想采样信号的幅频特性，在折叠频率以内和给定的理想幅频特性无明显差异，并做记录；（2）改变采样频率为300Hz，T=1/300，观察所得到的频谱特性曲线的变化，并做记录；（3）进一步减小采样频率为200Hz，T=1/200，观察频谱“混淆”现象是否明显存在，说明原因，并记录这时候的幅频特性曲线。 结果分析： （1）采样频率为1000Hz时，信号、幅频、相频分别为： 可以看出在折叠频率500Hz以内和理想幅频特性无明显差异。 （2）采样频率为300Hz时，信号、幅频、相频分别为： 此时和（1）中比较起来频谱有明显泄漏。 （3）采样频率为200Hz时，信号、幅频、相频分别为： 上图表明当采样频率进一步降低时，频谱泄漏更为严重。 2、离散信号、系统和系统响应的分析 （1）观察信号xb（n）和系统hb(n)的时域和幅频特性；利用线性卷积求信号通过系统以后的响应。比较系统响应和系统h的时域及幅频特性。注意它们之间有无差异，绘出图形。 结果分