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]
B=[0; 1]
phi=inv(s*eye(2,2)-A)
X=[1;1]
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;
disp(X," output = ")
|
