Control Design for Electronic Power Converters

PhD Thesis Control Design for Electronic Power Converters By Carolina Albea Sanchez´ Supervised by Carlos Canudas de Wit and Francisco Gordillo...

0 downloads 173 Views 2MB Size
Control Design for Electronic Power Converters Carolina Albea

To cite this version: Carolina Albea. Control Design for Electronic Power Converters. Automatic. Institut National Polytechnique de Grenoble - INPG; Universidad de Sevilla, 2010. English.

HAL Id: tel-00539077 https://tel.archives-ouvertes.fr/tel-00539077 Submitted on 24 Nov 2010

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L’archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d’enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

UNIVERSITE DE GRENOBLE INSTITUT POLYTECHNIQUE DE GRENOBLE No. attribu´e par la biblioth`eque THESE EN COTUTELLE INTERNATIONALE pour obtenir le grade de DOCTEUR DE L’Universit´e de Grenoble d´elivr´e par l’Institut polytechnique de Grenoble et de L’Universit´e de Sevilla (Espagne) Sp´ecialit´e: AUTOMATIQUE-PRODUCTIQUE pr´epar´ee au D´epartement Automatique du GIPSA-lab dans le cadre de l’Ecole Doctorale Electronique, Electrotechnique, Automatique, Traitement du Signal et au laboratoire Departamento de Ingenier´ıa de Sistemas y Autim´atica pr´esent´ee et soutenue publiquement par

´ Carolina ALBEA-SANCHEZ le 27/09/2010 TITRE

Control Design for Electronic Power Converters DIRECTEURS DE THESE M. Carlos Canudas de Wit ´ M. Francisco Gordillo Alvarez

Directeur de Recherche CNRS Professeur, Universidad de Sevilla JURY

M. Javier Aracil Santoja M. Wilfrid Perruquetti M. Enric Fossas Colet M. Laurent Fesquet M. Luis, Martinez Salamero

Professeur, Universidad de Sevilla, Professeur, Ecole Centrale de Lille, Professeur, Universit´e Polit`ecnica de Catalunya, Maitre de Conf´erence de l’INPG, Professeur, Universit´e de Tarragona,

Pr´esident Rapporteur Rapporteur Examinateur Examinateur

PhD Thesis

Control Design for Electronic Power Converters

By Carolina Albea S´anchez Supervised by Carlos Canudas de Wit and Francisco Gordillo

A mi familia y a Alexandre

Contents

Aknowledgements

vii

Notation

ix

1

Introduction

1

1.1

Introduction to power converters . . . . . . . . . . . . . . . . . . . . . . .

1

1.1.1

Converters classification . . . . . . . . . . . . . . . . . . . . . . .

2

1.1.2

DC-DC converters . . . . . . . . . . . . . . . . . . . . . . . . . .

2

1.1.3

DC-AC converter . . . . . . . . . . . . . . . . . . . . . . . . . . .

3

1.1.4

AC-DC converter . . . . . . . . . . . . . . . . . . . . . . . . . . .

3

1.1.5

AC-AC converter . . . . . . . . . . . . . . . . . . . . . . . . . . .

4

Research motivation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

4

1.2.1

Boost inverter . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

4

1.2.2

DC-DC Vdd-Hopping converter . . . . . . . . . . . . . . . . . . .

5

Main objectives . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

6

1.3.1

Controlling a boost inverter . . . . . . . . . . . . . . . . . . . . .

7

1.3.2

Controlling a DC-DC Vdd-hopping converters . . . . . . . . . . .

7

Thesis structure . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

8

1.2

1.3

1.4

i

I Controlling a DC-AC Boost Converter 2 Introduction 2.1

2.2

9 11

Boost inverter . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

12

2.1.1

System description . . . . . . . . . . . . . . . . . . . . . . . . . .

13

Control problem objectives . . . . . . . . . . . . . . . . . . . . . . . . . .

16

2.2.1

Control law . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

17

2.2.2

Adaptive control . . . . . . . . . . . . . . . . . . . . . . . . . . .

20

2.2.3

Attraction domain . . . . . . . . . . . . . . . . . . . . . . . . . .

21

3 Control of the DC-AC boost converter by energy shaping

23

3.1

Normalized average model . . . . . . . . . . . . . . . . . . . . . . . . . .

23

3.2

Energy shaping control for generation of oscillations . . . . . . . . . . . .

25

3.2.1

Approach overview . . . . . . . . . . . . . . . . . . . . . . . . . .

25

3.2.2

Controller design . . . . . . . . . . . . . . . . . . . . . . . . . . .

27

3.2.3

Control law for the full system . . . . . . . . . . . . . . . . . . . .

29

3.2.4

Simulation results . . . . . . . . . . . . . . . . . . . . . . . . . .

30

Synchronization problem . . . . . . . . . . . . . . . . . . . . . . . . . . .

33

3.3.1

Boost inverter synchronization . . . . . . . . . . . . . . . . . . . .

33

3.3.2

Synchronization with the electrical grid . . . . . . . . . . . . . . .

38

Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

39

3.3

3.4

4 Adaptive control 4.1

41

Design of an adaptive control . . . . . . . . . . . . . . . . . . . . . . . . .

41

4.1.1

43

Adaptation law . . . . . . . . . . . . . . . . . . . . . . . . . . . .

4.2

4.3

5

4.1.2

Error equation . . . . . . . . . . . . . . . . . . . . . . . . . . . .

43

4.1.3

Stability properties . . . . . . . . . . . . . . . . . . . . . . . . . .

44

Stability considerations of the full closed-loop system . . . . . . . . . . . .

45

4.2.1

Tuned system . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

45

4.2.2

Closed-loop system . . . . . . . . . . . . . . . . . . . . . . . . . .

46

4.2.3

Singular perturbation form . . . . . . . . . . . . . . . . . . . . . .

47

4.2.4

Slow sub-system . . . . . . . . . . . . . . . . . . . . . . . . . . .

49

4.2.5

Boundary layer fast subsystem . . . . . . . . . . . . . . . . . . . .

49

4.2.6

Simulations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

52

Conclusions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

52

Estimation of the attraction domain

57

5.1

Problem statement . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

58

5.2

An approach of estimation of the attraction domain. . . . . . . . . . . . . .

59

5.2.1

Sum of squares optimization . . . . . . . . . . . . . . . . . . . . .

62

5.3

Application to the boost inverter . . . . . . . . . . . . . . . . . . . . . . .

63

5.4

Conclusions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

65

II Controlling a DC-DC Vdd-Hopping converter

67

6

Introduction

69

6.1

Optimization of the energy consumption in SoCs . . . . . . . . . . . . . .

71

6.2

Vdd-Hopping DC-DC converter . . . . . . . . . . . . . . . . . . . . . . .

72

6.3

Non-linear control application to Vdd-Hopping DC-DC converter . . . . .

74

6.3.1

High-performance controllers . . . . . . . . . . . . . . . . . . . .

75

6.3.2

Energy-aware controller . . . . . . . . . . . . . . . . . . . . . . .

76

6.3.3

Approximate stability analysis for the energy-aware controller . . .

77

6.3.4

Advanced energy-aware controller . . . . . . . . . . . . . . . . . .

78

7 High-performance control for the DC-DC Vdd-Hopping converter 7.1

81

Mathematical model of Vdd-Hopping mechanism . . . . . . . . . . . . . .

82

7.1.1

Mathematical model for control design . . . . . . . . . . . . . . .

82

Control laws . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

85

7.2.1

Control proposed in [99] . . . . . . . . . . . . . . . . . . . . . . .

85

7.2.2

Controller No. 1: linear controller . . . . . . . . . . . . . . . . . .

86

7.2.3

Controller No.2: feedback linearization . . . . . . . . . . . . . . .

89

7.2.4

Controller No.3: Lyapunov-based design . . . . . . . . . . . . . .

90

Performance evaluation . . . . . . . . . . . . . . . . . . . . . . . . . . . .

93

7.3.1

Energy evaluation . . . . . . . . . . . . . . . . . . . . . . . . . . .

94

7.3.2

Summary . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

95

Advanced Lyapunov’s controller . . . . . . . . . . . . . . . . . . . . . . .

97

7.4.1

Optimal voltage reference computation . . . . . . . . . . . . . . .

98

7.4.2

Adaptive feedback control design . . . . . . . . . . . . . . . . . .

103

7.5

Simulation of the advanced Lyapunov’s controller . . . . . . . . . . . . . .

104

7.6

Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

106

7.2

7.3

7.4

8 Energy-aware controller for the Vdd-Hopping converter 8.1

Control design without current-peak managing . . . . . . . . . . . . . . .

109 110

8.2

8.3 9

Control redesign with current-peak managing . . . . . . . . . . . . . . . .

110

8.2.1

Time discretization . . . . . . . . . . . . . . . . . . . . . . . . . .

112

8.2.2

Simulation of ENARC controller in the Vdd-Hopping system . . .

113

Conclusions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

115

Approximate stability analysis of the DC-DC Vdd-Hopping converter

117

9.1

Stability with control (7.10) . . . . . . . . . . . . . . . . . . . . . . . . . .

118

9.2

Stability with control (8.6) . . . . . . . . . . . . . . . . . . . . . . . . . .

121

9.3

Conclusions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

127

10 Sub-optimal control considering delays and parameter uncertainties

129

10.1 Problem statement . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

131

10.1.1 Alternative representation for the saturated control (10.2) and the error equation (10.1) . . . . . . . . . . . . . . . . . . . . . . . . .

132

10.1.2 State-space representation . . . . . . . . . . . . . . . . . . . . . .

133

10.1.3 Stability and disturbance rejection problem . . . . . . . . . . . . .

133

10.2 H∞ control design . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

134

10.2.1 Descriptor model transformation . . . . . . . . . . . . . . . . . . .

134

10.2.2 Condition for state-space representation . . . . . . . . . . . . . . .

135

10.2.3 Control design . . . . . . . . . . . . . . . . . . . . . . . . . . . .

136

10.3 Robust control tuning . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

137

10.4 Sub-optimal control result . . . . . . . . . . . . . . . . . . . . . . . . . .

139

10.5 Simulation Results . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

140

10.5.1 Uncertain clock frequency. . . . . . . . . . . . . . . . . . . . . . .

140

10.5.2 Uncertain PMOS resistance . . . . . . . . . . . . . . . . . . . . .

140

10.5.3 Uncertain load parameter . . . . . . . . . . . . . . . . . . . . . . .

140

10.6 Evaluation of the tuning methods . . . . . . . . . . . . . . . . . . . . . . .

143

10.7 Conclusions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

145

12 Conclusions and future work

147

12.1 Conclusions and contribution summary . . . . . . . . . . . . . . . . . . .

147

12.2 Future work . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

150

List of Publications

152

Aknowledgements The period of this thesis brought me valuable new experiences in both a scientific and personal perspective. This thesis would not have been possible without the constant support and motivation of many people. Firstly, I would like to express my sincere gratitude to my two supervisors, the Professor Francisco Gordillo at the ‘Universidad de Sevilla’ (Spain), in the Nonlinear Control group from the ‘Departmento de Ingen´ıa de Sistemas y Autom´atica’; and the Director of Research Carlos Canudas de Wit at the CNRS (France), in the NeCS team from the ‘D´epartement d’Automatique de GIPSA-Lab’ at the ‘Institut Polytechnique de Grenoble’ and at the ‘Institut National de Recherche en Informatique et en Automatique de Grenoble’ (France). Their wide knowledge and their experience have been of great value for me. Their encouraging and personal guidance have provided a good background for the present thesis. Many thanks to the Professor Wildfrid Perruquetti at the ‘l’Ecole Central de Lille’ (France), and the Professor Enric Fossas Colet at the ‘Universit Polit`ecnica de Catalunya ’ (Spain), for accepting to review my work as well as for participating in my jury committee. Their constructive suggestions on the thesis are really appreciated for me. I am grateful that in the midst of their activities, they accepted these tasks. I am also grateful the Professor Javier Aracil at the ‘Universidad de Sevilla’ (Spain), the Maitre de Conference Laurent Fesquet at the ‘Institut Polytechnique de Grenoble’ (France) and the Professor Luis Martinez Salamero at the ‘Universit´e de Tarragona’ (Spain), for their commitment to take part in my jury committee. During this period, I have collaborated with many colleagues, from PhD students, teachers to professors for whom I have great regard, and I wish to extend my warmest thanks to all those who have helped me from a professional and personal point of view. Especially to Alicia, Mirko, Antonio, Isa, Pablo, Alejandro, J¨orn, Guilherme, Vicente, Fabio, Manolo ´ Gil, Manolo Vargas, Manolo Ruiz, Jose Angel, Manolo Lopez, Javier, Eduardo, Paco Rubio, Paco Salas, Carlos Vivas, Fernando, Asun, Amparo and David in the ‘Departmento de Ingen´ıa de Sistemas y Autom´atica’, at the ‘Universidad de Sevilla’. And thanks to Luc, Lara, Sylvain, Riccardo, Nicolas, Emilie, Gabriel, Denis, Wenjuan, Brandon, Katerina, John-Jairo, vii

Gregory, Daniel, Nicolas and Olivier in the ‘D´epartement d’Automatique de GIPSA-Lab’ at the ‘lInstitut Polytechnique de Grenoble’ and at the ‘Institut National de Recherche en Informatique et en Automatique de Grenoble’. I am grateful to the secretaries and librarians in the automatic departments from Spain, (Manolo, Antonio, Silvia) and France, (Virginie, Marie-Rose, Marie-Therese, Patricia, Myriam, Elodie) for helping me to run smoothly and for assisting me in many different ways. I wish to thank my best friends from Spain, JuanMa, Maribel, Olga, compa˜neros de ISF, Manolo, Sonia, and from France, Johana, Mikel, Luis, Eduardo, for helping me through difficult periods, and for providing to me their friendship, entertainment and attention. My warm special thanks to my family, especially to my mother Marilo, my sister Myrian, Domingo and my father Pep´ın, for their consideration of my research abroad. And my deep lovely thanks to Alexandre for its emotional supports. Without their encouragement and understanding it would have been impossible for me to finish this work. Lastly, I offer my regards and blessings to all of those who supported me in any respect during the completion of the project.

Notation The following symbols and conventions will be used consistently throughout the thesis.

≡ 6 ≡ ≈ , < (>) ≤ (≥) ≪ (≫) ± ∀ ∈ ⊂ ∩ : → xT ∞ ∑ |x| kxk kxk p [a, b] {1, 2, ..., N} xk q−1 Rn R+ N Z

identically equal not identically equal approximately equal defined as less (greater) than less (greater) than or equal to much less (greater) than plus and minus for all belongs to subset of intersected with such as tend to the transpose of a vector x infinity summation the absolute value of x the norm of a vector x the p-norm of a vector x closed interval from a to b the set from 1 to N discrete variable discrete-time delay or parameter, i.e., xk−1 = q−1 xk the n-dimensional Euclidean space the semi-positive-dimensional Euclidean space the set of natural numbers the set of integer numbers

ix

Re z (x, y) In diag[a1 , ..., an] O(·) f : S1 → S2 ∂f ∂x



d dt λmin(P)

P>0 P≥ e s sign max min exp sin cos dist(p, M) limx→c f (x) satM m (x) C o is the convex hull of a set round(x) ζ+ ζ− ∆ζ , ζ + − ζ − L2 H∞   [xx]

the real part of a complex variable z metric space n-dimensional identity matrix a diagonal matrix with diagonal elements a1 to an order of magnitude notation a function f mapping a set S1 into a set S2 the Jacobian matrix the first derivative of y with respect to time the first derivative of y with respect to time the minimum eigenvalue of symmetric matrix P a positive definite matrix a positive semidefinite matrix P neperian number the Laplace variable sign function maximum minimum exponential function sine cosine the distance from a point p to a set M limit of f (x)  as x approaches c  M if x > M x if m ≤ x ≤ M is defined as  m if x < m.

is the nearest integer to x denotes ζ (k + 1) denotes ζ (k) the value of ζ in two consecutive sampling time T is the space of {x} with the norm: kxk22 , ∑∞ k=0 x x < ∞ H-infinite method designation of the end of definition designation of the end of proof see reference number xx in the bibliography

Acronyms AC BVP DC DVS ENARC GALS HPF IVP LDVS LMI LPF PC PHC PI PLL PMOS PSS QoS SDP SMPC SoC SOS THD VLSI VHDL-AMS

Alterning current Boundary Value Problem Direct current Dynamic Voltage Scaling ENergy Aware Controller Globally-Asynchronous and Locally Synchronous Systems High Pass Filter Initial Value Problem Local Dynamic Voltage Scaling Linear Matrix Inequalities Low Pass Filter Personal Computer PHase Controller Proportional-Integral Phase-Lock Loop P-channel Metal-Oxide-Semiconductor field-effect transistor Power Supply Selector Quality of Service Semidefinite Program Switched-Mode Power Converter System on Chip Sum Of Squares Total Harmonic Distortion Very Large Scaling Integration Verilog Hardware Description Language-Analog and Mixed-Signal

Chapter 1 Introduction 1.1 Introduction to power converters Power converters are electronic circuits associated to the conversion, control, and conditioning of electric power. The power range can be from milliwatts, mobile phone, for example, to megawatts, in electric power transmission systems. Reliability of the power converters become a key industrial focus. Electronic devices and control circuit must be highly robust in order to achieve a high useful life. A special accent must be set on the total efficiency of the power electronic circuits. Firstly, because of the economic and environmental value of wasted power and, secondly, because of the cost of energy dissipated that it can generate. Even a small improvement in converter power efficiency translates to improved profitability of the investment in the electronic market [33, 100]. Among all electronic converters, the most common technology is switched-mode power converters (SMPC) [118]. They convert the voltage input to another voltage signal, by storing the input energy temporarily and then releasing that energy to the output at a different voltage. This switched-mode conversion has a particular interest due to the fact that it can switch at high frequency in a very efficient way. Power is controlled (even modified) by controlling the timing that the electronic switches are “on” and “off”. A much greater emphasis is required on achieving high-power efficiency in low-power level electronic technology, since few low-power circuits can tolerate a power efficiency less than 85%. Converters are used in these circuits in order to change the supply voltage in the blocks of the System on Chips (SoCs) according to performance requirements, for power efficiency reasons. Research have been focused on developing electronic circuits that can be employed as switches. e.g. approximating ideal closed or open switches, as the Vdd-hopping converter [98]. 1

1.1. Introduction to power converters

2

1.1.1 Converters classification Power converters control the flow of power between two systems by changing the character of electrical energy: from direct current to alternating current, from one voltage level to another voltage, or in some other way. Here, some important way to classify the power converters are described. The aim of this section is not to make a rigorous converter classification, either to make a state of the art, because it is not the purpose of this thesis. It is only desired to understand some properties of these kind of circuits. The most common classification of power conversion systems is based on the waveform of the input and output signals, in the case whether they are alternating current (AC) or direct current (DC) [33], thus: • DC to DC. • DC to AC. Inverter. • AC to DC. Rectifier. • AC to AC. Transformer. At the same time, the devices within converters can be switched in different ways [72, 79, 100]. If the devices switch at the line frequency (normally, 50Hz or 60HZ), they can be line frequency converters (naturally commutated converters) or high-frequency switching (forced-commutated converters). Depending on the character of the input source, they may be voltage-source converters or current-source converters. Moreover, converters may be of low, medium or high voltage and/or current level. Another sort of classification may be performed according to the size of the output signal obtained from the input signal; if the converter accomplishes a lower output signal it is well known as step-down, and if it obtains a larger signal, it is known as step-up [2].

1.1.2 DC-DC converters DC-DC converters are electronic circuits that change the DC operating voltage or current. They have recently aroused the interest in the current market due to its wide range of applicability. Normally, they are designed in order to transfer power from the input to the output

Chapter 1. Introduction

3

in one direction. However, in the case of switches topologies, the power moving may be also bidirectional, being very useful to develop new converter topologies for other applications, as can be an inverter topology [25]. They have a particular interest in low-power circuits, as cellular phone and personal computers (PCs). This sort of technology are composed of many sub-circuits that require an own voltage level from an external supply (higher, lower or even negative) or battery. DC-DC converters have a special role in these kind of systems, since they can be employed to change the voltage from a partial lowered battery voltage thereby. This is based on the Dynamic Voltage Scaling technic (DVS) [27, 89, 136]. The main idea of DVS is to vary the supply voltage in order to consume a minimal amount of energy. This fact improves the power efficiency and saves space in spite of using multiple batteries to accomplish the same voltage level [82].

1.1.3 DC-AC converter DC-AC converters, or commonly named inverters, can obtain a certain amplitude and frequency of the AC voltage and/or current without using normally an intermediate DC stage. This electrical device is a power electronic oscillator [118]. An electronic oscillator is just an electric circuit that produces a repetitive signal, as a sine-wave output signal. Generally, they are SMPCs. These kind of circuits require an efficient control for the switches devices that, in many occasions, can be quite complex due to system structure. Therefore, to design a suitable control law currently is a subject of much research [22, 107].

1.1.4 AC-DC converter The process that converts AC to DC is known as rectification, hence, these converters are also called rectifier. Among others applications, they are used in power supplies and detector of radio signals. The rectification can be half-wave or full-wave. In the first case (half-wave rectification), only one half of the input waveform can be employed to reach the desired output. Therefore, only this half AC wave (positive or negative) is converted. The efficiency will depend of the kind of application. It is clear, that it is not useful for power transfer. The full-wave rectification can convert the whole of the input waveform to achieve the constant output signal. It becomes more efficient [59, 137].

4

1.2. Research motivation

1.1.5 AC-AC converter AC-AC converters are employed to transform an AC input signal to another AC output signal with an arbitrary amplitude. Likewise, depending on the converter complexity, the frequency can be changed as well. The efficiency of these kind of systems depends on the type of circuit employed. It is clear that a higher power density and reliability will be obtained with a conversion in one single stage [139].

1.2 Research motivation A lot of research has recently been focused on converters due to the increasing deal of interest in power electronics. This is mainly caused by their broad applicability domain that includes battery-operating portable equipment, computers, appliances, vehicles, industrial electronic equipment, uninterruptible power supplies, telecommunication systems and much more. This current research is specially focused on finding highly-efficient converter topologies for every system application and, on designing control mechanisms that accomplishes the converter objectives. On this way, one or more electrical parameters can be regulated with a high reliability and efficiency, e.g., the supply voltage of an appliance, the temperature of an oven, the speed of a motor, the supply voltage within calcul node of a SoCs [41, 141]. Tackling the control problem in detail for every converter is out of the scope of a thesis. That is why, among all variety of converters, this thesis is focused on providing a control solution for two converters topologies, which have some interesting properties and applications. The converters that will be dealt with are: firstly, a switch inverter topology; and secondly, a DC-DC converter for low power application.

1.2.1 Boost inverter As was said before, inverters are devices that obtain a current output signal capables from passing through zero. The inverters are generally SMPCs, and their topologies are derived from coupling one or more basic switch topologies. Among them, it can be found the boostbuck inverter [95], the buck-boost inverter [90], the buck inverter [127], the boost inverter [25]. The first part of this thesis is focused on a boost inverter. Its interest is due to its step-up property, which is achieved through a signal stage. In this case, two DC-DC boost converters are connected with a load between them, thus it has a bidirectional current [25].

Chapter 1. Introduction

5

In this system, four switches have to be controlled by two control signals, in order to control the two output voltages of each DC-DC boost converter. Not only the voltage amplitude must be controlled, but the phase of both signal must also be controlled to achieve the specified output voltage. Several control laws have been designed for this converter from other authors [25, 126, 149], applying different control strategies. In this thesis, a novel control strategy based on energy shaping for generation of oscillations is employed [16,58], being the control objective a limit cycle. The novelty of this method is that it does not need to track any reference signal to achieve an oscillatory character in the output signal.

1.2.2 DC-DC Vdd-Hopping converter The second part of the thesis is focused on a DC-DC converter employed in low-power applications. As mentioned, the demand for high efficiency DC-DC converters is increasing dramatically, especially in battery-operated devices such as cellular phones and personal computers. In SoCs, to extend battery life has a particular role. By employing DC-DC converters based on power-saving, power efficiency in SoCs can be significantly increased, thereby extending battery life. The goal of these efficiency DC-DC converters is to adapt dynamically the supply voltage of the chip according to the required performance level. This is the DVS idea mentioned before. Numerous DC-DC converters employed for this aim have been proposed over the years to increase the power efficiency of an SoCs. The most commonly used topologies in DC-DC converters in low-power electronics are: continuous buck converters [119, 143, 153], boost converters [36], buck-boost [125, 142] converters and charge pump [125], among others. However, while converters may decrease conduction losses, additional losses can be added if switched devices are employed. In low-power applications where a high-efficiency is required, other different topologies far from switched-mode are employed. In [98], a discrete DC-DC converter was proposed based on the ‘Vdd-Hopping’ technique. This method is expired in scaling the voltage supply Vdd in a discrete way, delivering two small voltage levels according to the optimum Vdd required for every performance [75, 106, 128]. Hence its name of ‘DC-DC Vdd-Hopping converter’. Therefore, this technique replaces the continuous adjustable voltage, just to two set-points [83], so that it reaches a high-efficiency and reduced size. Likewise, it is a very simple system, becoming easily controllable [55]. This converter is employed in a French gouvernement project, with a very ambitious objective: ‘to reduce the size of the SoCs to 32nm’. For this, a new technology must be

1.3. Main objectives

6

Converter Boost Inverter

Vdd-Hopping Converter

Power level More suitable for medium and high power Low power

Conversion DC-AC

DC-DC

Scales normal

Model order 4th

microor nano-scales

1st

Table 1.1: Main differences between boost inverter and Vdd-Hopping converter

developed, since the currently technology applied to 45nm can not be employed for physical reasons. Here, the control may take a particular role since, a suitable control law can achieve the equilibrium and the demanded requirements, as well. For instance, it has to achieve the highest efficiency (among other goals) to achieve the global project objective.

1.3 Main objectives

The two selected converters have different natures and applications, and hence, they may have different control objectives. They covers a wide range of the power converter domain. The boost inverter normally is applied to medium and high power level for normal scales; and DC-DC Vdd-hopping converter is used in low-power technology for micro-scales or nano-scales. Likewise, conversions are DC-to-DC, and DC-to-AC. The DC-AC converter is based on the switched-mode classical topology, as is the boost inverter; and the DC-DC converter has a topology far from the common structures. The complexity of the systems are quite different, from a 1st -order model in the DC-DC converter to a 4th -order. Table 1.1 summarizes these differences. That does not mean that these two applications completely cover all power converter domain. In fact, there are other features that have not been taken into account, as is the different natural- or forced-commuted characteristic, the input sources, the level of the output signal, among others. These two converter applications, as for its work context as for its different characteristics, have some different control objectives.

Chapter 1. Introduction

7

1.3.1 Controlling a boost inverter The first application is focused on controlling an SMPC boost inverter. This converter is particularly interesting because it does not only allow to generate an alternating current, but it can also obtain an output voltage larger than the input signal. It has a high efficiency due to its switching character. Nevertheless, it has a non-minimum phase, 4th -order model. In addition, the desired behavior is not an equilibrium point but a limit cycle. Due to all the mentioned boost inverter characteristics, the main objective is to design a control law that guarantees not only the convergence to the desired limit cycle, but also the stability of it, with the particularity that no external reference is applied to the system. Likewise, the system has to accomplish right performance not only for known loads, but also for unknown loads. Another important aim is to estimate a set of initial voltage and current values, for which the system variables tend to the desired limit cycles when the control law before is applied to the boost inverter. If all these objectives are achieved, a control system guarantees a stable and robust behavior from an initial condition, which is within an estimated attraction region. And, in addition, the system is autonomous in the sense that no reference signals are needed.

1.3.2 Controlling a DC-DC Vdd-hopping converters The second application deals with the control of a discrete DC-DC Vdd-Hooping converter. This is a low-power converter with a high-efficiency. Furthermore, it has nice properties, for instance, it has a 1st -order and its control objective is an equilibrium. Nevertheless, in lowpower technology, this level of efficiency may not become enough if certain requirements are demanded (e.g. high energy-efficiency, small current peaks and reduced space) to achieve a certain objective. For this, to design a control law focused on achieving an optimal energyefficiency may be an attractive control problem in order to reach this objective. Indeed, the control problem of the Vdd-Hopping converter in this thesis comes directly demanded by the industry. Concretely, it is included in a French national project called ARAVIS, sponsored by the global competitive cluster Minalogic1. The main control objective of this converter is to guarantee that the system reaches the desired equilibrium point, achieving certain required features as: high-efficiency, stability, low computational cost, robustness with respect to parameter uncertainties and robustness with respect to delays due to synchronization and computation issues [45]. In this way, the control law must be designed taken these objectives into account. 1 http://www.minalogic.com/

8

1.4. Thesis structure

1.4 Thesis structure This thesis, as is noted above, is composed of two parts. Part I covers Chapters 2 to 5 while Part II covers Chapter 6 to 10. Conclusions are drawn in Chapter 11. The first part deals with controlling the boost inverter. In Chapter 2, the model of the double boost converter (boost inverter) is presented. Likewise, the objectives are specified in details, just as a particular solution is proposed in order to resolve the raised problem. Chapter 3 shows the general idea of producing oscillating behavior by means of the generation of a limit cycle through energy shaping. This idea yields a controller for the boost inverter, but it is shown that the behavior is not acceptable due to a lack of synchronization. Therefore, a phase controller is added to achieve the synchronization of an isolated boost inverter as well as the synchronization of the boost inverter with a pre-specified signal. Chapter 4 deals with the unknown-load case, which is solved by means of an adaptation mechanism design. A stability analysis for the full-system is also studied by using singular perturbation analysis. Chapter 5 is devoted to develop a method of estimating an attraction domain. This method deals with control and state constraints. It is employed to provide an estimated attraction domain for the boost inverter. The second part of the thesis is focused on controlling the DC-DC Vdd-Hopping converter. In Chapter 6, a summary of the ARAVIS project work context, where this research is included, is performed. Likewise, the control objectives required for this DC-DC converter in the ARAVIS project are defined. In Chapter 7, a set of controllers are presented and discussed. From the control solution that offers the best performance, a controller is developed in order to achieve the control objectives. For this, optimal control theory as well as adaptation methods are applied. Nevertheless, it has an important drawback, its implementation is not simple, thus it is not suitable in the ARAVIS project. Next, another controller is developed in Chapter 8. This proposed control solution is developed from the simplest control implementation of the set of controllers presented in Chapter 7. This controller presents good properties for the project. In Chapter 9, a rigorous stability analysis is developed for the closed-loop system with this last controller. Chapter 10 presents an optimal tuning mechanism for the control constants in order to deal with delays and parameter uncertainties. This development copes with resolving a H∞ problem, proposing some Linear Matrix Inequalities (LMIs) developed from Lyapunov Krasovskii method.

Part I Controlling a DC-AC Boost Converter

9

Chapter 2 Introduction DC-DC power converters have a very large presence in all kind of electronic circuits, from industrial applications (spacecraft power systems, DC motor drives, telecommunication equipment) to personal applications (PCs, office equipment, electrical appliance). These systems provide a regulated DC voltage level (Vo ) from an unregulated DC voltage level (Vin ). High efficiency is the most important requirement for DC-DC converters in a wide range of load power, since it directly affects the battery lifetime [42]. It can be achieved using ‘switched-mode’. A switched-mode power converter (SMPC) is characterized by rapidly switching on and off some devices, transferring a rate of energy from the input to the output. This rate of energy is controlled by a duty cycle1 to minimize the dissipated energy. The switching effect is achieved by transistors, which dissipates little power when it is outside of its active region. In addition, SMPCs have an inductor, whose main function is to limit the current slew rate through the power switch. This action help to limit the otherwise high peak current. Moreover, the inductor stores the energy, which can be recovered in the discharge phase [43, 118]. This approach is also used in alternating current (AC) applications. The basic components of the switching circuit can be rearranged to form a: 1. Buck converter. It is a step-down: the output voltage is lower than the input voltage. 2. Boost converter. It is a step-up: the output voltage is larger than the input voltage. 3. Buck-boost converter. It can be a step-down or a step-up. Its main characteristic is that it inverts the polarity of the voltage. 4. Cuk converter. It has the same features that a buck-boost converter, with other different configuration. 1 Duty

cycle is the fraction of time that a system is operated.

11

2.1. Boost inverter

12

5. SEPIC converter. It can be a step-down or a step-up, but it does not invert the voltage polarity.

From these topologies other converters can be obtained [96, 100]. A buck or boost topology, by oneself, can not achieve alternating current. Physical reasons prevent the output current signal from passing through zero. Hence, some topologies have been proposed in order to obtain the alternating current condition. Traditionally, DC-AC converters (or inverters) are based on the buck topology. Nevertheless, this kind of configurations obtain an AC output voltage lower in amplitude than the input voltage [90]. In applications that require a boosting output, this problem is solved by using two-stages. One-stage to change the signal from DC to AC, and the other stage, to raise the amplitude [121]. These topologies have the drawback of needing more space and dissipating more energy since they use more components. In [25], a new inverter was proposed composed of two boost converter. It is known as boost inverter. This inverter has as main advantage that it generates an AC output voltage from a lower DC voltage in a single stage. As side effect, it has a higher efficiency and a better signal quality with respect to the traditional buck inverters [126, 147]. These nice properties are only achieved with a suitable controller. Hence, to design an appropriate control law has an important relevance for these kind of circuits. The boost inverter may be used in diverse applications, as for example in photovoltaic system market. The solar cells can charge a battery up with a DC voltage of 48V . When they are used in domestic installations, a standard domestic AC power is required as power supply [15, 29]. Therefore, a boost inverter provides in these kind of applications a better benefit. Its structure allows to isolate as well as to increase the voltage. Moreover, it ensures that the power conversion is done with reduced energy losses [3, 4]. Figure 2.1 represents a domestic photovoltaic installation.

2.1 Boost inverter A boost circuit is usually employed as a DC-DC converter, being especially interesting because it generates an output voltage larger than its input voltage, i.e., it is a voltage elevator. In [46, 48], there is a proposition of using this boost circuit as a way to convert DC voltage into an oscillating voltage. However, alternating current cannot be generated with this converter, since the output current cannot change its sign. For this, an inverter is yielded by duplicating the boost circuit [25].

Chapter 2. Introduction

13

PHOTOVOLTAIC ARRAY DC 48V

BOOST INVERTER

AC 220Vrms 50Hz ELECTRICAL GRID AC 220Vrms 50Hz

Figure 2.1: Domestic photovoltaic installation.

The boost inverter is made up of two DC-DC converters2 and a load connected differentially across them, having a bidirectional current (see Fig. 2.2). Each converter produces a DC-biased sine wave output, v1 and v2 , so that each source generates a unipolar voltage. Voltages v1 and v2 should present a phase shift equal to 180◦ , to maximize the voltage excursion across the load. In this way, to generate an oscillatory signal without bias is possible. The circuit implementation is shown in Fig. 2.3. In order to simplify the analysis, a part of the boost inverter is replaced by a constant voltage source as is shown in Fig. 2.4. Once the desired results are obtained, they are extrapolated to the full inverter. Note that, this replacement shows more clearly the bidirectional current of each boost DC-DC converter.

2.1.1 System description Now, some assumptions about the boost inverter are presented. Assume that: 2 Throughout

this part of the thesis, each part of the boost inverter will be referred as ‘boost DC-DC converter’ since each part is a normal boost converter that is commonly used as a DC-DC converter. Nevertheless, it should be taken into account that, in the boost inverter, each part does not act as a DC-DC converter.

2.1. Boost inverter

14

load

+ converter A

+

v1

converter B

v2





Figure 2.2: Basic representation of the boost inverter.

+ Vo Q2

R

+ v1



L1

Q4 L2

C1 −

C2 Vin

Q1

+ v2 −

Q3

Figure 2.3: Boost inverter model.

+ Vo Q2



R +

L1 C1

v1 −

Vin

+ v2 −

Q1

Figure 2.4: Boost inverter model with replacement of a voltage source.

Chapter 2. Introduction

15

• all the components are ideal and the currents of the converter are continuous, • the power supply is constant and known, • the converter operates at a high-switching frequency, • the inductances L1 = L2 , and the capacitances C1 = C2 , are known and symmetric, • v1 and v2 are positive and sinusoidal voltages. The circuit shown in Fig. 2.4 is driven by the transistor ON/OFF inputs, Qi . This yields two modes of operation illustrated in Fig. 2.5. Formally, this circuit generates a switched model. For control purposes, it is common to use an average model described in terms of the mean current and voltage levels [97]. This averaging process may reach an averaged, smooth, nonlinear, continuous-time ordinary differential equation (ODE), as will be seen below. + Vo

L1 iL1

+ v1 −

Vin

R C1

+ Vo

L1 iL1

+ v1 −

Vin





R C1

Q1 + Q2 v2 −

ON OFF

Q1 + Q v2 2 −

OFF ON

Figure 2.5: Operation modes.

If the control variable, q, is defined as q = 0 when Q1 = OFF and Q2 = ON, and q = 1, when Q1 = ON and Q2 = OFF, the converter dynamic equations are diL1 = −v1 + qv1 +Vin dt dv1 v1 v2 C1 = iL1 − qiL1 − + . dt R R

L1

(2.1) (2.2)

Now, u1 = 1 − q is taken as the control action in equations (2.1)–(2.2); obtaining diL1 = −u1 v1 +Vin dt dv1 v1 v2 C1 = u1 i L 1 − + , dt R R

L1

(2.3) (2.4)

2.2. Control problem objectives

16

where u1 is the control variable, which can only take two values u1 ∈ {0, 1}. However, 1 R T +t it is usual to consider its average value u1 (t) = T t u1 (s)ds where T is the switching period [97]. Therefore, u1 is a continuous variable defined as u1 ∈ [0, 1]. The full inverter structure according to Fig. 2.3 is diL1 = −u1 v1 +Vin (2.5) dt dv1 v1 v2 C1 = u1 i L 1 − + (2.6) dt R R diL L2 2 = −u2 v2 +Vin (2.7) dt v1 v2 dv2 C2 = u2 i L 2 + − , (2.8) dt R R where u2 controls the other part of the full system (remember that this part has been replaced by a constant voltage source). L1

The main difficulty of system (2.5)–(2.8) copes with its control, because of: • system nonlinearities. The control signals multiply the state variables. This kind of system are more difficult to study [76]. • The linear part of system (2.5)–(2.8) is nonminimum phase because it has poles in the positive semiplane. Therefore, it is not stable in open-loop [131]. • It is 4th order, which is a relatively high order. • The current signal is indirectly controlled. An alternating current signal can be achieved by a suited control-law that controls the voltage output [48]. • A phase shift of 180◦ is not necessarily achieved. In order to reach a right performance with this inverter, both voltages signals must present this phase shift [25]. • Boost inverter is a double oscillator, thus it does not present two equilibrium points but two limit cycles. • The control law variables are saturated because of duty-cycle signals [25]. • Loads in this kind of systems are unknown or/and slowly variable [63].

2.2 Control problem objectives As mentioned before, the main objective for the boost inverter is to generate alternating current. Hence, it can produce an oscillating voltage centered around zero, and thus, it can

Chapter 2. Introduction

17

achieve negative voltages. Its main drawback, however, is its control due to the complexity of the system structure. The general control objectives for the boost inverter, which are common in switching electronics converters, are: 1. to generate a stable output voltage with an amplitude equal to the desired voltage. Furthermore, in certain applications, it is required that the output voltage has a prespecified phase; 2. to ensure the performance for unknown or/and slowly variable loads; 3. in the case that the control law does not guarantee global stability, to study an attraction domain composed of all initial conditions that ensure a convergency to the system right performance. This estimation of the region of attraction is important for the design of the starting phase. In this thesis, these general objectives can be achieved for a particular solution made up of some proposed specific objectives: 1. to design a suitable control law for the duty cycle by using energy shaping, without introducing reference signals; 2. to achieve an anti-synchronization3 between the voltage signals of each side of the circuit; 3. to propose an adaptive control to deal with unknown and/or slowly varying loads and, 4. to estimate an attraction domain for the resulting system. Figure 2.6 shows a block diagram of the solution proposed in this thesis for the boost inverter control problem. Note that the user has to specify the desired amplitude and frequency of the output voltage, as well as that the initial conditions belong to an estimated attraction domain to ensure the system convergency.

2.2.1 Control law The control of switched-mode inverter is usually accomplished by tracking a reference (sinusoidal) signal [20,35,126,150]. The use of this external signal makes the closed-loop system 3 In

this thesis it is said that two sinusoidal signals of the same frequency are anti-synchronized, or in antiphase, when the phase shift between them is equal to 180◦.

2.2. Control problem objectives

18

a ˆ

xˆ2

ADAPTIVE CONTROL

OBSERVER

x¯ V t CONTROL 1

u1

DC-DC BOOST CONVERTER 1

−V

+

V

V PHASE CONTROLLER

t −V

t

LOAD −V

CONTROL 2

DC-DC BOOST CONVERTER 2

− V t −V

Figure 2.6: General control problem.

non-autonomous in such a way that its analysis is more involved than if it were autonomous. In these kind of systems, the control objective can be seen as the generation of a stable limit cycle defined by a given amplitude and frequency. Figure 2.7 shows the control objective, which corresponds to the valley of a certain surface with a ‘Mexican-hat’ shape [115]. If a control law is able to produce such a limit cycle, alternating current will be generated without the need to introduce a time-dependent reference signal. The generation of limit cycles to produce self-oscillations has been successfully applied to electro-mechanical systems [56, 57]. Applications to electronic devices are [16, 58], where a three-phase UPS and a boost converter are controlled using this method. The idea behind [19, 73] is similar but there a sliding mode controller is proposed. In Chapter 3, this approach is applied to a nonlinear boost inverter [11]. Several approaches have been applied to control this topology of inverter. For instance, in [25, 126], sliding mode method is applied and, in [148], the control is based on passivity. Nevertheless, these methods need a reference signal. Therefore, the main contribution in this part of the thesis is to control the boost inverter without using any reference signal, i.e., the system becomes autonomous [24, 150]. Figure 2.8 shows the autonomous structure that replaces the standard feedback control loop. This is a sub-control problem, since both current and voltage signals are controlled by the same control variable.

Chapter 2. Introduction

19

Figure 2.7: Desired energy function: Mexican-hat shape.

Ref erence

+ CONTROL

SYSTEM



Figure 2.8: Autonomous feedback control loop.

Output

2.2. Control problem objectives

20

It is shown that the direct application of the approach proposed in [16] and [58] does not fulfill the objective due to the lack of anti-synchronization between both parts of the circuit. In order to achieve anti-synchronization, a phase controller (PHC) in an external loop is added to the previous control law. This approach is also applied to synchronize4 the output with a given signal. An example of such a configuration is the synchronization of the boost inverter with the electrical grid (as in the photovoltaic case) in order to achieve a satisfactory power factor, which is shown in Fig. 2.6. The circuit performance is validated in the simulation of a practical case presented in this chapter. Previous results were extended to the case that the load is not purely resistive but it is inductive, as is usually the case in industrial applications. In [10], the extension to the controller based on energy shaping method considering an inductive load was presented, taking also into account the PHC.

2.2.2 Adaptive control Previously, a control law satisfying previous requirements was designed in Chapter 3. It has been supposed that the load is known and constant. However, it is well known in industrial applications that the load can be unknown or suffer perturbations. This problem in switched-mode converters is usually dealt with by using adaptation mechanisms along with other techniques such as feedback stabilization [63], input-output feedback linearization [64], backstepping [123, 135], grid-point modeled [102], sliding modes [28, 134, 145], predictive control [94] or fuzzy logic control [40]. In [108], an adaptive control is obtained for a part of the DC-DC boost converter, which is controlled using the oscillation generation approach mentioned before. This adaptive controller is computed using passivity arguments. This approach is not easily applicable to the boost inverter because its model is more involved than the converter of [108]. In Chapter 4, the goal is to design a load-adaptation mechanism for the boost inverter controlled by energy shaping methodology. In order to estimate the load, a state observer is designed for any system variable even when the state variables are measured (Fig. 2.6). This provides a fast, successful adaptation of the load parameter [7]. This approach is applied by simulation to a real industrial case. The stability of the full system is analyzed by singular perturbation analysis, [76,78]. For the sake of simplicity, the phase controller is not considered in this analysis. The resulting adaptive control is tested by simulations. 4 In

this thesis it is said that two sinusoidal signals of the same frequency are synchronized when the phase shift between them is equal to 0◦ .

Chapter 2. Introduction

21

The extension to case of unknown and non-purely resistive load using the adaptation mechanism before is published in [12]. In this application two parameters of the load are adapted at the same time.

2.2.3 Attraction domain In Chapter 3, a control law for the boost inverter is designed satisfying its main objectives. Ideally, the designed control law guarantees global stability by means of a Lyapunov function. In practice, however, the control law does not achieve global stability due to two reasons: firstly, the ideal control signal cannot be implemented globally due to control signal saturation; and secondly, the circuit imposes physical constraints on certain state variables: the capacitor voltages, for example, cannot be negative. Consequently, it is necessary to estimate an attraction region for the boost inverter. This attraction region is composed of all initial conditions of the system that guarantee the convergency to the right behavior. There is a starting phase, which is very common in this kind of systems [13, 91, 156], that must bring the state of the system into a point inside this region. The attraction domain estimation problem presents several difficulties. The main drawback is the complexity of the control law, which is a rational function with a high degree polynomial numerator. Moreover, it is necessary to highlight that the desired behavior does not correspond to an equilibrium point but to a limit cycle. Therefore, to obtain an estimated attraction domain for the boost inverter can be quite involved. There exist many published methods to estimate the region of attraction (see, for example [54, 76] and the references therein). One example of this kind of methods is based on Lyapunov theory, in which closed Lyapunov-function level surfaces are employed to determine approximate sort of ‘conservative’ estimations for the region of attraction [76] (see Fig. 2.9). These methods often employ polynomial systems [85, 124, 146]. Chapter 5 presents a method of estimating an attraction domain, considering state and control-signal constraints. This approach can be applied to a class of system, whose local system stability was previously guaranteed by a Lyapunov function, as is the problem proposed here. By means of employing this Lyapunov function to estimate a ‘conservative’ attraction domain, a simple computational approach can be generated, although the model and/or control law have a relative high degree and complexity. In order to apply the method, the closed-loop system must be in a polynomial form, in such a way that the problem is transformed in a sum of squares (SOS) optimization problem [117]. This method is applied to estimate an attraction domain for the boost inverter. It is remarked, that the computed attraction region obtained from this Lyapunov function considers physical system constraints, containing control law saturations. The application of this method is very simple and satisfactory results are obtained.

2.2. Control problem objectives

22

ESTIMATED ATTRACTION Attraction DOMAIN

ATTRACTION DOMAIN

Domain

SATURATIONS DESIRED LIMIT CYCLE

CONSTRAINT

Figure 2.9: System attraction domain with constraints.

Chapter 3 Control of the DC-AC boost converter by energy shaping This chapter exposes a novel control strategy for nonlinear boost inverter. Boost inverter is nonminimum phase 4th order nonlinear system, which has not an equilibrium point but a limit cycle. The control objective is not only to obtain a right system performance, but also to guaranty the system stability. In addition, it is necessary to mention that it is a sub-control problem. The control law has to control voltage as well as current signal. The idea behind is based on generating an autonomous stable oscillator. The interesting advantage of this method is that an external reference signal is not needed. This aim is achieved by using energy-shaping methodology with a suitable Hamiltonian function which defines the desired system behavior [44]. This approach guaranties the system stability. The only missed thing in the developed controller is to synchronize the voltage signals with a phase shift of 180◦ . This is important in order to obtain the desired response. For this, a phase controller is added to the control law in order to achieve 180◦ -synchronization between both parts of the circuit, as is shown in Fig. 3.1. In addition, this idea is used to synchronize the voltage output with a pre-specified signal, e.g. synchronization with the electrical grid. The resulting control is tested by means of simulations.

3.1 Normalized average model Assume system (2.3)–(2.4) is only subject to a resistive load. In order to simplify the control study, a known change of variable is employed [18, 133], in order to achieve a normalize 23

3.1. Normalized average model

24 V

t CONTROL 1

u1

−V DC-DC BOOST CONVERTER 1

x2 +



V

Vin

V

ω ˜

t

PHASE CONTROLLER

−V

Vout t

LOAD −V

CONTROL 2

u2

DC-DC BOOST CONVERTER 2

− V

x4 t

−V

Figure 3.1: Controlled boost inverter with PHC.

model: r 1 L1 x1 = iL Vin C1 1 v1 x2 = Vin r 1 L1 x3 = iL Vin C1 2 v2 x4 = Vin

(3.1) (3.2) (3.3) (3.4)

where x1 and x3 are the averaged currents and x2 and x4 are the averaged voltages. The normalized time scale is t˜ = ωnt (3.5) with

ωn = √

1 , L1C1

(3.6)

which yields x˙1 = −u1 x2 + 1 x˙2 = u1 x1 − ax2 + ax4 , x˙3 = −u2 x4 + 1 x˙4 = u2 x3 + ax2 − ax4 , where a =

1 R

q

L1 C1 .

Note that ωn is the natural frequency and a is twice the damping.

(3.7) (3.8) (3.9) (3.10)

Chapter 3. Control of the DC-AC boost converter by energy shaping

25

Remark 3.1 x˙1 , x˙2 , x˙3 and x˙4 are time derivatives of x1 , x2 , x3 and x4 , respectively, with respect to t˜.

As mentioned in the chapitre before, for simplicity, the simplified boost inverter (2.3)– (2.4) is dealt with and, later, the results are extrapolated to the full system. Focussing on the simplified system. If u1 is eliminated in (3.7)–(3.8), next equation is obtained x1 (1 − x˙1 ) = x2 (x˙2 + ax2 − ax4 ). (3.11) This equation is an implicit equation, which relates the state variables (x2 , x4 ) and their time derivatives and does not depend on the control signal u. Note that x4 can be considered an exogenous input in system (3.7)–(3.8). Equation (3.11) can be understood as the internal dynamic of the system. If x˙1 = 0 and x˙2 = 0 is performed, the equilibrium manifold is x1 = ax2 (x2 − x4 ). In this way, the internal dynamic of system (3.7)–(3.8) given by (3.11) acts as a constraint on the system states. From Eq. (3.11), it is possible to see that given x4 , and only controlling x1 , variable x2 can be indirectly controlled1 . Moreover, the stability of the system is maintained [48].

3.2 Energy shaping control for generation of oscillations 3.2.1 Approach overview The generation of alternating current in electronic converters can be achieved by generating a stable limit cycle without the need to introduce a reference signal. To do this, an oscillatory target system may be defined and by matching its equations and system equations (3.7)–(3.8) a control law can be obtained. In order to define the target system, consider the following energy-like function 1 H0 (η1 , η2 ) = Γ21 (η1 , η2 ), (3.12) 4 where η1 and η2 are state variables and Γ1 (η1 , η2 ) , ω 2 (η1 − η10 )2 + (η2 − η20 )2 − µ . Parameters ω , η10 , η20 and µ > 0 should be chosen so that the closed curve Γ1 = 0 defines the desired behavior. This curve is an ellipse centered at point (η10 , η20 ). A dynamical system can be defined such that this closed curve is its limit set. This can be reached by 1 For

the full system, it is had: x1 (1 − x˙1 ) = x2 (x˙2 + a(x2 − x4 )) and x3 (1 − x˙3 ) = x4 (x˙4 + a(x4 − x2 )). Thus, controlling x1 and x3 , the desired behaviors for x2 and x4 can be obtained.

3.2. Energy shaping control for generation of oscillations

26

adopting H0 as a Hamiltonian function [16, 108], and defining the Hamiltonian dynamical system # " ∂H #   " 1 0 −k η˙ 1 a1 Γ1 ∂ η1 = (3.13) ∂ H0 , η˙ 2 − Γ11 −ka2 ∂η 2

which, after using (3.12), results in

η˙ 1 = (η2 − η20 ) − ka1 ω 2 (η1 − η10 )Γ1 η˙ 2 = −ω 2 (η1 − η10 ) − ka2 (η2 − η20 )Γ1 .

(3.14) (3.15)

 H˙0 = −Γ21 ka1 ω 4 (η1 − η10 )2 + ka2 (η2 − η20 )2 ≤ 0,

(3.16)

Taking into account that

by using the LaSalle invariance principle it can be seen that, for all initial conditions except the center of the ellipse, the trajectories of the system tend to the curve Γ1 = 0. Figure 3.2 shows this energy-like function.

Figure 3.2: Desired energy function: mexican-hat.

The behavior of the target system (Γ1 = 0) corresponds to the desired sinusoidal behavior for the DC-AC converter. Constants ω , η10 , η20 and µ are design parameters for the frequency, bias and amplitude of the desired behavior, while ka1 and ka2 define the speed of the transient response. Note that η˙ 1 and η˙ 2 are in this case time derivatives of η1 and η2 with respect to t˜, in order to work with the normalized averaged model (3.7)–(3.8).

Chapter 3. Control of the DC-AC boost converter by energy shaping

27

3.2.2 Controller design System (3.7)–(3.8) can not be directly transformed to the form of system (3.14)–(3.15), but this can be done using the new change of coordinates given below: x21 + x22 2 ζ2 = x1 − ax22 + ax2 x4 + ζ20

ζ1 =

(3.17) (3.18)

where ζ20 is an offset term that will be a tuning parameter. From (3.17)–(3.18), it is easy to see that

ζ˙1 = ζ2 − ζ20 ζ˙2 = 1 + 2a2 x22 − 3a2 x4 x2 + a2 x24 + ax2 x˙4 − u1 (x2 + 2ax1 x2 − ax4 x1 ).

(3.19) (3.20)

It is not easy to obtain simple relationships x1 = f (ζ1 , ζ2 ) and x2 = f (ζ3 , ζ4 ) from (3.17)– (3.18) due to the quadratic terms. Nevertheless, this change of variables is a diffeomorphism if and only if x2 + 2ax1 x2 − ax4 x1 6= 0, as follows from the inverse function theorem. In Chapter 5 it will be seen that this constraint restricts the domain of attraction of the desired limit cycle when the controller obtained below is applied. Looking at target system structure (3.14)–(3.15) and comparing it with (3.19)–(3.20) the choice ka1 = 0 is obvious, resulting in the target system

ζ˙1 = ζ2 − ζ20 ζ˙2 = −ω 2 (ζ1 − ζ10 ) − k1 Γ1 (ζ2 − ζ20 ),

(3.21) (3.22)

where, for sake of simplicity, ka2 has been denoted as k1 . The attraction of curve Γ = 0 can still be proved by the LaSalle invariance principle. The control law, u, that matches (3.19)–(3.20) and (3.21)–(3.22) is u1 =

1 + 2a2 x22 − 3a2 x4 x2 + a2 x24 + ax2 x˙4 + k1 Γ1 (ζ2 − ζ20 ) + ω 2 (ζ1 − ζ10 ) . x2 + 2ax1 x2 − ax4 x1

(3.23)

Indeed, u1 varies dependently on x, as can be noted from Eq. (3.17)– (3.18). This controller has several problems. First, the denominator in (3.23) may be zero (this is the same necessary condition for (3.19)–(3.20) to be a diffeomorphism). Furthermore, in other cases, the resulting u1 can violate the constraint 0 ≤ u1 ≤ 1. In Chapter 5, an estimation for the region of attraction of the desired limit cycle will be obtained by taking these problems into account. It is assumed that a starting strategy will bring the state of the system into this region of attraction [13, 91, 156].

28

3.2. Energy shaping control for generation of oscillations

Parameters η10 , η20 and µ have to be defined as a function of the desired behavior. For this, it is necessary to obtain an analytical expression of the desired objective curve in plane x1 − x2 . Assume that the desired time evolutions for x2 and x4 are x∗2 = A sin ω t + B x∗4 = −A sin ω t + B,

(3.24) (3.25)

x∗1 = aα0 + α1 cos ω t + β1 sin ω t

(3.26)

where A, B and ω take pre-specified values to obtain the desired evolution for v1 and iL1 using (3.1)–(3.4), (3.5) and (3.6). In addition, note that these desired evolutions allow us to remove the bias in the output. The origin of time in (3.24)–(3.25) is arbitrary in such a way that no phase shift value is imposed (signal synchronization will be achieved below). Assume that the desired steady state for x1 can be approximated by

This assumption is very common in the field of electronics [20, 37, 48, 60]. By substituting (3.24)–(3.25) and (3.26) in (3.11) 1 aα0 + (β1 + aα0 α1 ω ) sin ω t + (α1 − aα0 β1 ω ) cos ω t + ω (α12 − β12 ) sin 2ω t 2 1 −α1 β1 ω cos 2ω t = aA2 − 2aAB sin ω t − ω AB cos ω t + ω A2 sin 2ω t − aA2 cos 2ω t. 2 If the second order harmonics are neglected, the corresponding coefficients can be equated: aα0 = aA2 β1 + aα0 α1 ω = −2aAB α1 − aα0 β1 ω = −ω AB.

When this system is resolved for α0 , α1 and β1 ,

α0 = A2 ω AB(2a2 A2 + 1) α1 = 1 + a 2 A4 ω 2 aAB(ω 2 A2 − 2) β1 = − . 1 + a 2 A4 ω 2

(3.27) (3.28) (3.29)

The next problem is to show that the desired behavior for ζ1 and ζ2 is an ellipse and defining the ellipse parameters (ω , ζ10 , ζ20 and µ ) in terms of the desired behavior for x2 . For this, it is necessary to obtain the desired evolution for ζ1 and ζ2 by applying the change of variables (3.17)–(3.18) to (3.24)–(3.25) and (3.26). 1 ζ1 = [(aα0 + α1 cos ω t + β1 sin ω t)2 + (A sin ω t + B)2 ] (3.30) 2 ζ2 = aα0 + α1 cos ω t + β1 sin ω t − a(A sin ω t + B)2 + a(−A2 sin2 ω t + B2 ) + ζ20(3.31)

Chapter 3. Control of the DC-AC boost converter by energy shaping

29

Expanding these expressions in Fourier terms yields (0)

(11)

cos ω t + ζ1

(0)

(11)

cos ω t + ζ2

ζ1 = ζ1 + ζ1 ζ2 = ζ2 + ζ2

(12)

sin ω t + ζ1

(21)

cos 2ω t + ζ1

(12)

sin ω t + ζ2

(22)

sin 2ω t

(3.32)

(21)

cos 2ω t + ζ2

(22)

sin 2ω t.

(3.33)

By equating (3.30)–(3.31) with (3.32)–(3.33) the following Fourier coefficients, are obtained 2a2 α02 + α12 + β12 + A2 + 2B2 4 (11) ζ1 = aα0 α1 (0)

ζ1 = (12)

ζ1

= aα0 β1 + AB

α12 − β12 − A2 4 α β (22) 1 1 ζ1 = 2 (0) ζ2 = ζ20 (21)

=

(11)

= α1

(12)

= β1 − 2aAB

ζ1

ζ2 ζ2

(21)

ζ2

(22)

ζ2

= aA2 = 0.

(21)

(22)

(21)

(22)

Assuming that the double frequency terms ζ1 , ζ1 , ζ2 and ζ2 can be neglected, these expressions can be approximated by an ellipse in the plane ζ1 , ζ2 since (3.30)–(3.31) yields (11)

ωζ1

(12)

ωζ1

(12)

= −ζ2

(11)

= ζ2

.

(3.34) (3.35)

The parameters of this ellipse are given by (0)

ζ10 = ζ1

(0)

ζ20 = ζ2 µ

(11) (12) = ω 2 ((ζ1 )2 + (ζ1 )2 ).

(3.36) (3.37) (3.38)

3.2.3 Control law for the full system The boost inverter is composed of two DC-DC converters. Therefore, it has two control signals. For that, the previous control law, which is used for a part of the system, is used in

3.2. Energy shaping control for generation of oscillations

30

order to obtain the control law for the other part of the system. By comparing the normalized full model Eqs. (3.7)–(3.10) and (3.7)–(3.8) there is a similar structure for the pairs of current and voltage of both boost DC-DC converters. Therefore, the two control laws are easily obtained. Control law u2 is obtained by using symmetry. The control laws are 1 + a2 (2x22 − 3x2 x4 + x24 + x2 x˙4 ) + k1 Γ1 (ζ2 − ζ20 ) + ω 2 (ζ1 − ζ10 ) , κ1 (x) (3.39) x2 + 2ax1 x2 − ax4 x1 1 + a2 (2x24 − 3x2 x4 + x22 + x4 x˙2 ) + k2 Γ2 (ζ4 − ζ40 ) + ω 2 (ζ3 − ζ30 ) u2 = , κ2 (x),(3.40) x4 + 2ax3 x4 − ax2 x3 u1 =

where Γ1 (ζ1 , ζ2 ) = ω 2 (ζ1 − ζ10 )2 + (ζ2 − ζ20 )2 − µ Γ2 (ζ3 , ζ4 ) = ω 2 (ζ3 − ζ30 )2 + (ζ4 − ζ40 )2 − µ .

(3.41) (3.42)

The expressions for the time derivatives x˙2 and x˙4 are taken directly from the normalized equations of the boost inverter. The stability is proved taking: 1 H = (Γ21 + Γ22 ), 4 whose differetation is: H˙ = −Γ21 k(ζ2 − ζ20 )2 − Γ22 k(ζ4 − ζ40 )2 ≤ 0.

3.2.4 Simulation results The following simulation shows how the controller is applied in a practical case. It is desired to obtain an output voltage Vo = 220 √2 sin(50 · 2π t) from an input voltage Vin = 48V . 2

These simulations are performed considering, R = 100Ω, L = 250µ H, C = 250µ F. The desired frequency and voltage amplitude are f = T1 = 50Hz and 220Vrms, respectively.

a 2

Note that, ωn = 4 · 103 rad s = 0.01 ≪ 1 .

and

ϖ = 2π f = 3.14 · 102 rad s ,

i.e.,

ωn > 5ϖ

and

Parameter A in Eqs. (3.24)–(3.25) has to be the half of the desired output voltage amplitude, and B is chosen so that x2 and x4 are always positive. In order to obtain this voltage, the parameters are A = 3.33 and B = 9.37 with ω = 0.078 in the normalized variables (x1 , x2 ).

Chapter 3. Control of the DC-AC boost converter by energy shaping

31

Figure 3.3 shows the results of a simulation using a commutation frequency of 50KHz and employing a sample time of 0.1T s (remind that T is the commutation frequency period). Both DC-DC converters achieve the desired limit cycle.

a)

b)

14

3 2

12

1 x2 , x4

ζ2 , ζ4

10 0

8 −1 6

4 −5

−2

0 x1 , x3

5

−3 0

50

100

ζ1 , ζ3

Figure 3.3: a) Evolution of (x1 , x2 ) (solid) and (x3 , x4 ) (dashed); b) evolution of (ζ1 , ζ2 ) (solid) and (ζ3 , ζ4 ) (dashed).

Figure 3.4 shows the boost inverter output voltage. Note that the system does not show overshoot, converging to the desired behaviour very fast. It can be seen that the desired amplitude is not achieved. The reason is that the previous design does not force the phase shift between signals v1 and v2 to be in anti-phase (180◦ phase shift). Figure 3.5 shows that, as a result, this goal is not achieved. The next section deals with this problem. The control signals are shown in Fig. 3.6. For the parameter chosen in the application, the control law signal oscillates between 0.08 and 0.17. As pointed out by the jury member this signal has values smaller than 0.1, what is not good for the implementation. This is a circuit design problem. If the circuit is designed for more suitable values of the duty cycle, the proposed controller would also lead to satisfactory results.

3.2. Energy shaping control for generation of oscillations

32

300 200

Vout(V)

100 0 −100 −200 −300 0

0.02

0.04 0.06 Time(s)

0.08

0.1

Figure 3.4: Output voltage of the boost inverter.

600 550

V1V2(V)

500 450 400 350 300 0

0.02

0.04 0.06 Time(s)

0.08

0.1

Figure 3.5: Output voltages of the first (solid) and second (dashed) boost DC-DC converters.

Chapter 3. Control of the DC-AC boost converter by energy shaping

33

1

u1 , u2

0.8

0.6

0.4

0.2

0 0

0.02

0.04

0.06

0.08

0.1

Time(sec) Figure 3.6: Time evolution of the control laws.

3.3 Synchronization problem The controllers developed above for boost inverters do not synchronize the two parts of the circuit with a phase shift of 180◦ since each one controls independently a DC-DC converter. Therefore, in the above design, the voltage signal did not present the phase shift mentioned before. In order to get the desired output voltage, it is necessary to synchronize these signals, in such a way that they present a phase shift to 180◦ . In this section, a phase controller (PHC), inspired by the configuration of a phase-lock loop (PLL) [66]– [1], is added. The PHC allows us to achieve the desired phase shift between the output of the two DC-DC converters as well as to synchronize the boost inverter output with respect to a specified voltage signal, as in the case of synchronization with the electrical grid.

3.3.1 Boost inverter synchronization The objective is to synchronize voltage signals x2 and x4 in anti-phase. The method is illustrated in Fig. 3.7. The normalized voltage of the second DC-DC converter, x4 , is taken as a reference signal and the normalized voltage of the first DC-DC converter, x2 , is the signal to be synchronized with x4 in anti-phase. These are the inputs to the PHC. The output is a frequency variation, ω˜ , which is added to the nominal frequency, ω , in the Control 1 block and the resultant frequency is entered in (3.39). The output of the converter is a sinusoidal signal of that resulting frequency.

3.3. Synchronization problem

34

PHC x2

ω˜ Control 1

Converter 1

x4

R

Control 2

Converter 2

u1

u2

ω Figure 3.7: Block diagram of boost DC-AC converter with output voltages in anti-phase by the PHC.

The PHC block diagram appears in Fig. 3.8. The multiplier obtains the product, x′2 × x′4 , in such way that its output, once filtered by a low pass filter (LPF), is a measure of the deviation of the phase shift with respect to 90◦ , [66]. For this reason, one of the inputs of the multiplier, e.g. x′2 , is obtained by passing voltage x2 through a high pass filter (HPF) in order to eliminate its continuous component. Likewise, x′4 is obtained after passing x4 through another HPF, changing its sign and integrating it. In the following, an intuitive explanation of the correct behavior of the full system is presented. It is assumed that ω˜ (t) is small and varies slowly. Under this assumption it can be expected that the introduction of the PHC does not affect the normal behavior of the controllers (3.39) and (3.40) (apart from the phase shift between x2 and x4 as desired). In this way, it can be assumed that, after a transient period, x2 and x4 evolve as sinusoidal signals: x2 = A sin((ω + ω˜ (t))t + ϕi′ ) + B ≈ A sin(ω t + ϕi (t)) + B x4 = A sin(ω t + π ) + B

(3.43) (3.44)

where the origin of time has been chosen in such a way that ϕi (t) represents the phase shift between x2 and the desired behavior for x2 . Note, ϕi′ is constant. Likewise, it is desired that ϕi = 2nπ with n ∈ Z. Assuming that both HPFs eliminate the bias terms, the signals that enter into the multiplier in Fig. 3.8 are x′2 ≈ A sin((ω t + ϕi (t)) x′4 = A cos(ω t + π ) +CPHC where CPHC is generated by the integrator, being eliminated by the LPF. Assuming that the LPF filters out every sinusoidal signal of frequency greater or equal than ω , then A2 Kvd A2 Kvd ω˜ (t) = sin(ϕi (t) − π ) = − sin(ϕi (t)). 2 2

Chapter 3. Control of the DC-AC boost converter by energy shaping

x2

x′2

HPF

× x4

HPF

35

ω s

−1

Kvd

LPF

ω˜

x′4

Figure 3.8: Conceptual block diagram of the PHC.

where Kvd is a positive design parameter. In this expression, sin(α ) cos(β ) = 12 (sin(α + β ) + sin(α − β )) has been employed. On the other side, from Eq. (3.43), it is easy to obtain ϕ˙ i ≈ ω˜ , and thus, A2 Kvd ϕ˙ i ≈ − sin(ϕi (t)). 2 Obviously ϕi (t) converges to 2nπ with n ∈ Z, which corresponds to the desired behavior. Simulation results The previous values are used for the boost inverter parameters. The high pass filter applied is: 1.4s . s+ω The LPF is a second order Butterworth filter [69]: √ (s + 0.008( 22



√ 2 2

1 j))(s + 0.008(

√ 2 2

+

√ 2 2

.

(3.45)

j))

The value of the PHC gain is Kvd = 5 · 10−4. The results of the PHC application are shown in Fig. 3.9. Voltages v1 and v2 in anti-phase. In Fig. 3.10 the boost inverter output voltage is represented. Figure 3.11 shows the ripple in the inductance currents, which is quite acceptable. The output signal has a total harmonic distortion (THD) below 0.22% for a 50-Hz output voltage. Figure 3.12 shows the signal spectrum of the signals v2 and v4 . As can be seen, the harmonics of the fundamental frequency wave of the obtained output is quite satisfactory. This result justifies the first harmonic approximations carried out during the design of the control law. Of course, this is only valid for the chosen parameters and it does not prove the usefulness of the law in a general sense. In fact, the approximation does not work when the circuit parameters are not chosen adequately, but it is thought that, when

3.3. Synchronization problem

36

b) 600

550

550

500

500 V1V2(V)

V1V2(V)

a) 600

450

450

400

400

350

350

300 0

300 0.02

0.04 0.06 Time(s)

0.08

0.1

0.6

0.62

0.64 0.66 Time(s)

0.68

0.7

Figure 3.9: Output voltage of the first (solid) and second (dashed) boost DC-DC converters. In a) the transient time is shown and in b) the steady-state is shown.

b) 300

200

200

100

100 Vout(V)

Vout(V)

a) 300

0

0

−100

−100

−200

−200

−300 0

−300 0.02

0.04 0.06 Time(s)

0.08

0.1

0.6

0.62

0.64 0.66 Time(s)

0.68

0.7

Figure 3.10: Output voltage with PHC. In a) the transient time is shown and in b) the steady-state is shown.

Chapter 3. Control of the DC-AC boost converter by energy shaping

i1

37

i2

100

i(A)

50

0

−50

−100

0.6

0.62

0.64

t(s)

0.66

0.68

0.7

Figure 3.11: Inductance currents.

25

Magnitude (dB)

20

15

10

5

0 0

50

100 150 Frequency (Hz)

200

Figure 3.12: Output voltage signal spectrum.

250

3.3. Synchronization problem

38

the circuit is designed properly (taking into account the voltage, current and power levels), the approximations will yield good results. Notice once more that these approximations are common in the literature [47, 48].

3.3.2 Synchronization with the electrical grid In some applications, such as renewable energy plants, an inverter is necessary to inject energy from production plants into the electrical grid. In this case, the problem is to synchronize the voltage output signals with the pre-specified signal. For this, the normalized voltage signals of both DC-DC converters (x2 , x4 ) are treated with two PHCs, as is shown in Fig. 3.13. x2 is synchronized with the grid voltage using a phase shift of 180◦ by using PHC1, whose structure is shown in Fig. 3.14 and which is similar to the PHC in Fig. 3.8. Signal x4 is synchronized with the grid voltage using  a zero phase shift by means of PHC2 π A ′ shown in Fig. 3.15. In this case, g = ω sin ω t − 2 .

For simplicity, the stability analysis of the full system is not delivered, which is similar to the previous subsection. Grid PHC1

PHC2

ω˜ 1

x2

Control 1

Converter 1

ω˜ 2

x4 R

Control 2

Converter 2

u1

u2

ω Figure 3.13: Block diagram of boost DC-AC converter with output voltage synchronized with the electrical grid.

Simulation results

The grid voltage is 2 Vgrid = 220 √ sin(100π t) 2

(3.46)

The values used in this simulation for the filter parameters and gain, Kvd , are the same ones used previously.

Chapter 3. Control of the DC-AC boost converter by energy shaping

x′2

HPF

x2

× g

HPF

39

−1

1 s

LPF

Kvd

ω˜

g′

Figure 3.14: Conceptual block diagram of PHC1.

x′4

HPF

x4

× g

HPF

+1

1 s

LPF

Kvd

ω˜

g′

Figure 3.15: Conceptual block diagram of PHC2.

The performance of the synchronization of the boost DC-AC converter with the electrical supply voltage is represented in Fig.3.16 showing satisfactory behavior.

3.4 Conclusion A control strategy for the complex nonlinear boost inverter was presented. The method is based on energy-shaping methodology, which generates a limit cycle guaranteeing the system stability. The obtained control law is a complex expression. Nevertheless, it has an important relevance: the system does not require any external reference signals. The resulting controller achieves the objective by adding a phase controller. The same idea is used in order to solve the problem of grid electrical synchronization. The control laws designed in this chapter depends on the value of the resistive load, which is not necessarily known. Next chapter will deal with unknown and/or slowly varying loads, which are supposed to be constant.

3.4. Conclusion

40

300 200

Vo (V )

100 0 −100 −200 −300 0

0.05

0.1

0.15

0.2

t(s) Figure 3.16: Electrical supply voltage (solid) and simulated output voltage synchronized with PHCs (dashed).

Chapter 4 Adaptive control In the previous chapter, a control law for the nonlinear boost inverter has been designed, which guaranties the system stability. Likewise, a phase controller has been proposed in order to achieve the desired phase in the output voltage signal. This development has been performed assuming known load. Nevertheless, it is usual, that the load is unknown or/and change slowly. In this chapter an adaptive control is designed for the boost inverter in order to cope with unknown and/or varying resistive load (see Fig. 4.1). This is a common problem in the field of electronics. Different control strategies have been applied to provide a solution to this standard problem in switched-mode converters [28, 63, 64, 94, 135]. Adaptation mechanism for similar controllers, as the one presented in previous chapter, were used in [108] for the case of the boost converter. The fact that the boost converter model is a 4th -order system makes the design of the adaptation law more involved. A state observer for some of the converter variables is designed even when the state variables are measured. In order to analyze the stability of the full system singular perturbation analysis is used [76]. For simplicity, the phase-lock system is not considered in this analysis. The resultant adaptive control is tested by means of simulations.

4.1 Design of an adaptive control An adaptive law (or a load observer) is proposed to cope with load variations and/or load uncertainties. This observer is designed based only on a one-sided circuit, which contains enough information to make this parameter observable. Therefore, the study of the full twosided circuit is avoided due to symmetry considerations. 41

4.1. Design of an adaptive control

42

a ˆ

xˆ2

ADAPTIVE CONTROL

OBSERVER

x¯ V t CONTROL 1

−V

u1

DC-DC BOOST CONVERTER 1

x2 +



V

Vin

ω ˜

t −V

V PHASE CONTROLLER

Vout t

LOAD −V

CONTROL 2

u2

DC-DC BOOST CONVERTER 2

− V

x4 t

−V

Figure 4.1: Controlled boost inverter with observer.

The model problem for one-sided circuit (left part of the Fig. 2.3) (3.7)-(3.8), can be rewritten compactly as: x˙l = Ul xl + aDl y + El a˙ = 0 y = x2 − x4 ym = Mxl

(4.1) (4.2) (4.3) (4.4)

with xl = [x1 , x2 ]T ; x4 can be considered as an input, and 

     0 −u1 0 1 Ul = , Dl = , El = , M = I2×2 . u1 0 −1 0 Note, that y ∈ R1 and ym ∈ R2 .

Remark 4.1 In what follows, it is assumed that both voltage and current, ym , are measurable, and thus accessible for control use.

Remark 4.2 From remark 4.1, note that xl and y are measurable variables.

Chapter 4. Adaptive control

43

4.1.1 Adaptation law The proposed adaptation law is comprised of a state observer for one side of the inverter boost, plus an adaptation law for parameter a. From Eq. (4.4): K0 (ym − yˆm ) = K0 (Mxl − M xˆl ) = K(xl − xˆl ) where K0 , K ∈ R2×2 are constant design matrix. Therefore, the adaptation law has the following structure: xˆ˙l = Ul xl + aD ˆ l y + El + K(xl − xˆl ) a˙ˆ = β (xl , xˆl ),

(4.5) (4.6)

where β (xl , xˆl ) is the adaptation law to be designed, xˆl is the estimated state of xl and aˆ is the estimated value of a. From remark 4.2, the real state of xl and y in Eq. (4.5) can be used. Note that even if xl is accessible, the adaptation law designed here requires the additional (or extended) state observer. This will become clear during the analysis of the error equation system, as will be shown below.

4.1.2 Error equation Assume that a is a constant parameter (a˙ = 0) or that it changes slowly (a˙ ≈ 0) and define the following error variables: x˜l = xl − xˆl ,

a˜ = a − a, ˆ

˙ˆ a˙˜ = −a.

Error equations are now derived from (4.1)–(4.4) together with (4.5)–(4.6) x˙˜l = −K x˜l + aD ˜ ly a˙˜ = −β (xl , xˆl ).

(4.7) (4.8)

Let K be of the form, K , α I,

α >0

and P = I be the trivial solution of PK T + KP = −Q, with Q = 2α I. Now, introducing V = x˜Tl Px˜l +

a˜2 γ

(4.9)

4.1. Design of an adaptive control

44

it follows that   a˙˜ T T ˙ V = −x˜l Qx˜l + 2a˜ x˜l PDl y + γ   a˙ˆ T T = −x˜l Qx˜l + 2a˜ x˜l PDl y − . γ The adaptation law is now designed by canceling the terms in parentheses, i.e. a˙ˆ = γ (DTl Px˜l )y.

(4.10)

4.1.3 Stability properties The observer and the adaptive law error equations are now fully defined. These equations are: x˙˜l = −K x˜l + aD ˜ ly T a˙˜ = −γ (Dl Px˜l )y.

(4.11) (4.12)

The stability properties of these equations follow from the Lyapunov function, V , defined above. Note that with choice (4.10), it follows V˙ = −x˜Tl Qx˜l From standard Lyapunov arguments [76], it is proved that error variables x˜l and a˜ are bounded. Moreover, asymptotic stability is established by LaSalle’s invariance principle [76]. For this, consider the level set Vc = V (x˜l , a, ˜ y) ≤ c0 for sufficiently large c0 > 0 and where V˙ (x˜l , a, ˜ y) ≤ 0. From Eqs. (4.11)–(4.12), note that V˙ (x˜l , a, ˜ y) is negative everywhere, except on the line x˜l = 0. Note that x˜l ≡ 0 that implies x˙˜l ≡ 0, is only obtained if a(t)D ˜ l y(t) ≡ 0.

(4.13)

In addition, observe that if y behaves as a sinusoidal (as is expected from the control problem formulation) the unique asymptotic solution for a˜ is a˜ = 0, as long as y 6≡ 0, ∀t ≥ 0. Consequently, the maximum invariant set in V˙c (x˜l , a, ˜ y) = 0 corresponds to the single point (x˜l = 0, a˜ = 0) Therefore, every solution starting in Vc approaches this point as t → ∞. The following lemma summarizes the above results, assuming that y 6≡ 0. Note, that the manifold y(t) ≡ 0 has to be carefully analyzed since can cause problems (e.g. in Eq. 4.13).

Chapter 4. Adaptive control

45

An intuitive explanation about the right system performance (as will be seen by simulation below) is the PHC introduction, see Chapter 3. The PHC objective is to achieve x2 + x4 = 0, thus avoiding y(t) = x2 − x4 = 0. A rigorous explanation has not been possible to provide it a cause of the system and control nature.

Lemma 4.3 Consider the open-loop system (4.1)–(4.4), and the observer (4.5)–(4.6) with K = α I such that K is a solution for PK T + KP = −Q, then the observer states have the following properties: i) The estimated states xˆl , aˆ are bounded. ii) limt→∞ xˆl (t) = x(t). iii) limt→∞ a(t) ˆ = a, if y(t) 6≡ 0, ∀t ≥ 0.

4.2 Stability considerations of the full closed-loop system In the previous section, the stability properties have been presented for the extended observer. These properties are independent of the evolution of the system state variables. The stability of the complete system is analyzed in this section. The open-loop two-sided inverter (3.7)-(3.8) plus (2.7)-(2.8) normalized, can be compactly rewritten as: x˙ = U (u1, u2 )x + aDy + E y = x2 − x4

(4.14) (4.15)

with x = [x1 , x2 , x3 , x4 ]T , and 

     1 0 −u1 0 0 0  u1 0 0 0   −1  0      U =  0 0 0 −u2  , D =  0  , E =  1 . 0 0 0 u2 0 1

4.2.1 Tuned system The tuned system is defined as the ideal closed-loop system controlled by the tuned feedback laws u∗1 = κ1 (x, a∗ ) and u∗2 = κ2 (x, a∗ ), (from Eq. (3.39) and Eq.(3.40), respectively ), where a∗ is the exact value of a.

4.2. Stability considerations of the full closed-loop system

46

The tuned system given in (3.24)–(3.25) and (3.26) states x˙ = U (u∗1, u∗2 )x + aDy + E = U (κ1(x, a∗ ), κ2 (x, a∗ ))x + aDy + E , f (x)

(4.16) (4.17) (4.18)

and, it achieves an asymptotically orbitally stable periodic solution, i.e. x∗ (t) = x∗ (t + T ). In Section 3, it has been shown that functions Γ1 and Γ2 defined in (3.41)–(3.42) tend to zero. They correspond to periodic sinusoidal solutions of period T = 2π /ω . Consequently, y∗ = x∗2 − x∗4 is also sinusoidal.

4.2.2 Closed-loop system In practice, the control laws that are effectively applied depend on the estimation, a, ˆ of ˆ and uˆ2 = κ2 (x, a), ˆ respectively. parameter a. This control laws are denoted as uˆ1 = κ1 (x, a) Note that these control laws depend on state x and not on their estimations, xˆl , because state x is directly measured. The role of xˆl is to make possible the design of the adaptation law for parameter a. The closed-loop equation resulting from the use of uˆ1 = κ1 (x, a), ˆ uˆ2 = κ2 (x, a), ˆ u∗1 = κ1 (x, a∗ ) and u∗2 = κ2 (x, a∗ ) is written as x˙ = U (uˆ1 , uˆ2 )x + aDy + E +U (u∗1, u∗2 )x −U (u∗1 , u∗2 )x = f (x) + [U (uˆ1, uˆ2 ) −U (u∗1, u∗2 )] x,

(4.19) (4.20)

Let us assume that aˆ ∈ [aˆm , aˆM ] and denote a˜ , a∗ − a. ˆ Applying mean-value theorem [76] yields U (u∗1 , u∗2 ) −U (uˆ1, uˆ2 ) = T (x)a, ˜ being

T (x) ,

"

I ∂∂ua1 |a=aˆ¯ 0

#

0 , I ∂∂ua2 |a=aˆ¯

where, aˆ¯ takes any value belonging to the interval A , [min{a∗ , aˆm }, max{a∗ , aˆM }], and   0 −1 I, . 1 0

Chapter 4. Adaptive control

47

The term T (x)a˜ captures the mismatch between the estimated and the true value of the load. In view of the discussion above, this term has the following property:

Property 4.4 For small enough constants εx > 0 and εa > 0. Let M = {(x, a) ˜ : dist(x, x∗ ) ≤ εx , |a| ˜ ≤ εa } be a compact domain that includes the asymptotic periodic solutions from the tuned system and the exact value of a i.e. a∗ . Then, T (x)a˜ has ∀(x, a) ˜ ∈ M, the following properties:

i) it is continuous, analytic, and free of singularities ii) lima→0 T (x)a˜ = 0. ˜

Putting (4.20) together with the observer error system yields the complete set of closedloop equations, with y = y(x) x˙ = f (x) − T (x)ax ˜ x˙˜l = −α x˜l + aD ˜ ly T a˙˜ = −γ (Dl Px˜l )y,

(4.21) (4.22) (4.23)

remain that K = α I. The stability consideration discussed here will be based on the timescale separation. The main idea is that with the suitable choice of gains (as discussed below), observer equation (4.22)-(4.23) can be seen as the fast subsystem and equation (4.21) as the slow subsystem. Note again that this time-scale separation should be enforced by a particular choice of the observer and adaption constants α and γ , respectively.

4.2.3 Singular perturbation form In order to rewrite the system above in the standard singular perturbation form, it is necessary to follow the next steps:

• introduce a¯ = αa˜ , • select γ = α • define ε =

1 α

48

4.2. Stability considerations of the full closed-loop system

With these considerations, Eqs. (4.21)–(4.23) are rewritten: x˙ = f (x) − T (x)α a( ¯ a)x ˜ ε x˙˜l = −x˜l + a( ¯ a)D ˜ l y, T ε a˙¯ = −(Dl Px˜l )y, where ε > 0 is small parameter (for a larger value of α ). Note that, this particular selection of gains imposes relative gains for the adaptation, γ , and defines precisely, how the observer gain is related to γ . Remark 4.5 The perturbed variable parameter, ε = ς (a), ˜ and for a side effect, the adaptation gain, γ , must fulfill   1 1 , ε ≪ min ω2 k  γ ≫ max ω 4 , k2 . These relationships with respect to the tuning parameter, k, and desired frequency, ω , ensures the convergency of the observer and adaptation parameter, a. Letting z = [x˜l , a( ¯ a)] ˜ T yields the general form x˙ = f (x) − T (x)α a( ¯ a)x ˜ ε z˙ = g(x, z)

(4.24) (4.25)

with x(t0 ) = x0 , x ∈ R4 , z(t0 ) = z0 , z ∈ R3 , and   −x˜l + a( ¯ a)D ˜ ly g(z, x) = −(DTl Px˜l )y According to the singular perturbation analysis, the next steps must be followed:

1. Find a stationary solution of the fast subsystem (4.25) by finding the roots of the equation g(x, z) = 0, i.e. z = φ (x). 2. Replace this solution in the slow subsystem (4.24), and find the resulting slow system x˙ = f (x) − T (x)φ (x)x. 3. Check the boundary layer properties of the fast subsystem along one particular solution of x˙ = f (x) − T (x)φ (x)x.

Chapter 4. Adaptive control

49

4.2.4 Slow sub-system The previous step 1 requires to find the roots of x˜l = a( ¯ a)D ˜ ly 0 = −a( ¯ a)D ˜ Tl PDl y2 . Note that DTl PDl = 2, and that the above equation has multiple solutions, i.e x˜l = 0 a( ¯ a)y ˜ (x) = 0 2

which means that if y ≡ 0, there is one solution for x˜l = 0, and infinite solutions for a. ¯ However, if y 6≡ 0, for instance, in the particular tuned solution y∗ = A cos(ω t), then   x˜l z = φ (x) = =0 a( ¯ a) ˜ becomes an isolated root. Now, step 2 is considered. For the previous particular solution, and taking into account aˆ that a¯ = a− α = 0, i.e. aˆ = a, the slow model is written as: x˙ = f (x) − T (x) · 0 · x = f (x),

(4.26)

which is nothing other than the tuned system whose solutions x(t) = x∗ (t) are sinusoidal.

4.2.5 Boundary layer fast subsystem The next step is to evaluate the stability of the boundary layer system in the finite time interval t ∈ [t0,t1 ]. This is obtained by evaluating the fast subsystem (4.25) along one particular solution of the quasi-steady-state x p (t), and by re-scaling time t to the stretched time coordinates τ = (t − t0 )/ε . The fast subsystem (4.25) evaluated along this trajectory is d xˆ˜l = −xˆ˜l1 dτ 1 d ˆ ˆ¯ a)y x˜l = −x˜ˆl2 − a( ˜ p dτ 2 d aˆ¯ = x˜ˆl2 y p , dτ which can be rewritten as: d zˆ = J(y p )ˆz = J(τ , ω , ε )ˆz. dτ

(4.27)

4.2. Stability considerations of the full closed-loop system

50

with



 −1 0 0 J =  0 −1 −y p  0 yp 0

(4.28)

Under these conditions, system (4.27) is reduced to the autonomous linear system d zˆ = J(τ , ω , 0)ˆz = J(y p 0 )ˆz. dτ

(4.29)

Consider the y p0 ∈ Dx , with Dx , {x : |y| = |x2 − x4 | > ε0 }. The above system has the following properties. Property 4.6 The eigenvalues of J(y p ), for [t, x p, z] ∈ [t0 ,t1] × Dx × R3 , are all strictly negative, i.e.

λ1 = −1 (

(4.30) p

)

1 − 4y p 2 <0 2 ) ( p −1 − 1 − 4y p 2 λ3 = Re < 0, 2

λ2 = Re

−1 +

(4.31) (4.32)

where ε0 > 0 is a constant. Therefore J(y p ) is Hurwitz in the considered domain. As a result, there exists a matrix P(y p ) = P(y p )T > 0 and a Q(y p ) > 0 such that the standard Lyapunov equation holds: P(y p )J(y p) + J(y p )T P(y p ) = −Q(y p ). From standard Lyapunov arguments, it follows that for all t ∈ [t0,t1 ],    t − t0 ||ˆz(t, ε )|| ≤ c1 exp −λmin (Q(y p )) . ε Tikhonov’s theorem [76] can now be used to summarize the previous result. Theorem 4.7 There exists a positive constant ε ∗ such that for all y p0 ∈ Dx , and 0 < ε < ε ∗ , the singular perturbation problem of (4.24)-(4.25) has a unique solution, x(t, ε ), z(t, ε ) on [t0,t1 ], and x(t, ε ) − x p(t) = O(ε ) z(t, ε ) − zˆ p(t/ε ) = O(ε )

(4.33) (4.34)

Chapter 4. Adaptive control

51

hold uniformly for t ∈ [t0,t1], where zˆ p (τ ) is the solution of the boundary layer model (4.29). Moreover, given any tb > t0 , there is ε ∗∗ ≤ ε ∗ such that z(t, ε ) = O(ε ) holds uniformly for t ∈ [tb,t1 ] whenever, ε < ε ∗∗ .

In order to extend this result to an infinite time interval, it is necessary that the boundary layer system is exponentially stable in a neighborhood of the tuned slow solution x p (t) for all t ≥ t0 . This may not be a simple proof, and it will be left for future investigation. Instead, the effectiveness of this approach is shown below using simulation. An intuitive yet not completely rigorous explanation for the resulting good behavior in the infinite time interval can be given with the help of Fig. 4.2. Notice that the Hurwitz nature of Jacobian (10.13) is only lost when y = x2 − x4 = 0. Since the fast motion, z, evolves with almost constant y (vertical lines in Fig. 4.2), y will not reach the value zero during this motion provided that y is initially far enough from zero. Once the slow manifold is reached, the slow variable will evolve in the domain z = 0. This domain corresponds to the case when the adaptation mechanism has reached its objective and parameter a is correctly estimated. In this domain y may reach the value zero but, intuitively, it is assumed that the system, once the adaptation law has reached the correct value, will present a behavior that is similar to the known load case, whose stability is proved in Chapter 3.

0.015

kzk

0.01

0.005

0 12 10 8 6

x2

4

−3

−2

0

−1

1

2

3

x1

Figure 4.2: Evolution of two trajectories in the state subspace (x1 , x2 , kzk). The last part of the trajectory is in the plane kzk = 0.

4.3. Conclusions

52

4.2.6 Simulations The inverter parameter values are the same as those for the known load case given in the Subsection 3.2.4, where the load resistance has been 50Ω and, therefore, a = 0.01. The desired output is 2 Vo = 220 √ sin(100π t). 2 At time t = 0s, the chosen value for the adaptation parameter is aˆ = 0.001 (which corresponds to R0 = 1000Ω, i.e. the relative error is 90%). Later, at time t = 0.5s, a load variation is produced from R = 100Ω to R0 = 1000Ω, such that, parameter a is again equal to 0.01 Once again, a commutation time of 50KHz is employed resulting a sample time of 0.1T s. Figures 4.3 and 4.4 show the evolution of the output voltage and the voltages of every boost DC-DC after t = 0s, respectively. Note that the circuit tends to the desired behavior. During this period the adaptation mechanism does not destabilize the system. Figures 4.5 and 4.6 show the evolution of the output voltage and the voltages of every boost DC-DC around t = 0.5s, respectively. Note that the circuit continues with the desired behavior. The time scale is the real time scale before the change of variable. Note that during this time, when the perturbation in the resistance and, thus, the corresponding adaptation mechanism is activated, the output signal does not undergo any significant variation. The adaptation of parameter a is represented in Figs. 4.7 and 4.8 where the load-change instants in the transition and steady-state are zoomed respectively. Note that the adaptation is very fast relative to the time scale of the system. In each of these graphs two evolutions are presented for two different values of ε . Note the smaller ε is, the faster the adaptation is.

4.3 Conclusions In this chapter an adaptation law is added to the previous control law in order to deal with a common and important problem in the field of electronics: slowly varying and/or unknown loads. For that, an adaptive control for unknown load is developed, which adapts the load very fast with respect to the time evolution of the system. The method is based on using a state observer on one-sided inverter and assuming that the state variables are measured. The stability of the complete system is proved by rewriting the system in the standard sin-

Chapter 4. Adaptive control

53

300 200

Vo (V )

100 0 −100 −200 −300 0

0.02

0.04

t(s)

0.06

0.08

0.1

Figure 4.3: Output voltage with the adaptation of an unknown load in t = 0s.

600

v1 , v2 (V )

550 500 450 400 350 300 0

0.02

0.04

t(s)

0.06

0.08

0.1

Figure 4.4: Output voltages of the first (solid) and second (dashed) boost DC-DC converters with the adaptation for an unknown load in t = 0s.

4.3. Conclusions

54

300 200

Vo (V )

100 0 −100 −200 −300 0.45

0.5

t(s)

0.55

Figure 4.5: Output voltage with the adaptation of an unknown load in t = 0.5s.

600

v1 , v2 (V )

550 500 450 400 350 300 0.45

0.5

t(s)

0.55

Figure 4.6: Output voltages of the first (solid) and second (dashed) boost DC-DC converters with the adaptation for an unknown load in t = 0.5s.

Chapter 4. Adaptive control

55

b) 0.01

0.005

0.005

0

0





a) 0.01

−0.005

−0.005

−0.01 0

2

−0.01 0

4

2

t(s) x 10−5

4

t(s) x 10−5

Figure 4.7: Time-evolution of the fast variable a˜ with ε = 0.01 in a) and ε = 0.001 in b). The perturbation is at t = 0s.

b) 0.01

0.005

0.005

0

0





a) 0.01

−0.005

−0.01

−0.005

0.5

0.5

t(s)

0.5

−0.01

0.5

0.5

0.5

t(s)

Figure 4.8: Time-evolution of the fast variable a˜ with ε = 0.01 in c) and ε = 0.001 in d). The perturbation is at t = 0.5s.

56

4.3. Conclusions

gular perturbation form; hence some relationship between the adaptation gain, γ , the observer matrix parameter, α , and the perturbed variable parameter, ε , are achieved. Another important relationship between the perturbed variable parameter, ε , and the system frequency, ω , is achieved in the analysis of the boundary layer fast subsystem. Finally, the stability is established by means of Tikhonov’s theorem [76]. The stability proof extension for an infinite time interval, will be an objective for further investigation. The assumption that voltage and current are measurable simplifies the observation problem. No persistant signals are required to prove the stability and the no noise is included in the measurable signals. As future work, an extension of this development could be done assuming that only the currents are measurable and including experimental results in order to validate all assumptions established in this work.

Chapter 5 Estimation of the attraction domain In the previous chapter, a controller has been developed for the boost inverter. Global stability of the closed-loop system has been guaranteed by a Lyapunov function. However, state and control-signal constraints including saturations have been disregarded. In this chapter, the objective is to estimate a satisfactory attraction domain for the boost inverter taking into account physical system constraints, containing control signal saturations. That is, to estimate a region composed of all initial conditions corresponding to trajectories that converge towards the right system behavior. The system will be driven into this attraction domain by a starting phase; this is common on the field of electronics [13,91,156]. This attraction region is estimated by using a novel method developed in this chapter. Estimating an attraction domain may be involved if there are physical system constraints. This problem may present a high degree of difficulty due to the system and control-law nonlinearities, including saturation-like constraints. The term saturation-like constraints is used for non-linear functions γ (u) that appear in the system model and they become the identity, γ (u) = u, in certain regions of the state space that include the desired behavior, (those regions are referred to regions in which such constraints are not active). Functions of this sort include typical control signal saturation as well as others, such as rate limiters, for example. Other constraints on the state variables can be considered as well. There exist many published methods to estimate the region of attraction (see, for example [54, 76] and the references therein). A kind of such methods is based on Lyapunov theory: closed Lyapunov-function level surfaces enclose (conservative) estimations for the region of attraction [76]. These methods often employ polynomial systems [85, 124, 146]. There exist powerful mathematical tools that can be used in the computation of the maximum acceptable level for polynomial systems [32, 68, 111, 151]. Some of these tools could be further developed for application to non-polynomial systems as well [31]. Application of

57

58

5.1. Problem statement

these methods would imply to seek for Lyapunov functions in order to be able to deal with the constraints. The search for a Lyapunov function by means of the numerical estimation method may be seen as an advantage as the user would not be required to propose a Lyapunov function. In cases like these, however, the computational method has to solve a much more difficult problem and it will be hard to tackle with systems of moderate complexity. Furthermore, saturation-like functions, which are one of the most common nonlinearities in practice, are usually out of the scope of these techniques. Usually saturation-like functions are only considered in the case of linear systems [5, 21, 110, 114]. A simple idea that solves the problem of estimation of the attraction domain for polynomial non-linear systems (and some others) with saturation-like constraints and state constraints is presented. The idea is to take advantage of the unconstrained global stability analysis and use this result to obtain a ‘conservative’ estimation of the region of attraction for the constrained case. Using the Lyapunov function of this analysis, there is no need to look for a Lyapunov function while estimating the domain of attraction and, thus, this problem becomes much simpler. On the other hand, the estimated attraction region is included in the domain where the saturation-like constraints are not active and, therefore, the method introduces a new source of conservatism. This obvious idea may be successful in hard problems when all other methods fail. In this chapter, this approach is employed to the boost inverter to estimate an attraction region. The application of the method produces good, albeit conservative, results. However, certain difficulties should be mentioned here: 1) the system equations are quite involved and can even present non-polynomial terms, namely rational functions; 2) the desired attractor is not an equilibrium point but a limit cycle. Both complications make the use of any other analysis method a formidable task. For sake of simplicity, in the application of this method to the boost inverter, the phase controller dynamic is not taken into account either the adaptation law.

5.1 Problem statement In Section 3.2 a control law (Eqs. (3.39)–(3.40)) for boost inverter not has been only designed, but has been also proved that for all initial conditions except the origin, the trajectories of the system tend to the curves Γi = 0 for i = 1, 2. Nevertheless, there are several constraints in the state variable that make this analysis useless from the practical point of view. These constraints are of several types: C1. Saturation-like constraint: sat10 (ui ) for i = 1, 2, makes control law (3.39)–(3.40) not to be feasible in the full state space. Therefore, this constraint is ‘soft’ in the sense that if

Chapter 5. Estimation of the attraction domain

59

the system arrives at a point where the constraints are violated, the analysis of Section 3.2 is no longer valid for the system with constraints. The point may nonetheless still lie inside the attraction domain of the desired limit cycle. Therefore, these constraints are saturation-like. C2. Capacitor voltages must be strictly positive in this circuit, which implies xi > 0; 2, 4. This is a ‘hard’ constraint since this situation should be avoided.

i=

C3. Finally, the control law is not feasible when any of the denominators in (3.39)–(3.40) are zero. This constraint is actually contained in C1, since denominators close to zero would imply large (positive or negative) values for ui ; i = 2, 4.

The objective of this work is to obtain a (possibly conservative) estimation for the region of attraction of the resultant system taking these physical constraints into account.

5.2 An approach of estimation of the attraction domain. Normally, in every control system, the control signal is subject to physical constraints such as saturations, rate-limiters, etc. As for control designs, such constraints are typically disregarded and the resulting control law is applied to the actuator. In this way, if the designed control law is ud = α (x), where x is the state variable, the actual control signal is u = γ (ud ) = γ (α (x)), where γ (·) is a saturation-like function. This approach is valid when the actual expression for u is used in the stability analysis of the resultant system. It is however quite common to neglect the actuator constraints in the stability analysis so as to simplify the analysis. In fact, the local stability property is not usually affected by these constraints, since in a neighborhood D of the desired attractor they are not active, that is, γ (α (x)) = α (x) ∀x ∈ D. However, the resulting attraction domain may be affected by constraints γ . This study deals with the estimation of this attraction domain based on a stability analysis that neglects the constraints. The analysis can also take into account state variable constraints in the following sense. Assume that there exists a ‘forbidden’ region in the state space. This means that the system state must remain within the boundaries of a pre-specified admissible (‘safe’) region. The estimation of the domain of attraction should take into account these constraints. Formally, the problem can be stated as follows: Actual system. Consider a control system type defined as such: x˙ = fa (x, u),

(5.1)

60

5.2. An approach of estimation of the attraction domain.

where x ∈ S ⊂ Rn , u ∈ Rnu . Function fa may include saturation-like functions. Furthermore, due to physical considerations, the state of the system must not go out of an admissible region T.  Unconstrained system Assume that an approximate model of the system is x˙ = f (x, u),

(5.2)

where function f : Rn → Rn and, besides, in D, f (x, u) ≡ fa (x, u). The reader can consider that this approximate model includes neither the saturation-like functions nor the state constraints.  Assume that a control law u = α (x) has been designed for the unconstrained model (5.2) for a given control objective. Remark 5.1 The control objective is not necessarily the stabilization of an equilibrium point, but perhaps the stabilization of limit cycles, for instance, as seen in the examples in Section 3.2, can be considered. Assumption 5.2 There is a widely known radially unbounded Lyapunov function V (x), in which a compact positively invariant set Ω, ∂∂Vx f (x, α (x)) ≤ 0. Let M be the largest invariant subset of the set for which V˙ = 0 in Ω. By the LaSalle invariance principle, assumption 5.2 guarantees that the trajectories of the unconstrained model tend to M . It is implicitly assumed that this is the desired behavior. Notice that if the original Lyapunov theorem is used to prove global stability, the previous assumption is also fulfilled. Assumption 5.2 also guarantees local stability for system (5.1). The problem lies in the estimation of the domain of attraction. The key is clearly to ensure that the system state remain within the boundaries of the region where saturations are not active, thus introducing new constraints. A conservative estimation of the region of attraction can then be easily obtained. The advantage of the relative ease with which it is obtained, however, is compromised by the fact that it may be far too conservative. Nevertheless, in many problems this simple idea may give satisfactory results. △

Assumption 5.3 Consider system (5.1) with control law u = α (x). Let be A = D ∩ T , that is, the intersection between the safe region and the region where the saturation-like functions are not active. It is assumed that this set can be estimated by a set of inequalities g(x) ≥ 0, where g : Rn → Rng .

Chapter 5. Estimation of the attraction domain

61

Now the problem can be transformed as follows:

Given a control system x˙ = f (x, u) with constraints in both the state variables and the control input g(x) ≥ 0, assume that a control law u = α (x) has been designed such that global stability is confirmed when no constraints are taken into account. The problem lies in estimating a region of attraction for the real system with constraints when this control law is applied.

A ‘conservative’ estimation for the attraction domain of the system with constraints is given by the following theorem:

Theorem 5.4 Under assumptions 5.2 and 5.3, assume that there exists a constant c > 0 such that in the set Ωc = {x : V (x) ≤ c}, the constraints g(x) ≥ 0 are satisfied. Then, all trajectories of the system with constraints starting at Ωc tend to M ∩ Ωc . Proof: Since in Ωc the saturation-like constraints are satisfied, the results for the unconstrained system are valid in Ωc . Therefore, V˙ ≤ 0 in Ωc and Ωc is positively invariant. Furthermore, since V (x) is radially unbounded Ωc is compact. The statement can be validated by applying LaSalle’s invariance principle.

Remark 5.5 Since M ∩ Ωc ⊂ M , the theorem guarantees the desired asymptotic behavior for the system with constraints.

Remark 5.6 As with other techniques for estimation of attraction domain, the present method is conservative. In this case the conservativeness is mainly due to two facts:

• The estimation of the region of attraction is restricted to domains in which V ≤ c. • The method considers the saturation-like functions as hard constraints. Nevertheless, there may be points where the saturations are active in the actual attraction domain.

Using Theorem 5.4, the problem is reduced to finding a value c > 0 such that g(x) ≥ 0 at the points where V (x) ≤ c. In order to use numerical tools for the determination of c, as will be seen in the next sections, the optimization problem can be stated as follows:

5.2. An approach of estimation of the attraction domain.

62

Problem 5.7 Maximize c subject to:

(V (x) − c) + pi (x)gi (x) − εi ≥ 0 i = 1, . . ., N,

(5.3)

where pi (x) are unknown semi-definite positive functions and εi > 0; i = 1, . . ., N. The purpose of constraint (5.3) is to validate the Theorem 1 hypothesis. To observe this, notice that at the boundary of the set Ωc , V (x) = c and, thus, the above constraints are reduced to pi (x)gi (x) ≥ εi > 0. As functions pi ≥ 0, the constraints gi (x) ≥ 0 are satisfied at the points on the boundary of Ωc . Furthermore, in the interior of this set, V (x) − c < 0 and, thus, these constraints are also satisfied. The pi functions lend even more degrees of freedom, thereby increasing problem feasibility. The small εi constants are pre-specified and are necessary in order to avoid problems at the points where pi (x) = 0. The introduction of εi parameters constitute a new source of conservatism. Remark 5.8 The region of attraction is estimated without the necessity of computing V˙ . In this work, sum-of-squares optimization is used in order to solve this problem. For this, a new assumption is needed. Assumption 5.9 Functions f (x, u) and g(x) are polynomial.

5.2.1 Sum of squares optimization Sum of squares optimization is an optimization technique based on the Sum Of Squares (SOS) decomposition for multivariate polynomials. A multivariate polynomial p(x) is said to be a SOS, if there exist polynomials f1 (x), ..., fm(x), such that: m

p(x) = ∑ fi2 (x) i=1

and therefore, p(x) ≥ 0 [117]. A SOS program has the following form [117]: Minimize the linear objective function rT c, where c is a vector formed from the (unknown) coefficients of:

Chapter 5. Estimation of the attraction domain

63

• polynomials pi (x), for i = 1, 2, ..., N1 • sum of squares pi (x), for i = N1 + 1, ..., N2 such that N

g0, j (x) + ∑ pi (x)gi, j (x) = 0

for j = 1, 2, , . . ., M1 ,

g0, j (x) + ∑ pi (x)gi, j (x) are SOS,

for j = M1 + 1, . . . , M2 ,

i=1 N

i=1

where w is the linear objective function weighting coefficients vector, and gi, j (x) represent certain scalar constant coefficient polynomials. Currently, SOS programs are solved by reformulating them as semi-definite programs (SDPs), which in turn are solved efficiently, e.g., using interior point methods. Several commercial as well as non-commercial software packages are available to solve SDPs. SOSTOOLS [116] is a Matlab toolbox that performs this conversion automatically, calls the SDP solver, and converts the SDP solution back to that of the solution of the original problem. The problem stated in the previous section can be addressed solving the following SOS problem: Problem 5.10 Maximize c subject to: (V (x) − c) + pi (x)gi (x) − εi are SOS;

i = 1, . . . , N,

(5.4)

where pi are unknown SOS polynomials. This problem is more restricted than that presented in the previous section. Nevertheless, any solution to SOS problem 5.10 is a solution to problem 5.7. Remark 5.11 Assumption 5.9 can be relaxed since other types of functions, such as trigonometric functions [109] or rational functions (e.g., the application examples in the next section) can be considered.

5.3 Application to the boost inverter The method developed before is used in order to obtain an estimation of the attraction domain for the boost inverter when system physical constraints are taken into account. This method

5.3. Application to the boost inverter

64

is useful because the model and control law have a relative high degree and complexity. Furthermore, the global stability proof for the unconstrained problem is available by means a Lyapunov approach. In Chapter 5.1 it has been proved that, under no constraints, all trajectories (except the one starting at the origin) of system (2.5)–(2.8) with control laws (3.39)–(3.40) tend to the desired limit cycle. Remind that the Lyapunov function used is: V=

Γ21 Γ22 + . 2 2

(5.5)

The constraints are (only constraints C1 and C2 are presented here; constraint C3 will be discussed later): • ui (x) ≤ 1

i = 1, 2

• ui (x) ≥ 0

i = 1, 2

• x2 > 0 • x4 > 0. The expressions for u1 and u2 , which are given by (3.39) and (3.40) are not polynomial but rational functions. Nevertheless, writing them as quotient of polynomials ui (x) = ni (x)/di (x) all the constraints can be formulated in standard form. For this, it can be assumed that polynomial d(x) does not vanish at any point of the objective curve Γi (x) = 0 i = 1, 2. Otherwise, control laws (3.39)–(3.40) are not valid for this problem. Therefore, the sign of d(x) is constant along Γi (x) = 0 and, by continuity, in a neighborhood of this curve. By numerical inspection, it can be checked that, for the circuit parameters given below, d(x) > 0 on Γi (x) = 0. With this consideration in mind, constraints (1)–(3) can be written as polynomial constraints: • di (x) − ni (x) ≥ 0 i = 1, 2 • ni ≥ 0

i = 1, 2

• x2 > 0 • x4 > 0 Thus, the problem to solve is minimize (−c)

(5.6)

Chapter 5. Estimation of the attraction domain

65

subject to: (V (x) − c) + p1 (x)(d1 (x) − n1 (x)) − ε1 ≥ 0

(5.7)

(V (x) − c) + p3 (x)n1 (x) − ε3 ≥ 0

(5.9)

(V (x) − c) + p2 (x)(d2 (x) − n2 (x)) − ε2 ≥ 0

(5.8)

(V (x) − c) + p4 (x)n2 (x) − ε4 ≥ 0

(5.10)

(V (x) − c) + p6 (x)x4 − ε6 ≥ 0

(5.12)

(V (x) − c) + p5 (x)x2 − ε5 ≥ 0

(5.11)

Notice that constraints C3 are considered in the previous set of constraints. Indeed, constraints (5.9)–(5.10) implies n1 (x) ≥ ε3 > 0 and n2 (x) ≥ ε4 > 0, respectively, for V (x) ≤ c, while constraints (10.15)–(10.16) implies d(x) ≥ n1 (x) + ε1 and d(x) ≥ n2 (x) + ε2 , respectively, for V (x) ≤ c. This implies that d(x) > 0 in V (x) ≤ c. The following analysis can be directly modified for the case when d(x) < 0. Results The values of the circuit parameters are taken from Subsection 3.2.4. Software SeDuMi [144] was used as the SDP solver under SOSTOOLS. The values for parameters εi are chosen equal to 10−6 while the chosen order for the unknown polynomials pi is 3. The solution is obtained in approximately ten minutes on a PC (1.66 GHz Intel Core2): 23.26. This result is probably conservative as has been pointed out in Remark 5.6. As a way to corroborate this result, by numerical inspection, it has been found that for x(1) = (0 − 0.1 0.2 5.8)⊤, which corresponds to V (x(1) ) = 33.02, the constraint x2 ≥ 0 is violated.

5.4 Conclusions The problem considered in this chapter is the estimation of the attraction domain for the boost inverter with the control law proposed in Chapter 3, which does not present a global stability due to certain physical constraints. For this a method for the estimation of the attraction region considering general physical constraints is presented. This approach can be applied to systems with a global Lyapunov stability achieved without considering saturations and other kind of limitations. This is a common situation since, saturations are neglected in many stability analysis. The idea is to take advantage of the Lyapunov level sets, finding the maximum Lyapunov level in such a way that constraints are fulfilled inside it. This makes that the computed attraction domain is a ‘conservative’ estimation. This method is useful even when the degree and complexity of the equations is high.

66

5.4. Conclusions

For application of the method, powerful computational tools exist when the system is polynomial, such as SOS. Therefore, the closed-loop system needs to be polynomial or rational (however, there exist cases where SOS programming have been applied to trigonometrical and other terms [109]). Consequently, the problem is transformed in a sum of squares optimization problem. Conservativeness of the method has also been discussed. In the application to the boost inverter the system as the constraints are rewritten in a polynomial form, as the method requires. For simplicity reasons neither phase controller nor adaptive control has been taken into account.

Part II Controlling a DC-DC Vdd-Hopping converter

67

Chapter 6 Introduction The development of low-power electronic devices has raised up in recent years. Very-LargeScale Integration (VLSI) is mostly used in products related to information technology, such as PCs, mobile devices and digital consumer equipments. Motivated by the Moore’s law and market evolutions [132], ARAVIS project (Architecture avanc´ee Re-configurable and Asynchrone pour Video et radio logicielle Int´egr´ee Sur puce) sponsored by Minalogic1 , looks for architecture and design solutions that allow the production of embedded computational platforms in its scalability limit. It proposes a generalization of certain techniques in order to obtain a solution to the technology variability problem in 32nm, what will represent an input toward the development of a new paradigm. This part of the thesis is included in the ARAVIS project context. Currently System on Chips (SoCs) technology: 90nm, 65nm and, even, 45nm can not be applied any more to the technology of 32nm due to the semiconductor material behavior in small scale. The main problem in this kind of scale is the occurrence of technology variability phenomenon [154], that generates quite disparate performances in a same chip. Consequently, a new architecture must be developed in order to answer to this issue. In Fig. 6.1, an example of this problem is shown. It presents a fault or low performance of a computational node2 in high speed. The ARAVIS project is focused on three technology keys:

• re-configurable structure with respect to applicability requirements. It can be accomplished by programming the flexible interconnections between the clustered nodes of the SoC computational unit [113]. 1

http://www.minalogic.com/ common computational unit of a SoC is composed of clustered nodes [74]

2 The

69

70

Figure 6.1: Technology variability problem in a chip in a computational node working in high speed.

• Globally asynchronous locally synchronous method [92], in order to release the communication constraints between remote points, and • dynamic management of the power consumption and activity with respect to constraints are achieved by control theory application [30, 65, 122]. The last key looks for achieving a good trade-off between activity, power consumption and Quality of Service (QoS), what is one of the challenges in future embedded architecture. This dynamic management is especially difficult for 45nm and 32nm, which are at the limit of the scalability. Figure 6.2 shows chips integration in 32nm scale.

Figure 6.2: Integration of future 32 nm chips.

Advanced control strategies have, therefore, to be designed in order to make the performance fits the requirement, minimizing the energy losses on a chip. This second part of the

Chapter 6. Introduction

71

thesis focuses on designing control strategies for certain loop of the chip that optimizes the energy consumption.

6.1 Optimization of the energy consumption in SoCs As has been mentioned before, one of the challenges in the ARAVIS project is the energyconsumption reduction in SoCs, which can be got by means of automatic control methods. Control loops can be applied in different architecture levels: cluster frequency and voltage supply, cluster computational power and management of the quality of service provided by the application [30]. Figure 6.3 shows these different control loops.

Cluster

Cluster

f, Vdd

f, Vdd Processing Nodes

Processing Nodes

Voltage Controller

Voltage Controller

Vdd - Hopping

Vdd - Hopping

Programmable Ring

Programmable Ring

Energy Controller

Energy Controller

Speed1, No. of Instructions1, Deadlines1

Speed2, No. of Instructions2, Deadlines2

QoS Controller

Figure 6.3: Sketch of the three control loops.

Generally, the power consumption can be reduced if the local core voltage or/and the clock frequency are decreased. For this dynamic control loop in frequency and voltage is more and more important in the embedded systems [122]. Thus, these control loops allow the adaptation of computational power in the cluster level, and hence, the power-saving.

6.2. Vdd-Hopping DC-DC converter

72

Another control loop is used in order to achieve an energy-performance trade-off. The system must not always work in the maximum power level if it can work in other lower power levels. This is possible if each task is performed before to certain deadline [45]. In [39], a control solution has been presented, applying ‘robust control’. This solution reaches a control that can reduce or, even, reject the influence of the variability problem. The last control loop employed in this kind of architecture looks for a trade-off between QoS, computational limitations and global energy consumption.

6.2 Vdd-Hopping DC-DC converter Microprocessors in SoC have a computational load that requires a time-varying performance. Consequently, the SoCs can achieve a substantial energy efficiency, if they reach to operate at the minimum performance level required by the active software processes. Dynamic Voltage Scaling (DVS) is a known technique that manages the system power consumption [27, 89, 136]. The operation principle is to adjust the processor supply voltage to the minimum level of performance required by the system application. Therefore, DC-DC converters are a key element in a DVS mechanism, since they can adapt this supply voltage. However, these converters have a different structure than the standards ones, because they must change the operating voltage in a dynamic way [23]. A dynamic continuous buck converter for DVS systems was presented in [88], which provides good performance. It, however, limits SoC scaling properties due to the size of inductive component. In the framework of SoC miniaturization, a Vdd-Hopping DC-DC converter was developed, whose basic structure is showed in Section 6.4. Note that this converter is composed by two supply sources and a Power Supply Selector (PSS) [99]. In this structure, the inductive element is replaced by a set of PMOS transistors3, reducing the converter required size. Vdd-Hopping converter is specially interesting, because it may deal with Local Dynamic Voltage Scaling (LDVS) [23, 157] adapted to Globally Asynchronous and Locally Synchronous Systems (GALS) systems [92]. The main idea for GALS systems is to replace the global clock by several independent synchronous blocks. Every synchronous block operates with an own internal clock and they are communicated asynchronously by each other. This mode of operation provides additional flexibility. This fact allows to use energy-aware converter structures such as DVS architectures. DVS is applied in every synchronous block, that is why, it is called LDVS. 3 P-channel

Metal-Oxide-Semiconductor field-effect transistor.

Chapter 6. Introduction

73

The Vdd-Hopping converter, as has been mentioned before, is basically made up of a PSS and two external supply voltages. A high voltage supply, Vh , for a unit running at nominal speed and a low voltage supply, Vl , for a unit running at reduced speed. vc is the output voltage of the system.

Vh

Vl

Set of PMOS

uk e CLK CONTROL

vc − +

vr

LPM

Figure 6.4: DC-DC Vdd-Hopping converter structure.

A PSS is constituted by a group of PMOS transistors connected in parallel with common drain, source and bulk, but separated gates in order to scale the output voltage from a low voltage level to a high voltage level (rising transient-period) and from a high voltage level to a low voltage level (falling transient-period). Furthermore, the PMOS transistor that connects the Vl to the voltage output vc is switched on when Vl is the selected power supply. This reduces the dissipated energy when the unit running is at low speed. Other component in the PSS is a control block that provides a control signal uk for the PMOS transistors. Besides it generates a reference signal vr . Likewise, this control block has as inputs: a clock signal (CLK), a local power manager signal (LPM), that orders to the PSS to start the hopping sequence and the error voltage signal from a comparator e. Load model The load model taken in this work is an impedance which depends on the core voltage, vc , and sometimes, also on the clock frequency, ωn [98]. In this thesis, the load model provided in [98] is employed. This model approaches a low frequency in function of the core voltage, thus the load only depends on the voltage vc . It is composed of a current supply, Ileak , a capacitance, C, and a dynamic resistance, rL , representing the dynamic and short-circuit consumption. It is shown in Fig. 6.5.

6.3. Non-linear control application to Vdd-Hopping DC-DC converter

74

C

rL

Ileak

Figure 6.5: Load model.

6.3 Non-linear control application to Vdd-Hopping DC-DC converter Control objective for this kind of discrete DC-DC converter is to achieve the target voltage providing a correct and reliable operation during the switching transitions, allowing to accomplish the main ARAVIS projet objective: SoCs miniaturization. Therefore, the control must achieve: • high energy efficiency, • system stability, • small current peaks, • fast transient periods, • robustness with respect to parameter uncertainty, • robustness with respect to delays and • easy implementation. Figure 6.6 shows the control problem for the DC-DC Vdd-Hopping converter. A simple discrete controller was proposed in [99] to handle the two voltage levels of the Vdd-Hopping technique. In this control structure only one transistor can be switched at each sampling time. This limits the ability of the converter to make fast transient-periods, and hence the possibility to optimize the energy consumption. In addition, the employed voltage reference was a ramp with a computed slope to obtain the smallest possible current peaks.

Chapter 6. Introduction

75

LDVS

GALS

Local Dynamic Voltage Scaling Architecture

Globally Asynchronous and Locally Synchronous Systems

Vdd-Hopping Converter PSS

+ Two voltages sources

CONTROL!!!!!

Stability of the closed-loop system

high energy efficiency, small current peaks, fast transient periods, robustness with respec to parameter uncertainty, robustness with respec to delays and easy implementation.

Figure 6.6: DC-DC Vdd-Hopping converter control problem.

6.3.1 High-performance controllers In Chapter 7, a set of high-performance control laws for the DC-DC Vdd-Hopping converter is developed, without the constraint that only one transistor can be switched at each sampling time. [8]. This allows to obtain a richer set of control sequences and, thus, a better expected performance with respect to the issues previously mentioned. The controllers are developed in order to improve the tracking capability and its regulation characteristic. As a side effect, it is also observed that the transient current peaks are reduced. However, the computational cost is increased. The different controllers are compared in terms of: transient response, quality of the induced load current and power consumption. It will be seen that the most suitable controller is the one based on Lyapunov theory. However, this controller is enhanced in order to achieve a high energy saving, minimum current peaks and a suitable performance with unknown load resistive parameter, which are very common properties in the field of microelectronics [70, 103]. Firstly, an optimal evolution of the voltage reference will be computed applying optimal control theory [80, 86, 155]. This optimal voltage reference is computed looking for minimizing the current peaks as well as energy losses. Secondly, an adaptive controller is proposed for an unknown load resistive component, offering a suitable system performance by simulation even if it is time-varying. Figure 6.7 shows a sketch of this high performance controller.

6.3. Non-linear control application to Vdd-Hopping DC-DC converter

76

In spite of the good advantages that this enhanced controller offers, it presents much more computational blocks (more complex implementation) than the controller proposed in [99] and thus, although, it reduces notably energy consumption, it does not corresponds to a feasible solution for the ARAVIS project context. For this, another controller is developed in Chapter 8 based on the one that presents the easiest implementation from the set of controllers previously presented. This new developed controller must be easy to implement and maintain all the good characteristics obtained by the controller mentioned previously in Chapter 7. ADAPTIVE CONTROL

OPTIMAL REFERENCE

v∗

+ CONTROL

-

Vdd-HOPPING CONVERTER

vc

Figure 6.7: Vdd-Hopping converter closed-loop with optimal reference and adaptation parameters.

6.3.2 Energy-aware controller The previous proposed controller provides very suited properties. The only drawback is that it has a complex implementation, what can exclude it from the ARAVIS project scope. Among the set of controllers presented before, there is one with a simple implementation. In Chapter 8, this controller is selected and developed taking into account the objectives given before to obtain a suitable controller for the ARAVIS project. The controller with a feasible structure is a ‘linear controller, what makes it to present a simple implementation. Hence, it takes a relevant interest in the industrial applications [17, 84, 87, 120]. Nevertheless, this controller does not provide the best system performance. Now, from this structure, an ‘advanced linear controller’ is developed in Chapter 8. An innovative approach based on saturations with time-varying limits that manages the current peaks during the transient periods is proposed for this controller. Furthermore, energy-efficiency is improved when a step voltage reference is used instead of a ramp voltage reference. Consequently, it does not only reduces the current peaks but achieves also a fast transient response and reduces energy dissipation. In summary, this proposed ‘advanced linear controller’ is focused on limiting the current peaks while transient periods are reduced, being energy-aware.

Chapter 6. Introduction

77

Generally, controllers applied in the industry are implemented in discrete-time. Hence that, the previous advanced controller is discretized. This control law is patent pending under the name of ENergy-AwaRe Control (ENARC) [6]. In order to show the consumption energy saving that this controller may achieve, a comparison with a previous controller published in [99] is performed. It is shown that this controller can diminish the energy consumption 96% with respect to the previous published controller. In addition, other control objectives are achieved: faster transients, small current peaks, and easy implementation. For the last characteristic, it is remarked that this controller needs less computational blocks than the other controllers presented above. The stability property is analyzed in Chapter 9. Thanks mainly to its innovative current-peak aware the ENARC control becomes attractive for industrial applications. These results will be validated in the ARAVIS project by using VHDL-AMS4 simulator, since the Vdd-Hopping system with the load is an hybrid system between analog and digital elements (the controller will be implemented digitally). It is expected that the ENARC controller is physically implemented in the new generation of SoCs of 45nm and/or 36nm developed in the project context.

6.3.3 Approximate stability analysis for the energy-aware controller In Chapter 9, an approximate global stability analysis of the equilibrium is developed for the Vdd-Hopping system with the continuous-time version of the ENARC controller [9]. This controller coincides with the ‘advanced linear controller’ mentioned before. For sake of simplicity, this analysis is performed in continuous-time and it is not rigorous. This analysis allows to have an intuition of the closed-loop system behaviour with the discrete ENARC controller. This assumption is very common [71, 152]. As the Vdd-Hopping model as the continuous-time version of the ENARC controller are nonlinear. Among their nonlinearities, there is a saturation with limits depending on system state. This makes the analysis involve. For simplicity, a preliminary stability analysis is performed for the Vdd-Hopping system with the ‘linear controller’ proposed in the set of high-performance controllers developed in Chapter 7. This analysis is based on LaSalle’s invariance principle [76]. Then, it is extended to the continuous-time version of the ENARC. This continuous-time controller has a saturation mechanism, the system can work in three 4 It

is a hardware description language used in electronic design automation. It is capable to define the behavior of mixed-signal systems, because it describes digital systems as well as analog systems. This is a very used powerful software.

6.3. Non-linear control application to Vdd-Hopping DC-DC converter

78

operating-modes: non-saturated system (Region I), saturated system in the upper limit (Region II) and saturated system in the lower limit (Region III). It is proved that the equilibrium is placed in Region I. The stability proof is performed in two parts. Firstly, it is proved that the system in saturated mode converges to the non-saturated mode in finite time, and that, it crosses the saturation limit to the Region I without being able to return towards the saturated mode. Secondly, it is also proved that the system converges to the desired point in Region I. This analysis applies LaSalle’s invariance principle for a bounded domain limited by a Lyapunov level curve and the saturation limit lines. Figure 6.8 represents this idea. Once the global stability analysis is ensured for the ENARC controller, the control objective that must be dealt with are: the robustness under delay presence and parameter uncertainties.

e Ω2

Region II

Region I

P2

Region III

σ Figure 6.8: Representation of the system operating regions

6.3.4 Advanced energy-aware controller Previously, an energy-aware controller has been presented, which not only satisfies certain requirements for SoCs, but their stability properties have been also proved. Now, in Chapter 10, this controller is improved in order to satisfy the control objective of achieving robustness with respect to uncertain parameters and delays. The delays are presents in the input and output of the control block for computational and synchronization issues, respectively. So far, delays and parameter uncertainties have not been considered in the Vdd-Hopping

Chapter 6. Introduction

79

system. In SoC, however, these issues are very usual [112, 140]. The controller presents delays due mainly to two reasons: synchronization issues and power-performance tradeoff [45]. In addition, another point to take into account is the parameter uncertainties, which can cause an unpredictable impact on the power, performance and reliability of the system [34, 93]. In Chapter 10, the ENARC controller is improved dealing with delays and parameter uncertainties. A sub-optimal ‘conservative’ control tuning approach for the control gains is presented based on linear control theory. This method is sub-optimal because, for simplicity, it is developed for an approximate ENARC controller version that does not consider the current-peaks management. The sub-optimal control gains are obtained solving a H∞ control problem [26, 53, 101]. This problem is dealt with Lyapunov Krasovskii theory [53, 138], which provides some stability conditions through Linear Matrix Inequalities (LMIs). Consequently, a robust equilibrium stability as well as a robust disturbance rejection under parameter uncertainties are ensured for a delay-time system. This kind of approach are used in industrial applications [14, 129]. In this chapter, the robustness properties achieved by the closed-loop system with the sub-optimal constants applied to the ENARC controller are shown by simulations. In summary, the ENARC controller with the sub-optimal constants applied to the VddHopping DC-DC converter can satisfy all the control objectives mentioned before for lowpower technology. Therefore, it can have an important impact in the recent trends in the miniaturization of electronic devices.

80

6.3. Non-linear control application to Vdd-Hopping DC-DC converter

Chapter 7 High-performance control for the DC-DC Vdd-Hopping converter The development of low-power electronic devices has raised up in recent years. Very-LargeScale Integration (VLSI) is mostly used in information technology related products, such as PCs, mobile devices and digital consumer equipments. In a System on Chip (SoC), several levels of supply voltages are required to reduce power consumption. It can be reached applying the Dynamic Voltage Scaling (DVS) concept, which is an interesting method that manages dynamically the microprocessor supply voltage Vdd according to various loading conditions. DC-DC converters may be used in order to apply the DVS concept. Among this kind of converters, a high efficiency discrete converter is found: the DC-DC Vdd-Hopping converter [99], which is implemented applying, as its name references, the Vdd-Hopping approach [75, 106, 128]. This technic delivers two distinct small voltage levels with a very small current, according to the required performance level. Consequently, it achieves a high energy-efficiency. Therefore, its operation principle is to vary the voltage from a low voltage level to a high voltage level, and reciprocally. A controller for this DC-DC converter must be developed taking into account the context where it will be implemented. One of the main control problem in low-power DC-DC converters is to achieve a high energy-efficiency. Furthermore, DC-DC Vdd-Hopping converters must be able to adapt to various loading conditions and achieve high efficiency over a wide load-current range, which is critical for extended battery life. Moreover, to keep the rate of change of the device voltage providing a correct and reliable operation during the switch transition is also important. In this chapter, a set of controllers for the Vdd-Hopping converter is developed, obtaining 81

82

7.1. Mathematical model of Vdd-Hopping mechanism

a high-performance. Among these proposed controllers, the one that offers a best transient performance is selected and enhanced with the aim to deal with the unknown resistive component of the load as well as to minimize the dissipated energy and current peaks, what is very important in the field of microelectronics. Current peaks and power consumption are minimized by computing an optimal evolution for the voltage reference. Likewise, an adaptive controller is proposed to deal with the unknown load resistive parameter. Consequently, the obtained high performance controller can acquires a high consideration on electronic devices. Generally, the power consumption in a SoC can be reduced if the local core voltage or/and the clock frequency are decreased, that is why a GALS systems is developed in the ARAVIS project where the Vdd-Hopping converter is embedded. Likewise, the clock frequency has to satisfy that the task (for instance, the execution of the control laws designed here) is performed before a deadline and that the minimum required local clock frequency that guaranties the critical path (longest path delay) of the whole chip is fulfilled [45]. Normally, in order to take into account all of these issues in a GALS system when the frequency and voltage have to rise, firstly the voltage is rising and later the frequency is rising. On the contrary, when the frequency and voltage have to fall, firstly the frequency is fallen, and later the voltage is fallen.

7.1 Mathematical model of Vdd-Hopping mechanism The discrete DC-DC Vdd-Hopping converter presented in [99] for SoCs shows several advantages: high efficiency and reduced size, since it does not need passive components. This converter uses the Vdd-Hopping technique [75, 106, 128] in order to obtain a LDVS architecture for a GALS system. Figure. 6.4 shows this connected Vdd-Hopping structure.

7.1.1 Mathematical model for control design For simplicity the low voltage supply, Vl , is disregarded for control design purposes (see Fig. 6.4). The main objective is that the core voltage vc achieves the high and low voltage levels by switching the PMOS transistors. In this configuration,, at least, one transistor must always be switched on. Figure 7.1 shows an electrical representation of the Vdd-Hopping converter without the low voltage supply, Vl , connected to the load that has been described in Section 6.2. Assumption 7.1 PMOS transistors are modeled as ideal resistors when they are switched

Chapter 7. High-performance control for the DC-DC Vdd-Hopping converter

83

CONTROL

uk M1

R1

M2

R2

Mn−1

Rn−1

Mn

RN

+

Il (vc)

IMPEDANCE

Vh

vc

Z(vc)



Figure 7.1: Vdd-Hopping, voltage supply and load.

on and as resistors with infinite resistance when they are switched off. They are considered to have the same electrical characteristic.

The voltage loop equation yields the relationship Il (vc ) =

Vh − vc , Ruk

where

Ruk ,

R0 . uk

(7.1)

uk is the number of transistors switched on, thus, uk ∈ U = {1, 2, ..N} and it is the control variable. Likewise, R0 is the PMOS transistor resistance. In this kind of system, it can be assumed that all transistors have the same transistor resistance R0 = R1 = R2 = . . . = RN . The current through the set of PMOS transistors, Il , depends on voltage vc and control signal uk . Thus, Il varies during the hopping transients. In this work, as mentioned in Section 6.2, the load model presented in [98] is employed: Il = f (vc , fclk ) = Idyn + Ishort + Istat + Icap Idyn = Kdyn fclk vc Ishort = Kshort fclk (vc − 2Vth )3 Ileak = Kleak dvc Icap = C , dt

(7.2) (7.3) (7.4) (7.5) (7.6)

where C, Kdyn , Kstat and Kleak depending on the real consumption estimation. Vth is the

7.1. Mathematical model of Vdd-Hopping mechanism

84

threshold voltage, which must fulfill vc > 2Vth . And fclk , which represents the clock frequency, may be approximated as a function of the core voltage, vc . Figure 7.2 shows the representation of the load model used in this research. rL represents the dynamic resistance. For simplicity reasons, firstly, an constant average value of rL is taken in order to design controllers. Later, the real time-varying parameter, rL , will be taken into account. The averaged load resistance RL is given by 1 RL , t f − t0

Z tf t0

rL dt.

where t0 and t f are the initial and final time, respectively, in the rising transient period. Assume that RL has the same value in the falling transient period. Il

Icap

Idyn + Ishort

rL

C

IMPEDANCE

Ileak

Figure 7.2: Load model.

Let us combine the specific form of the load, Eqs. (7.2)–(7.6), with system (7.1). The voltage equation can be expressed as v˙c = −β vc + b(Vh − vc )uk − δ ,

(7.7)

where • β , CR1 L > 0 and δ ,

Ileak C

> 0 depend on the load.

• b , CR1 0 > 0 depends on PMOS resistance, R0 , and on load parameter C. Define the voltage error as: e , vr − vc , where vr is a voltage reference. Thus, the associated error voltage equation is e˙ = −β e + b(vr −Vh )uk − buk e + β vr + δ + v˙r .

(7.8)

Chapter 7. High-performance control for the DC-DC Vdd-Hopping converter

85

7.2 Control laws The objective of this section is to present some high-performance control laws for the VddHopping converter, fulfilling the requirements mentioned before. All these controllers are designed to provide stable behaviors, and different control methodologies are used in each case. Some simulations are performed, such that the behavior of the closed-loop system with the different controllers are shown. In these simulations, N = 24 is taken as the total number of PMOS transistors in the model shown in Fig. 7.1. Note that, at least, one active transistor must be always switched on. The voltage supply is Vh = 1.2V . The reference signal, vr , follows a linear time evolution between the low voltage level Vl = 0.8V and the high voltage level Vch = 1.2 − ∆h . This signal has a slope specified by the designer, which is inspired by [99]. Remark 7.2 For physical reasons, the maximum voltage achieved, vc , must be Vch = Vh − ∆h where ∆h ∈ R and is small. ∆h depends on the voltage supply, PMOS resistance and load parameter. The system resistances are RL = 27.7Ω and R0 = 31.41Ω, the capacitance is C = 9nF while Kleak = 1.67 · 10−3, the threshold voltage is Vth = 0.4V , and system clock frequency is ωn = 500MHz. The sampling frequency has the same value that the clock frequency. The difference between the high voltage supply, Vh , and the high core voltage, Vc h , is ∆h = 0.08V and the slope of the reference signal, vr , is 1.067 · 106V /s.

7.2.1 Control proposed in [99] The development of the set of the high performance controllers for the Vdd-Hopping converter is inspired by the ‘intuitive control’ used in [99], under the form: uk = satN 1 {uk−1 + sign(e)}

(7.9)

In this law, no more than one transistor switches at each sampling time according to the sign of the voltage error. Previously to the works developed in this thesis [8, 9], control (7.9) was the only published controller. Therefore, this controller has the limitation that one only transistor can be switched on or off at every sampling time. Figure 7.3 shows the implementation of this controller. Likewise, Fig. 7.4 shows a simulation for this controller by using Matlab. Note that the performance presents an oscillatory behavior, with important current peaks.

7.2. Control laws

86

CONTROL

vr +

1

e

vc −

−1

+

uk−1

satN 1

uk

+

z −1

Figure 7.3: Intuitive control from [99].

In what follows, other control alternatives are proposed without the constraint that only one transistor can be switched on or off, as long as the number of transistor is limited by 1 and N. This introduces a saturation in every control law for the Vdd-Hopping converter. Remark 7.3 Every control law is designed in such a way that the desired output voltage corresponds to one of the saturation bounds 1 or N, since they corresponds to the lower or higher voltage level, respectively. Assumption 7.4 The sampling time is chosen in such a way that the controllability and observability properties are preserved. Control laws will be designed using directly the nonlinear continuous-time equation (7.8). This will lead to a continuous-time controller expression that will be approximately discretized. This approach is very common in the field of automatic control [71, 81]. The implementation of these discrete-time controllers are shown by block diagrams. The time evolution for the reference signal employed in [99] is maintained in the simulations of each developed controller in this section. However, later, it will be seen that, by means of choosing an suitable reference, the closed-loop system performance can be enhanced.

7.2.2 Controller No. 1: linear controller The first proposed controller is based on a linear structure, namely a PI (ProportionalIntegral) controller. This controller is the most common industrial control solution [17, 84, 120]. It has to be designed to cope with possible steady-state errors.

Chapter 7. High-performance control for the DC-DC Vdd-Hopping converter

87

a) NTrans

Nmax=24 20 10

Nmin=1 0 0

0.1

0.2

0.3

0.4

0.5 t(s) b)

0.6

0.5 t(s) c)

0.6

0.5 t(s)

0.6

0.7

0.8

0.9

1 −6

x 10

V(V)

Vh= 1.2 Vch=1.12 1

Vl= 0.8

I(A)

0

0.08 0.06 0.04 0.02 0 0

0.1

0.1

0.2

0.3

0.2

0.3

0.4

0.4

0.7

0.8

0.9

1 −6

x 10

0.7

0.8

0.9

1 −6

x 10

Figure 7.4: Intuitive control. Evolution of: a) number of PMOS transistors switched on, b) vr (dashed) and vc (solid), c) current. Il .

The proposed control law is: uk = satN 1 round {K1 e + K2 σ } , where σ corresponds to

R t+Ts t

edt and being Ts =

1 ωn

(7.10)

the sampling time.

The constants K1 and K2 are tuned off line. The tuning process must take into account the sampling frequency, low and high voltage level, number of transistors, load parameters and PMOS resistance. A tuning method is proposed to ensure the right system performance. For this, the closedloop system is linearized around a set point, such that, K = [K1 , K2 ] are defined by ensuring that A + BK is Hurwitz. This tuning mechanism is 2ξ ωn − (ukl b + β ) b(Vh −Vl ) ωn2 K2 = , b(Vh −Vl ) K1 =

(7.11) (7.12)

being ukl the minimum value of transistors that have to remain switched on, i.e., the minimum value of uk , uk l = 1, and ξ is a design parameter. In Chapter 9, constraints on ξ are obtained, such that, the closed-loop system is asymptotically stable to the equilibrium point [e, σ ] = [0, σ¯ ]. The closed-loop system with control (7.10) can suffer wind-up. This phenomenon has

7.2. Control laws

88

not been studied and it is a future work. It is highlighted that in the performed simulations this effect has not displayed. The discretization of this controller is: uk = sat1N {uk−1 + round(K¯ 1 ∆ek + K¯ 2 ek )} ,

(7.13)

where K2 K¯ 1 , K1 − 2 ¯ K2 , K2 Ts .

(7.14) (7.15)

Equations (7.14)–(7.15) are common relationships between continuous- and discretetime system [104]. Next up, Fig. 7.5 shows the approximate discretization of this controller, for digital implementation. CONTROL

¯2 K

vr +

ek

+

vc −

z

¯1 + K

+

− −1ek−1

+

uk−1

round +

satN 1

uk

Figure 7.5: Digital PI controller.

Figure 7.6 shows a simulation of the system (7.8) with control (7.13). In this simulation is chosen as design parameter value ξ = 0.05. The interval I defined in Lemma 9.1 (in Section 9.1) results [0.01, 0.08]. Therefore, the equilibrium stability will be able to be guaranteed according to Chapter 9. This controller could be the most suitable for physical implementation not only because it requires a reduced number of computational blocks, but also because it does not require model information. Moreover, as it can be seen in Fig. 7.6, it provides better performance in both voltage and current variables with respect to the intuitive controller. Nevertheless, it still presents some peaks in the current signal.

Chapter 7. High-performance control for the DC-DC Vdd-Hopping converter

89

In summary, this control law can make that the system achieve the suited stability properties, as will be seen in Chapter 9, with the characteristic to requires a reduced number of computational blocks.

NTrans

a) 20 10 0 0

0.2

0.4

0.6

0.8

V(V)

t(s) b) 1.2 1 0.8 0

0.2

0.4

x 10

0.6

0.8

I(A)

t(s) c) 0.05 0 0

0.2

0.4

1 −6

1 −6

x 10

0.6

0.8

t(s)

1 −6

x 10

Figure 7.6: Control No.1. Evolution of: a) number of PMOS transistors switched on, b) the vr (dashed) and vc (solid). c) the current Il .

7.2.3 Controller No.2: feedback linearization The second proposed controller is designed by using feedback linearization technique. This leads to a continuous-time linear system in closed-loop. The aim of this method is that the closed-loop system becomes e˙ = −K3 e − K4 σ ,

which achieves the suited properties of a stable linear system. Remind that σ = The controller has the following form:

uk =

K3 e + K4

where K3 and K4 are positive constant.

edt + β (vr − e) + v˙r + δ b(Vh + e − vr )

R t+Ts t

R t+Ts t

edt.

(7.16)

7.2. Control laws

90

Note that, if K3 and K4 are positive and the next Lyapunov function is chosen: ( e2 Vlin = + K4 2

R t+Ts t

edt)2 , 2

which differentiation is V˙lin = −K3 e2 ≤ 0, then the equilibrium stability with control (7.16) is guaranteed. For physical implementation, the controller (7.16) has the next discrete-time approximation, ¯  K3 (ek − ek−1 ) + K¯ 4 ek + β Ts (vrk − ek ) + vrk − vrk−1 + δ Ts N uk = sat1 round (7.17) bTs (Vh + ek − vrk ) where K¯ 3 and K¯ 4 follow the similar change of parameter given for K¯ 1 and K¯ 2 in Controller No.1 (Eqs (7.14)–(7.15)), i.e.: K4 K¯ 3 , K3 − 2 ¯ K4 , K4 Ts .

(7.18) (7.19)

Note that, the saturation and rounding functions are necessarily considered in the discretetime controller. Figure 7.7 shows the implementation of the approximate discrete-time controller No. 2, Eq. (7.17). The performance of the control (7.17) with K3 = 2.4 and K4 = 1.44 is displayed in Fig. 7.8. These constants have been tuned by ensuring that the closed-loop system is Hurwitz. Observe that the current peaks have been reduced with respect to the ‘intuitive controller’, obtaining a smoother current signal. Moreover, the voltage evolves towards the voltage reference with hardly oscillations. As a consequence, the dissipated energy will be reduced, as will be shown below. Nevertheless, this controller has two drawbacks, it directly uses model parameters and needs a larger number of computational blocks.

7.2.4 Controller No.3: Lyapunov-based design The last controller is designed guarantying closed-loop Lyapunov stability conditions for the equilibrium, e = 0. Once again the design is performed employing the continuous-time error equation (7.8).

Chapter 7. High-performance control for the DC-DC Vdd-Hopping converter

91

CONTROL z −1

vrk−1 + − +

+

βTs

+



δTs

ek

+

vc

+

¯4 K

vr

+

+

¯3 K





z −1

× ÷

+

+

round

satN 1

uk

ek−1

Vh

+ + −

bTs

Figure 7.7: Digital feedback linearization control.

NTrans

a) 20 10 0 0

0.2

0.4

0.6

0.8

V(V)

t(s) b) 1.2 1 0.8 0

0.2

0.4

x 10

0.6

0.8

I(A)

t(s) c) 0.05 0 0

0.2

0.4

1 −6

x 10

0.6 t(s)

1 −6

0.8

1 −6

x 10

Figure 7.8: Control No. 2. Evolution of: a) number of PMOS transistors switched on, b) the vr (dashed) and vc (solid). c) the current Il .

7.2. Control laws

92

Consider the following Lyapunov function candidate Vlyap =

e2 . 2

Its time derivative is V˙lyap = −β e2 + (b(vr −Vh )uk − buk e + β vr + δ + v˙r )e.

(7.20)

The negativeness of V˙ can be assured canceling the undesired terms. This can be done by choosing β vr + v˙r + δ uk = , (7.21) b(Vh + e − vr ) then Eq. (7.20) will be V˙lyap = −β e2 ≤ 0.

Therefore, e = 0 is asymptotically stable.

The approximate discrete-time version of Eq. (7.21) considering the saturation and rounding function for physical implementation purposes is:   β Ts vrk + vrk − vrk−1 + δ Ts N uk = sat1 round , (7.22) bTs (Vh + ek − vrk ) Figure 7.9 shows a block diagram of this discrete-time controller. CONTROL

z −1

vrk−1 + − +

δTs

+

vr +

vc

+

βTs

ek

+

× ÷

round

satN 1

uk



Vh

+ + −

bTs

Figure 7.9: Digital Lyapunov control.

The performance of this controller is displayed by simulation in Fig. 7.10. Note that, the application of this controller to the Vdd-Hopping converter reduces the current peaks,

Chapter 7. High-performance control for the DC-DC Vdd-Hopping converter

93

obtaining smoother voltage and current evolutions. Note that the obtained performance is similar to the one obtained with the feedback linearization controller (see Fig. 7.8). However, this controller presents less computational blocks. Compare Fig 7.7 and Fig. 7.9.

NTrans

a) 20 10 0 0

0.2

0.4

0.6

0.8

V(V)

t(s) b) 1.2 1 0.8 0

0.2

0.4

x 10

0.6

0.8

I(A)

t(s) c) 0.05 0 0

0.2

0.4

1 −6

x 10

0.6 t(s)

1 −6

0.8

1 −6

x 10

Figure 7.10: Control No. 3. Evolution of: a) number of PMOS transistors switched on, b) the vr (dashed) and vc (solid). c) the current Il .

7.3 Performance evaluation In this section a performance evaluation is performed for the resulting voltage and current signals, after applying the previous controllers. The voltage signal performance is evaluated by computing the mean and variance of the voltage error. Likewise, the current signal performance is evaluated by computing the maximum current peaks produced as well as its Power Spectral Density (PSD). This PSD is computed using all the recorded data, since this decomposition is computed after the simulation. Table 7.1 presents the mean and variance of the voltage error signal and maximum peak of the current signal.

Note that, all the new proposed controllers improve the system performance with respect

7.3. Performance evaluation

94

Intuitive Contr. No. 1 Contr. No. 2 Contr. No. 3

Mean Error 3.32 · 10−3 2.4 · 10−3 2.52 · 10−3 2.24 · 10−3

Var. Error 6.59 · 10−5 2.54 · 10−5 4.66 · 10−5 3.37 · 10−5

Max. Curr. Peak 4.0 · 10−2 2.5 · 10−2 0.5 · 10−2 0.5 · 10−2

Table 7.1: Performance evaluation.

to the solution given in [99]. Furthermore, equilibrium stability has been guarantied for Controller No.2 and No.3 in the previous section. Likewise, equilibrium stability of Controller No.1 will be proved in a chapter dedicated to such purpose (Chapter 9) because of its complexity. From this point of view, among these new proposed controllers, the most interesting one is Controller No.3, since it provides the best voltage and current performance. This can be observed in Table 7.1 and in Fig. 7.11. Observe that controller No.3 ((d) in Fig. 7.11) provides a PSD smaller than the other controllers (see (a), (b) and (c) in Fig. 7.11).

7.3.1 Energy evaluation In the set of PMOS, the accumulated dissipated energy in the transient period depends on the control law employed, i.e., on the switching sequence. For instance, undesirable oscillatory current profile can be obtained with certain controllers. This non-smooth behavior of the transient current may result in a higher energy consumption. The purpose of this subsection is to evaluate the energy cost associated with each one of the controllers presented in previous section. The estimation of the dissipated energy in the PMOS transistors during the transientperiod is Ed =

Z tf t0

(Vh − vc )Il dt

where t0 is the initial time and t f is the final time in such transient period. Figure 7.12 and Tab. 7.2 show the dissipated energy during the rising transient period. Note that the energy consumption for all controllers presented above is improved with respect to the intuitive controller. This is due to the smoother behavior of voltage and current signals obtained with these controllers. Likewise, note that, the smallest dissipated energy is achieved with Controller No.3, since it provides the best performance.

−80

−80

−90

−90 Power/frequency(dB/Hz)

Power/frequency(dB/Hz)

Chapter 7. High-performance control for the DC-DC Vdd-Hopping converter

−100 −110 −120 −130

−100 −110 −120 −130

−140 −150 0

5

10

15 20 25 Frequency (MHz)

30

−140 0

35

5

10

15 20 25 Frequency (MHz)

30

35

(b) PI controller

−80

−80

−90

−90

Power/frequency(dB/Hz)

Power/frequency(dB/Hz)

(a) Intuitive controller

−100 −110 −120 −130 −140 0

95

−100 −110 −120 −130

5

10

15 20 25 Frequency (MHz)

30

35

(c) Feedback linearization controller

−140 0

5

10

15 20 25 Frequency (MHz)

30

35

(d) Lyapunov Controller

Figure 7.11: Power spectral density.

DISSIPATED TOTAL ENERGY (µ J) Intuitive control 7.2 Controller No.1 6.8 Controller No.2 6.2 Controller No.3 4.8

Table 7.2: Total energy dissipated in rising transient period.

7.3.2 Summary The intuitive control proposed in [99] provides a reasonable tracking at the expense of an oscillatory behavior due to its own limitation. This involves that the current signal time profile presents a high frequency behavior with some substantial peaks, in particular when the total PMOS parallel resistances are larger. This seems to be the main cause of a larger dissipated energy.

7.3. Performance evaluation

96 −6

x 10 7

Intuitive

6

Control No.1

Energy

5 Control No.2

4 3

Control No.3

2 1 0 0.5

1

1.5

2 t(s)

2.5

3 −7

x 10

Figure 7.12: Energy dissipated during the rising transient period.

Linear control (Controller No.1) does not need model knowledge. This controller also reduces the current peaks with respect to the intuitive controller. The dissipated energy reduction according to this intuitive controller is 5%. Control by linearization (Controller No.2) yields a smother current and voltage timeprofiles, reducing the current peaks. However, it directly needs system knowledge and presents more computational blocks. In terms of energy consumption, this controller improves the ‘intuitive control’ by 14%. Lyapunov’s controller (Controller No.3) requires also model knowledge and a certain number of computational blocks. The highlight of this controller is its energy consumption reduction, which is due to the smoother behavior of the voltage and current time profiles. This involves that the controller reduces by 32% the energy consumption with respect to the ‘intuitive control’. Although the Lyapunovs controller presents very nice characteristics, they may be enhanced by changing the voltage signal reference and adapting the resistive load parameter. Firstly, a signal reference can be computed looking for minimizing the dissipated energy and the current peaks as well. And secondly, note that, the controller depends directly on the resistive load parameter. Therefore, an adaptive controller can be designed to cope with variations and/or uncertainties on this load parameter. These two issues will be seen in the next section. Generally, the power consumption can be reduced if the local core voltage or/and the

Chapter 7. High-performance control for the DC-DC Vdd-Hopping converter

97

clock frequency are decreased, that is why a GALS systems is developed in the ARAVIS project where the Vdd-Hopping converter is embedded. The frequency is chosen in such a way that the task (for instance, the execution of the control laws designed here) is performed before a deadline and that the minimum required local clock frequency that guaranties the critical path (longest path delay) on the corresponding clock domain [45]. Normally, in order to take care of this issue when the frequency and voltage have to rise, firstly the voltage is rising and later the frequency is rising. On the contrary, when the frequency and voltage have to fall, firstly the frequency is fallen, and later the voltage is fallen.

7.4 Advanced Lyapunov’s controller The Lyapunov’s controller (Eq. (7.21)) presents very suited properties for the Vdd-Hopping converter. However, this controller can be improved. On the one hand, minimization of energy consumption and current peaks are desired. This can be achieved finding an appropriate evolution for the voltage reference, v∗c (t), by applying optimal control theory [80, 86, 155]. On the other hand, note that the Lyapunov’s controller depends on the resistance load parameter, β . However, this parameter is, in many occasions, difficult to estimate and may change with time, as mentioned in Section 7.1. Therefore, a second objective is to design an adaptation law in order to obtain an estimation βˆ for the unknown parameter. The proposed control architecture including the optimal reference and the adaptation mechanism is shown in Fig. 7.13. ADAPTIVE CONTROL

βˆ

OPTIMAL REFERENCE

v∗

+ CONTROL

uk

SYSTEM

vc

-

Figure 7.13: Vdd-Hopping closed-loop with optimal evolution of the reference and adaptation parameter.

7.4. Advanced Lyapunov’s controller

98

7.4.1 Optimal voltage reference computation Assume that the desired voltage is constant. The problem may be formulated as to find a continuous-time voltage reference trajectory from a voltage initial value vc (t0) to set-point vr , minimizing current peaks, ∆I, and the dissipated energy. This problem will be addressed applying continuous-time optimal control theory [80, 86, 155]. In order to optimize the current peaks, time derivative of the current I˙l is included in the performance index. In every sampling time, a certain number of transistors will be switched on. The total number of PMOS transistors switched on at the previous sampling time is denoted by u− k , and the total number of PMOS transistors switched on at the current sampling time is denoted by u+ k . Consequently, the number of PMOS transistors switched on or off in − every sampling time is given by ∆uk = u+ k − uk . The current peaks are due to the sudden change of the PMOS resistance at the sampling times. These peaks ∆Il = Il+ − Il− are given by ∆Il =

Vh − vc + Vh − vc (uk − u− )= ∆uk . k R0 R0

− The same notation given above for u+ k and uk is used here for variable Il . Therefore, the continuous-time approximation for the current peaks is

Vh − vc u. ˙ I˙l ≈ R0 Another way to achieve this same expression is taking time derivative of Il given by Eq. (7.1). Rigorously, the time derivative of this current is Vh − vc v˙c I˙l = u˙ − u. R0 R0 Nevertheless, it can be shown by simulation that during a typical transient-period, the last term is very small (see Fig. 7.14). This simulation is performed using the same parameters given in Section 7.2. This graph supports the previous argument. Take the following performance index J=

Z τ

L(e, u,t)dt,

(7.23)

0

where the final time τ is free and the Lagrangian L(e, u,t) is chosen in order to penalize: • voltage error e,

Chapter 7. High-performance control for the DC-DC Vdd-Hopping converter

99

9

Current time−derivative terms

9

x 10

8 7 6 5 4 3

Vh −vc R0 u˙

2

Vh −vc v˙c R0 u˙ − R0 u

− Rv˙c0 u

1 0 0

1

2

3 t(s)

4

5

6 −9

x 10

Figure 7.14: Current time-derivative terms.

• dissipated power P = (Vh − vc )Il and • current peaks I˙l . For this, the following Lagrangian is chosen 2

L = q1 e + q2



(Vh − vr + e)2 u R0

2



Vh − vr + e + u˙ R0

2

,

(7.24)

where q1 and q2 are positive weighting constants. The first term of Eq. (7.24) penalizes the voltage error, the second one penalizes the dissipated power in the set of transistors and 2 c) the last one, the current peaks. This dissipated power is modeled as Pw = (Vh −v u, where R0 vc = vr − e. Let us consider a 2-dimensional optimal control problem x = [e, u] with x˙ = [e, ˙ ν ], where ν , u. ˙ Thus, the Hamiltonian function is   (Vh − vr + e)4 u2 Vh − vr + e 2 H = q1 e + q2 + ν + λ1 [b(−Vh + vr − e)u + β + δ ] + λ2 ν . R0 R20 (7.25) 2

7.4. Advanced Lyapunov’s controller

100

Solving the algebraic equation ∂ H(e, ν , λ1 , λ2 ) ∗ = 0, ∂ν ν =ν

the optimal ν ∗ (x, λ ) is

−λ2 ν = 2 ∗



R0 Vh − vr + e

2

,

which gives the optimal Hamiltonian expression H ∗ (e, λ1, λ2 ) = q1 e2 + q2

λ22 R20 (Vh − vr + e)4 u2 − + λ1 [b(vr −Vh − e)u + β + δ ]. 2(Vh − vr + e)2 R20 (7.26)

The optimal solution is associated with the set of differential equations:

∂ H∗ = b(vr −Vh − e)u + β + δ = e˙ (7.27) ∂ λ1  2 R0 ∂ H ∗ −λ2 = u˙ = ν (7.28) = ∂ λ2 2 Vh − vr + e ∂ H∗ ( λ 2 R0 ) 2 (Vh − vr + e)3 2 = 2q1 e + 4q2 u + − buλ1 − β λ1 = −λ˙ 1 (7.29) ∂e 2(Vh − vr + e)3 R20 ∂ H ∗ 2q2 u(Vh − vr + e)4 = + bλ1 (vr −Vh − e) = −λ˙ 2 ∂u R20

(7.30)

with the boundary conditions, e(0) = vr − vc (0) e(τ ) = 0 u(0) = number of transistors switched on in t = 0. u(τ ) = number of transistors switched on in t = τ .

(7.31) (7.32) (7.33) (7.34)

and, the transversality condition H ∗ (τ ) = 0.

(7.35)

Note that, this is a nonlinear Boundary Value Problem (BVP) with a transversality condition, since the final time τ is unknown. Solving (7.27)–(7.35), yields e∗ , from which, the optimal voltage evolution v∗c = vr − e∗ can be derived. This evolution can be employed as reference for the controllers developed in Section 7.2. Numerical solution

Chapter 7. High-performance control for the DC-DC Vdd-Hopping converter

101

The problem raised before: finding a solution for (7.27)–(7.35) with (7.31)–(7.34) and (7.35), is a complex problem because it is a nonlinear BVP with a four dimensional character and it has a transversality condition. For this, a numerical solution is proposed. Nevertheless, finding this numerical solution is also an involved task. There is not so many tools that cope with this kind of problems. In this case the Matlab function‘bvp4c’ has been employed. Function bvp4c [130] combines the solution of Initial Value Problem (IVP) for Ordinary Differential Equations (ODEs) and the solution of algebraic equations, being a non-shooting code. The nonlinear algebraic equations are solved iteratively by linearization, providing an initial guess over a mesh and taking into account the boundary conditions. This is due to the fact that can have more than one solution and, thus, a guess for the desired solution must be provided by designers, which includes an initial mesh for this desired solution. Function bvp4c controls the error of the numerical solution and adapts the mesh in every iteration to obtain an accurate numerical solution with a modest number of mesh points. Thus, obtaining an ‘residual’ error is common. If the residual error is small, then the solution provided by function bvp4c is a good solution. Function bvp4c is not directly applicable for the present problem since it cannot handle the transversality condition. Thus, this function has been used iteratively in order to obtain a solution that fulfills condition (7.35). The system parameters given in Section 7.2 are reported. Furthermore, it is considered the rising transient period, i.e, when output voltage goes from the low voltage level to the high voltage level. Therefore, the next boundary conditions are selected: e(0) = vr − vc (0) e(τ ) = 0 u(0) = 1 u(τ ) = N

(7.36) (7.37) (7.38) (7.39)

The following values for the weighting constants are chosen1 , q1 = 0.64 q2 = 0.32.

(7.40) (7.41)

Using as initial guess 8t

e(t) = 0.3e−10

8

u(t) = 24 − 23e−10 t λ1 (t) = 106t + 105 λ2 (t) = −3.2 · 107t 2 − 3.2 · 107t + 100, 1 As

usual in optimal control problems, they have been chosen in a trial-and-error procedure, checking by simulations the solutions obtained.

7.4. Advanced Lyapunov’s controller

102

Note that this initial guess has a complex form and, because it has been difficult to obtain. As mentioned before, this problem is a complex problem, and finding a solution has been very involved. However, with the future numerical methods, it is expected that new tools for this kind of problem will be researched and developed. The nonlinear BVP (7.27)–(7.30) with the specified boundary conditions (7.36)–(7.39) reaches the numerical solution shown in Fig. 7.15. Note that, the boundary conditions in e∗ and u∗ are satisfied when τ = 23.3 · 10−9 s. From e∗ , the optimal evolution of the voltage reference v∗ can be obtained. Figure 7.16 shows the evolution of H ∗ , whose value at τ = 23.3 · 10−9s is close to zero, fulfilling the transversality condition. Notice that, this voltage reference has been computed for the rising transient period. For the falling transient period, a similar procedure can be applied.

a)

b)

u∗ (NTrans)

e∗ (V )

0.4 0.3 0.2

25 20 15 10

0.1

5

0 0

0.5

1.5

2

0 0

2.5

0.5

−8

x 10

c)

4

8

1 t(s)

x 10

1 t(s)

1.5

2

2.5 −8

x 10

d) 1000 0

6

−1000

λ2∗

λ1∗

4

−2000

2 0 0

−3000

0.5

1 t(s)

1.5

2

2.5 −8

x 10

−4000 0

0.5

1 t(s)

1.5

2

2.5 −8

x 10

Figure 7.15: Optimal numerical solution. a) error evolution, b) control evolution, c) λ1 evolution and d) λ2 evolution.

Chapter 7. High-performance control for the DC-DC Vdd-Hopping converter

103

12

3.5

x 10

3 2.5

H∗

2 1.5 1 0.5 0 −0.5 0

0.5

1

1.5

t(s)

2

2.5 −8

x 10

Figure 7.16: H ∗ evolution.

7.4.2 Adaptive feedback control design The Lyapunov’s controller (Section 7.2) has been designed under the assumption that the parameter β is known. In this section, an adaptive law is proposed in order to cope with the case when the load parameter β is unknown. Let us denote βˆ as the estimated value for the load parameter. This estimated parameter will be used in control law (7.21) instead of its real value. The application of this law to system (7.8) yields e˙ = −β e + β vr − βˆ vr . (7.42) Let us assume that β is a constant parameter which involves β˙ = 0 (the case when β is time-varying will be discussed in next section) and define

β˜ = β − βˆ ,

˙ β˙˜ = −βˆ .

For the adaptive control system, the next Lyapunov function candidate is proposed W=

e2 β˜ 2 + , 2 2 γ1

where γ1 is a positive design parameter that may define the adaptation speed.

(7.43)

7.5. Simulation of the advanced Lyapunov’s controller

104

Differentiating W with respect to time, yields ˙ = −β e2 + β˜ W

β˙˜ vr e + γ1

!

.

Note that β > 0, as has been seen above. The adaptive law is designed by canceling the term in brackets, i.e.: ˙ βˆ = −β˙˜ = γ1 vr e.

(7.44)

˙ = −β e2 . This achieves W Asymptotic stability is established by LaSalle’s invariance principle [76]. For this, consider the level set Wc = W (e, β˜ ) ≤ c0 for sufficiently large c0 > 0, where W˙ ≤ 0. This set is compact and positively invariant. ˙ = 0 on e = 0. Furthermore, note from Eq. (7.42), that Note that W e(t) ≡ 0



e(t) ˙ ≡0



β˜ (t) ≡ 0.

˙ = 0 corresponds to the single point Therefore, the maximum invariant set in Wc with W ˜ P1 = (e = 0, β = 0), thus, every solution starting in Wc approaches the desired point P1 as t → ∞.

7.5 Simulation of the advanced Lyapunov’s controller In this section, some simulations using the parameter presented in Section 7.4 are performed. The resulting controller (7.22) is

βˆ Ts v∗rk + v∗rk − v∗rk−1 + δ Ts N uk = sat1 round bTs (Vh + ek − v∗rk ) (

)

(7.45)

where v∗rk and v∗rk−1 comes from the discretization of the optimal voltage reference, which has been previously obtained. For implementation, the values of v∗rk can be stored in a table. In the same way, βˆ is adapted by the discrete-time approximation of the adaptation law (7.44):

β˜k = β˜k−1 − Ts γ1 v∗rk ek

Chapter 7. High-performance control for the DC-DC Vdd-Hopping converter

105

The data reported in Section 7.2 are used in the simulations. In order to perform more realistic tests, a more precise model for the load is considered in such a way that β depends on rL , i.e., it is time-varying. The bounds on β are: βmin = 1.38 · 107 for vc = Vl and βmax = 5.9 · 107 for vc = Vh . As initial estimated values is taken: βˆ = 0. Figure 7.17 shows the closed-loop performance by employing the optimal voltage reference and the adaptation mechanism. Note that when the adaptation mechanism is implemented the system can achieve a similar performance to the case of known load. Although, the adaptive control introduces a delay in the system response, small current peaks and faster transient periods are obtained. In addition, thanks to the obtained voltage reference the energy consumption is reduced. The accumulated dissipated energy in the rising transient period is 4.8µ J using the Lyapunov controller (Eq. (7.22)) and 0.5µ J using the advanced Lyapunov controller (Eq. (7.45)). That means 90% energy saving. Likewise, the transient period is reduced to 23.3ns.

NTrans

a) 20 10 0 0

0.2

0.4

0.6

0.8

V(V)

t(s) b) 1.2 1 0.8 0

0.2

0.4

x 10

0.6

0.8

I(A)

t(s) c) 0.05 0 0

0.2

0.4

1 −6

x 10

0.6 t(s)

1 −6

0.8

1 −6

x 10

Figure 7.17: Vdd-Hopping with control (7.45) and adaptation βˆ . Evolution of a) number of PMOS transistors switched on, b) vr (dashed) and vc (solid) and c) current Il .

The adaptation of the load resistive component β is shown in Fig. 7.18. Note that β approaches its real value, in spite of the fact that β is time-varying. Observe that the timeevolutions of βˆ and β are superimposed. Consequently, the reliability and efficiency of the controller (7.45), which uses the optimal voltage reference and the adaptation mechanism, has been validated. Besides, the fact,

7.6. Conclusion

106 7

6

x 10

β ,βˆ

4

2

0 0

1

2

3

t(s)

4

5

6 −7

x 10

Figure 7.18: Time-evolution of βˆ (solid) and β (dashed).

that there exists a time-varying load resistive component is not relevant for the right system performance. In addition, this controller has a small energy consumption as well as small current peaks and faster transient-periods.

7.6 Conclusion In this work, a set of controllers has been designed for a Vdd-Hopping converter. Most of these controllers improve performance over the one used in [99] in terms of transient responses, as has been seen in Section 10.5. This controller is a very simple controller with a strong limitation: only one transistor can be switched on or off in every sampling time. The good results obtained with the set of controllers developed in this chapter come from applying control theory as well as the possibility to let such controllers to switch more than one transistor at once. In a performance evaluation presented in Section 10.5 to the set of controllers, it has been concluded that the best one seems to be Lyapunov’s controller. This is not only for a better signal performance but also for a better energy consumption. Nevertheless, this controller can be enhanced, if both optimal and adaptive control are developed. These control approaches allow to diminish energy consumption and current peaks and deal with unknown load parameters, respectively. A method to obtain an optimal reference has been developed applying optimal control theory [80, 86, 155]. Nevertheless, the problem stated for this method presents a high complexity, because it is a BVP with transversally condition for a 4th -order optimal problem. Hence, that an optimal reference has been computed from of this problem. This numerical solution has been obtained using the Matlab function bvp4c. It has been a involved task, and

Chapter 7. High-performance control for the DC-DC Vdd-Hopping converter

107

it is expected that new mathematic tools will be developed to make easier to compute this voltage reference. This result achieves a reduction of current peaks and 90% energy saving with respect to the previous Lyapunov controller. This fact makes that the total energy saving with respect to the intuitive controller used in [99] is 93% reduction. In addition, an adaptive strategy is developed in order to deal with the load modeling error. Moreover, in order to prove the reliability of this adaptive controller, it has been introduced by simulation that parameter β is time-varying, as is common in practice. The suited performance of the results in the Vdd-Hopping converter has been shown by simulations. In summary, an advanced controller which does not need knowledge of the load resistive parameter has been obtained. This controller reduces energy consumption as well as current peaks and transient-periods. Nevertheless, it presents a more complex physic implementation, since it requires more computational blocks. Next chapter, a new controller will be developed to cover this drawback, at the same time that all nice achieved properties are maintained.

108

7.6. Conclusion

Chapter 8 Energy-aware controller for the Vdd-Hopping converter In the previous chapter a set of high-performance controller has been designed for the DCDC Vdd-Hopping converter. From the set of controllers, the one that provides a better performance has been selected and some developments have been done to fulfill with low-power technology requirements. Nevertheless, it presents some number of computational blocks, what can be translated in a complex implementation in certain industrial applications, as in the ARAVIS project. That is why the controller with presents the less number the computational blocks from the set of controller proposed before is chosen and enhanced in order to accomplish all control objectives. It is based on a PI structure [84, 120]. From this controller, an optimal nonlinear energy-aware controller is obtained. The proposed solution is a discrete-time control mechanism, which does not need to track any timeindexed voltage reference. This control law only needs to know the set-point. As an important innovation, the proposed a control introduces a saturation with time-varying limits, which reduces the current peaks. These facts involve an important diminishing of energy consumption. Moreover, its computing cost have been reduced. It is patented under the name of ENergy-AwaRe Control (ENARC) [6]. In the ARAVIS project, this controller will be simulated in VHDL-AMS. In addition, it is expected that it is implemented in the innovated 45nm or/and 36nm SoC developed in ARAVIS.

109

8.1. Control design without current-peak managing

110

8.1 Control design without current-peak managing In this section, a controller based on linear control theory is designed to cope with possible steady-state errors [105]. The relevant characteristic of this controller is that presents a relative low number of computational blocks. This makes it interesting in industry. Previously, it has been enhanced assuming that more than one transistor is switched at once and, dealing with a tracking problem. Now, let us consider a stabilization point problem. Therefore, the reference will be constant, i.e., a step. This looks for making faster transient periods. Nevertheless, important current peaks can be generated being able to damage the system. The ‘linear control’ proposed in Section 7.2 (Eq. (7.10)) has been: uk = satN 1 {round(K1 e + K2 σ )} ,

(8.1)

where σ corresponds to tt+Ts edt, which can be considered as a new variable that augmentes the system dimension. Constants K1 and K2 are chosen with the tunning mechanism (7.11)– (7.12). R

Figure 8.1 shows a simulation to display the stability properties of this control law. The data reported in Section 7.2 are employed in this simulation. The time-varying reference is changed by a step reference, what makes to obtain faster transient periods. Nevertheless, the constant reference generates a no-desired important current-peak (as has been predicted before), which may increase the dissipated energy as well as can damage the physical system. Note that this current peak is not symmetric when system is falling down. It is due to the current definition (Eq. (7.1) in Section 7.1). Observe that current variable directly depends on the number of PMOS transistors switched on, uk , that is, smaller current for lower voltage.

8.2 Control redesign with current-peak managing The controller presented before is modified in order to achieve a high-performance from a point of view of current-peaks. Current-peaks can be managed by introducing a pre-specified maximum admissible current-peak constraint, e.g., introducing an on-line saturation mechanism. This current-peak constraint is defined in the Vdd-Hopping system by ∆Ilmax . In Section 7.4, an expression that relates the current variation with the number of transistor switched on or off in every sampling time has been given. In order to make easier the reading, this expression is here advocated ∆Il =

Vh − vc ∆uk . R0

(8.2)

Chapter 8. Energy-aware controller for the Vdd-Hopping converter

111

NTrans

a) 20 10 0 0

0.2

0.4

0.6

0.8

V(V)

t(s) b) 1.2 1 0.8 0

0.2

0.4

x 10

0.6

0.8

I(A)

t(s) c) 0.1 0 0

0.2

0.4

1 −6

1 −6

x 10

0.6

0.8

t(s)

1 −6

x 10

Figure 8.1: Control (8.1) with step reference. Evolution of a) number of PMOS transistors switched on, b) vr (dashed) and vc (solid) and c) current Il .

It comes from the discrete property of the Vdd-Hopping converter, while that its core voltage is a continuous-time variable, as is seen in [99]. Assume that a maximum admissible current variation for Eq. (8.2) is defined as ∆Ilmax , in such a way that Vh − vc ∆Ilmax ≤ ∆uk ≤ ∆Ilmax . (8.3) R0 This constraint can be introduced in the system by saturating the maximum and minimum PMOS transistors that can be switched on or off and every sampling time, i.e., the maximum in minimum admissible ∆uk are R0 ∆Ilmax , αkM Vh − vc R0 ∆um ∆Il , αkm , k =− Vh − vc max

∆uM k =

(8.4) (8.5)

being αkM > 0 and αkm < 0. Note that the saturation limits depends on the output voltage. Thus, control (8.1) is modified in order to consider current-peaks uk = satN 1



α¯ M round(satα¯ km (K1 e + K2 σ )) k



,

(8.6)

8.2. Control redesign with current-peak managing

112

where the current-peak constraints are employed according to next expressions:

α¯ kM , uk−1 + αkM α¯ km , uk−1 + αkm .

(8.7) (8.8)

Remark 8.1 The current peak constraint must allow to switch, at least, one transistor in order to guarantee that the system achieves the desired voltage level. Therefore, the admissible current peak constraint has to be larger than a minimum bound guaranteeing this condition. From Eq. (8.3), it can be obtained ∆Ilmax ≥

Vc −Vlow Vch −Vlow min(∆uk ) = h 1 R0 R0

where min(∆uk ) is the minimum number of PMOS transistors that can be switched on or off in every sampling time. As uk is determined by a rounding function, the minimum of (∆uk ) must be larger or equal than 1 in order to have that the constraint ∆Ilmax allows switching at least one transistor.

The aware management of the current-peaks makes Control (8.6) to be an innovated controller for the Vdd-Hopping converter. As a side effect, an important reduction of the dissipated energy is achieved by means of a trade-off between faster transient period and small current-peak as will be shown below. The optimization of the energy dissipation is a crucial point in the miniaturization of microsystems.

8.2.1 Time discretization Usually, the controllers are implemented in discrete-time. That is why control (8.6) is approximately discretized, yielding the structure   αkM N ¯ ¯ uk = sat1 uk−1 + round(satα m (K1 ∆ek + K2 ek )) , (8.9) k

where K¯ 1 and K¯ 2 follows the transformation given in (7.14)–(7.15) in Section 7.2, and where αkm and αkM are given in Eqs. (8.4)–(8.5). This controller is patent pending under the name of ENergy-AwaRe Controller (ENARC) [6]. The structure of the ENARC controller is shown in Fig. 8.2. This controller is composed of:

Chapter 8. Energy-aware controller for the Vdd-Hopping converter

113

• two gains, K¯ 1 and K¯ 2 , adjusted by a tuning method given by Eqs (7.11)–(7.12), in such a way that asymptotic stability of the equilibrium of the closed-loop system is guaranteed. The stability will be analyzed in Chapter 9. • A current limit mechanism, that computes on-line the maximum and minimum number m of PMOS transistors switched in every sampling time, ∆uM k and ∆uk , respectively. The limits over the switched-transistor variation are necessary in order to ensure that a maximum admissible current peak ∆Ilmax is respected. • A saturation mechanism that limits the maximum and minimum switched-transistor variation computed before. • A rounding for digital control signal and • an output saturation mechanism, that limits the minimum and maximum PMOS transistor number switched on at every sampling time.

1. Error Tracking Filter vr ADC vc

ADC

vrk vck

¯2 K ek

¯1 K Memory

ek−1

∆uk

3. On-line saturation mechanisms

∆usk 4. Rounding mechanism ∆¯ usk

∆uM k

∆um k

5.Output saturation mechanism u ¯sk−1

∆Imax

2. On-line current limits mechanism

Memory

u ¯sk u ¯sk

satN 1

Figure 8.2: ENARC structure patented in [6].

8.2.2 Simulation of ENARC controller in the Vdd-Hopping system Now, some simulations are performed to display system signal evolutions. The reference signal is a step, as has been chosen in Section 8.1. For implementation, some values of αkm and αkM can be stored in a table. For this, the reported parameter values from Section 7.2 are chosen. Consequently, from Eqs. (7.14)–(7.15) and Eqs. (7.11)–(7.12) K¯ 1 = −19.3 K¯ 2 = 39.27.

(8.10) (8.11)

8.2. Control redesign with current-peak managing

114

The maximum admissible current peak constraint is chosen according to the equation given in Remark 8.1 Vc −Vl ∆Il = h 0.6. R0 Figure 8.3 shows a simulation of the resultant closed-loop system with Eq.(8.9). Note, that the system response obtains faster transient period and reduces considerably the current peak with respect to the previous controller (see Fig. 8.1).

NTrans

a) 20 10 0 0

0.2

0.4

0.6

0.8

V(V)

t(s) b) 1.2 1 0.8 0

0.2

0.4

x 10

0.6

0.8

I(A)

t(s) d) 0.05 0 0

0.2

0.4

1 −6

x 10

0.6 t(s)

1 −6

0.8

1 −6

x 10

Figure 8.3: ENARC with step reference. Evolution of the: a) number of PMOS transistors switched on, b) vr (dashed) and vc (solid), c) current. Il .

Energy-consumption evaluation. Now, let us discuss the main objective concerning to the controller: the energy consumption. The cumulated dissipated energy in the set of PMOS during the rising transient time depends on the type of control law employed, and thus, on the switching sequence. Remember that in the control structure published in [99] only one transistor can be switched at each sampling time and a ramp reference has been employed. However, with the ENARC controller more than one transistor may be switched in every sampling-time and a step reference

Chapter 8. Energy-aware controller for the Vdd-Hopping converter

115

has been employed. The energy consumption during the rising transient time has been reduced from 7.2µ J (with the control published in [99]) to 0.26µ J (with ENARC), i.e., 96% reduction. The energy consumption is much higher using the ‘intuitive’ controller than using the ENARC controller. Note, thus that, a non-smooth behavior of the current signal and a larger transient time may result in a higher energy consumption.

8.3 Conclusions In this chapter, some important advances have been developed for the ‘linear controller’ (Controller No.1, in Section 7.2) proposed for the Vdd-Hopping system. The ‘linear controller’ has a relevant interest due to its simple implementation, since it requires a relative low number of computational blocks. A nonlinear discrete-time controller has been designed for the Vdd-Hopping system with the aim of reducing the dissipated energy. This controller has an energy-aware management of current-peaks in the set of PMOS transistors. In addition, a step reference has been used, thus, it only needs to know the two desired voltage levels. This result comes from the possibility to control more than one transistor at once, i.e., to switch more than one transistor in a same sampling time. As a side effect the transient-periods are diminished. This controller has been compared with the ‘intuitive’ control used in [99]. In this context, the ENARC controller reduces the energy consumption a 96%. The ENARC controller does not only have the same properties that the high-performance Lyapunov controller developed in Chapter 7; but it also has a simple implementation. And thus, it may present special interest for industrial applications in the fields of microelectronics. This mechanism is an innovative controller for the discrete Vdd-Hopping converter applied in the ARAVIS project. It is focused on achieving the project global objective. The features that have not yet been studied are the presence of delay as well as parameter uncertainty, which will be taken into account in Chapter 10. In addition, the closed-loop stability of Vdd-Hopping with ENARC will be studied in next chapter.

116

8.3. Conclusions

Chapter 9 Approximate stability analysis of the DC-DC Vdd-Hopping converter In Chapter 8, a controller has been developed based on a linear structure. This is a nonlinear discrete-time energy-aware controller called ENARC, which may fulfill most of the specific control objectives for the Vdd-Hopping converter, mentioned in Chapter 6. Among them its low computational cost can be highlighted. This characteristic makes it attractive for industrial applications [17, 84, 120]. However, in that chapter, the equilibrium stability of the system with the ENARC has been not studied. Therefore, the convergence and stability of the equilibrium of the closed-loop Vdd-Hopping is not reliable. This chapter focuses on studying the stability of the nonlinear system (7.7) with the ENARC controller developed in Chapter 8. As Vdd-Hopping model is continuous, the stability analysis is performed in continuous-time, assuming that the ENARC stability property is ensured through its continuous-time version (this continuous-time version of the ENARC has been presented in control (8.6)). This is a very common assumption [71, 152]. For simplicity, a preliminary equilibrium stability analysis of the closed-loop system is studied when the controller does not consider the current peak issues. This controller has been the ‘linear controller’ presented in Section 7.2 (Eq. (7.10)). Then, a stability analysis is performed when control (8.6) is employed. This controller introduces a type of nonlinearity: a saturation with dynamic limits. This fact makes that the nonlinear closed-loop system works in three operating-modes: non-saturated system, saturated system in the upper limit and saturated system in the lower limit. It is seen that the equilibrium is in non-saturated mode. The stability analysis is based on LaSalle’s invariance principle [76]. The stability analysis presented here is not rigorous because it is based on some system approximations in continuous time. The complexity of the discrete system with the ENARC 117

9.1. Stability with control (7.10)

118

does not allow to find another simple way to prove the global stability.

9.1 Stability with control (7.10) In this section, a stability analysis is performed for the nonlinear model (7.7) of the VddHopping converter with control (7.10) developed in Section 8.1. For simplicity, the saturation that limits the number of the PMOS transistors switched on as well as the rounding function are disregarded. This analysis is based on LaSalle’s invariance principle [76]. To help the reading of the thesis, the equation of the control is recalled, uk = K1 e + K2 σ

(9.1)

Remember that σ = t edt and the control parameters K1 and K2 are constant and chosen according to the tuning mechanism (7.11)–(7.12) described in Chapter 7. These tuning equations depend on a design parameter ξ . The following lemma deals with tuning the parameter ξ , such that, K1 and K2 have certain desired properties. R t+Ts

Lemma 9.1 Consider the interval  ukl b + β I= , 2ωn

 ukl b + β ωn . + 2ωn 2(bukh + β )

If ξ is chosen in I, then the following inequalities are satisfied K1 > 0

and

K1 (bukh + β ) − K2 < 0,

where ukh is the upper-bound of uk , i.e., ukh = N. Notice that K2 > 0 from Eq. (7.12). Remark 9.2 The term

ωn 2(bukh + β ) is always positive from the parameter properties given in Section 7.1, thus it is clear that the interval I is not empty. Then, there is always a possibility to find the suitable values for the control parameters K1 and K2 . Assume that the tuning of these control parameters is achieved according to the previous rules. The next step concerns the stability analysis of the closed-loop system. The following theorem establishes the global stability of system (7.7) under the control law given in (9.1).

Chapter 9. Approximate stability analysis of the DC-DC Vdd-Hopping converter

119

Theorem 9.3 Consider system (7.7) with the controller (9.1). Then, if K1 and K2 are positive, the equilibrium of the system is globally stable.

Proof: From the dynamics of system (7.7), its equilibrium satisfies the equation 0 = (vr −Vh )bu¯k + β vr + δ . Thus, an expression of the control input at the equilibrium is straightforwardly obtained: u¯k =

β vr + δ = K2 σ¯ , (Vh − vr )b

(9.2)

where v¯c and σ¯ are the equilibrium values of vc and σ , respectively. Note that u¯k corresponds to the saturation bounds 1 and N when the set point is the lower and higher values, respectively, as has been defined in Remark 7.3. Now, the new variable wk , which represents the difference between the control uk considered at any position and at the equilibrium, is introduced: wk , uk − u¯k . (9.3) Rewriting the dynamic of system (7.7) using this new variable, yields: e˙ = −(β + buk )e + (vr −Vh )bwk + (vr −Vh )bu¯k + β vr + δ ,

(9.4)

Substituting the expression of u¯k from (9.2) into the right-hand side of (9.4), the following equation is obtained e˙ = −(β + buk )e + (vr −Vh )bwk . According to the assumption of the theorem, the control constant K2 is positive. Moreover, from Section 7.1, the positivity of β and b is guaranteed. Since the control uk is equal to satN 1 {round(K1 e+K2 σ )}, it is clear that uk is positive and consequently so is β +buk +K2 . The next step of the proof is based on the Lyapunov’s theorem. Consider the Lyapunov function candidate of the form: V (e, σ ) =

e2 (σ − σ¯ )2 + K2 , 2b(Vh − vr ) 2

Notice that, the reference voltage vr is constant. From the system properties, the high voltage Vh is greater than the reference voltage vr . As has been seen just before, b is positive. Thus, V is indeed a positive definite function of e and σ . Thanks to (9.1) and (9.2), V can be expressed as follows e2 (uk − u¯k − K1 e)2 V (e, σ ) = + . 2b(Vh − vr ) 2K2

(9.5)

9.1. Stability with control (7.10)

120

The differentiation of the function V along the trajectories of (9.1) leads to (β + b(u¯k + wk ))e V˙ = − − K1 e2 + b(Vh − vr ) 1 [(K1 e˙ + K2 e) − (K1 e˙ + K2 e)] ·[(K1 e + K2 σ ) − K2 σ¯ − K1 e] . K2 2

Introducing (−K1 e2 + K1 e2 ) in the previous equality and assuming that the differentiation of u˙k from Eq. (9.1) is u˙k = K1 e˙ + K2 e. Then, the differatation of the function can be expressed as:   (β + b(u¯k + wk )) ˙ V =− + K1 e2 . b(Vh − vr )

(9.6)

Note that K1 and K2 have been assumed positive, and uk = u¯k +wk is assumed be positive, thus, V (e, σ ) ≥ 0 V˙ (e, σ ) ≤ 0. Asymptotic stability is established by LaSalle’s invariance principle [76]. For this, consider the level set Ω1 = V (e, σ ) ≤ c1 for sufficiently large c1 > 0. This set is compact and positively invariant. This is represented in Fig, 9.1. From Eq. (9.6), notice that V˙ (e, σ ) is negative everywhere, except on the line e = 0, where V˙ (e, σ ) = 0. Unless σ = σ¯ , this is impossible from the closed-loop system (7.7) with control (9.1), since e(t) ≡ 0



e(t) ˙ ≡0



0 ≡ −(vr −Vh )bK2 σ + β vr + δ .

The last equation is satisfied in

σ≡

vr + δ = σ¯ . (Vh − vr )bK2

Consequently, the maximum invariant set in Ω1 with V˙ (e, σ ) = 0 corresponds to the single point P2 = (e = 0, σ = σ¯ ). Therefore, every solution starting in Ω1 approaches P2 as t → ∞.

In summary, the global stability of the Vd-Hopping with Eq. (9.1) has been established. In next section, this result will be employed to analyze the global stability with control (8.6).

Chapter 9. Approximate stability analysis of the DC-DC Vdd-Hopping converter

121

e Ω1

P2

σ Figure 9.1: Representation of the compact set Ω1 .

9.2 Stability with control (8.6) Based on the previous study, the objective of this section is to extend the proof before to control (8.6). Note that, it is the continuous-time version of the ENARC controller, and it is nonlinear due to control saturations. Here, the saturation that limits the number of the PMOS transistors switched on as well as the rounding function are also disregarded for simplicity reasons. An added difficulty is that the limits of the saturations depend on the state e. The equation of control (8.6) is expressed here, to make easier the reading: α¯ M

uk = satα¯ km (K1 e + K2 σ ), k

(9.7)

with R0 ∆Il Vh − vc max R0 α¯ km = uk−1 + αkm = uk−1 − ∆Il , Vh − vc max

α¯ kM = uk−1 + αkM = uk−1 +

(9.8) (9.9)

as has been defined in Eqs. (8.4)–(8.5) and Eqs. (8.7)–(8.8). The stability analysis comes from dividing the space (e, σ ) in three regions (see Fig 9.2): • Region I, where the system does not saturates.

9.2. Stability with control (8.6)

122

• Region II, where the system saturates in the upper limit. • Region III, where the system saturates in the lower limit. Remark 9.4 The equilibrium of system (7.7), i.e., uk = u¯k , is in Region I.

This is easy to see from Eq.(9.2) and Eqs. (9.8)–(9.9): R0 ∆Imax h −vc R0 K2 σ¯ − V −vc ∆Imax h

K2 σ¯ + V

u¯k = sat

(K2 σ¯ ),

Notice that Vh − vc and R0 are positive as it has been defined in Section 7.2 and that ∆Imax is 0 also positive according to Section 8.2. Then, VhR−v ∆Imax > 0. c The proof will be defined in two parts. In the first part, convergence to the non-saturated region in finite time will be proven. In the second part, convergence to the desired point, once the system is in Region I is proven. For the first part, two properties are performed for the saturated cases (Region II and III), which define the convergence to the non-saturated region (Region I). These properties are based on the variable u˙k , since it is directly affected by the saturation limits: αkm (e) and αkM (e).

Assumption 9.5 For system (7.8) with a suited sampling time, next expression in continuoustime can be taken into account α¯ M

u˙k = satα¯ km (K1 e˙ + K2 e). k

Equation before comes from Eq. (9.7) and Eqs. (9.8)–(9.9) Firstly, next property is defined for Region II:

Property 9.6 In Region II, if K1 is constant and positive, system (7.7) satisfies • u˙k > αkM (e) > 0 • u¨k ≤ εM < 0 being εm , • e>0

K1 ∆Imax C

= constant

(9.10)

Chapter 9. Approximate stability analysis of the DC-DC Vdd-Hopping converter

123

e

Region II

Region I

Region III

σ Figure 9.2: Representation of the system operating regions.

where αkM (e) is defined by Eq. (9.8) and vc , vr − e (as has been seen in Section 7.1). Besides, εM comes from differentiating u˙k and employing the definition of αkM (e) from Eq. (9.8). Parameters C and ∆Imax are positive according to Section 7.2 and Section 8.2, respectively.

Remark 9.7 u˙k is decreasing with time derivative that is bounded away from zero, as is shown in Property 9.6, and hence it will reach αkM (e) (see Eq. 9.8) in finite time.

Secondly, the following property is defined for Region III:

Property 9.8 In Region III, if K1 is constant and positive, system (7.7) satisfies • u˙k < αkm (e) < 0 • u¨k ≥ εm > 0 being εm ,

K1 ∆Imax C

= constant

• e<0 where αkm (e) is defined by Eq. (9.9) and vc , vr −e (as has been seen in Section 7.1). Besides, εm comes from differentiating u˙k and employing the definition of αkm (e) from Eq. (9.9). Parameters C and ∆Imax are positive according to Section 7.2 and Section 8.2, respectively.

9.2. Stability with control (8.6)

124

Remark 9.9 u˙k is creasing, with time derivative that is bounded away from zero, as is shown in Property 9.8, and hence it will reach αkm (e) (see Eq. 9.9) in finite time.

These properties and remarks about the different regions allow to prove the convergence of variable u˙k to the saturation limits. For this, next lemma can be stated

Lemma 9.10 If the next conditions are satisfied: K1 > 0,

K2 > 0

and K1 (bukh + β ) − K2 < 0,

and taking into account the Assumption 9.5 then, system (7.7) with controller (9.7) saturated in the upper or lower saturation limit converges to the non-saturation region in a finite time.

Proof: Firstly, Region II (see Fig.9.3) is studied, i.e., the case when system saturates in the upper saturation limit. Employing u˙k which is affected directly by αkM (e), as has been seen in Eq. (9.10). It is desired to prove that: u˙k → αkM (e)

e

Region II

Region I

Region III

σ Figure 9.3: Region II in the representation of the system operating regions.

Chapter 9. Approximate stability analysis of the DC-DC Vdd-Hopping converter

125

For this, the following Lyapunov function candidate is selected: W = u˙k − αkM (e) = K1 e˙ + K2 e − αkM (e) > 0. Differentiating ˙ = K1 e¨ + K2 e˙ + W

αkM (e) e. ˙ (Vh − vr + e)

Next up, from Eq. (8.4), αkM (e) is differentiated

 αkM (e) W˙ = −b(Vh − vc )K1 u˙k + −K1 (buk + β ) + K2 + e. ˙ Vh − vr + e 

(9.11)

Note that b and (Vh − vc ) are positive from Section 7.2. Besides, from the statement the positivity of K1 is also established. And Property 9.6 maintains u˙k > 0 in this region. Therefore, the first term on the right-hand side is negative for every time instant. For simplicity, this first term is defined as ς (e) , b(Vh − vc )K1 u˙k > 0. (9.12) Rewriting Eq. (9.11), such that the next inequalities are satisfied   (e) α k W˙ = −ς (e) + −K1 (buk + β ) + K2 + e˙ Vh − vr + e   αk (e) K2 e < −ς (e) + K1 (buk + β ) − K2 − Vh − vr + e K1   αk (e) K2 e < −ς (e) + K1 (bukh + β ) − K2 − Vh − vr + e K1 < 0. The first inequality comes from applying u˙k = K1 e˙ + K2 e > 0, thus −e˙ < KK21e . Furthermore, form Property 9.6, it is known that e < 0. A maximum bound of uk , ukh , is taken in the second inequality. From Lemma 9.10, K1 (bukh + β ) − K2 < 0 is fulfilled, then the last inequality is satisfied W˙ < 0.

In addition, from Property 9.6, it is ensured that u˙k converges to the boundary. The proof for Region III (see Fig. 9.4) is symmetric to the one developed before for Region II, applying Property 9.8.

9.2. Stability with control (8.6)

126

e

Region II

Region I Region III

σ

Figure 9.4: Region III in the representation of the system operating regions.

Finally, from the lemma stated before about the convergence of the saturated region to the non-saturated region, the equilibrium localization given in Remark 9.4 and the global stability property of the Region I established in Theorem 9.3, next theorem ensures the global stability property of the system (7.7) with controller (9.7). Theorem 9.11 If K1 and K2 are positive, and K1 (bukh + β ) − K2 < 0, then the equilibrium of system (7.7) with controller (9.7) is globally stable. Proof: There exist two state space regions corresponding to the case when system saturates in the upper limit (Region II), and when the system saturates in the lower limit (Region III). By Lemma 9.10, the system operating in Region II or III converges to the nonsaturated region (Region I) in finite time. In addition, Properties 9.6 and 9.8 guarantee that the system once is in Region I can not cross this saturation lines towards Regions II or III. Now, the proof is concluded by using the assumptions K1 and K2 that are constant and positive, K1 (bukh + β ) − K2 < 0 and by advocating the La Salle’s invariance principle. There exists a set Ω2 limited by the level curve V (e, σ ) = c1 for sufficiently large c1 and the saturation limits, which is compact and positively invariant. Figure 9.5 represents the set Ω2 . From Lemma 9.10 the state of the system reaches Ω2 in finite time. As has been seen in Section 9.1, the maximum invariant set with V˙ (e, σ ) = 0 corresponds to the single point P2 = (e = 0, σ = σ¯ ). LaSalle principle establishes that every system evolution in Ω2 approaches P2 as t → ∞.

Chapter 9. Approximate stability analysis of the DC-DC Vdd-Hopping converter

127

e Ω2

Region II

Region I

P2

Region III

σ Figure 9.5: Representation of the invariant set Ω2 .

9.3 Conclusions In this chapter the equilibrium stability of the nonlinear Vdd-Hopping system when control (9.1) has been analyzed . This controller does not manage the occurrence of current peak in the system. The proof of this stability analysis has been based on LaSalle’s invariance principle. Next up, the analysis has been extended to prove equilibrium stability with the innovative nonlinear controller (9.7). This controller manages current peaks through saturations. The saturation limits depends on the system state, making difficult to prove global stability. The rigorous global stability analysis deals with three operating modes: no saturated, saturated in the upper limit and saturated in the lower limit. It has been ensured that the equilibrium is located in the non-saturated mode. When system is saturated, it will converge to the nonsaturated mode in finite time without being able to return to the non-saturated case. This analysis follows LaSalle’s invariance principle for a domain bounded by Lyapunov level and the saturation lines. The analysis presented here has been performed for a continuous-time version of the ENARC controller. It has been assumed that the stability property is maintained for the discrete-time ENARC controller. This kind of assumptions is very common in control theory

128

9.3. Conclusions

[71, 152]. The only control objectives that have not yet been considered is the equilibrium robustness in view of delays presence and parameter uncertainties that the system can suffer. These issues will be seen in Chapter 10.

Chapter 10 Sub-optimal control considering delays and parameter uncertainties In Chapter 8, an energy-aware controller has been designed for the Vdd-Hopping converter. This converter satisfies control requirements for this low-power converter implemented in the ARAVIS project context. Furthermore, in Chapter 9 a stability analysis of the equilibrium has been performed. In order to cover all control requirements specified in Chapter 6, this chapter focuses on system robustness with respect to delays and parameter uncertainties in the Vdd-Hopping system at the same time that the stability is guaranteed. They are common issues in SoC [112, 140]. Figure 10.1 shows Vdd-Hopping mechanism including delays. The system has a computational h2 -sample-period delay at the control block input required to ensure that the system is synchronized with the cluster clock [77]. Likewise, there is a computational h1 -sampleperiod delay associated with computational issues in the control block output. The size of this last delay depends on a trade-off between power consumption and performance. Generally, the power consumption can be reduced if the local core voltage or/and the clock frequency are decreased. However, this fact produces that the computational speed diminishes, in such a way that the size of the existing delay decreases. Therefore, the h1 -sample-period delay depends on the local clock frequency1 [38]. In many cases, applications do not require the full computational power at any time. The performance requirement is that the task is performed before a deadline. Therefore, it is possible to have a low frequency, which ensures system performance, and hence, to allow reducing the power consumption. On the other hand, there is a minimum required local clock frequency that guaranties the critical path (longest path delay) on the corresponding clock domain circuit [45]. That is why, a low frequency of 200MHz is taken here, for the local clock. This frequency introduces an one-sample-period delay in the control block output. In some cases, it is considered that the delays are fixed and 1 The

higher the clock frequency is, the longer the delay is.

129

130

known. The other relevant issue in low-power technology are the parameter uncertainties, which can generate a non-desirable performance and lack of reliability of the system [34, 93]. In summary, delays presence and parameter uncertainties must be considered in the design of the controller. The control design procedure given in Lemma 9.1 provides a range of values for parameter ξ that guarantees equilibrium asymptotic stability for system (7.7) for fixed parameters. This stability analysis has been performed in Chapter 9. Although this tuning method deals with the current peak management, it does not consider that some system parameters may be uncertain. Likewise, perturbation rejection and delays were not considered either. In order to ensure the robustness of the Vdd-Hopping converter with respect to uncertain parameters and delays presence, a ‘conservative’ sub-optimal tuning approach for K¯ 1 and K¯ 2 is now developed for an approximate ENARC controller version that, for simplicity, disregards the current peak management. For this, the system is rewritten into a suitable state-space form to formulate a robust H∞ problem that can be solved by using Lyapunov-Krasovskii theory [53, 62]. In this process, the saturation in the controller is considered [26, 50]. The designed controller guarantees asymptotic stability, disturbance rejection as well as robustness of the system with respect to delays and uncertain parameters. The problem is expressed in terms of Linear Matrix Inequalities (LMIs). Likewise, an attraction domain is estimated in such a way that a regional stabilization for the saturated control is guaranteed. The robustness properties of the closed-loop system with the sub-optimal control tuning design applied to the ENARC controller are tested by simulations. In these simulations is taken into account that the load dynamic parameter is not constant RL , but is time-varying rL , as mentioned in Chapter 7. Remember that for control purposes, it has been taken as constant in Section 7.1.

vr + -

e

z −h2

CONTROL

K1 K2

uk

z −h1

Vdd-HOPPING

R0

LOAD

+ R ,C L

vc

Figure 10.1: Sub-optimal control tuning for the ENARC controller.

An evaluation of the two control tuning approaches developed in this part of the thesis for the ENARC controller is performed.

Chapter 10. Sub-optimal control considering delays and parameter uncertainties

131

10.1 Problem statement In this section, we are going to formulate the problem rewriting the system (7.8) in an approximate discrete-time form considering both delays and uncertain parameters. This approximate time-discretization is performed by using forward euler method [104], by assuming that the sampling time is small enough to the system evolution. Likewise, for simplicity, an approximate version of the ENARC controller is taken. From Fig. 10.1, note that h1 and h2 are the number of sampling periods in the input and in the output of the control block, respectively. Taking this notation and considering a small sampling period, Ts , the associated approximated discrete-time voltage error equation after employing forward Euler method is ek+1 = (1 − Tsβ )ek + Ts b(vr −Vh )uk−h1 + Ts (β vr + δ ) − bTs uk−h1 ek . Remember that the sampling frequency takes the clock frequency value, Ts = Considering the state xk = [ek

1 ωn .

ek−1 ]T , the applied control law is

uk−h1 = satN 1 {uk−1−h1 + Kxk−h }, where h , h1 + h2 and K = [K¯ 1

(10.1)

(10.2)

K¯ 2 ].

The parameters RL , R0 and C that correspond to β , b and δ and ωn that defines the Ts in model (10.1), can be considered uncertain. Each uncertain parameter is within an uncertainty interval, whose corresponding extremes are • uk ∈ [ukl = 1, ukh = N], M • RL ∈ [Rm L , RL ] , M • R0 ∈ [Rm 0 , R0 ],

• Ts ∈ [Tsm , TsM ]. Remark 10.1 The asymptotic stability of system (10.1) is guaranteed in a point within an uncertainty interval for the low level voltage, Il , and within an uncertainty interval for the high level voltage, Ih , bounded by   m m Rm RL (ukl Vh − RM L (ukl Vh − R0 Kleak ) 0 Kleak ) , (10.3) Il , m M ukl Rm ukl Rm L + R0 L + R0  M  m RM RL (ukhVh − RM L (ukh Vh − R0 Kleak ) 0 Kleak ) Ih , , . (10.4) M m ukh RM ukh RM L + R0 L + R0

10.1. Problem statement

132

Other added problem is that the system can suffer some exogenous disturbances. Consequently, the main objective is to find the optimal gain K in such away that Control (10.2) is robust with respect to delays as well as parameter uncertainties. Likewise, this optimal K must guarantee asymptotic stability and disturbance rejection for the known constant delays, h1 and h2 .

10.1.1 Alternative representation for the saturated control (10.2) and the error equation (10.1) Firstly, some lemmas are given to rewrite the saturated control (10.2) and an alternative form. Lemma 10.2 [67], for xk in R n , assume that there exist K, G ∈ R 1×n such that 1 < uk−1−h1 + Gxk−h < N, then n satN1 {uk−1−h1 + Kxk−h } ∈ Co α (m) (uk−1−h1 + Kxk−h ) + (1 − α (m))(uk−1−h1 + Gxk−h ),

o m = 1, 2, .

Lemma 10.3 Assume that there exists G ∈ R 1×n , P1 2 > 0 ∈ R n×n and c > 0 such that for any xk ∈ X, where  X = xk : xTk P1 xk ≤ c−1 ,

then, 1 < uk−1−h1 + Gxk−h < N, and Control (10.2) admits the following representation 2

uk−h1 =



m=1

h i λm k (α (m)(uk−1−h1 + Kxk−h ) + (1 − α (m))(uk−1−h1 + Gxk−h )

2

=



m=1

h i λm k uk−1−h1 + α (m) Kxk−h + (1 − α (m) )Gxk−h ) ,

being ∑2m=1 λmk = 1, with λmk ≥ 0, for all k > 0. Then, defining Ωα ,

2



m=1

λmk α (m) ,

for

all 0 ≤ λmk ≤ 1,

2

∑ λmk = 1

m=1

where the vertices of the polytope are given by α (m) , Eq. (10.1) can be rewritten ek+1 = (1 − Ts β )ek + Ts b(vr −Vh )(uk−1−h1 + α (m) Kxk−h

+(1 − α (m) )Gxk−h ) + Ts (β vr + δ ) − bTs uk−1−h1 ek

for m = 1, 2. 2P 1

is a positive matrix defined to guarante system stability.

(10.5)

Chapter 10. Sub-optimal control considering delays and parameter uncertainties

133

10.1.2 State-space representation The saturated control and error equation, as redefined before, allow to system (10.1) rewrite it in a state-space form. For this, an expression for the dynamic part of the controller, uk−1−h1 , is obtained. Next expression is achieved from (10.1) uk−h1 =

ek+1 − (1 − Ts β )ek − Ts (β vr + δ ) + bTs uk−h1 ek , Ts b(vr −Vh )

and it is delayed one sampling period uk−1−h1 =

ek − (1 − Ts β )ek−1 − Ts (β vr + δ ) + bTs uk−1−h1 ek−1 . Ts b(vr −Vh )

(10.6)

Now, Eq. (10.6) is applied to (10.5) obtaining ek+1 = (2 − Tsβ )ek − (1 − Ts β )ek−1 + Ts b(vr −Vh )(α (m)K + (1 − α (m) )G)xk−h −Ts b(uk−h1 ek − uk−1−h1 ek−1 ), m = 1, 2. (10.7) This can be rewritten in the following matrix form: (m)

xk+1 = A(uk−h1 , uk−1−h1 )xk + Bu¯k−h where

m = 1, 2,

(10.8)

 2 − Ts β − Ts buk−h1 Ts β − 1 + Ts buk−1−h1 A= , 1 0   Ts b(vr −Vh ) B= , 0 

(m)

being u¯k−h = (α (m) K + (1 − α (m) )G)xk−h for m = 1, 2, with α (m) = [0, 1]. uk−h1 and uk−1−h1 are treated as uncertain parameters, whose values will be inside the uncertainty interval [1, N].

10.1.3 Stability and disturbance rejection problem Equation (10.8) can be rewritten in the following explicit closed-loop form with an L2 perturbation added, in such a way that a H∞ problem can be formulated. xk+1 = Axk + B(α (m) K + (1 − α (m) )G)xk−h + Bw wk xl = φl , ∀l ∈ [−h, 0] zk = I2 xk ,

m = 1, 2,

(10.9) (10.10) (10.11)

10.2. H∞ control design

134

with 

 bw1 1 bw1 2 Bw = , bw2 1 bw2 2 where xk , zk , wk ∈ R n are the state vector, controlled output and exogenous disturbance input, respectively. φk is the initial condition and h ≥ 0 ∈ R is a fixed and known delay. Moreover, Problem 10.4 The problem is to find a X(P1, c), a vector G and K such that,

a) Lemma 10.3 holds and, hence, the closed-loop system (10.8) and b) there exists a Lyapunov-Krasovskii functional Vk > 0, such that Vk+1 − Vk along the solution of (10.9) fulfills Vk+1 −Vk < 0,

(10.12)

when the system is not perturbed, and for any perturbation input, there exists a minimum disturbance attenuation, γ ∗ ≥ 0, such that, for all γ ≥ γ ∗ the L2 gain between the perturbation vector wk , and the output vector zk is less or equal to γ . i.e. kzk k22 − γ 2 kwk k22 < 0, ∀wk ∈ L2 for φl = 0, −h ≤ l ≤ 0.

(10.13)

The solution to this problem guarantees the system stability as well as the disturbance rejection for the time-delay system (10.9)–(10.11).

10.2 H∞ control design In order to cope with this problem a mathematical manipulation of Eq. (10.9) is performed via a descriptor model transformation [51]. The descriptor approach is just a variable change, which makes easier to work with Lyapunov-Krasovskii functional [52].

10.2.1 Descriptor model transformation Equation (10.9) is manipulated in order to achieve the previous objectives. For this, a descriptor model transformation is applied.

Chapter 10. Sub-optimal control considering delays and parameter uncertainties

135

Considering: k−1



ψk ,

yk , xk+1 − xk ,

yi .

i=k−h

Next, Eq. (10.9) is rewritten in the descriptor form [51]:     yk + xk xk+1 = , 0 −yk + Axk − xk + B(α (m) K + (1 − α (m) )G)xk + Bw wk From xk−h = xk − ψk , this system can be compactly written as:     0 0 ¯ E x¯k+1 = Ax¯k − ψk + w, (m) (m) Bw k B(α K + (1 − α )G)

m = 1, 2.

m = 1, 2,

(10.14)

where



 I2 , −I2

I2 A¯ , (m) A + B(α K + (1 − α (m) )G) − I2

m = 1, 2.

E , diag{I2 , 0},

  x x¯k , k , yk

10.2.2 Condition for state-space representation Firstly, the condition a) of the Problem 10.4 is dealt with. From [49], it is seen that, in order (i) (i) to guarantee 1 < uk−1−h1 + Gxk−h < N, for i = 1, 2, where uk−1−h1 = {1, N} ∀x ∈ X given in (10.3), it is necessary that next equations are satisfied (i)

2N ≥ N(1 + cxTk−h P1 xk−h ) ≥ 2(uk−1−h1 + Gxk−h )

(10.15)

2 ≤ (1 + cxTk−h P1 xk−h ) ≤ 2(uk−1−h1 + Gxk−h ),

(10.16)

(i)

which correspond to 

for (10.15) and

for (10.16).



1

" (i)  N − 2uk−1−h1 xTk−h ∗

" (i)  −1 + 2uk−1−h1 1 xTk−h ∗

G cP1

−G cNP1 #

#

1 xk−h



1 xk−h

≥ 0,



≥ 0,

(10.17)

i = 1, 2.

(10.18)

These matrices can be rewritten in a suitable form by employing the Schur’s complement. Likewise defining Y = GQ1 , applying P¯1 = Q1 P1 Q1 and pre and post-multiplying by diag{1, Q1}, next LMIs from (10.17) and (10.18) are obtained " # " # c Y c Y ≥ 0, ≥ 0, i = 1, 2, (10.19) (i) (i) ∗ (N 2 − 2Nuk−1−h1 )P¯1 ∗ (−1 + 2uk−1−h1 )P¯1 respectively.

10.2. H∞ control design

136

10.2.3 Control design Now, the condition b) of the Problem 10.4 can be formulated in terms of Linear Matrix Inequalties (LMIs) [53]. Fulfillment of condition (10.12) plus condition (10.13) is looked for. Take ζ , [x¯k

ψk

wk ]T , then objective (10.13) is satisfied if Vk+1 −Vk + xTk xk − γ 2 wTk wk ≤ ζ T Γ0 ζ < 0.

(10.20)



 P1 P2 Define P , T , being P2 Hermitian. P2 0 For this purpose, as Lyapunov-Krasovskii candidate is considered Vk = V1,k +V2,k +V3,k ,

(10.21)

being V1,k = x¯Tk EPE x¯k , h

V2,k =

k−1

∑ ∑

yTi Ryi ,

n=1 i=k−n k−1 V3,k = xTi Sxi , i=k−h



P1 > 0 R>0

S > 0,

(10.22) (10.23) (10.24)

where V1,k guaranties asymptotic stability of system (10.14) without delays. Delay-dependent as well as delay-independent criteria are considered in V2,k and V3,k , respectively [53, 62]. Next, a sufficient condition for asymptotic stability and disturbance rejection is derived.

Theorem 10.5 Consider system (10.9)–(10.11) with energy-bounded wk and control law u¯k−h = α (m) Kxk−h + (1 − α (m) )Gxk−h for m = 1, 2, where h > 0 ∈ R is a known constant delay and K, G ∈ R 1×n . If there exist S, R, P1 > 0 ∈ R n×n such that the LMIs (10.19) plus the following LMIs are satisfied:  ¯T ¯ A PA − EPE + diag{In , hR} (m) Γ ,  ∗ ∗

−A¯ T P



   0 S + (m) (m) 0 B(α K + (1 − α )G) − 1h R − S ∗

  0 T ¯ A P Bw   <0  0 −γ 2 In (10.25)

for m = 1, 2,, then the equilibrium of the closed-loop system (10.9)–(10.11) is asymptotically stable and there is a value γ ∗ such that for γ < γ ∗ condition (10.13) is fulfilled.

Chapter 10. Sub-optimal control considering delays and parameter uncertainties

137

Proof: The goal is to satisfy Vk+1 − Vk + zTk zk − γ 2 wk T wk < 0 for both disturbance rejection and asymptotic stability of the equilibrium for system (10.14). Lyapunov-Krasovskii method yields: V1,k+1 −V1,k = x¯Tk+1 EPE x¯k+1 − x¯Tk EPE x¯k o n T (m) T T (m) T T T T ¯T = x¯k A − ψk [0 α K B + (1 − α )G B ] + wk [0 Bw ]       0 0 ¯ ×P Ax¯k − (m) ψk + w − x¯Tk EPE x¯k (m) Bw k α BK + (1 − α )BG = x¯Tk [A¯ T PA¯ − EPE]x¯k + ηk + νk , m = 1, 2. where

ηk =

¯ −x¯Tk AP

m = 1, 2.

 0 ψ − ψkT [0 α (m) K T BT + (1 − α (m) )GT BT ]PA¯ x¯k α (m) BK + (1 − α (m) )BG k



νk =

wTk [0

V2,k+1 −V2,k =

BTw ]P



   0 0 T ¯ wk + x¯k AP w + wk [0 BTw ]PA¯ x¯k , Bw Bw k

hyTk Ryk −

h



yTk−n Ryk−n

n=1



x¯Tk

  1 0 0 x¯k − ψkT Rψk 0 hR h

This inequality is obtained developing Eq. (10.23) and applying Jensen Inequality [61]. V3,k+1 −V3,k = xTk Sxk − xTk−h Sxk−h = xTk Sψk + ψkT Sxk − ψkT Sψk These developed expressions are applied to inequality (10.20), in such a way that the LMIs (10.25) are obtained.

10.3 Robust control tuning Now, in this section, the uncertain parameters given in Section 10.1 are taken into account, guarantying the properties achieved above, stability as well as disturbance rejection for the time-delay Vdd-hopping system. That means, to obtain a robust saturated control under

10.3. Robust control tuning

138

parameter uncertainties, satisfying those properties. For this, Theorem 10.5 is extended in the case of polytopic uncertainties. Denote Ω , [A BK and assume that Ω ∈ C o{Ω j , Ω=

Bw

α

uk−h1

uk−1−h1 ]

j = 1, ..., 128}, namely

n

∑ λ jΩ j,

j=1

for

some,

0 ≤ λ j ≤ 1,

and being the vertices of the polytope described by Ω j = [A( j)

n

∑ λj = 1

j=1

B( j) K

( j)

Bw

α ( j)

( j)

( j)

uk−h1

uk−1−h1 ] for j = 1, 2, ..., 128. Pre and post-multiplying LMI (10.25) by Q = diag{Q1 , Q1 , Q1 , Q1 , In } and taking Q1 = P2−1 > 0 and P¯1 = Q1 P1 Q1 , R¯ = Q1 RQ1 , S¯ = Q1 SQ1 , the following sufficient condition is achieved:

Theorem 10.6 Consider system (10.9)–(10.11) with energy-bounded wk and the control law uk−h = Kxk−h where h ≥ 0 ∈ R is a known constant delay and K, G ∈ R 1×n . If there exist −1 n×n for j = ¯ ¯ ¯ T,Y ∈ R n×1 and Q1 ∈ R n×n with K = T Q−1 1 , G = Y Q1 and R, P1 , S > 0 ∈ R 1, ..., 128 such that the LMIs (10.19) and

 ( j) ( j) Γ¯ 1 Γ¯ 2   ∗ P¯1 − 2Q1 + hR¯ Γ¯ ( j) ,   ∗ ∗ ∗ ∗ j = 1, ..., 128.

−α ( j) B( j) T − (1 − α ( j))B( j)Y + S¯ −α ( j) B( j) T − (1 − α ( j))B( j)Y ¯ − Rh − S¯( j) ∗

 Bw Q1  Bw Q1   < 0, (10.26) 0  −γ 2 Q 1

where T T T ( j) Γ¯ 1 , Q1 A( j) + A( j) Q1 − 2Q1 + α ( j) T T B( j) + (1 − α ( j) )Y T B( j) + α B( j) T

+(1 − α )B( j)Y + In

T T T ( j) Γ¯ 2 , P¯1 + Q1 A( j) − 2Q1 + α ( j) B( j) + (1 − α ( j) )Y T B( j) ,

are satisfied. Then, in the vertices j and i, the equilibrium is asymptotically stable as well as the disturbances are rejected in the entire polytope.

Proof: This is an extension of Theorem 10.5 for a polytopic uncertainties with some mathematical manipulations. Therefore, this theorem follows Theorem 10.5 proof.

Chapter 10. Sub-optimal control considering delays and parameter uncertainties

139

Remark 10.7 This robust control tuning method is conservative due to the definition of the matrix P, as well as, the attraction domain, X.

Corollary 10.8 Gain K, obtained from T and Q1 in Theorem 10.6, fulfills Theorem 10.5 and consequently guaranties both robust stability and robust disturbance rejection for a fixed delay.

The extension of this approach considering the saturation mechanism of the ENARC controller is open for future work.

10.4 Sub-optimal control result In this section, the previous sub-optimal control tuning is applied with the data reported in Section 7.2. Now, the clock frequency is taken fclk = 200MHz. This frequency introduces an one-sample-period delay (h1 = 1) in the control block output due to a power-performance trade-off. Likewise h2 = 2, thus h = 3. The uncertain parameters take the following ranges: • transistor characteristic, R0 , from 25Ω to 38Ω, • load dynamic resistance, RL , from 55.53Ω to 72.46Ω, • load capacitance, C, from 1pF to 1nF and • clock frequency, ωn , from 125MHz to 600MHz. Then, LMIs (10.6) are resolved, obtaining K¯ 1 = −7179 K¯ 2 = 12114. 

0.0004 This was obtained for c = 7.56, P1 = 0.0008

  0.0008 ·1011 and G = −35.09 1.931

(10.27) (10.28)  736.45 .

In this computation, any perturbation was not taken into account. However, this method can be applied for any L2 exogenous disturbance. Note that even if the control constant tuning is conservative, there is a feasible solution.

140

10.5. Simulation Results

10.5 Simulation Results In this section some simulations are performed in order to show the properties that the closed-loop system can achieve when the sub-optimal control gains are used in the ENARC controller. Likewise, a comparison between the performance achieved with respect to the control gains obtained by the previous control tuning given in Section 7.2 (K¯ 1 = −19.3 and K¯ 2 = 39.27), and the gains got by the sub-optimal control tuning (K¯ 1 = −7179 and K¯ 2 = 12114) is made. For these simulations, the parameters values given in Section 7.2 and the data given above are taken. We want to remark that in the following simulations, the delay is h = 3 and the load dynamic resistance, RL , will be time-varying, rL .

10.5.1 Uncertain clock frequency. In this kind of systems, clock frequency can be changes. The closed-loop system robustness (with K¯ 1 = −7179 and K¯ 2 = 12114) is displayed when the sampling frequency is ωn = 200MHz in Fig. 10.2 and ωn = 400MHz in Fig. 10.3. Observe that the effect of the delay is shown in system response. Note that the equilibrium is robust with respect to parameter uncertainties and delay.

10.5.2 Uncertain PMOS resistance In this first evaluation, it is assumed that the electrical characteristic of the PMOS can suffer changes. For this, 0.8R0% and 1.2R0 % is changed. This is shown in Fig. 10.4 and 10.5, respectively. Note that the system in the low voltage level converges to 0.677V and 0.797, which are inside the interval given by (10.3), Il = [0.675V, 0.799V ]. Likewise, the high voltage level converges to 1.133V , which is inside the interval given by (10.4), Ih = [1.132V, 1.155V ]. Therefore, system converges to the uncertain intervals. These tests show once again the system robustness.

10.5.3 Uncertain load parameter Finally, an example shows that system performance is sensitive to K¯ 1 and K¯ 2 .

Chapter 10. Sub-optimal control considering delays and parameter uncertainties

141

NTrans

a) 20 10 0 0

0.2

0.4

0.6

0.8

V(V)

t(s) b) 1.2 1 0.8 0

0.2

0.4

x 10

0.6

0.8

I(A)

t(s) d) 0.02 0 0

0.2

0.4

1 −6

1 −6

x 10

0.6

0.8

t(s)

1 −6

x 10

Figure 10.2: ωn = 200MHz and K¯ 1 = −7179 and K¯ 2 = 12114. Evolution of the: a) number of PMOS transistors switched on, b) vr (dashed) and vc (solid), c) current Il .

NTrans

a) 20 10 0 0

0.2

0.4

0.6

0.8

V(V)

t(s) b) 1.2 1 0.8 0

0.2

0.4

x 10

0.6

0.8

I(A)

t(s) d) 0.05 0 0

0.2

0.4

1 −6

x 10

0.6 t(s)

1 −6

0.8

1 −6

x 10

Figure 10.3: ωn = 400MHz and K¯ 1 = −7179 and K¯ 2 = 12114. Evolution of the: a) number of PMOS transistors switched on, b) vr (dashed) and vc (solid), c) current Il .

10.5. Simulation Results

142

NTrans

a) 20 10 0 0

0.2

0.4

0.6

0.8

V(V)

t(s) b) 1.2 1 0.8 0

0.2

0.4

x 10

0.6

0.8

I(A)

t(s) d) 0.02 0 0

0.2

0.4

1 −6

1 −6

x 10

0.6

0.8

t(s)

1 −6

x 10

Figure 10.4: 0.8R0 % and K¯ 1 = −7179 and K¯ 2 = 12114. Evolution of the: a) number of PMOS transistors switched on, b) vr (dashed) and vc (solid), c) current Il .

NTrans

a) 20 10 0 0

0.2

0.4

0.6

0.8

V(V)

t(s) b) 1.2 1 0.8 0

0.2

0.4

x 10

0.6

0.8

I(A)

t(s) d) 0.02 0 0

0.2

0.4

1 −6

x 10

0.6 t(s)

1 −6

0.8

1 −6

x 10

Figure 10.5: 1.2R0 % and K¯ 1 = −7179 and K¯ 2 = 12114. Evolution of the: a) number of PMOS transistors switched on, b) vr (dashed) and vc (solid), c) current Il .

Chapter 10. Sub-optimal control considering delays and parameter uncertainties

143

The capacitance employed in the previous simulations have been C = 1nF. In the following, it is desired to validate the system robustness when C = 1pF, i.e., 1000 times smaller. The lack of knowledge of the load in the real applications of these systems may achieve this change of three order of magnitude. Some simulations using both the original and suboptimal control tuning are made. In Fig. 10.6, it can be seen a wrong behavior of the controller when the previous control tuning given in Section 7.2 K¯ 1 = −19.3 and K¯ 2 = 39.27 are used. Nevertheless, Fig. 10.7 shows the simulation employing the sub-optimal control tuning, K¯ 1 = −19.3 and K¯ 2 = 39.27. Note that in Fig 10.6, the system does not respond to voltage variation. However, in Fig. 10.7 the system performance is satisfactory. This example shows the great robustness of the system when the sub-optimal control tuning is employed.

NTrans

a) 20 10 0 0

0.2

0.4

0.6

0.8

V(V)

t(s) b) 1.2 1 0.8 0

0.2

0.4

x 10

0.6

0.8

I(A)

t(s) d) 0.02 0 0

0.2

0.4

1 −6

1 −6

x 10

0.6 t(s)

0.8

1 −6

x 10

Figure 10.6: C = 1pF and K¯ 1 = −19.3 and K¯ 2 = 39.27. Evolution of the: a) number of PMOS transistors switched on, b) vr (dashed) and vc (solid), c) current Il .

10.6 Evaluation of the tuning methods In this section, we want to perform an evaluation of the two control tuning approaches presented in this part of the thesis.

10.6. Evaluation of the tuning methods

144

NTrans

a) 20 10 0 0

0.2

0.4

0.6

0.8

V(V)

t(s) b) 1.2 1 0.8 0

0.2

0.4

x 10

0.6

0.8

I(A)

t(s) d) 0.02 0 0

0.2

0.4

1 −6

x 10

0.6 t(s)

1 −6

0.8

1 −6

x 10

Figure 10.7: C = 1pF and K¯ 1 = −7179 and K¯ 2 = 12114. Evolution of the: a) number of PMOS transistors switched on, b) vr (dashed) and vc (solid), c) current Il .

The first tuning method for the ENARC controller was presented in Chapter 7. It was developed focussed on ensuring the equilibrium convergence around of a set point for the linearized closed-loop system. Later, in Chapter 9, an approximate stability analysis for the nonlinear system and the approximate continuous-time ENARC controller with this first tuning approach was performed. Consequently, this tuning mechanism takes into account the management of the current peaks and the fast transient periods in such a way that a high energy efficiency is achieved. However, the saturation of the total number of PMOS transistors, system delays and parameter uncertainties have not been studied . The second tuning method was developed in this chapter. It was developed taking into account the saturation of the total number of transistors, fast transient periods, delays, perturbation rejection and uncertain parameter as well as stability issues. Moreover, the closedloop system with this method also offers a high energy efficiency. For this, it was called sub-optimal control mechanism. Nevertheless, in this approach was not considered the management of the current peaks. An estimation of an attraction domain, that ensures a maximum variation of the switched transistors in every sampling time, could guarantee that these current peaks are small. As future work, it is desired to extend this result regarding the saturation mechanism that manages the current peaks. This last control tuning takes into account more control objectives, and thus, it consideres a closed-loop system closer of the real one. In order to test the achieved properties, some

Chapter 10. Sub-optimal control considering delays and parameter uncertainties

145

simulations are done employing the ENARC controller (i.e., the controller with the current peak management) with this sub-optimal tuning approach. It can be seen, the right system behavior.

10.7 Conclusions An energy-aware control has been developed for the Vdd-Hopping system. This controller achieves almost of the control requirements for SoCs technology. Nevertheless, parameter uncertainties and delays have been disregarded. Generally, this kind of systems have delays due to synchronization issues and performance constraints. Furthermore, parameters can change due to task requirements or can be varying on time. In this chapter, these important issues have been dealt with. An sub-optimal ‘conservative’ control tuning approach has been developed for the ENARC controller in order to achieve a robust closed-loop system with respect to the parameter uncertainties and delays. For this, the system is rewritten in a state-space form. The control has been based on H∞ theory applied to time-delay systems [26, 53, 101]. For this, some LMIs have been developed following Lyapunov-Krasovskii method. Conservativeness of the method for the Vdd-Hopping system has been discussed. This kind of method to tune a linear controller are employed in industrial applications [14, 129]. System robustness has been showed by means of some simulations. An evaluation of the two tunning methods presented in this part of the thesis have been performed taking into account the approximations employed in both approaches. A future research will be performed, in order to extend this result considering the current peak management in the control signal. Therefore, with this development, all control requirements specified in Chapter 6 have been achieved. A controller, pending patent under the name ENARC, was designed for the Vdd-Hopping system with the aim of reducing the dissipated energy. This controller has an energy-aware management of current-peaks in the set of PMOS transistors. In addition, a step reference is used, thus, it only needs to know the two set-points. This result comes from the possibility to control more than one transistor at once, i.e., to switch more than one transistor in a same sampling time. As a side effect the transient-periods are diminished. These improvements make that system is more energetically efficient. In a comparison performed with an ‘intuitive’ controller published in [99], it has been showed that energy-consumption is reduced a 96%. Furthermore, it presents a relative low number of computational blocks, what makes this controller feasible for industry applications. Finally, the closed-loop system has a robust equilibrium stability with respect to parameter uncertainties and delays.

10.7. Conclusions

146

In summary, in this work a controller for the Vdd-Hopping system has been obtained. This controller has the next properties: • high energy efficiency, • system stability, • small current peaks, • fast transient periods, • robustness with respect to parameter uncertainty, • robustness with respect to delays and • easy implementation. Hence, it achieves an interesting relevance for SoCs applications and, thus, in ARAVIS project implementation. An implementation of the ENARC controller in VHDL-AMS will be performed in the project context, in order to validate its performed.

Chapter 12 Conclusions and future work 12.1 Conclusions and contribution summary This thesis contributes to provide nonlinear control problem solutions to several classes of power converters. Specifically, this thesis deals with a DC-AC converter for applications of medium or high power, and a DC-DC converter for low-power applications. The first problem was developed in the ‘Departamento de Ingenier´ıa de Sistemas y Autom´atica’ at the ‘Universidad de Sevilla’ (Spain). And the second electronic application was raised in a project of the French government, sponsored by the international competitiveness pole Minalogic, called ARAVIS. The control problem was tackled in the ‘D´epartement d’Automatique de GIPSA-Lab’ at the ‘Institut Polytechnique de Grenoble’ and at the ‘Institut National de Recherche en Informatique et en Automatique de Grenoble’ (France). Both converters are regarded with respect to their work contexts, and thus, the associated control objectives are different. Due to the dissimilar natures of both applications mentioned before, the thesis is composed of two parts. The first part is focused on the control problem of the DC-AC converter. Its structure is based on a double DC-DC boost converter, obtaining a boost inverter. It is no-minimum phase 4th order nonlinear system. Its main objective is to achieve the desired voltage with a suitable control law, that does not require any reference signal. Furthermore, other problems, as no purely resistive and known loads, and the proposition of an estimated attraction region have been coped with. The second part of the thesis is devoted to a discrete DC-DC converter for low-voltage application in SoC. Its structure is the result of employing DC-DC converters in this kind of technology. It is based on the Vdd-Hopping technic in order to fulfill Dynamic Voltage Scaling, hence, it is called ‘Vdd-Hopping converter’. This system is a 1st order nonlinear system. Although, its simple model may not be attractive for control applications, it has relevant interest in this field due to the project context, where 147

148

12.1. Conclusions and contribution summary

it is dealt with. The control objectives come from the ARAVIS project. In this context, the relevant aim for this DC-DC converter is to achieve a high-efficiency, while the current peaks and the transient periods are minimized. This must be reached with a controller that ensures the convergence to the desired equilibrium points, global stability and robust behaviour with respect to delays and parameter uncertainties. Next up, the most important contributions of this thesis are highlighted. For the first part, i.e., controlling the boost inverter circuit, the most important contributions have been:

1. A control law for the boost inverter has been designed based on an energy shaping approach for oscillation generation. This method provides a relevant property, the system is autonomous, and hence, its analysis and implementation are easier, since the system needs no reference signal. Furthermore, the energy shaping approach employed to obtain the control structure ensures global stability. 2. A phase controller inspired by a phase-lock loop is presented in order to synchronize the two output voltages of both parts of the system. This idea is extended to synchronize the circuit with an external signal, as for example, with the electrical grid. 3. By means of developing a control adaptive for the unknown or/and slowing varying load connected to the boost inverter, the desired output voltage is always achieved. This control adaptive needs of a state observer for some variables, although all variables are measured. Global stability of the full system is proved by using singular perturbation method, for this, the system is rewritten in the suitable form by using time-scale separation. 4. The previous problems dealt with before have been extended to a load that is not purely resistive but also has an inductive component. These developed works, have been not considered in this thesis to make a simple reading, since they are just an extension. However, as the inverter control law as the adaptive controller developed for the boost inverter with an inductive load were published in [10] and in [12], respectively. 5. The last issue considered has been the estimation of an attraction region for the boost inverter. It provides a set of initial conditions corresponding to trajectories that converge towards the desired system behavior. This problem comes from the real nature of the boost inverter, which has several constraints, including saturations. This makes that the system has not global stability with the Lyapunov function obtained from energy shaping approach. Inspired by this estimation problem of an attraction region for the boost inverter, a general estimation method for this class of problem has been proposed. It takes advantage of certain defined Lyapunov function, that ensures global stability, raising a

Chapter 12. Conclusions and future work

149

simple optimization problem. This fact makes that the estimation is ‘conservative’. This problem is resolved rewriting it as a sum of squares optimization. This method is applied to the boost inverter. For this, the previous Lyapunov function, has been employed. The method provides a good ‘conservative’ estimation for the inverter.

The second part of the thesis is focused on the control problem of the DC-DC VddHopping converter. Here, the provided contributions have been:

1. A set of high-performance controllers has been proposed. These controllers have been developed by applying several control theories focused on reaching the desired equilibrium point by the closed-loop system. From the set of controllers, the one that provides a best performance is selected. In order to satisfy the objectives, some developments have been made to this controller, applying optimal and adaptive control theory. This control solution in spite of providing nice properties to the Vdd-Hopping converter behaviour, has a relevant drawback: its computational cost is very large. Therefore, it is not a suitable solution for the ARAVIS project. 2. An innovative controller for the DC-DC Vdd-Hopping converter has been developed based on the control structure with the smallest computational cost from the set of control solutions proposed before. This controller has been designed based on energyaware concept. Its originality is due to its current-peak managing through saturations with dynamic limits depending on the state of the system. It achieves the desired equilibrium points increasing the energy saving, reducing the current peaks and diminishing the transient periods. It covers almost all requirements in low-power technology. In addition, its simple structure is remarkable for industrial applications. 3. A global stability analysis of the nonlinear model of the Vdd-Hopping converter with the nonlinear controller presented before has been developed. The stability analysis of the closed-loop system has been involved due mainly to the saturation limits depending on the system state. The global analysis is ensured by LaSalle’s invariance principle. For simplicity, this analysis has been performed in continuous-time. It is assumed that the stability properties in discrete-time are conserved, as is common in control. 4. The last contribution is performed with respect to delays in parameter uncertainties. On the one hand, some delays can be presented due to the regarded work context of the Vdd-Hopping converter. They may be caused by synchronization issues, as well as, by providing an energy-performance trade-off. These delays can be considered as one only constant delay. On the other hand, depending on the specific application, the system parameters can be diverse and time-varying during the transient periods, as well. Consequently, the problem is to find the optimal gains for the controller. For this, an optimal tuning mechanism for these control gains based on H∞ theory is proposed. They are obtained resolving some Linear Matrix Inequalities (LMIs),

12.2. Future work

150

which have been developed by Lyapunov Krasovskii method. These LMIs ensure the equilibrium robustness with respect to delays and parameter uncertainties.

Finally, it can be concluded that nonlinear control theories can be powerful tools to provide solutions of several natures in industrial applications, in such a way that global stability of these systems is guaranteed.

12.2 Future work Following the investigations described in this thesis, the next future work will be taken up: In the boost inverter: 1. A physic implementation of the boost inverter will be made in order to test the performance of the control law proposed from applying energy shaping. In addition, this will allow to compare this control law with other controllers that have been already published for the same inverter. 2. An extension of the proposed adaptive control can be performed considering that not all the states are measured. 3. In the case of the global stability analysis to the system with the adaptive controller, to extend this analysis (which has been performed employing singular perturbation analysis) for a infinite time. 4. To find a less conservative solution to estimate the attraction region by employing another advanced tool to solve the sum of squares optimization problem. 5. . In the Vdd-Hopping converter: 1. An implementation of the controller proposed in this thesis will be performed in VHDL-AMS, in order to test its real approximate behaviour before to implement it in SoCs. 2. A better numerical solution to obtain an optimal voltage reference for the Lyapunov controller will be performed with another advanced mathematical tool. 3. Extension of the stability analysis of the closed-loop system in discrete-time.

Chapter 12. Conclusions and future work

151

4. Extension of the sub-optimal tuning approach for the control gains considering the saturation of the current peak management.

152

12.2. Future work

List of Publications In conferences: 1. C. Albea, C. Canudas-de-wit. and F. Gordillo, Control and Stability Analysis for the Vdd-hopping Mechanism, 3rd IEEE Multi-conference on Systems and Control (MSC’09), Saint Petersburg, Russia. 8-10 July 2009. 2. C. Albea, C. Canudas-de-wit. and F. Gordillo, Advanced Control Design for Voltage Scaling Converters, 34th IEEE Conference of the IEEE Industrial Electronics Society (IECON’08), Orlando, Florida, LA, USA.10-13 Nov. 2008 3. C. Albea, F. Gordillo and C. Canudas De Wit, Adaptive Control of the Boost Inverter with Load RL, 17th IFAC World Conference, Korea Republic. 6-11 July 2008. 4. C. Albea, C. Canudas-de-wit. and F. Gordillo, Diseo de Controladores para Convertidores con Escalado de Tensin, XXIX Jornadas de Automtica, Tarragona, Spain. 3-5 Sept. 2008. 5. C. Albea and F. Gordillo., Control of the boost DC-AC converter with RL load by energy shaping, 46th IEEE Conference on Decision and Control (CDC’07), New Orleans, LA, USA. 12-14 Dec. 2007. 6. M. Fiacchini, T. Alamo, C. Albea and E. Fernandez Camacho, ”Control Applications, Adaptive model predictive control of the hybrid dynamics of a fuel cell system, 1er IEEE Multi-conference on Systems and Control (MSC’07), Singapore. 1-3 Oct. 2007. 7. C. Albea, C. Canudas De Wit and F. Gordillo, Adaptive Control of the Boost DC-AC Converter, 1er IEEE Multi-conference on Systems and Control (MSC’07), Singapore. 1-3 Oct. 2007 8. C. Albea and F. Gordillo, C. Canudas De Wit, Control Adaptativo del Inversor Boost, XXVIII Jornadas de Automtica, 5-8 Huelva, Spain. 5-8 Sept. 2007. 9. C. Albea and F. Gordillo, Estimation of the Region of Attraction for a Boost DCAC Converter Control Law, 7th IFAC Symposium on Nonlinear Control Systems (Nolcos’07), Pretoria, South Africa. 22-24 August 2007. 153

154

12.2. Future work

10. C. Albea C, F. Gordillo and J. Aracil Santoja, Control of the boost DC-AC converter, 32th IEEE Conference of the IEEE Industrial Electronics Society (IECON’06), Paris, France. 6-10 Nov. 2006 11. C. Albea C., M. G. Ortega, F. Salas and F. R. Rubio, Aplicacion del control Hinf al PPCAR, XXVII Jornadas de Automtica, Almeria, Spain. 6-9 Sept. 2006. 12. C. Albea and F. Gordillo, Control del Convertidor Boost DC-AC por moldeo de Energa, XXVII Jornadas de Automtica, Almeria, Spain. 6-9 Sept. 2006

In journals: 1. C. Albea, F. Gordillo and C. Canudas De Wit., Adaptive Control Design for a Boost Inverter, in Engineering Control Technology, submitted 2009, accepted August 2010, published September 2010. 2. C. Albea and F. Gordillo, On the estimation of attraction domains for polynomial systems with constraints. Application to electronic converters” (submitted). 3. C. Albea, C. Canudas De Wit. and F. Gordillo, High Performance Control Design for Dynamic Voltage Scaling Devices, (under preparation). 4. C. Albea, F. Gordillo and C. Canudas De Wit. , Robust Control for Low-Power Converters (under preparation).

Patent: 1. C. Albea and C. Canudas de Wit, Dispositif de commande numrique pour un tableau de transistors PMOS en parallele Patent No: 08/07342. Filing date: 22 Dec. 2008 Publ. date: 22 January 2009.

Bibliography [1] D. Abramovitch. Phase-Locked Loops: a control centric tutorial. In Proc. IEEE American Control Conference (ACC), 1:1 –15, 2002. [2] Enrique Acha. Switching Power Supply Design. McGraw-Hill, 1997. [3] E. Achille, T. Martir´e, C. Glaize, and C. Joubert. Optimized DC-AC boost converters for modular photovoltaic grid-connected generators. 2, 2004. [4] R. Akhter. A new technique of PWM boost inverter for solar home application. BRAC University Journal, 4(1):39–45, 2007. [5] T. Alamo, A. Cepeda, D. Limon, and E.F. Camacho. Estimation of the domain of attraction for saturated discrete-time systems. International Journal of Systems Science, 37(8):575–583, 2006. [6] C. Albea and C. Canudas-de Wit. Dispositif de commande num´erique pour un tableau de transistors pmos en parallele. Patent Number:08/07342. [7] C. Albea, C. Canudas-de Wit, and F. Gordillo. Adaptive Control of the Boost DC-AC Converter. In Proc. IEEE Conference on Control Applications (CCA), pages 611–616, 2007. [8] C. Albea, C. Canudas-de Wit, and F. Gordillo. Advanced control design for voltage scaling converters. In Proc. IEEE Industrial Electronics (IECON), 2008. [9] C. Albea, C. Canudas-de Wit, and F. Gordillo. Control and stability analysis for the vdd-hopping mechanism. In Proc. IEEE Conference on Control Applications (CCA), 2009. [10] C. Albea and F. Gordillo. Control of the boost DC-AC converter with RL load by energy shaping. In Proc. IEEE Decision and Control (CDC), pages 2417–2422, 2007. [11] C. Albea, F. Gordillo, and J. Aracil. Control of the Boost DC-AC Converter by Energy Shaping. In Proc. IEEE Industrial Electronics (IECON), pages 754–759, 2006. [12] C. Albea, F. Gordillo, and C. Canudas-de Wit. Adaptive Control of the Boost Inverter with Load RL. In Proc. IFAC World Congress, pages 3316–3321, 2008. 155

156

BIBLIOGRAPHY

[13] J.M. Alonso, J. Ribas, J.J.D. Coz, A.J. Calleja, E.L. Corominas, and M. Rico-Secades. Development of a distributive control scheme for fluorescent lighting based on LonWorks technology. IEEE Trans. on Industrial Electronics, 47(6):1253–1262, 2000. [14] J. Anthonis, A. Seuret, J.P. Richard, and H. Ramon. Design of a pressure control system with dead band and time delay. IEEE Trans. on Control Systems Technology, 15(6):1103–1111, 2007. [15] T. Antonakopoulos, S. Pressas, and V. Makios. Single chip controller for a highswitching frequency DC/AC power inverter. International Journal of Electronics, 73(1):229–240, 1992. [16] J. Aracil and F. Gordillo. On the control of oscillations in DC-AC converters. In Proc. IEEE Industrial Electronics (IECON), 4, 2002. [17] K.J. Astrom, T. H¨agglund, C.C. Hang, and W.K. Ho. Automatic tuning and adaptation for PID controllers-a survey. Control Engineering Practice, 1(4):699–714, 1993. [18] I. Barbi and S.A.O. da Silva. Sinusoidal line current rectification at unity power factor with boost quasiresonant converters. pages 553–562, 1990. [19] D. Biel, E. Fossas, F. Guinjoan, E. Alarc´on, and A. Poveda. Application of slidingmode control to the design of a buck-basedsinusoidal generator. IEEE Trans. on Industrial Electronics, 48(3):563–571, 2001. [20] D. Biel, F. Guinjoan, E. Fossas, and J. Chavarr´ıa. Sliding-mode control design of a boost-buck switching converter for AC signal generation. IEEE Trans. on Circuits and Systems I: Regular Papers, 51(8):1539–1551, 2004. [21] F. Blanchini and S. Miani. Any domain of attraction for a linear constrained system is a tracking domain of attraction. In Proc. IEEE in the Conference on Decision and Control (CDC), 5:5065–5070, 2000. [22] S.R. Bowes and D. Holliday. Optimal regular-sampled PWM inverter control techniques. IEEE Trans. on Industrial Electronics, 54(3):1547–1559, 2007. [23] T.D. Burd and R.W. Brodersen. Design issues for dynamic voltage scaling. In Proc. IEEE International Symposium on Low Power Electronics and Design (ISLPED), pages 9–14, 2000. [24] R. C´aceres and I. Barbi. Sliding mode controller for the boost inverter. pages 247–252, 1996. [25] R C´aceres and I Barbi. A boost DC-AC converter: analysis, design, and experimentation. IEEE Trans. on Power Electronics, 14(1):134–141, 1999. [26] Y.Y. Cao, Z. Lin, and T. Hu. Stability analysis of linear time-delay systems subject to input saturation. IEEE Trans. on Circuits and Systems I: Fundamental Theory and Applications, 49(2):233–240, 2002.

BIBLIOGRAPHY

157

[27] Z. Cao, B. Foo, L. He, and M. Van der Schaar. Optimality and improvement of dynamic voltage scaling algorithms for multimedia applications. pages 681 –690, 2010. [28] J.M. Carrasco, J.M. Quero, F.P. Ridao, M.A. Perales, and L.G. Franquelo. Sliding mode control of a DC/DC PWM converter with PFC implementedby neural networks. IEEE Trans. on Circuits and Systems I: Fundamental Theory and Applications, 44(8):743–749, 1997. [29] L. Castaner and S. Silvestre. Modelling photovoltaic systems using PSpice. Wiley, 2003. ˚ en. Feedback–feedforward schedul[30] A. Cervin, J. Eker, B. Bernhardsson, and K.E. Arz´ ing of control tasks. Real Time Systems, 23(1):25–53, 2002. [31] G. Chesi. Estimating the domain of attraction for non-polynomial systems via LMI optimizations. Automatica, 45(6):1536–1541, 2009. [32] G. Chesi, A. Garulli, A. Tesi, and A. Vicino. Solving quadratic distance problems: an LMI-based approach. IEEE Trans. on Automatic Control, 48(2):200–212, 2003. [33] D. Christiansen. Electronics engineers’ handbook. 1996. [34] A. Das, S. Ozdemir, G. Memik, and A. Choudhary. Evaluating voltage islands in CMPs under process variations. 2007. [35] K.C.A. de Souza, M.R. de Castro, and F. Antunes. A DC/AC converter for singlephase grid-connected photovoltaic systems. In Proc. IEEE Industrial Electronics (IECON), 4:3268 – 3273, 2000. [36] H. Deng, X. Duan, N. Sun, Y. Ma, A.Q. Huang, and D. Chen. Monolithically Integrated Boost Converter Based on 0.5-mm CMOS Process. IEEE Trans. on Power electronics, 20(3):628–638, 2005. [37] J.W. Dixon and B.T. Ooi. Indirect current control of a unity power factor sinusoidal currentboost type three-phase rectifier. IEEE Trans. on industrial electronics, 35(4):508–515, 1988. [38] R. Dobkin, R. Ginosar, and C.P. Sotiriou. Data synchronization issues in GALS SoCs. pages 170–179, 2004. [39] S. Durand and N. Marchand. Fast Predictive Control of Micro Controller’s EnergyPerformance Tradeoff. In Proc. IEEE Conference on Control Applications (CCA), 2009. [40] C. Elmas, O. Deperlioglu, and H.H. Sayan. Adaptive fuzzy logic controller for DCDC converters. Expert Systems with Applications, 36(2):1540–1548, 2009.

158

BIBLIOGRAPHY

[41] J.H.R. Enslin, M.S. Wolf, D.B. Snyman, and W. Swiegers. Integrated photovoltaic maximum power point tracking converter. IEEE Trans. on industrial electronics, 44(6):769–773, 1997. [42] R. Erickson and D. Maksimovic. High efficiency DC-DC converters for batteryoperated systems with energy management. Worldwide Wireless Communications, 1995. [43] R.W. Erickson and D. Maksimovi´c. Fundamentals of power electronics. Kluwer Academic Pub, 2001. [44] G. Escobar, A.J. Van der Schaft, and R. Ortega. A Hamiltonian viewpoint in the modeling of switching power converters. Automatica-Oxford-, 35:445–452, 1999. [45] L. Fesquet and H. Zakaria. Controlling Energy and Process Variability in System-onChips: needs for control theory. In Proc. IEEE Conference on Control Applications (CCA), 27(1):302–308, 2009. [46] E. Fossas, J.M. Olm, and H. Sira-Ram´ırez. Iterative approximation of limit cycles for a class of Abel equations. Physica D: Nonlinear Phenomena, 237(23):3159–3164, 2008. [47] E. Fossas, J.M. Olm, and H. Sira-Ram´ırez. Iterative approximation of limit cycles for a class of Abel equations. Physica D: Nonlinear Phenomena, 237(23):3159–3164, 2008. [48] E. Fossas-Colet and J.M. Olm-Miras. Asymptotic tracking in dc-to-dc nonlinear power converters. Discrete and Continuous Dynamical Systems Series B, 2(2):295– 311, 2002. [49] E. Fridman, A. Pila, and U. Shaked. Regional stabilization and Hinf control of timedelay systems with saturating actuators. International Journal of Robust and Nonlinear Control, 13(9):885–907, 2003. [50] E. Fridman, A. Seuret, and J.P. Richard. Robust sampled-data stabilization of linear systems: an input delay approach. Automatica, 40(8):1441–1446, 2004. [51] E. Fridman and U. Shaked. A Descriptor System Approach to H∞ Control of Linear Time-Delay Systems. IEEE Trans. on Automatic Control, 47(2):253–270, 2002. [52] E. Fridman and U. Shaked. An lmi approach to stability of discrete delay systems. In Proc. of the European Control Conference, 2003. [53] E. Fridman and U. Shaked. Delay-dependent Hinf control of uncertain discrete delay systems. European Journal of Control, 11(1):29–39, 2005. [54] R. Genesio, M. Tartaglia, and A. Vicino. On the estimation of asymptotic stability regions: State of the art and new proposals. in Trans. on IEEE Automatic Control, 30(8):747–755, 1985.

BIBLIOGRAPHY

159

[55] M. Gildersleeve, H.P. Forghani-zadeh, and G.A. Rinc´on-Mora. A comprehensive power analysis and a highly efficient, mode-hopping DC-DC converter. pages 153– 156. [56] F. Gomez-Estern, J. Aracil, and F. Gordillo. The Hopf bifurcation and controlled oscillations in electromechanical systems. In Proc. Mediterranean Conference on Control and Automation (MED), 2002. [57] F. Gordillo, J. Aracil, and F. Gomez-Estern. Stabilization of autonomous oscillations and the Hopf bifurcation in the ball and beam. In Proc. IEEE Decision and Control (CDC), 4:3924 – 3925, 2002. [58] F. Gordillo, DJ Pagano, and J. Aracil. Autonomous oscillation generation in electronic converters. In Proc. IEEE International Workshop on Electronics and System Analysis (IWESA), pages 57–65, 2004. [59] K. Graubard and D.K. Hartline. Full-wave rectification from a mixed electricalchemical synapse. Science, 237(4814):535, 1987. [60] R. Griø´ , E. Fossas, and D. Biel. Sliding mode control of a full-bridge unity power factor rectifier. Lecture Notes in Control and Information Sciences, pages 139–148, 2002. [61] K. Gu, V. Kharitonov, and J. Chen. Stability of time-delay systems. Birkhauser, 2003. [62] K. Gu and S.I. Niculescu. Stability Analysis of Time-Delay Systems: A Lyapunov Approach. Advanced Topics in Control Systems Theory, Lecture Notes from FAP, pages 139–170, 2005. [63] Sira-Ram´ırez H., R. Tarantino-Alvarado, and O. Llanes-Santiago. Adaptive feedback stabilization in PWM-controlled DC-to-DC power supplies. International Journal of Control, 57(3):599–625, 1993. [64] A. Hadri-Hamida, A. Allag, M.Y. Hammoudi, S.M. Mimoune, S. Zerouali, M.Y. Ayad, M. Becherif, E. Miliani, and A. Miraoui. A nonlinear adaptive backstepping approach applied to a three phase PWM AC-DC converter feeding induction heating. Communications in Nonlinear Science and Numerical Simulation, 14(4):1515–1525, 2009. [65] J. Hellerstein and S. Parekh. Feedback control of computing systems. Wiley-IEEE Press, 2004. [66] G.C. Hsieh and J.C. Hung. Phase-Locked Loop techniques. A survey. IEEE Trans. on Industrial Electronics, 43(6):609–615, 1996. [67] T. Hu and Z. Lin. Control systems with actuator saturation: analysis and design. Birkhauser, 2001.

160

BIBLIOGRAPHY

[68] B. Huber and B. Sturmfels. A polyhedral method for solving sparse polynomial systems. in Trans Mathematics of Computation, 64(212):1541–1555, 1995. [69] E.C. Ifeachor and B.W. Jervis. Digital Signal Processing: A Practical Approach. Pearson Education, 2002. [70] K.S. Jeong, W.Y. Lee, and C.S. Kim. Energy management strategies of a fuel cell/battery hybrid system using fuzzy logics. Journal of Power Sources, 145(2):319– 326, 2005. [71] M. Johansson and A. Rantzer. Computation of piecewise quadratic Lyapunov functions for hybridsystems. IEEE Trans. on automatic control, 43(4):555–559, 1998. [72] J.G. Kassakian and M.F. Schlecht. High-frequency high-density converters for distributed power supplysystems. Proceedings of the IEEE, 76(4):362–376, 1988. [73] M.S. Kaster and D.J. Pagano. Control of Autonomous Oscillations in Buck-Based Inverters. pages 899–903, 2005. [74] N. Kavaldjiev and G.J.M. Smit. A survey of efficient on-chip communications for SoC. 2003. [75] H. Kawaguchi, G. Zhang, S. Lee, and T. Sakurai. An LSI for VDD-hopping and MPEG4 system based on the chip. pages 918–921, 2001. [76] H.K. Khalil. Nonlinear Systems. Prentice Hall, third edition edition, 2002. [77] M.H. Klein, J.P. Lehoczky, and R. Rajkumar. Rate-monotonic analysis for real-time industrial computing. IEEE Computer, 27(1):24–33, 1994. [78] P.V. Kokotovi´c, H.K. Khalil, and J. O’reilly. Singular Perturbation Methods in control Analysis and Design. SIAM, 1999. [79] A. Kusko and M.T. Thompson. Power quality in electrical systems. McGraw-Hill Professional, 2007. [80] H. Kwakernaak and R. Sivan. Linear optimal control systems. Wiley-Interscience New York, 1972. [81] R.G. Langlois, D.M. Hanna, and R.J. Anderson. Implementing preview control on an off-road vehicle with active suspension. page 340, 1992. [82] S. Lee and T. Sakurai. Run-time voltage hopping for low-power real-time systems. In Proc. IEEE conference on Design automation (CDA), pages 806–809, 2000. [83] S. Lee and T. Sakurai. Run-time voltage hopping for low-power real-time systems. pages 806–809, 2000. [84] W. Leonhard. Control of electrical drives. Springer-Verlag, 2001.

BIBLIOGRAPHY

161

[85] A. Levin. An analytical method for estimation the domain of attraction for polynomial differential equations. in Trans. on IEEE Automatic Control, pages 2471–2475, 1994. [86] F.L. Lewis and V.L. Syrmos. Optimal control. Wiley-Interscience, 1995. [87] G.P. Liu and S. Daley. Optimal-tuning PID control for industrial systems. Control Engineering Practice, 9(11):1185–1194, 2001. [88] P. Liu, J. Liu, and L. Geng. A dynamic buck converter with ultra fast response and low voltage ripples designed for DVS systems. IEICE Electronics Express, 6(21):1490– 1496, 2009. [89] D. Ma. Automatic Substrate Switching Circuit for On-Chip Adaptive Power-Supply System. IEEE Trans.on Circuits and Systems II: Express Briefs, 54(7):641–645, 2007. [90] H. Ma and S. Han. Buck inverter with sliding mode control strategy. 1, 2004. [91] M. Malinowski, G. Marqu´es, M. Cichowlas, and M.P. Kazmierkowski. New direct power control of three-phase PWM boost rectifiers under distorted and imbalanced line voltage conditions. 1, 2003. [92] D. Marculescu and A. Iyer. Power and performance evaluation of globally asynchronous locally synchronous processors. In Proc. International Symposium on Computer Architecture (ISCA), 2002, 30:2. [93] F. Martorell, S.D. Cotofana, and A. Rubio. An analysis of internal parameter variations effects on nanoscaled gates. IEEE Trans. on Nanotechnology, 7(1):24, 2008. [94] A.C.H. MeVay and R. Sarpeshkar. Predictive comparators with adaptive control. IEEE Trans. on Circuits and Systems II: Analog and Digital Signal Processing, 50(9):579– 588, 2003. [95] C. Meza, D. Biel, J. Mart´ınez, and F. Guinjoan. Boost-buck inverter variable structure control for grid-connected photovoltaic systems. pages 1318–1321, 2005. [96] R.D. Middlebrook. Power Electronics: Topologies, Modelling, and Measurement. Advances in switched-mode power conversion, page 17, 1983. ´ [97] R.D. Middlebrook and S. Cuk. A general unified approach to modelling switchingconverter power stages. International Journal of Electronics, 42(6):521–550, 1977. [98] S. Miermont. Contrle distribu de la tension d’alimentation dans les architectures GALS et proposition d’un slecteur dynamique d’alimentation. PhD thesis, Institute National Polytechnique de Grenoble, December 2008. [99] S. Miermont, P. Vivet, and M. Renaudin. A Power Supply Selector for Energy-and Area-Efficient Local Dynamic Voltage Scaling. Lectures Notes in Computer Science, 4644:556, 2007.

162

BIBLIOGRAPHY

[100] N. Mohan, T.M. Undeland, and W.P. Robbins. Power electronics. Wiley New York, 1995. [101] Y.S. Moon, P. Park, W.H. Kwon, and Y.S. Lee. Delay-dependent robust stabilization of uncertain state-delayed systems. International Journal of Control, 74(14):1447– 1455, 2001. [102] T.C.T. Ng, F.H.F. Leung, and P.K.S. Tam. A simple adaptive control strategy for regulated switching DC-DC converter based on grid-point concept. In Proc. IEEE Power Electronics Specialists Conference (PESC), 2:1494 –1498, 1996. [103] K.J. Nowka, G.D. Carpenter, E.W. MacDonald, H.C. Ngo, B.C. Brock, K.I. Ishii, T.Y. Nguyen, and J.L. Burns. A 32-bit PowerPC system-on-a-chip with support for dynamic voltage scaling and dynamic frequency scaling. IEEE Journal of Solid-State Circuits, 37(11):1441–1447, 2002. [104] K. Ogata. Discrete-time control systems. Prentice-Hall Englewood Cliffs, NJ, 1987. [105] K. Ogata. Modern control engineering. Prentice Hall, 2009. [106] K. Onizuka and T. Sakurai. VDD-Hopping Accelerator for On-Chip Power Supplies Achieving Nano-Second Order Transient Time. pages 145–148, 2005. [107] M. Osama and T.A. Lipo. A new inverter control scheme for induction motor drives requiringwide speed range. IEEE Trans. on Industry Applications, 32(4):938–944, 1996. [108] D. J. Pagano, J. Aracil, and F. Gordillo. Autonomous oscillation generation in the boost converter. In Proc. IFAC World Congress, 2005. [109] A. Papachristodoulou and S. Prajna. On the construction of Lyapunov functions using the sum of squares decomposition. In Proc. IEEE Decision and Control (CDC), 3:3482–3487, 2002. [110] T.E. Pare, H. Hindi, J.P. How, and S.P. Boyd. Synthesizing stability regions for systems with saturating actuators. In Proc. IEEE Conference Decision and Control (CDC), 2:1981–1982, 1998. [111] P.A. Parrilo. Structured Semidefinite Programs and Semialgebraic Geometry Methods in Robustness and Optimization. PhD thesis, California Institute of Technology, Pasadena, California, 2000. [112] A. Pavlov and M. Sachdev. CMOS SRAM circuit design and parametric test in nanoscaled technologies: process-aware SRAM design and test. Springer-Verlag, 2008. [113] E. Pinheiro, R. Bianchini, E.V. Carrera, and T. Heath. Dynamic cluster reconfiguration for power and performance. Compilers and operating systems for low power, 2001.

BIBLIOGRAPHY

163

[114] C. Pittet, S. Tarbouriech, and C. Burgat. Stability regions for linear systems with saturating controls viacircle and Popov criteria. 5:4518–4523, 1997. [115] L. Pontryagin, A. Andronov, and A. Vitt. On the statistical treatment of dynamical systems. Noise in nonlinear dynamical systems, 1:329, 1989. [116] S. Prajna, A. Papachristodoulou, and P.A. Parrilo. Introducing SOSTOOLS: a general purpose sum of squares programming solver. In Proc. IEEE Conference on Decision and Control (CDC), 1:741–746, 2001. [117] S. Prajna, A. Papachristodoulou, P. Seiler, and P.A. Parrilo. SOSTOOLS, Sum of Squares Optimization Toolbox for MATLAB. available from http://www.cds.caltech.edu/sostools, 2004. [118] A.I. Pressman, K. Billings, and T. Morey. Switching power supply design. McGrawHill Professional, 2009. [119] S.K. Reynolds. A DC-DC converter for short-channel CMOS technologies. IEEE Trans. on Solid-State Circuits, 32(1):111–113, 1997. [120] L. Reznik, O. Ghanayem, and A. Bourmistrov. PID plus fuzzy controller structures as a design base for industrial applications. Engineering Applications of Artificial Intelligence, 13(4):419–430, 2000. [121] M. Rico-Secades, E.L. Corominas, J.M. Alonso, J. Ribas, J. Cardes´ın, A.J. Calleja, and J. Garc´ıa-Garc´ıa. Complete low-cost two-stage electronic ballast for 70-W highpressure sodium vapor lamp based on current-mode-controlled buck-boost inverter. IEEE Transac. on Industry applications, 41(3):728–734, 2005. [122] D. R´ıos, A. Buhrig, and M. Renaudin. Power Consumption Reduction using dynamic control of Microprocessor performance. Integrated Circuit and System Design. International Workshop on Power and Timing Modeling, Optimization and Simulation (PATMOS), pages 10–18, 2005. [123] M. R´ıos-Bol´ıvar, A.S.I. Zinober, and H. Sira-Ram´ırez. Dynamical Sliding Mode Control via Adaptive Input-Output Linearization: A Backstepping Approach. Robust control via variable structure and Lyapunov techniques, page 15, 1996. [124] B.G. Romanchuk. Computing regions of attraction with polytopes: planar case. Automatica, 32(12):1727–1732, 1996. [125] T. Saito, H. Torikai, and W. Schwarz. Switched dynamical systems with double periodic inputs: An analysis tool and its application to the buck-boost converter. IEEE Trans. on Circuits and Systems I: Fundamental Theory and Applications, 47(7):1038– 1046, 2000. [126] P. Sanch´ıs, A. Ursaea, E. Gubia, and L. Marroyo. Boost DC–AC Inverter: A New Control Strategy. IEEE Trans. on Power Electronics, 20(2):343–353, 2005.

164

BIBLIOGRAPHY

[127] P. Sanch´ıs, A. Ursua, E. Gubia, and L. Marroyo. Design and experimental operation of a control strategy for the buck-boost DC-AC Inverter. In Proc. IEEE Electric Power Applications, 152(3):660–668, 2005. [128] M. Seno, Nakai, S. Akuii, K. Seno, T. Meguro, T. Seki, Hashiguchi A. Kondo, T, H. Kumano, and M. Shimura. Dynamic Voltage and Frequency Management for a Low-Power Embedded Microprocessor. IEEE Journal of Solid-State Circuits, 40(1), 2005. [129] A. Seuret, F. Michaut, J.P. Richard, and T. Divoux. Networked control using GPS synchronization. page 6, 2006. [130] L.F. Shampine, J. Kierzenka, and M.W. Reichelt. Solving boundary value problems for ordinary differential equations in Matlab with bvp4c. Manuscript available at http://200.13.98.241/ martin/irq/tareas1, 2000. [131] Y.B. Shtessel, A.S.I. Zinober, and I.A. Shkolnikov. Sliding mode control of boost and buck-boost power converters using method of stable system centre. Automatica, 39(6):1061–1067, 2003. [132] O. Silven, T. Rintaluoma, and K. Jyrkka. Implementing energy efficient embedded multimedia. Future, 2000(3000):700–1000, 2010. [133] H. Sira-Ram´ırez. DC-to-AC power conversion on a’boost’converter. International Journal of Robust and Nonlinear Control, 11(6):589–600, 2001. [134] H. Sira-Ram´ırez and M. R´ıos-Bol´ıvar. Sliding mode control of DC-to-DC power converters via extended linearization. IEEE Trans. on Circuits and Systems I: Fundamental Theory and Applications, 41(10):652–661, 1994. [135] H. Sira-Ram´ırez, M. R´ıos-Bol´ıvar, and A.S.I. Zinober. Adaptive Dynamical InputOutput Linearization of DC to DC Power Converters: A Backstepping Approach. Robust and Nonlinear Control, 7(3):279–296, 1997. [136] D.C. Snowdon, S. Ruocco, and G. Heiser. Power management and dynamic voltage scaling: Myths and facts. 2005. [137] J.A. Solomon and G. Sperling. Full-wave and half-wave rectification in second-order motion perception. Vision Research, 34(17):2239–2257, 1994. [138] S. Song, J.K. Kim, C.H. Yim, and H.C. Kim. H∞ control of discrete-time linear systems with time-varying delays in state. Automatica, 35:1587–1591, 1999. [139] S. Srinivasan and G. Venkataramanan. Comparative evaluation of PWM AC-AC converters. 1, 1995. [140] A. Srivastava, S. Shah, K. Agarwal, D. Sylvester, D. Blaauw, and S. Director. Accurate and efficient gate-level parametric yield estimation considering power/performance correlation. In Proc. of Austin Center for Advanced Studies, 2005.

BIBLIOGRAPHY

165

[141] A.J. Stratakos, R.W. Brodersen, and S.R. Sanders. High-efficiency low-voltage DCDC conversion for portable applications. pages 105–110, 1994. [142] A.J. Stratakos, S.R. Sanders, and R.W. Brodersen. A low-voltage CMOS DC-DC converter for a portable battery-operated system. pages 619–619, 1994. [143] A.J. Stratakos, S.R. Sanders, and R.W. Brodersen. A low-voltage CMOS DC-DC converter for a portable battery-operatedsystem. In Proc. IEEE Power Electronics Specialists Conference (PESC), pages 619–626, 1994. [144] J.F. Sturm. Using SeDuMi 1.02, a MATLAB toolbox for optimization over symmetric cones. Optimization Methods and Software, 11–12:625–653, 1999. [145] S.C. Tan, Y.M. Lai, and C.K. Tse. Adaptive control schemes for stabilizing switching frequency of sliding mode controlled power converters. In Proc. IEEE Power Electronics and Applications, 2005. [146] B. Tibken and K.F. Dilaver. Computation of subsets of the domain of attraction for polynomial systems. In Proc. of the IEEEConference Decision and Control (CDC), 3:2651–2656, 2002. ´ [147] N. V´azquez, J. Alvarez, C. Aguilar, and J. Arau. Some critical aspects in sliding mode control design for the boostinverter. pages 76–81, 1998. ´ [148] N. V´azquez, J. Alvarez, and J. Arau. Passivity based control for the boost inverter: details of implementation. 2002. ´ [149] N. V´azquez, D. Cort´es, C. Hern´andez, J. Alvarez, and J. Arau. A new nonlinear control strategy for the boost inverter. 3, 2003. ´ [150] N. V´azquez, D. Cort´es, C. Hernandez, J. Alvarez, J. Arau, and J. Alvarez. A new nonlinear control strategy for the boost inverter. In Proc. IEEE Power Electronics Specialist Conference (PESC), 3:1403–1407, 2003. [151] J. Verschelde. Algorithm 795: Phcpack: a general-purpose solver for polynomial systems by homotopy continuation. in ACM Trans. on Mathematical Software (TOMS), 25(2):251–276, 1999. [152] H.O. Wang, K. Tanaka, and M.F. Griffin. An approach to fuzzy control of nonlinear systems: Stability and design issues. IEEE Trans. on Fuzzy Systems, 4(1):14–23, 1996. [153] G.Y. Wei and M. Horowitz. A fully digital, energy-efficient, adaptive power-supply regulator. IEEE Trans. on Solid-State Circuits, 34(4):520–528, 1999. [154] M. Yap-San-Min. Analyse statistique de l’impact des variations locales sur les courses de signaux dans une m´emoire sram embarqu´ee. Master’s thesis, Universit´e Montpellier II, 2008.

166

BIBLIOGRAPHY

[155] K. Zhou, J.C. Doyle, and K. Glover. Robust and optimal control. Prentice Hall Englewood Cliffs, NJ, 1996. [156] L. Zhu, K. Wang, F.C. Lee, and J.S. Lai. New start-up schemes for isolated full-bridge boost converters. IEEE Transactions on Power Electronics, 18(4):946–951, 2003. [157] Y. Zhu and F. Mueller. Feedback EDF scheduling exploiting dynamic voltage scaling. In Proc. IEEE Real-Time and Embedded Technology and Applications Symposium (RTAS), pages 84–93, 2004.