数字信号处理实验94781new.docVIP

实验一 实验结果 2 调用filter解差分方程，由系统对u n 的响应判断稳定性 程序如下： A [1,-0.9];B [0.05,0.05]; %系统差分方程系数向量B和A x1n [1 1 1 1 1 1 1 1 zeros 1,50 ]; %产生信号x1 n R8 n x2n ones 1,128 ; %产生信号x2 n u n hn impz B,A,58 ; %求系统单位脉冲响应h n subplot 2,2,1 ;y h n ;tstem hn,y ; %调用函数tstem绘图 title a 系统单位脉冲响应h n ;box on y1n filter B,A,x1n ; %求系统对x1 n 的响应y1 n subplot 2,2,2 ;y y1 n ;tstem y1n,y ; title b 系统对R8 n 的响应y1 n ;box on y2n filter B,A,x2n ; %求系统对x2 n 的响应y2 n subplot 2,2,4 ;y y2 n ;tstem y2n,y ; title c 系统对u n 的响应y2 n ;box on 运行结果： （3）调用conv函数计算卷积 程序如下： x1n [1 1 1 1 1 1 1 1 ]; %产生信号x1 n R8 n h1n [ones 1,10 zeros 1,10 ]; h2n [1 2.5 2.5 1 zeros 1,10 ]; y21n conv h1n,x1n ; y22n conv h2n,x1n ; figure 2 subplot 2,2,1 ;y h1 n ;tstem h1n,y ; %调用函数tstem绘图 title d 系统单位脉冲响应h1 n ;box on subplot 2,2,2 ;y y21 n ;tstem y21n,y ; title e h1 n 与R8 n 的卷积y21 n ;box on subplot 2,2,3 ;y h2 n ;tstem h2n,y ; %调用函数tstem绘图 title f 系统单位脉冲响应h2 n ;box on subplot 2,2,4 ;y y22 n ;tstem y22n,y ; title g h2 n 与R8 n 的卷积y22 n ;box on 运行结果： （4）谐振器分析 程序如下： un ones 1,256 ; %产生信号u n n 0:255; xsin sin 0.014*n +sin 0.4*n ; %产生正弦信号 A [1,-1.8237,0.9801];B [1/100.49,0,-1/100.49]; %系统差分方程系数向量B和A y31n filter B,A,un ; %谐振器对u n 的响应y31 n y32n filter B,A,xsin ; %谐振器对u n 的响应y31 n figure 3 subplot 2,1,1 ;y y31 n ;tstem y31n,y ; title h 谐振器对u n 的响应y31 n ;box on subplot 2,1,2 ;y y32 n ;tstem y32n,y ; title i 谐振器对正弦信号的响应y32 n ;box on 运行结果： 思考题 （1）①对输入信号序列分段；②求单位脉冲响应h n 与各段的卷积；③将各段卷积结果相加。 （2）如果信号经过低通滤波器，把信号的高频分量滤掉，时域信号的剧烈变化将被平滑，由实验内容（1）结果图 a 、 b 和 c 可见，经过系统低通滤波使输入信号、和的阶跃变化变得缓慢上升与下降。 实验三 实验程序及实验结果 实验程序（1）： x1n [ones 1,4 ]; %产生序列向量x1 n R4 n M 8;xa 1: M/2 ; xb M/2 :-1:1; x2n [xa,xb]; %产生长度为8的三角波序列x2 n x3n [xb,xa]; X1k8 fft x1n,8 ; %计算x1n的8点DFT X1k16 fft x1n,16 ; %计算x1n的16点DFT X2k8 fft x2n,8 ; %计算x1n的8点DFT X2k16 fft x2n,16 ; %计算x1n的16点DFT X3k8 fft x3n,8 ; %计算x1n的8点DFT X3k16 fft x3n,16 ; %计算x1n的16点DFT %以下绘制幅频特性曲线 subplot 2,2,1 ;mstem X1k8 ; %绘制8点DFT的幅频特性图 title 1a 8点DFT[x_