数字信号处理实验三时域及频域采样定理.docVIP

学生实验报告 开课学院及实验室： 电子楼317 2013 年 4 月 8 日 学院 机械与电气工程学院 年级、专业、班 姓名 学号 实验课程名称 数字信号处理实验 成绩 实验项目名称 实验三　时域及频域采样定理 指导老师 一、实验目的 熟悉并加深对采样定理的理解，了解采样信号的频谱和模拟信号频谱之间的关系。 二、实验原理 模拟信号经过理想采样，形成采样信号。采样信号的频谱和模拟信号频谱之间的关系如下： 此式告诉我们，采样信号的频谱是由模拟信号的频谱按照采样角频率周期性的延拓形成的，由此得到结论：采样频率必须大于等于模拟信号最高频率的两倍以上，才不会引起频率混叠。但用此式在计算机上进行计算不方便，下面我们将导出另外一个公式，以便在计算机上进行实验。 对模拟信号进行理想采样的公式如下式： 　　　　　　（3.1） 对上式进行傅立叶变换，得到： 将上式的积分号和求和号交换次序，得到： 在上式的积分号内，只有当时，才有非零值，因此： 式中，在树值上等于由采样得到的时域离散信号，如果再将代入，得到： 　　　　　　　　　（3.2） 上式的右边就是序列的傅立叶变换，即： 　　　　　　　　　（3.3） 上式说明采样信号的傅立叶变换可用相应序列的傅立叶变换得到，只要将自变量用代替即可。 这里有一个问题要解释，采样信号的频谱是将模拟信号的频谱按照采样角频频率为周期，进行周期性延拓形成的，而序列的傅立叶变换是以为周期，这里是否一致？答案是肯定的。因为按照公式，当时，，因此序列的傅立叶变换以为周期，转换到模拟域就是以采样频率为周期。另外，是的折叠频率，如果产生频率混叠，就是在该处附近发生，在数字域中，就是在附近易产生频谱混叠。有了以上的公式和概念，就可以用计算机研究对模拟信号的采样定理。 下面分析频域采样定理。对信号x(n)的频谱函数，在[0，2π]上等间隔采样N点，得到 （3.4） 则N点IDFT[]得到的序列就是原序列x(n)以N为周期进行周期延拓后的主值区序列，公式为： (3.5) 由上式可知，频域采样点数N必须大于等于时域离散信号的长度M(即N≥M)，才能使时域不产生混叠，则N点IDFT[]得到的序列就是原序列x(n),即=x(n)。如果NM，比原序列尾部多N-M个零点；如果NM，z则=IDFT[]发生了时域混叠失真，而且的长度N也比x(n)的长度M短，因此。与x(n)不相同。 在数字信号处理的应用中，只要涉及时域或者频域采样，都必须服从这两个采样理论的要点。对比上面叙述的时域采样原理和频域采样原理，得到一个有用的结论，这两个采样理论具有对偶性：“时域采样频谱周期延拓，频域采样时域信号周期延拓”。因此放在一起进行实验。 三、使用仪器、材料 1、硬件：计算机 2、软件：Matlab 四、实验步骤 (一) 时域采样定理实验 1.　给定模拟信号如下： 假设式中A=444.128, , rad/s，将这些参数代入上式中，对进行傅立叶变换，得到，画出它的幅频特性，如图3.1所示。根据该曲线可以选择采样频率。 图3.1　的幅频特性曲线 2.　按照选定的采样频率对模拟信号进行采样，得到时域离散信号： 这里给定采样频率如下：，300Hz，200Hz。分别用这些采样频率形成时域离散信号，按顺序分别用、、表示。选择观测时间。 3.　计算的傅立叶变换： 　　（3.6） 式中，，分别对应三种采样频率的情况。采样点数用下式计算： 　　　　　　　　　（3.7） （3.6）式中，是连续变量。为用计算机进行数值计算，改用下式计算： 　　　（3.8） 式中，，，；。可以调用MATLAB函数fft计算3.8式。 4.　打印三种采样频率的幅度曲线，，；。 （二）频域采样定理实验 给定信号如下： 编写程序分别对频谱函数在区间上等间隔采样32和16点，得到： 再分别对进行32点和16点IFFT，得到： 分别画出、的幅度谱，并绘图显示x(n)、的波形，进行对比和分析，验证总结频域采样理论。 提示：频域采样用以下方法容易变程序实现。 ① 直接调用MATLAB函数fft计算就得到在的32点频率域采样 ② 抽取的偶数点即可得到在的16点频率域采样，即。 ③ 当然也可以按照频域采样理论，先将信号x(n)以16为周期进行周期延拓，取其主值区（16点），再对其进行16点DFT(FFT),得到的就是在的16点频率域采样。 五、实验过程原始记录（数据、图表、计算等） (一) 时域采样定理实验MATLAB代码 A=444.128; a=50*pi*sqrt(2); w0=50*pi*sqrt(2); fs1=1000; fs2=300; fs3=2