Scilab Manual for

Control Systems & Electrical Machines

by Dr Lochan Jolly
Electronics Engineering
Thakur College of Engineering & Technology1

Solutions provided by
Dr Lochan Jolly
Electronics Engineering
Mumbai/TCET

June 17, 2019

1 Funded by a grant from the National Mission on Education through ICT,

http://spoken-tutorial.org/NMEICT-Intro. This Scilab Manual and Scilab codes
written in it can be downloaded from the ”Migrated Labs” section at the website
http://scilab.in
1
Contents

List of Scilab Solutions 3

1 To plot pole zero plot of given system and comment on

stability 5

2 To study transient response for given system 7

3 To sketch root locus, bode plot and nyquist plot for given
system 10

4 To compare open loop and closed loop system. 12

5 To use block reduction technique to get the equivalent sys-

tem 14

6 To implement state model reduction technique for the given

system 16

2
List of Experiments

Solution 1.01 polezeroplot . . . . . . . . . . . . . . . . . . . . . 5

Solution 2.02 transientresponse . . . . . . . . . . . . . . . . . . 7
Solution 3.03 stabilityanalysis . . . . . . . . . . . . . . . . . . . 10
Solution 4.04 systemcomparison . . . . . . . . . . . . . . . . . . 12
Solution 5.05 blockreduction . . . . . . . . . . . . . . . . . . . . 14
Solution 6.06 statemodel . . . . . . . . . . . . . . . . . . . . . . 16

3
List of Figures

1.1 polezeroplot . . . . . . . . . . . . . . . . . . . . . . . . . . . 6

2.1 transientresponse . . . . . . . . . . . . . . . . . . . . . . . . 9

3.1 stabilityanalysis . . . . . . . . . . . . . . . . . . . . . . . . . 11

4.1 systemcomparison . . . . . . . . . . . . . . . . . . . . . . . . 13

5.1 blockreduction . . . . . . . . . . . . . . . . . . . . . . . . . . 15

6.1 statemodel . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17

4
 Experiment: 1

To plot pole zero plot of given

system and comment on
stability

Scilab code Solution 1.01 polezeroplot

1 // exp1 f i n d p o l e z e r o p l o t o f g i v e n t r a n s f e r
functions
2 //OS=Windows XP s p 3
3 // S c i l a b v e r s i o n 5 . 4 . 0
4 // s a m p l e v a l u e s
5 s = poly (0 , ’ s ’ ) ;
6 n =[2+3* s +4* s ^2]; // i n p u t n u m e r a t o r o f t r a n s f e r
function
7 d =[1+ s + s ^2]; // i n p u t d e n o m e n a t o r o f t r a n s f e r f u n c t i o n
8 h = syslin ( ’ c ’ ,n / d ) ;
9 plzr ( h ) ;

5
Figure 1.1: polezeroplot

6
 Experiment: 2

To study transient response for

given system

Scilab code Solution 2.02 transientresponse

1 // 2 To s t u d y t r a n s i e n t r e s p o n s e f o r g i v e n
system1
2 //OS=Windows XP s p 3
3 // S c i l a b v e r s i o n 5 . 4 . 0
4 // s a m p l e v a l u e s
5 s = %s
6 // ( ’ E n t e r t h e v a l u e o f damping z f o r u n d e r damped
system ’ )
7 z =0.9
8 num =1; den =36+2* z * s + s ^2; // i n p u t n u m e r a t o r and
denomenator o f t r a n s f e r f u n c t i o n
9 TF = syslin ( ’ c ’ ,num , den ) // t r a n s f e r f u n c t i o n
10 // d i s p (TF)
11 subplot (131)
12 t =0:0.1:10;
13 y1 = csim ( ’ s t e p ’ , t , TF ) ; // t i m e r e s p o n s e
14 title ( ’ Underdamped s y s t e m ’ ) ;
15 plot (t , y1 )
16 // ( ’ E n t e r t h e v a l u e o f damping z f o r c r i t i c a l l y

7
 damped system ’ )
17 z =6
18 num =1; den =36+2* z * s + s ^2;
19 TF = syslin ( ’ c ’ ,num , den )
20 subplot (132)
21 t =0:0.1:10;
22 y1 = csim ( ’ s t e p ’ , t , TF ) ;
23 title ( ’ C r i t i c a l l y damped s y s t e m ’ ) ;
24 plot (t , y1 )
25 // ( ’ E n t e r t h e v a l u e o f damping z f o r o v e r damped
system ’ )
26 z =20
27 num =1; den =36+2* z * s + s ^2;
28 TF = syslin ( ’ c ’ ,num , den )
29 subplot (133)
30 t =0:0.1:10;
31 y1 = csim ( ’ s t e p ’ , t , TF ) ;
32 title ( ’ Overdamped s y s t e m ’ ) ;
33 plot (t , y1 )

8
Figure 2.1: transientresponse

9
 Experiment: 3

To sketch root locus, bode plot

and nyquist plot for given
system

Scilab code Solution 3.03 stabilityanalysis

1 // 3 To s k e t c h r o o t l o c u s , bode p l o t and n y q u i s t
p l o t f o r given system
2 //OS=Windows XP s p 3
3 // S c i l a b v e r s i o n 5 . 4 . 0
4 // s a m p l e v a l u e s
5 s = %s ;
6 T = syslin ( ’ c ’ ,25+30* s +5* s ^2 ,168+206* s +89* s ^2+16* s ^3+ s
^4) ; // t r a n s f e r f u n c t i o n
7 subplot (131)
8 title ( ’ Bode p l o t ’ )
9 bode ( T ) // bode p l o t
10 subplot (132)
11 nyquist ( T ) // n y q u i s t p l o t
12 subplot (133)
13 evans ( T ) ; // r o o t l o c u s

10
Figure 3.1: stabilityanalysis

11
 Experiment: 4

To compare open loop and

closed loop system.

Scilab code Solution 4.04 systemcomparison

1 // 4 To compare open l o o p and c l o s e d l o o p s y s t e m .

2 //OS=Windows XP s p 3
3 // S c i l a b v e r s i o n 5 . 4 . 0
4 // s a m p l e v a l u e s
5 s = %s
6 z =0.1
7 num =1; den =36+2* z * s + s ^2; // i n p u t n u m e r a t o r and
denomenator o f t r a n s f e r f u n c t i o n
8 TF = syslin ( ’ c ’ ,num , den ) // t r a n s f e r f u n c t i o n o f
system
9 disp ( TF )
10 subplot (121)
11 t =0:0.1:30;
12 y1 = csim ( ’ s t e p ’ , t , TF ) ; // t r a n s i s e n t r e s p o n s e
13 title ( ’ Open Loop s y s t e m ’ ) ;
14 plot (t , y1 )
15 s = %s
16 num1 = s ; den1 =1+ s ;
17 TF1 = syslin ( ’ c ’ , num1 , den1 )

12
 Figure 4.1: systemcomparison

18 num2 =1; den2 =3;

19 TF2 = syslin ( ’ c ’ , num2 , den2 )
20 S1 = TF1 /. TF2
21 disp ( S1 )
22 subplot (122)
23 t =0:0.1:30;
24 y1 = csim ( ’ s t e p ’ , t , S1 ) ;
25 title ( ’ C l o s e d Loop s y s t e m ’ ) ;
26 plot (t , y1 )

This code can be downloaded from the website wwww.scilab.in

13
 Experiment: 5

To use block reduction

technique to get the equivalent
system

Scilab code Solution 5.05 blockreduction

1 // 5 To u s e b l o c k r e d u c t i o n t e c h n i q u e t o g e t
the e q u i v a l e n t system
2 // t h i s i s i s a t y p i c a l c a s e where f o r g i v e n b l o c k
d a i g r a m where h3 and h4 form a l o o p which i s i n
s e r i e s t o h1 and h3 i s n e g a t i v e f e e d b a c k t o t h e
combination
3 //OS=Windows XP s p 3
4 // S c i l a b v e r s i o n 5 . 4 . 0
5 // s a m p l e v a l u e s
6 clc
7 s = poly (0 , ’ s ’ ) ;
8 n1 =[2* s ]; // i n p u t n u m e r a t o r and d e n o m e n a t o r o f
transfer function
9 d1 =[3+ s ^2];
10 h1 = syslin ( ’ c ’ , n1 / d1 ) ; // t r a n s i e n t r e s p o n s e
11 n2 =[2];
12 d2 =[ s +4];

14
 Figure 5.1: blockreduction

13 h2 = syslin ( ’ c ’ , n2 / d2 ) ;
14 n3 =[4];
15 d3 =[ s +6];
16 h3 = syslin ( ’ c ’ , n3 / d3 ) ;
17 n4 =[ s ];
18 d4 =[1+2* s ];
19 h4 = syslin ( ’ c ’ , n4 / d4 ) ;
20 h5 = h3 /. h4
21 h6 = h1 * h5
22 h7 = h6 /. h3
23 disp ( h7 , ” Reduced t r a n s f e r f u n c t i o n =” )

15
 Experiment: 6

To implement state model

reduction technique for the
given system

Scilab code Solution 6.06 statemodel

1 // 6 To i m p l e m e n t s t a t e model r e d u c t i o n
t e c h n i q u e f o r the given system
2 //OS=Windows XP s p 3
3 // S c i l a b v e r s i o n 5 . 4 . 0
4 // s a m p l e v a l u e s
5 clc
6 clear all
7 s = poly (0 , ’ s ’ ) ;
8 X1 =[0 1; -6 -5]; //X1 , X2 X3 r e p r e s e n t A B C m a t r i x o f
state space repesentation
9 X2 =[0 ; 1];
10 X3 =[8 1 ]
11 [ n1 , n2 ] = size ( X1 )
12 I = eye ( n1 , n2 ) // i d e n t i t y m a t r i x
13 X = s *I - X1
14 phi = inv ( X ) // i n v e r s e o f m a t r i x
15 Y = X3 * phi

16
 Figure 6.1: statemodel

16 Z = Y * X2
17 // s y s= t f 2 s s ( Z )
18 disp (Z , ” The t r a n f e r function representation of
s y s t e m i s =” )
19 // d i s p ( s y s )

17