Robust permanent magnet synchronous motor dead-beat model prediction method

Abstract

The invention relates to a robust dead-beat model prediction method for a permanent magnet synchronous motor, which relates to the field of motor control and comprises the following steps of: sampling stator current and stator voltage of the permanent magnet synchronous motor, and outputting the stator current and the stator voltage to an MRAS observer; the MRAS observer calculates a rotating speed estimation value and a resistance estimation value according to the obtained parameters; a dead-beat prediction controller is adopted to receive corresponding parameters of the motor and parameters output by the MRAS observer, and an ideal voltage vector is calculated in each sampling period and is output to a discrete space vector modulation module; and a standardized triangular area where the ideal voltage vector is located is determined through processing of the discrete space vector modulation module, three vertex vectors of the standardized triangular area are used as candidate voltage vectors respectively, and an optimal control vector is selected to be output to the inverter. The method has the advantages of high efficiency, low pulse output and small calculation amount, can realize high-precision voltage prediction and control under the condition of parameter uncertainty, and improves the robustness.

Description

Technical Field

[0001] This invention relates to the field of motor control, and more specifically to a robust method for predicting deadbeat models of permanent magnet synchronous motors. Background Technology

[0002] In traditional permanent magnet synchronous motor control strategies, finite control set model predictive control (FCIMM) is used. FCS-MPC This method has attracted widespread attention due to its rapid dynamic response and excellent steady-state performance. However, its inherent drawbacks include high torque ripple, large computational load, and its control performance is easily affected by the accuracy of motor model parameters.

[0003] Another common control scheme is single-vector model predictive torque control (MPC). MPTC This scheme applies only one voltage vector to the system in each sampling period. Although it can maintain a fast dynamic response speed, it will result in large torque ripple and current harmonics, making it difficult for the system to achieve ideal control accuracy and stability during steady-state operation.

[0004] Furthermore, deadbeat control theoretically possesses superior dynamic response characteristics and stronger steady-state performance; however, this method is more sensitive to changes in motor parameters, and its system robustness is inferior to traditional finite set model predictive control. FCS-MPC )method.

[0005] In view of this, this application has conducted in-depth research on this basis, resulting in this case. Summary of the Invention

[0006] The purpose of this invention is to provide a robust deadbeat model prediction method for permanent magnet synchronous motors, which can achieve high efficiency, low ripple output, and low computational load. At the same time, it can realize high-precision voltage prediction and control under parameter uncertainty conditions, while improving robustness.

[0007] To achieve the above objectives, the solution of the present invention is: a robust method for predicting deadbeat models of permanent magnet synchronous motors, comprising the following steps: Step 1, sampling the permanent magnet synchronous motor in... The stator current and stator voltage at time t are denoted as . and And output it to the MRAS observer; Step 2: The MRAS observer calculates the rotor speed estimate and stator resistance estimate for each sampling period based on the acquired parameters; Step 3: Receive the reference torque and flux linkage amplitude of the permanent magnet synchronous motor using a deadbeat predictive controller. The stator current, the rotor speed estimation value and the stator resistance estimation value output by the MRAS observer are determined at each sampling period, and an ideal voltage vector is calculated and output to a discrete space vector modulation module; Step 4, in the discrete space vector modulation module, a voltage plane in the stationary reference frame is divided into a plurality of normalized triangular regions according to the direction of the mother boundary line, a normalized triangular region is determined by calculating the geometric distance between the ideal voltage vector and three selected mother boundary lines, and the voltage vectors corresponding to the three vertices of the normalized triangular region are respectively taken as candidate voltage vectors, and then the optimal control vector is selected according to the minimization criterion; α - β The voltage plane in the stationary reference frame is divided into a plurality of normalized triangular regions according to the direction of the mother boundary line, a normalized triangular region is determined by calculating the geometric distance between the ideal voltage vector and three selected mother boundary lines, and the voltage vectors corresponding to the three vertices of the normalized triangular region are respectively taken as candidate voltage vectors, and then the optimal control vector is selected according to the minimization criterion; Step 5, the switching duration corresponding to the optimal control vector is calculated, and a signal is output to an inverter; PWM The inverter outputs the stator voltage and the stator current at the moment to the MRAS observer, and then step 2 is repeated. The inverter outputs the stator voltage and the stator current at the moment to the MRAS observer, and then step 2 is repeated. The inverter outputs the stator voltage and the stator current at the moment to the MRAS observer, and then step 2 is repeated.

[0008] In step 2, the following steps are included; Step 2-1, a reference model is established, and a reference flux linkage value at the moment is calculated according to the reference model and transmitted to an error calculation module; wherein the reference model is , wherein represents the flux linkage value at the moment calculated by the reference model, represents the stator voltage vector at the moment, represents the nominal resistance, represents the stator current vector at the moment, represents the sampling period of the control system; Step 2-2, an adaptive prediction model is established, and a stator flux linkage estimation value at the moment and a stator current estimation value at the moment are calculated according to the adaptive prediction model and transmitted to the error calculation module described below, wherein the adaptive prediction model is , , , ; , wherein respectively represent the stator flux linkage estimation value at the moment and the stator flux linkage estimation value at the moment, represents the stator resistance estimation value at the moment, represents the stator resistance estimation value at the moment, represents the stator resistance estimation value at the moment, represents the stator resistance estimation value at the moment, ​​​​This represents the rotor speed estimate output by the MRAS observer during the current sampling period. express Stator current estimate at time [time] Represents the imaginary unit. This represents the mapping function used to calculate the stator current estimate from the stator flux estimate and the rotor speed estimate; Steps 2-3: Obtain the parameter error using the error calculation formula, which is as follows: , ; In the formula, express Momentary current error, express Momentary flux error express The stator reference current value at that moment; Steps 2-4, using Pi The adaptive adjustment model adaptively adjusts the parameter error and then adjusts the adaptively adjusted parameter. Rotor speed estimate at time and The estimated stator resistance values ​​at each time step are fed back to the adaptive prediction model for state prediction in the next sampling period. Pi The adaptive adjustment model is as follows: , ; In the formula, , This represents the adjustment coefficient. Indicates passage Pi The adaptive adjustment model was calculated to obtain Estimated rotor speed at any given time , This represents the integral adjustment coefficient. Indicates a discrete-time index; Steps 2-5: The MRAS observer outputs the updated current sampling period. and , This represents the estimated rotor speed. This represents the estimated value of the stator resistance.

[0009] In step 3, the discretized mathematical model of the permanent magnet synchronous motor is established as follows: In the formula, Indicates the first Stator current vector within each sampling period, This represents the ideal voltage vector. This represents the output of the MRAS observer. Estimated stator resistance at time [time] This indicates the value estimated by the MRAS observer based on the stator flux linkage during the current sampling period. and rotor speed estimate The calculated back electromotive force, Indicates stator inductance, Indicates the sampling period of the control system; Next, in the The stator current vector is obtained through the current sampling device within each sampling period. A prediction model is established based on the stator current vector. To achieve deadbeat tracking of the reference current by the stator current vector in the next sampling period, the predicted current value in the discretized mathematical model is set to... Compared with reference current value Equal to each other, and substituted into the discretized mathematical model, the ideal voltage vector is obtained as follows: .

[0010] In step 4, α - β In the coordinate system, the voltage plane formed by the inverter output voltage is hexagonal in shape. The voltage plane includes three sets of parallel parent boundary lines, with three non-parallel parent boundary lines serving as directional baselines.

[0011] Includes the following steps: Step 4-1: Determine the standardized triangular region. Let the coordinates of the endpoints of the ideal voltage vector be ( ). u α , u β The equations for the three non-parallel parent boundary lines in the voltage plane are established as follows: , , ; In the formula, , , These represent the three non-parallel parent boundary lines of the voltage plane. u α , u β They represent respectively Clarke After coordinate transformation α - β Two-phase equivalent voltage components in a coordinate system This indicates the DC bus voltage of the inverter; Then, the perpendicular distance from the end point of the ideal voltage vector to the three mother boundary lines is calculated, and the calculation formula is , , ; Next, the normalized distance integers are obtained by quantizing respectively , , , and the distance integers are , wherein , denotes the normalized coefficient, denotes the number of equal parts in which the hexagon is divided in one sampling period; Then, the standardized triangular region is determined, and the number is , wherein , , denote the normalized integer coordinates of the three sides respectively; Step 4-2, determine three candidate voltage vectors; First, the reference point coordinates of the standardized triangular region are calculated , , and the calculation formula is , , , wherein denotes the coordinate of the point on the α axis, denotes the coordinate of the point on the β axis, denotes the basic voltage step; Then, according to the reference point coordinates , and the basic voltage step Δ, the candidate voltage vectors corresponding to the three vertices of the standardized triangular region are calculated, and the complex form is as follows , , , denotes the imaginary unit; Step 4-3, the future response of each candidate voltage vector is predicted and calculated by using a cost function, and the cost function is , , wherein denotes the reference electromagnetic torque, denotes the reference stator flux amplitude; denotes the next sampling period electromagnetic torque estimate value obtained under the condition of applying the candidate voltage vector, respectively represent the next sampling period stator flux linkage estimation values obtained under the condition of applying the candidate voltage vector, 、 、 、 respectively represent the weight coefficients of torque, flux linkage, switching loss and current constraint, represent the number of switching actions required when the inverter switches from the current switching state to the switching state corresponding to the candidate voltage vector V ( i ) corresponding to the candidate voltage vector, which is used to measure the switching loss; represent the current constraint penalty term; represent the index number of the candidate voltage vector, which is used to distinguish different candidate voltage vectors, represent the first candidate voltage vector; Step 4-4, establish an optimal control vector selection formula, select the candidate voltage vector with the minimum comprehensive performance index as the optimal control vector, and the optimal control vector selection formula is, .

[0012] The inverter is a three-phase two-level inverter, and the inverter has six effective switching states.

[0013] After step 4-2, according to the sign correction according to the parity relationship of and + , if the parity of and + is the same, the negative sign is taken when calculating the vertex coordinates corresponding to the candidate voltage vector, if the parity of and + is different, the positive sign is taken when calculating the vertex coordinates corresponding to the candidate voltage vector, and the negative sign and the positive sign are used to distinguish the directionality of adjacent regions in the standardization triangular region numbering process.

[0014] After adopting the above method, the present invention has the following beneficial effects: 1. The dual-layer control structure is adopted: the first layer realizes the prediction and control of voltage under the condition of parameter uncertainty through the deadbeat prediction model that does not rely on velocity estimation; the second layer adopts an improved discrete space vector modulation module, which transforms the traditional method of enumerating multiple virtual voltage vectors into only evaluating three candidate voltage vectors, thereby completing the synthesis of predicted voltage in a finite set of candidate voltage vectors, reducing the amount of computation, and helping to reduce output ripple. The two parts are functionally independent and cooperate with each other during operation. The former provides an ideal voltage reference, and the latter ensures the continuity of physical realization, thereby improving the real-time performance and robustness of the system.

[0015] 2. This invention employs an MRAS observer based on Lyapunov theory, which can simultaneously acquire real-time stator resistance estimates and rotor speed estimates, thereby effectively reducing drift errors and ensuring system stability. Attached Figure Description

[0016] Fig. 1 This is a system flowchart of the present invention.

[0017] Fig. 2 This is a control block diagram of the MRAS observer in this invention.

[0018] Fig. 3 This is the voltage plane containing the ideal voltage vector endpoints in this invention. Detailed Implementation

[0019] To further explain the technical solution of the present invention, the present invention will be described in detail below through specific embodiments.

[0020] A robust deadbeat-free model prediction method for permanent magnet synchronous motors (PMSMs) is proposed. This method is based on common control systems used in PMSM control, such as... Fig. 1 As shown, the control system includes an inverter, which is a conventional two-level three-phase inverter. A conventional two-level three-phase inverter has eight switching states, and in each of these eight switching states, it generates six non-zero voltage vectors. U 1~ U 6) and two zero-voltage vectors. These non-zero voltage vectors are the basic voltage vectors. In other words, the two-level three-phase inverter has six effective switching states. These six effective switching states are all the inverter's conventional switching states, and therefore will not be described further. In this embodiment, the six basic voltage vectors generated by the three-phase two-level inverter under six different switching states are as follows: Fig. 3As shown, they are 60° apart in space, and each end point together constitutes a regular hexagon, which represents the limit range of the inverter output voltage, and the space formed is the voltage plane for subsequent space vector modulation.

[0021] In this embodiment, as shown in the figure, the deadbeat model prediction method of the permanent magnet synchronous motor includes the following steps. Fig. 1-3

[0022] Step 1, sample the permanent magnet synchronous motor to obtain the current stator current and current stator voltage, in this embodiment, define the current sampling time as , the next sampling time as , so the current stator current and current stator voltage correspond to and respectively, and then output the sampled and to the MRAS observer.

[0023] Further, the permanent magnet synchronous motor is provided with a current sampling device (such as a current sensor) and a voltage sensor, the current sensor is used to sample the current stator current, and the voltage sensor is used to sample the current stator voltage. It is worth mentioning that the current sensor and the voltage sensor are configured on the conventional permanent magnet synchronous motor, which is a conventional operation in the motor field, so it will not be described.

[0024] Step 2, the MRAS observer calculates the rotor speed estimation value ( ) and the stator resistance estimation value ( ) according to the obtained and .

[0025] Specifically, step 2-1, establish a reference model, calculate the reference flux value at according to the reference model, and transmit it to the error calculation module.

[0026] In detail, the reference model is based on the motor stator voltage equation, which is: , wherein represents the stator voltage vector, represents the stator current vector, represents the stator resistance, represents the stator flux vector, represents the rotor speed, represents the imaginary unit, which is used to express the voltage vector in the α - β coordinate system as a complex number.

[0027] ​Due to For the speed coupling term, it will cause the flux linkage estimation sensitive to the rotor speed variation, therefore, in order to eliminate the above influence, the components related to the speed and its high frequency disturbance are grouped into the same function , which is defined as: , where, represents the high frequency and nonlinear error term. Thus, the reference equation independent of the speed can be obtained as: , where, represents the stator voltage vector at the current time, represents the nominal resistance, represents the stator current vector at the current time, represents the reference flux linkage value at the current time, represents the sampling period of the control system.

[0028] And since the influence of the above function can be neglected in modeling, the above reference model is: . Where, represents the sampling period, represents the calculated reference flux linkage value at the next time (i.e. time).

[0029] It is worth mentioning that the above reference model relies on the measurable voltage and current, and does not contain unknown speed terms. With this reference model, the reference flux linkage value at the time can be calculated, and then combined with the flux equation to convert the reference current value at the next time (i.e. time) . Finally, the obtained reference flux linkage value and reference current value are transmitted to the error calculation module Step 2-2, the adaptive prediction model is used to calculate the stator flux linkage estimation value at the time and the stator current estimation value at the time , and transmitted to the error calculation module.

[0030] In detail, the adaptive prediction model is an online observer based on the mathematical model of the motor, which aims to use the estimated parameters at the current time (i.e. time) to predict the flux linkage and current estimation values of the motor at the next time (i.e. time), which is based on the state prediction of the discretized voltage equation, as follows.

[0031] The flux linkage prediction equation is established as: , where, ​They represent Time and The estimated value of the stator flux linkage at time t. express Estimated stator resistance at time [time] This represents the rotor speed estimate output by the MRAS observer during the current sampling period.

[0032] The current prediction equation is established as follows: In the formula, express Stator current estimate at time, function f This represents the flux linkage-current relationship for solving the current from the flux linkage and rotor speed, i.e., the mapping function for calculating the stator current estimate from the stator flux linkage estimate and the rotor speed estimate.

[0033] Steps 2-3: In the error calculation module, the parameters transmitted from the reference model and the adaptive prediction model are processed using the error calculation formula to obtain the parameter acquisition error. This parameter error is then sent to... Pi In the adaptive adjustment model.

[0034] To elaborate, the error calculation formula is as follows: , ; In the formula, express Momentary current error, express Momentary flux error express The stator reference current value at a given time.

[0035] Steps 2-4: Employing Lyapunov functions as a basis Pi The adaptive adjustment model adaptively adjusts the parameter error and then adjusts the result accordingly. Rotor speed estimate at time and The estimated stator resistance values ​​at each time point are fed back to the adaptive prediction model for state prediction in the next sampling period.

[0036] To elaborate, let's first define a Lyapunov candidate function to characterize the energy change of the system. This Lyapunov candidate function is: In the formula, Let Lyapunov function be denoted as and let .... Total error energy at any given time; express flux linkage error at time , express Current error at any given time.

[0037] To ensure the monotonic decrease of Lyapunov candidate function, it is required to make its difference satisfy the following formula: , wherein represents the change of Lyapunov function, and thus the basic form of the parameter update rate is: , ; , wherein , represents an adjustment coefficient, represents the rotor speed estimation value at the moment calculated by the adaptive adjustment model. Pi

[0038] To improve the steady-state accuracy and response speed, an integral compensation term is introduced on the basis of the basic form of the parameter update rate to form the formula update structure Pi , and the formula update structure is the adaptive adjustment model Pi , and the adaptive adjustment model is: Pi Pi , ; , wherein , represents an integral adjustment coefficient, and the summation term represents the error accumulation from the initial moment to the current moment (i.e., the moment ). represents a discrete time index, which is used to represent the historical moment from the initial sampling period to the kth sampling period.

[0039] Step 2-5, the MRAS observer calculates the latest parameter estimation value according to the adaptive adjustment model above, and outputs the updated rotor speed estimation value Pi and the stator resistance estimation value at the current sampling period, which are transmitted to the subsequent module for deadbeat predictive control.

[0040] Step 3, the deadbeat predictive controller receives the reference torque of the permanent magnet synchronous motor , the flux amplitude , the stator current at the moment, and the real-time estimated motor rotor speed and the real-time estimated stator resistance output by the MRAS observer, and calculates the ideal voltage vector in each sampling period and outputs it to the discrete space vector modulation module.

[0041] ​​​​​To elaborate, the discretized mathematical model of the permanent magnet synchronous motor is established as follows: ; Will and Forced equality yields the ideal voltage vector. The formula for the ideal voltage vector is: .

[0042] In the formula, Represents the ideal voltage vector. This represents the value estimated by the MRAS observer. Estimated stator resistance at time [time] Indicates stator inductance, This represents the back electromotive force (EMF) obtained by the MRAS observer based on the flux linkage estimate and the rotor speed estimate. This back EMF needs to be combined with the rotor electric angular velocity obtained from the MRAS observation. The calculation is performed using the following formula: , In the formula, This is the estimated value of the stator flux linkage. j It is the imaginary unit.

[0043] Step 4: Forming a discrete space vector modulation module in the discrete space vector modulation module. α - β The voltage plane in the coordinate system is divided into multiple standardized triangular regions according to the direction of the parent boundary line. The standardized triangular regions are determined by calculating the geometric distance between the ideal voltage vector and the three selected parent boundary lines. The voltage vectors corresponding to the three vertices of the standardized triangular regions are used as candidate voltage vectors. Then, the optimal control vector is selected according to the minimization criterion.

[0044] To elaborate, step 4-1 is to determine the standardized triangular region.

[0045] Because in α - β In the coordinate system, the voltage plane of the inverter output voltage is distributed in a regular hexagon. Its geometric boundary contains three pairs of parallel line segments, that is, the voltage plane includes three sets of parallel parent boundary lines. Therefore, only three non-parallel parent boundary lines need to be defined as directional baselines, and the boundary structure of the entire voltage limit space can be generated by translating along these three parent boundary lines. The regular hexagon can be represented by three non-parallel parent boundary lines. In this embodiment, the three non-parallel parent boundary lines are respectively... , , Their equations are: , , ; In the formula, u α , u β They represent respectively Clarke After coordinate transformation α - β Two-phase equivalent voltage components in a coordinate system This indicates the DC bus voltage of the inverter.

[0046] Let the coordinates of the endpoints of the ideal voltage vector (i.e., DB-VV) be ( u α , u β The perpendicular distances from the endpoints of the ideal voltage vector to the three selected parent boundary lines are denoted as follows: d 1, d 2, d 3. The calculation formula is as follows: , , Thus, through three sets of vertical distances ( , , Together, they determined the position of the ideal voltage vector in the regular hexagonal grid formed by the periodic translation of the parent boundary line.

[0047] To facilitate digital indexing and fast addressing, and to facilitate digital implementation, , , The normalized distance integers are obtained by quantization respectively, and these distance integers are respectively , , and with The unique identifier of the standardized triangle region containing the ideal voltage vector is represented in the form of , as follows.

[0048] The quantized distance integer is: , ,in, , This represents the normalization coefficient, used to map continuous geometric distances to integer indices; This represents the number of equal divisions of the hexagon within a sampling period.

[0049] Then, the unique number of the above-mentioned standardized triangular region is: , h m Represents the normalized integer distance, and represents the ideal voltage vector at the _th ...m The quantization level number along the boundary line is a key parameter for achieving rapid geometric positioning. Among these parameters, the triplet ( , , It is a location index used to uniquely identify the standardized triangular region. 。

[0050] It is worth mentioning that, through the geometric positioning method, it is possible to... α - β In the coordinate system, directly based on the quantized distance integer , , Precisely locate the standardized triangular region containing the ideal voltage vector (DB-VV), a location process based on the parent boundary line. , , The periodic translational characteristics divide the regular hexagonal voltage limit region into several standardized small triangular regions, i.e., standardized triangular regions. The position index of each standardized triangular region can be determined by three quantized index integers ( , , The only certainty is that among them These represent the relative distances of the ideal voltage vector to the three parent boundary lines, and also reflect its relative geometric position within the triangle.

[0051] Step 4-2: Determine three candidate voltage vectors.

[0052] Based on step 4-1, the standardized triangular region containing the ideal voltage vector was determined, and the coordinates of the reference point and the basic voltage step size of this region were calculated. The calculation formula is as follows: , , In the formula, This indicates that the point is at α Coordinates on the axis This indicates that the point is at β Coordinates on the axis; This represents the basic voltage step size, corresponding to the minimum voltage resolution after discrete space vector modulation (DSVM) subdivision.

[0053] Then, based on the reference point coordinates of this standardized triangular region ( a , b Given the basic voltage step size, calculate the candidate voltage vectors corresponding to the three vertices of the normalized triangular region. Their complex form is as follows: , , , It represents the imaginary unit. It's worth mentioning that... , , The geometric relationship forms an equilateral triangular unit, whose interior points can be described by a linear combination of the three vertices, which is used to approximate a continuous ideal voltage vector in discrete space.

[0054] Furthermore, when the sampling period is divided into even numbers, to ensure that the geometric direction of the candidate voltage vector is consistent with the numbering direction, it is necessary to... and + The parity correction sign is as follows: if and + If the parity of the candidate voltage vectors is the same, then a negative sign is taken when calculating the vertex coordinates of the candidate voltage vectors. Otherwise, use "+". The positive and negative signs here are used to distinguish the directionality of adjacent areas during triangle numbering, preventing reverse deviations in vector indexing during calculations, and have no effect on voltage amplitude.

[0055] In this way, by following the steps above, the three candidate voltage vectors closest to the ideal voltage can be quickly determined without traversing all virtual voltage vectors.

[0056] Step 4-3: Calculate the future response of each candidate voltage vector using a cost function, which is: , In the formula, Indicates the reference electromagnetic torque. Indicates the reference stator flux linkage amplitude; This represents the electromagnetic torque estimate obtained in the next sampling period under the applied candidate voltage vector. These represent the stator flux linkage estimates obtained in the next sampling period under the applied candidate voltage vector; , , , These represent the weighting coefficients for torque, flux linkage, switching losses, and current constraints, respectively. This indicates that the inverter is switching from its current switching state to the candidate voltage vector. The number of switching actions in the corresponding switching state is used to suppress excessively high switching frequency; This represents the current constraint penalty term, which takes a larger value when the estimated stator current exceeds the set limit, and zero otherwise, and is used to prevent overcurrent.

[0057] Step 4-4: Establish the optimal control vector selection formula, and select the candidate voltage vector with the smallest comprehensive performance index as the optimal control vector. The optimal control vector selection formula is as follows: .

[0058] It should be noted that the optimal control vector , which is the actual output voltage of the inverter during the current sampling period.

[0059] Thus, by employing the above method, while ensuring the accuracy of torque and flux regulation, switching frequency fluctuations and overcurrent risks can be effectively suppressed, achieving the desired voltage vector. High-precision approximation and robust control.

[0060] Step 5: The discrete space vector modulation module calculates the switching duration corresponding to the optimal control vector and outputs it. PWM Signal (i.e.) Sa , Sb , Sc (This is given to the inverter, which is based on...) PWM The signal is switched on and off, and output is also performed. Stator voltage at time and The stator current is fed into the MRAS observer at a certain time, and then step 2 is repeated.

[0061] This invention presents a robust deadbeat model prediction method for permanent magnet synchronous motors (PMSMs). It organically integrates a predicted MRAS observer with algebraic geometric discrete space vector modulation into a deadbeat model predictive control structure, forming a self-correcting predictive control method with low computational complexity, high robustness, and no speed measurement requirement. Specifically, it employs an MRAS observer based on Lyapunov theory to simultaneously predict stator resistance and rotor speed, effectively reducing drift error and ensuring stability. Furthermore, it constructs a deadbeat predictive controller that eliminates the speed term, directly avoiding control drift caused by speed estimation errors. Moreover, it uses an improved discrete space vector modulation method to reduce the number of candidate voltage vectors to three, significantly lowering computational complexity.

[0062] The above description is only a preferred embodiment of this invention. Any equivalent changes and modifications made within the scope of the claims of this invention shall fall within the scope of the claims of this invention.

Claims

1. A robust method for predicting deadbeat models of permanent magnet synchronous motors, characterized in that, The steps include: Step 1, sampling the permanent magnet synchronous motor. The stator current and stator voltage at time t are denoted as . and And output it to the MRAS observer; Step 2: The MRAS observer calculates the rotor speed estimate and stator resistance estimate for each sampling period based on the acquired parameters; Step 3: Receive the reference torque and flux linkage amplitude of the permanent magnet synchronous motor using a deadbeat predictive controller. The stator current at each time point, along with the rotor speed estimate and stator resistance estimate output by the MRAS observer, are used to calculate the ideal voltage vector in each sampling period and output it to the discrete space vector modulation module. Step 4: In the discrete space vector modulation module, a process is formed in... α - β The voltage plane in the coordinate system is divided into multiple standardized triangular regions according to the direction of the parent boundary line. The standardized triangular regions are determined by calculating the geometric distance between the ideal voltage vector and the three selected parent boundary lines. The voltage vectors corresponding to the three vertices of the standardized triangular regions are used as candidate voltage vectors. Then, the optimal control vector is selected according to the minimization criterion. Step 5: Calculate the switching duration corresponding to the optimal control vector and output it. PWM The signal is sent to the inverter, and the inverter outputs... Stator voltage at time and The stator current is fed into the MRAS observer at a certain time, and then step 2 is repeated.

2. The robust prediction method for a deadbeat-free model of a permanent magnet synchronous motor according to claim 1, characterized in that: Step 2 includes the following steps; Step 2-1: Establish a reference model and calculate based on the reference model. The reference flux linkage value is constantly monitored and transmitted to the error calculation module; wherein, the reference model is... In the formula, express The flux linkage value calculated from the reference model at each time step. express Stator voltage vector at time , Indicates nominal resistance. express Stator current vector at any given moment Indicates the sampling period of the control system; Step 2-2: Establish an adaptive prediction model and calculate the results based on the adaptive prediction model. Stator flux linkage estimate at time and The estimated stator current value at time t is obtained and transmitted to the error calculation module described below, wherein the adaptive prediction model is: ， ； In the formula, They represent Time and The estimated value of the stator flux linkage at time t. express Estimated stator resistance at time [time] This represents the rotor speed estimate output by the MRAS observer during the current sampling period. express Stator current estimate at time [time] Represents the imaginary unit. This represents the mapping function used to calculate the stator current estimate from the stator flux estimate and the rotor speed estimate; Steps 2-3: Obtain the parameter error using the error calculation formula, which is as follows: ， ； In the formula, express Momentary current error, express Momentary flux error express The stator reference current value at that moment; Steps 2-4, using Pi The adaptive adjustment model adaptively adjusts the parameter error and then adjusts the adaptively adjusted parameter. Rotor speed estimate at time and The estimated stator resistance values ​​at each time step are fed back to the adaptive prediction model for state prediction in the next sampling period. Pi The adaptive adjustment model is as follows: ， ； In the formula, , This represents the adjustment coefficient. Indicates passage Pi The adaptive adjustment model was calculated to obtain Estimated rotor speed at any given time , This represents the integral adjustment coefficient. Indicates a discrete-time index; Steps 2-5: The MRAS observer outputs the updated current sampling period. and , This represents the estimated rotor speed. This represents the estimated value of the stator resistance.

3. The robust prediction method for a deadbeat-free model of a permanent magnet synchronous motor according to claim 1, characterized in that: In step 3, the discretized mathematical model of the permanent magnet synchronous motor is established as follows: In the formula, Indicates the first Stator current vector within each sampling period, This represents the ideal voltage vector. This represents the output of the MRAS observer. Estimated stator resistance at time [time] This indicates the value estimated by the MRAS observer based on the stator flux linkage during the current sampling period. and rotor speed estimate The calculated back electromotive force, Indicates stator inductance, Indicates the sampling period of the control system; Next, in the The stator current vector is obtained through the current sampling device within each sampling period. A prediction model is established based on the stator current vector. To achieve deadbeat tracking of the reference current by the stator current vector in the next sampling period, the predicted current value in the discretized mathematical model is set to... Compared with reference current value Equal to each other, and substituted into the discretized mathematical model, the ideal voltage vector is obtained as follows: 。 4. A robust method for predicting deadbeat models of permanent magnet synchronous motors according to claim 1, 2, or 3, characterized in that: In step 4, α - β In the coordinate system, the voltage plane formed by the inverter output voltage is hexagonal in shape. The voltage plane includes three sets of parallel parent boundary lines, with three non-parallel parent boundary lines serving as directional baselines.

5. The robust prediction method for a deadbeat-free model of a permanent magnet synchronous motor according to claim 4, characterized in that: Includes the following steps: Step 4-1: Determine the standardized triangular region. Let the coordinates of the endpoints of the ideal voltage vector be ( ). u α , u β The equations for the three non-parallel parent boundary lines in the voltage plane are established as follows: ， ， ； In the formula, , , These represent the three non-parallel parent boundary lines of the voltage plane. u α , u β They represent respectively Clarke After coordinate transformation α - β Two-phase equivalent voltage components in a coordinate system This indicates the DC bus voltage of the inverter; Then, the perpendicular distance from the endpoint of the ideal voltage vector to the three parent boundary lines is calculated using the following formula: , , ; Next, , , The normalized distance integers are obtained by quantization, and the distance integers are... ,in, , Represents the normalization coefficient. This represents the number of equal divisions of the hexagon within a sampling period; Then, the standardized triangular region is determined, and its number is... ,in, , , These represent the normalized integer coordinates of the three sides, respectively. Step 4-2: Determine three candidate voltage vectors; First, calculate the coordinates of the reference point in the standardized triangular region. , The calculation formula is as follows: , , In the formula, This indicates that the point is at α Coordinates on the axis This indicates that the point is at β Coordinates on the axis Indicates the basic voltage step size; Then, based on the coordinates of the reference point ( , Given the basic voltage step size Δ, calculate the candidate voltage vectors corresponding to the three vertices of the standardized triangular region. The complex form of these vectors is as follows: , , , Represents the imaginary unit; Step 4-3: Calculate the future response of each candidate voltage vector using a cost function, which is: ， In the formula, Indicates the reference electromagnetic torque. Indicates the reference stator flux linkage amplitude; This represents the electromagnetic torque estimate obtained in the next sampling period under the applied candidate voltage vector. These represent the stator flux linkage estimates obtained in the next sampling period under the applied candidate voltage vector; , , , These represent the weighting coefficients for torque, flux linkage, switching losses, and current constraints, respectively. This indicates that the inverter has switched from its current switching state to the candidate voltage vector. V ( i The number of switching actions required for the corresponding switching state is used to measure switching losses. Indicates the current constraint penalty term; The index number represents the candidate voltage vector, used to distinguish different candidate voltage vectors. Indicates the first One candidate voltage vector; Step 4-4: Establish the optimal control vector selection formula, and select the candidate voltage vector with the smallest comprehensive performance index as the optimal control vector. The optimal control vector selection formula is as follows: 。 6. The robust prediction method for a deadbeat-free model of a permanent magnet synchronous motor according to claim 5, characterized in that: The inverter is a three-phase two-level inverter, and the inverter has six effective switching states.

7. The robust prediction method for a deadbeat-free model of a permanent magnet synchronous motor according to claim 5, characterized in that: After step 4-2, according to and + The parity correction sign, if and + If the parity of the candidate voltage vectors is the same, then a negative sign is taken when calculating the vertex coordinates corresponding to the candidate voltage vector. and + If the parity of the candidate voltage vector is different, a positive sign is taken when calculating the vertex coordinates corresponding to the candidate voltage vector. The negative sign and the negative sign are used to distinguish the directionality of adjacent regions during the numbering of the standardized triangular region.

Images