数字信号处理实验报告-五个实验.docx

《数字信号处理实验报告-五个实验.docx》由会员分享，可在线阅读，更多相关《数字信号处理实验报告-五个实验.docx（36页珍藏版）》请在优知文库上搜索。

1、教字信号处理实验报告实验一信号、系统及系统响应一、实验目的1、熟悉连续信号经理想采样前后的频谱变化关系，加深对时域采样定理的理解;2、熟悉时域离散系统的时域特性；3、利用卷积方法观察分析系统的时域特性；4、掌握序列傅立叶变换的计算机实现方法，利用序列的傅立叶变换对连续信号、离散信号及系统响应进行频域分析。二、实验原理及方法采样是连续信号数字处理的第一个关键环节。对采样过程的研究不仅可以了解采样前后信号时域和频域特性发生变化以及信号信息不丢失的条件，而且可以加深对傅立叶变换、Z变换和序列傅立叶变换之间关系式的理解。对一个连续信号与进行理想采样的过程可用下式表示：JCa)=JGa)P其中K为的理想

2、采样，伙。为周期脉冲，即P=(t-nT)W=-OO)的傅立叶变换为1X4(4)=7Xg(QTmJm=-上式表明XgC)为XtS网的周期延拓。其延拓周期为采样角频率g=2T).只有满足采样定理时，才不会发生频率混叠失真。在实验时可以用序列的傅立叶变换来计算X“C/Q)。公式如下：.g)=X()离散信号和系统在时域均可用序列来表示。为了在实验中观察分析各种序列的频域特性，通常对X(W)在0,21上进行M点采样来观察分析。对长度为N的有限长序列x(n),有：NTX(e)=Zx(m)e-M=O其中，k=一k,k=0,l,M-IM时域离散线性非移变系统的输入/输出关系为y(n)=x(n)*h(n)=Zx

3、(m)h(n-ni)/I=-OO上述卷积运算也可在频域实现Y(ej)=X(ej)H(ej)三、实验程序s=yesinput(PleaseSelectTheStepOfExperiment:n一.(1时域采样序列分析s=str2num(s);closeall;Xb=impseq(O,0,1);Ha=stepseq(l,1,10)；Hb=impseq(O,O,3)+2.5*impseq(l,O,3)+2.2*impseq(2,O,3)+impseq(3,0,3);i=0;while(s);%时域采样序列分析if(s=l)1=1；k=0;while(1)if(k=0)A=yesinput(,plea

4、seinputtheAmplitude:n,.444.128,100,100O);a=yesinput(,pleaseinputtheAttenuationCoefficientAn222.144,100,600);w=yesinput(,pleaseinputtheAngleFrequence(rads):n,.222.144,100,600);endk=k+l;fs=yesinput(,pleaseinputthesamplefrequence:n,.1000,100,1200);Xa=FF(A,a,w,fs);i=i+l;string+,fs三,num2str(fs);figure(i)

5、DFT(Xa,50,string);l=yesinputl=str2num(l);end%系统和响应分析elseif(s=2)kk=str2num(kk);while(kk)if(kk=l)m=conv(Xb,Hb);N=5;i=i+l;figure(i)string=(,hb(n),);Hs=DFT(Hb,4,string);i=i+l;figure(i)stringCxb(n),);DFT(Xb,2,string);string=(,y(n)=xb(n)*hb(n),);elseif(kk=2)m=conv(Ha,Ha);N=19;string=(,y(n)=ha(n)*(ha(n),);

6、elseif(kk=3)Xc=stepseq(l,1,5);m=conv(Xe,Ha);N=14;string=(,y(n)=xc(n)*ha(n),);endendendi=i+l;figure(i)DFT(m,N,string);kk=yesinputkk=str2num(kk);end卷积定理的验证elseif(s=3)A=I;a=0.5;w=2,0734;fs=l;Xal=FF(A,a,w,fs);i=i+l;figure(i)string=(,Thexal(n)(A=l,a=0.4,T=l),);Xa,wDFT(Xal,50,string);i=i+l;figure(i)string

7、=(,hb(n),);Hs=DFT(Hb,4,string);Ys=Xs.*Hs;y=conv(Xal,Hb);N=53;i=i+l;figure(i)string=(,y(n)=xa(n)*hb(n),);yy,w=DFT(y,N,string);i=i+l;figure(i)subplot(2,2,1)plot(wpi,abs(yy);axis(-2202);xlabel(,wpi,);ylabel(,Ys(jw),);title(FTx(n)*h(n)J,);subplot2,3)plot(wpi,abs(Ys);axis(-2202);xlabel(,wpi,);ylabelCYs(j

8、w),);title(,FTxs(n).FTh(n)J,);endendend子函数:离散傅立叶变换及X(n),FTx(n)的绘图函数functionc,l=DFT(x,N,str)n=0:N-l;k=-200:200;w=(pi100)*k;l=w;c=x*Xc=stepseq(l,1,5);子函数：产生信号functionc=FF(A,a,w,fs)n=o:50-1;c=A*exp(-a)*nfs).*sin(w*nfs).*stepseq(0,0,49);子函数：产生脉冲信号functionx,n=impseq(n,nl,n2)n=nl:n2;x=(11-nO)=O;子函数：产生矩形框信

9、号functionx,n=Stepseq(n,nl,n2)n=nl:n2;x=(n-nO=O);四、实验内容及步骤1、认真复习采样理论，离散信号与系统，线性卷积，序列的傅立叶变换及性质等有关内容，阅读本实验原理与方法。2、编制实验用主程序及相应子程序。3、调通并运行实验程序，完成下列实验内容：分析采样序列的特性，产生采样信号序列乙，使A=444.128,a=502,0=502-o(xa(t)的无失真采样频率约为1000Hz)。a. 取采样频率f=lkHz,即T=Ims。观察所得采样儿()的幅频特性X(ftl)|和原图中的幅频特性曲线在折叠频率附近有无明显差别。应当注意，实验中所得频谱是用序列的

10、傅立叶变换公式求得的，所以在频率度量上存在关系：=,为数字频率，Q为模拟频率。b. 改变采样频率，=300Hz,观察X(e加)|的变化，并做记录(打印曲线)；进一步降低采样频率，s=200Hz,观察频谱混叠是否明显存在，说明原因，并记录(打印)这时的X()曲线。由图形可知，当采样频率为IOooHZ时，采样序列在折叠频率附近处，即=乃处无明显频谱混叠。b步实验结果如下图所示:信号的原形由图可知，当采样频率进一步降低时，主瓣宽度逐渐变宽，频率混叠现象也逐渐严重。存在较明显的失真现象。时域离散信号、系统和系统响应分析。a.观察信号项,()和系统/5)的时域和频域特性；利用线性卷积求信号45)通过系统

11、45)的响应y(n),比较所求响应y(n)和儿()的时域及频域特性，注意它们之间有无差别。绘图说明，并用所学理论解释所得结果。实验结果如下图所示:b.观察系统4()对信号z5)的响应特性。利用线性卷积求系统响应y(n),并判断y(n)图形及其非零值序列长度是否与理论结果一致，对乙()=4()=Rn)()，说出一种定性判断y(n)图形正确与否的方法。调用序列傅立叶变换数值计算子程序，求得丫(0网)，观察|丫(/叫)|特性曲线，定性判断结果的正确性。改变/()的长度，取N=5,重复该实验。注意参数变化的影响,说明变化前后的差异，并解释所得结果。实验结果如下图所示:N=IO:N=5：欲判断结果正确与

12、否，可以先对其进行运算，算出其卷积，再与图形对照。当N=Io时，峰值较高，且峰值很窄，变换之后图形频带主值部分比较集中；N=5时情况与之相反。卷积定理的验证。将实验中的信号换为5）,使a=0.4,。0=2,0734,A=I,T=I,重复实验a,打印|丫（，四）|曲线；对主程序做简单修改，按式（1039）计算y（/W）=X.（eM并绘出y（e网）|曲线，与前面直接对y（n）进行傅立叶变换所得幅频特性曲线进行比较，验证时域卷积定理。实验所得结果如下:FTx(n)h(n)四、思考题1、在分析理想采样序列特性的实验中，采样频率不同，相应理想采样序列的傅立叶变换频谱的数字频率度量是否都相同？它们所对应的

13、模拟频率是否相同？为什么？答：由0=T可知，若采样频率不同，则其周期T不同，相应的数字频率。也不相同；而因为是同一信号，故其模拟频率。保持不变。2、在卷积定理验证的实验中，如果选用不同的频域采样点数M值，例如，选M=IO和M=20,分别做序列的傅立叶变换，求得Y(ejt)=Xa(ejt)Hh(ejt),k=O,l,M-l所得结果之间有无差异？为什么？答：有差异。因为所得y(m)图形由其采样点数唯一确定，由频域采样定理可知，若M小于采样序列的长度N,则恢复原序列时会发生时域混叠现象。实验二用FFT作谱分析一、实验目的1、进一步加深DFT算法原理和基本性质的理解(因为FFT只是DFT的一种快速算法

14、，所以FFT的运算结果必然满足DFT的基本性质)。2、熟悉FFT算法原理和FFT子程序的应用。3、学习用FFT对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析误差及其原因，以便在实际中正确应用FFT。二、实验步骤1、复习DFT的定义、性质和用DFT作谱分析的有关内容。2、复习FFT算法原理与编程思想，并对照DITFFT运算流图和程序框图,读懂本实验提供的FFT子程序。3、编制信号产生子程序，产生以下典型信号供谱分析：+1,On3S6)=Y84n7O,其它n4一，0h3X3O2)=Y3-,4h70,其它nX.七()=CoSAx5(?)=sin24(r)=cos8+cos16加+cos20万应当注意，如果给出的是连续信号儿(。，则首先要根据其最高频率确定采样速率工以及由频率分辨率选择采样点数N,然后对其进行软件采样(即计算K)=%(7j,0N-1),产生对应序列对信号乙，频率分辨率的选择要以能分辨开其中的三个频率对应的谱线为准则。对周期序列，最好截取周期的整数倍进行谱分析