以此题为切入点,我们分析一下如何用matlap来解决这个问题。
我们知道求零状态响应,要利用激励卷积冲激响应得到零状态响应,所以要求零状态响应就先得掌握冲激响应和卷积。
而冲激响应如何得到呢?
借用他人的ppt来描述一下解题思路:
所以,我们把冲激响应转变成为求拉氏逆变换。
而拉氏逆变换有一个很重要的方法就是留数法。
在matlap中 有一个相应的语法为
[r,p,k] = residue(b,a)
[b,a] = residue(r,p,k)
因此,我们求得冲激响应之后,零状态响应=冲激响应*x(t)即可。
a=[1 4 3];
b=[1];
[r,p]=residue(b,a);
dt=0.05;
t=0:dt:10;
h=r(1)*exp(p(1)*t)+r(2)*exp(p(2)*t);
figure
plot(t,h)
N=length(t);
x1=ones(1,N);
y1=conv(x1,h)*dt;
figure
plot(t,y1(1:N))
x2=t;
y2=conv(x2,h)*dt;
figure
plot(t,y2(1:N))
打印结果:
(1)
(2)
(3)