Power Stage Design of Fourth-Order DC-DC Converters by Means of Principal Components Analysis
Power Stage Design of Fourth-Order DC-DC Converters by Means of Principal Components Analysis
Power Stage Design of Fourth-Order DC-DC Converters by Means of Principal Components Analysis
Fig. 2. Voltage mode control scheme. Fig. 3. Power losses budget distribution.
where and represent the phase margin and gain margin if the system is asymptotically stable, the following matrices can
respectively, as in Fig. 5. The third solution seems to be the be defined:
suitable in terms of dynamic performances but, as can be de-
duced from Fig. 5(b), there is a dangerous approaching of the
loop gain to the 0 dB line which can cause multiple crossover
(2)
points. Since the frequency response is susceptible of parame-
ters values variations due to tolerances and drifts and also de-
pends on the varying working conditions through the duty cycle
value, the presence of multiple crossover points can lead to the (3)
system instability.
The goal of this paper is to show how to select an optimal as the Controllability and Observabilty Gramian matrices
combination of passive components in order to avoid the which respectively generate the subspaces of controllable and
problem of double crossover and to make the voltage mode observable states [29], [30].
control design easier and effective. Any Gramian matrix defined as
(4)
III. PCA IN SYSTEM ANALYSIS
has the following properties:
The PCA [18] is a mathematical theory inherited from statis- 1) is symmetric and positive semidefinite;
tics and used to obtain low-dimensional approximate descrip- 2) the eigenvelues are real, nonnegative and sorted in de-
tions for multidimensional systems. PCA is mathematically de- creasing order ;
fined as an orthogonal linear transformation that transforms the 3) the eigenvectors are orthonormal ;
data to a new coordinate system such that the greatest variance 4) the eigenvalues are given by
by any projection of the data comes to lie on the first coordi-
nate (called the first principal component), the second greatest
variance on the second coordinate, and so on. It involves the
where is the -row of the matrix .
computation of the eigenvalue decomposition or Singular value
Since the Gramian is symmetric it is also diagonalizable
decomposition of a data set. Such a theory can be adapted to the
Linear Systems (LS) [18] in order to find a new representation (5)
useful for the order reduction. Let it be considered a -order
system Input-State-Output (I-S-O) model: and
(1) (6)
If, for example, the new coordinate system based on the the system shows the same behavior of a second-order converter,
transformation matrix , which diagonalizes the controlla- thus eliminating all the undesired effects of a double resonance
bility gramian : peak. Let it be considered the model of a SEPIC converter ob-
(8) tained through the Averaging State Space technique
(9) (14)
where
where
(10)
verter, one could understand for which values of the parameters ters considering the losses of all passive components as well as the losses of
switching devices according to [32].
3The existence of the internally balanced representation is demonstrated. The 5Damping factor was obtained directly extracting and inverting the value of
transformation matrix P can be obtained through different algorithms, see [31]. the Q factor from the control-to-output transfer function evaluated in Matlab.
DE NARDO et al.: POWER STAGE DESIGN OF FOURTH-ORDER DC–DC CONVERTERS BY MEANS OF PCA 2871
Fig. 7. Optimal solutions set: (a) V and (b) V . Sidebar gives the efficiency value.
TABLE I
DATABASE OF THE COMMERCIALLY AVAILABLE PARAMETERS solutions at , together with all the 4096 possible solu-
tions have been plotted in Fig. 6. In order to have a simpler inter-
pretation, the solution sets shown in Fig. 6 in the three-dimen-
sional space have been represented by using two bidimensional
plots, shown in Fig. 7, wherein the “Efficiency” dimension was
compressed and shown in the gray-scale palette.
The Pareto front is [17] the set of design solutions that are
nondominated, in the sense that no other solutions exhibit an
improvement of at least one performance figure without wors-
ening the others.
In Fig. 7 the ripple-based design solution introduced in
Section II is compared with the Pareto-optimal ones at both
The solutions can be laid out in the three-dimensional the input voltage values. The Pareto front obtained at
space of the performances or into the four-dimensional space consists of four solutions, while eight solutions give the Pareto
of parameters. front at . As can be deduced from Fig. 7 and Table II,
The design procedure has been conducted at both the extreme solutions in the upper left sides of the fronts are suitable for
values of the input voltage fixed at the beginning of Section II small Q factor, and therefore maximum damping factor, while
for the SEPIC case study, namely, 2.7 and 4.2 V, in the sequel solutions in the lower right sides show the best dynamic be-
referred to as and . havior because their higher Th factor make their corresponding
After the exhaustive search, the two Pareto-optimal fronts Bode diagram of the open loop transfer function closer to
[17], composed of a total of four solutions at and eight that one typical of a second-order circuit. It is worth noting
2872 IEEE TRANSACTIONS ON POWER ELECTRONICS, VOL. 23, NO. 6, NOVEMBER 2008
Fig. 9. (a) PCA Threshold (sidebar) vs. ripples ratios @ V , (b) damping Factor (sidebar) vs. ripples ratios @ V .
TABLE IV
(a) RIPPLES VALUES CORRESPONDING TO THE SOLUTIONS BELONGING TO THE
PARETO FRONT @ V . (b) RIPPLES VALUES CORRESPONDING TO THE
SOLUTIONS BELONGING TO THE PARETO FRONT @ V
are better than the one obtained by the ripple based design and
this is obtained with lower values of inductances (Table III). Fig. 10. (a) Open-loop frequency response of the four selected solutions and
An important issue is that in order to get low values of and (b) magnification to highlight the differences.
high values of (lower order behavior) the values of C1 and
C2 tend to reach the minimum and the maximum capacitance
respectively, meanwhile the inductance values L1 and L2 tend , it can be put in evidence that equalization of the values
to be equalized. Looking at the pictures in Figs. 8, 9, referred to of the current ripples leads the system to
2874 IEEE TRANSACTIONS ON POWER ELECTRONICS, VOL. 23, NO. 6, NOVEMBER 2008
Fig. 11. (a) Open-loop frequency response of four (#1,#4,#6,#8) selected solutions @ V and (b) magnification to highlight the differences.
Fig. 12. Bode diagrams of the best solutions in terms of Th value: (a) @ V , (b) V .
Fig. 13. Bode diagrams corresponding to the solution #3: (a) @ V , (b) @ V .
approach better Th values, while an increased value of the ratio The previous analysis helps the designer in understanding the
gives better values of the damping factor. behavioral dependence from the ripples values as well as from
the components values: this is usually not possible by a simple
DE NARDO et al.: POWER STAGE DESIGN OF FOURTH-ORDER DC–DC CONVERTERS BY MEANS OF PCA 2875
Fig. 14. Bode diagrams corresponding to the solution #3 in presence of a 610% tolerance affecting the four reactive parameters: black = upper and lower
boundaries, gray = nominal solution (a) @ V , (b) @ V .
analysis of the transfer function in terms of its poles and zeroes Fig. 16. Load step-up transient Vo response, comparison among selected solu-
especially because of its involved symbolic form. tions @ V .
The features of the best solutions in frequency domain can be
also analyzed by looking at Figs. 10 and 11 for both values.
In order to distinguish the different curves corresponding to the , even if such solutions show a higher quality factor (see
solutions of the Pareto front, only four out of eight of them have Fig. 7). The best solutions in terms of values obtained at
been depicted in Fig. 11 for . and at , thus the solutions #8 and #4 respec-
As it can be seen in Fig. 10, the transfer functions corre- tively, have been also compared in the frequency domain with
sponding to the four solutions at show a much better the ripple based design solutions.
shape with respect to those ones depicted in Fig. 4 and refer- Diagrams shown in Fig. 12 confirm that each solution opti-
ring to ripple-based designed solutions. There are no high res- mized for a given input voltage value (#4 @ and #8 @
onance peaks and the dynamics is very similar to that one of a ) exhibits the best behavior at that input voltage value
second-order system. By looking at the sequence of solutions, (#8 in Fig. 12(a) and #4 in Fig. 12(b)). Nevertheless, even if
from #1 to #4, a slight increase in the first resonance peak is they deteriorate at the value different with respect to that
observed, which is related to the increasing efficiency (as evi- one for which they are optimized (#8 in Fig. 12(b) and #4 in
denced in Table II), while the second resonance peak tends to Fig. 12(a)), they keep better than those ones obtained by means
disappear, as it was expected due to the higher PCA threshold of the ripple based approach. Figs. 13 and 14 show that solu-
value. The same comments can be referred to the solutions at tion #3 in both Pareto fronts (see Fig. 7), corresponding to the
2876 IEEE TRANSACTIONS ON POWER ELECTRONICS, VOL. 23, NO. 6, NOVEMBER 2008
same set of parameters at both and , represents [4] J. Betten and R. Kollman, Underutilized SEPIC Outperforms the Fly-
a good compromise and exhibits a small sensitivity with respect back Topology Courtesy of PlanetAnalog 07/05/2005.
[5] L. H. Dixon, High Power Factor Preregulator Using the SEPIC Con-
to large uncertainties, e.g., 10%, affecting inductance and ca- verter Unitrode seminar SEM900, Topic 6, 1993.
pacitance values. [6] A. Huang, The Power Management of PDA—The Application of
SEPIC Circuit Analog Integrations Corporation Application Note,
The four solutions at have been also compared in the AN020, Apr. 2001.
time domain by using PSIM software: a control loop ensuring a [7] S. Rahman and F. C. Lee, “Nonlinear program based optimization
30-kHz crossover frequency and 52 of phase margin for each of boost and buck-boost converter designs,” in Proc. PESC 81-IEEE
Power Elec. Spec. Conf., Boulder, CO, 1981, pp. 180–191.
one of them [28] has been designed. A load step from 25 to [8] S. Balachandran and F. C. Lee, “Algorithms for power converter design
250 mA has been applied and in Fig. 15 the transient response optimization,” IEEE Trans. Aerosp. Electron. Syst., vol. AES-17, no. 3,
given by the Pareto-optimal combinations is compared that one pp. 422–432, 1981.
[9] C. J. Wu, F. C. Lee, S. Balachandran, and H. L. Goin, “Design opti-
corresponding to the traditional design. The higher cross-over mization for a half-bridge dc-dc converter,” IEEE Trans. Aerosp. Elec-
frequency shown by the four solutions improves the settling tron. Syst., vol. AES-18, no. 4, pp. 497–508, 1982.
time during the transient, while the higher resonance peak of [10] S. Busquets-Monge, G. Sormekun, and E. Hertz, “Power converter de-
sign optimization: A GA-based design approach to optimization of
the ripple-based solution has a detrimental effect on the corre- power electronics circuits,” IEEE Ind. Applicat. Mag., Jan./Feb. 2004.
sponding response even if has a 12 dB attenuation (Fig. 5(a)). [11] J. M. Kwon, W. Y. Choi, J. J. Lee, E. H. Kimand, and B. H. Kwon,
A final look to solutions #1 to #4 at evidences a “Continuous-conduction-mode SEPIC converter with low reverse-re-
covery loss for power factor correction,” IEE Proc. Electr. Power Appl.,
reduction of the voltage ripple due to the increasing value of vol. 153, no. 5, September 2006.
the inductance and the reduced effect of the second resonance [12] G. Spiazzi and P. Mattavelli, “Design criteria for power factor prereg-
peak in terms of the pseudo-oscillation inspired by the transient ulators based on SEPIC and Cuk converters in continuous conduction
mode,” in Proc. IEEE IAS Conf, 1994, pp. 1084–1089.
(Fig. 16), while the total voltage undershoot remains almost the [13] N. D. Muhamad, M. R. Sahid, A. H. M. Yatim, N. R. N. Idris, and M.
same. S. Ayob, Design of Power Stage and Controller for DC-DC Converter
In conclusion the fourth solution shows a better frequency Systems Using PSPICE 0-7803-9296-5/05 2005 © IEEE, pp. 903–908.
[14] N. Femia, M. L. Cruoglio, G. Lisi, A. DeNardo, M. Slocchi, and W.
response and a higher efficiency than the others, but this is ob- Berns, “Novel design approach for POL buck DC-DC converters,” in
tained with higher inductance values of the two inductors, PCIM 2006.
and , whilst the first solution still guarantee good perfor- [15] N. Femia and M. Slocchi, “Novel mathematical approach for exploring
point of load converter design,” in Proc. PCIM 2006.
mances by using the lowest values of inductance. The solutions [16] N. Femia and G. Spagnuolo, “State-space models and order reduction
between the first one and the fourth one represent different trade- for DC-DC switching converters in discontinuous modes,” IEEE Trans.
offs. Power Electronics, vol. 10, no. 6, Nov. 1995.
[17] K. Deb, Multi-Objective Optimization Using Evolutionary Algo-
rithms. New York: Wiley, 2001.
VI. CONCLUSION [18] B. C. Moore, “Principal component analysis in linear systems: Control-
lability, observability, and model reduction,” IEEE Trans. Automatic
In this paper, the problem of the selection of passive compo- Contr., vol. AC-26, 1981.
nents of the power stage of fourth-order dc-to-dc converters is [19] D. Schelle and J. Castorena, Buck-Converter Design Demystified. :
Power Electronics Technology, June 2006, pp. 46–53.
addressed. A novel design procedure based on Principal Com- [20] N. Mohan, T. Undeland, and W. Robbins, Power Electronics: Con-
ponents Analysis has been introduced, which improves the con- verters, Applications and Design. New York: Wiley, 2003, pp.
trol-to-output transfer function of the converter. Passive com- 195–196.
[21] R. W. Erickson and D. Maksimović, Fundamentals of Power Elec-
ponents are selected by PCA-based design so that the frequency tronics, 2nd ed. New York: Kluwer Academic, 2001.
responses of the fourth-order converter, which exhibit couples of [22] R. B. Ridley, Analyzing the Sepic Converter. : Power Systems Design
complex poles and zeroes and complicated dependencies on cir- Europe, Nov. 2006, pp. 14–18.
[23] S. B. Yaakov, D. Adar, and G. Rahav, A SPICE Compatible Behav-
cuit parameters, is made similar to a second-order system. The ioral Model of SEPIC Converters 0-7803-3500-7/96, 1996 © IEEE,
PCA-based design method is supported by a multiobjective op- pp. 1668–1674.
timization algorithm which works on a database of real compo- [24] D. Adar, G. Rahav, and S. B. Yaakov, “Behavioral average model of
SEPIC converters with coupled inductors,” Electron. Lett., vol. 32, no.
nents taken from manufacturers’ components lists. PCA-based 17, pp. 1525–1526, Aug. 15, 1996.
design solutions obtained for SEPIC converter result to be much [25] R. B. Ridley, “A new, continuous-time model for current-mode con-
better compared with those ones given by the usual ripple-based trol,” IEEE Trans. Power Electron., vol. 8, pp. 271–280, Apr. 1991.
[26] L. H. Dixon, Control Loop Design, SEPIC Preregulator Example Uni-
design approach, both in time and frequency domains. trode Seminar SEM900, Topic 7, 1993.
[27] W. Gu, Small Signal Modeling for Current Mode Controlled Cuk and
SEPIC Converters 0-7803-8975-1/05, IEEE PEDC 2005, pp. 906–910.
REFERENCES [28] Venable and H. Dean, “The K factor: A new mathematical tool for
[1] L. C. Gomes de Freitas, M. G. Simoes, C. A. Canesin, and L. C. de stability analysis and synthesis,” in Proc. Powercon 10, San Diego, CA,
Freitas, “Programmable PFC based hybrid multipulse power rectifier 1983, pp. H1-1–H1-12.
for ultra clean power application,” IEEE Trans. Power Electron., vol. [29] R. E. Kalman, “Irreducible realizations and the degree of a rational
21, no. 4, pp. 959–966, July 2006. matrix,” SIAM J. Appl. Math., vol. 13, no. 2, pp. 520–544, 1965.
[2] M. Chen and J. Sun, “Reduced-order averaged modeling of active- [30] E. Kreindler and P. E. Sarachik, “On the concept of controllability
clamp converters,” IEEE Trans. Power Electron., vol. 21, no. 2, pp. and observability of linear systems,” IEEE Trans. Automat. Contr., vol.
487–494, March 2006. AC-9, pp. 129–136, 1964.
[3] S. Chakraborty, A. K. Jain, and N. Mohan, “A novel converter topology [31] A. J. Laub, “Computation of system balancing transformations and
for multiple individually regulated outputs,” IEEE Trans. Power Elec- other applications of simultaneous diagonalization algorithms,” IEEE
tron., vol. 21, no. 2, pp. 361–369, March 2006. Trans. Automatic Contr., vol. AC-32, no. 2, Feb. 1987.
DE NARDO et al.: POWER STAGE DESIGN OF FOURTH-ORDER DC–DC CONVERTERS BY MEANS OF PCA 2877
[32] J. Klein, Synchronous Buck MOSFET Loss Calculation With Excel Marco Nicolò received the M.E. degree in electrical
Model Fairchild Semiconductor Application Note, AN-6005, 2006. engineering from the University of Salerno, Italy, in
the 2005.
In 2007 he joined National Semiconductor Corpo-
ration, Fuerstenfeldbruck, Germany, as Design Engi-
neer. Currently, he is working in the NSC Application
Design Center and his main interest is in power elec-
tronics.