PREDICTIVE CONTROL METHOD OF CURRENT INCREMENT FOR PERMANENT MAGNET SYNCHRONOUS MOTOR UNDER HIGH-SPEED OPERATION

Abstract

The present disclosure provides a predictive control method of current increment for a permanent magnet synchronous motor includes: substituting a mathematical expression of a stator voltage during one control period into a continuous time domain current model to obtain a discrete current prediction model and a predicted current at the next time point; obtaining a predicted current increment from a current increment prediction model by subtracting a predictive current at a present time point from a predictive current at a next time point; establishing a cost function according to a preset reference current increment and the predicted current increment; obtaining an optimal voltage increment by minimizing the cost function; superposing the optimal voltage increment on a stator voltage of a present control period to obtain an optimal stator voltage of a next control period for controlling control the permanent magnet synchronous motor.

Claims

1. A predictive control method of current increment for a permanent magnet synchronous motor under a high-speed operation, comprising the following steps: (1) establishing a mathematical expression of a stator voltage during one control period according to a position change of a rotor of the permanent magnet synchronous motor during the one control period; (2) obtaining a continuous time domain current model of the permanent magnet synchronous motor by solving a continuous time domain equation of the permanent magnet synchronous motor; (3) substituting the mathematical expression of the stator voltage during the one control period into the continuous time domain current model of the permanent magnet synchronous motor for solving solutions to obtain a discrete current prediction model suitable for the permanent magnet synchronous motor under the high-speed operation, and then obtaining a predicted current at a next time point by using the discrete current prediction model; (4) subtracting a predictive current at a present time point from a predictive current at a next time point to obtain a current increment prediction model suitable for the permanent magnet synchronous motor under the high-speed region, and then obtaining a predicted current increment calculated from the current increment prediction model; (5) establishing a cost function by taking a squared error at an end of each control period between a preset reference current increment and a predicted current increment as an evaluation criterion, and then minimizing the cost function by solving a convex optimization problem for the cost function to obtain an optimal voltage increment; (6) superposing the optimal voltage increment on a stator voltage of a present control period to obtain an optimal stator voltage of a next control period, and then apply the optimal stator voltage of the next control period to the permanent magnet synchronous motor.

2. The predictive control method of the current increment for the permanent magnet synchronous motor under the high-speed operation according to claim 1, wherein in step (1), the mathematical expression of the stator voltage during the one control period is: $\left[\begin{array}{c}{u}_d\left(t\right)\\ u_q\left(t\right)\end{array}\right]=\left[\begin{array}{cc}\cos \left[\left(t-k{T}_s\right){\omega }_r\right]& \sin \left[\left(t-k{T}_s\right){\omega }_r\right]\\ -s\mathrm{in}\left[\left(t-k{T}_s\right){\omega }_r\right]& \cos \left[\left(t-k{T}_s\right){\omega }_r\right]\end{array}\right]\left[\begin{array}{c}{u}_{d,k}\\ u_{q,k}\end{array}\right]$ where u.sub.d(t) and u.sub.q(t) are a d-axis stator voltage and a q-axis stator voltage, respectively; T.sub.s is the one control period; u.sub.d,k is a d-axis component of a stator voltage vector at time point kT.sub.s, and u.sub.q,k is a q-axis component of the stator voltage vector at time point kT.sub.s, where the subscript d represents d-axis, the subscript q represents q-axis, and the subscript k represents an ordinal number of a control period; ω.sub.r is an electrical angular velocity of the rotor of the permanent magnet synchronous motor; k represents the ordinal number of the control period; t represents a present time point.

3. The predictive control method of the current increment for the permanent magnet synchronous motor under the high-speed operation according to claim 1, wherein in step (3), the discrete current prediction model is: $i_s\left(k+1\right)=A_0\left(k\right)i_s\left(k\right)+B_0\left(k\right)u_s\left(k\right)+D_0\left(k\right)$ $A_0\left(k\right)=\left[\begin{array}{cc}\cos \left({\omega }_{r,k}{T}_s\right)& \frac{L_{q0}}{L_{d0}}\sin \left({\omega }_{r,k}{T}_s\right)\\ -\frac{L_{d0}}{L_{q0}}\sin \left({\omega }_{r,k}{T}_s\right)& \cos \left({\omega }_{r,k}{T}_s\right)\end{array}\right]$ $B_0\left(k\right)=\left[\begin{array}{cc}\frac{T_s}{L_{d0}}\cos \left({\omega }_{r,k}{T}_s\right)& \frac{T_s}{L_{d0}}\sin \left({\omega }_{r,k}{T}_s\right)\\ -\frac{T_s}{L_{q0}}\sin \left({\omega }_{r,k}{T}_s\right)& \frac{T_s}{L_{q0}}\cos \left({\omega }_{r,k}{T}_s\right)\end{array}\right]$ $D_0\left(k\right)={\left[\frac{{\psi }_{f0}\left[\cos \left({\omega }_{r,k}{T}_s\right)-1\right]}{L_{d0}}-\frac{{\psi }_{f0}\sin \left({\omega }_{r,k}{T}_s\right)}{L_{q0}}\right]}^T$ where i.sub.s(k+1) represents a predicted current vector at time point (k+1)T.sub.s; i.sub.s(k) represents a stator current vector at time point kT.sub.s; u.sub.s(k) represents a stator voltage vector at time point kT.sub.s; A.sub.0(k) represents a coefficient matrix of i.sub.s(k); B.sub.0(k) represents a coefficient matrix of u.sub.s(k); D.sub.0(k) represents a coefficient matrix of a back electromotive force; ω.sub.r,k is an electrical angular velocity of the rotor of the permanent magnet synchronous motor at time point kT.sub.s, where r represents a symbol related to the rotor; k represents an ordinal number of a control period; L.sub.d0, L.sub.q0, and ψ.sub.f0 are nominal values of a d-axis stator inductance, a q-axis stator inductance, and a permanent magnet flux linkage, respectively.

4. The predictive control method of the current increment for the permanent magnet synchronous motor under the high-speed operation according to claim 1, wherein in step (4), the current increment prediction model is:
Δi.sub.s(k+1)=A.sub.0(k)Δi.sub.s(k)+B.sub.0(k)Δu.sub.s(k)
Δi.sub.s(k)=i.sub.s(k)−i.sub.s(k−1)
Δu.sub.s(k)=u.sub.s(k)−u.sub.s(k−1) where, Δi.sub.s(k+1) represents the predicted current increment calculated from the current increment prediction model; Δi.sub.s(k) represents a stator current increment between a stator current at time point kT.sub.s and a stator current at time point (k−1)T.sub.s; Δu.sub.s(k) represents a stator voltage increment between a stator voltage at time point kT.sub.s and a stator voltage at time point (k−1)T.sub.s; i.sub.s(k−1) represents a stator current vector at time point; (k−1)T.sub.s; u.sub.s(k−1) represents a stator voltage vector at time point (k−1)T.sub.s; A.sub.0(k) represents a coefficient matrix of Δi.sub.s(k); B.sub.0(k) represents a coefficient matrix of Δu.sub.s(k).

5. The predictive control method of the current increment for the permanent magnet synchronous motor under the high-speed operation according to claim 1, wherein in step (5), the cost function is established as: $\underset{\Delta u_s}{\min }J={\left[\Delta i_s^{\mathrm{ref}}-\Delta i_s\left(k+2\right)\right]}^T\left[\Delta i_s^{\mathrm{ref}}-\Delta i_s\left(k+2\right)\right]+\Delta {u_s\left(k+1\right)}^{T}P\Delta u_s\left(k+1\right)$ $\mathrm{Satisfy}:\{\begin{array}{c}.Math.u_s\left(k\right)+\Delta u_s\left(k+1\right).Math.\le U_{\max }\\ .Math.i_s\left(k+1\right)+\Delta i_s\left(k+2\right).Math.\le I_{\max }\end{array}$ where Δi.sub.s.sup.ref represents the preset reference current increment; P is a weight factor matrix; U.sub.max and I.sub.max are a maximal voltage and a maximal current of a permanent magnet synchronous motor driving system; the superscript T represents a matrix transpose operation; “Satisfy” represents constraint conditions; Δi.sub.s(k+2) represents the predicted current increment calculated from the current increment prediction model; Δu.sub.s(k+1) represents a stator voltage increment between a stator voltage at time point (k+1)T.sub.s and a stator voltage at time point kT.sub.s; J is a value of the cost function.

6. The predictive control method of the current increment for the permanent magnet synchronous motor under the high-speed operation according to claim 1, wherein in step (6), the optimal voltage increment is added to the stator voltage of the present control period to obtain the optimal stator voltage of the next control period, and the optimal stator voltage is:
u.sub.s.sup.opt(k+1)=u.sub.s(k)+Δu.sub.s.sup.opt(k+1) where u.sub.s(k) represents a stator voltage at time point kT.sub.s; u.sub.s.sup.opt(k+1) represents an optimal stator voltage at time point (k+1)T.sub.s; Δu.sub.s.sup.opt(k+1) represents an optimal voltage increment from time point kT.sub.s to time point (k+1)T.sub.s.

Description

BRIEF DESCRIPTION OF DRAWINGS

[0025] FIG. 1 is the control block diagram which illustrates the predictive control method of current increment suitable for the permanent magnet synchronous motor under the high-speed operation.

[0026] FIG. 2 is the simulation current waveforms of the conventional current predictive model when the deadtime is set to 3 μs and the simulation current waveforms of the current increment prediction model.

[0027] FIG. 3 is the current prediction error waveforms of the current increment prediction model with inductance mismatch under different motor operating conditions.

[0028] FIG. 4 is the experimental waveforms of d-axis current, q-axis current and a phase current under different motor operating conditions for the predictive current control method based on the conventional current predictive model and for the current increment prediction model.

DESCRIPTION OF EMBODIMENTS

[0029] Embodiments of the predictive control method of the current increment suitable for the permanent magnet synchronous motor under high-speed operation are explained with reference to the drawings.

[0030] In the following, the method of the disclosure is further introduced based on the detailed principle and situation:

[0031] 1. Establish the Model of the Permanent Magnet Synchronous Motor:

[0032] The rotating coordinate is established and the d-axis aligns on the rotor flux. The continuous time domain model of the permanent magnet synchronous motor is

[00004]$\begin{array}{cc}\frac{{\mathrm{di}}_s}{\mathrm{dt}}=A_s{i}_s+B_s{u}_s+D_s& \left(1\right)\end{array}$

[0033] In the equation (1),

[00005]$\begin{array}{ccc}{A}_s=\left[\begin{array}{cc}-\frac{r_{s0}}{L_{d0}}& {\omega }_r\frac{L_{q0}}{L_{d0}}\\ -{\omega }_r\frac{L_{d0}}{L_{q0}}& -\frac{r_{s0}}{L_{q0}}\end{array}\right];& B_s\left[\begin{array}{cc}\frac{1}{L_{d0}}& 0\\ 0& -\frac{r_{s0}}{L_{q0}}\end{array}\right];& D_s=\left[\begin{array}{c}0\\ -\frac{{\omega }_r{\psi }_{f0}}{L_{q0}}\end{array}\right]\end{array};$

[0034] where r.sub.s0, L.sub.d0, L.sub.q0, and ψ.sub.f0 are the nominal values of stator resistance, d-axis inductance, q-axis inductance, and permanent magnet flux linkage, respectively; ω.sub.r is the electrical angular velocity; A.sub.s represents the coefficient matrix of the current item; B.sub.s represents the coefficient matrix of the voltage item; D.sub.s represents the coefficient matrix related to the back electromotive force; i.sub.s(t)=[i.sub.d(t), i.sub.q(t)].sup.T, where i.sub.d(t), and i.sub.q(t) are the d-axis stator current and the q-axis stator current, respectively; u.sub.s(t)=[u.sub.d(t), u.sub.q(t)].sup.T, where u.sub.d(t) and u.sub.q(t) are the d-axis stator voltage and the q-axis stator voltage, respectively; t represents the time.

[0035] By solving (1), the current model of the permanent magnet synchronous motor in continuous time domain is expressed as

[00006]$\begin{array}{cc}{i}_s\left(t\right)=e^{\left(t-{\mathrm{kT}}_s\right)A_s}{i}_s\left(k\right)+\left[e^{\left(t-{\mathrm{kT}}_s\right)A_s}-I\right]A_s^{-1}{D}_s+{\int }_{{\mathrm{kT}}_s}^t{e}^{\left(t-\tau \right)A_s}{B}_s{u}_s\left(\tau \right)d\tau & \left(2\right)\end{array}$

[0036] In equation (2), T.sub.s is the control period; I is the identity matrix.

[0037] 2. Establish the Discrete Current Prediction Model of the Permanent Magnet Synchronous Motor Under High-Speed Operation Condition:

[0038] The existing predictive current control methods commonly employ the prediction model derived from the first-order forward Euler approximation method which assumes that the value of ω.sub.rT.sub.s is small enough so that the variation of rotor position angle during one control period can be ignored. The exponential term e.sup.(t-kTs)As is equivalently simplified as (t−kT.sub.s)A.sub.s+I. Substituting the above assumptions into equation (2) and discretizing equation (2), the conventional current prediction model obtained by using the first-order forward Euler approximation method is discretized as

[00007]$\begin{array}{cc}{i}_s\left(k+1\right)=A_{c0}\left(k\right)i_s\left(k\right)+B_{c0}\left(k\right)u_s\left(k\right)+D_{c0}\left(k\right)& \left(3\right)\end{array}$$\begin{array}{cc}{A}_{c0}\left(k\right)=\left[\begin{array}{cc}1-\frac{T_s{R}_{s0}}{L_{d0}}& {\omega }_{r,k}\frac{T_s{L}_{q0}}{L_{q0}}\\ -{\omega }_{r,k}\frac{T_s{L}_{d0}}{L_{q0}}& 1-\frac{T_s{R}_{s0}}{L_{q0}}\end{array}\right]\text{⁠};& B_{c0}\left(k\right)=\left[\begin{array}{cc}\frac{T_s}{L_{d0}}& 0\\ 0& \frac{T_s}{L_{q0}}\end{array}\right];\end{array}$$D_c\left(k\right)={\left[\begin{array}{cc}0& -\frac{T_s{\omega }_{r,k}{\psi }_{f0}}{L_{q0}}\end{array}\right]}^T;$

[0039] In equation (3), i.sub.s(k+1)=[i.sub.d,k+1, i.sub.q,k+1].sup.T represents the predicted current vector at time point (k+1)T.sub.s, and i.sub.s(k+1)=[i.sub.d,k+1, i.sub.q,k+1].sup.T, where i.sub.d,k+1 and i.sub.q,k+1 are the d-axis predictive current and q-axis predicted current at time point (k+1)T.sub.s, where the subscript d represents d-axis, the subscript q represents q-axis, and the subscript (k+1) represents the (k+1)th control period; i.sub.s(k) represents the stator current vector at kT.sub.s time point, and i.sub.s(k)=[i.sub.d,k, i.sub.q,k].sup.T, where i.sub.d,k and i.sub.q,k are d-axis current and q-axis current at time point kT.sub.s, where the subscript k represents the kth control period t; u.sub.s(k) represents the stator voltage vector at time point kT.sub.s, and u.sub.s(k)=[u.sub.d,k, u.sub.q,k].sup.T, where u.sub.d,k and u.sub.q,k are d-axis stator voltage and q-axis stator voltage at time point kT.sub.s; A.sub.c0(k) represents the coefficient matrix of the stator current item at time point kT.sub.s; B.sub.c0(k) represents the coefficient matrix of the stator voltage item at time point kT.sub.s; D.sub.c0(k) represents the coefficient matrix related to the back electromotive force at time point kT.sub.s; ω.sub.r,k is the electrical angular velocity at time point kT.sub.s, where r indicates a symbol related to the rotor; k represents the ordinal number of the control period.

[0040] However, the assumption of e.sup.(t-kTs)As≈(t−kT.sub.s)A.sub.s+I is invalid when the motor works under the high-speed operation, and the variation of rotor position angle during one control period cannot be ignored. This disclosure considers the variation of rotor position angle during one control period, and the stator voltage u.sub.s(t)=[u.sub.d(t),u.sub.q(t)].sup.T in equation (2) during one control period can be expressed as

[00008]$\begin{array}{cc}\left[\begin{array}{c}{u}_d\left(t\right)\\ u_q\left(t\right)\end{array}\right]=\left[\begin{array}{cc}\cos \left[\left(t-{\mathrm{kT}}_s\right){\omega }_r\right]& \sin \left[\left(t-{\mathrm{kT}}_s\right){\omega }_r\right]\\ -\sin \left[\left(t-{\mathrm{kT}}_s\right){\omega }_r\right]& \cos \left[\left(t-{\mathrm{kT}}_s\right){\omega }_r\right]\end{array}\right]\left[\begin{array}{c}{u}_{d,k}\\ u_{q,k}\end{array}\right]& \left(4\right)\end{array}$

[0041] In equation (4), u.sub.d,k and u.sub.q,k are d-axis voltage and q-axis voltage at time point kT.sub.s, respectively; kT.sub.s≤t≤(k+1)T.sub.s.

[0042] Ignoring the stator resistance voltage drop and substituting equation (4) into equation (2), the discrete current prediction model is obtained

[00009]$\begin{array}{cc}{i}_s\left(k+1\right)=A_0\left(k\right)i_s\left(k\right)+B_0\left(k\right)u_s\left(k\right)+D_0\left(k\right)& \left(5\right)\end{array}$$A_0\left(k\right)=\left[\begin{array}{cc}\cos \left({\omega }_{r,k}{T}_s\right)& \frac{L_{q0}}{L_{d0}}\sin \left({\omega }_{r,k}{T}_s\right)\\ -\frac{L_{d0}}{L_{q0}}\sin \left({\omega }_{r,k}{T}_s\right)& \cos \left({\omega }_{r,k}{T}_s\right)\end{array}\right];B_0\left(k\right)=\left[\begin{array}{cc}\frac{T_s}{L_{d0}}\cos \left({\omega }_{r,k}{T}_s\right)& \frac{T_s}{L_{d0}}\sin \left({\omega }_{r,k}{T}_s\right)\\ -\frac{T_s}{L_{q0}}\sin \left({\omega }_{r,k}{T}_s\right)& \frac{T_s}{L_{q0}}\cos \left({\omega }_{r,k}{T}_s\right)\end{array}\right]$$D_0\left(k\right)=\left[\begin{array}{c}\frac{{\psi }_{f0}}{L_{d0}}[\cos \left({\omega }_{r,k}{T}_s\right)-1\\ -\frac{{\psi }_{f0}}{L_{q0}}\sin \left({\omega }_{r,k}{T}_s\right)\end{array}\right].$

[0043] In equation (5), A.sub.0(k) represents the coefficient matrix of the stator current item at time point kT.sub.s; B.sub.0(k) represents the coefficient matrix of the stator voltage item at time point kT.sub.s; D.sub.0(k) represents the coefficient matrix related to the back electromotive force at time point kT.sub.s.

[0044] Compared with equation (3), equation (5) considers the influence of rotor movement in each control period on the actual operation trajectories of stator current and voltage, so that it can reflect the change of the stator current in one control period more accurately. However, the dead-time effect and motor parameter mismatch still cause the prediction error.

[0045] 3. Establish the Current Increment Prediction Model of the Permanent Magnet Synchronous Motor Under the High-Speed Operation.

[0046] The inverter output voltage error caused by the dead-time effect is related to the three-phase switching states of the inverter and the directions of three-phase currents. Because the three-phase switching mode of the inverter is fixed, and the directions of the three-phase currents do not change frequently, so the voltage errors between two adjacent control periods caused by the dead-time effect can be seemed to be equal. Therefore, the voltage error caused by the dead-time effect can be eliminated to a certain extent by subtracting the stator voltages from another one stator voltage in two adjacent control period. In the motor drives, ω.sub.r can be seemed to be constant during two adjacent control periods since the control period is short enough, so A.sub.0(k), B.sub.0(k), and D.sub.0(k) can be seemed to be constant during two adjacent control periods. Subtracting the predicted current at time point (k−1)T.sub.s from the predictive current at time point kT.sub.s based on equation (5), the current increment predictive model appropriate for permanent magnet synchronous motor under the high-speed operation is obtained as:

Δi.sub.s(k+1)=A.sub.0(k)Δi.sub.s(k)+B.sub.0(k)Δu.sub.s(k)  (6)

[0047] In equation (6), Δi.sub.s(k+1)=[Δi.sub.d,k+1, Δ.sub.q,k+1].sup.T represents the predicted current increment calculated from the current increment prediction model, where Δi.sub.d,k+1 and Δi.sub.q,k+1 are d-axis predictive current increment and q-axis predictive current increment, respectively; Δi.sub.s(k) represents the stator current increment between the stator current at time point kT.sub.s and the stator current at time point (k−1)T.sub.s, i.e., Δi.sub.s(k)=i.sub.s(k)−i.sub.s(k−1), and Δi.sub.s(k)=[Δi.sub.d,k, Δi.sub.q,k].sup.T, where Δi.sub.d,k is d-axis stator current increment between the d-axis stator current at time point kT.sub.s and the d-axis stator current at time point (k−1)T.sub.s, and Δi.sub.q,k is q-axis stator current increment between the q-axis stator current at time point kT.sub.s and the q-axis stator current at time point (k−1)T.sub.s; Δu.sub.s(k) represents the stator voltage increment between the stator voltage at time point kT.sub.s and the stator voltage at time point (k−1)T.sub.s, i.e., Δu.sub.s(k)=u.sub.s(k)−u.sub.s(k−1), and Δu.sub.s(k)=[Δu.sub.d,k, Δu.sub.q,k].sup.T where Δu.sub.d,k is d-axis stator voltage increment between the d-axis stator voltage at time point kT.sub.s and the d-axis stator voltage at time point (k−1)T.sub.s, and Δu.sub.q,k is q-axis stator voltage increment between the q-axis stator voltage at time point kT.sub.s and the q-axis stator voltage at time point (k−1)T.sub.s; u.sub.s(k−1) represents the stator voltage at time point (k−1)T.sub.s.

[0048] The control variable in equation (6) is the stator voltage increment Δu.sub.s(k), which indicates that the current increment prediction model can reduce the output voltage error caused by the inverter dead-time effect. Comparing equation (6) with equation (5), it can be seen that the coefficients Δ.sub.0(k) and B.sub.0(k) in equation (5) and equation (6) are equal, but the back electromotive force item is eliminated in equation (6), i.e., the current increment prediction model is independent of the permanent magnet flux linkage, and is only affected by the stator inductance.

[0049] 4. Establish the Cost Function

[0050] The cost function is established by taking the squared error between the preset reference current increment and the predictive current increment as the evaluation criterion. The cost function is applied to evaluate the error of the stator current increment at the end of each control period corresponding to the stator voltage increment. Considering the delay compensation problem of the predictive current control, the cost function is established as:

[00010]$\begin{array}{cc}\underset{\Delta u_s}{\min }J={\left[\Delta i_s^{\mathrm{ref}}-\Delta i_s\left(k+2\right)\right]}^T\left[\Delta i_s^{\mathrm{ref}}-\Delta i_s\left(k+2\right)\right]+\Delta {u_s\left(k+1\right)}^{T}P\Delta u_s\left(k+1\right)& \left(7\right)\end{array}$$\mathrm{Satisfy}:\{\begin{array}{c}.Math.u_s\left(k\right)+\Delta u_s\left(k+1\right).Math.\le U_{\max }\\ .Math.i_s\left(k+1\right)+\Delta i_s\left(k+2\right).Math.\le I_{\max }\end{array}$

[0051] where Δi.sub.s.sup.ref represents the reference current increment; P is the weight factor which is used to determine the importance of voltage increment; U.sub.max and I.sub.max are the maximal voltage and the maximal current of the permanent magnet synchronous motor driving system; the superscript T represents the matrix transpose operation; “Satisfy” represents the constraint conditions; Δi.sub.s(k+2) represents the predicted current increment calculated from the current increment prediction model; Δu.sub.s(k+1) represents the stator voltage increment between the stator voltage at time point (k+1)T.sub.s and the stator voltage at time point kT.sub.s; J is the value of the cost function. The voltage increment item in the cost function is used to reduce the dynamic overshoot of the motor, and to prevent the motor and power switching suffering from voltage surge and current surge.

[0052] 5. Obtain the Optimal Stator Voltage Vector

[0053] Substituting equation (6) into equation (7), the cost function is established as:

[00011]$\begin{array}{cc}\underset{\Delta u_s}{\min }J=\Delta {u_s\left(k\right)}^T\left(B_0^T{B}_0+P\right)\Delta u_s\left(k\right)-2\Delta {u_s\left(k\right)}^T{B}_0^T\left[\Delta i_s^{\mathrm{ref}}-A_0\left(k\right)\Delta i_s\left(k\right)\right]+{\left(\Delta i_s^{\mathrm{ref}}\right)}^T\Delta i_s^{\mathrm{ref}}-2{\left(\Delta i_s^{\mathrm{ref}}\right)}^T{A}_0\left(k\right)\Delta i_s\left(k\right)+\Delta {i_s\left(k\right)}^T{A}_0^T{A}_0\left(k\right)\Delta i_s\left(k\right)& \left(8\right)\end{array}$

[0054] According to the convex optimization theory, the extreme value of the cost function J can be obtained by calculating the partial derivative of formula (8) with respect to Δu.sub.s(k) and make it zero:

[00012]$\begin{array}{cc}\frac{\partial J}{\partial \Delta u\left(k\right)}=0& \left(9\right)\end{array}$

[0055] By solving equation (9), the optimal voltage increment minimizing the value of J can be derived as

Δ.sub.s.sup.opt=[B.sub.0.sup.TB.sub.0+P.sup.−].sup.−1B.sub.0.sup.T[ΔAi.sub.s.sup.ref−A.sub.0Δi.sub.s(k+1)]  (10)

[0056] By superposing the optimal voltage increment on the stator voltage at the present control period, the optimal stator voltage is obtained as:

u.sub.s.sup.opt(k+1)=u.sub.s(k)+Δ.sub.s.sup.opt(k+1)  (11)

[0057] In equation (11), Δu.sub.s.sup.opt(k+1) represents the optimal voltage increment.

[0058] The detailed implementation process of the disclosure is shown in FIG. 1 in the form of a control block diagram. This method improves the current prediction accuracy under high-speed operation and eliminates the use of the permanent magnet flux during implementing the predictive control method of the current increment.

[0059] The feasibility of the proposed method is verified by combining detailed simulation and experimental data and FIG. 2 to FIG. 4.

[0060] To verify the practicability and validity of the proposed predictive current control method based on the current increment prediction model, the simulations and experiments are carried out on a 20-kW permanent-magnet synchronous motor system. The parameters of the tested motor are shown in TABLE I. In the experimental platform, the controller taking the DSP (TMS320F28335) as the core is employed for algorithm implementation, and the dynamometer is an induction motor controlled by S120 produced by Siemens.

TABLE-US-00001 TABLE I parameter Symbols Values Units Number of pole pairs n.sub.p 4 Rated speed n.sub.N 3000 r/min Rated torque T.sub.N 64 Nm Rotor flux linkage Ψ.sub.f0 0.07574 Wb Stator resistance r.sub.s0 0.00114 Ω d-axis stator inductance L.sub.d0 0.2 mH q-axis stator inductance L.sub.q0 0.555 mH

[0061] 1. The Influence of the Dead-Time Effect on the Current Control Performance

[0062] To eliminate the effect of parameter mismatch, this disclosure checks the effect of the dead-time on the current control performances of the predictive current control method based on the conventional current prediction model and predictive current increment control method in this disclosure by simulation. In the simulation, the conventional current prediction model is shown in (3), the control period T.sub.s is set to 200 μs, and the dead-time t.sub.d is set to 3 μs. FIG. 2 shows the simulation waveforms of preset d-axis reference current i.sub.d.sup.ref, d-axis current measured value i.sub.d, preset q-axis current i.sub.q.sup.ref and q-axis current measured value i.sub.q of the predictive current control method based on the conventional current prediction model and predictive control method of the current increment, respectively when the permanent magnet synchronous motor works under different rotating speed n and different load torque T.sub.L.

[0063] From FIG. 2, it can be seen that there is tracking error of between i.sub.d.sup.ref and i.sub.d of the predictive current control method based on the conventional current predictive model when the motor works at 300 r/min with 64 Nm load torque, and the i.sub.d can track i.sub.d.sup.ref in other operation conditions. However, the tracking error between i.sub.q and its reference of the predictive current control method based on the conventional current predictive model is significant in different operation conditions. From FIG. 2, it can also be seen that both i.sub.d and i.sub.q of the predictive current increment control method can track their references well.

[0064] 2. The Analysis of Current Prediction Error Under Parameter Mismatch

[0065] In the practical motor drives, there is the errors between the nominal inductances (L.sub.d0 and L.sub.d0) and the actual inductances (L.sub.d and L.sub.d). Define ΔL.sub.d and ΔL.sub.q as the perturbation values of the inductances. Substituting L.sub.d=L.sub.d0+ΔL.sub.d and L.sub.q=L.sub.q0+ΔL.sub.q into equation (6), the predictive current increment of current increment prediction model considering inductance mismatch is obtained as:

[00013]$\begin{array}{cc}\Delta i_{\mathrm{sp}}\left(k+1\right)=\Delta i_s\left(k+1\right)+e\left(k+1\right)& \left(12\right)\end{array}$$\begin{array}{cc}& \left(13\right)\end{array}$$\{\begin{array}{c}e\left(k+1\right)=\left[\begin{array}{c}{e}_{d,k+1}\\ e_{q,k+1}\end{array}\right]=\Delta A\left(k\right)\left[\begin{array}{c}\Delta i_{d,k}\\ \Delta i_{q,k}\end{array}\right]+\Delta B\left(k\right)\left[\begin{array}{c}\Delta u_{d,k}\\ \Delta u_{q,k}\end{array}\right]\\ \Delta A\left(k\right)=\left[\begin{array}{cc}0& \frac{\Delta L_q{L}_{d0}-\Delta L_d{L}_{q0}}{L_{d0}\left(L_{d0}+\Delta L_d\right)}\sin \left({\omega }_{r,k}{T}_s\right)\\ \frac{\Delta L_q{L}_{d0}-\Delta L_d{L}_{q0}}{L_{q0}\left(L_{q0}+\Delta L_q\right)}\sin \left({\omega }_{r,k}{T}_s\right)& 0\end{array}\right]\\ \Delta B\left(k\right)=\left[\begin{array}{cc}-\frac{T_s\Delta L_d}{L_{d0}\left(L_{d0}+\Delta L_d\right)}\cos \left({\omega }_{r,k}{T}_s\right)& -\frac{T_s\Delta L_d}{L_{d0}\left(L_{d0}+\Delta L_d\right)}\sin \left({\omega }_{r,k}{T}_s\right)\\ \frac{T_s\Delta L_q}{L_{q0}\left(L_{q0}+\Delta L_q\right)}\sin \left({\omega }_{r,k}{T}_s\right)& -\frac{T_s\Delta L_q}{L_{q0}\left(L_{q0}+\Delta L_q\right)}\cos \left({\omega }_{r,k}{T}_s\right)\end{array}\right]\end{array}$

[0066] In equation (12), Δi.sub.sp(k+1) represents the predicted current increment considering inductance mismatch; e(k+1)=[e.sub.d,k+1, e.sub.q,k+1].sup.T represents the current prediction error vector caused by inductance mismatch, where e.sub.d,k+1 and e.sub.q,k+1 are the d-axis current prediction error and q-axis current prediction error, respectively; Δi.sub.s(k+1) has been shown in equation (6).

[0067] In equation (13), ΔL.sub.d represents the error between d-axis actual inductance La and d-axis nominal inductance L.sub.d0; ΔL.sub.q represents the error between q-axis actual inductance L.sub.q and q-axis nominal inductance L.sub.q0; ΔA(k) represents the coefficient matrix of the stator current increment item; ΔB(k) represents the coefficient matrix of the stator voltage increment item.

[0068] From equation (13), it can be seen that the amplitude and sign of current prediction error (e.sub.d,k+1 and e.sub.q,k+1) are affected by the state variables (Δi.sub.d,k and Δi.sub.q,k), the control variables (Δu.sub.d,k and Δu.sub.q,k) and the angular velocity ω.sub.r. To intuitively illustrate the influence of inductance mismatch on the current prediction error of the current increment prediction model, waveforms of e.sub.d,k+1 and e.sub.d,k+1 are shown in FIG. 3 for the current increment prediction model under different operation conditions (motor rotational speed n and load torque T.sub.L) when L.sub.d and L.sub.q increase up 20% and 50%, respectively.

[0069] For the current increment prediction model, e.sub.d,k+1 and e.sub.q,k+1 always fluctuate around zero as shown in FIG. 3. Such current prediction errors mainly impact the current ripples, but have almost no effect on the current tracking error. In addition, the larger the fluctuation range of e.sub.d,k+1 (or e.sub.q,k+1) is, the larger the effect of inductance mismatch on the d-axis current ripple (or q-axis current ripple) is. From FIG. 3, it can also be seen the fluctuation ranges of e.sub.d,k+1 and e.sub.q,k+1 change very little under different motor operation conditions. When the load torque is light, the amplitudes of d-axis current and q-axis current are not large enough compared with the fluctuation ranges of e.sub.d,k+1 and e.sub.q,k+1, so the effects of inductance mismatches on the d-axis current ripple and q-axis current ripple are great. However, as the torque increases, the amplitudes of d- and q-axis currents become large, and the fluctuation ranges of e.sub.d and e.sub.q become small compared with amplitudes of d- and q-axis currents. So the effects of inductance mismatch on the current ripples reduce. The changes of the inductances under the conditions with different load torque in above analysis are set the same (ΔL.sub.d=20% L.sub.d0 and ΔL.sub.q=50% L.sub.q0). However, there are a positive correlation between the inductances and stator currents. In other words, the changes of the inductances in light load are much smaller than the setting values. Therefore, the effect of the inductance mismatch on the current ripples of predictive current increment control is not great in the practical motor operating.

[0070] In summary, inductance mismatch has the impact on the current ripple of the predictive control method of the current increment, but the effect of the inductances mismatch is small in the practical motor operating. In addition, the inductance mismatch has almost no effect on the current tracking error.

[0071] 3. The Comparison of Steady-State Performance

[0072] This disclosure compares the steady-state performance of predictive current increment control with the steady-state performance of predictive current control based on conventional current predictive model in a 20 kW PMSM drives. The parameter of the tested PMSM is shown in TABLE I. In the experiment, the motor works at 300 r/min and 7500 r/min, and the output power of the motor is 20 kW.

[0073] FIG. 4 shows the steady-state waveforms of the presuppose d-axis reference current i.sub.d.sup.ref, the d-axis current i.sub.d, the presuppose q-axis reference current i.sub.q.sup.ref the q-axis current i.sub.q and a-phase current i.sub.a of the predictive current control method based on the conventional current prediction model and current increment prediction model, respectively under different operation condition. From the experimental results, it can be seen that the current ripples of the predictive current increment control do not increase compared with the predictive current control method based on the conventional current predictive model when the motor works at 300 r/min, although current increment prediction model ignores the stator resistance. However, both the d-axis current ripple and q-axis current ripple of predictive current increment control are significantly less than the d-axis current ripple and q-axis current ripple of predictive current control method based on the conventional current prediction model when the motor works at 7500 r/min. In addition, it can be seen that i.sub.d and i.sub.q of predictive current increment control can track their references well, but there are obvious d-axis current tracking error and q-axis current tracking error of predictive current control method based on the conventional current prediction model caused by the dead-time affect and parameter mismatch.

[0074] The disclosure does not limit the type of each device except for a special description. so long as the device can complete the above functions.

[0075] Technical personnel in this field can understand that the attached figure is only a schematic diagram, and the serial number of the above disclosure implementation cases is only for description, which does not represent the advantages and disadvantages of the implementation cases.

[0076] The above is only a better implementation case of the invention, which is not used to limit the invention. Any modification, equivalent replacement, improvement, etc. within the spirit and principle of the invention should be included in the protection scope of the invention.