CN115133802B - Inverter model prediction control method - Google Patents

Inverter model prediction control method Download PDF

Info

Publication number
CN115133802B
CN115133802B CN202210634289.7A CN202210634289A CN115133802B CN 115133802 B CN115133802 B CN 115133802B CN 202210634289 A CN202210634289 A CN 202210634289A CN 115133802 B CN115133802 B CN 115133802B
Authority
CN
China
Prior art keywords
current
inverter
switching
voltage vector
voltage
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202210634289.7A
Other languages
Chinese (zh)
Other versions
CN115133802A (en
Inventor
戴瑜兴
彭子舜
赵振兴
廖石波
朱勇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
SHENZHEN JINGQUANHUA ELECTRONICS CO LTD
Wenzhou University
Original Assignee
SHENZHEN JINGQUANHUA ELECTRONICS CO LTD
Wenzhou University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by SHENZHEN JINGQUANHUA ELECTRONICS CO LTD, Wenzhou University filed Critical SHENZHEN JINGQUANHUA ELECTRONICS CO LTD
Priority to CN202210634289.7A priority Critical patent/CN115133802B/en
Publication of CN115133802A publication Critical patent/CN115133802A/en
Application granted granted Critical
Publication of CN115133802B publication Critical patent/CN115133802B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02MAPPARATUS FOR CONVERSION BETWEEN AC AND AC, BETWEEN AC AND DC, OR BETWEEN DC AND DC, AND FOR USE WITH MAINS OR SIMILAR POWER SUPPLY SYSTEMS; CONVERSION OF DC OR AC INPUT POWER INTO SURGE OUTPUT POWER; CONTROL OR REGULATION THEREOF
    • H02M7/00Conversion of ac power input into dc power output; Conversion of dc power input into ac power output
    • H02M7/42Conversion of dc power input into ac power output without possibility of reversal
    • H02M7/44Conversion of dc power input into ac power output without possibility of reversal by static converters
    • H02M7/48Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode
    • H02M7/53Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal
    • H02M7/537Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal using semiconductor devices only, e.g. single switched pulse inverters
    • H02M7/5387Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal using semiconductor devices only, e.g. single switched pulse inverters in a bridge configuration
    • H02M7/53871Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal using semiconductor devices only, e.g. single switched pulse inverters in a bridge configuration with automatic control of output voltage or current
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/381Dispersed generators

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Inverter Devices (AREA)

Abstract

The invention discloses a predictive control method of an inverter model, which comprises a switch state selection unit, a model predictive control modeling unit and a cost function optimizing unit; according to the invention, only the cost function corresponding to the selected switching state is calculated in the rolling optimization process, and the cost function combines multiple aspects of inverter switching loss, specific order output current content and the like, so that compared with the traditional model prediction control method which adopts a traversal optimization method to obtain the optimal switching state and the cost function which only contains current control, the method provided by the invention has the advantages of reducing the cost function calculation amount, reducing the output current harmonic content, reducing the inverter switching loss and the like.

Description

Inverter model prediction control method
Technical Field
The invention relates to the technical field of automatic control, in particular to a predictive control method for an inverter model.
Background
The model predictive control has good applicability to complex systems, and for nonlinear systems, the nonlinear systems do not need to be solved by adopting an average method in a mode of combining a limited number of switch states with a dynamic model of the system. When the control algorithm is designed, the transfer function of the system does not need to be deduced, and only a prediction model of the system needs to be obtained, so that the systematic control can be realized. The model predictive control can design a cost function containing a plurality of constraint conditions according to an optimization target, and multi-dimensional control is realized. The model predictive control is gradually applied to the field of industrial power electronics by virtue of the advantages of adaptability to complex systems, solving of the problem of power electronics nonlinearity and the like, but the requirements of the control algorithm on a control chip are higher and higher with the improvement of control precision.
The three-phase inverter model predictive controls voltage vectors corresponding to a limited number of switching states as shown in fig. 4, for a total of 8 switching states, 7 different voltage vectors, where V0= V7. In the traditional model predictive control, a traversal method is adopted to calculate cost functions corresponding to all switch states in a period, and the switch state when the cost function is the minimum value is selected to act on the inverter, so that the control purpose is achieved.
Along with the gradual trend of high frequency of power electronic equipment, the sampling period is compressed to a shorter time, and a cheap chip has limited computing capability and is difficult to meet the requirement; although the high performance chip can solve this problem, it also causes problems such as an increase in cost. The traditional model prediction control needs to calculate cost functions corresponding to all switching vectors, and the chip calculation burden is large, so that the working frequency of the converter is difficult to further improve, and even the control effect is influenced because of time delay.
The prior art discloses a model predictive control method for a grid-connected inverter CN112383237A, which can implement and allocate each weight factor of a cost function according to a specific working requirement, so that the switching state combination is more reasonable, but the problem of too high calculation amount of the cost function is not solved substantially.
Disclosure of Invention
In order to solve the problems of the conventional model predictive control, the invention provides a predictive control method of an inverter model, aiming at solving the influence caused by the traversal optimization mode adopted by the conventional model predictive control, reducing the calculation amount of a cost function, reducing the harmonic content of output current, improving the quality of the output waveform of the inverter and reducing the switching loss of the inverter, and the specific technical scheme is as follows:
the model predictive control method of the inverter is model predictive control for reducing cost function calculation amount based on a grid-connected inverter, and comprises the following specific steps:
s1, establishing a discrete model of grid-connected inverter current through a difference method, and converting the discrete model into a model under a static coordinate system;
s2, sampling at the moment k to obtain a network side current i COM (k) And the network side voltage u COM (k) And substituting the voltage into a model under the static coordinate system, and simultaneously detecting the network side voltage u COM (k) Phase, and combining with reference current amplitude, constructing a grid-side voltage u according to the voltage loop COM (k) Reference current i of the same phase ref (k);
S3, judging a voltage vector V corresponding to the optimal switching state at the k-1 moment COMn (k-1) whether the voltage vector is zero or not and calculating the predicted current value of the voltage vector corresponding to the time of k-1
Figure BDA0003681393890000021
S4, designing a cost function according to a switching loss term and a specific order output current error term, and predicting a current value according to the k +1 moment
Figure BDA0003681393890000022
And a reference current i ref Calculating to obtain a cost function g n Comparing the cost functions to obtain a minimum cost function g min To obtain the optimal switching function S with the minimum cost function i According to an optimum switching function S i And controlling the on-off of a switching tube of the inverter to ensure that the error between the grid-test current and the reference current at the moment k +1 is minimum.
Further, the specific steps of S1 are as follows:
establishing a discrete model of the grid-connected inverter current by a difference method:
Figure BDA0003681393890000023
in the formula u a 、u b 、u c Is the network side voltage, i a 、i b 、i c Is a three-phase current, L a 、L b 、L c Filter inductor, R, for connection of an inverter to a power supply system g As a lineEquivalent resistance, U dc Is a direct current side voltage;
introduction of a switching function S for expressing the switching states of the individual legs of an inverter i
Figure BDA0003681393890000024
Using Euler' S formula to correct said switching function S i Discretizing to obtain a discrete model of the grid-connected inverter:
Figure BDA0003681393890000031
in the formula v a 、v b 、v c For each phase of the inverter leg output voltage v a =S a *U dc ,v b =S b *U dc ,v c =S c *U dc ,T s Is a sampling period, due to L a =L b =L c Therefore, L is uniformly used for substitution;
converting the discrete model into a model under a static coordinate system:
Figure BDA0003681393890000032
in the formula
Figure BDA0003681393890000033
To predict the current vector of the current, i COM (k) Is a net side vector current, u COM (k) Is the grid side voltage vector, v COMj (k) And outputting a voltage vector for the bridge arm of the inverter.
Further, the specific steps of S3 are as follows:
if V COMn (k-1) is a non-zero voltage vector, the voltage vector and two voltage vectors V close to this voltage vector are calculated COMx (k)、V COMy (k)、V COMz (k) And zero voltage vector V COM7 (k)(V COM0 (k) Predicted current value corresponding to time k +1
Figure BDA0003681393890000034
And &>
Figure BDA0003681393890000035
/>
If V COMn (k-1) is a zero voltage vector V COM7 (k-1)(V COM0 (k-1)), the predicted current values at the time k +1 corresponding to all the voltage vectors are calculated
Figure BDA0003681393890000036
Further, the specific steps of S4 are as follows:
the content of the m-order output current at the k moment is as follows:
Figure BDA0003681393890000037
the content of the m-order output current at the k +1 th moment is as follows:
Figure BDA0003681393890000038
in the formula I m Is the m-th harmonic amplitude, N is the number of samples in one sampling period,
Figure BDA0003681393890000039
in order to be a factor of rotation,
Figure BDA0003681393890000041
the difference between the m-order output current content at the k moment and the k +1 moment can be obtained as follows:
Figure BDA0003681393890000042
will be provided with
Figure BDA0003681393890000043
And a reference current i ref (k + 1) the difference gives the error current:
Figure BDA0003681393890000044
in the formula i errα (k + 1) and i errβ (k + 1) is the real and imaginary parts of the error current, i refα (k + 1) and i refβ (k + 1) are the real and imaginary parts of the reference current respectively,
Figure BDA0003681393890000045
and &>
Figure BDA0003681393890000046
Respectively the real part and the imaginary part of the predicted current;
calculating the specific subharmonic content of the error current, and when m is 1, obtaining a mathematical expression of the fundamental amplitude component:
Figure BDA0003681393890000047
when m is a specific subharmonic order, a mathematical expression for the harmonic magnitude component can be obtained:
Figure BDA0003681393890000048
taking fundamental wave and harmonic amplitude component as one constraint condition of the cost function, and taking the action times q of the switching device as another constraint condition to obtain the cost function:
Figure BDA0003681393890000049
in the formula, the first term is iterative discrete Fourier transform of fundamental current error, and the second term is iterative discrete Fourier transform of weighted harmonic current errorSum of leaf transforms, m being a specific harmonic order, λ m A weight coefficient for each order harmonic, said weight coefficient lambda m The third term is the action times of the switching device, lambda is a weight coefficient, and q is the total switching action number of the three-way bridge arm.
Further, the number of switching operations in the switching state corresponding to the zero voltage vector is 1, and the number of switching operations in the switching state corresponding to the non-zero voltage vector is 0 or 1.
Preferably, the switch states corresponding to the zero voltage vector include two switch states, 000 and 111.
Preferably, the switch states include 8 states, and 8 voltage vectors, i.e. v, can be taken COMj (k) Wherein j =0,1, 2.
Compared with the prior art, the invention has the beneficial effects that:
according to the invention, the calculation time of the controller is reduced by adopting model prediction control for reducing the calculation amount of the cost function, the cost function is designed according to the switching loss and the specific order output current content, and compared with the traditional model prediction control, the grid-connected inverter has the characteristics of higher response speed, higher working efficiency, less output waveform harmonic content and the like.
Drawings
FIG. 1 is a flowchart of a prediction method according to an embodiment of the present invention.
Fig. 2 is a block diagram illustrating a structure and control of a grid-connected inverter according to an embodiment of the present invention.
Fig. 3 is a flowchart of the embodiment of the present invention applied to a grid-connected inverter.
Fig. 4 is a voltage vector diagram corresponding to each switching state of the grid-connected inverter under the complex plane.
Fig. 5 is a vector diagram of a grid-connected inverter reference current and a grid-side current in a complex plane according to an embodiment of the present invention.
Detailed Description
The present invention is further illustrated by the following description of specific embodiments, which are not intended to limit the invention, and various modifications and improvements can be made by those skilled in the art based on the basic idea of the invention, but the invention is within the protection scope of the invention.
Referring to fig. 1 to 5, an embodiment of the present invention is as follows:
fig. 2 shows a topology and a control block diagram of the grid-connected inverter. In FIG. 2, U dc Is a DC source voltage, L a 、L b 、L c Is an AC side filter inductor, R g Is a line equivalent resistance, i s Is a net side current, wherein i a 、i b 、i c Is a three-phase current of u a 、u b 、u c Is the network side voltage, omega s 、θ S For the grid voltage angular frequency and phase angle, i ref For reference current, n (k + 1) is the serial number corresponding to the optimal switch state at the last moment, V COM0 (k)~V COM7 (k) Is 8 voltage vectors, V COMx (k)、V COMy (k)、V COMz (k) Is three successive voltage vectors, S a 、S b 、S c As a function of inverter leg switching.
The system mainly comprises a grid-connected inverter topology structure chart, a switching state selection unit, a modeling unit and a cost function optimizing unit. In the switching state selection unit, the optimal switching state at the moment is deduced from the switching state value of the inverter at the last moment; in the model prediction control modeling unit and the cost function optimizing unit, network side current and voltage are detected, an inverter prediction model is established, cost functions corresponding to a plurality of switching states obtained from the switching state selection unit are calculated, the switching state when the cost function is the minimum value is selected and acts on the inverter, and the output current of the inverter can better track the reference current.
As shown in fig. 5, a vector diagram of a grid-connected inverter reference current and a grid-side current in a complex plane according to an embodiment of the present invention. The diagram (a) is a partially enlarged diagram of the diagram (b), and in general, the converter operating frequency is high, and the reference current can be considered to be constant in two preceding and succeeding sampling periods. It can be seen from fig. 5 that the optimal switching states of the converter mainly occur in several switching states in several sampling periods, and by calculating the cost functions corresponding to the respective switching states at a time, it can be found that when the voltage vector corresponding to the optimal switching state is a non-zero vector, the corresponding cost function value is minimum, and the cost functions corresponding to the voltage vectors closer to the voltage vector (the included angle with the optimal voltage vector is smaller) are all smaller than the other cost functions, and the cost functions corresponding to the voltage vectors farther away from the optimal voltage vector (the included angle with the optimal voltage vector is larger) are larger. According to the characteristic, only cost functions corresponding to a plurality of cost functions close to the optimal switching state and a zero vector at the first moment are calculated at each moment, and the included angle between the switching vectors corresponding to the rest switching states and the optimal voltage vector at the last moment is large, so that the cost functions corresponding to the voltage vectors are not calculated, and the calculation amount is greatly reduced.
As shown in fig. 1, the implementation steps of the model predictive control for reducing the calculation amount of the cost function based on the grid-connected inverter provided by the invention are as follows:
s1, establishing a discrete model of grid-connected inverter current through a difference method. The discrete model current expression is as follows:
Figure BDA0003681393890000061
in the formula u a 、u b 、u c Is the network side voltage, i a 、i b 、i c Is a three-phase current, L a 、L b 、L c Filter inductor, R, for connection of an inverter to a power supply system g Is a line equivalent resistance, U dc Is the dc side voltage. For expressing the switching states of the individual legs of the converter, the switching function Si is defined as
Figure BDA0003681393890000071
Discretizing the formula (1) by using an Euler formula to obtain a discrete model of the grid-connected inverter, wherein the discrete model is as follows:
Figure BDA0003681393890000072
in the formula v a 、v b 、v c For each phase of the inverter leg output voltage v a =S a *U dc ,v b =S b *U dc ,v c =S c *U dc ,T s Is a sampling period, due to L a =L b =L c Therefore, the model is uniformly replaced by L, and due to the characteristic of the nonlinearity of the switching function, in order to simplify the model, the discrete model is converted into a model in a static coordinate system as follows:
Figure BDA0003681393890000073
in the formula
Figure BDA0003681393890000074
To predict the current vector of the current, i COM (k) Is a net side vector current, u COM (k) Is the grid side voltage vector, v COMj (k) And outputting a voltage vector for the bridge arm of the inverter.
As shown in FIG. 4, the inverter legs have 8 different switching states, v COMj (k) 8 voltage vectors may be taken, where j =0,1, 2.
S2, sampling at the moment k to obtain a network side current i COM (k) And the network side voltage u COM (k) And substituting it into formula (4); and detecting the phase of the network side voltage, and constructing a reference current i in the same phase with the network side voltage through a voltage loop ref (k)。
S3, judging a voltage vector V corresponding to the optimal switching state at the k-1 moment COMn (k-1) is a zero voltage vector. If it is a non-zero voltage vector, the voltage vector and two voltage vectors V close to the voltage vector are calculated COMx (k)、V COMy (k)、V COMz (k) And zero voltage vector V COM7 (k)(V COM0 (k) Predicted current value corresponding to time k +1
Figure BDA0003681393890000075
And &>
Figure BDA0003681393890000076
When the voltage vector corresponding to the optimal switching state at the moment of k-1 is V COM3 (k-1), calculating a voltage vector V COM2 (k),V COM3 (k),V COM4 (k) And V COM7 (k) Corresponding predicted current value
Figure BDA0003681393890000081
And &>
Figure BDA0003681393890000082
When the voltage vector corresponding to the optimal switching state at the moment of k-1 is V COM6 (k-1), then the voltage vector is calculated as V COM5 (k),V COM6 (k),V COM1 (k) And V COM7 (k) Time-corresponding predicted current value
Figure BDA0003681393890000083
And &>
Figure BDA0003681393890000084
If V COMn (k-1) is a zero voltage vector V COM7 (k-1)(V COM0 (k-1)), the predicted current value ^ at the time k +1 corresponding to all the voltage vectors is calculated>
Figure BDA0003681393890000085
Step 4, reference current i in two moments before and after ref (k) And i ref (k + 1) is approximately the same, so i can be considered to be ref (k+1)=i ref (k) In that respect Designing a cost function according to the switching loss and the specific order output current content, and calculating the specific order output current content by iterative discrete Fourier calculation, wherein the specific method comprises the following steps: the content of the m-order output current at the k moment is as follows:
Figure BDA0003681393890000086
the content of the m-order output current at the k +1 th moment is as follows:
Figure BDA0003681393890000087
in the formula I m Is the m-th harmonic amplitude, N is the number of samples in one sampling period,
Figure BDA0003681393890000088
in order to be a factor of rotation,
Figure BDA0003681393890000089
the difference between the formula (5) and the formula (6) can be obtained,
Figure BDA00036813938900000810
the formula (7) is an iterative discrete Fourier transform expression, and the amplitude component of the specific secondary output current is obtained through a few calculation processes. Will be provided with
Figure BDA00036813938900000811
And a reference current i ref (k + 1) the difference gives the error current:
Figure BDA00036813938900000812
in the formula i errα (k + 1) and i errβ (k + 1) is the real and imaginary parts of the error current, i refα (k + 1) and i refβ (k + 1) are the real and imaginary parts of the reference current respectively,
Figure BDA00036813938900000813
and &>
Figure BDA00036813938900000814
The real and imaginary parts of the predicted current, respectively. MeterCalculating the specific subharmonic content of the error current, and when m is 1, obtaining a mathematical expression of the fundamental amplitude component as follows:
Figure BDA0003681393890000091
when m is a specific subharmonic order, the mathematical expression for the harmonic amplitude component is given as:
Figure BDA0003681393890000092
/>
the equations (9) and (10) are mathematical expressions of fundamental wave and harmonic amplitude component, which are taken as one constraint condition of the cost function, and the number q of switching device actions is taken as another constraint condition:
Figure BDA0003681393890000093
wherein the first term is the iterative discrete Fourier transform of the fundamental current error, the second term is the sum of the iterative discrete Fourier transforms of the weighted harmonic current error, m is a specific harmonic order, and λ m A weight coefficient for each order harmonic, said weight coefficient lambda m The third term is the action times of the switching device, lambda is a weight coefficient, and q is the total switching action number of the three-way bridge arm. Predicting current value according to k +1 time
Figure BDA0003681393890000094
And a reference current i ref Calculating to obtain a cost function g n Comparing the cost functions to obtain a minimum cost function g min The detailed comparison process is shown in FIG. 3: when n is 2, 3, 4 and 5, g is calculated and compared n-1 、g n 、g n+1 And g 0 (ii) a When n is 7, calculating and comparing the remaining four cost functions; when n is 1, g is calculated and compared 6 、g 1 、g 2 And g 0 (ii) a When n is 6, calculate g 5 、g 6 、g 1 And g 0 Thereby obtaining a cost function g corresponding to the output minimum voltage vector min And then obtaining an optimal switching function S which minimizes the cost function i According to an optimum switching function S i And controlling the on-off of a switching tube of the inverter to ensure that the error between the grid-measured current and the reference current at the moment k +1 is minimum.
The method for reducing the calculated amount of the cost function is designed according to the distribution condition of each cost function in each moment and the characteristics of output voltage at two adjacent moments, the cost function is designed according to the switching loss and the content of the specific-order output current, the accuracy of the output current is ensured while the calculated amount of the cost function is reduced, and the method has great application significance for improving the working frequency of a converter and reducing the calculated amount of a chip.
The foregoing embodiments are merely illustrative of the principles and utilities of the present invention and are not intended to limit the invention. It will be appreciated by those skilled in the art that changes may be made in these embodiments without departing from the principles and spirit of the invention, the scope of which is defined in the claims and their equivalents.

Claims (7)

1. The inverter model prediction control method is characterized by being based on the model prediction control of a grid-connected inverter for reducing cost function calculation amount, and specifically comprising the following steps of:
s1, establishing a discrete model of grid-connected inverter current through a difference method, and converting the discrete model into a model under a static coordinate system;
s2, sampling at the moment k to obtain a network side current i COM (k) And the network side voltage u COM (k) And substituting the voltage into a model under the static coordinate system, and simultaneously detecting the network side voltage u COM (k) Phase, and combining with reference current amplitude, constructing a grid-side voltage u according to the voltage loop COM (k) Reference current i of the same phase ref (k);
S3, judging a voltage vector V corresponding to the optimal switching state at the k-1 moment COMn (k-1) whether the voltage vector is zero or not and calculating the predicted current value of the voltage vector corresponding to the time of k-1
Figure FDA0003681393880000011
S4, designing a cost function according to a switching loss term and a specific order output current error term, and predicting a current value according to the k +1 moment
Figure FDA0003681393880000012
And a reference current i ref Calculating to obtain a cost function g n Comparing the cost functions to obtain a minimum cost function g min And then obtaining an optimal switching function S which minimizes the cost function i According to an optimum switching function S i And controlling the on-off of a switching tube of the inverter to ensure that the error between the grid-test current and the reference current at the moment k +1 is minimum.
2. The inverter model predictive control method of claim 1, wherein the specific steps of S1 are as follows:
establishing a discrete model of the grid-connected inverter current by a difference method:
Figure FDA0003681393880000013
in the formula u a 、u b 、u c Is the network side voltage, i a 、i b 、i c Is a three-phase current, L a 、L b 、L c Filter inductance, R, for connection of an inverter to a power grid g Is a line equivalent resistance, U dc Is a direct current side voltage;
introduction of a switching function S for expressing the switching states of the individual legs of an inverter i
Figure FDA0003681393880000021
Using Euler' S formula to correct said switching function S i Discretizing to obtain a discrete model of the grid-connected inverter:
Figure FDA0003681393880000022
in the formula v a 、v b 、v c For each phase of the inverter leg output voltage v a =S a *U dc ,v b =S b *U dc ,v c =S c *U dc ,T s Is a sampling period, due to L a =L b =L c Therefore, L is uniformly used for substitution;
converting the discrete model into a model under a static coordinate system:
Figure FDA0003681393880000023
in the formula
Figure FDA0003681393880000024
To predict the current vector of the current, i COM (k) Is a net side vector current, u COM (k) Is the grid side voltage vector, v COMj (k) And outputting a voltage vector for the bridge arm of the inverter.
3. The inverter model predictive control method of claim 1, wherein the specific steps of S3 are as follows:
if V COMn (k-1) is a non-zero voltage vector, the voltage vector and two voltage vectors V close to this voltage vector are calculated COMx (k)、V COMy (k)、V COMz (k) And zero voltage vector V COM7 (k)(V COM0 (k) Predicted current value corresponding to time k +1
Figure FDA0003681393880000025
And
Figure FDA0003681393880000026
if V COMn (k-1) is a zero voltage vector V COM7 (k-1)(V COM0 (k-1)), the predicted current values at the time k +1 corresponding to all the voltage vectors are calculated
Figure FDA0003681393880000027
4. The inverter model predictive control method of claim 1, wherein the specific steps of S4 are as follows:
the content of the m-order output current at the k moment is as follows:
Figure FDA0003681393880000028
the content of the m-order output current at the k +1 th moment is as follows:
Figure FDA0003681393880000031
in the formula I m Is the m-th harmonic amplitude, N is the number of samples in one sampling period,
Figure FDA0003681393880000032
is a factor of the rotation of the optical fiber,
Figure FDA0003681393880000033
the difference between the m-order output current content at the k moment and the k +1 moment can be obtained as follows:
Figure FDA0003681393880000034
will be provided with
Figure FDA0003681393880000035
And a reference current i ref (k + 1) the difference gives the error current:
Figure FDA0003681393880000036
in the formula i errα (k + 1) and i errβ (k + 1) is the real and imaginary parts of the error current, i refα (k + 1) and i refβ (k + 1) are the real and imaginary parts of the reference current respectively,
Figure FDA0003681393880000037
and
Figure FDA0003681393880000038
respectively the real part and the imaginary part of the predicted current;
calculating the specific subharmonic content of the error current, and when m is 1, obtaining a mathematical expression of the fundamental amplitude component:
Figure FDA0003681393880000039
when m is a specific subharmonic order, a mathematical expression for the harmonic magnitude component can be obtained:
Figure FDA00036813938800000310
taking fundamental wave and harmonic amplitude component as one constraint condition of the cost function, and taking the action times q of the switching device as another constraint condition to obtain the cost function:
Figure FDA00036813938800000311
wherein the first term is the iterative discrete Fourier transform of the fundamental current error, the second term is the sum of the iterative discrete Fourier transforms of the weighted harmonic current error, m is a specific harmonic order, and λ m A weight coefficient for each order harmonic, said weight coefficient lambda m The third term is the action times of the switching device, lambda is a weight coefficient, and q is the total switching action number of the three-way bridge arm.
5. The inverter model predictive control method of claim 1, wherein the number of switching operations of the switching state of the zero voltage vector is 1, and the number of switching operations of the switching state of the non-zero voltage vector corresponding to the zero voltage vector is 0 or 1.
6. The inverter model predictive control method of claim 5, wherein the zero voltage vector corresponds to switching states comprising 000 and 111 switching states.
7. The inverter model predictive control method of claim 5, wherein the switching states include 8 states, and 8 voltage vectors, i.e., v, can be taken COMj (k) Wherein j =0,1, 2.
CN202210634289.7A 2022-06-07 2022-06-07 Inverter model prediction control method Active CN115133802B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210634289.7A CN115133802B (en) 2022-06-07 2022-06-07 Inverter model prediction control method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210634289.7A CN115133802B (en) 2022-06-07 2022-06-07 Inverter model prediction control method

Publications (2)

Publication Number Publication Date
CN115133802A CN115133802A (en) 2022-09-30
CN115133802B true CN115133802B (en) 2023-03-24

Family

ID=83378952

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210634289.7A Active CN115133802B (en) 2022-06-07 2022-06-07 Inverter model prediction control method

Country Status (1)

Country Link
CN (1) CN115133802B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115765508A (en) * 2022-11-25 2023-03-07 天津大学 Prediction control method for equivalent space vector model of modular multilevel converter
CN116780929B (en) * 2023-05-04 2024-01-30 东莞市海柯电子有限公司 Variable weight prediction control method based on SiWBG cascade H-bridge converter
CN116707342B (en) * 2023-08-08 2023-12-19 安徽大学 Data-driven energy storage inverter voltage prediction control method

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110045610A (en) * 2019-04-18 2019-07-23 中国地质大学(武汉) Inverter modified multistep model predictive control method, equipment and storage equipment
CN113991739A (en) * 2021-11-05 2022-01-28 福州大学 Simplified vector fixed-frequency prediction current control method for grid-connected inverter
WO2022105068A1 (en) * 2020-11-17 2022-05-27 江苏时代新能源科技有限公司 Control method and device for grid-connected inverter

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104680424B (en) * 2015-01-30 2018-09-14 国家电网公司 The voltage power situation predictor method in large-sized photovoltaic power station

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110045610A (en) * 2019-04-18 2019-07-23 中国地质大学(武汉) Inverter modified multistep model predictive control method, equipment and storage equipment
WO2022105068A1 (en) * 2020-11-17 2022-05-27 江苏时代新能源科技有限公司 Control method and device for grid-connected inverter
CN113991739A (en) * 2021-11-05 2022-01-28 福州大学 Simplified vector fixed-frequency prediction current control method for grid-connected inverter

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
模型预测控制在光伏并网逆变器中的应用;潘健等;《湖北工业大学学报》;20171015;第32卷(第05期);82-84,92 *

Also Published As

Publication number Publication date
CN115133802A (en) 2022-09-30

Similar Documents

Publication Publication Date Title
CN115133802B (en) Inverter model prediction control method
Machado et al. A neural network-based dynamic cost function for the implementation of a predictive current controller
CN103595069B (en) Photovoltaic generating system grid side converter model predictive control method under unbalance voltage
Sujith et al. Optimization of harmonics with active power filter based on ADALINE neural network
CN112994482B (en) Model prediction power dynamic compensation control method of double-PWM frequency conversion system
Canciello et al. Multi-objective adaptive sliding manifold control for More Electric Aircraft
WO2011095947A2 (en) Method and system for facilitating design of a high voltage (hvdc) control system, an hvdc system and a method for optimising an hvdc system
CN108347190B (en) Rapid prediction method based on modular multilevel converter structure
CN108959671B (en) Real-time simulation modeling method for half-bridge and full-bridge modular multilevel converter
Merabet et al. Direct neural method for harmonic currents estimation using adaptive linear element
Balouji et al. Deep learning based predictive compensation of flicker, voltage dips, harmonics and interharmonics in electric arc furnaces
CN108347057B (en) LCC-MMC mixed direct-current power grid alternating iteration load flow calculation method
CN109004852B (en) Model prediction control strategy of modular multilevel converter
Maciążek Power theories applications to control active compensators
Soodi et al. Design, Optimization and Experimental Verification of a Low Cost Two-Microcontroller Based Single-Phase STATCOM
WO2024082730A1 (en) Finite control set model predictive control method for llcl-type battery energy storage converter
Thirumoorthi et al. Adaptive method for power quality improvement through minimization of harmonics using artificial intelligence
Noda et al. A practical steady-state initialization method for electromagnetic transient simulations
CN108767887B (en) Sub-module capacitor voltage prediction correction method for high-voltage direct-current transmission
Weaver et al. Real-time Hardware-in-the-Loop simulation for optimal Dc microgrid control development
Dai et al. An improved dragonfly algorithm with higher exploitation capability to optimize the design of hybrid power active filter
CN114597955A (en) Three-phase LCL grid-connected NPC inversion system based on rapid model predictive control
CN113452069A (en) Multi-step model prediction control method for inverter station supplying power to passive network
CN107276123A (en) A kind of variable coefficient frequency reducing model prediction direct Power Control method and device
CN113904578A (en) Non-weight coefficient model prediction control method for single-phase cascaded H-bridge converter

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20230208

Address after: 518000 jingquanhua Industrial Park, No. 325, Guiyue Road, Guanlan street, Longhua District, Shenzhen, Guangdong

Applicant after: Shenzhen JingQuanHua Electronics Co.,Ltd.

Applicant after: Wenzhou University

Address before: 518000 jingquanhua Industrial Park, No. 325, Guiyue Road, Guanlan street, Longhua District, Shenzhen, Guangdong

Applicant before: Shenzhen JingQuanHua Electronics Co.,Ltd.

GR01 Patent grant
GR01 Patent grant