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

2实验一信号、系统及系统响应一、实验目的1、熟悉连续信号经理想采样前后的频谱变化关系，加深对时域采样定理的理解；2、熟悉时域离散系统的时域特性；3、利用卷积方法观察分析系统的时域特性；4、掌握序列傅立叶变换的计算机实现方法，利用序列的傅立叶变换对连续信号、离散信号及系统响应进行频域分析。二、实验原理及方法采样是连续信号数字处理的第一个关键环节。对采样过程的研究不仅可以了解采样前后信号时域和频域特性发生变化以及信号信息不丢失的条件，而且可以加深对傅立叶变换、Z变换和序列傅立叶变换之间关系式的理解。对一个连续信号)(txa进行理想采样的过程可用下式表示：)()()(^tpttxxaa其中)(^txa为)(txa的理想采样，p(t)为周期脉冲，即mnTttp)()()(^txa的傅立叶变换为)]([1)(^smamjXTjaX上式表明^)(jXa为)(jXa的周期延拓。其延拓周期为采样角频率（T/2）。只有满足采样定理时，才不会发生频率混叠失真。在实验时可以用序列的傅立叶变换来计算^)(jXa。公式如下：TwjwaeXjX|)()(^离散信号和系统在时域均可用序列来表示。为了在实验中观察分析各种序列的频域特性，通常对)(jweX在[0，2]上进行M点采样来观察分析。对长度为N的有限长序列x(n)，有：njwNnjwkkemxeX)()(103其中，kMk2，k=0,1,……M-1时域离散线性非移变系统的输入/输出关系为mmnhmxnhnxny)()()(*)()(上述卷积运算也可在频域实现)()()(jjjeHeXeY三、实验程序s=yesinput(PleaseSelectTheStepOfExperiment:\n一.(1时域采样序列分析s=str2num(s);closeall;Xb=impseq(0,0,1);Ha=stepseq(1,1,10);Hb=impseq(0,0,3)+2.5*impseq(1,0,3)+2.2*impseq(2,0,3)+impseq(3,0,3);i=0;while(s);%时域采样序列分析if(s==1)l=1;k=0;while(1)if(k==0)A=yesinput('pleaseinputtheAmplitude:\n',...444.128,[100,1000]);a=yesinput('pleaseinputtheAttenuationCoefficient:\n',...222.144,[100,600]);w=yesinput('pleaseinputtheAngleFrequence(rad/s):\n',...222.144,[100,600]);endk=k+1;fs=yesinput('pleaseinputthesamplefrequence:\n',...1000,[100,1200]);Xa=FF(A,a,w,fs);i=i+1;string+['fs=',num2str(fs)];figure(i)DFT(Xa,50,string);1=yesinput1=str2num(1);end%系统和响应分析4elseif(s==2)kk=str2num(kk);while(kk)if(kk==1)m=conv(Xb,Hb);N=5;i=i+1;figure(i)string=('hb(n)');Hs=DFT(Hb,4,string);i=i+1;figure(i)string('xb(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)');elseif(kk==3)Xc=stepseq(1,1,5);m=conv(Xc,Ha);N=14;string=('y(n)=xc(n)*ha(n)');endendendi=i+1;figure(i)DFT(m,N,string);kk=yesinputkk=str2num(kk);end卷积定理的验证elseif(s==3)A=1;a=0.5;w=2,0734;fs=1;Xal=FF(A,a,w,fs);i=i+1;figure(i)string=('Thexal(n)(A=1,a=0.4,T=1)');[Xa,w]DFT(Xal,50,string);i=i+1;figure(i)5string=('hb(n)');Hs=DFT(Hb,4,string);Ys=Xs.*Hs;y=conv(Xal,Hb);N=53;i=i+1;figure(i)string=('y(n)=xa(n)*hb(n)');[yy,w]=DFT(y,N,string);i=i+1;figure(i)subplot(2,2,1)plot(w/pi,abs(yy));axis([-2202]);xlabel('w/pi');ylabel('|Ys(jw)|');title(FT[x(n)*h(n)]');subplot(2,2,3)plot(w/pi,abs(Ys));axis([-2202]);xlabel('w/pi');ylabel('|Ys(jw)|');title('FT[xs(n)].FT[h(n)]');endendend子函数:离散傅立叶变换及X(n),FT[x(n)]的绘图函数function[c,l]=DFT(x,N,str)n=0:N-1;k=-200:200;w=(pi/100)*k;l=w;c=x*Xc=stepseq(1,1,5);子函数：产生信号functionc=FF(A,a,w,fs)n=o:50-1;c=A*exp((-a)*n/fs).*sin(w*n/fs).*stepseq(0,0,49);子函数：产生脉冲信号function[x,n]=impseq(n0,n1,n2)n=[n1:n2];6x=[(n-n0)==0];子函数：产生矩形框信号function[x,n]=stepseq(n0,n1,n2)n=[n1:n2];x=[(n-n0=0)];四、实验内容及步骤1、认真复习采样理论，离散信号与系统，线性卷积，序列的傅立叶变换及性质等有关内容，阅读本实验原理与方法。2、编制实验用主程序及相应子程序。3、调通并运行实验程序，完成下列实验内容：①分析采样序列的特性，产生采样信号序列)(nxa，使A=444.128，a=502,2500。（)(txa的无失真采样频率约为1000Hz）。a.取采样频率sf=1kHz，即T=1ms。观察所得采样)(nxa的幅频特性|)(jeX|和原图中的幅频特性曲线在折叠频率附近有无明显差别。应当注意，实验中所得频谱是用序列)(nxa的傅立叶变换公式求得的，所以在频率度量上存在关系：T，为数字频率，为模拟频率。b.改变采样频率，sf=300Hz，观察|)(jeX|的变化，并做记录（打印曲线）；进一步降低采样频率，sf=200Hz，观察频谱混叠是否明显存在，说明原因，并记录（打印）这时的|)(jeX|曲线。a步实验结果如下图所示：由图形可知，当采样频率为1000Hz时，采样序列在折叠频率附近处，即=处无明显频谱混叠。b步实验结果如下图所示：7由图可知，当采样频率进一步降低时，主瓣宽度逐渐变宽，频率混叠现象也逐渐严重。存在较明显的失真现象。②时域离散信号、系统和系统响应分析。a.观察信号)(nxb和系统)(nhb的时域和频域特性；利用线性卷积求信号)(nxb通过系统)(nhb的响应y(n)，比较所求响应y(n)和)(nhb的时域及频域特性，注意它们之间有无差别。绘图说明，并用所学理论解释所得结果。实验结果如下图所示：8b.观察系统)(nha对信号)(nxc的响应特性。利用线性卷积求系统响应y(n)，并判断y(n)图形及其非零值序列长度是否与理论结果一致，对)()()(10nRnhnxac，说出一种定性判断y(n)图形正确与否的方法。调用序列傅立叶变换数值计算子程序，求得)(kjeY，观察|)(kjeY|特性曲线，定性判断结果的正确性。改变)(nxc的长度，取N=5，重复该实验。注意参数变化的影响，说明变化前后的差异，并解释所得结果。9实验结果如下图所示：N=10：N=5：欲判断结果正确与否，可以先对其进行运算，算出其卷积，再与图形对照。当N=10时，峰值较高，且峰值很窄，变换之后图形频带主值部分比较集中；N=5时情况与之相反。③卷积定理的验证。将实验②中的信号换为)(nxa，使a=0.4，,0734,20A=1，T=1，重复实验②a，打印|)(kjeY|曲线；对主程序做简单修改，按式（10.3.9）计算)()()(kkkjbjajeHeXeY，并绘出|)(kjeY|曲线，与前面直接对y(n)进行傅立叶变换所得幅频特性曲线进行比较，验证时域卷积定理。10实验所得结果如下：11四、思考题1、在分析理想采样序列特性的实验中，采样频率不同，相应理想采样序列的傅立叶变换频谱的数字频率度量是否都相同？它们所对应的模拟频率是否相同？为什么？答：由T可知，若采样频率不同，则其周期T不同，相应的数字频率也不相同；而因为是同一信号，故其模拟频率保持不变。2、在卷积定理验证的实验中，如果选用不同的频域采样点数M值，例如，选M=10和M=20，分别做序列的傅立叶变换，求得)()()(kkkjbjajeHeXeY，k=0,1,…,M-1所得结果之间有无差异？为什么？答：有差异。因为所得)(kjeY图形由其采样点数唯一确定，由频域采样定理可知，若M小于采样序列的长度N，则恢复原序列时会发生时域混叠现象。12实验二用FFT作谱分析一、实验目的1、进一步加深DFT算法原理和基本性质的理解(因为FFT只是DFT的一种快速算法，所以FFT的运算结果必然满足DFT的基本性质)。2、熟悉FFT算法原理和FFT子程序的应用。3、学习用FFT对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析误差及其原因，以便在实际中正确应用FFT。二、实验步骤1、复习DFT的定义、性质和用DFT作谱分析的有关内容。2、复习FFT算法原理与编程思想，并对照DIT—FFT运算流图和程序框图，读懂本实验提供的FFT子程序。3、编制信号产生子程序，产生以下典型信号供谱分析：nRnx411n，30nnx2n8，74n0，其它nn4，30nnx3n3，74n0，其它nnnx4cos4nnx8sin5ttttx20cos16cos8cos6应当注意，如果给出的是连续信号txa，则首先要根据其最高频率确定采样速率sf以及由频率分辨率选择采样点数N，然后对其进行软件采样(即计算nTxnxa，10Nn)，产生对应序列nx。对信号tx6，频率分辨率的选择要以能分辨开其中的三个频率对应的谱线为准则。对周期序列，最好截取周期的整数倍进行谱分析，否则有可能产生较大的分析误差。134、编写主程序下图给出了主程序框图，供参考。本实验提供FFT子程序和通用绘图子程序。主程序框图开始读入长度N调用信号产生子程序产生实验信号调用绘图子程序(函数)绘制时间序列波形图调用FFT子程序(函数)计算信号的DFT调用绘图子程序(函数)绘制kX曲线结束14三、实验结果直接运行程序，按照实验内容及程序提示键入1~8，分别对nx1~nx6及nxnxnx547、njxnxnx548进行谱分析。输出nxnx51~的波形及其8点DFT和16点DFT，nx6的16点、32点和64点采样序列及其DFT。1、nx1及其8点和16点DFT2、nx2及其8点和16点DFT153、nx