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 /46/CH6/EX6.1/Example6_1.sce | |
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 '46/CH6/EX6.1/Example6_1.sce')
-rwxr-xr-x | 46/CH6/EX6.1/Example6_1.sce | 26 |
1 files changed, 26 insertions, 0 deletions
diff --git a/46/CH6/EX6.1/Example6_1.sce b/46/CH6/EX6.1/Example6_1.sce new file mode 100755 index 000000000..cc04a958a --- /dev/null +++ b/46/CH6/EX6.1/Example6_1.sce @@ -0,0 +1,26 @@ +//Example 6.1
+clc;
+syms s t;
+tau=1;//min
+R=1/9;//ft/cfm
+A=9;
+//from Equation 6.8
+g=R/(tau*s+1);
+disp(g,'H(s)/Q(s)=')
+//from Example 4.5
+disp('Q(t)=90[u(t)-u(t-0.1)')
+//where u(t) is a unit step function, the laplace transform of it gives
+Qs=90*(1-exp(-0.1*s))/s
+disp(Qs,'Q(s)=')
+Hs=Qs*g;
+disp(Hs,'H(s)=')
+//taking first term for t<0.1, the second term goes equals to zero
+Ht=ilaplace('10*(1/(s*(s+1)))',s,t);//t<0.1
+disp(Ht,'H(t)=')
+disp('H(t)=10(1-yexp(-(-t-0.1)))')//t>0.1
+Ht=10*((1-exp(-t))-(1-exp(-(-t-0.1))));
+disp(Ht,'H(t)=')
+//from Eq.(5.16)
+Ht=R*A*exp(-(t/tau));//impulse
+disp(Ht,'H(t)=')
+
|