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

《数字信号处理--实验三 时域及频域采样定理》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、学生实验报告开课学院及实验室：电子楼3172013年4月8日学院机械与电气工程学院年级、专业、班姓名学号实验课程名称数字信号处理实验成绩实验项目名称实验三　时域及频域采样定理指导老师一、实验目的熟悉并加深对采样定理的理解，了解采样信号的频谱和模拟信号频谱之间的关系。二、实验原理模拟信号经过理想采样，形成采样信号。采样信号的频谱和模拟信号频谱之间的关系如下：此式告诉我们，采样信号的频谱是由模拟信号的频谱按照采样角频率周期性的延拓形成的，由此得到结论：采样频率必须大于等于模拟信号最高频率的两倍以上，才不会引起频率混叠。但用此式在计算机上进

2、行计算不方便，下面我们将导出另外一个公式，以便在计算机上进行实验。对模拟信号进行理想采样的公式如下式：　　　　　　（3.1）对上式进行傅立叶变换，得到：将上式的积分号和求和号交换次序，得到：在上式的积分号内，只有当时，才有非零值，因此：式中，在树值上等于由采样得到的时域离散信号，如果再将代入，得到：　　　　　　　　　（3.2）上式的右边就是序列的傅立叶变换，即：　　　　　　　　　（3.3）上式说明采样信号的傅立叶变换可用相应序列的傅立叶变换得到，只要将自变量用代替即可。这里有一个问题要解释，采样信号的频谱是将模拟信号的频谱按照采样角频

3、频率为周期，进行周期性延拓形成的，而序列的傅立叶变换是以为周期，这里是否一致？答案是肯定的。因为按照公式，当时，，因此序列的傅立叶变换以为周期，转换到模拟域就是以采样频率为周期。另外，是的折叠频率，如果产生频率混叠，就是在该处附近发生，在数字域中，就是在附近易产生频谱混叠。有了以上的公式和概念，就可以用计算机研究对模拟信号的采样定理。下面分析频域采样定理。对信号x(n)的频谱函数，在[0，2π]上等间隔采样N点，得到（3.4）则N点IDFT[]得到的序列就是原序列x(n)以N为周期进行周期延拓后的主值区序列，公式为：(3.5)由上式可

4、知，频域采样点数N必须大于等于时域离散信号的长度M(即N≥M)，才能使时域不产生混叠，则N点IDFT[]得到的序列就是原序列x(n),即=x(n)。如果N>M，比原序列尾部多N-M个零点；如果N

5、。三、使用仪器、材料1、硬件：计算机2、软件：Matlab四、实验步骤(一)时域采样定理实验1.　给定模拟信号如下：假设式中A=444.128,,rad/s，将这些参数代入上式中，对进行傅立叶变换，得到，画出它的幅频特性，如图3.1所示。根据该曲线可以选择采样频率。图3.1　的幅频特性曲线2.　按照选定的采样频率对模拟信号进行采样，得到时域离散信号：这里给定采样频率如下：，300Hz，200Hz。分别用这些采样频率形成时域离散信号，按顺序分别用、、表示。选择观测时间。3.　计算的傅立叶变换：　　（3.6）式中，，分别对应三种采样频率的

6、情况。采样点数用下式计算：　　　　　　　　　（3.7）（3.6）式中，是连续变量。为用计算机进行数值计算，改用下式计算：　　　（3.8）式中，，，；。可以调用MATLAB函数fft计算3.8式。4.　打印三种采样频率的幅度曲线，，；。（二）频域采样定理实验给定信号如下：编写程序分别对频谱函数在区间上等间隔采样32和16点，得到：再分别对进行32点和16点IFFT，得到：分别画出、的幅度谱，并绘图显示x(n)、的波形，进行对比和分析，验证总结频域采样理论。提示：频域采样用以下方法容易变程序实现。①直接调用MATLAB函数fft计算就得到

7、在的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=200;Tp=0.05;%观察时间Tp=50msn1=0:Tp*fs1-1;%产生不同的长度区间n1,n

8、2,n3n2=0:Tp*fs2-1;n3=0:Tp*fs3-1;T1=1/fs1;%不同的采样频率对应的采样周期T2=1/fs2;T3=1/fs3;x1=A*exp(-a*n1*T1).*sin(w0*n1*T1);%产