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 /29/CH7/EX7.5.8/exa7_5_8.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 '29/CH7/EX7.5.8/exa7_5_8.sce')
-rwxr-xr-x | 29/CH7/EX7.5.8/exa7_5_8.sce | 27 |
1 files changed, 27 insertions, 0 deletions
diff --git a/29/CH7/EX7.5.8/exa7_5_8.sce b/29/CH7/EX7.5.8/exa7_5_8.sce new file mode 100755 index 000000000..b3b883b12 --- /dev/null +++ b/29/CH7/EX7.5.8/exa7_5_8.sce @@ -0,0 +1,27 @@ +//caption:value_of_K_in_terms_of_T1_and_T2
+//example 7.5.8
+//page 207
+s=%s;
+syms K T1 T2;
+m=s*(s*T1+1)*(s*T2+1)
+G=K/m;
+CH=1+G;
+disp("on simplyfying CH")
+CH=m+K;
+CH=simple(CH);
+disp('=0',CH,"characterstics_eq,CH=");
+c0=coeffs(CH,'s',0);
+c1=coeffs(CH,'s',1);
+c2=coeffs(CH,'s',2);
+c3=coeffs(CH,'s',3);
+b=[c0 c1 c2 c3]
+n=4;
+routh=[b([4,2]);b([3,1])];
+routh=[routh;-det(routh)/routh(2,1),0]
+t=routh(2:3,1:2)
+routh=[routh;-det(t)/t(2,1),0]
+disp(routh,"routh=")
+disp("for given system to be stable:");
+disp("K>0 and (-(K*T1*T2-T2-T1)/(T2+T1))>0");
+disp("which gives:");
+disp("0<K<((1/T1)+(1/T2))");
\ No newline at end of file |