时域采样与频域采样 实验报告.doc

《时域采样与频域采样 实验报告.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、实验二时域采样与频域采样学校：西南大学班级：通信工程班一、实验目的时域采样理论与频域采样理论是数字信号处理中的重要理论。要求掌握模拟信号采样前后频谱的变化，以及如何选择采样频率才能使采样后的信号不丢失信息；要求掌握频率域采样会引起时域周期化的概念，以及频率域采样定理及其对频域采样点数选择的指导作用。二、实验原理时域采样定理的要点是采样频率必须大于等于模拟信号最高频率的两倍以上，才能使采样信号的频谱不产生频谱混叠。频域采样定理的要点是：a)对信号x(n)的频谱函数X(ejω)在[0，2π]上等间隔采样N点，得到则N点IDFT[]得到的序列就是原序列x(n)以N为周期进行周期延拓后的主值区序列

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

3、*T);Xk=T*fft(xnt,M);yn='xa(nT)';subplot(3,2,1);tstem(xnt,yn);boxon;title('(a)Fs=1000Hz');k=0:M-1;fk=k/Tp;subplot(3,2,2);plot(fk,abs(Xk));title('(a)T*FT[xa(nT)],Fs=1000Hz');xlabel('f(Hz)');ylabel('幅度');axis([0,Fs,0,1.2*max(abs(Xk))])（Fs=300Hz和Fs=200Hz的程序与上面Fs=1000Hz完全相同。）（2）频域采样理论的验证。M=27;N=32;n=0:

4、M;xa=0:floor(M/2);xb=ceil(M/2)-1:-1:0;xn=[xa,xb];Xk=fft(xn,1024);X32k=fft(xn,32);x32n=ifft(X32k);X16k=X32k(1:2:N);x16n=ifft(X16k,N/2);subplot(3,2,2);stem(n,xn,'.');boxontitle('(b)三角波序列x(n)');xlabel('n');ylabel('x(n)');axis([0,32,0,20])k=0:1023;wk=2*k/1024;subplot(3,2,1);plot(wk,abs(Xk));title('(a)

5、FT[x(n)]');xlabel('omega/pi');ylabel('

6、X(e^j^omega)

7、');axis([0,1,0,200])k=0:N/2-1;subplot(3,2,3);stem(k,abs(X16k),'.');boxontitle('(c)16点频域采样');xlabel('k');ylabel('

8、X_1_6(k)

9、');axis([0,8,0,200])n1=0:N/2-1;subplot(3,2,4);stem(n1,x16n,'.');boxontitle('(d)16点IDFT[X_1_6(k)]');xlabel('n');ylabel('x_

10、1_6(n)');axis([0,32,0,20])k=0:N-1;subplot(3,2,5);stem(k,abs(X32k),'.');boxontitle('(e)32点频域采样');xlabel('k');ylabel('

11、X_3_2(k)

12、');axis([0,16,0,200])n1=0:N-1;subplot(3,2,6);stem(n1,x32n,'.');boxontitle('(f)32点IDFT[X_3_2(k)]');xlabel('n');ylabel('x_3_2(n)');axis([0,32,0,20])四、实验结果①、时域采样理论的验证程序运行结果②、频

13、域采样理论的验证程序五、思考题先对原序列x(n)以N为周期进行周期延拓后取主值区序列，再计算N点DFT则得到N点频域采样：