blob: 30645d251c77016ae7210502af47553037e638b0 (
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
27
28
29
30
31
32
33
34
|
//function//
A=[0 1;-3 -4]
B=[0;1]
C=[1 0]
x=[0;0]
disp(x,"x(t)=')
s=poly(0,'s');
[Row Col]=size(A) //Size of a matrix A
m=s*eye(Row,Col)-A //sI-A
n=det(m) //To Find The Determinant of si-A
p=inv(m) ; // To Find The Inverse Of sI-A
syms t s m;
disp(p,"phi(s)=") //Resolvent Matrix
for i=1:Row
for j=1:Col
//Taking Inverse Laplace of each element of Matrix(phi(s))
q(i,j)=ilaplace(p(i,j),s,t);
end;
end;
disp(q,"phi(t)=") //State Transition Matrix
t=(t-m)
q=eval(q) //At t=t-m ,evaluating q i.e phi(t-m)
r=integ(q*B,m)//Integrate q w.r.t m (Indefinite Integration)
m=0 //Upper limit is t
g=eval(r) //Putting the value of upper limit in q
m=t //Lower Limit is 0
h=eval(r) //Putting the value of lower limit in q
y=(h-g);
printf("x(t)= phi(t)*x(0) + integ(phi(t-m)*B) w.r.t m from 0 t0 t \n")
//x(t)=phi(t)*x(0) + integ(phi(t-m)*B) w.r.t m from 0 t0 t
y1=(q*x)+y;
disp(y1,"x(t)=")
y2=C*y1;
disp(y2,"Output Response=")
|
