实验二_时域采样和频域采样.docxVIP

学习参考 学习参考 实验目的 时域采样理论与频域采样理论是数字信号处理中的重要理论 。要求掌握模 拟信号采样前后频谱的变化，以及如何选择采样频率才能使采样后的信号不丢 失信息；要求掌握频率域采样会引起时域周期化的概念 ，以及频率域采样定理 及其对频域采样点数选择的指导作用。 二、实验原理及方法 1、时域采样定理的要点： 对模拟信号Xa(t)以间隔T进行时域等间隔理想采样，形成的采样信号的 频谱刃(2)是原模拟信号频谱 Xa(j「)以采样角频率-s (门s=2-jy)为周 期进行周期延拓。公式为： 1 °0 涯(2)二 FT[?a(t)]二匚 Xa(j- - jn-s) 采样频率s必须大于等于模拟信号最高频率的两倍以上 ，才能使采样信 号的频谱不产生频谱混叠。利用计算机计算上式并不方便，下面我们导出另外 一个公式，以便用计算机上进行实验。 理想采样信号Xa(t)和模拟信号Xa(t)之间的关系为： Xa(t)二 Xa(t) 、(t - nT) n =joo 对上式进行傅立叶变换，得到： 乳(“)=.；[Xa(t)、(t- nT)]ejtdt - n oO「 oO 「Xa(t)、(t - nT)e jtdt 在上式的积分号内只有当t二nT时，才有非零值，因此： O0 . ?a（jC） =送 Xa （nT）e°—T n^joci 上式中，在数值上Xa（nT） = x（n）,再将」二T代入，得到： O0 兌a（j」）二 X（n）e-j n n 二.：: 上式的右边就是序列的傅立叶变换 X（eP ），即 刃a（j「）= x（ej） “ 上式说明理想采样信号的傅立叶变换可用相应的采样序列的傅立叶变换得到 只要将自变量①用T代替即可。 2、频域采样定理的要点： a）对信号x（n）的频谱函数X（ej?）在[0，2 n上等间隔采样N点，得到 XN（k）=X（e 气匚嗦,k=0,1,2川 |,N —1 则N点IDFT[XN（k）]得到的序列就是原序列x（n）以N为周期进行周期延拓后的 主值区序列，公式为： QO Xn（ n）= IDFTXn K n）P [ x n（ iN NR 口 i =q b）由上式可知，频域采样点数N必须大于等于时域离散信号的长度 M（即 N M），才能使时域不产生混叠，则N点IDFT[XN（k）]得到的序列XN（n）就是原 序列x（n）,即XN（n）=x（n）。如果nM，XN（n）比原序列尾部多N-M个零点；如果 NM，z则XN（n）=IDFT[ XN（k）]发生了时域混叠失真，而且XN（n）的长度n也比 x（n）的长度M短，因此。XN（n）与x（n）不相同。 在数字信号处理的应用中，只要涉及时域或者频域采样，都必须服从这两 个采样理论的要点。 对比上面叙述的时域采样原理和频域采样原理，得到一个有用的结论，这 两个采样理论具有对偶性：时域采样频谱周期延拓，频域采样时域信号周期延 拓”。因此放在一起进行实验。 三、实验内容及步骤 1、时域采样理论的验证 给定模拟信号，Xa(t)=Ae%in(0ot)u(t) 式中A=444.128，: =50 2 n，10=50 2冗rad/s，它的幅频特性曲线如图 2.1 图2.1 Xa(t)的幅频特性曲线 现用DFT(FFT求该模拟信号的幅频特性，以验证时域采样理论 按照Xa(t)的幅频特性曲线，选取三种采样频率，即Fs=1kHz，300Hz， 200Hz。观测时间选T 200Hz。观测时间选 Tp =50ms 。 为使用DFT，首先用下面公式产生时域离散信号，对三种采样频率，采样 % %观察时间Tp=50ms %不同的采样频率 %产生不同的长度区间n1,n2,n3 %产生采样序列x1(n) %产生采样序列x2(n) %产生采样序列x3(n) %采样序列x1(n)的FFT变换 %采样序列x2(n)的FFT变换 %采样序列x3(n)的FFT变换 序列按顺序用Xl(n) , x2(n) , x3(n)表示 x(n) = Xa(nT) = Ae七T sinQonT)u(nT) 因为采样频率不同，得到的Xl(n) , X2(n) , X3(n)的长度不同,长度(点 数)用公式N =Tp Fs计算。选FFT的变换点数为M=64，序列长度不够64的 尾部加零。 X(k)=FFT[x(n)]， k=0,1,2,3,-----,M-1 co k = — k 式中k代表的频率为 M 。 要求：编写实验程序，计算X1(n)、X2(n)和X3(n)的幅度特性，并绘图显 示 析频谱混叠失真。 Matlab源代码： A=444.128;a=50*sqrt(2)*pi;w0=50*sqrt(2)*pi; 邛=50/1000;F1=1000;F2=300;F3=200; T1=1/F1;T2=1/F2;T3=1