实验三时域及频域采样定理

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

2、对模拟信号进行理想采样的公式如下式：（3.1）对上式进行傅立叶变换，得到：将上式的积分号和求和号交换次序，得到：在上式的积分号内，只有当时，才有非零值，因此：式中，在树值上等于由采样得到的时域离散信号，如果再将代入，得到：（3.2）上式的右边就是序列的傅立叶变换，即：（3.3）上式说明采样信号的傅立叶变换可用相应序列的傅立叶变换得到，只要将自变量用代替即可。这里有一个问题要解释，采样信号的频谱是将模拟信号的频谱按照采样角频频率为周期，进行周期性延拓形成的，而序列的傅立叶变换是以为周期，这里是否一致？答案是肯定的。因为按照公式，当时，因此序列的傅立叶变换以为周期，转换到模拟域就是以采样频率为周期

3、。另外，是的折叠频率，如果产生频率混叠，就是在该处附近发生，在数字域中，就是在附近易产生频谱混叠。有了以上的公式和概念，就可以用计算机研究对模拟信号的采样定理。下面分析频域采样定理。对信号x(n)的频谱函数，在0，2上等间隔采样N点，得到 （3.4）则N点IDFT得到的序列就是原序列x(n)以N为周期进行周期延拓后的主值区序列，公式为： (3.5)由上式可知，频域采样点数N必须大于等于时域离散信号的长度M(即NM)，才能使时域不产生混叠，则N点IDFT得到的序列就是原序列x(n),即=x(n)。如果N>M，比原序列尾部多N-M个零点；如果N<M，z则=IDFT发生了时域混叠失真，而

4、且的长度N也比x(n)的长度M短，因此。与x(n)不相同。 在数字信号处理的应用中，只要涉及时域或者频域采样，都必须服从这两个采样理论的要点。对比上面叙述的时域采样原理和频域采样原理，得到一个有用的结论，这两个采样理论具有对偶性：“时域采样频谱周期延拓，频域采样时域信号周期延拓”。因此放在一起进行实验。三、实验用MATLAB函数介绍1.fft功能：一维快速傅立叶变换（FFT）。Xk=fft(xn,N)：采用FFT算法计算序列向量xn的N点DFT，缺省N时，fft函数自动按xn的长度计算xn的DFT。当N为2的整数次幂时，fft按基2算法计算，否则用混合基算法。2.ifft功能：一维快速逆傅立叶

5、变换（IFFT）。调用格式：与fft相同。四、实验内容和步骤(一) 时域采样定理实验1.给定模拟信号如下：假设式中A=444.128, , rad/s，将这些参数代入上式中，对进行傅立叶变换，得到，画出它的幅频特性，如图3.1所示。根据该曲线可以选择采样频率。图3.1的幅频特性曲线实验过程及原始数据：clf;A=444.128;a=50*pi*sqrt(2);w0=50*pi*sqrt(2);fs=1000; %采样频率1000HZT=1/fs; n=0:0.05*fs-1; %产生的长度区间nxt=A*exp(-a*n*T).*sin(w0*n*T); %产生采样序列xt(n)f=fft(x

6、t,length(n); %采样序列xt(n)的FFT变换k1=0:length(f)-1;fk1=k1/0.05; %设置xt(n)的频谱的横坐标的取值subplot(1,2,1)stem(n,xt,'.') %xt离散图title('(a)fs=1000Hz');xlabel('n');ylabel('x(n)');subplot(1,2,2)plot(fk1,abs(f)title('(a) FTx(nT),Fs=1000Hz');xlabel('f(Hz)');ylabel('幅度&

7、#39;)2.按照选定的采样频率对模拟信号进行采样，得到时域离散信号：这里给定采样频率如下：，300Hz，200Hz。分别用这些采样频率形成时域离散信号，按顺序分别用、表示。选择观测时间。实验过程及原始数据：clf;A=444.128;a=50*sqrt(2)*pi;w0=50*sqrt(2)*pi;Tp=50/1000;F1=1000;F2=300;F3=200; T1=1/F1;T2=1/F2;T3=1/F3; n1=0:0.5: Tp*F1-1;n2=0:0.5:Tp*F2-1;n3=0:0.5:Tp*F3-1; x1=A*exp(-a*n1*T1).*sin(w0*n1*T1); x2

8、=A*exp(-a*n2*T2).*sin(w0*n2*T2); x3=A*exp(-a*n3*T3).*sin(w0*n3*T3); f1=fft(x1,length(n1); f2=fft(x2,length(n2); f3=fft(x3,length(n3); k1=0:length(f1)-1;fk1=k1/Tp; k2=0:length(f2)-1;fk2=k2/Tp; k3=0:length(f3)-1;fk3=k3/Tp; subplot(3,2,1)stem(n1,x1,'.')title('(a)Fs=1000Hz');xlabel('

9、n');ylabel('x1(n)');subplot(3,2,3)stem(n2,x2,'.')title('(b)Fs=300Hz');xlabel('n');ylabel('x2(n)');subplot(3,2,5)stem(n3,x3,'.')title('(c)Fs=200Hz');xlabel('n');ylabel('x3(n)');subplot(3,2,2)plot(fk1,abs(f1)title('(a) FTx

10、a(nT)的频谱特性,Fs=1000Hz');xlabel('f(Hz)');ylabel('幅度')subplot(3,2,4)plot(fk2,abs(f2)title('(b) FTxa(nT) 的频谱特性,Fs=300Hz');xlabel('f(Hz)');ylabel('幅度')subplot(3,2,6)plot(fk3,abs(f3)title('(c) FTxa(nT) 的频谱特性,Fs=200Hz');xlabel('f(Hz)');ylabel('

11、;幅度') 3.计算的傅立叶变换：（3.6）式中，分别对应三种采样频率的情况。采样点数用下式计算：（3.7）（3.6）式中，是连续变量。为用计算机进行数值计算，改用下式计算：（3.8）式中，；。可以调用MATLAB函数fft计算3.8式。4.打印三种采样频率的幅度曲线，；。f=input('f= '); %输入取样频率f=1000T=1/f;n=0:50/(1000*T);M=64;%设置信号有关参数A=444.128; a=50*sqrt(2.0);w0=50*sqrt(2.0)*pi;xn=A*exp(-a*n*T).*sin(w0*n*T); % 产生x(n)Xk

12、=fft(xn,1024); %1024 点FFTx(n),用于近似序列 x(n)的FT X64k=fft(xn,M); %FFTx(n)x64n=ifft(X64k); %64 点IFFTX64(k) 得到64(n)subplot(2,2,1);stem(n,xn,'.');grid;box ontitle('(a) x1(n)');xlabel('n');ylabel('x1(n)');k=0:M-1;subplot(2,2,2);stem(k,abs(X64k),'.');grid;title('(b

13、) 64 点频域采样 输入取样频率f=1000');xlabel('k');ylabel('|X64(k)|');k=0:1023;wk=2*k/1024;subplot(2,2,3);plot(wk,abs(Xk);grid;title('(c)FTx1(n)');xlabel('omega/pi');ylabel('|X(ejomega)|');n1=0:M-1;subplot(2,2,4);stem(n1,x64n,'.');grid;box ontitle('(d) 64 点

14、')ifft(X64(k);xlabel('n');ylabel('x64(n)');f=input('f= '); % 输入取样频率f=300T=1/f;n=0:50/(1000*T);M=64;%设置信号有关参数A=444.128; a=50*sqrt(2.0);w0=50*sqrt(2.0)*pi;xn=A*exp(-a*n*T).*sin(w0*n*T); % 产生x(n)Xk=fft(xn,1024); %1024 点FFTx(n), 用于近似序列x(n)的FT X64k=fft(xn,M); %FFTx(n)x64n=ifft

15、(X64k); %64点 IFFTX64(k)得到 x64(n)subplot(2,2,1);stem(n,xn,'.');grid;title('(a) x2(n)');xlabel('n');ylabel('x2(n)');k=0:M-1;subplot(2,2,2);stem(k,abs(X64k),'.');grid;box ontitle('(b) 64 点频域采样输入取样频率f=300');xlabel('k');ylabel('|X64(k)|');k=

16、0:1023;wk=2*k/1024;subplot(2,2,3);plot(wk,abs(Xk);grid;title('(c)FTx2(n)');xlabel('omega/pi');ylabel('|X(ejomega)|');n1=0:M-1;subplot(2,2,4);stem(n1,x64n,'.');grid;box ontitle('(d) 64 点IFFTX64(k)');xlabel('n');ylabel('x64(n)');f=input('f= &

17、#39;); % 输入取样频率f=200T=1/f;n=0:50/(1000*T);M=64; %设置信号有关参数A=444.128; a=50*sqrt(2.0);w0=50*sqrt(2.0)*pi;xn=A*exp(-a*n*T).*sin(w0*n*T); %产生 x(n)Xk=fft(xn,1024); %1024 点FFTx(n), 用于近似序列x(n)的FT X64k=fft(xn,M); %FFTx(n)x64n=ifft(X64k); %64 点IFFTX64(k) 得到x64(n)subplot(2,2,1);stem(n,xn,'.');grid;box

18、ontitle('(a) x3(n)');xlabel('n');ylabel('x3(n)');k=0:M-1;subplot(2,2,2);stem(k,abs(X64k),'.');grid; title('(b) 64点频域采样 输入取样频率f=200');xlabel('k');ylabel('|X64(k)|');k=0:1023;wk=2*k/1024;subplot(2,2,3);plot(wk,abs(Xk);grid;title('(c)FTx3(n)&#

19、39;);xlabel('omega/pi');ylabel('|X(ejomega)|');n1=0:M-1;subplot(2,2,4);stem(n1,x64n,'.');grid;box ontitle('(d) 64 µã')ifft(X64(k);xlabel('n');ylabel('x64(n)'); 采样频率fs=1000采样频率fs=300采样频率fs=200（二）频域采样定理实验给定信号如下： 编写程序分别对频谱函数在区间上等间隔采样32和16点，得到： 再

20、分别对进行32点和16点IFFT，得到： 分别画出、的幅度谱，并绘图显示x(n)、的波形，进行对比和分析，验证总结频域采样理论。提示：频域采样用以下方法容易变程序实现。 直接调用MATLAB函数fft计算就得到在的32点频率域采样 抽取的偶数点即可得到在的16点频率域采样，即。 当然也可以按照频域采样理论，先将信号x(n)以16为周期进行周期延拓，取其主值区（16点），再对其进行16点DFT(FFT),得到的就是在的16点频率域采样。xn1=1:14;xn2=13:-1:1;xn=xn1,xn2;n=0:26;subplot(3,2,1);stem(n,xn,'.');grid

21、 ontitle('xn')xlabel('n');ylabel('xn');%绘制x(n)的频谱函数.xk=fft(xn,1024);%用1024点近似xn的谱函数N=0:1023;subplot(322);plot(N,abs(xk);grid ontitle('x(k)');xlabel('n');ylabel('x(k)');%取样32点xk32=fft(xn,32);xk16=xk32(1:2:32);%取样16个点xn32=ifft(xk32);xn16=ifft(xk16);%绘制xn32subplot(323);N=0:31;stem(N,xn32,'.');grid ontitle('xn32');xlabel('n');ylabel('xn32');%绘制xn32的谱函数subplot(324);plot(N,abs(xk32);grid ontitle('xk32');xlabel('n');ylabel('xk3