Model Answer Control 2 - 2018 - 2019
Model Answer Control 2 - 2018 - 2019
Model Answer Control 2 - 2018 - 2019
( )
( )
Design a lead compensator such that the dominant closed-loop poles are located at
√ . (14 marks)
Answer 2-a):
The angle deficiency at the closed loop poles √ is
Page 1 of 10
Electrical Engineering Department
]txep TpyT[ ECOMP 423: Automatic Control (2)
4th year Computers And Systems Section
Final Examination 1st Term 2018/2019
Date: 02/01/2019 Time: 3.0 Hours
Model Answer (100 Marks) Faculty of Engineering
Aswan University
– First, draw a horizontal line passing through point P, A P
the desired location for one of the dominant closed-
loop poles. This is shown as line PA in following figure. 15 15
– Draw also a line connecting point P and the origin O.
C
– Bisect the angle between the lines PA and PO,
x B D
– Draw two lines PC and PD that make angles Pc Zc O
with the the bisector PB.
– The intersections of PC and PD with the negative real axis give the necessary
locations for the pole and zero of the lead network.
– Also →
– Determine the value of Kc of the lead compensator from the magnitude
( )
condition. ( ) ( ) ( )
( )( )
( )
| |
( )( ) √
– The Kc is calculated as
– Hence, the lead compensator Gc(s) just designed is given by
( )
b) A unity feedback control system has an open loop transfer function G(s)
( )
( )( )
Design a lag compensator to achieve the following specification: Velocity error constant Kv = 8
and phase margin PM = 40o. (16 marks)
Answer 2-b):
We shall use a lag compensator of the form
( ) ̂ = ̂ , ( )
• Define ̂
• Define also ( )= ( ) ( )( )
• The first step in the design is to adjust the gain K to meet the required static
velocity error constant. Thus,
( ) ( ) ( ) ( )
( )( )
Or K= 8
Page 2 of 10
Electrical Engineering Department
]txep TpyT[ ECOMP 423: Automatic Control (2)
4th year Computers And Systems Section
Final Examination 1st Term 2018/2019
Date: 02/01/2019 Time: 3.0 Hours
Model Answer (100 Marks) Faculty of Engineering
Aswan University
• With K=8, the compensated system satisfies the steady-state performance
requirement.
• We shall next plot the Bode diagram of ( ) ( )( )
• The magnitude curve and phase-angle curve of G1(jω) are shown in the
following figure.
Bode Diagram
100
20 dB
50
Magnitude (dB)
-50
-100
-90
-128o
-135
Phase (deg)
-140o
-180
-186.02o
-225
-270
-2 -1 0 1 2
10 10 10 10 10
Frequency (rad/s)
• From this plot, the phase margin is found to be –-6.02°, which means that
the gain-adjusted but uncompensated system is unstable.
• Noting that the addition of a lag compensator modifies the phase curve of the
Bode diagram, we must allow 5° to 12° to the specified phase margin to
compensate for the modification of the phase curve.
• Since the frequency corresponding to a phase margin of 40° is 0.856 rad/sec,
the new gain crossover frequency (of the compensated system) must be chosen
near this value.
• To avoid overly large time constants for the lag compensator, we shall choose
the corner frequency ω = 1/T (which corresponds to the zero of the lag
compensator) to be 0.1 rad/sec.
• Since this corner frequency is not too far below the new gain crossover
frequency, the modification in the phase curve may not be small. Hence, we
Page 3 of 10
Electrical Engineering Department
]txep TpyT[ ECOMP 423: Automatic Control (2)
4th year Computers And Systems Section
Final Examination 1st Term 2018/2019
Date: 02/01/2019 Time: 3.0 Hours
Model Answer (100 Marks) Faculty of Engineering
Aswan University
• add about 12° to the given phase margin as an allowance to account for the
lag angle introduced by the lag compensator.
• The required phase margin is now 52°.The phase angle of the uncompensated
open-loop transfer function is –128° at about ω =0.603 rad/sec. So we
choose the new gain crossover frequency to be 0.603 rad/sec.
• To bring the magnitude curve down to 0 dB at this new gain crossover
frequency, the lag compensator must give the necessary attenuation, which in
this case is –20 dB. Hence,
• Or β = 10
• The other corner frequency ω=1/(βT), which corresponds to the pole of the lag
compensator, is then determined as
• Thus, the transfer function of the lag compensator is
( ) ̂( ) = ̂
( ) ( )
Since the plant has an integrator, we use the second method of Ziegler-Nichols
tuning rules. By setting Ti = ∞ and Td = 0, we obtain the closed-loop transfer
function as follows: C ( s) Kp
R( s) s( s 1)(s 5) K p
The value of Kp that makes the system marginally stable so that sustained oscillation
occurs can be obtained by use of Routh's stability criterion. Since the characteristic
equation for the closed-loop system is
Page 4 of 10
Electrical Engineering Department
]txep TpyT[ ECOMP 423: Automatic Control (2)
4th year Computers And Systems Section
Final Examination 1st Term 2018/2019
Date: 02/01/2019 Time: 3.0 Hours
Model Answer (100 Marks) Faculty of Engineering
Aswan University
Examining the coefficients of the first column of the Routh table, we find that
sustained oscillation will occur if Kp = 30. Thus, the critical gain Kcr is
With gain K, set equal to , the characteristic equation becomes
√
Referring to Table 1, we determine Kp, Ti and Td as follows:
The PID controller has a pole at the origin and double zero at s = -1.4235
( )
b) Consider the following transfer function system: ( )
. Obtain the state-space
representation of this system in Controllable canonical, observable canonical and diagonal
forms. (8 marks)
Answer 3-b):
To obtain the state-space representation in controllable canonical form we rewrite
the previous equation in the following form:
( ) ( )
( ) ( )
From this last equation, the following two equations may be obtained:
( ) ( ) ( )
( ) ( )
Now define state variables as follows:
( ) ( )
( ) ( )
Page 5 of 10
Electrical Engineering Department
]txep TpyT[ ECOMP 423: Automatic Control (2)
4th year Computers And Systems Section
Final Examination 1st Term 2018/2019
Date: 02/01/2019 Time: 3.0 Hours
Model Answer (100 Marks) Faculty of Engineering
Aswan University
Then, clearly, ( ) ( )
Which may be rewritten as ̇ ( ) ( )
Noting that ( ) ( ) ( ) ( )
Or ̇ ( ) ( ) ()
̇ ( ) ( )
[ ] [ ][ ] [ ] ( )
̇ ( ) ( )
( )
( ) [ ][ ]
( )
To obtain the state-space representation in observable canonical form we rewrite
the previous equation in the following form: ( ) ( ) ( )
By dividing the entire equation by s2 and rearranging, we obtain
( ) ( ) ( )
Now define state variables as follows
( ) [ ( ) ( )]
( ) ( )
( ) ( )
By substituting the last equation into previous equations and multiplying both
sides of the equations by s, we obtain
( ) ( ) ( )
( ) ( )
Taking the inverse Laplace transforms of the preceding n equations and writing
them in the reverse order, we get
̇ ( ) ( ) ( )
̇ ( ) ()
( ) ( )
̇ ( ) ( )
Or [ ] [ ][ ] [ ] ( )
̇ ( ) ( )
( )
( ) [ ][ ]
( )
To obtain the state-space representation in diagonal form we use partial fraction
expansion then the previous equation will be in the following form:
( )
( ) ( ) ( )
Last equation may be written as ( ) ( ) ( )
( ) ( )
Define the state variables as follows:
( ) ( )
Page 6 of 10
Electrical Engineering Department
]txep TpyT[ ECOMP 423: Automatic Control (2)
4th year Computers And Systems Section
Final Examination 1st Term 2018/2019
Date: 02/01/2019 Time: 3.0 Hours
Model Answer (100 Marks) Faculty of Engineering
Aswan University
( ) ( )
( )
Which may be rewritten as
( ) ( )
( ) ( ) ( )
The inverse Laplace transforms of these equations give
̇ ( ) ()
̇ ( ) ( ) ()
( ) ( ) ( )
̇ ( ) ( )
Or [ ] [ ][ ] [ ] ( )
̇ ( ) ( )
( )
( ) [ ][ ]
( )
Question 4: (18 marks)
Consider a system defined by:
Thus, ,
For the determination of the state-feedback gain matrix, we use
[ ]
Where [ ] [ ]
Hence s the state-feedback gain matrix K is
[ ][ ][ ] [ ][ ] [ ]
Page 7 of 10
Electrical Engineering Department
]txep TpyT[ ECOMP 423: Automatic Control (2)
4th year Computers And Systems Section
Final Examination 1st Term 2018/2019
Date: 02/01/2019 Time: 3.0 Hours
Model Answer (100 Marks) Faculty of Engineering
Aswan University
Using this state-feedback gain matrix K, the control signal u is given by
( )
[ ][ ]
( )
Suppose that we use the observed-state feedback control instead of the actual-state
feedback control, or
̃ ( )
[ ][ ]
̃ ( )
b) Design a full order state observer, the desired eigenvalues for the observer matrix are
by applying Ackerman’s formula method. (9 marks)
Answer 4-b):
Ackerman's formula ( )[ ] [ ]
( ) [ ] [ ] [ ] [ ]
[ ][ ] [ ] [ ]
The observer equation is given by:
̇
( ̃) ( ̇) ̃
Since [ ]
̃
Then ̇
( ̃) ( ̇ )̃
Or
( ̃)̇ ̃
[ ] {[ ] [ ] [ ]} [ ] [ ]
̇
( ̃) ̃
( ̃)̇ ̃
[ ] [ ][ ] [ ]
( ̃)̇ ̃
Page 8 of 10
Electrical Engineering Department
]txep TpyT[ ECOMP 423: Automatic Control (2)
4th year Computers And Systems Section
Final Examination 1st Term 2018/2019
Date: 02/01/2019 Time: 3.0 Hours
Model Answer (100 Marks) Faculty of Engineering
Aswan University
Question 5: (24 marks)
a) Consider the following system: (12 marks)
̇
̇
̇
Let Liapunov function is:
( )
I) Show that ̇ ( ) ( )
Answer 5-a-I):
̇( ) ̇ ̇ ̇ ̇ ̇ ̇ ̇
( )
II) For what range of is this system stable?
Answer 5-a-II):
This system will be stable if ( ) positive definite ̇ ( ) is negative definite, then
( ) [ ][ ][ ]
[ ] [ ]
√ √
[ ]
√ √ √
[ ]
Then and
To be matrix P positive definite real symmetric matrix must
Page 9 of 10
Electrical Engineering Department
]txep TpyT[ ECOMP 423: Automatic Control (2)
4th year Computers And Systems Section
Final Examination 1st Term 2018/2019
Date: 02/01/2019 Time: 3.0 Hours
Model Answer (100 Marks) Faculty of Engineering
Aswan University
b) The state equations for the harmonic oscillator have the form (12 marks)
̇ () () ̇ () () ( ).
It is desired to find the optimal control signal u such that the performance index is minimized:
∫ ( )
Answer 5-b):
( ) ( ) ( ) ( )
[ ][ ] [ ][ ] [ ] [ ] [ ][ ][ ] [ ]
[ ]
[ ] [ ] [ ] [ ] [ ]
[ ] [ ]
[ ]
Substituting this P matrix into the following equation gives the optimal
K matrix
[ ][ ][ ] [ ]
Page 10 of 10