From 8ac15bc5efafa2afc053c293152605b0e6ae60ff Mon Sep 17 00:00:00 2001 From: Siddharth Agarwal Date: Tue, 3 Sep 2019 18:27:40 +0530 Subject: Xcos examples from textbooks and for blocks --- Working_Examples/3885/CH5/EX5.6/Ex5_6.sci | 27 +++++++++++++++++++++++++++ 1 file changed, 27 insertions(+) create mode 100755 Working_Examples/3885/CH5/EX5.6/Ex5_6.sci (limited to 'Working_Examples/3885/CH5/EX5.6') diff --git a/Working_Examples/3885/CH5/EX5.6/Ex5_6.sci b/Working_Examples/3885/CH5/EX5.6/Ex5_6.sci new file mode 100755 index 0000000..d513434 --- /dev/null +++ b/Working_Examples/3885/CH5/EX5.6/Ex5_6.sci @@ -0,0 +1,27 @@ +//control systems by Nagoor Kani A +//Edition 3 +//Year of publication 2015 +//Scilab version 6.0.0 +//operating systems windows 10 +// Example 5.6 + +clc; +clear; +s=poly(0,'s') +a=(s^7)+(5*s^6)+(9*s^5)+(9*s^4)+(4*s^3)+(20*s^2)+(36*s)+36 +b=coeff(a) +n=length(b) +R=routh_t(a) +disp(R,'the routh array is;') +//divide characterstic equation by auxilary polynomial to get quotient polynomial +//routh table for quotient poly nomial +a1=(s^3)+(5*s^2)+(9*s)+9 //quotient poly nomial +b1=coeff(a1) +n1=length(b1) +R1=routh_t(a1) +disp(R1,'the routh array for quotient poly nomial is;') +ap=s^4+4 //auxillary polynomial +r=roots(ap) +disp(r,'the roots are') +disp('the system is unstable') +disp('two roots on right half of s plane and five roots lie in left half of s plane') -- cgit