diff options
| author | prashantsinalkar | 2018-02-03 11:01:52 +0530 |
|---|---|---|
| committer | prashantsinalkar | 2018-02-03 11:01:52 +0530 |
| commit | 7bc77cb1ed33745c720952c92b3b2747c5cbf2df (patch) | |
| tree | 449d555969bfd7befe906877abab098c6e63a0e8 /3889/CH5/EX5.5/Ex5_5.sce | |
| parent | d1e070fe2d77c8e7f6ba4b0c57b1b42e26349059 (diff) | |
| download | Scilab-TBC-Uploads-master.tar.gz Scilab-TBC-Uploads-master.tar.bz2 Scilab-TBC-Uploads-master.zip | |
Diffstat (limited to '3889/CH5/EX5.5/Ex5_5.sce')
| -rw-r--r-- | 3889/CH5/EX5.5/Ex5_5.sce | 26 |
1 files changed, 26 insertions, 0 deletions
diff --git a/3889/CH5/EX5.5/Ex5_5.sce b/3889/CH5/EX5.5/Ex5_5.sce new file mode 100644 index 000000000..e94a985f3 --- /dev/null +++ b/3889/CH5/EX5.5/Ex5_5.sce @@ -0,0 +1,26 @@ +//Example 5.5
+//page 321
+//Control Systems: Principles and Design
+//M Gopal, Second Edition, Tata McGraw-Hill
+//Chapter:Concepts of Stability: Routh Stability Array
+
+xdel(winsid())//close all graphics Windows
+clear;
+clc;
+//transfer function
+s=%s;
+P=s^6+2*s^5+8*s^4 + 12* s^3 +20*s^2 +16*s +16;
+//routh array
+disp('Routh Array')
+r=routh_t(P)
+disp(r)
+[r,num]=routh_t(P)
+disp('As there are no sign changes,there is no root on right half of s plane')
+disp('However as two rows are having same values, we decide stability on basis of auxiliary equation')
+//Extracting and finding roots of auxiliary polynomial
+a=flipdim(r(2,:),2)
+b=kron(a(1:$), [1 zeros(1,1)])
+P1=roots(poly(coeff(b),"x","coeff"))
+disp(P1,"Roots of auxiliary equation=")
+disp("As roots of auxiliary equation are purely imaginary, system is marginally stable")
+
|