收藏切换
A Cascade-Formed Accurate Quasi-PR Controller Realization by Pole-Zero Placement
收藏切换
PDF
Chao GAO1, Shuyu ZHANG2, Bo LU3, Wenlong DING4, Ka Nang LEUNG1, Poh Chiang LOH1
CPSS Transactions on Power Electronics and Applications | 2024, 9(1) : 51 - 62
Less
收藏切换
CPSS Transactions on Power Electronics and Applications | 2024, 9(1): 51-62
Papers
A Cascade-Formed Accurate Quasi-PR Controller Realization by Pole-Zero Placement
Full
Chao GAO1, Shuyu ZHANG2, Bo LU3, Wenlong DING4, Ka Nang LEUNG1, Poh Chiang LOH1
Affiliations
  • 1 The Chinese University of Hong Kong Department of Electronic Engineering Hong Kong China
  • 2 Stanford University School of Engineering Stanford 94305 United States
  • 3 Fudan University Academy for Engineering and Technology Shanghai 200433 China
  • 4 Shandong University School of Control Science and Engineering Jinan 250100 China
  • Chao Gao received the B.S. degree from Wuhan University, Wuhan, China, in 2018, and the M.S. degree from Huazhong University of Science and Technology, Wuhan, China, in 2021, both in electrical engineering. He is currently working toward the Ph.D. degree from the Chinese University of Hong Kong, Hong Kong, China. His research interests include modeling and control of power electronic converters.

    Shuyu Zhang received her B.Eng. degree in Electrical Engineering from Huazhong University of Science and Technology, Wuhan, China, in 2023. She is currently pursuing a Master’s degree in Electrical Engineering at Stanford University. Her research focuses on multilevel converters and the study of electromagnetic interference in Power Electronic Systems.

    Bo Lu received the B.S. and M.S. degrees from Huazhong University of Science and Technology, Wuhan, China, in 2018 and 2021, respectively, both in electrical engineering. He is currently working toward the Ph.D. degree from Fudan University, Shanghai, China. His research interests include high-precision and high-efficient power amplifiers.

    Wenlong Ding received the B.S. and M.S. degrees from Qufu Normal University, Rizhao, China, in 2011 and 2014, respectively, and the Ph.D. degree from Shandong University, Jinan, China, in 2019, all in electrical engineering. From November 2019 to August 2020, he was a Post-Doctoral Research Fellow with the Department of Electronic Engineering, The Chinese University of Hong Kong, Hong Kong, China. He was with the Hong Kong Applied Science and Technology Research Institute (ASTRI), Hong Kong, from September 2020 to November 2021. Since 2021, he has been with the School of Control Science and Engineering, Shandong University. His current research interests include applications of wide band gap (WBG) devices, multilevel converters, and battery charging/testing technology.

    Ka Nang Leung received the B.Eng., M.Phil., and Ph.D. degrees in Electrical and Electronic Engineering from Hong Kong University of Science and Technology (HKUST), Clear Water Bay, Hong Kong, in 1996, 1998, and 2002, respectively. In 2002, he was a Visiting Assistant Professor at HKUST. In 2005, he joined with the Department of Electronic Engineering, Chinese University of Hong Kong, Hong Kong, where he is currently an Associate Professor. His research interests include power-management integrated circuits and low-voltage low-power analog integrated circuits. Dr. Leung is the Chairman of the IEEE (Hong Kong) Electron Device/Solid-State Circuit Joint Chapter in 2012. He is a member in the Editorial Boards and a guest editor of a special issue of Energies. He serves as a Paper Reviewer in numerous IEEE journals and IEEE international conferences. Moreover, he is actively involved in the organization of IEEE international conferences. He is a co-recipient of the Best Paper Awards of TENCON in 2015 and the IEEE Student Symposium ED/SSC in 2011, 2014 and 2019.

    Poh Chiang Loh received the B.Eng. (hons.) and M.Eng. degrees from the National University of Singapore, Singapore, in 1998 and 2000, respectively, and the Ph.D. degree from Monash University, Melbourne, Vic, Australia, in 2002, all in electrical engineering. From 2013 to 2015, he was a Professor with Aalborg University, Aalborg, Denmark. Since 2015, he has been a Tenured Full Professor with the Chinese University of Hong Kong, Hong Kong, China. His research interests include power converters and their grid applications.

Published: 2024-03-10 doi: 10.24295/CPSSTPEA.2023.00046
Outline
收藏切换

PR controller has been widely researched in various control systems for its robustness and simplicity. However, a traditional PR controller with relatively small integral gain, used for higherorder harmonics to keep stability, will cause increases in magnitudes, and decreases in phase around resonant frequency, and jeopardize stability. These all call for a more precise realization of PR controller. This paper proposes a cascadeformed PR controller realization method, which proves to realize PR controller more accurately even with a relatively small integral gain. The method is to decompose a PR controller into multiple independent PR units, and each PR unit is realized by mapping PR controller's parameters to its pole and zero positions. The distance between a polezero pair is found related to frequency characteristic error and is restricted accordingly to limit the error. Comprehensive comparisons of PR controllers realized by cascade form and the traditional parallel form have been conducted theoretically and experimentally, verifying that the cascade realization method is more accurate.

Cascade form  /  pole-zero placement  /  PR controller  /  realization
Chao GAO, Shuyu ZHANG, Bo LU, Wenlong DING, Ka Nang LEUNG, Poh Chiang LOH. A Cascade-Formed Accurate Quasi-PR Controller Realization by Pole-Zero Placement[J]. CPSS Transactions on Power Electronics and Applications, 2024 , 9 (1) : 51 -62 . DOI: 10.24295/CPSSTPEA.2023.00046
I $\mathrm{N}$ recent decades, proportional resonant (PR) controllers have Lattracted increasing research attention [1], especially in power electronics such as motor drives [2], active power filters [3], grid-tied inverters with LCL filter [4], and so on, for its excellent robustness and simplicity [5],[6].
Traditionally, PR controllers are realized by parallel form, i.e., the parallel of a proportional (P) controller and one or multiple resonant (R) controllers [7]. R controller is derived initially from transforming synchronous frame proportional integrator (PI) controller into the stationary frame for integrating sinusoidal signal [2],[8]. Thus, the ideal R controller is also called general integrator (GI)[9] or sinusoidal signal integrator (SSI)[10] with infinite gain at its resonant frequency and a relatively small gain at other frequencies. The parallel of multiple R controllers provides infinite gain at multiple resonant frequencies, and thus zero tracking error for those frequencies is guaranteed [11]. Further, with the P controller added to R controllers, higher gain at frequencies other than resonant frequencies is obtained, broadening bandwidth [12]. Since any finite value is trivial when added to an infinite value, a parallel-formed ideal PR controller can maintain its frequency characteristics at a resonant frequency the same as that of the corresponding $\mathrm{R}$ controller. Therefore, the infinite gain of the $\mathrm{R}$ controller is the prerequisite of the accurate realization of the parallel-formed PR controller.
In practice, the usage of a damped PR (quasi-PR) controller, whose resonant peak is finite, is recommended for several reasons. First, an ideal PR controller acts like an infinite quality factor network and is unrealizable [13],[14]. Besides,[8],[15] find quasi-PR controller is less sensitive to fundamental frequency deviation because of the broader bandwidth of its resonant peaks. Meanwhile, its phase changes slower around resonant frequencies [16] alleviating the anomalous peaks of the closed-loop transfer function that incorporates a PR controller [17],[18]. Thus, a quasi-PR controller is widely used.
There are multiple ways to design a PR controller [19]-[21]. However, none of these design methods consider the deviation of frequency characteristic caused by replacing ideal R controller with damped R controller. The situation exacerbates with resonant frequency increasing. Multiple application needs PR controller to be tuned at relatively high frequency. For example, an active power filter (APF) usually needs a PR controller to be tuned at up to 13th order harmonic [22] to compensate the corresponding harmonic. More critically, when APF is used to filter switching harmonics [23], the resonant frequency of PR controller should be tuned at thousands of Hertz. Similarly, when a PR controller is applied to the motor drive to suppress torque ripple caused by distorted back electromotive force, it also needs to be tuned at ${6np}\left({n \in \mathbf{Z}}\right)$ times mechanical angle speed where $p$ is the number of pole pairs [24]. It is found that the integral gain must decrease with resonant frequency increasing to keep stability [25],[26], and lead to more error on frequency characteristics. Reference [26] proposes an iteration algorithm to correct the errors of phases, which takes a large amount of calculation and does not consider the magnitude errors.
In this paper, the inaccuracy of the parallel realization of the quasi-PR controller is first demonstrated and will be tackled by a new cascade realization form, which more accurately realizes the magnitudes and phases of quasi-PR controller at resonant frequencies. Since this paper deals with only quasi-PR controllers rather than ideal PR controllers, in the following parts, the prefix “quasi-” will be saved for conciseness, and the words “PR or R controller” refer to the quasi-PR or quasi-R controller.
This paper is organized as follows. In Section II, the inaccuracy of the PR controller’s parallel realization is demonstrated. In Section III, the PR controller’s pole-zero distribution characteristics and resonant decoupling property are formulated, based on which a cascade realization method for PR controller is proposed. Section IV presents the relevant experimental results and analysis, followed by Section V, the conclusion.
Since this paper only focuses on the new PR controller realization method, we use a commonly used L-filtered grid-tied converter to illustrate the performance of the controller without excluding the possibility of other application situations. The schematic diagram of the grid-tied converter with its controller is shown in Fig. 1, where $L$ and $R$ are the inductance and parasitic resistance of the L filter, respectively, ${i}_{\mathrm{o}}\left( s\right)$ is the output current, ${i}^{* }\left( s\right)$ is the reference current, $e\left( s\right)$ is the error signal, ${G}_{\mathrm{{PR}}}\left( s\right)$ is the PR controller which is traditionally realized by the parallel form. The transfer function from reference voltage, ${u}_{\text{ref }}\left( s\right)$, to output current, ${i}_{\mathrm{o}}\left( s\right)$, is:
$ P\left( s\right)= \frac{{i}_{\mathrm{o}}\left( s\right)}{{u}_{\text{ref }}\left( s\right)} =\frac{{e}^{-{1.5sT}}}{{sL}+ R}$
where $T$ is the sampling period, and ${e}^{-{1.5}\mathrm{\;s}T}$ represents the 1.5 sampling delay introduced by digital computation and PWM. A conventional parallel PR controller can be expressed as:
${G}_{\mathrm{{PR}}}\left( s\right)= {K}_{\mathrm{P}}+ \mathop{\sum }\limits_{h}{G}_{\mathrm{R}h}\left( s\right)$
where ${K}_{\mathrm{P}}$ is the proportional gain, and ${G}_{\mathrm{R}h}\left( s\right)$ is the $\mathrm{R}$ controller tuned at $h$ order harmonic, which can be expressed as:
${G}_{\mathrm{R}h}\left( s\right)= {K}_{\mathrm{I}h}\frac{2{\omega }_{\mathrm{c}h}\left\lbrack {s\cos \left({\varphi }_{h}\right)- h{\omega }_{1}\sin \left({\varphi }_{h}\right)}\right\rbrack }{{s}^{2}+ 2{\omega }_{\mathrm{c}h}s +{h}^{2}{\omega }_{1}^{2}}$
where ${K}_{\mathrm{I}h}$ is the integral gain, i.e., the resonant magnitude of ${G}_{\mathrm{R}h}\left( s\right)$, and ideally it should also be the resonant magnitude of ${G}_{\mathrm{{PR}}}\left( s\right),{\omega }_{\mathrm{c}h}$ is the resonant cutoff frequency, ${\varphi }_{h}$ is the phase lead of ${G}_{\mathrm{R}h}\left( s\right)$, and is also expected to be the phase of ${G}_{\mathrm{{PR}}}\left( s\right)$ at the corresponding resonant frequency, $h{\omega }_{1}$. By solving the denominator of (3), the poles of the PR controller are given by:
${s}_{1,2}= -{\omega }_{ch}\pm j\sqrt{{h}^{2}{\omega }_{1}^{2}- {\omega }_{ch}^{2}}$
At frequencies other than resonant frequencies, the PR controller can be regarded as a proportional controller, and all $\mathrm{R}$ controllers can be neglected. Thus, ${K}_{\mathrm{P}}$ decides the bandwidth of the controller with the following relationship:
${K}_{\mathrm{p}}= {\alpha }_{\mathrm{e}}L $
where ${\alpha }_{\mathrm{c}}$ is the bandwidth of the controller and usually set to ${\omega }_{\mathrm{s}}/{10}$(${\omega }_{\mathrm{s}}$ represents the sampling angle frequency). With parameters listed in Table I, ${K}_{\mathrm{P}}$ is calculated as 15.7. Phase lead is used to compensate for the 1.5 sampling delay [26], and improve phase margin, i.e.,
${\varphi }_{h}= {1.5}\times {2\pi }{f}_{1}{hT}$
where ${f}_{1}$ is the grid fundamental frequency.
${K}_{\mathrm{I}h}$ can be approximated as the resonant magnitude [16]. Thus, it should be set high enough for a small tracking error [16]. Meanwhile, at high frequencies, a high ${K}_{\mathrm{I}h}$ may critically reduce the stability margin and response speed [25],[26]. Thus, the tuning of ${K}_{\mathrm{I}h}$ differs in resonant frequency and should strike a balance between small steady-state error and good stability.
A parallel PR controller tuned at $h \in \{ 1,3,5,\ldots,{19}\}$ order harmonics is realized aiming to resist grid voltage harmonics at those frequencies. ${K}_{\mathrm{I}h}$ is set to 100, and other parameters can be found in Table I. The Bode plot of the parallel PR controller is depicted in Fig. 2. It can be obviously found that at lower resonant frequencies, the resonant magnitudes exceed ${K}_{\mathrm{I}h}({100}$ or ${40}\mathrm{\;{dB}}$), which is circled in red. For clearer investigation Table II lists the expected resonant magnitude $\left({K}_{1h}\right)$, and phase lead $\left({\varphi }_{h}\right)$, and actual values of the resultant parallel PR controller $\left({{M}_{\mathrm{{PR}}h},{\varphi }_{\mathrm{{PR}}h}}\right)$ at selected 1st,3rd,17th, and 19th order harmonics to see its characteristics at low and high frequency. It can be found that the resonant magnitude and phase of parallel PR controller deviates from the expected value a lot. PR controller’s magnitude, ${M}_{\mathrm{{PR}}h}$, exceeds ${K}_{\mathrm{I}h}$ by up to 15.2 at $h = 1$. PR controller’s phase, ${\varphi }_{\mathrm{{PR}}h}$, lags ${\varphi }_{h}$ by about ${9.2}^{\circ }$ at $h =$ 19, and this may jeopardize stability.
Fig. 3 illustrates the magnitude error and phase error of a PR controller at resonant frequency. Each vector represents a complex value of either an $\mathrm{R}$ controller, or a $\mathrm{P}$ controller, or their summation, the PR controller. The ${K}_{\mathrm{I}}$ of ${G}_{\mathrm{R}}\left( s\right)$ or ${G}_{\mathrm{{PR}}}\left( s\right)$ in dashed red line is higher than that of ${G}_{\mathrm{R}}\left( s\right)$ or ${G}_{\mathrm{{PR}}}\left( s\right)$ in solid black line. Phase error, ${E}_{\varphi }= \angle {G}_{\mathrm{{PR}}}\left( s\right)- \angle {G}_{\mathrm{R}}\left( s\right)$, is the angle between ${G}_{\mathrm{{PR}}}\left( s\right)$ and ${G}_{\mathrm{R}}\left( s\right)$, as marked in Fig. 3. It can be seen that $\left|{E}_{\varphi }\right|$ increase as ${K}_{\mathrm{I}}$ decreases. On the other hand, magnitude error, ${E}_{\mathrm{M}}$, is the difference between $\left|{G}_{\mathrm{{PR}}}\right|$ and $\left|{G}_{R}\right|$, which is obtained graphically by drawing a circle (dotted line) centered at the endpoint of ${G}_{\mathrm{{PR}}}\left( s\right)$ with a radius of ${K}_{\mathrm{I}}$. Obviously, $\left|{E}_{\mathrm{M}}\right|$ also increase as ${K}_{\mathrm{I}}$ decreases. Thus, for high-order harmonics where ${K}_{1h}$ is set small, the traditional parallel realization will result in a relatively significant increase in magnitude and decrease in phase. A new realization form is proposed to tackle this problem in the following section.
In the analysis of the parallel PR controller (addition form), we always try to ignore the smaller part in the summation, which needs to be more accurate. In the analysis of the transfer function in cascade form (multiplication form), the synthesized frequency characteristics will be the exact the summation of all frequency characteristics of multiplier transfer functions since $\angle {G}_{1}{G}_{2}= \angle {G}_{1}+ \angle {G}_{2}$, and ${20}\log \left|{{G}_{1}{G}_{2}}\right|= {20}\log \left|{G}_{1}\right|+$ ${20}\log \left|{G}_{2}\right|$. Thus, in this section, we try to analyze the PR controller in cascade form and find ways to realize a PR controller directly in cascade form.
Fig. 4 shows a typical pole-zero map of a PR controller tuned at only two resonant frequencies, ${\omega }_{1}$ and $3{\omega }_{1}$. It can be found that poles are located around $\pm {jh}{\omega }_{1}$ as specified by (4). Zeros are extremely close to their corresponding poles with respect to other poles (note the scale of the real and imaginary axes). In fact, it is a common characteristic of PR controller’s pole-zero map that poles and zeros are located closely in pairs around resonant frequencies. Thus, a PR controller can be expressed as:
${G}_{\mathrm{{PR}}}\left( s\right)= {K}_{\mathrm{P}}\mathop{\prod }\limits_{i}\frac{s -{z}_{i}}{s -{p}_{i}}$
where ${z}_{i}$ and ${p}_{i}$ are a pair of zero and pole tuned at the same frequency.
According to (7), the frequency characteristics of the PR controller at $\omega$ will be:
${G}_{\mathrm{{PR}}}\left({j\omega }\right)= {K}_{\mathrm{P}}\frac{\left({{j\omega }- {z}_{1}}\right)\left({{j\omega }- {z}_{3}}\right)\left({{j\omega }- {\bar{z}}_{1}}\right)\left({{j\omega }- {\bar{z}}_{3}}\right)}{\left({{j\omega }- {p}_{1}}\right)\left({{j\omega }- {p}_{3}}\right)\left({{j\omega }- {\bar{p}}_{1}}\right)\left({{j\omega }- {\bar{p}}_{3}}\right)} $
where $\bar{x}$ means the conjugate value of $x$. All poles and zeros can be calculated by simply setting the denominator or numerator of (2) to 0 . The difference between $s ={j\omega }$ and a pole or zero is denoted as a vector, ${v}_{x}, x \in \{$ poles, zeros $\}$. For clarity, $-{v}_{x}$ is shown in Fig. 4 with the reverse direction of ${v}_{x}$. Therefore,(8) becomes:
$ G\left({j\omega }\right)= {K}_{\mathrm{P}}\frac{{v}_{z1}{v}_{z3}{v}_{\bar{z}1}{v}_{\bar{z}3}}{{v}_{p1}{v}_{p3}{v}_{\bar{p}1}{v}_{\bar{p}3}}$
When $s ={j\omega }$ is at the position (see Fig. 4 ) far away from all poles and zeros (with respect to the distance within a pole-zero pair), a pair of pole and zero can be regarded as overlapping, and ${v}_{zx}/{v}_{px}$ is approximated to 1, where ${z}_{x}$ and ${p}_{x}$ are a pair of zero and pole. Then, ${G}_{\mathrm{{PR}}}\left({j\omega }\right)$ is approximated to:
${G}_{\mathrm{{PR}}}\left({j\omega }\right)\approx {K}_{\mathrm{P}}$
For example, ${j2}{\omega }_{1}$ is far away from the resonant frequency, $j{\omega }_{1}$ and ${j3}{\omega }_{1}$ where poles and zeros are located. And ${K}_{\mathrm{P}}= {15.7}$ is highly close to the actual value ${G}_{\mathrm{{PR}}}\left({{j2}{\omega }_{1}}\right)= {15.62}\angle -{0.65}^{\circ }$ confirming the accuracy of the approximation (10).
Moreover, if $s$ is close to a pole-zero pair, the assumption of ${v}_{zx}/{v}_{px}\approx 1$ for that pole-zero pair will not hold, and (9) will be approximated to:
${G}_{\mathrm{{PR}}}\left({j\omega }\right)\approx {K}_{\mathrm{P}}\frac{{v}_{zx}}{{v}_{px}}= {K}_{\mathrm{P}}\frac{{j\omega }- {z}_{x}}{{j\omega }- {p}_{x}}$
For example, when $s = j{\omega }_{1},\left({11}\right)$ becomes:
${G}_{\mathrm{{PR}}}\left({j{\omega }_{1}}\right)\approx {K}_{\mathrm{p}}\frac{j{\omega }_{1}- {z}_{1}}{j{\omega }_{1}- {p}_{1}}= {116.2}\angle {5.01}^{\circ }$
which is very close to the actual value, ${G}_{\mathrm{{PR}}}\left({j{\omega }_{1}}\right)= {115.58}$ $\angle {4.7}^{\circ }$ confirming the accuracy of the approximation (11).
To summarize,(10) and (11) effectively approximate PR controller’s frequency characteristics. When considering frequency far from resonant frequencies, poles and zeros align, PR controller is simplified to ${K}_{\mathrm{P}}$, i.e.,(10). When considering frequency near a resonant frequency, the corresponding pole-zero pair dominate, simplifying the PR controller to ${K}_{\mathrm{P}}\left({s - z}\right)/\left({s - p}\right)$, i.e.,(11). Thus, $\left({s - z}\right)/\left({s - p}\right)$ can be regarded as a PR controller unit (PRU) determining a resonant peak, and (7) can be rewritten as:
${G}_{\mathrm{{PR}}}\left( s\right)= {K}_{\mathrm{P}}\prod \mathrm{{PRU}}$
In real-coefficient PR controller, PRUs always appear in conjugate pairs. Each PRU, represented as $\left({s -{z}_{h}}\right)/\left({s -{p}_{h}}\right)$, and its conjugate, $\left({s -{\bar{z}}_{h}}\right)/\left({s -{\bar{p}}_{h}}\right)$, regulate the positive and negative sequences, respectively. If the conjugate PRU is not used, only the positive or negative sequence will be regulated resulting in a reduced order generalized integrator (ROGI)[27].
This property of the PR controller, termed “resonant decoupling” in this paper, offers insights, and allows one to realize the frequency characteristics of a PR controller by cascading individual PRUs, as depicted in Fig. 1. This concept will be detailed in the following subsections. In the next subsection, a detailed examination of the resonant decoupling property will be conducted, and the errors associated with its application are assessed.
The reason that omitting a pole-zero pair will cause difference on frequency characteristics is that the pole and zero do not align precisely. Consider a pole-zero pair: the zero, ${z}_{k}= {r}_{z}+ j{\omega }_{z}$, and the pole, ${p}_{k}= {r}_{p}+ j{\omega }_{p}$, resonant at ${\omega }_{p}= k{\omega }_{1}\approx {\omega }_{z}$. According to (7), the phase at the frequency $\omega$, which is far away from ${\omega }_{\mathrm{p}}$, can be expressed as
$\angle {G}_{\mathrm{{PR}}}\left({j\omega }\right)= \mathop{\sum }\limits_{{i \neq k}}\angle \frac{{j\omega }- {z}_{i}}{{j\omega }- {p}_{i}}+ \angle \frac{{j\omega }- {z}_{k}}{{j\omega }- {p}_{k}}= {\varphi }^{* }+ \bigtriangleup {\varphi }_{k}$
where ${\varphi }^{* }$ denotes the phase $\angle {G}_{\mathrm{{PR}}}\left({j\omega }\right)$ with pole-zero pair(${p}_{k}$, ${z}_{k}$), omitted, and $\Delta {\varphi }_{k}$ denotes the phase variation caused by $\left({{p}_{k},{z}_{k}}\right)$. When resonant decoupling property is implemented, the pole-zero pair $\left({{p}_{k},{z}_{k}}\right)$, is omitted, and therefore, a phase error, $-\Delta {\varphi }_{k}$, will be caused on frequency characteristics at $\omega$. As shown in Fig. 5, $\Delta {\varphi }_{k}$ is the angle formed by two vectors $\left({s -{z}_{k}}\right)$ and $\left({s -{p}_{k}}\right).{\gamma }_{p}$ and ${\gamma }_{z}$ are relatively small and can be approximated by their tangent values, which yields
$\bigtriangleup {\varphi }_{k}= {\gamma }_{z}- {\gamma }_{p}\approx \tan \left({\gamma }_{z}\right)- \tan \left({\gamma }_{p}\right)\\= \frac{\left|{r}_{z}\right|}{\omega -{\omega }_{z}}- \frac{\left|{r}_{p}\right|}{\omega -{\omega }_{p}}\approx \frac{{r}_{p}- {r}_{z}}{\omega -{\omega }_{p}}$
(15) indicates that $\Delta {\varphi }_{k}$ is in proportion to $\left({{r}_{p}- {r}_{z}}\right)$ and in inverse proportion to $\left({\omega -{\omega }_{p}}\right)$. Fig. 6(a) shows the relationship between the absolute phase variation, $\left|{\Delta {\varphi }_{k}}\right|$, and $\left|{\operatorname{Real}\left({{z}_{k}- {p}_{k}}\right)}\right|$ and $\left|{\omega -{\omega }_{p}}\right|$ graphically.
So, the difference between the real parts of a pole-zero pair, $\operatorname{Real}\left({{z}_{k}- {p}_{k}}\right)= \left({{r}_{p}- {r}_{z}}\right)$, should be restricted so that the phase error, $-\Delta {\varphi }_{k}$, is acceptable. The maximum allowable phase draft can be uniformly set to $\pm {0.02}$ rad $\left({\pm {1.15}^{\circ }}\right)$ at resonant frequency, $h{\omega }_{1}$, for all $h$, i.e.
$\left|{\Delta {\varphi }_{k}}\right|= \left|\frac{\operatorname{Real}\left({{z}_{k}- {p}_{k}}\right)}{h{\omega }_{1}- k{\omega }_{1}}\right|< {0.02}$
Consequently
$- {0.02}\min \bigtriangleup \omega <\operatorname{Real}\left({{z}_{k}- {p}_{k}}\right)< {0.02}\min \bigtriangleup \omega $
where $\min {\Delta \omega }= \min \left|{h - k}\right|{\omega }_{1}$. Usually, resonant frequencies are set first, i.e., the $\min {\Delta \omega }$ is known. Thus,(17) restricts the real part of the zeros.
On the other hand, the magnitude of ${G}_{\mathrm{{PR}}}\left(\mathrm{s}\right)$ at $\omega$ is
$\left|{{G}_{\mathrm{{PR}}}\left({j\omega }\right)}\right|= {K}_{\mathrm{p}}\mathop{\prod }\limits_{{j \neq k}}\left|\frac{{j\omega }- {z}_{\mathrm{j}}}{{j\omega }- {pj}}\right|\cdot \left|\frac{{j\omega }- {z}_{k}}{{j\omega }- {p}_{k}}\right|$
where the factor ${K}_{\mathrm{p}}\mathop{\prod }\limits_{{j \neq k}}\left|\frac{{j\omega }- {z}_{\mathrm{j}}}{{j\omega }- {pj}}\right|$ denoted as ${M}^{* }$ is the magnitude $\left|{{G}_{\mathrm{{PR}}}\left({j\omega }\right)}\right|$ with pole-zero pair, $\left({{p}_{k},{z}_{k}}\right)$, omitted. As for another part
$\left|\frac{{j\omega }- {z}_{k}}{{j\omega }- {p}_{k}}\right|= \left|\frac{{r}_{z}+ j\left({{\omega }_{z}- \omega }\right)}{{r}_{p}+ j\left({{\omega }_{p}- \omega }\right)}\right|$
when $\omega$ is far away from ${\omega }_{p},{r}_{z}$ and ${r}_{p}$ are far smaller than $\left|{{\omega }_{p}- \omega }\right|$, or $\left|{{\omega }_{z}- \omega }\right|$, they can be neglected, and the magnitude of ${G}_{\mathrm{{PR}}}\left({j\omega }\right)$ is approximated to
$\left|{{G}_{\mathrm{{PR}}}\left(\omega \right)}\right|\approx {M}^{* }\left|\frac{{\omega }_{z}- \omega }{{\omega }_{p}- \omega }\right|= {M}^{* }\left({1 +\frac{{\omega }_{z}- {\omega }_{p}}{{\omega }_{p}- \omega }}\right)$
Then the variation in the magnitude of ${G}_{\mathrm{{PR}}}\left({j\omega }\right)$ is
$\bigtriangleup \left|{{G}_{\mathrm{{PR}}}\left({j\omega }\right)}\right|= \left|{{G}_{\mathrm{{PR}}}\left({j\omega }\right)}\right|- {M}^{* }= {M}^{* }\frac{{\omega }_{z}- {\omega }_{p}}{{\omega }_{p}- \omega }$
which means the variation of the magnitude, $\Delta \left|{{G}_{\mathrm{{PR}}}\left({j\omega }\right)}\right|$, is in proportion to ${M}^{* }$ and a factor $\left({{\omega }_{z}- {\omega }_{p}}\right)/\left({{\omega }_{p}- \omega }\right)$, which is introduced by a pole-zero pair, $\left({{p}_{k},{z}_{k}}\right)$, tuned at $k{\omega }_{1}$. Fig. 6(b) illustrates the relationship between the absolute magnitude variation in per unit, $\left|\frac{\bigtriangleup \left|{{G}_{\mathrm{{PR}}}\left({{jh}{\omega }_{1}}\right)}\right|}{{M}^{* }}\right|$, and the variations in $\left|{\operatorname{Imag}\left({{z}_{k}- {p}_{k}}\right)}\right|$ and $\left|{\omega -{\omega }_{p}}\right|$.
The bias of the imaginary part of a zero, Imag$\left(z_{k}-p_{k}\right)=\left(\omega_{z}-\omega_{p}\right)$, should be restricted such that the omitted pole-zero pair does not affect the magnitudes of other resonant frequencies too much. The maximum allowable deviation ratio of $\left|{{G}_{\mathrm{{PR}}}\left({j\omega }\right)}\right|$ can be set uniformly to $2\%$ at resonant frequency, $h{\omega }_{1}$, for all $h$, i.e.
$\left|\frac{\bigtriangleup \left|{{G}_{\mathrm{{PR}}}\left({{jh}{\omega }_{1}}\right)}\right|}{M}\right|= \left|\frac{\operatorname{Imag}\left({{z}_{k}- {p}_{k}}\right)}{k{\omega }_{1}- h{\omega }_{1}}\right|< 2\%$
Then
$- {0.02}\min \bigtriangleup \omega <\operatorname{Imag}\left({{z}_{k}- {p}_{k}}\right)< {0.02}\min \bigtriangleup \omega $
(17) and (23) restrict zeros within squares centered at their paired poles. In this way, the maximum effect by omitting a pole-zero pair to frequency characteristics at other resonant frequencies is limited to $\pm 2\%$ magnitude error and $\pm {0.02}$ rad phase error. However, note that when omitting multiple pole-zero pairs, the error will be superimposed, and may exceed the abovementioned limit. But this should not be a worry since (17) and (23) use $\min {\Delta \omega }$, which corresponds to the adjacent pole-zero pair to be omitted. As for further pole-zero pair, its distance to the frequency of concern is much larger, and the influence will be small enough to be omitted.
The problem that paralleling scheme fails to realize the PR controller accurately will be tackled by cascade realization in this subsection, which can realize the PR controller more accurately.
A PR controller can be solely determined by ${K}_{\mathrm{P}},{K}_{\mathrm{I}h},{\varphi }_{h}$, and ${\omega }_{\mathrm{c}h}$, provided that $h \in H$ and resonant frequencies $h{\omega }_{1}$ are known. Those parameters appear explicitly in the expressions of parallel-formed PR controllers (2) and (3), and should be tuned directly by traditional design methods as summarized in Section II. However, these parameters do not explicitly show in the expression of cascade-formed PR controller (7), except for ${K}_{\mathrm{P}}$. Thus, by mapping all those parameters to the locations of poles and zeros, a cascade realization of a PR controller can be done.
The location of a pole can be obtained directly since their imaginary part is the resonant frequency $h{\omega }_{1}\left\lbrack {18}\right\rbrack$, and its real part is $-{\omega }_{\mathrm{{ch}}}$ according to (4).
The locations of zeros are determined by desired frequency characteristics around resonant frequencies. According to the resonance-decoupling property of PR controller, the frequency response around a particular resonant frequency $k{\omega }_{1}$ can be approximated by
${G}_{\mathrm{{PR}}}\left( s\right)\approx {K}_{\mathrm{p}}\frac{s -{z}_{k}}{s -{p}_{k}}$
where ${p}_{k}= -{\omega }_{\mathrm{c}k}+ {jk}{\omega }_{1},{z}_{k}= {r}_{k}+ j{\omega }_{k}$. The magnitude of ${G}_{\mathrm{{PR}}}\left( s\right)$ approaches its maximum at resonant point, $s ={jk}{\omega }_{1}$, which should be the integral gain ${K}_{\mathrm{I}k}$, i.e.
${K}_{\mathrm{I}k}= {K}_{\mathrm{p}}\left|\frac{{jk}{\omega }_{1}- {z}_{k}}{{jk}{\omega }_{1}- {p}_{k}}\right|$
namely
$\left|{{jk}{\omega }_{1}- {z}_{\mathrm{k}}}\right|= \frac{{K}_{\mathrm{I}k}{\omega }_{\mathrm{c}k}}{{K}_{\mathrm{p}}}$
which defines a trajectory of ${z}_{k}$, a circle centered at ${jk}{\omega }_{1}$ with a radius of ${K}_{\mathrm{I}k}{\omega }_{\mathrm{c}k}/{K}_{\mathrm{P}}$ (see blue line in Fig. 7 ). In addition, the phase at the resonant frequency, $k{\omega }_{1}$, should be ${\varphi }_{k}$:
${\varphi }_{k}= \angle {G}_{\mathrm{{PR}}}\left({{jk}{\omega }_{1}}\right)\approx \angle \left(\frac{{z}_{k}- {jk}{\omega }_{1}}{{p}_{k}- {jk}{\omega }_{1}}\right)$
$\angle \left(\frac{{z}_{k}- {jk}{\omega }_{1}}{{p}_{k}- {jk}{\omega }_{1}}\right)$ is the angle $\angle \mathrm{{ZRP}}$ as shown in Fig. 7. Thus, ${z}_{k}$
is located on the half-line starting from ${jk}{\omega }_{1}$ and form an angle ${\varphi }_{k}$ counterclockwise with $\left({{p}_{k}- {jk}{\omega }_{1}}\right)$ (see green line in Fig. 7 ). Combining the magnitude and phase conditions together, ${z}_{k}$ is at the intersection of the circle (26) and the half-line (27). For now, a pair of pole and zero are located, namely a PRU, $\left({s -{z}_{k}}\right)/\left({s -{p}_{k}}\right)$ is determined. The conjugate PRU, $\left({s -{\bar{z}}_{k}}\right)/$ $\left({s -{\bar{p}}_{k}}\right)$, can be obtained immediately if needed. By repeating the above steps, other PRUs tuned at other frequencies can be realized. Finally, the PR controller is obtained by cascading all PRUs and a ${K}_{\mathrm{P}}$ in the form of (7) or (13). From the standpoint of how subordinate controllers are linked, this proposed implementation method connects PRUs in a cascading manner differing from the conventional approach where P controllers and $\mathrm{R}$ controllers are connected in parallel. Therefore, this implementation can be seen as a dual form of the traditional parallel setup. It’s worth noting that the PRU doesn’t necessarily have to function as a PR controller. It can also be designed as a notch filter to block signals at a specific frequency. An ideal notch filter exhibits zero magnitude at its resonant frequency, meaning ${K}_{\mathrm{I}k}= 0$ and ${z}_{k}= {jk}{\omega }_{1}$.
In the previous subsection, how a zero should stay within a square region centered on its paired pole, with both sides measuring ${0.04}\min {\Delta \omega }$ (as illustrated by the red line in Fig. 7 ) is discussed. To maintain both magnitude and phase unchanged, this condition can be met by adjusting ${\omega }_{\mathrm{c}h}$. According to the magnitude condition (26), the location of the zero is also related to ${\omega }_{\mathrm{c}h}$. With ${K}_{\mathrm{I}h}$ set to ${100}, h \in \{ 1,3,\ldots,{19}\}$ and other parameters listed in Table I, the maximal values for ${\omega }_{\mathrm{c}h}$ for all $h$ are obtained (the minimal values are zero). The minimum of these maximum values of ${\omega }_{\mathrm{c}h}$ is calculated to be 0.95 justifying the suitability of the value 1 assigned to ${\omega }_{\mathrm{c}h}$ in Section II.
Since controllers are all discrete for digital implementation, we ultimately need to obtain a PR controller in the $z$ domain. Conventional discretization method for PR controller is Tustin with prewarping, which discretizes PR controllers by discretizing each $\mathrm{R}$ controller first and then paralleling them. This inherently suits parallel-formed PR controller. If cascade-formed PR controller adopts Tustin with prewarping to discretize, it has to be decomposed into partial-fraction form or parallel form. That would be laborious. Considering this aspect, the cascade realization of PR controller in $\mathrm{z}$ domain will be elaborated in this subsection, which can be done directly without the prior $s$ domain realization.
As concluded in the $s$ domain, the pole-zero-gain form of a discrete PR controller is:
${G}_{\mathrm{{PR}}}\left( z\right)= {K}_{\mathrm{p}}\mathop{\prod }\limits_{i}\frac{z -{z}_{i}}{z -{p}_{i}}$
where ${z}_{i}$ and ${p}_{i}$ are a pair of zero and pole tuned at the same frequency. The frequency response of (28) can be obtained by substituting $z$ with ${e}^{j\omega T}$, symmetrically to the $s$ domain where $s$ is substituted by ${j\omega }$. Geometrically this means $z$ varies along the unit circle versus $s$ along the imaginary axis. Thus, all the above analyses in the $s$ domain can be applied to the $z$ domain with a bit of adaptation.
Firstly, the pole should be placed at ${e}^{-{\omega }_{ch}T +{jh}{\omega }_{t}T}$, corresponding to ${p}_{h}= -{\omega }_{\mathrm{c}h}+ {jh}{\omega }_{1}$ in the $s$ domain. Secondly, when $z$ is located around ${e}^{{jh}{\omega }_{1}T}$, the circumference can be approximated by a straight line, and similarly, the angle $\angle \mathrm{{ZRP}}$ should be ${\varphi }_{h}$ (see Fig. 8 ), to make sure accurate phase lead. Thirdly, the ratio $\left|{z -{z}_{h}}\right|/\left|{z -{p}_{h}}\right|$ should be equal to ${K}_{1h}/{K}_{P}$, which determines a circular trajectory of ${z}_{h}$ centered at ${e}^{{jh}{\omega }_{1}T}$ with a radius of ${K}_{\mathrm{I}h}$ $\left({1 -{e}^{{jh}{\omega }_{1}T}}\right)/{K}_{\mathrm{P}}$ as shown in Fig. 8. Since ${z}_{h}$ is located on the circle and the half line at the same time, ${z}_{h}$ is obtained. The conjugate PRU can be obtained immediately, by $\left(z-\bar{z}_{h}\right) /\left(z-\bar{p}_{h}\right)$ if needed. Repeating the above steps gives other PRUs. And the final PR controller is realized by cascading all PRUs and ${K}_{\mathrm{P}}$ as (28).
Similarly, if the PRU is designed as a notch filter, ${K}_{1h}$ should be 0, and ${z}_{h}$ should be ${e}^{{jh}{\omega }_{1}T}$. Again, the distance between the pole and the zero should be restricted so as not to influence frequency characteristics in the distance. In $s$ domain, all poles and zeros are distributed along the imaginary axis, and the effect of zero’s movement can be separated into tangential and normal directions. Whereas in $z$ domain, such a situation no longer exists, tangential direction toward a pair of pole and zero may not be tangential anymore for another pole-zero pair. Therefore, the limit of zero is set instead as a circle centered at the pole with a radius of ${0.02}\min {\Delta \omega T}$ (see Fig. 8 ). With parameters in Table I, ${K}_{\mathrm{I}h}$ set to 100 and $h \in \{ 1,3,\ldots,{19}\}$, the maximal values for ${\omega }_{\mathrm{c}h}$ for all $h$ are obtained (the minimal values are zero). The minimum of these maximum values of ${\omega }_{\mathrm{c}h}$ is calculated to be 0.95 justifying setting ${\omega }_{\mathrm{c}h}$ to 1 .
Fig. 2 depicts the frequency responses of proposed cascade PR controllers realized in both $z$ domain and $s$ domain. As can be seen, they highly overlap indicating the $z$-domain realization method is able to accurately discretize the cascade PR controller, even without a preceding $s$-domain realization.
Fig. 2 displays the frequency responses of an illustrative PR controller implemented using both the parallel method and the proposed cascade method in the s domain. It can be seen that the resonant magnitudes of parallel $\mathrm{{PR}}$ controller exceed ${\mathrm{K}}_{\mathrm{I}h}$, in low frequency range. In addition, Fig. 2(b) is the zoom-in view of Fig. 2(a) around 19th order harmonic frequency. The parallel PR controller lags cascade PR controller by ${8.8}^{\circ }$ at 19th order harmonic frequency.
For a clearer comparison, Table III lists the expected resonant magnitudes $\left({K}_{\mathrm{I}h}\right)$ and phase leads $\left({\varphi }_{h}\right)$, and the actual value of the cascade PR controller at selected 1st, 3rd, 17th, and 19th order harmonics to see its characteristics at low frequency and high frequency. In contrast to parallel PR controller (see Table III) where magnitude error is as high as 15.2, the magnitude errors of cascade realization are at most -2.9. Regarding phase lead, the cascade realization manages to restrict phase error within $\pm {1}^{\circ }$ (see Table III) much smaller than the ${9.2}^{\circ }$ error of parallel PR controller (see Table III).
Phase lag in high frequency range is unfavorable to stability. This can be verified by Fig. 9 showcasing the Nyquist plots of open loop transfer functions, i.e.,
${G}_{\mathrm{o}}\left( s\right)= P\left( s\right){G}_{\mathrm{{PR}}}\left( s\right)$
with different PR controllers. In Fig. 9(a), the conventional parallel PR controller is used with no phase compensation. The system is stable with a highest order, $h$, of 11 (solid blue line), whereas it becomes unstable with $h$ reaching 13, when the Nyquist curve (red dashed line) encompasses the point,-1. For comparison, Fig. 9(b) shows the Nyquist plots of the open loop transfer function using parallel PR controller with compensation angle. The stability is significantly enhanced. The system is stable with PR controller tuned at up to 17th order harmonic, but unstable with 19th order harmonic. Further, when cascade PR controller with compensation angle is used (see Fig. 9(c)), the system can be stabilized with PR controller tuned at up to 19th order harmonic.
Table IV summarizes the highest order of harmonics different PR controllers can be tuned at while stabilizing the system of Fig. 1. Other parameters are listed in Table I. As can be seen, the highest order is related to two major factors, phase lead and integral gain. From a vertical view, the highest order the PR controller can be tuned at decreases as integral gain increases. From a horizontal view, a PR controller with phase lead can stabilize in a higher order than that without phase lead, justifying that phase lead can greatly increase stable margin. Moreover, cascade PR controller stabilizes the system with higher order harmonics than parallel PR controller because of more accurate phase lead.
To further verify the theoretical analysis, experiments are carried out on a single-phase grid-tied converter with an $\mathrm{L}$ filter, as shown in Fig. 10. The grid is emulated with a high-fidelity linear amplifier APS 15000. The applied half bridge module and the control platform are a PEB-SiC-8024 module and a B-BOX RCP control platform from Imperix, respectively. Sampling frequency and switching frequency are set to the same value of $5\mathrm{{kHz}}$. Other relevant system parameters of the single-phase grid-tied converter with L-Filter is presented in Table I.
To illustrate the harmonic resisting ability of multi-resonance PR controller, grid voltage is added with 10%3rd, 10% 9th and 10% 15th order harmonics. In experiment, DC link connects to only a capacitor without a DC voltage source. Thus, $d$-axis current reference is given by a DC voltage regulator to maintain a constant DC-link voltage. In steady state, the $d$-axis current is negligible. Meanwhile, the $q$-axis current reference, ${I}^{* }$, is set to constant ${10}\mathrm{\;A}.{K}_{1h}$ for $\mathrm{R}$ controllers is set to 180 for all $h$.
Fig. 11 shows the experimental result when a parallel PR controller tuned at only fundamental frequency is used. As can be seen, the output current is distorted heavily since no $\mathrm{R}$ controller tuned at harmonic frequencies is added. The FFT analysis conducted by the oscilloscope indicates that the magnitudes of 3rd, 9th, and 15th order harmonic components are about 1.6A, 2A and 2.8A in RMS, respectively.
Fig. 12 shows the experimental results with either a parallel PR controller (Fig. 12(a)) or a cascade PR controller (Fig. 12 (b)) implemented. Both PR controllers are tuned at all odd order harmonics up to 15th order. As it can be seen, both controllers can stabilize the system in consistency with Table IV. Besides, the output currents contain less harmonics compared with Fig. 11. Regarding their FFT analysis, the 3rd, 9th, and 15th order harmonic components are significantly reduced and can be neglected. This is because $\mathrm{R}$ controllers tuned at harmonics are added.
Further, Fig. 13 shows waveforms with the same setup as Fig. 12 except that an additional $\mathrm{R}$ controller tuned at 17th order harmonic is added to the PR controllers. It is obvious that the output current of the parallel PR controller (see Fig. 13(a)) contains an amount of 17th order harmonic. Whereas the output current of the cascade PR controller (see Fig. 13(b)) is rather sinusoidal with a negligible amount of harmonics. These are consistent with Table IV, that the parallel PR controller is unstable when tuned at 17th harmonic with a ${K}_{\mathrm{I}h}$ of 180, whereas the cascade PR controller is stable. It proves that the phase lag caused by inaccurate realization of parallel form jeopardizes the stability of the system.
This subsection investigates the PR controller’s harmonic-resisting ability with deviated fundamental frequency. The transfer function from grid voltage, ${v}_{\mathrm{{ac}}}$, to output current, ${i}_{\mathrm{o}}$, is derived as:
${G}_{\mathrm{v}}\left( s\right)= \frac{{i}_{\mathrm{o}}\left( s\right)}{{v}_{\mathrm{{ac}}}\left( s\right)} =- \frac{\frac{1}{{sL}+ R}}{1 + P\left( s\right){G}_{\mathrm{{PR}}}\left( s\right)} $
Fig. 14 presents the magnitude curves of (30) with cascade and parallel PR controllers. ${K}_{\mathrm{i}h}$ is set to 100 for all $h$. As can be seen, anomalous peaks appear around resonant frequencies [17], specifically 13th, 15th, and 17th order harmonics, indicating a larger current error will result from the corresponding frequency component of grid voltage. Nonetheless, the system with parallel PR controller always gives a larger anomalous peak compared with cascade PR controller, because of larger phase lag [17]. The comparison is stark from the zoom-in view around 17th harmonic. At 17.03 per unit, or 851.53Hz, the magnitudes of ${G}_{\mathrm{v}}\left( s\right)$ with parallel PR controller is 1.96 much larger than that of cascade PR controller, which is 0.48 .
To validate this comparison, relevant experiments were conducted. A 5% 17th harmonic was introduced into the grid voltage, with a fundamental frequency set to 50.09Hz. Consequently, the corresponding 17th harmonic frequency was 851.53 (17.03 per unit). Fig. 15 presents the experimental results for both cascade and parallel PR controllers. The current reference was set to zero to exclusively highlight the influence of the 17th order voltage harmonic. Upon conducting FFT analysis, it is evident that the harmonic current with the parallel PR controller is approximately 6A in RMS, almost three times the harmonic current observed with the cascade PR controller, which is approximately 2.2A in RMS. Note that the results may not be in strict quantitative alignment with Fig. 14, potentially due to challenges in precisely locating sharp anomalous peaks with the finite accuracy of the equipment used. Nevertheless, the outcome in Fig. 15 conclusively validates that the cascade PR controller exhibits smaller anomalous peaks compared to the parallel PR controller.
This paper proposes a cascade realization method for PR controllers to cope with the frequency characteristics error at resonant frequencies caused by small integral gain. Firstly, through pole-zero map analysis, the resonant decoupling property of the PR controller has been formulated, which enables PRUs to be designed or realized separately. Pole-zero placement method is used to realize PRU by mapping all parameters of the PR controller to the locations of the pole and zero of the PRU. The difference between a zero and its paired pole is restricted to restrict the frequency characteristics error. The proposed cascade realization is comprehensively compared with traditional parallel realization. It is found that cascade form realizes PR controller with much smaller magnitude and phase error than parallel form. Thus, cascade form can stabilize the system with $\mathrm{R}$ controller tuned at higher order compared to parallel form. Relevant experiments verify the correctness of the theoretical analysis and show the superiority of cascade realization.
  • National Natural Science Foundation of China(62203275)
  • Shandong Province Science Fund for Excellent Young Scientists Found Program(2022HYYQ-022)
  • National Key R&D Plan of China(2022YFF0712700)
[1]
W. Ding , L. Ming , C. Yin , M. Chen , P. C. Loh . "Differential-mode filters with common-mode neutral-point-balancing accelerators for sin-gle-phase symmetrical five-level converters". in IEEE Transactions on Power Electronics, 2021. 36 (8): 9209-9220.
[2]
D. N. Zmood D. G. Holmes . "Stationary frame current regulation of PWM inverters with zero steady-state error". in IEEE Transactions on Power Electronics, 2003. 18 (3): 814-822.
[3]
C. Gong , W. K. Sou , C. S. Lam . "Design and analysis of vector pro-portional-integral current controller for LC-coupling hybrid active power filter with minimum DC-link voltage". in IEEE Transactions on Power Electronics, 2021. 36 (8): 9041-9056.
[4]
W. Ding , L. Ming , C. Yin , P. C. Loh , B. Duan , C. Zhang . "An integrated common-mode fast-balancing mechanism for three-phase three-level converter with LCL filter". in IEEE Transactions on Power Electronics, 2021. 36 (11): 12694-12709.
[5]
Z. Li , Y. Li , P. Wang , H. Zhu , C. Liu , F. Gao . "Single-loop digital control of high-power 400-Hz ground power unit for airplanes". in IEEE Transactions on Industrial Electronics, 2010. 57 (2): 532-543.
[6]
M. Elkayam A. Kuperman . "Optimized design of multiresonant AC current regulators for single-phase grid-connected photovoltaic invert-ers". in IEEE Journal of Photovoltaics, 2019. 9 (6): 1815-1818.
[7]
L. Antonio De Souza Ribeiro , F. D. Freijedo , F. De Bosio , M. Soares Lima , J. M. Guerrero , M. Pastorelli . "Full discrete modeling, con-troller design, and sensitivity analysis for high-performance grid-forming converters in islanded microgrids". in IEEE Transactions on Industry Applications, 2018. 54 (6): 6267-6278.
[8]
R. Teodorescu , F. Blaabjerg , M. Liserre , P. C. Loh . "Proportion-al-resonant controllers and filters for grid-connected voltage-source con-verters". in IEE Proceedings: Electric Power Applications, 2006. 153 (5): 750-762.
[9]
Z. Xin , P. C. Loh , X. Wang , F. Blaabjerg , Y. Tang . "Highly accurate derivatives for LCL-filtered grid converter with capacitor voltage active damping". in IEEE Transactions on Power Electronics, 2016. 31 (5): 3612-3625.
[10]
R. I. Bojoi , G. Griva , V. Bostan , M. Guerriero , F. Farina , F. Profumo . "Current control strategy for power conditioners using sinusoidal signal integrators in synchronous reference frame". in IEEE Transactions on Power Electronics, 2005. 20 (6): 1402-1412.
[11]
B. Xie , K. Guo , M. Mao , L. Zhou , T. Liu , Q. Zhang , G. Hao . "Anal-ysis and improved design of phase compensated proportional resonant controllers for grid-connected inverters in weak grid". in IEEE Transac-tions on Energy Conversion, 2020. 35 (3): 1453-1464.
[12]
F. González-Espín , I. Patrao , E. Figueres , G. Garcerá . "An adaptive digital control technique for improved performance of grid connected inverters". in IEEE Transactions on Industrial Informatics, 2013. 9 (2): 708-718.
[13]
P. C. Loh , D. M. Vilathgamuwa , S. K. Tang , H. L. Long . "Multilevel dynamic voltage restorer". in IEEE Power Electronics Letters, 2004. 2 (4): 125-130.
[14]
D. N. Zmood , D. G. Holmes , G. H. Bode . "Frequency-domain anal-ysis of three-phase linear current regulators". in IEEE Transactions on Industry Applications, 2001. 37 (2): 601-610.
[15]
S. Zhou J. Liu . "Analysis and comparison of resonant-based cur-rent controllers implemented in stationary reference frame: A complex pole-zero placement perspective", in Proceedings of 2015 IEEE Energy Conversion Congress and Exposition (ECCE), Montreal, QC, Canada, 2015. 1624-1631.
[16]
M. Castilla , J. Miret , J. Matas , L. G. de Vicuña , J. M. Guerrero . "Control design guidelines for single-phase grid-connected photovol-taic inverters with damped resonant harmonic compensators". in IEEE Transactions on Industrial Electronics, 2009. 56 (11): 4492-4501.
[17]
A. G. Yepes , F. D. Freijedo , Ó. López , J. Doval-Gandoy . "Analysis and design of resonant current controllers for voltage-source converters by means of nyquist diagrams and sensitivity function". in IEEE Trans-actions on Industrial Electronics, 2011. 58 (11): 5231-5250.
[18]
A. G. Yepes , F. D. Freijedo , Ó. López , J. Doval-Gandoy . "High-per-formance digital resonant controllers implemented with two integrators". in IEEE Transactions on Power Electronics, 2011. 26 (2): 563-576.
[19]
C. Tang , K. Zhou , Y. Shu , Q. He , Q. Chen . "Analysis and design of multiple resonant current control for grid-connected converters". in IEEE Journal of Emerging and Selected Topics in Power Electronics, 2022. 10 (2): 2539-2546.
[20]
B. Lin , L. Peng , X. Liu . "Selective pole placement and cancellation for proportional-resonant control design used in voltage source inverter". in IEEE Transactions on Power Electronics, 2022. 37 (8): 8921-8934.
[21]
C. Lascu , L. Asiminoaei , I. Boldea , F. Blaabjerg . "Frequency response analysis of current controllers for selective harmonic com-pensation in active power filters". in IEEE Transactions on Industrial Electronics, 2009. 56 (2): 337-347.
[22]
C. Lascu , L. Asiminoaei , I. Boldea , F. Blaabjerg . "High performance current controller for selective harmonic compensation in active power filters". in IEEE Transactions on Power Electronics, 2007. 22 (5): 1826-1835.
[23]
H. Bai , X. Wang , P. C. Loh , F. Blaabjerg . "An active trap filter for switching harmonic attenuation of low-pulse-ratio inverters". in IEEE Transactions on Power Electronics, 2017. 32 (12): 9078-9092.
[24]
C. Xia , B. Ji , Y. Yan . "Smooth speed control for low-speed high-torque permanent-magnet synchronous motor using proportional-inte-gral-resonant controller". in IEEE Transactions on Industrial Electronics, 2015. 62 (4): 2123-2134.
[25]
C. Xie , X. Zhao , K. Li , J. Zou , J. M. Guerrero . "A new tuning meth-od of multiresonant current controllers for grid-connected voltage source converters". in IEEE Journal of Emerging and Selected Topics in Power Electronics, 2019. 7 (1): 458-466.
[26]
F. Hans , W. Schumacher , S. F. Chou , X. Wang . "Design of multi-frequency proportional-resonant current controllers for voltage-source converters". in IEEE Transactions on Power Electronics, 2020. 35 (12): 13573-13589.
[27]
C. A. Busada , S. Gómez Jorge , A. E. Leon , J. A. Solsona . "Current controller based on reduced order generalized integrators for distributed generation systems". in IEEE Transactions on Industrial Electronics, 2012. 59 (7): 2898-2909.
Year 2024 volume 9 Issue 1
PDF
260
152
Cite this Article
BibTeX
Article Info
doi: 10.24295/CPSSTPEA.2023.00046
  • Receive Date:2023-06-06
  • Online Date:2025-07-05
  • Published:2024-03-10
Article Data
Affiliations
History
  • Received:2023-06-06
  • Revised:2023-08-30
  • Accepted:2023-10-25
Funding
National Natural Science Foundation of China(62203275)
Shandong Province Science Fund for Excellent Young Scientists Found Program(2022HYYQ-022)
National Key R&D Plan of China(2022YFF0712700)
Affiliations
    1 The Chinese University of Hong Kong Department of Electronic Engineering Hong Kong China
    2 Stanford University School of Engineering Stanford 94305 United States
    3 Fudan University Academy for Engineering and Technology Shanghai 200433 China
    4 Shandong University School of Control Science and Engineering Jinan 250100 China

Corresponding:

Wenlong Ding.
References
Share
https://castjournals.cast.org.cn/joweb/dldzjs/EN/10.24295/CPSSTPEA.2023.00046
Share to
QR

Scan QR to access full text

Cite this article
BibTeX
Citations
表12种不同金属材料的力学参数

Family
属数
Number of
genus
种数
Number of
species
占总种数比例
Percentage of
total species (%)

Genus
种数
Number of
species
占总种数比例
Percentage of total
species (%)
鹅膏菌科Amanitaceae 2 11 5.26 鹅膏菌属 Amanita 10 4.78
小菇科 Mycenaceae 2 12 5.74 丝盖伞属 Inocybe 5 2.39
多孔菌科 Polyporaceae 8 14 6.70 蜡蘑属 Laccaria 5 2.39
红菇科 Russulaceae 3 23 11.00 小皮伞属 Marasmius 6 2.87
小菇属 Mycena 11 5.26
光柄菇属 Pluteus 5 2.39
红菇属 Russula 17 8.13
栓菌属 Trametes 5 2.39
关闭全屏
  • BibTeX
  • EndNote
  • RefWorks
  • TxT