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

实验一 用递推法解差分方程 以及用线性卷积法求线性时不变系统的输出 实验目的： 学习用递推法求解差分方程 学习用线性卷积法计算网络输出 实验原理 1．基于“递推法”的差分方程的时域求解 时域离散系统的输入输出关系可以用差分方程描述。 N阶常系数线性差分方程： （1） 需要预先给定条件：，才能进行差分方程的求解。 系统的因果性决定了递推法求解的顺序： 对于因果系统， 由当前时刻之前的N个状态，递推确定当前的系统输出； 对于非因果系统，由当前时刻之后的N个状态，递推确定当前的系统输出。 给定，利用差分方程，计算系统的输出 给定，利用差分方程，计算系统的单位脉冲响应 注：差分方程形式为 差分方程的时域求解方法中，较为典型的为“递推法”。 基于MATLAB的递推法差分方程的求解,往往采用filter函数。 其调用形式为：。 其中：y和x的起止位置一致，b [b0,…,bN],a [a0,…,aM] 2．线性卷积 对于线性时不变系统，若其单位脉冲响应为，则系统的输入输出关系还可以用线性卷积描述： （2） 其中： 若：已知有限长序列 则： 利用定义进行卷积的求取，关键就是在轴上寻找两个已知序列重叠的非零值区间。 实验内容 1．已知因果系统的差分方程如下： （1）输入信号，初始条件，试用递推法求解输出 （2）输入信号，初始条件，试用递推法求解输出 实现代码：xy100001.m（有关推导参见课本25页例1。5。2） % 利用差分方程，递推求解因果系统在不同初始状态的输出 % y n 0.9y n-1 +x n % n 0,…,150 % % Part 1,产生矩形序列，作为输入x n R10 n n [0:150] [x,nx] rectseq1 0,0,150,10 ; % Part 2,输入应为矩形序列，初始条件：y -1 1 y1 zeros size n ; y1 1 0.9+x 1 ; for k 1:150 y1 k+1 0.9*y1 k +x k+1 ; end % Part 3,输入应为矩形序列，初始条件：y -1 0 y2 zeros size n ; y2 1 x 1 ; for k 1:150 y2 k+1 0.9*y2 k +x k+1 ; end subplot 3,1,1 ; stem n,x ; axis [0,150,-0.1+min x ,0.1+max x ] ; title x n ; subplot 3,1,2 ; stem n,y1 ; axis [0,150,-0.1+min y1 ,0.1+max y1 ] ; title y1 n ; subplot 3,1,3 ; stem n,y2 ; axis [0,150,-0.1+min y2 ,0.1+max y2 ] ; title y2 n ; 运行结果： 2．已知系统差分方程为：, 用递推法求解因果系统的单位脉冲响应。要求：采用封闭形式表示；输出曲线。 解： 初始状态：； 所以： 综上： 基于MATLAB的递推法差分方程的求解,往往采用filter函数。 其调用形式为：。 其中：y和x的起止位置一致，b [b0,…,bN],a [a0,…,aM] ――――――――――――――――――――――――――――――――――――― 采用三个方法确定单位脉冲响应: xy100002.m % 利用差分方程求解因果系统的单位脉冲响应 % y n 0.9y n-1 +x n % n 0,…,150 % % Part 1,直接给出封闭形式的单位脉冲响应，记为h1（n） n [0:150] h1 0.9.^n; % Part 2,直接调用MATLAB函数，计算单位脉冲响应h2 n ,输入应为单位脉冲序列 a [1,-0.9]; b 1; x impseq 0,0,150 ; h2 filter b,a,x ; % Part 3,递推法，计算单位脉冲响应h3 n ,输入应为单位脉冲序列，初始条件：y -1 0 h3 zeros size n ; for k 0:150 if k 0 h3 k+1 1; else h3 k+1 0.9*h3 k ; end end subplot 3,1,1 ; stem n,h1 ; axis [0,150,-0.1+min h1 ,0.1+max h1 ] ; title 单位脉冲响应h1 n ; subplot 3,1,2 ; stem n,h2 ; axis [0,150,-0.1+min h2 ,0.1+max h2 ] ; title 单位脉冲响应h2 n ; subplot 3,1,3 ; stem n,h3 ; axis [0,150,-0.1+min h3 ,