initial commit / add all books

1 files changed, 32 insertions, 0 deletions

diff --git a/28/CH5/EX5.9/ex5_9.sce b/28/CH5/EX5.9/ex5_9.sce
new file mode 100755
index 000000000..e591f5a00
--- /dev/null
+++ b/28/CH5/EX5.9/ex5_9.sce
@@ -0,0 +1,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)=")