blob: e7c9302a64626ca7535370eb8d8b57a96d2d70a5 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
|
//caption:obtain_state_transistion_matrix
//example 12_56
//page 586
s=%s;
syms t
A=[0 1;0 -3]
[r c]=size(A);//size of matrix A
p=s*eye(r,c)-A;//s*I-A where I is identity matrix
q=det(p)//determinant of sI-A
r=inv(p)//inverse of sI-A
//for calculating state transistion matrix
ip=[0 0;0 0]
i=1;
j=1;
for i=1:2
for j=1:2
if(i==2 & j==1)
else
ip(i,j)=ilaplace(r(i,j),s,t);
j=j+1;
end
end
i=i+1;
end
r(2,1)=0
disp(ip,"state transistion matrix,ip(t)=");
|
