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/CH24/EX24.1/Example24_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/CH24/EX24.1/Example24_1.sce')
-rwxr-xr-x | 46/CH24/EX24.1/Example24_1.sce | 19 |
1 files changed, 19 insertions, 0 deletions
diff --git a/46/CH24/EX24.1/Example24_1.sce b/46/CH24/EX24.1/Example24_1.sce new file mode 100755 index 000000000..ef339d7e0 --- /dev/null +++ b/46/CH24/EX24.1/Example24_1.sce @@ -0,0 +1,19 @@ +//Example 24.1
+clc
+syms K b z w;
+Gz=K*(1-b)/(z-b)
+//where b=exp(-T/tau)
+//From Eq.(24.4)
+z=w+1/w-1;
+Gz=eval(Gz)
+disp('=0',1+Gz,'1+G(z)=')
+//which is equivalent to
+disp('(K+1)*(1-b)*w+(1+b)-K(1-b)=0')
+routh=[(K+1)*(1-b);(1+b)-K*(1-b)]
+//b is always positive and less than one and K is positive
+//The first element in the array is positive
+//For stability, the Routh test requires that all elements of the first column be positive
+//Therefore,
+disp('>0',routh(2,1))
+disp('K<(1+b)/(1-b)')
+
|