一、实验目的时域采样理论与频域采样理论是数字信号处理中的重要理论.要求掌握模拟信号采样前后频谱的变化,以与如何选择采样频率才能使采样后的信号不丢失信息；要求掌握频率域采样会引起时域周期化的概念,以与频率域采样定理与其对频域采样点数选择的指导作用.二、实验原理与方法时域采样定理的要点是：a.对模拟信号)(t x a 以间隔T 进行时域等间隔理想采样,形成的采样信号的频谱)(ˆΩj X是原模拟信号频谱()a X j Ω以采样角频率s Ω〔T s /2π=Ω〕为周期进行周期延拓.公式为： )](ˆ[)(ˆt x FT j X a a =Ω )(1∑∞-∞=Ω-Ω=n s a jn j X T b.采样频率s Ω必须大于等于模拟信号最高频率的两倍以上,才能使采样信号的频谱不产生频谱混叠.利用计算机计算上式并不方便,下面我们导出另外一个公式,以便用计算机上进行实验.理想采样信号)(ˆt xa 和模拟信号)(t x a 之间的关系为： 对上式进行傅立叶变换,得到：在上式的积分号内只有当nT t =时,才有非零值,因此：上式中,在数值上)(nT x a ＝)(n x ,再将T Ω=ω代入,得到：上式的右边就是序列的傅立叶变换)(ωj e X ,即 T j a e X j X Ω==Ωωω)()(ˆ上式说明理想采样信号的傅立叶变换可用相应的采样序列的傅立叶变换得到,只要将自变量ω用T Ω代替即可.频域采样定理的要点是：a) 对信号x<n>的频谱函数X<e j ω>在[0,2π]上等间隔采样N 点,得到则N 点IDFT[()N X k ]得到的序列就是原序列x<n>以N 为周期进行周期延拓后的主值区序列,公式为：b) 由上式可知,频域采样点数N 必须大于等于时域离散信号的长度M<即N ≥M>,才能使时域不产生混叠,则N 点IDFT[()N X k ]得到的序列()N x n 就是原序列x<n>,即()N x n =x<n>.如果N>M,()N x n 比原序列尾部多N-M 个零点；如果N<M,z 则()N x n =IDFT[()N X k ]发生了时域混叠失真,而且()N x n 的长度N 也比x<n>的长度M 短,因此.()N x n 与x<n>不相同.在数字信号处理的应用中,只要涉与时域或者频域采样,都必须服从这两个采样理论的要点. 对比上面叙述的时域采样原理和频域采样原理,得到一个有用的结论,这两个采样理论具有对偶性："时域采样频谱周期延拓,频域采样时域信号周期延拓〞.因此放在一起进行实验.三、 实验内容与步骤〔1〕时域采样理论的验证给定模拟信号,)()sin()(0t u t Ae t x t a Ω=-α式中A =444.128,α=502π,0Ω=502πrad/s,它的幅频特性曲线如图图)(t x a 的幅频特性曲线现用DFT<FFT>求该模拟信号的幅频特性,以验证时域采样理论.安照)(t x a 的幅频特性曲线,选取三种采样频率,即s F =1k Hz ,300Hz ,200Hz .观测时间选ms T p 50=.为使用DFT,首先用下面公式产生时域离散信号,对三种采样频率,采样序列按顺序用)(1n x ,)(2n x ,)(3n x 表示.因为采样频率不同,得到的)(1n x ,)(2n x ,)(3n x 的长度不同, 长度〔点数〕用公式s p F T N ⨯=计算.选FFT 的变换点数为M=64,序列长度不够64的尾部加零.X <k >=FFT[x <n >] , k =0,1,2,3,-----,M -1式中k 代表的频率为 k Mk πω2=. 要求： 编写实验程序,计算)(1n x 、)(2n x 和)(3n x 的幅度特性,并绘图显示.观察分析频谱混叠失真.〔2〕频域采样理论的验证给定信号如下：编写程序分别对频谱函数()FT[()]j X e x n ω=在区间]2,0[π上等间隔采样32和16点,得到)()(1632k X k X 和：再分别对)()(1632k X k X 和进行32点和16点IFFT,得到)()(1632n x n x 和：分别画出()j X e ω、)()(1632k X k X 和的幅度谱,并绘图显示x <n>、)()(1632n x n x 和的波形,进行对比和分析,验证总结频域采样理论.提示：频域采样用以下方法容易变程序实现.1〕直接调用MATLAB函数fft计算就得到在的32点频率域采样2〕抽取的偶数点即可得到在的16点频率域采样,即.3〕当然也可以按照频域采样理论,先将信号x<n>以16为周期进行周期延拓,取其主值区〔16点〕,再对其进行16点DFT<FFT>,得到的就是在的16点频率域采样 .四、实验程序1 时域采样理论的验证程序清单% 时域采样理论验证程序exp2a.mTp=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*2^0.5;omega=pi*50*2^0.5;xnt=A*exp<-alph*n*T>.*sin<omega*n*T>;Xk=T*fft<xnt,M>; %M点FFT[xnt>]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*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 频域采样理论的验证程序清单%频域采样理论验证程序exp2b.mM=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点FFT[x<n>], 用于近似序列x<n>的TFX32k=fft<xn,32> ;%32点FFT[x<n>]x32n=ifft<X32k>; %32点IFFT[X32<k>]得到x32<n>X16k=X32k<1:2:N>; %隔点抽取X32k得到X16<K>x16n=ifft<X16k,N/2>; %16点IFFT[X16<k>]得到x16<n>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>FT[x<n>]'>;xlabel<'\omega/\pi'>;ylabel<'|X<e^j^\omega>|'>;axis<[0,1,0,200]>k=0:N/2-1;subplot<3,2,3>;stem<k,abs<X16k>,'.'>;box ontitle<'<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 ontitle<'<d> 16点IDFT[X_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 ontitle<'<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 ontitle<'<f> 32点IDFT[X_3_2<k>]'>;xlabel<'n'>;ylabel<'x_3_2<n>'>;axis<[0,32,0,20]>五、实验程序运行结果与分析1、时域采样理论的验证程序运行结果exp2a.m如图所示.由图可见,采样序列的频谱的确是以采样频率为周期对模拟信号频谱的周期延拓.当采样频率为1000Hz时频谱混叠很小；当采样频率为300Hz时,在折叠频率150Hz附近频谱混叠很严重；当采样频率为200Hz时,在折叠频率110Hz附近频谱混叠更很严重.图2 时域采样理论的验证程序exp2b.m运行结果如图所示.图该图验证了频域采样理论和频域采样定理.对信号x<n>的频谱函数X<e jω>在[0,2π]上等间X k]得到的序列正是原序列x<n>以16为周期进行周期延拓隔采样N=16时,N点IDFT[()N后的主值区序列：x n与x<n>不相同,如图图<c>和<d>所示.当由于N<M,所以发生了时域混叠失真,因此.()NN=32时,如图图<c>和<d>所示,由于N>M,频域采样定理,所以不存在时域混叠失真,因x n与x<n>相同.此.()N。
