diff options
author | priyanka | 2015-06-24 15:03:17 +0530 |
---|---|---|
committer | priyanka | 2015-06-24 15:03:17 +0530 |
commit | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch) | |
tree | ab291cffc65280e58ac82470ba63fbcca7805165 /773/CH17/EX17.11 | |
download | Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2 Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip |
initial commit / add all books
Diffstat (limited to '773/CH17/EX17.11')
-rwxr-xr-x | 773/CH17/EX17.11/17_11.sci | 23 |
1 files changed, 23 insertions, 0 deletions
diff --git a/773/CH17/EX17.11/17_11.sci b/773/CH17/EX17.11/17_11.sci new file mode 100755 index 000000000..2a599a04d --- /dev/null +++ b/773/CH17/EX17.11/17_11.sci @@ -0,0 +1,23 @@ +//function//
+A=[0 1;-6 -5]
+x=[1;0];
+disp(x,"x(t)=')
+s=poly(0,'s');
+[Row Col]=size(A) //Size of a matrix
+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;
+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
+r=inv(q);
+r=simple(r); //To Find phi(-t)
+disp(r,"phi(-t)=")
+y=q*x; //x(t)=phi(t)*x(0)
+disp(y,"Solution To The given eq.=")
|