《时域采样与频域采样》由会员分享，可在线阅读，更多相关《时域采样与频域采样（7页珍藏版）》请在金锄头文库上搜索。
1、时域采样与频域采样1. 实验目的时域采样理论与频域采样理论是数字信号处理中的重要理论。要求掌握模拟信号采样前后频谱的变化，以及如何选择采样频率才能使采样后的信号不丢失信息；要求掌握频率域采样会引起时域周期化的概念，以及频率域采样定理及其对频域采样点数选择的指导作用。2. 实验原理与方法时域采样定理的要点是：a) 对模拟信号 以间隔 T 进行时域等间隔理想采样，形成的采样信号的频谱)(txa是原模拟信号频谱 以采样角频率 （ ）为周)(jX()aXjsT/2期进行周期延拓。公式为：)()(txFTjaa )(1nsajnb) 采样频率 必须大于等于模拟信号最高频率的两倍以上，才能使采样信号的s频
2、谱不产生频谱混叠。利用计算机计算上式并不方便，下面我们导出另外一个公式，以便用计算机上进行实验。理想采样信号 和模拟信号 之间的关系为：)(txa)(txanaTttt对上式进行傅立叶变换，得到：dtenttxjXjnaa )()(tTtjna)() 在上式的积分号内只有当 时，才有非零值，因此：TtnnTjaaexjX)()(上式中，在数值上 ，再将 代入，得到：xannjaexj)()(上式的右边就是序列的傅立叶变换 ，即jX TjaeXj)()(上式说明理想采样信号的傅立叶变换可用相应的采样序列的傅立叶变换得到，只要将自变量 用 代替即可。T频域采样定理的要点是：a) 对信号 x(n)的
3、频谱函数 X(ej )在0，2 上等间隔采样 N 点，得到2() ,0,1,jNkNXkeL则 N 点 IDFT 得到的序列就是原序列 x(n)以 N 为周期进行周期延拓后的主值()区序列，公式为：()IDFT()()(NNNixnXkxniRb) 由上式可知，频域采样点数 N 必须大于等于时域离散信号的长度 M(即 NM)，才能使时域不产生混叠，则 N 点 IDFT 得到的序列 就是原序列 x(n),即()()Nxn=x(n)。如果 NM， 比原序列尾部多 N-M 个零点；如果 NM，z 则()Nxn()xn=IDFT 发生了时域混叠失真，而且 的长度 N 也比 x(n)的长()NXk()N
4、x度 M 短，因此。 与 x(n)不相同。x在数字信号处理的应用中，只要涉及时域或者频域采样，都必须服从这两个采样理论的要点。对比上面叙述的时域采样原理和频域采样原理，得到一个有用的结论，这两个采样理论具有对偶性：“时域采样频谱周期延拓，频域采样时域信号周期延拓” 。因此放在一起进行实验。3. 实验内容及步骤（1）时域采样理论的验证。给定模拟信号， )(sin)(0tuAetxta式中 A=444.128， =50 ， =50 rad/s，它的幅频特性曲线如图2210.2.1 图 10.2.1 的幅频特性曲线)(txa 现用 DFT(FFT)求该模拟信号的幅频特性，以验证时域采样理论。安照 的
5、幅频特性曲线，选取三种采样频率，即 =1kHz，300 Hz，200 Hz。)(txa sF观测时间选 。msTp50为使用 DFT，首先用下面公式产生时域离散信号，对三种采样频率，采样序列按顺序用 ， ， 表示。)(1nx2)(3nx)(si(0nTuAeTTa因为采样频率不同，得到的 ， ， 的长度不同， 长度（点数）)1x23x用公式 计算。选 FFT 的变换点数为 M=64，序列长度不够 64 的尾部加零。spFNX(k)=FFTx(n) ， k=0,1,2,3,-,M-1式中 k 代表的频率为 。2要求： 编写实验程序，计算 、 和 的幅度特性，并绘图显示。观察分析)(12nx)(3
6、频谱混叠失真。程序清单：（1）时域采样理论的验证：定义函数：function tstem(xn,yn)n=0:length(xn)-1;stem(n,xn,.);box onxlabel(n);ylabel(yn);axis(0,n(end),min(xn),1.2*max(xn)程序：Tp=50/1000;Fs=1000;T=1/Fs;M=Tp*Fs;n=0:M-1;A=444.128;alph=pi*50*20.5;omega=pi*50*20.5;xnt=A*exp(-alph*n*T).*sin(omega*n*T);Xk=T*fft(xnt,M); yn=xa(nT);subplot
7、(2,1,1);tstem(xnt,yn);box on;title(a) Fs=1000Hz);k=0:M-1;fk=k/Tp;subplot(2,1,2);plot(fk,abs(Xk);title(a) T*FTxa(nT),Fs=1000Hz);xlabel(f(Hz);ylabel(幅度);axis(0,Fs,0,1.2*max(abs(Xk)figure(2)Tp=50/1000;Fs=300;T=1/Fs;M=Tp*Fs;n=0:M-1;A=444.128;alph=pi*50*20.5;omega=pi*50*20.5;xnt=A*exp(-alph*n*T).*sin(ome
8、ga*n*T);Xk=T*fft(xnt,M); yn=xb(nT);subplot(2,1,1);tstem(xnt,yn);box on;title(b) Fs=300Hz);k=0:M-1;fk=k/Tp;subplot(2,1,2);plot(fk,abs(Xk);title(b) T*FTxb(nT),Fs=300Hz);xlabel(f(Hz);ylabel(幅度);axis(0,Fs,0,1.2*max(abs(Xk)figure(3)Tp=50/1000;Fs=200;T=1/Fs;M=Tp*Fs;n=0:M-1;A=444.128;alph=pi*50*20.5;omega=
9、pi*50*20.5;xnt=A*exp(-alph*n*T).*sin(omega*n*T);Xk=T*fft(xnt,M); yn=xc(nT);subplot(2,1,1);tstem(xnt,yn);box on;title(c) Fs=200Hz);k=0:M-1;fk=k/Tp;subplot(2,1,2);plot(fk,abs(Xk);title(c) T*FTxc(nT),Fs=200Hz);xlabel(f(Hz);ylabel(幅度);axis(0,Fs,0,1.2*max(abs(Xk)频域采样定理M=27;N=32;n=0:M;xa=0:floor(M/2); xb=
10、 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,.);box ontitle(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)FTx(n);xlabel(omega/pi);ylabel(|X(ej
11、omega)|);axis(0,1,0,200)n=0:N/2-1;subplot(3,2,4);stem(n,x16n,.);box ontitle(d) 16 点的 IDFTX16(k);xlabel(n);ylabel(X16(n);axis(0,32,0,20)k=0:N/2-1;subplot(3,2,3);stem(k,abs(X16k),.);title(c)16 点频域采样);xlabel(k);ylabel(|X16(k)|);axis(0,8,0,200)n=0:N-1;subplot(3,2,6);stem(n,x32n,.);box ontitle(f) 32 点的 IDFTX32(k);xlabel(n);ylabel(x(n);axis(0,32,0,20)k=0:N-1;subplot(3,2,5);stem(k,abs(X32k),.);title(e)32 点频域采样);xlabel(k);ylabel(|X32(k)|);axis(0,15,0,200)
