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

学生实验报告开课学院及实验室：电子楼3172013年4月8日学院机械与电气工程学院年级、专业、班姓名学号实验课程名称数字信号处理实验成绩实验项目名称实验三时域及频域采样定理指导老师一、实验目的熟悉并加深对采样定理的理解，了解采样信号的频谱和模拟信号频谱之间的关系。二、实验原理模拟信号经过理想采样，形成采样信号。采样信号的频谱和模拟信号频谱之间的关系如下：1()(())aasmXjXjmT此式告诉我们，采样信号的频谱是由模拟信号的频谱按照采样角频率周期性的延拓形成的，由此得到结论：采样频率必须大于等于模拟信号最高频率的两倍以上，才不会引起频率混叠。但用此式在计算机上进行计算不方便，下面我们将导出另外一个公式，以便在计算机上进行实验。对模拟信号()axt进行理想采样的公式如下式：()()()aanxtxttnT（3.1）对上式进行傅立叶变换，得到：()()tjaaXjxtedt[()()]jtanxttnTedt将上式的积分号和求和号交换次序，得到：()()()jtaanXjxttnTedt在上式的积分号内，只有当tnT时，才有非零值，因此：()()njTaanXjxnTe式中，()axnT在树值上等于由采样得到的时域离散信号()xn，如果再将T代入，得到：()()jnanXjxne（3.2）上式的右边就是序列的傅立叶变换()jwXe，即：()()|jwaTXjXe（3.3）上式说明采样信号的傅立叶变换可用相应序列的傅立叶变换得到，只要将自变量用T代替即可。这里有一个问题要解释，采样信号的频谱()aXj是将模拟信号的频谱按照采样角频频率为周期，进行周期性延拓形成的，而序列的傅立叶变换是以2为周期，这里是否一致？答案是肯定的。因为按照公式sffT/2，当sff时，2，因此序列的傅立叶变换以2为周期，转换到模拟域就是以采样频率sf为周期。另外，/2sff是()aXj的折叠频率，如果产生频率混叠，就是在该处附近发生，在数字域中，就是在附近易产生频谱混叠。有了以上的公式和概念，就可以用计算机研究对模拟信号的采样定理。下面分析频域采样定理。对信号x(n)的频谱函数()jXe，在[0，2π]上等间隔采样N点，得到2()(),0,1,2,,1jNkNXkXekN（3.4）则N点IDFT[()NXk]得到的序列就是原序列x(n)以N为周期进行周期延拓后的主值区序列，公式为：()IDFT[()][()]()NNNNixnXkxniNRn(3.5)由上式可知，频域采样点数N必须大于等于时域离散信号的长度M(即N≥M)，才能使时域不产生混叠，则N点IDFT[()NXk]得到的序列()Nxn就是原序列x(n),即()Nxn=x(n)。如果NM，()Nxn比原序列尾部多N-M个零点；如果NM，z则()Nxn=IDFT[()NXk]发生了时域混叠失真，而且()Nxn的长度N也比x(n)的长度M短，因此。()Nxn与x(n)不相同。在数字信号处理的应用中，只要涉及时域或者频域采样，都必须服从这两个采样理论的要点。对比上面叙述的时域采样原理和频域采样原理，得到一个有用的结论，这两个采样理论具有对偶性：“时域采样频谱周期延拓，频域采样时域信号周期延拓”。因此放在一起进行实验。三、使用仪器、材料1、硬件：计算机2、软件：Matlab四、实验步骤(一)时域采样定理实验1.给定模拟信号如下：0()sin()()ataxtAetut假设式中A=444.128,250a,2500rad/s，将这些参数代入上式中，对()axt进行傅立叶变换，得到()aXj，画出它的幅频特性()~aXjff，如图3.1所示。根据该曲线可以选择采样频率。图3.1()axt的幅频特性曲线2.按照选定的采样频率对模拟信号进行采样，得到时域离散信号()xn：0()()sin()()anTaxnxnTAenTunT这里给定采样频率如下：1sfkHz，300Hz，200Hz。分别用这些采样频率形成时域离散信号，按顺序分别用1()xn、2()xn、3()xn表示。选择观测时间50pTms。3.计算()xn的傅立叶变换()jwXe：100()[()]sin()iinanTjwjninXeFTxnAenTe（3.6）式中，1,2,3i，分别对应三种采样频率的情况123111(,,)1000300200TsTsTs。采样点数用下式计算：piiTnT（3.7）（3.6）式中，是连续变量。为用计算机进行数值计算，改用下式计算：100()[()]sin()ikiknjwanTjwnMinXeDFTxnAenTe（3.8）式中，2kkM，0,1,2,3...k，1M；64M。可以调用MATLAB函数fft计算3.8式。4.打印三种采样频率的幅度曲线()~kjwkXew，0,1,2,3...k，1M；64M。（二）频域采样定理实验给定信号如下：其它02614271301)(nnnnnx编写程序分别对频谱函数()FT[()]jXexn在区间]2,0[上等间隔采样32和16点，得到)()(1632kXkX和：32232()(),0,1,2,31jkXkXek16216()(),0,1,2,15jkXkXek再分别对)()(1632kXkX和进行32点和16点IFFT，得到)()(1632nxnx和：323232()IFFT[()],0,1,2,,31xnXkn161616()IFFT[()],0,1,2,,15xnXkn分别画出()jXe、)()(1632kXkX和的幅度谱，并绘图显示x(n)、)()(1632nxnx和的波形，进行对比和分析，验证总结频域采样理论。提示：频域采样用以下方法容易变程序实现。①直接调用MATLAB函数fft计算3232()FFT[()]Xkxn就得到()jXe在]2,0[的32点频率域采样②抽取32()Xk的偶数点即可得到()jXe在]2,0[的16点频率域采样16()Xk，即1632()(2),0,1,2,,15XkXkk。③当然也可以按照频域采样理论，先将信号x(n)以16为周期进行周期延拓，取其主值区（16点），再对其进行16点DFT(FFT),得到的就是()jXe在]2,0[的16点频率域采样16()Xk。五、实验过程原始记录（数据、图表、计算等）(一)时域采样定理实验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,n2,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);%产生采样序列x2=A*exp(-a*n2*T2).*sin(w0*n2*T2);x3=A*exp(-a*n3*T3).*sin(w0*n3*T3);Xk1=fft(x1,length(n1));%采样序列x1(n)的FFT变换Xk2=fft(x2,length(n2));%采样序列x2(n)的FFT变换Xk3=fft(x3,length(n3));%采样序列x3(n)的FFT变换k1=0:length(Xk1)-1;fk1=k1/Tp;%x1(n)的频谱的横坐标的取值k2=0:length(Xk2)-1;fk2=k2/Tp;%x2(n)的频谱的横坐标的取值k3=0:length(Xk3)-1;fk3=k3/Tp;%x3(n)的频谱的横坐标的取值subplot(3,2,1)stem(n1,x1,'.')title('(a)Fs=1000Hz');xlabel('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(Xk1))title('(a)FT[xa(nT)],Fs=1000Hz');xlabel('f(Hz)');ylabel('幅度');subplot(3,2,4)plot(fk2,abs(Xk2))title('(b)FT[xa(nT)],Fs=300Hz');xlabel('f(Hz)');ylabel('幅度')subplot(3,2,6)plot(fk3,abs(Xk3))title('(c)FT[xa(nT)],Fs=200Hz');xlabel('f(Hz)');ylabel('幅度')实验图像：（二）频域采样定理实验MATLAB代码M=32;N=27;n=0:1:N-1;xn=(n=0&n=13).*(n+1)+(n=14&n=26).*(27-n);%产生x(n)Xk=fft(xn,1024);%1024点FFT[x(n)]X32k=fft(xn,32);%32点FFT[x(n)]x32n=ifft(X32k);%32点IFFT[X32(k)]得到x32(n)m=0:15;%x(n)以16为周期作周期延拓并取主值区间，得到序列x16nx16n=(m=0&m=13).*(m+1)+(m=14&m=26).*(27-m)+(m=0&m=10).*(11-m);X16k=fft(x16n,16);k=0:1023;wk=2*k/1024;%连续频谱图的横坐标取值subplot(3,2,1);plot(wk,abs(Xk));title('(a)FT[x(n)]');xlabel('\omega/\pi');ylabel('|X(e^j^\omega)|');axis([0,1,0,200])subplot(3,2,2);stem(n,xn,'.');title('(b)三角波序列x(n)');xlabel('n');ylabel('x(n)');axis([0,32,0,20])k=0:M/2-1;%离散频谱图的横坐标取值subplot(3,2,3);stem(k,abs(X16k),'.');title('(c)16点频域采样');xlabel('k');ylabel('|X_1_6(k)|');axis([0,8,0,200])n1=0:M/2-1;subplot(3,2,4);stem(n1,x16n,'.');title('(d)16IDFT[X_1_6(k)]');xlabel('n');ylabel('x_1_6(n)');axis([0,32,0,20])k=0:M-1;%离散频谱图的横坐标取值subplot(3,2,5);stem(k,abs(X32k),'.');title('(e)32点频域采样');xlabel('k');ylabel('|X_3_2(k)|');axis([0,16,0,200])n1=0:M-1;subplot(3,2,6);stem(n1,x32n,'.');title('(f)32IDFT[X_3_2(k)]');xlabel('n');ylabel('x_3_2(n)');axis([0,32,0,20])实验图像：六、实验结果及分析1.由时域采样定理实验的图像可见，采样序列的频谱是以采样频率为周期对模