PREDICTIVE CONTROL METHOD AND SYSTEM OF GRID-CONNECTED CONVERTER BASED ON STRUCTURALLY ADAPTIVE EXTENDED STATE OBSERVER

Abstract

A predictive control method of a grid-connected converter based on a structurally adaptive extended state observer (ESO) is provided. The method includes: obtaining data of a grid current and a grid voltage and converting the data to data in a dq coordinate system; calculating a selected voltage vector based on a voltage in the dq coordinate system; adaptively inputting the current in the dq coordinate system to a parallel ESO, a cascade ESO, or a hybrid cascade-parallel ESO, so as to obtain a current predicted current value and estimated total disturbance; carrying out a two-step grid current prediction based on the current predicted current value, the estimated total disturbance and the selected voltage vector; and taking minimizing a cost function as a control target, so as to obtain an optimal voltage vector based on a two-step grid current prediction result, and carrying out switching control of the grid-connected converter.

Claims

1. A predictive control method of a grid-connected converter based on a structurally adaptive extended state observer (ESO), comprising: obtaining data of a grid current and a grid voltage and converting the data to data in a dq coordinate system; calculating a selected voltage vector based on a voltage in the dq coordinate system; adaptively inputting a current in the dq coordinate system to a parallel ESO, a cascade ESO, or a hybrid cascade-parallel ESO, so as to obtain a current predicted current value and estimated total disturbance; carrying out a two-step grid current prediction based on the current predicted current value, the estimated total disturbance and the selected voltage vector; taking minimizing a cost function as a control target, so as to obtain an optimal voltage vector based on a two-step grid current prediction result, and carrying out switching control of the grid-connected converter; wherein adaptively inputting the current in the dq coordinate system to a parallel ESO, a cascade ESO, or a hybrid cascade-parallel ESO, so as to obtain a current predicted current value and estimated total disturbance includes: determining whether a d-axis current component in the dq coordinate system is a starting current or a disturbance current; if yes, firstly predicting the current predicted current value and the estimated total disturbance by using the parallel ESO according to the current in the dq coordinate system; after cycle, predicting the current predicted current value and the estimated total disturbance by using a hybrid cascade-parallel ESO with a first set bandwidth; after another cycle, predicting the current predicted current value and the estimated total disturbance by using a hybrid cascade-parallel ESO with a second set bandwidth; after another cycle, predicting the current predicted current value and the estimated total disturbance by using the cascade ESO; after setting a delay for a plurality of cycles, returning to obtaining the d-axis current component in the dq coordinate system again; if not, firstly predicting the current predicted current value and the estimated total disturbance by directly using the cascade ESO according to the current in the dq coordinate system; wherein carrying out a two-step grid current prediction based on the current predicted current value, the estimated total disturbance and the selected voltage vector includes: the current predicted current value is a result obtained after adaptively inputting to the parallel ESO, the cascade ESO or the hybrid cascade-parallel ESO, that is, a result obtained in a first prediction step; a result of a second current prediction step is as follows: ${\hat{i}}_{\mathrm{dq}}\left(k+2\right)={\hat{i}}_{\mathrm{dq}}\left(k+1\right)+T_S\left[\underset{j=1}{\overset{3}{.Math.}}{\hat{F}}_j\left(k+1\right)+u\left(k\right)\right]$ wherein k is sampling instant, T.sub.s is sampling time, {circumflex over (l)}.sub.dq(k+1) is the result obtained in the first prediction step and is a predicted estimated value of the grid current at the next sampling instant (k+1), is a constant control input gain, u(k) is a voltage of the converter caused by a switching state Sabc(k); {circumflex over (F)}.sub.j(k+1) represents a predicted disturbance value of the first step at a jth sub-frequency level .sub.0j, wherein j[1,2,3]; wherein taking minimizing a cost function as the control target, so as to obtain the optimal voltage vector, the cost function includes: $J={\left[{\hat{i}}_{\mathrm{dq}}\left(k+2\right)-i_{\mathrm{dq}}^{*}\left(k+2\right)\right]}^2$ wherein {circumflex over (l)}.sub.dq(k+2) is a predicted estimated value of the grid current at sampling instant (k+2), and i.sub.dq.sup.*(k+2) is a reference value of the grid current at the sampling instant (k+2).

2. The predictive control method of the grid-connected converter based on the structurally adaptive ESO according to claim 1, wherein determining whether a d-axis current component in the dq coordinate system is a starting current or a disturbance current includes: when the following formula is satisfied, the starting current or disturbance current appears: $i_{d\_4}=.Math.{}_0^{t\_4}\frac{d}{\mathrm{dt}}{i}_{d}d.Math.0.01i_{d\_\max }$ wherein i.sub.d_4 is an absolute change value of the d-axis current i.sub.d after 4 times of sampling of the d-axis current, t_4=4/T.sub.s and, T.sub.s is the sampling time, and i.sub.d_max is a maximum d-axis grid current.

3. The predictive control method of the grid-connected converter based on the structurally adaptive ESO according to claim 1, wherein the hybrid cascade-parallel ESO is as follows: $\stackrel{.}{{\hat{z}}_1}\left(t\right)={\hat{F}}_1\left(t\right)+u\left(t\right)-{}_{11}\left({\hat{z}}_1\left(t\right)-y\left(t\right)\right)\stackrel{.}{{\hat{F}}_1}\left(t\right)=-{}_{21}\left({\hat{z}}_1\left(t\right)-y\left(t\right)\right)\stackrel{.}{{\hat{z}}_2}\left(t\right)=\underset{j=1}{\overset{2}{.Math.}}{\hat{F}}_j\left(t\right)+u\left(t\right)-{}_{12}\left[{\hat{z}}_2\left(t\right)-{\hat{z}}_1\left(t\right)\right]\stackrel{.}{{\hat{F}}_2}\left(t\right)=-{}_{22}\left[{\hat{z}}_2\left(t\right)-{\hat{z}}_1\left(t\right)\right]\stackrel{.}{{\hat{z}}_3}\left(t\right)=\underset{j=1}{\overset{3}{.Math.}}{\hat{F}}_j\left(t\right)+u\left(t\right)-{}_{13}\left[{\hat{z}}_3\left(t\right)-{\hat{z}}_1\left(t\right)\right]\stackrel{.}{{\hat{F}}_3}=-{}_{23}\left[{\hat{z}}_3\left(t\right)-{\hat{z}}_1\left(t\right)\right]$ wherein {circumflex over ()}.sub.j(t) is a first order derivative of time {circumflex over (z)}.sub.j(t), {circumflex over (z)}.sub.j(t) is an ESO state variable of an estimated current predicted current {circumflex over (l)}.sub.dq,j(t), y(t) is a measured noisy output signal, {.sub.1j, .sub.2j} is a gain of an ESO of ESO.sub.k, wherein j[1,2,3]; .sub.0= is a total bandwidth of an ESO system, is a gain of a constant control input; u(t) is a controller input.

4. The predictive control method of the grid-connected converter based on the structurally adaptive ESO according to claim 1, wherein each switching state voltage of a three-phase two-level grid-connected converter is evaluated in the cost function; and a switching state corresponding to a voltage u.sub.dq with the minimum value of the cost function is applied as the switching state Sabc of the grid-connected converter.

5. A predictive control system of a grid-connected converter based on a structurally adaptive ESO, comprising: a data obtaining module, configured for obtaining data of a grid current and a grid voltage and converting the data to data in a dq coordinate system; a voltage vector selecting module, configured for calculating a selected voltage vector based on a voltage in the dq coordinate system; an adaptive state observing module, configured for adaptively inputting a current in the dq coordinate system to a parallel ESO, a cascade ESO, or a hybrid cascade-parallel ESO, so as to obtain a current predicted current value and estimated total disturbance; a two-step predicting module, configured for carrying out a two-step grid current prediction based on the current predicted current value, the estimated total disturbance and the selected voltage vector; a grid-connected converter controlling module, configured for taking minimizing a cost function as a control target, so as to obtain an optimal voltage vector based on a two-step grid current prediction result, and carrying out switching control of the grid-connected converter; wherein adaptively inputting the current in the dq coordinate system to a parallel ESO, a cascade ESO, or a hybrid cascade-parallel ESO, so as to obtain a current predicted current value and estimated total disturbance includes: determining whether a d-axis current component in the dq coordinate system is a starting current or a disturbance current; if yes, firstly predicting the current predicted current value and the estimated total disturbance by using the parallel ESO according to the current in the dq coordinate system; after cycle, predicting the current predicted current value and the estimated total disturbance by using a hybrid cascade-parallel ESO with a first set bandwidth; after another cycle, predicting the current predicted current value and the estimated total disturbance by using a hybrid cascade-parallel ESO with a second set bandwidth; after another cycle, predicting the current predicted current value and the estimated total disturbance by using the cascade ESO; after setting a delay for a plurality of cycles, returning to obtaining the d-axis current component in the dq coordinate system again; if not, firstly predicting the current predicted current value and the estimated total disturbance by directly using the cascade ESO according to the current in the dq coordinate system; wherein carrying out a two-step grid current prediction based on the current predicted current value, the estimated total disturbance and the selected voltage vector includes: the current predicted current value is a result obtained after adaptively inputting to the parallel ESO, the cascade ESO or the hybrid cascade-parallel ESO, that is, a result obtained in a first prediction step; a result of a second current prediction step is as follows: ${\hat{i}}_{\mathrm{dq}}\left(k+2\right)={\hat{i}}_{\mathrm{dq}}\left(k+1\right)+T_S\left[\underset{j=1}{\overset{3}{.Math.}}{\hat{F}}_j\left(k+1\right)+u\left(k\right)\right]$ wherein k is sampling instant, T.sub.s is sampling time, {circumflex over (l)}.sub.dq(k+1) is the result obtained in the first prediction step and is a predicted estimated value of the grid current at the next sampling instant (k+1), is a constant control input gain, u(k) is a voltage of the converter caused by a switching state Sabc(k); {circumflex over (F)}.sub.j(k+1) represents a predicted disturbance value of the first step at a jth sub-frequency level .sub.0j, wherein j[1,2,3]; wherein taking minimizing a cost function as the control target, so as to obtain the optimal voltage vector, the cost function includes: $J={\left[{\hat{i}}_{\mathrm{dq}}\left(k+2\right)-i_{\mathrm{dq}}^{*}\left(k+2\right)\right]}^2$ wherein {circumflex over (l)}.sub.dq(k+2) is a predicted estimated value of the grid current at sampling instant (k+2), and i.sub.dq.sup.*(k+2) is a reference value of the grid current at the sampling instant (k+2).

6. A terminal device, comprising a processor and a memory, wherein the processor is used to implement a plurality of instructions; and the memory is used to store the plurality of instructions, the plurality of instructions are suitable for being loaded by the processor and executing the predictive control method of the grid-connected converter based on the structurally adaptive ESO according to claim 1.

7. A computer-readable storage medium in which a plurality of instructions are stored, the plurality of instructions are suitable for being loaded by a processor of a terminal device and executing the predictive control method of the grid-connected converter based on the structurally adaptive ESO according to claim 1.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0029] FIG. 1 is a structurally schematic diagram of a three-phase grid-connected converter in an example of the present invention.

[0030] FIG. 2(a)FIG. 2(c) are respectively schematic diagrams of ESOs in parallel, cascade and hybrid cascade-parallel in examples of the present invention.

[0031] FIG. 3 is a flow chart of adaptive ESO structure transition in an example of the present invention.

[0032] FIG. 4 is a schematic diagram of a predictive control method of the grid-connected converter based on a structurally adaptive ESO in an example of the present invention.

[0033] FIG. 5 shows comparative performance of an adaptive ESO observer and a CESO observer during startup of a power converter in the example.

DESCRIPTION OF THE EMBODIMENTS

[0034] It should be noted that the following detailed explanations are illustrative and are intended to provide further explanations of this application. Unless otherwise indicated, all technical and scientific terms used in the present invention have the same meaning as that would normally be understood by a person of ordinary skill in the technique field to which this application belongs.

[0035] It should be noted that the terms used herein are intended only to describe specific embodiments and are not intended to limit exemplary embodiments under this application. As used herein, a singular form is also intended to include a plural form unless the context expressly indicates otherwise, and it should also be understood that when terms contain and/or include are used in this specification, they indicate the presence of features, steps, operations, devices, components and/or combinations thereof.

Example 1

[0036] In one or more embodiments, disclosed in combination with FIG. 4 is a predictive control method of a grid-connected converter based on a structurally adaptive ESO, specifically comprising the following process:

[0037] Step (1): obtaining data of a grid current and a grid voltage and converting the data to data in a dq coordinate system;

[0038] FIG. 1 shows a structurally schematic diagram of a three-phase grid-connected converter. A schematic diagram of a grid-connected power converter has the following dynamic model:

[00005]$\begin{array}{cc}L\frac{{\mathrm{di}}_{\mathrm{abc}}}{\mathrm{dt}}={\mathrm{Ri}}_{\mathrm{abc}}-e_{\mathrm{gabc}}+u_{\mathrm{abc}}& \left(1\right)\end{array}$

Wherein i.sub.abc is the grid current, e.sub.gabc is the grid voltage; u.sub.abc=(S.sub.abc) represents an output voltage of the power converter and is a function of a switching state S.sub.abc; L represents filter inductance, R represents filter resistance.

[0039] Formula (1) can be converted into:

[00006]$\begin{array}{cc}\frac{{\mathrm{di}}_{\mathrm{abc}}}{\mathrm{dt}}=F\left(t\right)+u_{\mathrm{abc}}\left(t\right)& \left(2\right)\end{array}$

Wherein u.sub.abc is a control input of the converter,

[00007]$=\frac{1}{L}$

[0040] is a constant control input gain,

[00008]$F=\frac{R}{L}\left(i_{\mathrm{abc}}-e_{\mathrm{gabc}}\right).$

[0041] In this example, obtaining the data of the grid current i.sub.abc and the grid voltage e.sub.gabc and converting the data to data in the dq coordinate system, so as to obtain a current i.sub.dq.sup.meas and a voltage e.sub.dq in the dq coordinate system.

[0042] Step (2): calculating, by using a phase locked loop (PLL) module, a voltage angle of the grid based on a voltage in the dq coordinate system. A voltage vector u.sub.dq is selected by applying the voltage phase angle of the grid.

[0043] Step (3): adaptively inputting the current in the dq coordinate system to a parallel ESO, a cascade ESO, or a hybrid cascade-parallel ESO, so as to obtain a current predicted current value .sub.dq(k+1) and estimated total disturbance {circumflex over (F)}.sub.dq;

[0044] In this example, the total disturbance F is adaptively estimated using the following three ESOs: the parallel ESO (PESO), the cascade ESO (CESO), and the hybrid cascade-parallel ESO (CP-ESO).

[0045] The three ESOs are described below:

[0046] (3-1) the parallel ESO (PESO)

[0047] As shown in FIG. 2(a), the PESO is described in a time domain as follows:

[00009]$\begin{array}{cc}\{\begin{array}{c}\stackrel{.}{{\hat{z}}_j}\left(t\right)={\hat{F}}_j\left(t\right)+u\left(t\right)-{}_{1j}\left[{\hat{z}}_j\left(t\right)-y\left(t\right)\right]\\ \stackrel{.}{{\hat{F}}_j}\left(t\right)=-{}_{2j}\left[{\hat{z}}_j\left(t\right)-y\left(t\right)\right]\end{array}& \left(3\right)\end{array}$

Wherein {circumflex over (z)}.sub.j(t).fwdarw..sub.dq,j(t) is an estimated state, y(t).fwdarw.i.sub.dq.sup.meas(t) is a measured output signal of noise, and .sub.1j, .sub.2j}j[1,2,3] is a gain of the observer of ESO.sub.j.square-solid.. The gain of the observer is further defined as follows:

[00010]${}_{11}=2{}_{01},{}_{21}={}_{01}^2,{}_{12}=2{}_{02},{}_{22}={}_{02}^2,{}_{13}=2{}_{03},{}_{23}={}_{03}^2;{}_{01}=\frac{{}_0}{M^2},{}_{02}=\frac{{}_0}{M},{}_{03}={}_0;{}_0$

is a total bandwidth of the ESO system.

[0048] It should be noted that z.sub.j.square-solid..square-solid.j{1,2,3} is an ESO state variable for estimating the current .sub.dq,j.

[0049] In addition, three sub-frequencies satisfy: .sub.01<.sub.02<.sub.03=.sub.0.

[0050] The PESO is implemented on a digital signal processor (or a microcontroller) as follows:

[00011]$\begin{array}{cc}\{\begin{array}{c}[{\hat{z}}_j\left(k+1\right)=[{\hat{z}}_j\left(k\right)+T_s\left[{\hat{F}}_j\left(k\right)+u\left(k\right)\right]-T_s{}_{1j}\left[{\hat{z}}_j\left(k\right)-y\left(k\right)\right]\\ {\hat{F}}_j\left(k+1\right)=-T_s{}_{2j}\left[{\hat{z}}_j\left(k\right)-y\left(k\right)\right]\end{array}& \left(4\right)\end{array}$

Wherein k is a discrete moment, T.sub.s is sampling time, {circumflex over (z)}.sub.j(tk).fwdarw..sub.dq,j(k) is an estimated state, y(k).fwdarw.i.sub.dq.sup.meas(k) is the measured output signal of the noise, and all other variables are described as above.

[0051] A corresponding current prediction in a first step:

[00012]$\begin{array}{cc}{\hat{i}}_{\mathrm{dq}}\left(k+1\right)={\hat{i}}_{\mathrm{dq}}\left(k\right)+T_S\left[{.Math.}_{j=1}^3{\hat{F}}_j\left(k\right)+u\left(k\right)\right]-T_s\left({}_{11}+{}_{12}+{}_{13}\right)\left[{\hat{i}}_{\mathrm{dq}}\left(k\right)-i_{\mathrm{dq}}\left(k\right)\right]& \left(5\right)\end{array}$

[0052] (3-2) the cascade ESO (CESO)

[0053] As shown in FIG. 2(b), the CESO is described in the time domain as follows:

[00013]$\begin{array}{cc}\{\begin{array}{c}\stackrel{.}{{\hat{z}}_1}\left(t\right)={\hat{F}}_1\left(t\right)+u\left(t\right)-{}_{11}\left[{\hat{z}}_1\left(t\right)-y\left(t\right)\right]\\ \stackrel{.}{{\hat{F}}_1}\left(t\right)=-{}_{21}\left[{\hat{z}}_1\left(t\right)-y\left(t\right)\right]\\ \stackrel{.}{{\hat{z}}_j}\left(t\right)={.Math.}_1^j{\hat{F}}_j\left(t\right)+u\left(t\right)-{}_{1j}\left[{\hat{z}}_j\left(t\right)-{\hat{z}}_{j-1}\left(t\right)\right]j\left\{2,3\right\}\\ \stackrel{.}{{\hat{F}}_j}\left(t\right)=-{}_{2j}\left[{\hat{z}}_j\left(t\right)-{\hat{z}}_{j-1}\left(t\right)\right]\end{array}& \left(6\right)\end{array}$

Wherein {circumflex over (z)}.sub.j(t) is an ESO state variable of an estimated current predicted current .sub.dq,j(t), {circumflex over ()}.sub.j(t) is a first order derivative with respect to {circumflex over (z)}.sub.j(t) time, and u(t) is a controller input. y(t).fwdarw.i.sub.dq.sup.meas(t) is the measured output signal of noise, and .sub.1j, .sub.2j}j[1,2,3] is the gain of the observer of ESO.sub.j.square-solid..

[0054] The gain of the observer is further defined as

[00014]${}_{11}=2{}_{01},{}_{21}={}_{01}^2,{}_{12}=2{}_{02},{}_{22}={}_{02}^2,{}_{13}=2{}_{03},{}_{23}={}_{03}^2;{}_{01}=\frac{{}_0}{M^2},{}_{02}=\frac{{}_0}{M},{}_{03}={}_0;{}_0$

is a total bandwidth of an ESO system.

[0055] It should be noted that z.sub.j.square-solid..square-solid.j{1,2,3} is an ESO state variable for estimating the current .sub.dq,j.

[0056] In addition, three sub-frequencies satisfy: .sub.01<.sub.02<.sub.03=.sub.0.

[0057] The CESO is implemented on the digital signal processor (or the microcontroller) as follows:

[00015]$\begin{array}{cc}{\hat{z}}_1\left(k+1\right)={\hat{F}}_1\left(k\right)+u\left(k\right)-{}_{11}\left[{\hat{z}}_1\left(k\right)-y\left(k\right)\right]{\hat{F}}_1\left(k+1\right)=-{}_{21}\left[{\hat{z}}_1\left(k\right)-y\left(k\right)\right]{\hat{z}}_j\left(k+1\right)={.Math.}_1^j{\hat{F}}_j\left(k\right)+u\left(k\right)-{}_{1j}\left[{\hat{z}}_j\left(k\right)-{\hat{z}}_{j-1}\left(k\right)\right]j\left\{2,3\right\}{\hat{F}}_j\left(k+1\right)=-T_s{}_{2j}\left[{\hat{z}}_j\left(k\right)-{\hat{z}}_{j-1}\left(k\right)\right]& \left(7\right)\end{array}$

Wherein k is the discrete moment, T.sub.s is the sampling time, {circumflex over (z)}.sub.j(t) is the ESO state variable of an estimated current predicted current .sub.dq,j(t), y(k).fwdarw.i.sub.dq.sup.meas(k) is the measured output signal of the noise, and all other variables are described as above. {circumflex over (F)}.sub.j(k+1) is a predicted disturbance value of a jth sub-frequency level .sub.0j in a first step.

[0058] A corresponding current prediction in the first step is as follows:

[00016]$\begin{array}{cc}{\hat{i}}_{\mathrm{dq}}\left(k+1\right)={\hat{i}}_{\mathrm{dq}}\left(k\right)+T_S\left[{.Math.}_{j=1}^3{\hat{F}}_j\left(k\right)+u\left(k\right)\right]-T_s{}_{13}\left[{\hat{i}}_{\mathrm{dq},3}\left(k\right)-{\hat{i}}_{\mathrm{dq},2}\left(k\right)\right]& \left(8\right)\end{array}$

[0059] (3-3) ESO in cascade parallel (CP-ESO)

[0060] As shown in FIG. 2(c), the CP-ESO is described in the time domain as follows:

[00017]$\begin{array}{cc}\{\begin{array}{c}\stackrel{.}{{\hat{z}}_1}\left(t\right)={\hat{F}}_1\left(t\right)+u\left(t\right)-{}_{11}\left({\hat{z}}_1\left(t\right)-y\left(t\right)\right)\\ \stackrel{.}{{\hat{F}}_1}\left(t\right)=-{}_{21}\left({\hat{z}}_1\left(t\right)-y\left(t\right)\right)\\ \stackrel{.}{{\hat{z}}_2}\left(t\right)={.Math.}_{j=1}^2{\hat{F}}_j\left(t\right)+u\left(t\right)-{}_{12}\left[{\hat{z}}_2\left(t\right)-{\hat{z}}_1\left(t\right)\right]\\ \stackrel{.}{{\hat{F}}_2}\left(t\right)=-{}_{22}\left[{\hat{z}}_2\left(t\right)-{\hat{z}}_1\left(t\right)\right]\\ \stackrel{.}{{\hat{z}}_3}\left(t\right)={.Math.}_{j=1}^3{\hat{F}}_j\left(t\right)+u\left(t\right)-{}_{13}\left[{\hat{z}}_3\left(t\right)-{\hat{z}}_1\left(t\right)\right]\\ \stackrel{.}{{\hat{F}}_3}=-{}_{23}\left[{\hat{z}}_3\left(t\right)-{\hat{z}}_1\left(t\right)\right]\end{array}& \left(9\right)\end{array}$

Wherein {circumflex over (z)}.sub.j(t) is an ESO state variable of an estimated current predicted current .sub.dq,j(t), y(t) is the measured output signal of noise, and .sub.1j, .sub.2j}j[1,2,3] is the gain of the observer of ESO.sub.j. The gain of the observer is further defined as follows:

[00018]${}_{11}=2{}_{01},{}_{21}={}_{01}^2,{}_{12}=2{}_{02},{}_{22}={}_{02}^2,{}_{13}=2{}_{03},{}_{23}={}_{03}^2;{}_{01}=\frac{{}_0}{M^2},{}_{02}=\frac{{}_0}{M},{}_{03}={}_0;$

.sub.0 is a total bandwidth of the ESO system.

[0061] It should be noted that z.sub.j.square-solid..square-solid.j{1,2,3} is an ESO state variable for estimating the current .sub.dq,j.

[0062] In addition, three sub-frequencies satisfy: .sub.01<.sub.02<.sub.03=.sub.0.

[0063] The CP-ESO is implemented on the digital signal processor (or the microcontroller) as follows:

[00019]$\begin{array}{cc}\{\begin{array}{c}{\hat{z}}_1\left(k+1\right)={\hat{z}}_1\left(k\right)+T_s\left\{{\hat{F}}_1\left(k\right)+u\left(k\right)-{}_{11}\left({\hat{z}}_1\left(k\right)-y\left(k\right)\right)\right\}\\ {\hat{F}}_1\left(k+1\right)={\hat{F}}_1\left(k\right)-T_s{}_{21}\left({\hat{z}}_1\left(t\right)-y\left(t\right)\right)\\ {\hat{z}}_2\left(k+1\right)={\hat{z}}_2\left(k\right)+T_s{.Math.}_{j=1}^2{\hat{F}}_j\left(k\right)+T_{s}u\left(t\right)-T_s{}_{12}\left[{\hat{z}}_2\left(k\right)-{\hat{z}}_1\left(k\right)\right]\\ {\hat{F}}_2\left(k+1\right)={\hat{F}}_2\left(k\right)-T_s{}_{22}\left[{\hat{z}}_2\left(k\right)-{\hat{z}}_1\left(k\right)\right]\\ {\hat{z}}_3\left(k+1\right)={\hat{z}}_3\left(k\right)+T_s{.Math.}_{j=1}^3{\hat{F}}_j\left(k\right)+T_{s}u\left(k\right)-T_s{}_{13}\left[{\hat{z}}_3\left(k\right)-{\hat{z}}_1\left(k\right)\right]\\ {\hat{F}}_3\left(k+1\right)={\hat{F}}_3\left(k\right)-T_s{}_{23}\left[{\hat{z}}_3\left(k\right)-{\hat{z}}_1\left(k\right)\right]\end{array}& \left(10\right)\end{array}$

Wherein k is the discrete moment, T.sub.s is the sampling time, {circumflex over (z)}.sub.j(tk).fwdarw..sub.dq,j(k) is the estimated state, y(k).fwdarw.i.sub.dq.sup.meas(k) is the measured output signal of the noise, and all other variables are described as above.

[0064] A corresponding current prediction in the first step:

[00020]$\begin{array}{cc}{\hat{i}}_{\mathrm{dq}}\left(k+1\right)={\hat{i}}_{\mathrm{dq}}\left(k\right)+T_S\left[{.Math.}_{j=1}^3{\hat{F}}_j\left(k\right)+u\left(k\right)\right]-T_s\left({}_{12}+{}_{13}\right)\left[{\hat{i}}_{\mathrm{dq},2}\left(k\right)-{\hat{i}}_{\mathrm{dq},1}\left(k\right)\right]& \left(11\right)\end{array}$

[0065] PESO has the best disturbance rejection capability. Therefore, it is most useful during device startup and when external disturbance is detected. It may minimize a peak current during transients. The CESO has the best noise suppression characteristics and would be applied during steady-state operations of the system. The CP-ESO would only be applied during transition from the PESO to the CESO.

[0066] Combined with FIG. 3, this example adaptively selects the observer by measuring a d-axis current i.sub.d, and a specific process is as follows:

[0067] Determining whether a d-axis current component in the dq coordinate system is a starting current or a disturbance current; if yes, firstly predicting the current predicted current value and the estimated total disturbance by using the parallel ESO according to the current in the dq coordinate system; after cycle, predicting the current predicted current value and the estimated total disturbance by using a hybrid cascade-parallel ESO (a bandwidth is 5.sub.0); after cycle, predicting the current predicted current value and the estimated total disturbance by using a hybrid cascade-parallel ESO (a bandwidth is 3.sub.0); after cycle, predicting the current predicted current value and the estimated total disturbance by using the cascade ESO; after preset cycles (such as 3 cycles), returning to obtaining the d-axis current component in the dq coordinate system again; if not, firstly predicting the current predicted current value and the estimated total disturbance by directly using the cascade ESO according to the current in the dq coordinate system.

[0068] In this example, when the following formula is true, the starting current or disturbance current appears:

[00021]$i_{d\_4}=.Math.{}_0^{t\_4}\frac{d}{\mathrm{dt}}{i}_{d}d.Math.0.01i_{d\_\max }$

Wherein i.sub.d_4 is an absolute change value of the d-axis gate current i.sub.d after 4 times of sampling of the d-axis gate current, t_4=4/T.sub.s; T.sub.s is the sampling time and i.sub.d_max is a maximum d-axis grid current.

[0069] For a microcontroller implementation, it is discretized to:

[00022]$i_{d\_4}=.Math.\underset{k=1}{\overset{4}{.Math.}}\left[i_d\left(k\right)-i_d\left(k-1\right)\right].Math.0.01i_{d\_\max }$

Wherein i.sub.d_4 is the absolute change of the current after four times of sampling, k is the discrete sampling instant, i.sub.d_max is the maximum d-axis grid current.

[0070] Step (4): carrying out a two-step grid current prediction based on the current predicted current value, the estimated total disturbance and the selected voltage vector;

[0071] In this example, a predicted current {circumflex over (l)}.sub.dq(k+2) in a second step is calculated using {circumflex over (l)}.sub.dq(k+1), {circumflex over (F)}.sub.dqand u.sub.dq for equation (6). The predicted current {circumflex over (l)}.sub.dq(k+2) and a reference current calculated using a PI controller are passed to a cost function stage.

[0072] In combination with FIG. 3, the predicted current at the current timing in the first step is obtained using step (3), which is obtained after discretizing the ESO (for PESO (5) or CESO (8) or CP-ESO (11)).

[0073] The currently active ESO (for PESO (5) or CESO (8) or CP-ESO (11)), {circumflex over (l)}.sub.dq(k+1) and .sub.j=1.sup.3{circumflex over (F)}.sub.j(k+1) are provided to calculate the prediction {circumflex over (l)}.sub.dq(k+2) in the second step in equation (12).

[0074] The predicted current in the second step is as follows:

[00023]$\begin{array}{cc}{\hat{i}}_{\mathrm{dq}}\left(k+2\right)={\hat{i}}_{\mathrm{dq}}\left(k+1\right)+T_S\left[{.Math.}_{j=1}^3{\hat{F}}_j\left(k+1\right)+u\left(k\right)\right]& \left(12\right)\end{array}$

Wherein k is the sampling instant, T.sub.s is the sampling time,

[00024]${}_{12}=2{}_{02},{}_{13}=2{}_{03},{}_{01}=\frac{{}_0}{M^2},{}_{02}=\frac{{}_0}{M},{}_{03}={}_0;$

.sub.0 is the bandwidth system of the entire ESO, {circumflex over (l)}.sub.dq(k+1) is the predicted estimated value of the grid current at the next sampling instant (k+1), {circumflex over (l)}.sub.dq(k) is the estimated current of the current discrete sample (k), T.sub.s is the sampling time, =1/L, u(k) is selected from Table 1.

[0075] A voltage S.sub.abc(k) of the converter caused by the switching state of i.sub.dq(k) is a measured current of the current sampling time and {circumflex over (F)}(k) is the estimated disturbance in the ESO bandwidth.

[0076] Step (5): taking minimizing a cost function as a control target, so as to obtain an optimal voltage vector based on a two-step grid current prediction result, and carrying out switching control of the grid-connected converter.

[0077] In this example, minimizing the cost function is specifically:

[00025]$\begin{array}{cc}J={\left[{\hat{i}}_{\mathrm{dq}}\left(k+2\right)-i_{\mathrm{dq}}^{*}\left(k+2\right)\right]}^2& \left(13\right)\end{array}$

Wherein i.sub.q.sup.*=0, and

[00026]$i_d^{*}=i_d+\left(V_{\mathrm{dc}}-V_{\mathrm{dc}}^{*}\right)\left(k_p+\frac{k_i}{s}\right);$

k.sub.p, k.sub.i are the adjusted V.sub.dc PI controller gain, a DC bus voltage, V.sub.dc.sup.* is a DC bus voltage reference.

[0078] For n={0,1, . . . ,7} in table 1, a voltage of each switching state corresponds to n={0,1, . . . ,7}, and the voltage u.sub.dq of each switching state is evaluated in the cost function. Among these 8 options, a voltage which obtains a minimum J is applied as the switching state S.sub.abc of the power converter. Evaluation is carried out in the cost function. Among these 8 options, the voltage which obtains the minimum J is applied as the switching state of the power converter.

TABLE-US-00001 TABLE 1 Switching states of a power converter n S abc u dq = T park. u abc 0 000 0 + j0 1 100 [00027]$\frac{V_{dc}}{3}\left(2+j0\right)$ 2 110 [00028]$\frac{V_{dc}}{3}\left(1+j\sqrt{3}\right)$ 3 010 [00029]$\frac{V_{dc}}{3}\left(-1+j\sqrt{3}\right)$ 4 011 [00030]$\frac{V_{dc}}{3}\left(-2+j0\right)$ 5 001 [00031]$\frac{V_{dc}}{3}\left(-1-j\sqrt{3}\right)$ 6 101 [00032]$\frac{V_{dc}}{3}\left(1-j\sqrt{3}\right)$ 7 111 0 + j0

[0079] FIG. 5 shows comparative performance of an adaptive ESO observer and a CESO observer during startup of a power converter in the example. The upper figure in FIG. 5 shows the d-axis current. It can be seen that the peak d-axis current produced by the traditional method (CESO) is 35.55A, while the peak current produced by the method in this example is lower, which is 24.55A. The lower figure in FIG. 5 shows an A phase current, and it can be seen that the method in this example significantly reduces a transient ripple and oscillation compared to the traditional method. This shows the effectiveness of the solution.

Example 2

[0080] In one or more embodiments, disclosed is a predictive control system of a grid-connected converter based on a structurally adaptive ESO, comprising: a data obtaining module, configured for obtaining data of a grid current and a grid voltage and converting the data to data in a dq coordinate system; a voltage vector selecting module, configured for calculating a selected voltage vector based on a voltage in the dq coordinate system; an adaptive state observing module, configured for adaptively inputting a current in the dq coordinate system to a parallel ESO, a cascade ESO, or a hybrid cascade-parallel ESO, so as to obtain a current predicted current value and estimated total disturbance; a two-step predicting module, configured for carrying out a two-step grid current prediction based on the current predicted current value, the estimated total disturbance and the selected voltage vector; a grid-connected converter controlling module, configured for taking minimizing a cost function as a control target, so as to obtain an optimal voltage vector based on a two-step grid current prediction result, and carrying out switching control of the grid-connected converter.

[0081] It should be noted that the specific implementation of the above modules has been explained in detail in example 1, and would not be detailed.

Example 3

[0082] In one or more embodiments, disclosed is a terminal device comprising a server, the server comprises a memory, a processor, and a computer program stored on the memory and capable of being executed by the processor, and the predictive control method of a grid-connected converter based on a structurally adaptive ESO in example 1 is implemented when the processor executes the program. For the sake of concise, it would not be repeated herein.

[0083] It should be understood that in this example, the processor may be a central processing unit (CPU), the processor may also be another general-purpose processor, a digital signal processor (DSP), an application-specific integrated circuit (ASIC), a field programmable gate array (FPGA) or another programmable logic device, a discrete gate or transistor logic device, a discrete hardware component, etc. The general-purpose processor may be a microprocessor or the processor may also be any conventional processor, etc.

[0084] The memory may include a read-only memory and a random access memory and provide instructions and data to the processor, and part of the memory may also include a non-volatile random access memory. For example, the memory may also store information of a device type.

[0085] In an embodiment process, the steps of the above method can be completed by an integrated logic circuit of a hardware in the processor or instructions in the form of a software.

Example 4

[0086] In one or more embodiments, disclosed is a computer-readable storage medium in which a plurality of instructions are stored, wherein the instructions are suitable for being loaded by a processor of a terminal device and executing the predictive control method of a grid-connected converter based on a structurally adaptive ESO in example 1.

[0087] Although the above describes the specific embodiments of the present invention in combination with the attached drawings, it is not a limitation of the protection scope of the present invention, and a person skilled in the art should understand that on the basis of the technical solutions of the present invention, various modifications or variants that can be made by a person skilled in the art without creative labor are still within the protection scope of the present invention.