实验4时域采样理论与频域采样定理验证文档格式.doc

1、 （b）采样频率必须大于等于模拟信号最高频率的两倍以上，才能使采样信号的 频谱不产生频谱混叠。 利用计算机计算上式并不方便，下面我们导出另外一个公式，以便用计算机上进行实验。 理想采样信号和模拟信号之间的关系为： 对上式进行傅立叶变换，得到：在上式的积分号内只有当时，才有非零值，因此： 上式中，在数值上，再将代入，得到： 上式的右边就是序列的傅立叶变换，即 上式说明理想采样信号的傅立叶变换可用相应的采样序列的傅立叶变换得到，只要将自变量用代替即可。频域采样定理的要点是：a） 对信号x（n）的频谱函数X（ej）在0，2上等间隔采样N点，得到则N点IDFT得到的序列就是原序列x（n）以N为周期进行

2、周期延拓后的主值区序列，公式为： （b）由上式可知，频域采样点数N必须大于等于时域离散信号的长度M（即NM），才能使时域不产生混叠，则N点IDFT得到的序列就是原序列x（n）,即=x（n）。如果NM，比原序列尾部多N-M个零点；如果NM，z则=IDFT发生了时域混叠失真，而且的长度N也比x（n）的长度M短，因此。与x（n）不相同。 在数字信号处理的应用中，只要涉及时域或者频域采样，都必须服从这两个采样理论的要点。 对比上面叙述的时域采样原理和频域采样原理，得到一个有用的结论，这两个采样理论具有对偶性：“时域采样频谱周期延拓，频域采样时域信号周期延拓”。因此放在一起进行实验。三、实验内容及步骤（

3、1）时域采样理论的验证。给定模拟信号， 式中A=444.128，=50，=50rad/s，它的幅频特性曲线如图10.2.1 图10.2.1 的幅频特性曲线现用DFT（FFT）求该模拟信号的幅频特性，以验证时域采样理论。安照的幅频特性曲线，选取三种采样频率，即=1kHz，300Hz，200Hz。观测时间选。为使用DFT，首先用下面公式产生时域离散信号，对三种采样频率，采样序列按顺序用，表示。 因为采样频率不同，得到的，的长度不同， 长度（点数）用公式计算。选FFT的变换点数为M=64，序列长度不够64的尾部加零。X（k）=FFTx（n） ， k=0,1,2,3,-,M-1 式中k代表的频率为 。

4、要求： 编写实验程序，计算、和的幅度特性，并绘图显示。观察分析频谱混叠失真。（2）频域采样理论的验证。给定信号如下：编写程序分别对频谱函数在区间上等间隔采样32和16点，得到： 再分别对进行32点和16点IFFT，得到：分别画出、的幅度谱，并绘图显示x（n）、的波形，进行对比和分析，验证总结频域采样理论。提示：频域采样用以下方法容易变程序实现。 直接调用MATLAB函数fft计算就得到在的32点频率域采样 抽取的偶数点即可得到在的16点频率域采样，即。 当然也可以按照频域采样理论，先将信号x（n）以16为周期进行周期延拓，取其主值区（16点），再对其进行16点DFT（FFT）,得到的就是在的1

5、6点频率域采样。四、思考题如果序列x（n）的长度为M，希望得到其频谱在上的N点等间隔采样，当NM时， 如何用一次最少点数的DFT得到该频谱采样？五、实验报告及要求（1） 运行程序,打印要求显示的图形。（2） 分析比较实验结果，简述由实验得到的主要结论。 （3） 简要回答思考题。（4） 附上程序清单和有关曲线。六、程序清单和信号波形1、时域采样理论的验证程序清单：% 时域采样理论验证程序Tp=64/1000;%观察时间Tp=64微秒%产生M长采样序列x（n）% Fs=1000;T=1/Fs;Fs=1000;M=Tp*Fs;n=0:M-1;f=n*Fs/M;A=444.128;alph=pi*50

6、*20.5;omega=pi*50*20.5;xn=A*exp（-alph*n*T）.*sin（omega*n*T）;Xk=T*fft（xn,M）;%M点FFTxnt）subplot（3,1,1）;plot（f,abs（Xk）;xlabel（f/Hz）;ylabel（|x1（jf）|title（x1（n）的幅度特性%=%Fs=300HzFs=300;%Fs=200HzFs=200;xnt=A*exp（-alph*n*T）.*sin（omega*n*T）;Xk=T*fft（xnt,M）;yn=xa（nT）;subplot（3,2,5）;tstem（xnt,yn）;box on; title（a）

7、 Fs=1000Hzk=0:fk=k/Tp;subplot（3,2,6）;plot（fk,abs（Xk）;（a） T*FTxa（nT）,Fs=1000Hzf（Hz）幅度axis（0,Fs,0,1.2*max（abs（Xk）;信号波形：2、 频域采样理论的验证M=27;N=32;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）的TFX32k=fft（xn,32）;%32点FFTx（n）x32n=ifft（X32k）;%32点IFFTX3

8、2（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（b） 三角波序列x（n）nx（n）axis（0,32,0,20）1023;wk=2*k/1024;%subplot（3,2,1）;plot（wk,abs（Xk）;（a）FTx（n）omega/pi|X（ejomega）|axis（0,1,0,200）N/2-1;subplot（3,2,3）;stem（k,abs（X16k）,（c） 16点频域采样k|X_1

9、_6（k）|axis（0,8,0,200）n1=0:subplot（3,2,4）;stem（n1,x16n,（d） 16点IDFTX_1_6（k）x_1_6（n）N-1;stem（k,abs（X32k）,（e） 32点频域采样|X_3_2（k）|axis（0,16,0,200）stem（n1,x32n,（f） 32点IDFTX_3_2（k）x_3_2（n）思考题简答 先对原序列x（n）以N为周期进行周期延拓后取主值区序列，再计算N点DFT则得到N点频域采样：七、实验总结1由图可见，采样序列的频谱的确是以采样频率为周期对模拟信号频谱的周期延拓。当采样频率为1000Hz时频谱混叠很小；当采样频率为300Hz时，在折叠频率150Hz附近频谱混叠很严重；当采样频率为200Hz时，在折叠频率110Hz附近频谱混叠更很严重。2频域采样定理的图验证了频域采样理论和频域采样定理。对信号x（n）的频谱函数X（ej）在0，2上等间隔采样N=16时， N点IDFT得到的序列正是原序列x（n）以16为周期进行周期延拓后的主值区序列：