summaryrefslogtreecommitdiff
path: root/3889/CH5/EX5.5/Ex5_5.sce
diff options
context:
space:
mode:
authorprashantsinalkar2018-02-03 11:01:52 +0530
committerprashantsinalkar2018-02-03 11:01:52 +0530
commit7bc77cb1ed33745c720952c92b3b2747c5cbf2df (patch)
tree449d555969bfd7befe906877abab098c6e63a0e8 /3889/CH5/EX5.5/Ex5_5.sce
parentd1e070fe2d77c8e7f6ba4b0c57b1b42e26349059 (diff)
downloadScilab-TBC-Uploads-7bc77cb1ed33745c720952c92b3b2747c5cbf2df.tar.gz
Scilab-TBC-Uploads-7bc77cb1ed33745c720952c92b3b2747c5cbf2df.tar.bz2
Scilab-TBC-Uploads-7bc77cb1ed33745c720952c92b3b2747c5cbf2df.zip
Added new codeHEADmaster
Diffstat (limited to '3889/CH5/EX5.5/Ex5_5.sce')
-rw-r--r--3889/CH5/EX5.5/Ex5_5.sce26
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")
+