summaryrefslogtreecommitdiff
path: root/62/CH7/EX7.48/ex_7_48.sce
diff options
context:
space:
mode:
authorpriyanka2015-06-24 15:03:17 +0530
committerpriyanka2015-06-24 15:03:17 +0530
commitb1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch)
treeab291cffc65280e58ac82470ba63fbcca7805165 /62/CH7/EX7.48/ex_7_48.sce
downloadScilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz
Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2
Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip
initial commit / add all books
Diffstat (limited to '62/CH7/EX7.48/ex_7_48.sce')
-rwxr-xr-x62/CH7/EX7.48/ex_7_48.sce22
1 files changed, 22 insertions, 0 deletions
diff --git a/62/CH7/EX7.48/ex_7_48.sce b/62/CH7/EX7.48/ex_7_48.sce
new file mode 100755
index 000000000..b2a7273f6
--- /dev/null
+++ b/62/CH7/EX7.48/ex_7_48.sce
@@ -0,0 +1,22 @@
+clc;
+A=[0 1;-6 -5];
+b=[0;1];
+c=[1 0];
+q0=[2;1];
+syms t T
+//y=c*e^(A*t)*q0+integ(c*e^A*(t-T)*b*x,0,t)
+e=spec(A);//eigen values
+b0=e(1)*%e^(e(2)*t)-e(2)*%e^(e(1)*t);
+b1=%e^(e(1)*t)-%e^(e(2)*t);
+eAt=b0*eye(A)+b1*A;
+b0=e(1)*%e^(e(2)*(t-T))-e(2)*%e^(e(1)*(t-T));
+b1=%e^(e(1)*(t-T))-%e^(e(2)*(t-T));
+eAtT=b0*eye(A)+b1*A;
+disp(eAt,"e^(At)");
+ceAtq0=c*eAt*q0;
+ceAtTq0=c*eAtT*q0;
+y=ceAtq0+integ(ceAtTq0*%e^-T,T);
+y0=limit(y,T,0);
+yt=limit(y,T,t);
+Y=yt-y0;
+disp(Y,"y(t)=") \ No newline at end of file