23 files changed, 430 insertions, 0 deletions
diff --git a/3885/CH5/EX5.1/Ex5_1.sci b/3885/CH5/EX5.1/Ex5_1.sci
new file mode 100644
index 000000000..15a526685
--- /dev/null
+++ b/3885/CH5/EX5.1/Ex5_1.sci
@@ -0,0 +1,17 @@
+//control systems by Nagoor Kani A
+//Edition 3
+//Year of publication 2015
+//Scilab version 6.0.0
+//operating systems windows 10
+// Example 5.1
+
+clc;
+clear; 
+s=poly(0,'s')
+a=(s^4)+(8*s^3)+(18*s^2)+(16*s)+5
+b=coeff(a)
+n=length(b)
+R=routh_t(a)
+disp(R,'the routh array is;')
+disp('since there is no sign change  the system is stable')
+disp('all the four roots lie left half of the s plane')
diff --git a/3885/CH5/EX5.12/Ex5_12.sci b/3885/CH5/EX5.12/Ex5_12.sci
new file mode 100644
index 000000000..c670d7cf1
--- /dev/null
+++ b/3885/CH5/EX5.12/Ex5_12.sci
@@ -0,0 +1,17 @@
+//control systems by Nagoor Kani A
+//Edition 3
+//Year of publication 2015
+//Scilab version 6.0.0
+//operating systems windows 10
+// Example 5.12
+
+clc;
+clear; 
+s=poly(0,'s')
+//the close loop tranfer function is K(1-s)/s*(s^2+5s+9)+K(1-s)
+a=(s^3)+(5*s^2)+(8*s)+1//the charater stic equation is assuming K=1
+b=coeff(a)
+n=length(b)
+R=routh_t(a)
+disp(R,'the routh array is;')
+disp('the valus of K lies between 0 to7.5 for the system to be stable')
diff --git a/3885/CH5/EX5.13/Ex5_13.sci b/3885/CH5/EX5.13/Ex5_13.sci
new file mode 100644
index 000000000..342823134
--- /dev/null
+++ b/3885/CH5/EX5.13/Ex5_13.sci
@@ -0,0 +1,17 @@
+//control systems by Nagoor Kani A
+//Edition 3
+//Year of publication 2015
+//Scilab version 6.0.0
+//operating systems windows 10
+// Example 5.13
+
+clc;
+clear;
+s=poly(0,'s')//defines s as  poly nomial variable
+h=syslin('c',(230/(s*(s+2)*(s+10))))//the value of K lies between 0 to 240 .the given transfer function assigned to variable h
+nyquist(h)
+show_margins(h,'nyquist')
+//calculation of K
+//when -0.00417K=-1the contour passes through -1+j0 and the corrospondig value of K is the limiting value of k for stabilty 
+K=1/0.00417
+disp(K,'the value of K is')
diff --git a/3885/CH5/EX5.14/Ex5_14.sci b/3885/CH5/EX5.14/Ex5_14.sci
new file mode 100644
index 000000000..bc9317727
--- /dev/null
+++ b/3885/CH5/EX5.14/Ex5_14.sci
@@ -0,0 +1,16 @@
+//control systems by Nagoor Kani A
+//Edition 3
+//Year of publication 2015
+//Scilab version 6.0.0
+//operating systems windows 10
+// Example 5.14
+
+clc;
+clear;
+s=poly(0,'s')//defines s as  poly nomial variable
+h=syslin('c',(0.6*(1+s)^2/(s^3)))//the system is stable for K vlaue greater than 0.5. the given transfer function assigned to variable h
+scf()
+nyquist(h)
+show_margins(h,'nyquist')
+//when K>0.5 -1+j0 is encircled in both clockwise and anticlockwise direction one time. so the system is stable
+disp('the system is stable for K vlaue greater than 0.5.')
diff --git a/3885/CH5/EX5.15/Ex5_15.sci b/3885/CH5/EX5.15/Ex5_15.sci
new file mode 100644
index 000000000..f4eef5eed
--- /dev/null
+++ b/3885/CH5/EX5.15/Ex5_15.sci
@@ -0,0 +1,16 @@
+//control systems by Nagoor Kani A
+//Edition 3
+//Year of publication 2015
+//Scilab version 6.0.0
+//operating systems windows 10
+// Example 5.15
+
+clc;
+clear;
+s=poly(0,'s')//defines s as  poly nomial variable
+h=syslin('c',(1+4*s)/(s*(1+s)*(1+2*s)))//the given transfer function assigned to variable h
+scf()
+nyquist(h)
+show_margins(h,'nyquist')
+disp('the closed loop system is unstable ')
+disp('two poles of closed loop system are lying on right half of s plane')
diff --git a/3885/CH5/EX5.16/Ex5_16.sci b/3885/CH5/EX5.16/Ex5_16.sci
new file mode 100644
index 000000000..25d98515e
--- /dev/null
+++ b/3885/CH5/EX5.16/Ex5_16.sci
@@ -0,0 +1,21 @@
+//control systems by Nagoor Kani A
+//Edition 3
+//Year of publication 2015
+//Scilab version 6.0.0
+//operating systems windows 10
+// Example 5.16
+
+clc;
+clear;
+s=poly(0,'s')//defines s as  poly nomial variable
+h=syslin('c',(7*(1+0.5*s)*(1+s)/(1+10*s)*(s-1)))//system is stable for K>6.....the given transfer function assigned to variable h
+scf()
+nyquist(h)
+show_margins(h,'nyquist')
+//calculation of K
+//when -0.01667K=-1the contur passes through -1=j0and the corrospong value of K is the limiting value of k for stabilty 
+K=1/0.1667
+disp(K,'the value of K is')
+disp('the open loop system is unstable')
+disp('for stability of closed loop system K>6')
+
diff --git a/3885/CH5/EX5.17/Ex5_17.sci b/3885/CH5/EX5.17/Ex5_17.sci
new file mode 100644
index 000000000..3eb81997a
--- /dev/null
+++ b/3885/CH5/EX5.17/Ex5_17.sci
@@ -0,0 +1,15 @@
+//control systems by Nagoor Kani A
+//Edition 3
+//Year of publication 2015
+//Scilab version 6.0.0
+//operating systems windows 10
+// Example 5.17
+
+clc;
+clear;
+s=poly(0,'s')//defines s as  poly nomial variable
+h=syslin('c',(5/(s*(1-s))))//the given transfer function assigned to variable h
+scf()
+nyquist(h)
+show_margins(h,'nyquist')
+disp('both open loop and closed loop system are unstable)
diff --git a/3885/CH5/EX5.18/Ex5_18.sci b/3885/CH5/EX5.18/Ex5_18.sci
new file mode 100644
index 000000000..d5327b765
--- /dev/null
+++ b/3885/CH5/EX5.18/Ex5_18.sci
@@ -0,0 +1,16 @@
+//control systems by Nagoor Kani A
+//Edition 3
+//Year of publication 2015
+//Scilab version 6.0.0
+//operating systems windows 10
+// Example 5.28
+
+clc;
+clear;
+s=poly(0,'s')//defines s as  poly nomial variable
+h=syslin('c',(s+2)/(s+1)*(s-1))//the given transfer function assigned to variable h
+scf()
+nyquist(h)
+show_margins(h,'nyquist')
+disp('open loop system is unstable')
+disp('closed loop system is stable;' )
diff --git a/3885/CH5/EX5.2/Ex5_2.sci b/3885/CH5/EX5.2/Ex5_2.sci
new file mode 100644
index 000000000..7e9f4d97c
--- /dev/null
+++ b/3885/CH5/EX5.2/Ex5_2.sci
@@ -0,0 +1,20 @@
+//control systems by Nagoor Kani A
+//Edition 3
+//Year of publication 2015
+//Scilab version 6.0.0
+//operating systems windows 10
+// Example 5.2
+
+clc;
+clear; 
+s=poly(0,'s')
+a=(s^6)+(2*s^5)+(8*s^4)+(12*s^3)+(20*s^2)+(16*s)+16
+b=coeff(a)
+n=length(b)
+R=routh_t(a)
+disp(R,'the routh array is;')
+ap=s^4+6*s^2+8
+r=roots(ap)
+disp(r,'the roots of auxilary polynomial;')
+disp('the system is marginally stable;')
+disp('four roots lying in imaginary axis ;')
diff --git a/3885/CH5/EX5.22/Ex5_22.sci b/3885/CH5/EX5.22/Ex5_22.sci
new file mode 100644
index 000000000..a411018c2
--- /dev/null
+++ b/3885/CH5/EX5.22/Ex5_22.sci
@@ -0,0 +1,19 @@
+//control systems by Nagoor Kani A
+//Edition 3
+//Year of publication 2015
+//Scilab version 6.0.0
+//operating systems windows 10
+// Example 5.22
+
+clc;
+clear;
+s=poly(0,'s')//defines s as  poly nomial variable
+G=syslin('c',(1/(s*(s^2+4*s+13))))//the given transfer function assigned to variable G assume K=1
+scf()
+evans(G)
+//calculation of K
+disp('the characterstic equation is given by : s^2+4*s^2+13*s+K')
+//put s=jw and equate real and imaginary parts
+//K=4*w^2
+K=4*13
+disp(K,'the value of K is ')
diff --git a/3885/CH5/EX5.23/Ex5_23.sci b/3885/CH5/EX5.23/Ex5_23.sci
new file mode 100644
index 000000000..feea17484
--- /dev/null
+++ b/3885/CH5/EX5.23/Ex5_23.sci
@@ -0,0 +1,19 @@
+//control systems by Nagoor Kani A
+//Edition 3
+//Year of publication 2015
+//Scilab version 6.0.0
+//operating systems windows 10
+// Example 5.23
+
+clc;
+clear;
+s=poly(0,'s')//defines s as  poly nomial variable
+h=syslin('c',(48/(s*(s+2)*(s+4))))//the given transfer function assigned to variable h assume K=1
+scf()
+evans(h)
+//calculation of K
+disp('the characterstic equation is given by : s^3+6*s^2+8*s+K')
+//put s=jw and equate real and imaginary parts
+//K=4*w^2
+K=6*8
+disp(K,'the value of K is ')
diff --git a/3885/CH5/EX5.24/Ex5_24.sci b/3885/CH5/EX5.24/Ex5_24.sci
new file mode 100644
index 000000000..91800d858
--- /dev/null
+++ b/3885/CH5/EX5.24/Ex5_24.sci
@@ -0,0 +1,22 @@
+//control systems by Nagoor Kani A
+//Edition 3
+//Year of publication 2015
+//Scilab version 6.0.0
+//operating systems windows 10
+// Example 5.24
+
+clc;
+clear;
+s=poly(0,'s')//defines s as  poly nomial variable
+h=syslin('c',(s+9)/(s*(s^2+4*s+11)))//the given transfer function assigned to variable h assume K=1
+scf()
+evans(h)
+//the characterstic equation is (s^3+4s^2+11s)+Ks+9K
+//put s=jw and equating real and imaginary parts to calculate K
+K=(4*(4.4)^2)/9//the value of w is 4.4
+disp(K,'the value of K is;')
+
+
+
+
+
diff --git a/3885/CH5/EX5.25/Ex5_25.sci b/3885/CH5/EX5.25/Ex5_25.sci
new file mode 100644
index 000000000..f3c7eb9d1
--- /dev/null
+++ b/3885/CH5/EX5.25/Ex5_25.sci
@@ -0,0 +1,17 @@
+//control systems by Nagoor Kani A
+//Edition 3
+//Year of publication 2015
+//Scilab version 6.0.0
+//operating systems windows 10
+// Example 5.25
+
+clc;
+clear;
+s=poly(0,'s')//defines s as  poly nomial variable
+h=syslin('c',(1/(s*(s+4)*(s^2+4*s+20))))// the given transfer function assigned to variable h assume K=1  
+scf()
+evans(h)
+//the characterstic equation is (s^4+8s^3+36s^2+80s)+K
+//put s=jw and equating real and imaginary parts to calculate K
+K=-(3.2)^4+36*(3.2)^2//the value of w is 3.2
+disp(K,'yhe value of K is;')
diff --git a/3885/CH5/EX5.26/Ex5_26.sci b/3885/CH5/EX5.26/Ex5_26.sci
new file mode 100644
index 000000000..d8e13f3b9
--- /dev/null
+++ b/3885/CH5/EX5.26/Ex5_26.sci
@@ -0,0 +1,17 @@
+//control systems by Nagoor Kani A
+//Edition 3
+//Year of publication 2015
+//Scilab version 6.0.0
+//operating systems windows 10
+// Example 5.26
+
+clc;
+clear;
+s=poly(0,'s')//defines s as  poly nomial variable
+h=syslin('c',(s+15)/(s*(s+1)*(s+5)))//the given transfer function assigned to variable h assume K=1
+scf()
+evans(h)
+//the characterstic equation is (s^3+6s^2+5s)+Ks+15K
+//put s=jw and equating real and imaginary parts to calculate K
+K=30/(-4.5)
+disp(K,'yhe value of K is;')
diff --git a/3885/CH5/EX5.27/Ex5_27.sci b/3885/CH5/EX5.27/Ex5_27.sci
new file mode 100644
index 000000000..36f588de1
--- /dev/null
+++ b/3885/CH5/EX5.27/Ex5_27.sci
@@ -0,0 +1,18 @@
+//control systems by Nagoor Kani A
+//Edition 3
+//Year of publication 2015
+//Scilab version 6.0.0
+//operating systems windows 10
+// Example 5.27
+
+clc;
+clear;
+s=poly(0,'s')//defines s as  poly nomial variable
+K=poly(0,'K')
+h=syslin('c',(s^2+6*s+25)/(s*(s+1)*(s+2)))//the given transfer function assigned to variable h assume K=1
+scf()
+evans(h)
+//the characterstic equation is (s^3+(3+K)s^2+(2+6K)s)25K
+//put s=jw and equating real and imaginary parts to calculate K
+a=roots((-6*K^2)+5*K-6)
+disp(a,'the value of K is')
diff --git a/3885/CH5/EX5.28/Ex5_28.sci b/3885/CH5/EX5.28/Ex5_28.sci
new file mode 100644
index 000000000..64e45a9dd
--- /dev/null
+++ b/3885/CH5/EX5.28/Ex5_28.sci
@@ -0,0 +1,17 @@
+//control systems by Nagoor Kani A
+//Edition 3
+//Year of publication 2015
+//Scilab version 6.0.0
+//operating systems windows 10
+// Example 5.28
+
+clc;
+clear;
+s=poly(0,'s')//defines s as  poly nomial variable
+h=syslin('c',(1/(s*(s^2+6*s+10))))//the given transfer function assigned to variable h assume K=1
+scf()
+evans(h)
+//the characterstic equation is (s^3+6s^2+10s)+K
+//put s=jw and equating real and imaginary parts to calculate K
+K=6*(3.2)^2//the valuw of w is 3.2
+disp(K,'the value of K is ')
diff --git a/3885/CH5/EX5.3/Ex5_3.sci b/3885/CH5/EX5.3/Ex5_3.sci
new file mode 100644
index 000000000..b295151fd
--- /dev/null
+++ b/3885/CH5/EX5.3/Ex5_3.sci
@@ -0,0 +1,16 @@
+//control systems by Nagoor Kani A
+//Edition 3
+//Year of publication 2015
+//Scilab version 6.0.0
+//operating systems windows 10
+// Example 5.3
+
+clc;
+clear; 
+s=poly(0,'s')
+a=(s^5)+(s^4)+(2*s^3)+(2*s^2)+(3*s)+5
+b=coeff(a)
+n=length(b)
+R=routh_t(a)
+disp(R,'the routh array is;')
+disp('the system is unstable')
diff --git a/3885/CH5/EX5.4/Ex5_4.sci b/3885/CH5/EX5.4/Ex5_4.sci
new file mode 100644
index 000000000..972e7e18d
--- /dev/null
+++ b/3885/CH5/EX5.4/Ex5_4.sci
@@ -0,0 +1,17 @@
+//control systems by Nagoor Kani A
+//Edition 3
+//Year of publication 2015
+//Scilab version 6.0.0
+//operating systems windows 10
+// Example 5.4
+
+clc;
+clear; 
+s=poly(0,'s')
+a=(9*s^5)-(20*s^4)+(10*s^3)-(s^2)-(9*s)-10
+b=coeff(a)
+n=length(b)
+R=routh_t(a)
+disp(R,'the routh array is;')
+disp('there are 3 sign changes in first column of routh array ')
+disp('three roots lie on right side of s plane so the system is unstable')
diff --git a/3885/CH5/EX5.5/Ex5_5.sci b/3885/CH5/EX5.5/Ex5_5.sci
new file mode 100644
index 000000000..78c85847f
--- /dev/null
+++ b/3885/CH5/EX5.5/Ex5_5.sci
@@ -0,0 +1,21 @@
+//control systems by Nagoor Kani A
+//Edition 3
+//Year of publication 2015
+//Scilab version 6.0.0
+//operating systems windows 10
+// Example 5.5
+
+clc;
+clear; 
+s=poly(0,'s')
+a=(s^7)+(9*s^6)+(24*s^5)+(24*s^4)+(24*s^3)+(24*s^2)+(23*s)+15
+b=coeff(a)
+n=length(b)
+R=routh_t(a)
+disp(R,'the routh array is;')
+disp('there is sign change in first column of routh array so ths system is unstable;')
+ae=s^4+s^2+1
+r=roots(ae)
+disp(r,'the roots of auxilary equation are')
+disp('two roots lie on right half of splane five roots lie on left half of s plane')
+
diff --git a/3885/CH5/EX5.6/Ex5_6.sci b/3885/CH5/EX5.6/Ex5_6.sci
new file mode 100644
index 000000000..d5134347f
--- /dev/null
+++ b/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')
diff --git a/3885/CH5/EX5.7/Ex5_7.sci b/3885/CH5/EX5.7/Ex5_7.sci
new file mode 100644
index 000000000..6111bfe21
--- /dev/null
+++ b/3885/CH5/EX5.7/Ex5_7.sci
@@ -0,0 +1,20 @@
+//control systems by Nagoor Kani A
+//Edition 3
+//Year of publication 2015
+//Scilab version 6.0.0
+//operating systems windows 10
+// Example 5.7
+
+clc;
+clear; 
+s=poly(0,'s')
+a=(s^5)+(4*s^4)+(8*s^3)+(8*s^2)+(7*s)+4
+b=coeff(a)
+n=length(b)
+R=routh_t(a)
+disp(R,'the routh array is;')
+ap=s^2+1
+r=roots(ap)
+disp(r,'the roots are')
+disp ('the roots of auxillary equation are in imagianry axis so the system is marginally stable')
+disp('three roots lie in left half of s plane')
diff --git a/3885/CH5/EX5.8/Ex5_8.sci b/3885/CH5/EX5.8/Ex5_8.sci
new file mode 100644
index 000000000..9de51e2a7
--- /dev/null
+++ b/3885/CH5/EX5.8/Ex5_8.sci
@@ -0,0 +1,28 @@
+//control systems by Nagoor Kani A
+//Edition 3
+//Year of publication 2015
+//Scilab version 6.0.0
+//operating systems windows 10
+// Example 5.8
+
+clc;
+clear; 
+s=poly(0,'s')
+a=(s^6)+(s^5)+(3*s^4)+(3*s^3)+(3*s^2)+(2*s)+1
+b=coeff(a)
+n=length(b)
+R=routh_t(a)
+disp(R,'the routh array is;') 
+//characterstic polynomial can be expressed as product of auxillary polynomial and quotient polynomial
+//divide characterstic equation by auxilary polynomial to get quotient polynomial
+//routh table for quotient poly nomial
+a1=(s^4)+(s^3)+(2*s^2)+(2*s)+1 //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^2+1
+r=roots(ap)
+disp(r,'the roots are')
+disp('the system is unstable')
+disp('two roota on imaginary axis ,two roots on right half of s plane and two roots lie in left half of s plane')
diff --git a/3885/CH5/EX5.9/Ex5_9.sci b/3885/CH5/EX5.9/Ex5_9.sci
new file mode 100644
index 000000000..d9b5660bd
--- /dev/null
+++ b/3885/CH5/EX5.9/Ex5_9.sci
@@ -0,0 +1,17 @@
+//control systems by Nagoor Kani A
+//Edition 3
+//Year of publication 2015
+//Scilab version 6.0.0
+//operating systems windows 10
+// Example 5.9
+
+clc;
+clear; 
+s=poly(0,'s')
+//the close loop tranfer function is K/s*(s+1)*(s+2)+K
+a=(s^3)+(3*s^2)+(2*s)+1//the characterstic equation assuming K=1
+b=coeff(a)
+n=length(b)
+R=routh_t(a)
+disp(R,'the routh array is;')
+disp('the valus of K lies between 0 to6 for the system to be stable')