《数字信号处理实验信号、系统及系统响应实验报告》由会员分享，可在线阅读，更多相关《数字信号处理实验信号、系统及系统响应实验报告（15页珍藏版）》请在金锄头文库上搜索。
1、 实验一 信号、系统及系统响应 1.实验目的（1） 熟悉连续信号经理想采样前后的频谱变化关系，加深对时域采样定理的理解。（2） 熟悉时域离散系统的时域特性。（3） 利用卷积方法观察分析系统的时域特性。（4） 掌握序列傅里叶变换的计算机实现方法，利用序列的傅里叶变换对连续信号、离散信号及系统响应进行频域分析。2.实验原理采样是连续信号数字处理的第一个关键环节。对采样过程的研究不仅可以了解采样前后信号时域和频域特性发生变化以及信号信息不丢失的条件，而且可以加深对傅立叶变换、Z变换和序列傅立叶变换之间关系式的理解。对一个连续信号进行理想采样的过程可用下式表示： 其中为的理想采样，p(t)为周期脉冲，
2、即 的傅立叶变换为 上式表明为的周期延拓。其延拓周期为采样角频率（）。只有满足采样定理时，才不会发生频率混叠失真。在实验时可以用序列的傅立叶变换来计算。公式如下：离散信号和系统在时域均可用序列来表示。为了在实验中观察分析各种序列的频域特性，通常对在0，2上进行M点采样来观察分析。对长度为N的有限长序列x(n)，有： 其中，k=0,1,M-1时域离散线性非移变系统的输入/输出关系为 上述卷积运算也可在频域实现 3.实验环境 应用MATLAB 6.5软件 操作系统：windows XP 开始4.主程序流程图 选择步骤1 产生信号序列xa(n)Xa(jw),并绘制xa(n)|Xa(jw)| 选择步骤
3、2,然后1 产生信号序列hb(n),并绘制hb(n) 调用序列傅氏变换数值计算子程序，求Hb(jw),并绘制|Hb(jw)| 求yb(n)=xb(n)*hb(n)，并绘制yb(n) 调用序列傅氏变换数值计算子程序，求Yb(jw),并绘制|Yb(jw)| 选择步骤2 求ya(n)=xc(n)*ha(n)，并绘制ya(n) 调用序列傅氏变换数值计算子程序，求Ya(jw),并绘制|Ya(jw)| 选择步骤3 求ya(n)=xc(n)*ha(n)，并绘制ya(n) 调用序列傅氏变换数值计算子程序，求Ya(jw),并绘制|Ya(jw)| 选择步骤3 产生信号序列xa(n),并绘制xa(n) 调用序列傅氏
4、变换数值计算子程序，求Xa(jw),并绘制|Xa(jw)| 求ya(n)=xc(n)*ha(n)，并绘制ya(n) 调用序列傅氏变换数值计算子程序，求Ya(jw),并绘制|Ya(jw)| 求|Y(jw)|=|Xa(jw)Hb(jw)|，并绘制|Y(jw)| 选择步骤 结束5.实验结果（1） 采样序列的特性。 一般称fs/2为折叠频率，只有当信号最高频率不超过该频率时才不会发生混叠现象，否则超过了fs/2的频率会折叠回来形成混叠现象，因此频率混叠均产生在fs/2附近。A. 采样频率fs=1000Hz 由图形可知，当采样频率为1000Hz时，采样序列在折叠频率附近处，即w=处无明显频谱混叠。B.采
5、样频率 fs=300Hz C.采样频率fs=200Hz 由图可知，当采样频率进一步降低时，主瓣宽度逐渐变宽，频率混叠现象也逐渐严重，存在较明显的失真现象。原因是采样频率太小，使最高频率fc超过了fs/2,超过了fs/2的频率会折叠回来而形成的混叠现象。（2）时域离散信号、系统和系统响应。 A. 理论值一个函数与单位脉冲序列的卷积等于函数本身，卷积得到的长度等于两个函数长度和减一。由图可知，yb(n)=xb(n），其长度13=4+10-1，所以理论与实际是一致的。B. 判断ya(n)是否正确的方法：ya(n)的长度L等于两个被卷积函数的长度和减去一，且ya(n)是关于n=(L-1)/2对称的，峰
6、值即为N值，对称轴左边由一逐渐按增一序列递增，右边按减一序列递减。由图知：19=10+10-1，且图形正确，所以做出的ya(n)是正确的。C. 当N=10时，峰值较高，且峰值很窄，变换之后图形频带主值部分比较集中，且峰值较高；当N=5时，峰值较矮，且峰值很宽，变换之后图形频带主值部分较为分散，且峰值较矮。（3）卷积定理的验证 a=0.4,=2.0734,A=1,T=1 Y(jw)=Xa(jw)*Hb(jw) 由图可知，由yb(n)=xa(n)*hb(n)经傅氏变换所得到的|Yb(jw)|和由|Yb(jw)|=|Xa(jw)Hb(jw)|所得到的|Yb(jw)|的图像是一样的，从而验证了时域卷积
7、定理。6.实验代码s = -1;while(s0) clc;s = input(*信号、系统及响应*nn选择实验步骤(默认1):n1:时域采样序列分析n2:系统和响应分析n3:卷积定理验证n0:退出n选择:,s);switch(s) case1,2,3,0 s = str2num(s); case s = 1; otherwise s = -1; endendclose all; while(s) %时域采样序列分析 if (s=1) A=444.128; a=50*sqrt(2)*pi; w=50*sqrt(2)*pi; n=0:50-1; fs = input(输入采样频率nfs=,s);
8、 %fs=1000,300,200 fs = str2num(fs); if isempty(fs) fs = 1000; disp(输入数据格式错误,使用默认值1000); else if(fs1) fs = 1000; disp(输入无效数据,使用默认值1000); end end c=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
9、,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):n1:容an2:容bn3:容b中xc(n)的长度改为5n0:退出n选择:,s); switch(l) case 1,2,3,0 l = str2num(l); otherwise l = 1; end while(l) if(l=1) %hb(n)的时域序列 hb=1,2.5,2.5,1; i=0:3; subplot(2,
10、2,1); stem(i,hb,.); axis(0 3 0 2.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,.
