实验三-时域及频域采样定理.docxVIP

广州大学学生实验报告 开课学院及实验室:电子楼3 1 72013年 月 日 开课学院及实验室:电子楼3 1 7 2013年 月 日 学院 机电 年级、专 业、班 姓名 学号 实验课程名称 数字信号处理实验 成绩 实验项目名称 实验三时域及频域釆样定理 指导 老师 一、实验目的 熟悉并加深对采样左理的理解，了解采样信号的频谱和模拟信号频谱之间的关系。 二.实验原理 模拟信号经过理想采样，形成采样信号。采样信号的频谱和模拟信号频谱之间的关系如下: A 1 Xa（jc）= 〒Y x山g-〃g）） 此式告诉我们，采样信号的频谱是由模拟信号的频谱按照采样角频率周期性的延拓形成的，由此 得到结论：采样频率必须大于等于模拟信号最髙频率的两倍以上.才不会引起频率混叠。但用此式在 计算机上进行计算不方便，下面我们将导出另外一个公式,以便在讣算机上进行实验。 对模拟信号进行理想采样的公式如下式: 屁（/）=乙（0 屁（/）=乙（0工犯一 M） n—x (3.1) 对上式进行傅立叶变换，得到： 立心0）=匸 2（/）小屁/ 11—00 将上式的积分号和求和号交换次序，得到: 力（丿g）=丘匸以声“-皿尢-巴〃 在上式的积分号内,只有当t = nT时才有非零值，因此: 式中，石(T)在树值上等于由采样得到的时域藹散信号班介),如果再将(o = QT代入,得到： TOC \o 1-5 \h \z A x n-x (3.2) 上式的右边就是序列的傅立叶变换x2),即： ；L0G) = X(/)Is. (33) 上式说明采样信号的傅立叶变换可用相应序列的傅立叶变换得到,只要将自变量°用GT代替即可。 A 这里有一个问题要解释，采样信号的频谱X“OG)是将模拟信号的频谱按照采样角频频率为周 期，进行周期性延拓形成的，而序列的傅立叶变换是以2兀为周期，这里是否一致?答案是肯上的。因 为按照公式co = GT = 2吋［人，当了 =仁时,少=2因此序列的傅立叶变换以2兀为周期，转换到模 A 拟域就是以采样频率人为周期。另外，D2是X“OG)的折叠频率，如果产生频率混叠，就是 在该处附近发生,在数字域中，就是在“二兀附近易产生频谱混叠。有了以上的公式和槪念，就可以 用计算机研究对模拟信号的采样立理。 下而分析频域采样左理。对信号x (“)的频谱函数X(^),在［0 , 2上等间隔采样N点，得到 X”伙) = X(严)“％，斤=0,l,2,..?,N — l (3.4) 则N点I DFT［XA,伙)］得到的序列就是原序列班/7)以N为周期进行周期延拓后的主值区序列，公 式为： 兀vS WIDFTIXw 伙)h =［工 x(〃 + N)］/\) (3. 5) /—X 由上式可知,频域采样点数N必须大于等于时域离散信号的长度M(即NM),才能使时域不产生 混叠，则N点IDFT［Xn伙)］得到的序列心何就是原序列x (n),即心(”)=x(n)°如果N M,xjv(n) 比原序列尾部多N-M个零点；如果NM, z则xv(n)=IDFT［ X.v伙)］发生了时域混叠失真，而且 .0()的长度N也比x(n)的长度M短，因此。xN(n)与x(n)不相同。 在数字信号处理的应用中，只要涉及时域或者频域采样，都必须服从这两个采样理论的要点。对 比上而叙述的时域采样原理和频域采样原理，得到一个有用的结论，这两个采样理论具有对偶性：“时 域采样频谱周期延拓，频域采样时域信号周期延拓”。因此放在一起进行实验。 三、实验用MATLAB函数介绍 1 . fft 功能:一维快速傅立叶变换(FFT)。 Xk=fft(xn,N)：采用FFT算法计算序列向量xn的N点DFT,缺省N时，ffl函数自动按xn的 长度计算xn的DFT。当N为2的整数次幕时，fft按基2算法计算,否则用混合基算法。 2. i f f t 功能:一维快速逆傅立叶变换(IFFT)。 调用格式:与fft相同。 四、实验内容和步骤 (一)时域采样定理实验 1.给左模拟信号如下： 石⑴=Ae-巾 sing』)“(f) 假设式中A=444. 1 2。= 50龙、总，Qo = 50^2 rad/s,将这些参数代入上式中，对兀°) 进行傅立叶变换，得到X。(丿°),画出它的幅频特性 进行傅立叶变换，得到X。(丿°),画出它的幅频特性 X) ，如图3?1所示。根据该曲线可以选 择采样频率。 图3.1 图3.1心⑴的幅频特性曲线 实验过程及原始数据： elf： A=4 4 4 ? 128; a=50* pi* sqrt (2); w 0=50* pi* sqrt(2); f s= 1 000; 务采样频率10 0 OHZ T=l/fs; n=0 ： O?05*fs-l; %产生的长度区间 x t =A* e ：p (— a *n*T)? *