数字信号处理实验-信号、系统及系统响应实验报告

实验一信号、系统及系统响应1.实验目的（1）熟悉连续信号经理想采样前后的频谱变化关系，加深对时域采样定理的理解。（2）熟悉时域离散系统的时域特性。（3）利用卷积方法观察分析系统的时域特性。（4）掌握序列傅里叶变换的计算机实现方法，利用序列的傅里叶变换对连续信号、离散信号及系统响应进行频域分析。2.实验原理采样是连续信号数字处理的第一个关键环节。对采样过程的研究不仅可以了解采样前后信号时域和频域特性发生变化以及信号信息不丢失的条件，而且可以加深对傅立叶变换、Z变换和序列傅立叶变换之间关系式的理解。对一个连续信号)(txa进行理想采样的过程可用下式表示：)()()(^tpttxxaa其中)(^txa为)(txa的理想采样，p(t)为周期脉冲，即mnTttp)()()(^txa的傅立叶变换为)]([1)(^smamjXTjaX上式表明^)(jXa为)(jXa的周期延拓。其延拓周期为采样角频率（T/2）。只有满足采样定理时，才不会发生频率混叠失真。在实验时可以用序列的傅立叶变换来计算^)(jXa。公式如下：^()()jnTaanXjxnTe()()jjnnXexne^()()|jaTXjXe《数字信号处理》实验——信号、系统及系统响应-1-离散信号和系统在时域均可用序列来表示。为了在实验中观察分析各种序列的频域特性，通常对)(jweX在[0，2]上进行M点采样来观察分析。对长度为N的有限长序列x(n)，有：njwNnjwkkemxeX)()(10其中，kMk2，k=0,1,……M-1时域离散线性非移变系统的输入/输出关系为mmnhmxnhnxny)()()(*)()(上述卷积运算也可在频域实现)()()(jjjeHeXeY3.实验环境应用MATLAB6.5软件操作系统：windowsXP4.主程序流程图开始《数字信号处理》实验——信号、系统及系统响应-2-5.实验结果产生信号序列xa(n),并绘制xa(n)选择步骤2,,然后1产生信号序列hb(n),并绘制hb(n)调用序列傅氏变换数值计算子程序，求Hb(jw),并绘制|Hb(jw)|求yb(n)=xb(n)*hb(n)，并绘制yb(n)调用序列傅氏变换数值计算子程序，求Yb(jw),并绘制|Yb(jw)|求ya(n)=xc(n)*ha(n)，并绘制ya(n)调用序列傅氏变换数值计算子程序，求Ya(jw),并绘制|Ya(jw)|选择步骤2选择步骤3求ya(n)=xc(n)*ha(n)，并绘制ya(n)调用序列傅氏变换数值计算子程序，求Ya(jw),并绘制|Ya(jw)|选择步骤3产生信号序列xa(n)Xa(jw),并绘制xa(n)|Xa(jw)|调用序列傅氏变换数值计算子程序，求Xa(jw),并绘制|Xa(jw)|求ya(n)=xc(n)*ha(n)，并绘制ya(n)调用序列傅氏变换数值计算子程序，求Ya(jw),并绘制|Ya(jw)|求|Y(jw)|=|Xa(jw)Hb(jw)|，并绘制|Y(jw)|选择步骤结束选择步骤1《数字信号处理》实验——信号、系统及系统响应-3-（1）采样序列的特性。一般称fs/2为折叠频率，只有当信号最高频率不超过该频率时才不会发生混叠现象，否则超过了fs/2的频率会折叠回来形成混叠现象，因此频率混叠均产生在fs/2附近。A.采样频率fs=1000Hz由图形可知，当采样频率为1000Hz时，采样序列在折叠频率附近处，即w=处无明显频谱混叠。B.采样频率fs=300HzC.采样频率fs=200Hz《数字信号处理》实验——信号、系统及系统响应-4-由图可知，当采样频率进一步降低时，主瓣宽度逐渐变宽，频率混叠现象也逐渐严重，存在较明显的失真现象。原因是采样频率太小，使最高频率fc超过了fs/2,超过了fs/2的频率会折叠回来而形成的混叠现象。（2）时域离散信号、系统和系统响应。A.)()(nnxb)3()2(5.2)1(5.2)()(nnnnnhb《数字信号处理》实验——信号、系统及系统响应-5-理论值一个函数与单位脉冲序列的卷积等于函数本身，卷积得到的长度等于两个函数长度和减一。由图可知，yb(n)=xb(n），其长度13=4+10-1，所以理论与实际是一致的。B.)()()(10nRnhnxac判断ya(n)是否正确的方法：ya(n)的长度L等于两个被卷积函数的长度和减去一，且ya(n)是关于n=(L-1)/2对称的，峰值即为N值，对称轴左边由一逐渐按增一序列递增，右边按减一序列递减。由图知：19=10+10-1，且图形正确，所以做出的ya(n)是正确的。《数字信号处理》实验——信号、系统及系统响应-6-C.)()(5nRnxc当N=10时，峰值较高，且峰值很窄，变换之后图形频带主值部分比较集中，且峰值较高；当N=5时，峰值较矮，且峰值很宽，变换之后图形频带主值部分较为分散，且峰值较矮。（3）卷积定理的验证a=0.4,Ω=2.0734,A=1,T=1《数字信号处理》实验——信号、系统及系统响应-7-Y(jw)=Xa(jw)*Hb(jw)由图可知，由yb(n)=xa(n)*hb(n)经傅氏变换所得到的|Yb(jw)|和由|Yb(jw)|=|Xa(jw)Hb(jw)|所得到的|Yb(jw)|的图像是一样的，从而验证了时域卷积定理。6.实验代码s=-1;while(s0)clc;s=input('******信号、系统及响应******\n\n选择实验步骤(默认1):\n[1]:时域采样序列分析\n[2]:系统和响应分析\n[3]:卷积定理验证\n[0]:退出\n选择:','s');switch(s)case{'1','2','3','0'}s=str2num(s);case{''}s=1;otherwises=-1;endendcloseall;while(s)%时域采样序列分析if(s==1)A=444.128;a=50*sqrt(2)*pi;w=50*sqrt(2)*pi;n=0:50-1;《数字信号处理》实验——信号、系统及系统响应-8-fs=input('输入采样频率\nfs=','s');%fs=1000,300,200fs=str2num(fs);ifisempty(fs)fs=1000;disp('输入数据格式错误,使用默认值1000');elseif(fs1)fs=1000;disp('输入无效数据,使用默认值1000');endendc=A*exp((-a)*n/fs).*sin(w*n/fs);subplot(2,2,1);stem(n,c,'.');xlabel('n');ylabel('xa(n)');title('xa(n)的时域序列');N=50;k=-200:200;w=k*pi/100;X=DFT(c,N);subplot(2,2,2);plot(w/pi,abs(X));xlabel('w/pi');ylabel('|X(jw)|');title('xa(n)的傅氏变换|X(jw)|');else%系统和响应分析if(s==2)l=input('系统和响应分析,请选择时域信号类型(默认1):\n[1]:内容②a\n[2]:内容②b\n[3]:内容②b中xc(n)的长度改为5\n[0]:退出\n选择:','s');switch(l)case{'1','2','3','0'}l=str2num(l);otherwisel=1;endwhile(l)if(l==1)%hb(n)的时域序列hb=[1,2.5,2.5,1];i=0:3;《数字信号处理》实验——信号、系统及系统响应-9-subplot(2,2,1);stem(i,hb,'.');axis([0302.5]);xlabel('n');ylabel('hb(n)');title('hb(n)的时域序列');%hb(n)的傅氏变换|Hb(jw)|N=4;k=-200:200;w=k*pi/100;Hb=DFT(hb,N);subplot(2,2,2);plot(w/pi,abs(Hb));xlabel('w/pi');ylabel('|Hb(jw)|');title('hb(n)的傅氏变换|Hb(jw)|');%xb(n)的时域序列xb=[1,0,0,0,0,0,0,0,0,0];i=0:9;subplot(2,2,3);stem(i,xb,'.');xlabel('n');ylabel('xb(n)');title('xb(n)的时域序列');%xb(n)的傅氏变换(Xb|jw|)N=10;k=-200:200;w=k*pi/100;Xb=DFT(xb,N);magXb=abs(Xb);subplot(2,2,4);plot(w/pi,magXb);xlabel('w/pi');ylabel('|Xb(jw)|');title('xb(n)的傅氏变换(Xb|jw|)');%yb(n)=xb(n)*hb(n)的时域序列yb=conv(xb,hb);figure;subplot(2,1,1);stem(0:12,yb,'.');xlabel('n');ylabel('yb(n)=xb(n)*hb(n)');title('yb(n)=xb(n)*hb(n)的时域序列');%yb(n)的傅氏变换(Yb|jw|)《数字信号处理》实验——信号、系统及系统响应-10-N=13;k=-200:200;w=k*pi/100;Yb=DFT(yb,N);subplot(2,1,2);plot(w/pi,abs(Yb));xlabel('w/pi');ylabel('|Yb(jw)|');title('yb(n)的傅氏变换|Yb(jw)|');elseif(l==2)%ya(n)=xc(n)*ha(n)的时域序列ha=[1,1,1,1,1,1,1,1,1,1];xc=ha;ya=conv(ha,xc);subplot(2,1,1);stem(0:18,ya,'.');xlabel('n');ylabel('ya(n)=xc(n)*ha(n)');title('ya(n)=xc(n)*ha(n)的时域序列');%ya(n)的傅氏变换(Ya|jw|)N=19;k=-200:200;w=k*pi/100;Ya=DFT(ya,N);subplot(2,1,2);plot(w/pi,abs(Ya));xlabel('w/pi');ylabel('|Ya(jw)|');title('ya(n)的傅氏变换|Ya(jw)|');elseif(l==3)ha=[1,1,1,1,1,1,1,1,1,1];xc=[1,1,1,1,1];%ya(n)=xc(n)*ha(n)的时域序列ya=conv(ha,xc);subplot(2,2,1);stem(0:13,ya,'.');xlabel('n');ylabel('ya(n)=xc(n)*ha(n)');title('ya(n)=xc(n)*ha(n)的时域序列');%ya(n)的傅氏变换(Ya|jw|)N=14;k=-200:200;《数字信号处理》实验——信号、系统及系统响应-11-w=k*pi/100;Ya=DFT(ya,N);subplot(2,2,2);plot(w/pi,abs(Ya));xlabel('w/pi');ylabel('|Ya(jw)|');title('ya(n)的傅氏变换(|Ya(jw)|');endendendl=input('请再选择信号类型(默认1):\n[1]:内容②a\n[2]:内容②b\n[3]:内容②b中xc(n)的长度改为5\n[0]:退出\n选择:','s');switch(l)case{'1','2','3','0'}l=str2num(l);otherwisel=1;endend%卷积定理验证elseif(s==3)A=1;a=0