discrete inversion formulas for the design of lead and lag...
TRANSCRIPT
![Page 1: Discrete Inversion Formulas for the Design of Lead and Lag ...zanasi/didattica/Veicolo_OLD/ECC_2009_Budapest_Formule.pdfCompensator design: the discrete time case The discrete block](https://reader035.vdocument.in/reader035/viewer/2022081620/610e99d9bf67680fc57462bc/html5/thumbnails/1.jpg)
Discrete Inversion Formulas for the Design of Lead and Lag Discrete Compensators
R. Zanasi
Computer Science Engineering Department (DII) University of Modena and Reggio Emilia, Italy E-mail: [email protected]
R. Morselli
CNH Italia S.p.A.,V.le delle Nazioni 55, Modena, ItalyE-mail: [email protected]
![Page 2: Discrete Inversion Formulas for the Design of Lead and Lag ...zanasi/didattica/Veicolo_OLD/ECC_2009_Budapest_Formule.pdfCompensator design: the discrete time case The discrete block](https://reader035.vdocument.in/reader035/viewer/2022081620/610e99d9bf67680fc57462bc/html5/thumbnails/2.jpg)
Outline
1) Inversion formulas for the continuous time case 2) Design of lead/lag compensators 3) Discrete Inversion formulas4) Remarks5) Numerical examples6) Conclusions
The design of first or second order compensators is a very classical topic in automatic control courses.
It can be done in a lot of different ways and usually the Bode diagrams are used. All these different ways are basically equivalent. Which is the simplest and clearest way? This is an open question!
This paper presents continuous and discrete inversion formulas which have a nice graphical interpretation on the Nyquist plane and seems to be useful for teaching.
R. Zanasi, R. Morselli Compensator design using discrete inversion formulas
![Page 3: Discrete Inversion Formulas for the Design of Lead and Lag ...zanasi/didattica/Veicolo_OLD/ECC_2009_Budapest_Formule.pdfCompensator design: the discrete time case The discrete block](https://reader035.vdocument.in/reader035/viewer/2022081620/610e99d9bf67680fc57462bc/html5/thumbnails/3.jpg)
Compensator design: the continuous time case
R. Zanasi, R. Morselli Compensator design using discrete inversion formulas
The considered block scheme:
The compensator:
Steady state and high frequency gains:
Lead compensator if : Lag compensator if :
Bode magnitude and phase plots:
![Page 4: Discrete Inversion Formulas for the Design of Lead and Lag ...zanasi/didattica/Veicolo_OLD/ECC_2009_Budapest_Formule.pdfCompensator design: the discrete time case The discrete block](https://reader035.vdocument.in/reader035/viewer/2022081620/610e99d9bf67680fc57462bc/html5/thumbnails/4.jpg)
Inversion Formulas - continuous time case
The design problem can often be formulated as follows. Find the parameters of compensator C(s) such that:
where and are the magnitude and the phase desired at frequency .
The problem is solved by using the following Inversion Formulas [1],[2],[3]:
[1] The M. Policastro, G. Zonta, “Un procedimento di calcolo delle reti correttrici nella sintesi in frequenza di sistemi di controllo”, University of Trieste, Internal Report no. 88, 1982.
[2] The G. Marro and R. Zanasi, “New Formulae and Graphics for Compensator Design'', 1998 IEEE International Conference On Control Applications, Trieste, Italy, 1998.
[3] The G. Marro, “Controlli Automatici‘”, Zanichelli, Bologna, 2004.
Admissible domains:
R. Zanasi, R. Morselli Compensator design using discrete inversion formulas
![Page 5: Discrete Inversion Formulas for the Design of Lead and Lag ...zanasi/didattica/Veicolo_OLD/ECC_2009_Budapest_Formule.pdfCompensator design: the discrete time case The discrete block](https://reader035.vdocument.in/reader035/viewer/2022081620/610e99d9bf67680fc57462bc/html5/thumbnails/5.jpg)
Without the Inversion Formulas
Without the Inversion Formulas the control problem can be solved reading the compensator parameters from lookup tables:
R. Zanasi, R. Morselli Compensator design using discrete inversion formulas
Lead compensator
![Page 6: Discrete Inversion Formulas for the Design of Lead and Lag ...zanasi/didattica/Veicolo_OLD/ECC_2009_Budapest_Formule.pdfCompensator design: the discrete time case The discrete block](https://reader035.vdocument.in/reader035/viewer/2022081620/610e99d9bf67680fc57462bc/html5/thumbnails/6.jpg)
Compensator design on the Nyquist plane
1) Chose a point B on the Nyquist Plane.
2) Red semicircle: are all the points that can be moved in B using a lead compensator.
4) Green quarter of plane: are all the points that can be moved in B using a lag compensator.
3) Lead compensators with bouded amplification:
5) Lag compensators with bounded attenuation:
R. Zanasi, R. Morselli Compensator design using discrete inversion formulas
Lead
Lag
![Page 7: Discrete Inversion Formulas for the Design of Lead and Lag ...zanasi/didattica/Veicolo_OLD/ECC_2009_Budapest_Formule.pdfCompensator design: the discrete time case The discrete block](https://reader035.vdocument.in/reader035/viewer/2022081620/610e99d9bf67680fc57462bc/html5/thumbnails/7.jpg)
Lead compensator design: numerical example
The system: The goal:
Using the Inversion Formulas one obtains the copensator:
Point B follows from the phase margin:
Point A is chosen within the admissible region:
Parameters:
Step responceswith and withoutthe compensator.
R. Zanasi, R. Morselli Compensator design using discrete inversion formulas
![Page 8: Discrete Inversion Formulas for the Design of Lead and Lag ...zanasi/didattica/Veicolo_OLD/ECC_2009_Budapest_Formule.pdfCompensator design: the discrete time case The discrete block](https://reader035.vdocument.in/reader035/viewer/2022081620/610e99d9bf67680fc57462bc/html5/thumbnails/8.jpg)
Lag compensator design: numerical example
The goal:
Point B follows from the phase margin:
Step responceswith and withoutthe compensator.
The system:
Point A is chosen within the admissible region:
Parameters:
Using the Inversion Formulas one obtains the copensator:
R. Zanasi, R. Morselli Compensator design using discrete inversion formulas
![Page 9: Discrete Inversion Formulas for the Design of Lead and Lag ...zanasi/didattica/Veicolo_OLD/ECC_2009_Budapest_Formule.pdfCompensator design: the discrete time case The discrete block](https://reader035.vdocument.in/reader035/viewer/2022081620/610e99d9bf67680fc57462bc/html5/thumbnails/9.jpg)
Lag-Lead compensator: the continuous time case
The compensator (3 parameters):
The constraints:
The frequency response:
where:
R. Zanasi, R. Morselli Compensator design using discrete inversion formulas
Admissible values:
Inversion formulas:
a) if you choose :
b) if you choose :
c) if you choose :
Three possible design choices:
![Page 10: Discrete Inversion Formulas for the Design of Lead and Lag ...zanasi/didattica/Veicolo_OLD/ECC_2009_Budapest_Formule.pdfCompensator design: the discrete time case The discrete block](https://reader035.vdocument.in/reader035/viewer/2022081620/610e99d9bf67680fc57462bc/html5/thumbnails/10.jpg)
Lag-Lead compensator: a numerical example
Step responcewith the compensator.
R. Zanasi, R. Morselli Compensator design using discrete inversion formulas
The system:
The lag-lead compensator:
Point A is chosen within the admissible region:
The goal (B):
Point B follows from the phase margin:
Parameters:
Compensator parameters:
![Page 11: Discrete Inversion Formulas for the Design of Lead and Lag ...zanasi/didattica/Veicolo_OLD/ECC_2009_Budapest_Formule.pdfCompensator design: the discrete time case The discrete block](https://reader035.vdocument.in/reader035/viewer/2022081620/610e99d9bf67680fc57462bc/html5/thumbnails/11.jpg)
Compensator design on the Nichols plane
R. Zanasi, R. Morselli Compensator design using discrete inversion formulas
B
A1
A2
Lead
Lag
0-180-360
0
The Inversion Formulas can be used also on the Nichols plane.
The graphical interpretation of the Lead and Lagdomains is still present: they are symmetric
The Lead and Lag domains can be plotted “by hand” less precisely.
G(jw)
![Page 12: Discrete Inversion Formulas for the Design of Lead and Lag ...zanasi/didattica/Veicolo_OLD/ECC_2009_Budapest_Formule.pdfCompensator design: the discrete time case The discrete block](https://reader035.vdocument.in/reader035/viewer/2022081620/610e99d9bf67680fc57462bc/html5/thumbnails/12.jpg)
Compensator design: the discrete time case
The discrete block scheme:
Discrete first-order compensator:
is a minimum phase system only if:
Steady state and high frequency gains:
Bode magnitude and phase plots:
Lead compensator if :
Lag compensator if :
R. Zanasi, R. Morselli Compensator design using discrete inversion formulas
![Page 13: Discrete Inversion Formulas for the Design of Lead and Lag ...zanasi/didattica/Veicolo_OLD/ECC_2009_Budapest_Formule.pdfCompensator design: the discrete time case The discrete block](https://reader035.vdocument.in/reader035/viewer/2022081620/610e99d9bf67680fc57462bc/html5/thumbnails/13.jpg)
Discrete Inversion Formulas
The admissible domains are equal to the ones given for the continuous time case:
The discrete design problem is solved by using the Discrete Inversion Formulas:
Cd(s) is minimum phase if and only if C(s) is minimum phase:
R. Zanasi, R. Morselli Compensator design using discrete inversion formulas
The discrete design problem can often be formulated as follows. Find the parameters and of compensator Cd(s) such that:
where and are the magnitude and the phase desired at frequency .
![Page 14: Discrete Inversion Formulas for the Design of Lead and Lag ...zanasi/didattica/Veicolo_OLD/ECC_2009_Budapest_Formule.pdfCompensator design: the discrete time case The discrete block](https://reader035.vdocument.in/reader035/viewer/2022081620/610e99d9bf67680fc57462bc/html5/thumbnails/14.jpg)
Discrete Inversion Formulas: details
R. Zanasi, R. Morselli Compensator design using discrete inversion formulas
2) Discrete design problem:
1) Auxiliary variables:
where:
and:
3) The design problem can be rewritten as:
4) After some manipulations:
5) At the end one obtains:
![Page 15: Discrete Inversion Formulas for the Design of Lead and Lag ...zanasi/didattica/Veicolo_OLD/ECC_2009_Budapest_Formule.pdfCompensator design: the discrete time case The discrete block](https://reader035.vdocument.in/reader035/viewer/2022081620/610e99d9bf67680fc57462bc/html5/thumbnails/15.jpg)
Discrete Inversion Formulas: remarks
The discrete inversion formulas can also be rewritten in the following form:
where and are the parameters given by continuous time inversion formulas.
C(s) is transformed in Cd(s) by using the “bilinear transformation with prewarping”:
+ =
Poles and zeros of the continuous and discrete inversion formulas:
R. Zanasi, R. Morselli Compensator design using discrete inversion formulas
![Page 16: Discrete Inversion Formulas for the Design of Lead and Lag ...zanasi/didattica/Veicolo_OLD/ECC_2009_Budapest_Formule.pdfCompensator design: the discrete time case The discrete block](https://reader035.vdocument.in/reader035/viewer/2022081620/610e99d9bf67680fc57462bc/html5/thumbnails/16.jpg)
Discrete inversion formulas: numerical example
The system: The goal (B):
Step responceswith continuous,
discrete and bilinear
compensators.
The discrete system:
Point A is chosen within the admissible region:
Parameters:
With discrete inversion formulas:
With bilinar transformaion:
R. Zanasi, R. Morselli Compensator design using discrete inversion formulas
![Page 17: Discrete Inversion Formulas for the Design of Lead and Lag ...zanasi/didattica/Veicolo_OLD/ECC_2009_Budapest_Formule.pdfCompensator design: the discrete time case The discrete block](https://reader035.vdocument.in/reader035/viewer/2022081620/610e99d9bf67680fc57462bc/html5/thumbnails/17.jpg)
Discrete inversion formulas: remarks
R. Zanasi, R. Morselli Compensator design using discrete inversion formulas
Frequency parameters to move A -> B :
Using the discrete inversion formulas: Using the continuous inversion formulas:
one obtains the discrete compensator:
one obtains:
Then applying the “bilinear transformation with prewarping”:
one obtains the same discrete compensator !!!
![Page 18: Discrete Inversion Formulas for the Design of Lead and Lag ...zanasi/didattica/Veicolo_OLD/ECC_2009_Budapest_Formule.pdfCompensator design: the discrete time case The discrete block](https://reader035.vdocument.in/reader035/viewer/2022081620/610e99d9bf67680fc57462bc/html5/thumbnails/18.jpg)
1) Inversion Formulas – continuous case
Easy to use (Lead, Lag and Lead/Lag compensators)
Simple graphical representations on Nyquist plane
Easy understandable by undergraduate students
2) Inversion Formulas – discrete case
Similar to the continuous time case
Same domains and same graphical representation
on the Nyquist plane
Useful for teaching and for discrete control design.
Conclusions
R. Zanasi, R. Morselli Compensator design using discrete inversion formulas