matlab求解系统响应filter,信号与系统MATLAB实验（教师版）

例2-6已知 y(k)?y(k?1)?0.9y(k?2))?f(k) 求单位响应。

a=[1,-1,0.9]; b=[1]; impz(b,a) impz(b,a,60) impz(b,a,-10:40)

4． 离散系统的零状态响应

MATLAB提供求离散系统零状态响应数值解函数filter，调用格式为filter(b,a,x),式中b和a是表示离散系统的向量，x是输入序列非零样值点行向量，输出向量序号同x一样。 例2-7 已知 y(k)?0.25y(k?1)?0.5y(k?2))?f(k)?f(k?1)， f(k)?()?(k) 求零状态响应， 范围0~20。 a=[1 -0.25 0.5]; b=[1 1]; t=0:20; x=(1/2).^t; y=filter(b,a,x) subplot(2,1,1) stem(t,x)

title('输入序列') subplot(2,1,2) stem(t,y)

title('响应序列') 三、上机实验内容 1．验证实验原理中程序

2．已知2y(k)?2y(k?1)?y(k?2)?f(k)?3f(k?1)?2f(k?2)，画单位响应波形。 3．已知y(k)?y(k?1)?0.25y(k?2)?f(k)，输入f(t)??(k)，画输出波形，范围0~15。

12k

13

实验三 连续时间LTI系统的时域分析

一、实验目的

1．学会用MATLAB求解连续系统的零状态响应； 2. 学会用MATLAB求解冲激响应及阶跃响应； 3．学会用MATLAB实现连续信号卷积的方法； 二、实验原理

1．连续时间系