《数字信号的处理》课后上机的题目.pdfVIP

实用标准文案 实验报告 第一章 : 时域离散信号和时域离散系统 *16. 已知两个系统的差分方程分别为 (1) y(n)=0.6y(n-1)-0.08y(n-2)+x(n) (2) y(n)=0.7y(n-1)-0.1y(n-2)+2x(n)-x(n-2) 分别求出所描述的系统的单位脉冲响应和单位阶跃响应 . 解 :( 可附程序 ) (1) 系统差分方程的系数向量为 B1=1,A1=[1,-0.6,0.08] (2) 系统差分方程的系数向量为 B2=[2,0,-1],A2=[1,-0.7,0.1] 调 用 MATLAB 函 数 filter 计 算 两 个 系 统 的 单 位 脉 冲 响 应 和 单 位 阶 跃 响 应 的 程 序%B1=1;A1=[1,-0.6,0.08]; %设差分方程 (1) 系数向量 B2=[2,0,-1];A2=[1,-0.7,0.1]; %设差分方程 (2) 系数向量 %系统 1 xn=[1,zeros(1,30)]; ys=0; %xn=单位脉冲序列 , 长度 N=31 xi=filtic(B1,A1,ys); %由初始条件计算等效初始条件输入序列 xi hn1=filter(B1,A1,xn,xi); %调用 filter 解差分方程 , 求系统输出信号 hn1 n=0:length(hn1)-1; subplot(2,2,1);stem(n,hn1,.) title((a) 系统 1 的系统单位脉冲响应 ); xlabel(n);ylabel(h(n)) xn=ones(1,30); %xn=单位阶跃序列 , 长度 N=31 sn1=filter(B1,A1,xn,xi); %调用 filter 解差分方程 , 求系统输出信号 sn1 n=0:length(sn1)-1; subplot(2,2,2);stem(n,sn1,.) title((b) 系统 1 的单位阶跃响应 ); xlabel(n);ylabel(s(n)) %系统 2 xn=[1,zeros(1,30)]; %xn=单位脉冲序列 , 长度 N=31 xi=filtic(B2,A2,ys); %由初始条件计算等效初始条件输入序列 xi hn2=filter(B2,A2,xn,xi); 精彩文档 实用标准文案 %调用 filter 解差分方程 , 求系统输出信号 hn2 n=0:length(hn2)-1; subplot(2,2,3);stem(n,hn2,.) title((a) 系统 2 的系统单位脉冲响应 ); xlabel(n);ylabel(h(n)) xn=ones(1,30); %xn=单位阶跃序列 , 长度 N=31 sn2=filter(B2,A2,xn,xi); %调用 filter 解差分方程 , 求系统输出信号 sn2 n=0:length(sn2)-1; subplot(2,2,4);stem(n,sn2,.) title((b) 系统 2 的单位阶跃响应 ); xlabel(n);ylabel(s(n)) ( 结果展示 ): (a)系 统 1的 系 统 单 位 脉 冲