2022年时域采样和频域采样文件 .pdf

一、实验目的时域采样理论与频域采样理论是数字信号处理中的重要理论。要求掌握模拟信号采样前后频谱的变化， 以及如何选择采样频率才能使采样后的信号不丢失信息；要求掌握频率域采样会引起时域周期化的概念，以及频率域采样定理及其对频域采样点数选择的指导作用。二、实验原理与方法时域采样定理的要点是：a. 对模拟信号)(txa以间隔 T进行时域等间隔理想采样，形成的采样信号的频谱)(?jX是原模拟信号频谱()aXj以采样角频率s（Ts/2）为周期进行周期延拓。公式为：)(?)(?txFTjXaa)(1nsajnjXTb. 采样频率s必须大于等于模拟信号最高频率的两倍以上，才能使采样信号的频谱不产生频谱混叠。利用计算机计算上式并不方便，下面我们导出另外一个公式，以便用计算机上进行实验。理想采样信号)(? txa和模拟信号)(txa之间的关系为：naanTttxtx)()()(?对上式进行傅立叶变换，得到：dtenTttxjXtjnaa )()()(?dtenTttxtjna)()(在上式的积分号内只有当nTt时，才有非零值，因此：nnTjaaenTxjX)()(?上式中，在数值上)(nTxa)(nx，再将T代入，得到：nnjaenxjX)()(?上式的右边就是序列的傅立叶变换)(jeX，即TjaeXjX)()(?名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 6 页 - - - - - - - - - 上式说明理想采样信号的傅立叶变换可用相应的采样序列的傅立叶变换得到，只要将自变量 用T代替即可。频域采样定理的要点是：a)对信号 x(n) 的频谱函数X(ej)在0， 2 上等间隔采样N 点，得到2( )() , 0,1,2,1jNkNXkX ekN则 N点 IDFT( )NXk 得到的序列就是原序列x(n)以 N 为周期进行周期延拓后的主值区序列，公式为：( )IDFT( )()( )NNNNixnXkx niNRnb)由上式可知，频域采样点数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）时域采样理论的验证给定模拟信号，)()sin()(0tutAetxta式中 A=444.128，=502，0=502rad/s ， 它的幅频特性曲线如图10.2.1 图 10.2.1 )(txa的幅频特性曲线名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 6 页 - - - - - - - - - 现用 DFT(FFT) 求该模拟信号的幅频特性，以验证时域采样理论。安照)(txa的幅频特性曲线，选取三种采样频率，即sF=1kHz，300Hz，200Hz。观测时间选msTp50。为使用 DFT，首先用下面公式产生时域离散信号，对三种采样频率，采样序列按顺序用)(1nx，)(2nx，)(3nx表示。)()sin()()(0nTunTAenTxnxnTa因为采样频率不同，得到的)(1nx，)(2nx，)(3nx的长度不同，长度（点数）用公式spFTN计算。选FFT 的变换点数为M=64 ，序列长度不够64 的尾部加零。X(k)=FFT x(n) ，k=0,1,2,3,-, M-1 式中 k代表的频率为kMk2。要求：编写实验程序，计算)(1nx、)(2nx和)(3nx的幅度特性，并绘图显示。观察分析频谱混叠失真。（2）频域采样理论的验证给定信号如下：其它02614271301)(nnnnnx编写程序分别对频谱函数()FT ( )jX ex n在区间2,0上等间隔采样32 和 16 点，得到)()(1632kXkX和：32232( )() , 0,1,2,31jkXkX ek1 6216( )() , 0,1,2,15jkXkX ek再分别对)()(1632kXkX和进行 32 点和 16 点 IFFT ，得到)()(1632nxnx和：3 23 23 2( )I F F T () ,0 , 1, 2 , 3 1xnXkn1 61 61 6( )I F F T () ,0 , 1, 2, 1 5xnXkn名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 6 页 - - - - - - - - - 分别画出()jX e、)()(1632kXkX和的幅度谱，并绘图显示x(n)、)()(1632nxnx和的波形，进行对比和分析，验证总结频域采样理论。提示：频域采样用以下方法容易变程序实现。1）直接调用MATLAB函数 fft 计算 就得到在 的 32 点频率域采样2）抽取的偶数点即可得到在 的 16 点频率域采样，即。3） 当然也可以按照频域采样理论，先将信号x(n)以 16 为周期进行周期延拓，取其主值区（16点） ，再对其进行16 点 DFT(FFT), 得到的就是在 的 16 点频率域采样。四、实验程序1 时域采样理论的验证程序清单% 时域采样理论验证程序exp2a.m Tp=64/1000; %观察时间Tp=64 微秒%产生 M 长采样序列x(n) % Fs=1000;T=1/Fs; 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); %M 点 FFTxnt) yn=xa(nT);subplot(3,2,1); tstem(xnt,yn); %调用自编绘图函数tstem 绘制序列图box on;title(a) Fs=1000Hz); k=0:M-1;fk=k/Tp; subplot(3,2,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) %= % Fs=300Hz 和 Fs=200Hz 的程序与上面Fs=1000Hz 完全相同。2 频域采样理论的验证程序清单%频域采样理论验证程序exp2b.m M=27;N=32;n=0:M; %产生 M 长三角波序列x(n) xa=0:floor(M/2); xb= ceil(M/2)-1:-1:0; xn=xa,xb; Xk=fft(xn,1024); %1024 点 FFTx(n), 用于近似序列x(n)的 TF X32k=fft(xn,32) ;%32 点 FFTx(n) x32n=ifft(X32k); %32 点 IFFTX32(k) 得到 x32(n) X16k=X32k(1:2:N); %隔点抽取X32k 得到 X16(K) x16n=ifft(X16k,N/2); %16 点 IFFTX16(k) 得到 x16(n) subplot(3,2,2);stem(n,xn,.);box on title(b) 三角波序列x(n);xlabel(n);ylabel(x(n);axis(0,32,0,20) 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 6 页 - - - - - - - - - k=0:1023;wk=2*k/1024; %subplot(3,2,1);plot(wk,abs(Xk);title(a)FTx(n); xlabel(omega/pi);ylabel(|X(ejomega)|);axis(0,1,0,200) k=0:N/2-1; subplot(3,2,3);stem(k,abs(X16k),.);box on title(c) 16 点频域采样 );xlabel(k);ylabel(|X_1_6(k)|);axis(0,8,0,200) n1=0:N/2-1; subplot(3,2,4);stem(n1,x16n,.);box on title(d) 16 点 IDFTX_1_6(k);xlabel(n);ylabel(x_1_6(n);axis(0,32,0,20) k=0:N-1; subplot(3,2,5);stem(k,abs(X32k),.);box on title(e) 32 点频域采样 );xlabel(k);ylabel(|X_3_2(k)|);axis(0,16,0,200) n1=0:N-1; subplot(3,2,6);stem(n1,x32n,.);box on title(f) 32 点 IDFTX_3_2(k);xlabel(n);ylabel(x_3_2(n);axis(0,32,0,20) 五、实验程序运行结果及分析1、 时域采样理论的验证程序运行结果exp2a.m 如图 10.3.2所示。由图可见，采样序列的频谱的确是以采样频率为周期对模拟信号频谱的周期延拓。当采样频率为1000Hz 时频谱混叠很小； 当采样频率为300Hz 时，在折叠频率150Hz附近频谱混叠很严重；当采样频率为 200Hz 时，在折叠频率110Hz附近频谱混叠更很严重。图 10.2.2 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 6 页 - - - - - - - - - 2 时域采样理论的验证程序exp2b.m 运行结果如图10.3.3 所示。图 10.3.3 该图验证了频域采样理论和频域采样定理。对信号 x(n)的频谱函数X(ej)在0，2 上等间隔采样 N=16 时， N 点 IDFT( )NXk 得到的序列正是原序列x(n)以 16 为周期进行周期延拓后的主值区序列：( )IDFT( )()( )NNNNixnXkx niNRn由于 NM ，频域采样定理，所以不存在时域混叠失真，因此。( )Nxn与 x(n)相同。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页，共 6 页 - - - - - - - - -