数字信号处理实验二时域采样和频域采样

1、实验二-时域采样和频域采样一、实验目的时域采样理论与频域采样理论是数字信号处理中的重要理论。要求掌握模拟 信号采样前后频谱的变化，以及如何选择采样频率才能使采样后的信号不丢失信 息；要求掌握频率域采样会引起时域周期化的概念，以及频率域采样定理及其对频域采样点数选择的指导作用。二、实验原理及方法1、时域采样定理的要点：a) 对模拟信号xa(t)以间隔t进行时域等间隔理想采样，形成的采样信号的频谱乂(j)是原模拟信号频谱xa(f 1 )以采样角频率"s (门s =2二/t)为周期进行周期延拓b) 采样频率门s必须大于等于模拟信号最高频率的两倍以上，才能使采样信号的 频谱不产生频谱混叠。利

2、用计算机计算上式并不方便，下面我们导出另外一个公 式，以便用计算机上进行实验。2、频域采样定理的要点：a)对信号x(n)的频谱函数x(ej)在0，2彳上等间隔采样n点则n点idft xn (k)得到的序列就是原序列x(n)以n为周期进行周期延拓后的 主值区序列。三、实验内容及步骤1、时域采样理论的验证 程序：clear;clca=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:tp*f1-1; n2=0:tp*f2-1; n3=0:tp*

3、f3-1;x1=a*exp(-a* n 1*t1).*si n(w0*in 1*t1);x2=a*exp(-a* n2*t2).*si n(w0* n2*t2);x3=a*exp(-a* n3*t3).*si n(w0* n3*t3);f1=fft(x1,le ngth( n1);f2=fft(x2,le ngth( n2);%f3=fft(x3,le ngth( n3);%k1=0:le ngth(f1)-1;fk1=k1/tp;%k2=0:le ngth(f2)-1;fk2=k2/tp;k3=0:le ngth(f3)-1;fk3=k3/tp;%subplot(3,2,1)stem( n1

4、,x1,'.')ti tle('(a)fs=1000hz');xlabel( 'n' );ylabel( 'x1( n)');subplot(3,2,3)stem( n2,x2,'.')ti tle('(b)fs=300hz');xlabel( 'n' );ylabel( 'x2( n)');subplot(3,2,5)stem( n3,x3,'.')ti tle('(c)fs=200hz');xlabel( 'n' )

5、;ylabel( 'x3( n)');subplot(3,2,2)plot(fk1,abs(f1)ti tle('(a) ftxa( nt),fs=1000hz');xlabel( 'f(hz)' );ylabel( ' u ? e)' subplot(3,2,4)plot(fk2,abs(f2)ti tle('(b) ftxa( nt),fs=300hz');xlabel( 'f(hz)' );ylabel( ' u ? e)' subplot(3,2,6)plot(fk3,abs

6、(f3)ti tle('(c) ftxa( nt),fs=200hz');xlabel( 'f(hz)' );ylabel( ' u ? e)'结果分析：由图 2.2 可见，采样序列的频谱的确是以采样频率为周期对模拟信号频谱的 周期延拓。当采样频率为1000hz时频谱混叠很小；当采样频率为 300hz时，在 折叠频率150hz附近频谱混叠很严重；当采样频率为200hz时，在折叠频率110hz 附近频谱混叠更很严重。由实验图像可以看出，时域非周期对应着频域连续。对连续时间函数对采样 使其离散化处理时，必须满足时域采样定理的要求， 否则，必将引起频域

7、的混叠。 要满足要求信号的最高频率fc不能采样频率的一半（fs/2），不满足时域采样定 理，频率将会在co = n附近，或者f=fs/2混叠而且混叠得最严重。 2、频域采样理论的验证 程序： clear;clcm=27;n=32;n=0: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,'.' );

8、boxontitle('(b)e y?2 - dod dx(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)ftx(n)');xlabel( 'omega/pi' );ylabel('|x(eajaomega)|');axis(0,1,0,200)k=0:n/2-1;subplot(3,2,3);stem(k,abs(x16k

9、),'.' );box ontitle( '(c)16 ? oo 2e ?)ix'bel('k' );ylabel( '|x_1_6(k)|');axis(0,8,0,200)n1=0:n/2-1;subplot(3,2,4);stem(n1,x16n,'.' );box ontitle( '(d)16?idftx_1_6(k)' );xlabel( 'n' );ylabel( 'x_1_6(n)');axis(0,32,0,20)k=0:n-1;subplot(3,

10、2,5);stem(k,abs(x32k),'.' );box ontitle( '(e)32 ? oo 2e ?)ix'bel('k' );ylabel( '|x_3_2(k)|');axis(0,16,0,200)n1=0:n-1;subplot(3,2,6);stem(n1,x32n,'.' );box ontitle( '(f)32 卩?idftx_3_2(k)' );xlabel( 'n' );ylabel( 'x_3_2(n)');axis(0,32,0,

11、20)结果分析：该图验证了频域采样理论和频域采样定理。对信号 x(n)的频谱函数x(ej co) 在0,2n止等间隔采样n=16时，n点idftxn(k)得到的序列正是原序列x( n) 以16为周期进行周期延拓后的主值区序列：0xn(n) =idft xn(k)h =' x(n 7n)rn(n)i由于n<m所以发生了时域混叠失真，因此，xn(n)与x(n)不相同，如图图3.3(c) 和(d)所示。当n=32时，如图图3.3(c)和(d)所示，由于n>m频域采样定理， 所以不存在时域混叠失真，因此，xn(n)与x(n)相同。由实验内容2的结果可知， 对一个信号的频谱进行采样处

12、理时，必须严格遵守频域采样定理，否则，用采样的离散频谱恢复原序列信号时，所得的时域离散序列是混叠失真，得不到原序列四、实验思考及解答如果序列x(n)的长度为m希望得到其频谱x(ej)在0,2二上的n点等间隔 采样，当n<m寸，如何用一次最少点数的dft得到该频谱采样？答：由实验内容2的结果可得：对于求频域采样点数n小于原时域序列长度 m的n点离散频谱时，可先对原序列x(n)以n为周期进行周期延拓后取主值区序 列0xn(n) = ' x(n in )rn(n)再计算n点dft则得到n点频域采样：xn(k) =dftxn(n)n =x(ej ) 2- , k =0,1,2, ,n -1xk但是，所求的n点离散频谱对应的时域离散序列是原序列x(n)以n为周期进行周期延拓后取主值区序列，而不是原序列 x(n)五、实验小结通过此次实验，对时域采样和频域采样的理论、定理的理解更加深入。采样