1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
|
A=[0 1; -2 -3]
U=[1/s;1/(s+2)]
B=[2 1; 0 1]
phi=inv(s*eye(2,2)-A)
X=[0;0]
a1=ilaplace(phi(1,1),s,t)
a2=ilaplace(phi(1,2),s,t)
a3=ilaplace(phi(2,1),s,t)
a4=ilaplace(phi(2,2),s,t)
S=[a1 a2;a3 a4]
disp(S,"%e^(A*t) = ")
ZIR=S*X
disp(ZIR," ZIR = ")
k=phi*B*U
b1=ilaplace(k(1,1),s,t)
b3=ilaplace(k(2,1),s,t)
ZSR=[b1;b3]
disp(ZSR,"ZSR = ")
X=ZIR+ZSR;
|
