数字信号处理实验.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:150y1(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:150y2(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:150if k==0h3(k+1)=1;elseh3(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);