信号与线性系统实验报告(16)

仿真结果为：

2、差分方程求解

% 修正原design5_2.m(现在的design5_2_bak.m) % y(k)-y(k-1)+0.35y(k-2)=2x(k)-x(k-1) % 初始条件为yzi(-1)=2,yzi(-2)=1 % 激励为单位阶跃序列

% 注意Matlab中序列的下标从1开始，不能取0或负数；下标与时间变量n的意义不同

clear;

N=2; % 差分方程阶数

a=[1 -0.95 0.9025]; % 差分方程分母系数向量 a(0)~a(N) b=[1/3 1/3 1/3]; % 差分方程分子系数向量 b(0)~b(N) k=20; % 输出样点数目

zi=[2 1]; % 初始状态 N 个,初始条件为yzi(-1)=2,yzi(-2)=1 yzi=[0*ones(1,k+N+1)]; % 初始化零输入响应 h=yzi;% 初始化单位函数响应 yzs=yzi;% 初始化零状态响应 for n=1:N

yzi(n)=zi(N-n+1); end

y=yzi;% 初始化全响应，并令y(-2)=yzi(-2)=1,y(-1)=yzi(-1)=2 n=[-N:k];

x=[n==0];%与x=impseq(0,-N,k);作用相同 %x=impseq(0,-N,k);

zic=filtic(b,a,zi); % 将零输入初始条件yzi(-1)=2,yzi(-2)=1转化为函数filter()所需的初始条件

h(N+1:end)=filter(b,a,x(N+1:end)); % 或由命令dimpulse实现