blob: e591f5a00257d58fef7722c7ab1cac7ffc51e59e (
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
|
s=%s;
syms t m
A=[0 1;-100 -20];
B=[0;100];
C=[1 0];
x=[0;0];
[r c]=size(A)
p=s*eye(r,c)-A
q=inv(p);
disp(q,"phi(s)=") // Resolvant matrix
for i=1:r;
for j=1:c;
q(i,j)=ilaplace(q(i,j),s,t)
end
end
disp(q,"phi(t)=") // State transition matrix
t=t-m;
q=eval(q)
// Integrate q w.r.t m
r=integrate(q*B,m)
m=0 // Upper limit is t
g=eval(r) // Puting upper limit in q
m=t // Lower limit is 0
h=eval(r) // Putting lower limit in q
y=(h-g);
disp(y,"y=")
printf("x(t)=phi(t)*x(0)+integrate(phi(t-m*B) w.r.t m from 0 to t)")
y1=(q*x)+y;
disp(y1,"x(t)=")
// transfer function
t=C*q*B;
disp(t,"T(s)=")
|
