Scilab Companion
Scilab Companion
Scilab Companion
Created by
Pooja Naik
B. E.
Instrumentation Engineering
Watumaull College of Electronics
College Teacher
Prof. Ashutosh Sharma
Cross-Checked by
Chaitanya
Author: M. N. Bandyopadhyay
Edition: 1
Year: 2009
ISBN: 9788120319547
1
Scilab numbering policy used in this document and the relation to the
above book.
For example, Exa 3.51 means solved example 3.51 of this book. Sec 2.3 means
a scilab code whose theory is explained in Section 2.3 of the book.
2
Contents
3
List of Scilab Codes
4
Exa 5.3 Routh array . . . . . . . . . . . . . . . . . . 27
Exa 5.4 Routh array . . . . . . . . . . . . . . . . . . 28
Exa 5.5 Routh array . . . . . . . . . . . . . . . . . . 28
Exa 5.6 Routh array . . . . . . . . . . . . . . . . . . 29
Exa 6.2 Root locus in scilab . . . . . . . . . . . . . . 31
Exa 6.2.2 location of the root locus between poles and
zeros . . . . . . . . . . . . . . . . . . . . . 32
Exa 6.3 Root locus . . . . . . . . . . . . . . . . . . . 32
Exa 6.4 root locus . . . . . . . . . . . . . . . . . . . 33
Exa 6.5 Root locus . . . . . . . . . . . . . . . . . . . 34
Exa 6.6 Root locus . . . . . . . . . . . . . . . . . . . 34
Exa 6.7 Root locus . . . . . . . . . . . . . . . . . . . 35
Exa 7.3.1a Bode plot . . . . . . . . . . . . . . . . . . . 36
Exa 7.3.1b Bode plot . . . . . . . . . . . . . . . . . . . 36
Exa 8.1 Nyquist plot . . . . . . . . . . . . . . . . . . 38
Exa 8.2 Nyquist plot . . . . . . . . . . . . . . . . . . 38
Exa 8.3 Nyquist plot . . . . . . . . . . . . . . . . . . 39
Exa 9.1 compensation in open loop control system . 41
Exa 10.1.1 Mass dashpot and spring arrangement . . . 44
Exa 10.3 determination of quadratic form . . . . . . . 45
Exa 10.4 Lipunovs method . . . . . . . . . . . . . . . 45
Exa 11.6 Jurys stability test . . . . . . . . . . . . . . 47
Exa 11.9.2a stability of linear continuous system . . . . . 48
Exa 11.9.2b stability of linear continuous system . . . . . 49
Exa 11.9.3 Schurcohn stability test . . . . . . . . . . . 50
Exa 15.2 Time domain specifications of second order
system . . . . . . . . . . . . . . . . . . . . . 51
Exa 15.4 transfer function of gyroscope . . . . . . . . 52
Exa 15.5 Transfer function of system . . . . . . . . . 52
Exa 15.6 comparison of sensitivities of two systems . 53
Exa 15.7 To find bandwidth of the transfer function . 53
Exa 15.8 Bandwidth of the transfer function . . . . . 54
Exa 15.9 Nyquist plot . . . . . . . . . . . . . . . . . . 55
Exa 15.12 solution of polynolynomial equation . . . . . 55
Exa 15.13 Time domain specifications . . . . . . . . . 56
Exa 15.14 Position servomotor . . . . . . . . . . . . . 57
Exa 15.15 steady output speed of DC motor . . . . . . 57
Exa 15.26 Routh array . . . . . . . . . . . . . . . . . . 58
5
Exa 15.27 To check reachability of the system . . . . . 59
Exa 15.28 Determine the stability of the system . . . . 60
Exa 15.31 Time domain specifications of second order
system . . . . . . . . . . . . . . . . . . . . . 60
Exa 15.33 Lipunovs method . . . . . . . . . . . . . . . 61
Exa 15.34 Find bandwidth of the transfer function . . 62
Exa 15.36 Impulse response of the transfer function . . 62
Exa 15.37 step response of the transfer function . . . . 63
Exa 15.38 Roots of characteristic equation . . . . . . . 63
Exa 15.39 Bode plot . . . . . . . . . . . . . . . . . . . 64
Exa 15.40 Nyquist plot . . . . . . . . . . . . . . . . . . 64
Exa 15.41 Nyquist plot . . . . . . . . . . . . . . . . . . 65
Exa 15.42 Bode plot . . . . . . . . . . . . . . . . . . . 65
Exa 15.43 Eigen values of matrix . . . . . . . . . . . . 66
Exa 15.44 State space representation of LTI system . . 66
Exa 15.45 Covariant matrix of A . . . . . . . . . . . . 67
Exa 15.49 Root locus of the transfer function . . . . . 67
Exa 15.50 Bode plot . . . . . . . . . . . . . . . . . . . 68
Exa 15.53 Statespace model of the differential equation 68
Exa 15.54 Nyquist plot . . . . . . . . . . . . . . . . . . 69
Exa 15.57 determination of zeta and wn . . . . . . . . 69
Exa 15.59 transfer function of signal flow graph . . . . 70
Exa 15.60 root locus . . . . . . . . . . . . . . . . . . . 71
6
Chapter 2
1 // Example 2 . 1
2 // s o l v i n g d i f f e r e n t i a l e q u a t i o n i n s c i l a b
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 function ydot = f (t , y )
7 ydot =(10/4) -(3* y /4)
8 endfunction
9 y0 =1;
10 t0 =0;
11 t =0:5:10;
12 y = ode ( y0 , t0 ,t , f )
13 // s i n c e ” t = 0 ; 5 ; 1 0 ”
14 // t h e a n s w e r i s c a l c u l a t e d f o r t = 0 : 5 : 1 0 ”
15 // t h u s t h e v a l u e o f ” y ” can be c a l c u l a t e d a t any
value of ” t ”.
7
Scilab code Exa 2.2 inverse of laplace transform using scilab
1 // Example 2 . 2
2 // I n v e r s e l a p l a c e t r a n s f o r m u s i n g s c i l a b
3 clear ; clc ;
4 xdel ( winsid () ) ;
5 s = %s ;
6 num =( s +6) ;
7 den =( s ^2+2* s +10) ;
8 F1 = syslin ( ’ c ’ ,num , den )
9 F = pfss ( F1 )
10 // s i n c e p f s s ( F1 ) i s n o t a b l e t o f a c t o r i s e F1 ,
therefore ,
11 // R e w r i t i n g n u m e r a t o r as , ( s +6)=( s +1+5) ;
12 // R e w r i t i n g t h e d e n o m i n a t o r as , ( s ˆ2+2∗ s +6)=( s +1)
ˆ2+3ˆ2;
13 disp ( ”F = [ ( ( s +1) / ( s +1) ˆ2+3ˆ2) + ( 5 / 3 ) ∗ ( 3 / ( s +1) ˆ2+3ˆ2) ] ”
)
14 // From t h e s t a n d a r d f o r m u l a o f i n v e r s e l a p l a c e
transform ;
15 // ( s +1) / ( s +1)ˆ2+3ˆ2=%eˆ− t ∗ ( c o s 3 t ) ;
16 // ( 5 / 3 ) ∗ ( 3 / ( s +1) ˆ2+3ˆ2) = ( 5 / 3 ) ∗%eˆ− t ∗ ( s i n 3 t ) ;
17 disp ( ” f ( t ) =(%eˆ− t ) ∗ ( c o s 3 t ) + ( 5 / 3 ) ∗ ( %eˆ− t ) ∗ ( s i n 3 t ) ” )
1 // Example 2 . 3
2 // c o m p u t i n g i n i t i a l v a l u e u s i n g s c i l a b
8
3 clear ; clc ;
4 xdel ( winsid () ) ;
5 s = %s
6 n3 =(3* s +2)
7 d3 = s *( s ^2+4* s +5)
8 F = n3 / d3
9 // A p p l y i n g i n i t i a l v a l u e t h e o r e m
10 // when l i m i t ” s ” t e n d s t o i n f i n i t y , f i n a l v a l u e o f ”
F” becomes ”0”
11
12 disp (0 , ” F i n a l v a l u e=” )
1 // Example s e c 2 . 4 . 2
2 // e i g e n v a l u e s
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 A =[0 6 -5;1 0 2;3 2 4]
7 B = spec ( A )
8 disp (B , ” E i g e n v a l u e s=” )
1 // Example 2 . 4
2 // c o m p u t i n g f ’ ( 0 + ) and f ’ ’ ( 0 + ) u s i n g s c i l a b
3 clear ; clc ;
4 xdel ( winsid () ) ;
9
5
6 s = %s ;
7 n4 =(4* s +1) ;
8 d4 = s *( s ^2+4* s +5) ;
9 F = n4 / d4
10 // As p e r i n i t i a l v a l u e theorem , l i m i t ” t ” t e n d s t o
z e r o and l i m i t ” s ” t e n d s t o i n f i n i t y
11
12 // f o r f ’ ( 0 + )
13 F1 = s * F +0
14
15 for s = %inf
16 disp ( ” f ’ ’ ( 0 + ) =4” )
17 end
18 // f o r f ’ ’ ( 0 + )
19 s = %s ;
20 F2 =(( s *( F1 ) ) -(4) )
21
22 for s = %inf
23 disp ( ” f ’ ’ ’ ’ ( 0 + ) =−15” )
24 end
1 // Example s e c 2 . 4 . a
2 // e i g e n v a l u e s
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 A =[1 -1;0 -1]
7 B = spec ( A )
8 disp (B , ” E i g e n v a l u e s=” )
10
Scilab code Exa 2.5 computing initial value of function F in scilab
1 // Example 2 . 5
2 // c o m p u t i n g f i n a l v a l u e o f f u n c t i o n F u s i n g s c i l a b
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 s = %s ;
7 n5 =(8* s +5) ;
8 d5 = s *( s +1) *( s ^2+4* s +5) ;
9 F = n5 / d5
10 F1 = s * F
11 // f o r f i n a l v a l u e l i m i t ” t ” t e n d s t o i n f i n i t y and
l i m i t ” s ” tends to zero .
12 //When s =0 , t h e v a l u e o f F1 w i l l be ” ( 5 / 5 ) =1”
13 for s =0
14 disp ( ” F i n a l v a l u e =1” )
15 end
Scilab code Exa 2.6 computing final value of function F using scilab
1 // Example 2 . 6
2 // c o m p u t i n g f i n a l value of function F using s c i l a b
3 clear ; clc ;
4 xdel ( winsid () ) ;
5 s = %s ;
6 n6 =(5) ;
7 d6 = s *( s ^2+49) ;
11
8 F = n6 / d6
9 F1 = s * F
10
11 disp ( ”THE SYSTEM HAS POLES ON IMAGINARY AXIS .
THEREFORE f ( t ) HAS NO FINAL VALUE” )
1 // Example 2 . 7
2 // I n v e r s e l a p l a c e t r a n s f o r m o f ” 2 / ( s ˆ 2 ∗ ( s +1) ) ” u s i n g
scilab
3 clear ; clc ;
4 xdel ( winsid () ) ;
5 s = %s ;
6 num =2;
7 den =( s ^2) *( s +1) ;
8 F1 = syslin ( ’ c ’ ,num , den )
9 F = pfss ( F1 )
10 // from t h e p a r t i a l f r a c t i o n d e c o m p o s i t i o n , t a k i n g
o u t 2 a s common term .
11 // The r e s u l t would be i n t h e form o f ”F ( s ) =2∗(1/ s
ˆ2 −1/ s +(1/ s +1) ) ”
12 disp ( ”F ( s ) = 2 ∗ ( ( 1 / s ˆ 2 ) −(1/ s ) +(1/( s +1) ) ) ” )
13 // From t h e s t a n d a r d f o r m u l a o f i n v e r s e l a p l a c e
transform ;
14 // ( 1 / s ˆ 2 )=t ; ( 1 / s ) = 1 ; ( 1 / ( s +1) )=%eˆ− t
15 disp ( ” f ( t ) =2∗( t −1+eˆ− t ) ” )
12
1 // Example 2 . 8
2 //Z t r a n s f o r m o f t h e s i g n a l x ( n ) = ( 0 . 5 ) ˆ n∗u ( n )
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 // u ( n ) i s u n i t s t e p i n p u t
7 n =2;
8 x =(0.5) ^ n ;
9 m =1;
10 w =1;
11 phi = tand (0) ;
12 a =1;
13 theta = tand (45) ;
14 [ cxz ]= czt (x ,m ,w , phi ,a , theta )
1 // Example 2 . 9
2 //Z t r a n s f o r m o f t h e s i g n a l x ( n ) =(a ) ˆ n∗u ( n ) +(b ) ˆ n∗u
(−n −1) .
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 // u ( n ) i s u n i t s t e p i n p u t
7 // a =0.5 and b =0.6
8 n1 =2;
9 x1 =(0.5) ^ n1 ;
10 m1 =1;
11 w1 =1;
12 phi1 = tand (0) ;
13 a1 =1;
14 theta1 = tand (45) ;
15 [ X1 ]= czt ( x1 , m1 , w1 , phi1 , a1 , theta1 )
13
16 n2 =2;
17 x2 =(0.6) ^ n2 ;
18 m2 =1;
19 w2 = -1;
20 phi2 = tand ( -45) ;
21 a2 =1;
22 theta2 = tand (45) ;
23 [ X2 ]= czt ( x2 , m2 , w2 , phi2 , a2 , theta2 )
24 X = X1 + X2 ;
25 disp (X , ” a n s=” )
1 // Example 2 . 1 0
2 // i n v e r s e Z t r a n s f o r m o f 1/(1 − a ∗ z ˆ −1)
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 // a=1
7 function y = f ( z ) ;
8 y = z /( z -1) // upon s i m p l i f i c a t i o n o f t h e g i v e n
equation
9 endfunction
10 intc (1+ %i ,2 - %i , f )
1 // Example 2 . 1 1
2 // i n v e r s e z t r a n s f o r m by power s e r i e s e x p a n s i o n
14
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 z = %z ;
7 num =2;
8 den =(2 -(3* z ^ -1) + z ^ -2) ;
9 X = syslin ( ’ c ’ , num / den )
10 // d i v i d i n g t h e n u m e r a t o r and d e n o m i n a t o r by 2
11 num1 =1;
12 den1 =1 -(1.5*( z ^ -1) ) +(0.5*( z ^ -2) ) ;
13 X1 = syslin ( ’ c ’ ,( num1 ) /( den1 ) )
14 // when mod ( z ) >1
15 // d e v e l o p i n g s e r i e s e x p a n s i o n i n n e g a t i v e power o f z
16 A1 =( num1 ) -( den1 )
17 // m u l t i p l y i n g t h e den2 by 1 and s u b t r a c t i n g i t
from num1
18 B1 =((1.5*( z ^ -1) ) -(0.5*( z ^ -2) ) ) -((1.5*( z ^ -1) ) * den1 )
19 // m u l t i p l y i n g t h e den2 by ( 1 . 5 ∗ z ˆ −1) and
s u b t r a c t i n g i t from r e m i n d e r o f A1
20 C1 =((1.75* z ^ -2) -(0.75* z ^ -3) ) -((1.75* z ^ -2) * den1 )
21 // m u l t i p l y i n g t h e den2 by ( 1 . 7 5 ∗ z ˆ −2) and
s u b t r a c t i n g i t from r e m i n d e r o f A1
22 D1 =((1.875* z ^ -3) -(0.875* z ^ -4) ) -((1.875* z ^ -3) * den1 )
23 // m u l t i p l y i n g t h e den2 by ( 1 . 8 7 5 ∗ z ˆ −3) and
s u b t r a c t i n g i t from r e m i n d e r o f A1
24 E1 =((1.9375* z ^ -4) -(0.9375* z ^ -5) ) -((1.9375* z ^ -4) * den1
)
25 // m u l t i p l y i n g t h e den2 by ( 1 . 9 3 7 5 ∗ z ˆ −4) and
s u b t r a c t i n g i t from r e m i n d e r o f A1
26 disp ( ” x1 ( n ) = 1 , 1 . 5 , 1 . 7 5 , 1 . 8 7 5 , 1 . 9 3 7 5 , . . . . . . . . ” )
27
28 // when mod ( z ) <0.5
29 // d e v e l o p i n g s e r i e s e x p a n s i o n i n p o s i t i v e power o f z
30 A2 =( num ) -((2* z ^2) * den ) // m u l t i p l y o n g t h e den by 2 ∗ ( z
ˆ 2 ) and s u b t r a c t i n g i t from num
31 B2 = A2 -(6* z ^3* den )
32 // m u l t i p l y o n g t h e den by 2 ∗ ( z ˆ 2 ) and s u b t r a c t i n g it
from A2
15
33 C2 = B2 -(14* z ^4* den )
34 // m u l t i p l y o n g t h e den by 2 ∗ ( z ˆ 2 ) and s u b t r a c t i n g it
from B2
35 D2 = C2 -(30* z ^5* den )
36 // m u l t i p l y o n g t h e den by 2 ∗ ( z ˆ 2 ) and s u b t r a c t i n g i t
from C2
37 E2 = D2 -(62* z ^6* den )
38 // m u l t i p l y o n g t h e den by 2 ∗ ( z ˆ 2 ) and s u b t r a c t i n g i t
from D2
39 disp ( ”X2 ( z ) =2∗ z ˆ2+6∗ z ˆ3+14∗ z ˆ4+30∗ z ˆ5+62∗ z
ˆ 6 + . . . . . . . . ”)
40 disp ( ” x2 ( n ) = . . . . . , 6 2 , 3 0 , 1 4 , 6 , 2 , 0 , 0 ” )
1 // Example 2 . 1 2
2 // i n v e r s e z t r a n s f o r m by p a r t i a l f r a c t i o n method
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 z = %z ;
7 num =1;
8 den =((1 - z ^ -1) ^2) *(1+ z ^ -1) ;
9 X = syslin ( ’ c ’ , num / den )
10 X1 = X / z
11 pfss ( X1 )
12 // by p a r t i a l f r a c t i o n t h e X1 w i l l be f a c t o r i s e d a s
( in terms of z )
13 disp ( ”X( z ) = ( 0 . 2 5 ∗ z / ( z +1) ) + ( 0 . 7 5 ∗ z / ( z −1) ) + ( 0 . 5 ∗ z / ( z
−1) ˆ 2 ) ” )
14 disp ( ”X( z ) = ( 0 . 2 5 / ( 1 + z ˆ −1) ) +(0.75/(1 − z ˆ −1) ) + ( 0 . 5 ∗ z / ( z
−1) ˆ 2 ) ” )
15 // 0 . 2 5 / ( 1 + z ˆ −1) i s t h e z t r a n s f o r m o f ” 0 . 2 5 ∗ ( − 1 ) ˆ n∗
16
u(n)”
16 // ( 0 . 7 5 / ( 1 − z ˆ −1) ) i s t h e z t r a n s f o r m o f ” 0 . 7 5 ∗ u ( n ) ”
17 // ( 0 . 5 ∗ z / ( z −1) ˆ 2 ) i s t h e z t r a n s f o r m o f ” 0 . 5 ∗ n∗u ( n ) ”
18 disp ( ” x ( n ) = 0 . 2 5 ∗ ( ( − 1 ) ˆ n ) ∗u ( n ) +0 .7 5∗ u ( n ) +0.5∗ n∗u ( n ) ” )
17
Chapter 3
1 // Example 3 . 1
2 // t y p e o f t h e s y s t e m
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 // f i g ( 3 . 1 4 )
7 s = %s ;
8 n1 =(2) ;
9 d1 =(( s ) *( s ^2+2* s +2) ) ;
10 A = n1 / d1
11 disp ( ” s i n c e one i n t e g r a t i o n i s b e i n g o b s e r v e d , i t is
TYPE 1 s y s t e m ” )
12
13 // f i g ( 3 . 1 5 )
14 s = %s ;
15 n2 =(5) ;
16 d2 =(( s +2) *( s ^2+2* s +3) ) ;
17 B = n2 / d2
18 disp ( ” s i n c e no i n t e g r a t i o n i s b e i n g o b s e r v e d , i t is
18
TYPE 0 s y s t e m ” )
19
20 // f i g ( 3 . 1 6 )
21 s = %s ;
22 n3 =( s +1) ;
23 d3 =(( s ^2) *( s +2) ) ;
24 C = n3 / d3
25 disp ( ” s i n c e two i n t e g r a t i o n i s b e i n g o b s e r v e d , i t is
TYPE 2 s y s t e m ” )
19
Chapter 4
1 // Example 4 . 1
2 // s t a t e e q u a t i o n
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 A =[0 1; -2 -3]
7 B =[0;1]
8 C =[0]
9 [ Ac Bc U ind ]= canon (A , B ) ;
10 disp ( clean ( Ac ) , ’ Ac= ’ ) ;
11 disp ( clean ( Bc ) , ’ Bc= ’ ) ;
12 disp (U , ’ t r a n s f o r m a t i o n m a t r i x U= ’ ) ;
1 // Example 4 . 3
20
2 // f o r g i v e n m a t r i x ”A” p r o v i n g e i g e n v a l u e s o f ”A”=”
t ˆ−1∗A∗T”
3 clear ; clc ;
4 xdel ( winsid () ) ;
5 A =[0 1 0;0 0 1; -6 -11 -6]
6 P = bdiag ( A ) // e i g e n v a l u e s o f ”A”
7
8 T =[1 1 1; -1 -2 -3; 1 4 9] // vandermode m a t r i x
9 inv ( T )
10
11 A1 = inv ( T ) * A * T // d i a g o n a l c a n o n i c a l form o f A
12
13 // t h u s ”P=A1” i s p r o v e d .
1 // Example s e c 4 . 6 a
2 // e x a m p l e o f c a n o n i c a l form
3 clear ; clc ;
4 xdel ( winsid () ) ;
5 A =[1 2 1;0 1 3;1 1 1];
6 B =[1;0;1];
7 C =[1 0 0;0 1 0;0 0 1]
8 S = cont_mat (A , B )
9 s = %s ;
10 D = s *C - A
11 det ( D )
12
13 // t h e c h a r a c t e r i s t i c e q u a t i o n i . e . d e t (D)=s ˆ3−3∗ s ˆ2−
s −3=0 i s o f t h e form o f
14 // s ˆ3+ a2 ∗Sˆ2+ a1 ∗ s+a0 =0. t h e r e f o r e c o m p a r i n g two
equation .
15
21
16 a2 = -3
17 a1 = -1
18 a0 = -3
19 M =[ a1 a2 1; a2 1 0;1 0 0]
20
21 P=S*M
22 A1 = inv ( P ) * A * P
23 B1 = inv ( P ) * B
1 // Example s e c 4 . 6 b
2 // e x a m p l e o f c a n o n i c a l form
3 clear ; clc ;
4 xdel ( winsid () ) ;
5 A =[1 2 1;0 1 3;1 1 1];
6 B =[1;0;1];
7 C =[1 1 0];
8 V =[ C ; C * A ; C * A ^2]
9
10 D = eye (3 ,3)
11 s = %s
12 E = s *D - A
13 det ( E )
14
15 // t h e c h a r a c t e r i s t i c e q u a t i o n i . e . d e t (E)=s ˆ3−3∗ s ˆ2−
s −3=0 i s o f t h e form o f
16 // s ˆ3+ a2 ∗Sˆ2+ a1 ∗ s+a0 =0. t h e r e f o r e c o m p a r i n g two
equation .
17
18 a2 = -3
19 a1 = -1
20 a0 = -3
22
21 M =[ a1 a2 1; a2 1 0;1 0 0]
22 F=M*V
23 Q = inv ( F )
24 A1 = inv ( Q ) * A * Q
25 B1 = inv ( Q ) * B
26 C1 = C * Q
1 // Example s e c 4 . 8
2 // J o r d a n c a n o n i c a l form
3 clear ; clc ;
4 xdel ( winsid () ) ;
5 A =[0 6 -5;1 0 2;3 2 4]
6 B = spec ( A )
7 // E i g e n v e c t o r s c o r r e s p o n d i n g t o e i g e n values of A
are
8 p1 =[2; -1; -2];
9 p2 =[1; -0.4285; -0.7142];
10 p3 =[1; -0.4489; -0.93877];
11 T =[ p1 p2 p3 ];
12 A1 = inv ( T ) * A * T
1 // Example s e c 4 . 1 0
2 // C o n t r o l l a b l e companion form
3 clear ; clc ;
4 xdel ( winsid () ) ;
23
5 A =[1 0 0;0 2 0;0 0 3]
6 B =[1 0;0 1;1 1]
7 b1 =[1;0;1]
8 b2 =[0;1;1]
9 u =[ B A * B A ^2 B ]
10 u1 =[1 0 1;0 1 0;1 1 3]
11 // u1 i s a r r a n g e d from [ b1 Aˆ ( v1 −1) ∗ b1 Aˆ ( v2 −1) ∗ b2 ]
12 // v1 and v2 a r e c o n t r o l l a b i l i t y i n d i c e s .
13 u1 =[ b1 A * b1 b2 ]
14 v1 =2;
15 v2 =1;
16 inv ( u1 )
17
18 p1 =[ -0.5 -0.5 0.5]
19 p2 =[0 1 0]
20 P =[ p1 ; p1 * A ; p2 ]
21 A1 = P * A * inv ( P )
22 B1 = P * B
23 C = eye (3 ,3)
24 s = %s
25 D = s *C - A1
26 E = det ( D )
27 routh_t ( E )
28 // t o g e t e q u a t i o n E , A must be e q u a l t o
29 A2 =[0 1 0;0 0 1; -6 -11 -6]
30 B2 =[0 0;1 0.5;0 1]
31 N1 =[6 1.5 4.5; -6 -11 8]
32 N = N1 * P
24
Chapter 5
1 // Example 5 . 1
2 // H u r w i t z s t a b i l i t y t e s t i n s c i l a b
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 s = %s
7 A = s ^4+8* s ^3+18* s ^2+16* s +4 // c h a r a c t e r i s t i c e q u a t i o n
8
9 // c o e f f i c i e n t s o f c h a r a c t e r i s t i c e q u a t i o n
10 a0 = det ( coeff (A ,4) )
11 a1 = det ( coeff (A ,3) )
12 a2 = det ( coeff (A ,2) )
13 a3 = det ( coeff (A ,1) )
14 a4 = det ( coeff (A ,0) )
15
16 D =[ a1 a0 0 0; a3 a2 a1 a0 ;0 a4 a3 a2 ;0 0 0 a4 ] //
Hurwitz determinant
17
18 // m i n o r s o f h u r w i t z d e t e r m i n a n t
25
19 D1 =[ a1 ]
20 det ( D1 )
21 D2 =[ a1 a0 ; a3 a2 ]
22 det ( D2 )
23 D3 =[ a1 a0 0; a3 a2 a1 ;0 a4 a3 ]
24 det ( D3 )
25 D4 =[ a1 a0 0 0; a3 a2 a1 a0 ;0 a4 a3 a2 ;0 0 0 a4 ]
26 det ( D2 )
1 // Example 5 . 2
2 // c o n s t r u c t i n g Routh a r r a y i n s c i l a b
3 clear ; clc ;
4 xdel ( winsid () ) ;
5 mode (0) ;
6
7 s = %s ;
8
9 A = s ^4+4* s ^3+4* s ^2+3* s ; // c h a r a c t e r i s t i c e q u a t i o n
10
11 k = poly (0 , ’ k ’ )
12
13 routh_t ((1) /A , poly (0 , ’ k ’ ) )
14 disp ( ”0<k < 2 . 4 3 7 5 ” )
15
16 // t h e f u n c t i o n w i l l a u t o m a t i c a l l y c o m p u t e s Routh
array
17 // from t h e Routh a r r a y t h e v a l u e o f ” k ” l i e s b e t w e e n
0 and 2 . 4 3 7 5
26
Scilab code Exa 5.2.2a Routh array
1 // Example s e c 5 . 2 . 2 a
2 // Routh a r r a y i n s c i l a b
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 s = poly (0 , ’ s ’ )
7 A = s ^5+ s ^4+2* s ^3+2* s ^2+4* s +6
8 routh_t ( A )
1 // Example s e c 5 . 2 . 2 b
2 // Routh a r r a y i n s c i l a b
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 s = poly (0 , ’ s ’ )
7 B = s ^5+2* s ^4+6* s ^3+12* s ^2+8* s +16
8 routh_t ( B )
9 // I n t h i s e x a m p l e a row o f z e r o f o r m s a t s ˆ 3 .
10 // The f u n c t i o n a u t o m a t i c a l l y t h e d e r i v a t i v e o f t h e
11 // a u x i l l a r y p o l y n o m i a l 2∗ s ˆ4+12∗ s ˆ2+16
12 // v i z =8∗ s ˆ3+24∗ s
27
Scilab code Exa 5.2.2c Routh array
1 // Example s e c 5 . 2 . 2 c
2 // Routh a r r a y i n s c i l a b
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 s = poly (0 , ’ s ’ )
7 p = poly (0 , ’ p ’ )
8 C = s ^5+ s ^4+2* s ^3+2* s ^2+3* s +5
9
10 // s u b s t i t u t i n g ” s =(1/ p ) ” i n B
11 // The r e s u l t i n g c h a r a c t e r i s t i c e q u a t i o n i s
12
13 C1 =5* p ^5+3* p ^4+2* p ^3+2* p ^2+ p +1
14 routh_t ( C1 )
1 // Example s e c 5 . 2 . 2 d
2 // Routh a r r a y i n s c i l a b
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 s = poly (0 , ’ s ’ )
7 D =2* s ^6+2* s ^5+3* s ^4+3* s ^3+2* s ^2+ s +1
8 routh_t ( D )
9
28
10 D1 = s ^2+1
11 // d i v i d i n g t h e main p o l y n o m i a l D by t h e auxillary
p o l y n o m i a l D1
12 D / D1
13 D2 =2* s ^4+2* s ^3+ s ^2+ s +1
14 routh_t ( D2 )
1 // Example 5 . 3
2 // C o n s t r u c t i n g Routh a r r a y i n s c i l a b
3
4 clear ; clc
5 xdel ( winsid () ) ;
6 mode (0) ;
7
8 s = %s ;
9
10 A = s ^4+4* s ^3+4* s ^2+3* s ; // c h a r a c t e r i s t i c equation
after simplification
11
12 k = poly (0 , ’ k ’ )
13
14 routh_t ((1) /A , poly (0 , ’ k ’ ) )
15
16 // s y s t e m w i l l c o n s t r u c t Routh a r r a y and
17 // from Routh a r r a y ” k ” must l i e b e t w e e n 0&39/16 i . e
(0<k < 2 . 4 3 7 5 )
18
19 disp ( ”0<k <39/16 ” )
29
Scilab code Exa 5.4 Routh array
1 // e x a m p l e 5 . 4
2 // c o n s t r u c t i n g Routh a r r a y i n s c i l a b
3 clear ; clc
4 xdel ( winsid () ) ; // c l o s e a l l windows
5 mode (0) ;
6 s = %s ;
7 A = s ^3+8* s ^2+26* s +40;
8
9 // c o n s i d e r p−p l a n e i s l o c a t e d t o t h e l e f t o f t h e s−
plane .
10 // d i s t a n c e b e t w e e n p−p l a n e and s−p l a n e i s 1 .
11 // i f t h e o r i g i n i s s h i f t e d from s−p l a n e t o t h e p−
p l a n e , then , s=p−1
12
13 z = %z
14 B = z ^3+5* z ^2+13* z +21; // s u b s t i t u t i n g s=p−1 i n t h e
e q u a t i o n o f A, t h e r e s u l t i n g e q u a t i o n w i l l be
15 routh_t ( B )
1 // Example 5 . 5 a
2 // c o n s t r u c t i n g Routh a r r a y i n s c i l a b
3 clear ; clc
4 xdel ( winsid () ) ; // c l o s e a l l windows
5 mode (0) ;
30
6 s = %s ;
7 A = s ^3+ s ^2 - s +1
8 routh_t ( A )
9
10 // Example 5 . 5 b
11
12 s = %s ;
13 B = s ^4 - s ^2 -2* s +2
14 routh_t ( B )
15 // i n t h i s e x a m p l e 0 o c c u r s i n t h e f i r s t column o f
the array
16 // f o r which s y s t e m a s s u m e s any s m a l l v a l u e ” e p s ”
and c o m p u t e s t h e a r r a y a u t o m a t i c a l l y .
1 // Example 5 . 6
2 // c o n s t r u c t i n g Routh a r r a y i n s c i l a b
3 clear ; clc ;
4 xdel ( winsid () ) ;
5 mode (0) ;
6
7 s = %s ;
8
9 A = s ^4+8* s ^3+24* s ^2+32* s ; // c h a r a c t e r i s t i c e q u a t i o n
10 k = poly (0 , ’ k ’ )
11
12 routh_t ((1) /A , poly (0 , ’ k ’ ) )
13 disp ( k =80)
14
15 // s i n c e from t h e f o u r t h row o f t h e Routh a r r a y
16 // t h e p o s i t i v e v a l u e ” k=80” w i l l g i v e r o o t s w i t h
zero r e a l part .
31
32
Chapter 6
1 // Example 6 . 2
2 // P l o t t i n g r o o t l o c u s
3 clear ; clc ;
4 xdel ( winsid () ) ;
5 s = %s ;
6 num =1;
7 den = s *( s +3) ^2;
8 G = syslin ( ’ c ’ , num / den ) ;
9 clf () ;
10 evans ( G ) ;
11 axes_handle . grid =[1 1]
12 mtlb_axis ([ -5 5 -5 5]) ;
13 // form t h e g r a p h i t can be s e e n t h a t t h e b r e a k away
p o i n t i s a t ”−1”
14 disp ( ” Break away p o i n t =−1” )
33
Scilab code Exa 6.2.2 location of the root locus between poles and zeros
1 // Example s e c 6 . 2 . 2
2 // l o c a t i o n o f r o o t l o c u s i n b e t w e e n p o l e s and z e r o s
.
3 clear ; clc ;
4 xdel ( winsid () ) ;
5 s = %s ;
6 num =(( s +1) *( s +2) ) ;
7 den =( s *( s +3) *( s +4) ) ;
8 G = syslin ( ’ c ’ , num / den ) ;
9 clf () ;
10 evans ( G ) ;
11 axes_handle . grid =[1 1]
12 mtlb_axis ([ -5 5 -5 5]) ;
1 // Example 6 . 3
2 // P l o t t i n g r o o t l o c u s
3 clear ; clc ;
4 xdel ( winsid () ) ;
5 s = %s ;
6 num =( s +2) ;
34
7 den1 =( s +1+( %i * sqrt (3) ) ) *( s +1+( %i * sqrt (3) ) ) ;
8 // upon s i m p l i f i c a t i o n t h e d e n o m i n a t o r becomes
9 den2 =( s ^2+2* s +4)
10 G = syslin ( ’ c ’ , num / den2 ) ;
11 clf () ;
12 evans ( G ) ;
13 axes_handle . grid =[1 1]
14 mtlb_axis ([ -5 5 -5 5]) ;
1 // Example 6 . 4
2 // P l o t t i n g r o o t l o c u s
3 clear ; clc ;
4 xdel ( winsid () ) ;
5 s = %s ;
6 num = -( s +2) ;
7 den1 =( s +1+( %i * sqrt (3) ) ) *( s +1+( %i * sqrt (3) ) ) ;
8 // upon s i m p l i f i c a t i o n t h e d e n o m i n a t o r becomes
9 den2 =( s ^2+2* s +4)
10 G = syslin ( ’ c ’ , num / den2 ) ;
11 clf () ;
12 evans ( G ) ;
13 axes_handle . grid =[1 1];
14 mtlb_axis ([ -3 3 -3 3]) ;
35
Scilab code Exa 6.5 Root locus
1 // Example 6 . 5
2 // P l o t t i n g r o o t l o c u s
3 clear ; clc ;
4 xdel ( winsid () ) ;
5 s = %s ;
6 num =1;
7 den = s *( s +4) *( s ^2+4* s +20) ;
8 G = syslin ( ’ c ’ , num / den ) ;
9 clf ;
10 evans ( G ) ;
11 axes_handle . grid =[1 1]
12 mtlb_axis ([ -5 5 -5 5]) ;
1 // Example 6 . 6
2 // P l o t t i n g r o o t l o c i in s c i l a b
3 clear ; clc ;
4 xdel ( winsid () ) ;
5 s = %s ;
6 num =( s +2) ;
7 den =( s +1) ^2;
36
8 t = syslin ( ’ c ’ , num / den ) ;
9 clf ;
10 evans ( t ) ;
11 axes_handle . grid =[1 1]
12 mtlb_axis ([ -4 4 -4 4]) ;
1 // Example 6 . 7
2 // P l o t t i n g r o o t l o c u s
3 clear ; clc ;
4 xdel ( winsid () ) ;
5 Beta =0
6 s = %s ;
7 num =1;
8 den = s *( s +1) *( s + Beta ) ;
9 G = syslin ( ’ c ’ , num / den ) ;
10 clf () ;
11 evans ( G ) ;
12 axes_handle . grid =[1 1]
13 mtlb_axis ([ -4 4 -4 4]) ;
37
Chapter 7
1 // Example 7 . 3 . 1 a
2 // Bode p l o t i n s c i l a b
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 s = poly (0 , ’ s ’ ) ;
7 H = syslin ( ’ c ’ ,(10*(1+ s ) ) ,s ^2*(1+.25* s +0.0625* s ^2) ) ;
8 clf () ;
9 bode (H ,0.1 ,1000)
1 // Example : ( i ) 7 . 3 . 1 b
2 // Bode p l o t i n s c i l a b
3 clear ; clc ;
38
4 xdel ( winsid () ) ;
5
6 s = poly (0 , ’ s ’ ) ;
7 G = syslin ( ’ c ’ ,(8*(1+0.5* s ) ) ,s *(1+2* s ) *(1+0.05* s
+0.0625* s ^2) ) ;
8 clf () ;
9 bode (G ,0.01 ,1000) ;
39
Chapter 8
1 // Example 8 . 1
2 // N y q u i s t p l o t
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 s = %s /2/ %pi ;
7 num =(1) ;
8 den = s *( s +1) ;
9 G = syslin ( ’ c ’ ,num , den )
10 clf () ;
11 nyquist ( G )
40
1 // Example 8 . 2
2 // N y q u i s t p l o t
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 s = %s /2/ %pi ;
7 // s i n c e t h e v a l u e o f ”K” and ” t a u ” i n t h e g i v e n
transfer function i s constant
8 // t h u s a s s u m i n g ”K=1” and ” t a u =1”
9 // t h e r e s u l t i n g t r a n s f e r f u n c t i o n i s ,
10 num2 =(1) ;
11 den2 =( s +1) ^2;
12 G = syslin ( ’ c ’ , num2 , den2 )
13 clf () ;
14 nyquist ( G )
1 // Example 8 . 3
2 // N y q u i s t p l o t
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 s = %s /2 / %pi ;
7 num =( s +3) ;
8 den =( s +1) *( s -1)
9 G = syslin ( ’ c ’ ,num , den )
10 clf () ;
11 nyquist ( G )
41
42
Chapter 9
1 // Example s e c 9 . 1
2 // c o m p e n s a t i o n i n open l o o p c o n t r o l s y s t e m
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 s = %s ;
7 disp ( ”G=(60∗ k ) / s ∗ ( s +1) ∗ ( s +6) ” )
8 // v e l o c i t y e r r o r c o n s t a n t ”Kv” when u n i t ramp i n p u t
i s a p p l i e d t o G i s ”5 k ” .
9 // I f ” k =1” , then , s t e a d y s t a t e e r r o r i s 0 . 2 .
10 // when ” k =35/60” G becomes
11
12 num =35;
13 den = s *( s +2) *( s +6) ;
14 G1 = syslin ( ’ c ’ ,num , den ) ;
15 subplot (1 ,2 ,1) ;
16 evans ( G1 )
17 // From t h e f i g u r e 9 . 1
18 OA = sqrt ((0.3) ^2+(2.8) ^2) ;
19 wn1 = OA
20 theta =84 // a n a l y t i c a l l y c a l c u l a t e d
43
21 zeta1 = cosd ( theta )
22 Ts1 =4/( zeta1 * wn1 ) // Ts1= s e t t l i n g time in seconds
23 // For z e t a t o be 0 . 6 and s e t t l i n g time l e s s than 0 . 4
sec
24 a = acosd (0.6)
25 //By d r a w i n g a n g l e ” a ” on t h e r o o t l o c u s
26 OB =1.26;
27 wn2 = OB ;
28 Ts2 =4/(0.6*1.26) // i n s e c o n d s
29 k =10.5/60;
30 // s u b s t i t u t i n g ” s =0” and ”60 k = 1 0 . 5 ” i n t h e e q u a t i o n
f o r G.
31 Kv1 =10.5/12 //Kv= v e l o c i t y e r r o r c o e f f i c i e n t
32 Ess1 = 1/ Kv1 // E s s= s t e a d y s t a t e e r r o r
33 //To g e t t h e r e q u i r e d v a l u e o f t h e z e t a , s t e a d y
s t a t e e r r o r i n c r e a s e s and s e t t l i n g t i m e i m p r o v e s .
34
35 // i n s e r t i n g one z e r o i n t h e e x p r e s s i o n f o r ”G”
36 disp ( ”G2=60∗k ∗ ( s +3) / s ∗ ( s +2) ∗ ( s +6) ” )
37 // c o n s i d e r i n g k=1
38 num1 =60*( s +3) ;
39 den1 =( s *( s +1) *( s +6) ) ;
40 G3 = syslin ( ’ c ’ , num1 , den1 ) ;
41 subplot (1 ,2 ,2) ;
42 evans ( G3 ) ;
43 // c o n s i d e r i n g ” z e t a =0.6” and d r a w i n g l i n e OA a t an
a n g l e 5 3 . 1 3 , on t h e r o o t l o c u s .
44 zeta =0.6;
45 OA1 =3.4;
46 wn = OA1
47 K =16/60
48 Kv =(60* K *3) /(2*6)
49 Ts =4/( zeta * wn ) // i n s e c o n d s
50 Ess =1/ Kv
44
45
Chapter 10
1 // Example s e c 1 0 . 1 . 1
2 // mass , d a s h p o t , s p r i n g a r r a n g e m e n t .
3 clear ; clc ;
4 xdel ( winsid () ) ;
5 M =1
6 K =2
7 F =2
8 A =[0 1; -2 -2]
9 C = eye ( A )
10 s = %s
11 D = s *C - A
12 X = inv ( D ) *[1;1]
13 // t a k i n g t h e l a p l a c e t r a n s f o r m o f X
14 disp ( ”X( t )=s q r t ( 5 ) ∗ s i n ( t+i n v ( t a n 0 . 5 ) ) ; s q r t ( 1 0 ) ∗ s i n (
t+i n v ( t a n −1/3) ) ” )
15 disp ( ” The s y s t e m i s a s y m p t o t i c a l l y s t a b l e ” )
46
Scilab code Exa 10.3 determination of quadratic form
1 // Example 1 0 . 3
2 // d e t e r m i n a t i o n o f q u a d r a t i c form
3 clear ; clc ;
4 xdel ( winsid () ) ;
5 // from t h e g i v e n e u a t i o n we g e t t h e f o l l o w i n g
6 A =[9 1 -2;1 4 -1; -2 -1 1]
7 det ( A )
8 A1 =[9 1;1 4]
9 det ( A1 )
10 // s i n c e d e t e r m i n a n t o f A and A1 i s p o s i t i v e
11 // t h e r e f o r e W i s p o s i t i v e d e f i n i t e .
12 disp ( ”W i s p o s i t i v e d e f i n i t e ” )
1 // Example 1 0 . 4
2 // Lipunov ’ s method
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 x1 = poly (0 , ’ x1 ’ ) ;
7 x2 = poly (0 , ’ x2 ’ ) ;
8 x11 = poly (0 , ’ x11 ’ ) ;
9 x22 = poly (0 , ’ x22 ’ ) ;
10 x2 = x11
11 // a s s u m i n g K1 and K2 e q u a l t o one .
12 disp ( ”W=x1 ˆ2+ x2 ˆ2 ” )
13 // ”W=x1 ˆ2+ x2 ˆ2” i s Liapunov ’ s f u n c t i o n
14 //W i s c h o s e n a r b i t r a r i l y , s i n c e t h e r e no s t a n d a r d
p r o c e d u r e f o r s e l e c t i n g W.
15 disp ( ”dW/ d t =2∗ x1 ∗ x11 +2∗ x2 ∗ x22 =−2∗( x2 ˆ2+ x2 ˆ 4 ) ” )
47
16 disp ( ” T h i s w i l l be n e g a t i v e s e m i d e f i n i t e and
t h e r e f o r e t h e s y s t e m w i l l be s t a b l e ” )
48
Chapter 11
1 // Example 1 1 . 6
2 // Jury ’ s s t a b i l i t y t e s t
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 z = %z ;
7 F =4* z ^4+6* z ^3+12* z ^2+5* z +1
8 // e q u a t i n g t h e e q u a t i o n F w i t h a4 ∗ z ˆ4+ a3 ∗ z ˆ3+ a2 ∗ z ˆ2+
a1 ∗ z 3+a0 .
9 a0 =1
10 a1 =5
11 a2 =12
12 a3 =6
13 a4 =4
14
15 b0 =[ a0 a4 ; a4 a0 ]
16 det ( b0 )
17 b1 =[ a0 a3 ; a4 a1 ]
18 det ( b1 )
19 b2 =[ a0 a2 ; a4 a2 ]
20 det ( b2 )
49
21 b3 =[ a0 a1 ; a4 a3 ]
22 det ( b3 )
23
24 c0 =[ det ( b0 ) det ( b3 ) ; det ( b3 ) det ( b0 ) ]
25 det ( c0 )
26 c1 =[ det ( b0 ) det ( b2 ) ; det ( b3 ) det ( b1 ) ]
27 det ( c1 )
28 c2 =[ det ( b0 ) det ( b1 ) ; det ( b3 ) det ( b2 ) ]
29 det ( c2 )
30
31 disp ( ” d e t ( a0 )<d e t ( a4 )= s a t i s f i e d ” )
32 disp ( ” d e t ( b0 )>d e t ( b3 )= s a t i s f i e d ” )
33 disp ( ” d e t ( c 0 )<d e t ( c 3 )=n o t s a t i s f i e d ” )
34
35 disp ( ” The s y s t e m i s u n s t a b l e ” )
1 // Example s e c 1 1 . 9 . 2 a
2 // s t a b i l i t y o f l i n e a r c o n t i n u o u s s y s t e m
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 s = %s ;
7 G =1/( s *( s +1) *( s +2) )
8 G1 = pfss ( G )
9 // t a k i n g Z t r a n s f o r m o f G1
10 z = %z ;
11 G2 =( z /(2*( z +1) ) ) -( z /( z + %e ^( -1) ) ) +( z /(2*( z + %e ^( -2) ) ) )
12 // upon s i m p l i f i c a t i o n we g e t t h e f o l l o w i n g
c h a r a c t e r i s t i c equation
13 B = z ^3 -(1.3* z ^2) +0.85* z -0.5
14 // s u b s t i t u t i n g ” z =(1+ r /1− r ) ” i n B
50
15 // t h e r e s u l t a n t e q u a t i o n i s B1
16 r = poly (0 , ’ r ’ ) ;
17 B1 =3.65* r ^3+1.95* r ^2+2.35* r +0.05
18 routh_t ( B1 )
19 disp ( ” The s y s t e m i s s t a b l e ” )
1 // Example s e c 1 1 . 9 . 2 b
2 // s t a b i l i t y o f l i n e a r c o n t i n u o u s s y s t e m
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 s = %s ;
7 G =5/( s *( s +1) *( s +2) )
8 G1 = pfss ( G )
9 // t a k i n g Z t r a n s f o r m o f G1
10 z = %z ;
11 G2 =5*(( z /(2*( z +1) ) ) -( z /( z + %e ^( -1) ) ) +( z /(2*( z + %e ^( -2)
))))
12 // upon s i m p l i f i c a t i o n we g e t t h e f o l l o w i n g
c h a r a c t e r i s t i c equation
13 B = z ^3 -(0.5* z ^2) +2.49* z -0.496
14 // s u b s t i t u t i n g ” z =(1+ r /1− r ) ” i n B
15 // t h e r e s u l t a n t e q u a t i o n i s B1
16 r = poly (0 , ’ r ’ )
17 B1 =3.5* r ^3 -2.5* r ^2+0.5* r +2.5
18 routh_t ( B1 )
19 disp ( ” The s y s t e m i s u n s t a b l e ” )
51
Scilab code Exa 11.9.3 Schurcohn stability test
1 // Example s e c 1 1 . 9 . 3
2 // S c h u r c o h n s t a b i l i t y t e s t
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 z = %z
7 G =1/(1 -((7/4) *( z ^ -1) ) -((1/2) *( z ^ -2) ) )
8 A2 =1 -((7/4) *( z ^ -1) ) -((1/2) *( z ^ -2) )
9 //K2= c o e f f i c i e n t o f z ˆ−2
10 K2 = -0.5
11 B2 = -0.5 -1.75*( z ^ -1) + z ^ -2
12
13 A1 =( A2 - K2 * B2 ) /(1 - K2 ^2)
14 //K1= c o e f f i c i e n t o f z ˆ−1
15 K1 = -3.5
16 //mod ( K1 ) >1 and mod ( K2 ) <1
17 disp ( ” The sytem i s u n s t a b l e ” )
52
Chapter 15
Scilab code Exa 15.2 Time domain specifications of second order system
1 // Example 1 5 . 2
2 // t i m e domain s p e c i f i c a t i o n s o f s e c o n d o r d e r s y s t e m
3 clear ; clc ;
4 xdel ( winsid () ) ;
5 mode (0) ;
6
7 // c o n v e r t i n g t h e g i v e n d i f f e r e n t i a l e q u a t i o n i n ” s ”
domain
8 // s i n c e x and y a r e c o n s t a n t s
9 // t h e r e f o e r e c o n s i d e r i n g ” x=y=1”
10
11 s = %s ;
12 g = s ^2+2* s ;
13 x = roots ( g )
14 wn = sqrt ( abs ( x (1) ) ) // undamped n a t u r a l f r e q u e n c y
15 zeta =(1/ wn ) // damping r a t i o
16 wd = wn * sqrt (1 - zeta ^2) // damped n a t u r a l f r e q u e n c y
17 Dc =( zeta * wn ) // Dc=damping c o e f f i c i e n t
18 Tc =1/( zeta * wn ) // Tc=t i m e c o n s t a n t o f t h e s y s t e m
53
Scilab code Exa 15.4 transfer function of gyroscope
1 // Example 1 5 . 4
2 // T r a n s f e r f u n c t i o n o f G y r o s c o p e
3 clear ; clc ;
4 xdel ( winsid () ) ;
5 // i n c a s e o f G y r o s c o p e t h e e q u a t i o n i s
6
7 disp ( ” ( J ∗ s ˆ2+B∗ s+K) t h e t a ( s )=H∗w( s ) ” )
8 // t h e r e f o r e
9 disp ( ” t h e t a ( s ) /w( s )=H/ J ∗ s ˆ2+B∗ s+K” )
1 // Example 1 5 . 5 ( f i g 1 5 . 4 )
2 // t r a n s f e r f u n c t i o n o f t h e s y s t e m
3 clear ; clc ;
4 xdel ( winsid () ) ;
5 mode (0) ;
6
7 s = poly (0 , ’ s ’ ) ;
8 //G1 and G2 a r e c o n n e c t e d i n s e r i e s
9 G1 = s ^2/( s +4) ^2
10 G2 =( s +1) /( s ^3*( s +3) )
11 //H1 i s f e e d b a c k l o o p
12 H1 =( s ^2+ s +1) /( s *( s +3) )
13 // Tf= t r a n s f e r f u n c t i o n
14 Tf =( G1 * G2 * H1 )
54
15 A = type ( s ) ;
16 disp (A , ’ Type o f t h e s y s t e m= ’ )
1 // Example 1 5 . 6
2 // c o m p a r i s o n o f s e n s i v i t y o f t h e two s y s t e m
3 clear ; clc ;
4 xdel ( winsid () ) ;
5 // k1 andk2 a r e s e r i e s b l o c k s o f t h e t r a n s f e r
function
6 k1 =100
7 k2 =100
8 // t r a n s f e r f u n c t i o n o f f i g . 1 5 . 5
9 T1 = k1 * k2 /(1+(0.0099* k1 * k2 ) )
10 // t r a n s f e r f u n c t i o n o f f i g . 1 5 . 6
11 T2 =( k1 /(1+(0.09* k1 ) ) ) *( k2 /(1+(0.09* k2 ) ) )
12 disp ( ” b o t h t r a n s f e r f u n c t i o n a r e e q u a l ” )
13 // s e n s i t i v i t y o f t h e t r a n s f e r f u n c t i o n T1 w i t h
r e s p e c t t o k1
14 T11 =1/(1+(0.0099* k1 * k2 ) )
15 // s e n s i t i v i t y o f t h e t r a n s f e r f u n c t i o n T2 w i t h
r e s p e c t t o k1
16 T12 =1/(1+(0.09* k1 ) )
17 disp ( ” The s y s t e m o f f i g 1 5 . 6 i s 10 t i m e s more
s e n s i t i v e than system o f f i g 1 5 . 5 with r e s p e c t to
v a r i a t i o n s i n k1 ” )
55
1 // Example 1 5 . 7
2 // f i n d bandwidth o f t h e t r a n s f e r f u n c t i o n
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 s = %s ;
7 O =1;
8 R =( s +1) ;
9 tf = O / R
10 disp ( ” when O/R( jw ) = 0 . 7 0 7 , w=wc” )
11
12 wc =(1/0.707) ^2 -1
13 // wc=bandwidth o f t h e t r a n s f e r f u n c t i o n
14
15 disp ( ” Hence t h e bandwidth i s 1 r a d / s e c ” )
1 // Example 1 5 . 8
2 // f i n d bandwidth o f t h e t r a n s f e r f u n c t i o n
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 s = %s ;
7 O =6;
8 R =( s ^2+2* s +6) ;
9 tf = O / R
10
11 disp ( ” when O/R( jw ) =6/ s q r t (wˆ4−8∗w+36) ” )
12
13 w =[+2 -2] // a f t e r d i f f e r e n t i a t i o n and s i m p l i f i c a t i o n
14
15 disp ( ” when O/R( jw ) =6/ s q r t (wˆ4−8∗w+36) , At w=+−2” )
56
16
17 peak =3/ sqrt (5)
1 // Example 1 5 . 9
2 // N y q u i s t p l o t
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 s = %s /2/ %pi ;
7 num =(1) ;
8 den = s ^3*( s +1) ;
9 G = syslin ( ’ c ’ ,num , den )
10 clf ()
11 nyquist ( G )
1 // Example 1 5 . 1 2
2 // p r o v e t h e s o l u t i o n of the equation
3 clear ; clc ;
4 xdel ( winsid () ) ;
5 // a s s u m i n g n=1
6 n =1;
7 z = %z ;
8 y(n)=z^n;
57
9 y ( n +1) = z ^( n +1) ;
10 y ( n +2) = z ^( n +2) ;
11 A = y ( n +2) +3* y ( n +1) +2* y ( n )
12 B=A/z
13 roots ( z ^2+3* z +2)
14 disp ( ” y ( n )=z ˆn i s s o l u t i o n o f p o l y n o m i a l e q u a t i o n ( z
+2) ∗ ( z +1)=0” )
1 // Example 1 5 . 1 3
2 // Time domain s p e c i f i c a t i o n s
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 J =5.5*10^ -2;
7 f =3.0*10^ -4;
8 disp ( ”wn=s q r t ( k / J ) =10ˆ3∗ s q r t ( k / 5 . 5 ) ” )
9 disp ( ” z e t a=s q r t ( 4 . 9 ∗ 1 0 ˆ − 3 / k ) ” )
10 // a t c r i t i c a l l y damped c o n d i t i o n ” z e t a =1” , t h e r e f o r e
11 k =4.09*10^ -3
12 // when k =1.5∗10ˆ −2
13 zeta = sqrt ((4.09*10^ -3) /(1.5*10^ -2) )
14 wn =10^3* sqrt (1.5*10^ -2/5.5)
15 wd =( wn /(2* %pi ) ) * sqrt (1 - zeta ^2)
16 //wd=f r e q u e n c y o f damped o s c i l l a t i o n
17 Pwd =1/ wd
18 //Pwd=p e r i o d o f damped o s c i l l a t i o n
58
Scilab code Exa 15.14 Position servomotor
1 // Example 1 5 . 1 4
2 // p o s i t i o n s e r v o m o t e r
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 // M i l= motor i n e r t i a r e f e r r e d t o t h e l o a d s i d e
7 Mil =20^2*0.45*10^ -6 // u n i t= kg .mˆ2
8
9 // Tr= T r a n s f o r m a t i o n r a t i o o f g e a r t r a i n b e t w e e n t h e
l o a d s h a f t and t h e t a c h o g e n e r a t o r
10 Tr =20*2
11
12 // t i l = t a c h o g e n e r a t o r i n e r t i a r e f e r r e d to the load
side
13 til =40^2*0.35*10^ -6 // u n i t= kg .mˆ2
14
15 // T i l= t o t a l i n e r t i a r e f e r r e d t o t h e l o a d s i d e
16 Til =(20*10^ -6) +(1.8*10^ -4) +(5.6*10^ -4) // u n i t= kg .m
ˆ2
17
18 // Mi= i n e r t i a r e f e r r e d t o t h e motor s i d e
19 Mi =(760*10^ -6) /400 // u n i t= kg .mˆ2
1 // Example 1 5 . 1 5
2 // s t e a d y o u t p u t s p e e d o f DC motor
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 //Jm= moment o f i n e r t i a o f motor
59
7 Jm =6.5*10^ -2;
8 //Fm= f r i c t i o n o f motor
9 Fm =3.5*10^ -3;
10 // a=g e a r r a t i o
11 a =1/100;
12 // J l= i n e r t i a o f l o a d
13 Jl =420;
14 // F l= f r i c t i o n o f l o a d
15 Fl =220;
16 // J= t o t a l moment o f i n e r t i a
17 J = Jm +( a ^2* Jl ) // u n i t=kg .mˆ2
18 //F= t o t a l f r i c t i o n
19 F = Fm +( a ^2* Fl ) // u n i t=kg .mˆ2
20 s = %s
21 //wm1=A n g u l a r v e l o c i t y i n f r e q u e n c y domain
22 wm1 =2/( s *(( J * s ) + F ) )
23 t =1;
24 //wm2=A n g u l a r v e l o c i t y i n t i m e domain
25 // s i n c e ” t =1” , wm2 i s i n i t i a l v a l u e o f a n g u l a r
velocity
26 wm2 =(2/ F ) *(1 -( %e ^(( -5.7*10^ -2) /(10.7*10^ -2) ) * t ) ) //
u n i t=r a d / s e c
27 //Nm1=motor s p e e d i n r p s ( i n i t i a l s p e e d )
28 Nm1 = wm2 /(2* %pi ) ;
29 //Nm2=motor s p e e d i n rpm
30 Nm2 =( wm2 /(2* %pi ) ) *60; // u n i t=rpm
31 // Nl=l o a d s p e e d
32 Nl =(1/100) *(( wm2 /(2* %pi ) ) *60) // u n i t=rpm
33 // Nos= s t e a d y o u t p u t s p e e d
34 // s i n c e Nos i s s t e a d y s p e e d , t h e e x p o n e n t i a l term o f
wn2 becomes 0 .
35 Nos =(1/100) *(60/(2* %pi ) ) *(2/(5.7*10^ -2) ) // u n i t=rpm
60
Scilab code Exa 15.26 Routh array
1 // Example 1 5 . 2 6
2 // C o n s t r u c t i n g Routh a r r a y i n s c i l a b
3
4 clear ; clc
5 xdel ( winsid () ) ;
6 mode (0) ;
7
8 A =[5 -6 -12; -1 1 2;5 -6 -11]
9 B = eye (3 ,3)
10 s = %s
11 C = s *B - A
12 D = s ^3+5* s ^2+5* s +1; // c h a r a c t e r i s t i c e q u a t i o n a f t e r
simplification
13 routh_t ( D )
14 disp ( ”No s i g n c h a n g e i n t h e f i r s t column , h e n c e t h e
system i s a s y m p t o t i c a l l y s t a b l e ”)
1 // Example 1 5 . 2 7
2 // To c h e c k w h e t h e r t h e s y s t e m i s r e a c h a b l e o r n o t
3
4 clear ; clc
5 xdel ( winsid () ) ;
6 mode (0) ;
7 A =[1 0;0 1]
8 B =[1;1]
9 Wc =[ A * B B ]
10 disp ( ” The r a n k o f Wc=(1∗1 −1∗1) =0 , and n o t e q u a l t o 2 .
Thus t h e g i v e n s y s t e m i s n o t r e a c h a b l e ” )
61
Scilab code Exa 15.28 Determine the stability of the system
1 // Example 1 5 . 2 8
2 // D e t e r m i n e t h e s t a b i l i t y o f t h e s y s t e m .
3
4 clear ; clc
5 xdel ( winsid () ) ;
6 mode (0) ;
7
8
9 z = %z
10
11 D = z ^3+6* z ^2+8* z -0.04; // c h a r a c t e r i s t i c e q u a t i o n
after simplification
12 routh_t ( D )
13 disp ( ” There i s s i g n c h a n g e i n t h e f i r s t column ,
hence the system i s u n s t a b l e ”)
Scilab code Exa 15.31 Time domain specifications of second order system
1 // Example 1 5 . 3 1
2 // t i m e domain s p e c i f i c a t i o n s o f s e c o n d o r d e r s y s t e m
3 clear ; clc ;
4 xdel ( winsid () ) ;
5 mode (0) ;
6
7 // c o n v e r t i n g t h e g i v e n d i f f e r e n t i a l e q u a t i o n i n ” s ”
domain
62
8 // s i n c e x and y a r e c o n s t a n t s
9 // t h e r e f o e r e c o n s i d e r i n g ” x=y=1”
10
11 s = %s ;
12 g = s ^2+5* s +7;
13 x = coeff ( g )
14 // c o m p a r i n g w i t h t h e s t a n d a r d e q u a t i o n o f s e c o n d
order system .
15 wn = sqrt ( x (: ,1) ) // undamped n a t u r a l f r e q u e n c y
16 zeta =(5/(2* wn ) ) // damping r a t i o
17 wd = wn * sqrt (1 - zeta ^2) // damped n a t u r a l f r e q u e n c y
18 Tc =1/( zeta * wn ) // Tc=t i m e c o n s t a n t o f t h e s y s t e m
1 // Example 1 5 . 3 3
2 // Lipunov ’ s method
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 x1 = poly (0 , ’ x1 ’ ) ;
7 x2 = poly (0 , ’ x2 ’ ) ;
8 x11 = poly (0 , ’ x11 ’ ) ;
9 x22 = poly (0 , ’ x22 ’ ) ;
10 x2 = x11
11 disp ( ” x22+x2+x2 ˆ3+ x1=0” )
12 // ( x1 , x2 ) h a s s i n g u l a r p o i n t a t ( 0 , 0 )
13 disp ( ”V=x1 ˆ2+ x2 ˆ2 ” )
14 // ”V=x1 ˆ2+ x2 ˆ2” i s Liapunov ’ s f u n c t i o n
15 //V i s p o s i t i v e f o r a l l v a l u e s o f x1 and x2 , e x c e p t
a t x1=x2=0
16 disp ( ”dV/ d t =2∗ x1 ∗ x2 −2∗x1 ∗ x2 −2∗x2 ˆ2−2∗ x2 ˆ4=−2∗ x2 ˆ2−2∗
x2 ˆ4 ” )
63
17 disp ( ”dV/ d t w i l l n e v e r be p o s i t i v e h e n c e o r i g i n i s
s t a b l e ”)
1 // Example 1 5 . 3 4
2 // f i n d bandwidth o f t h e t r a n s f e r f u n c t i o n
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 s = %s
7 A =1
8 B =( s +1)
9 tf = A / B
10
11 disp ( ” when A/B( jw ) =1/ s q r t ( 2 ) , w=w1” )
12
13 w1 =(1/0.707) ^2 -1
14 //w1=bandwidth o f t h e t r a n s f e r f u n c t i o n
15
16 disp ( ” Hence t h e bandwidth i s 1 r a d / s e c ” )
1 // Example 1 5 . 3 6
2 // i m p u l s e r e s p o n s e t r a n s f e r f u n c t i o n
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
64
6 s = %s ;
7 G = syslin ( ’ c ’ ,25 , s ^2+4* s +25) ;
8 t =0:0.05:5;
9 y = csim ( ’ i m p u l s ’ ,t , G ) ;
10 plot (t , y )
11 xtitle ( ’ I m p u l s e r e s p o n s e 2 5 / ( s ˆ2+4∗ s +25) ’ , ’ t s e c ’ , ’
Response ’ );
1 // Example 1 5 . 3 7
2 // s t e p r e s p o n s e t r a n s f e r f u n c t i o n
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 s = %s ;
7 G = syslin ( ’ c ’ ,25 , s ^2+4* s +25) ;
8 t =0:0.05:5;
9 y = csim ( ’ s t e p ’ ,t , G ) ;
10 plot (t , y )
11 xtitle ( ’ s t e p r e s p o n s e 2 5 / ( s ˆ2+4∗ s +25) ’ , ’ t s e c ’ , ’
Response ’ );
65
1 // Example 1 5 . 3 8
2 // f i n d r o o t s o f c h a r a c t e r i s t i c equation
3
4 clear ; clc ;
5 xdel ( winsid () ) ;
6 s = poly (0 , ’ s ’ )
7 G = s ^4+2* s ^3+ s ^2 -2* s -1
8 roots ( G )
1 // Example : 1 5 . 3 9
2 // Bode p l o t i n s c i l a b
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 s = poly (0 , ’ s ’ ) ;
7 G = syslin ( ’ c ’ ,(25) ,s ^2+4* s +25) ;
8 clf () ;
9 bode (G ,0.01 ,1000) ;
1 // Example 1 5 . 4 0
2 // N y q u i s t p l o t
3 clear ; clc ;
4 xdel ( winsid () ) ;
66
5
6 s = %s /2/ %pi ;
7 num =(1) ;
8 den =( s ^2+0.8* s +1) ;
9 G = syslin ( ’ c ’ ,num , den )
10 clf () ;
11 nyquist ( G )
1 // Example 1 5 . 4 1
2 // N y q u i s t p l o t
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 s = %s /2/ %pi ;
7 num =( s +2) ;
8 den =( s +1) *( s +1) ;
9 G = syslin ( ’ c ’ ,num , den )
10 clf () ;
11 nyquist ( G )
1 // Example : 1 5 . 4 2
67
2 // Bode p l o t i n s c i l a b
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 s = poly (0 , ’ s ’ ) ;
7 G = syslin ( ’ c ’ ,(64*( s +2) ) ,( s *( s +0.5) *( s ^2+3.2* s +64) ) ) ;
8 clf () ;
9 bode (G ,0.01 ,1000) ;
1 // Example : 1 5 . 4 3
2 // e i g e n v a l u e s o f m a t r i x A
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 A =[0 6 -5;1 0 2;3 2 4];
7 spec ( A )
1 // Example 1 5 . 4 4
2 // s t a t e s p a c e r e p r e s e n t a t i o n o f LTI s y s t e m
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 A =[0 1; -2 -3];
68
7 B =[0;1];
8 C =[1 1];
9 D =[0];
10 E =[0];
11
12 H = syslin ( ’ c ’ ,A ,B , C ) ;
13 s = %s ;
14 g = eye (2 ,2) ;
15 P =( - s * g ) -A
16 sm =[ P B ; C D ];
17 H1 = sm2ss ( sm )
1 // Example 1 5 . 4 5
2 // C o v a r i a n t m a t r i x o f ”A”
3 clear ; clc ;
4 xdel ( winsid () ) ;
5 A =[1 0 0;0 2 0;0 0 3]
6 mvvacov ( A )
1 // Example 1 5 . 4 9
2 // P l o t t i n g r o o t l o c i o f t h e t r a n s f e r f u n c t i o n k / s ∗ (
s +4) ∗ ( s ˆ2+4∗ s +20)
3 clear ; clc ;
4 xdel ( winsid () ) ;
5 s = %s ;
69
6 num =(1) ;
7 den = s *( s +3) *( s ^2+2* s +2) ;
8 G = syslin ( ’ c ’ , num / den ) ;
9 clf ;
10 evans ( G ) ;
11 mtlb_axis ([ -5 5 -5 5]) ;
1 // Example : 1 5 . 5 0
2 // Bode p l o t i n s c i l a b
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 s = poly (0 , ’ s ’ ) ;
7 G = syslin ( ’ c ’ ,(16*( s +2) ) ,( s *( s +0.5) *( s ^2+3.2* s +64) ) ) ;
8 clf () ;
9 bode (G ,0.01 ,1000) ;
1 // Example 15 53
2 // s t a t e s p a c e model o f d i f f e r e n t i a l e q u a t i o n .
3 clear ; clc ;
4 xdel ( winsid () ) ;
70
5
6 // c o n v e r t i n g t h e d i f f e r e n t i a l e q u a t i o n i n t e r m s o f
transfer function .
7 s = %s
8 // t r a n s f e r f u n c t i o n
9 A =1/( s ^3+6* s ^2+11* s +6)
10 B = tf2ss ( A )
1 // Example 1 5 . 5 4
2 // N y q u i s t p l o t
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 s = %s /2/ %pi ;
7 num =(2) ;
8 den = s *( s ^2+2* s +2) ;
9 G = syslin ( ’ c ’ ,num , den )
10 clf () ;
11 nyquist ( G )
1 // Example 1 5 . 5 7
2 // d e t e r m i n a t i o n o f z e t a & wn
3 clear ; clc ;
71
4 xdel ( winsid () ) ;
5
6 s = %s
7 num =10;
8 den = s ^2+2* s +10; // s i n c e k=0
9 G = num / den ;
10 B = coeff ( den )
11 //wn= undamped n a t u r a l f t r q u e n c y
12 wn = sqrt ( B (: ,1) )
13 // z e t a= damping r a t i o
14 zeta =2/(2* sqrt ( wn ) )
15 // when t i m e t t e n d s t o i n f i n i t y , s t a t i c e r r o r v i z .
e s s tends to 0 .
16 ess =0
17 // when ” z e t a = 0 . 6 5 ” i . e . ( z e t a 1 = 0 . 6 5 )
18 zeta1 =0.65
19 k0 =2* zeta1 * wn -2
1 // Example 1 5 . 5 9
2 // t r a n s f e r f u n c t i o n o f s i g n a l f l o w g r a p h
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 k1 =1;
7 k2 =5;
8 k3 =5;
9 s = %s ;
10 // From t h e g r a p h t h e t r a n s f e r f u n c t i o n i s
11 T =( k3 * k1 ) /( s ^3+ s ^2+( k3 * k1 ) +( k1 * k2 * s ^2) +5)
12 // s u b s t i t u t i n s ” s =0” i n t h e e q u a t i o n o f T
13 // and d i f f e r e n t i a t i n g and s i m p l i f y i n g t h e e q u a t i o n
72
14 // t h e f o l l o w i n g value of T w i l l appear
15 T1 =1/(1+ k1 )
1 // Example 1 5 . 6 0
2 // r o o t l o c u s
3 clear ; clc ;
4 xdel ( winsid () ) ;
5
6 s = %s ;
7 // s u b s t i t u t i n g ” a =15” i n t h e n u m e r a t o r
8 num =2*( s +15) ;
9 den = s *( s +2) *( s +10) ;
10 G = syslin ( ’ c ’ , num / den ) ;
11 evans ( G ) ;
12 axes_handle . grid =[1 1]
13 mtlb_axis ([ -5 5 -5 5]) ;
73