Luenberger Position Observer Based on Deadbeat-Current Predictive Control for Sensorless PMSM

Two problems can cause control performance degradation on permanent magnet synchronous motor (PMSM) systems, namely, fluctuation of PMSM parameters and the time delay between current sampling and command value update. In order to reduce the influence of these problems, a new current-predictive control strategy is proposed in this article for medium- and high-speed PMSM. This strategy is based on the discrete mathematical model of PMSM. This new control strategy consists of two main steps: First, an integrator is applied to calculate current compensation value; second, the predictive current value is obtained through deadbeat-current predictive method. The stability of predictive control system is also proved in the article. With this deadbeat-current predictive control scheme, the real current can reach the desired value within one control-step. Based on this new current control method, Luenberger observer and phase-locked loop position tracker is applied in this article. Experimental results for 0.4 kW surface-mounted PMSM confirm the validity and excellent performance for parameters fluctuation of new current predictive control.


Introduction
Nowadays, permanent magnet synchronous motors (PMSMs) have been widely used in industrial applications due to their high power and torque density, small size and wide-speed-range. In PMSM control system, the digital control is common applied for its high precision and excellent control performance. In particular, sensorless control and model predictive control are drawing increasing attention due to the concept of optimum control [1][2][3][4][5].
In sensorless control, position observer based on the voltage equations has been developed [6]. However, reliable operation of PMSM at high speeds still remains many challenges. For example, in medium-and high-speed, fundamental excitation model-based methods are applied to estimate position information [7][8][9][10][11][12][13][14][15][16][17]. Sliding model observer (SMO) is proposed in Reference [7] to estimate position, but signal chattering still exists in dynamic-state operation. A general phased-locked loop (PLL) position tracker is proposed to obtain position and speed information [8], but the gains in PLL are not easy to determine for its complex structure. A relationship between extended back electromotive force (extended back-EMF) and position error is built in [9], a model adaptive-reference system (MARS) is used in [10]. All of these methods are dependent on PMSM parameters excessively.
A nonlinear observer [11,12] is presented to adapt nonlinear characteristic of PMSM, but the relationship between system stability and feedback gains is not clear and the robustness of observer cannot reach to desired performance. The method in [13] is based on cubature Kalman filter (CKF), which needs a large amount of calculation within one PWM control step, hence the method is difficult to realize in engineering. A Luenberger observer is built to observe position information in static-coordinate in References [14][15][16][17], which causes high phase-lag on estimated signal.
In current controller, model predictive control has been widely applied to achieve various functions, such as speed regulation [18], flux control [19], torque control [20] and current control [21][22][23][24][25][26][27]. In general, current control is divided into three parts: hysteresis control [21], proportional PI controller [22] and current predictive control [23][24][25][26][27]. In recent years, an increasing amount of attention has been drawn to current predictive control for its concept of optimum control. The calculation for command voltage value in controller needs a large amount of time, which causes time delay on command value updating, hence the aim of deadbeat predictive control is to predict the command value. Finally, the command value can follow desired current value in PMSM. Compared with conventional current predictive, deadbeat current control can catch command value within one control step under prediction principle [28]. However, the deadbeat-current predictive control is dependent on exact PMSM model, which means the mismatch of system parameters will cause the hysteresis between actual current and desired value [29]. In order to solve this problem, several methods are proposed [30][31][32]. In Reference [30], stator current error is reduced by using a current-regulated delta modulator with current error-correction. A robust current predictive control method [31] is proposed to reduce the current error, caused by inductance fluctuation, between measured error and desired command value. A method based on generalized nonlinear model [32] is presented to predict current, which shows excellent performance on reduce error in dynamic-state operation.
In this paper, a sensorless Luenberger position observer based on the proposed deadbeat-current predictive control is applied. Compared with aforementioned sensorless control methods, the method in this paper can decrease signal phase-lag and few fluctuations on model parameters will not produce the large position error, because a current compensation method is proposed to reduce current ripple. Moreover, the current value can catch current desired value within one PWM control period. The contribution of this article can be concluded as follows: • In order to solve the phase-lag of estimated extended back-EMF, a Luenberger observer and two-order PLL position tracker are applied to observe the rotor position and speed. • Key PWM time-sequence logic relationship between current sampling and voltage command update is established. Based on this relationship, a deadbeat-current predictive method is proposed to calculate real-time voltage command value. • Furthermore, due to the unavoidable parameter fluctuation of PMSM during operation, a current compensation method is proposed to reduce current ripple.
The rest of this paper is organized as follows: Section 2 establish the discrete mathematical model of PMSM and applies Luenberger observer to estimate rotor position; Section 3 presents current prediction model by analyzing PWM time-sequence logic and proposes a current compensation method; Section 4 conducts comparative experiments, which verify the validity of the whole control system; Finally, the conclusions of this article are listed in Section 5.

PMSM Mathematical Model
Due to the error between real rotor position and estimated rotor position, two reference rotational frames are considered. One is the real rotational coordinate system, which includes real rotor position θ e and real rotor speed ω e , the other is the estimated rotational coordinate system, which includes estimated rotor positionθ e and estimated rotor speedω e . Figure 1 shows two different rotating frames and one static frame. Then position error and speed error are defined as ∆θ e and ∆ω e , respectively. Noticing that the real position and speed are unknown in operating control period, hence the PMSM model equations should be established in estimated rotor rotational frame. The voltage equation in d-q real rotational frame is established as follows: where u d and u q are stator voltage components in actual rotational frame, i d and i q are stator current components in actual rotational frame, R s is the stator resistance, L d and L q are the dand q-axis inductances, respectively.  Considering the position error, the sampling current can be transferred into estimated rotational frame by using estimated position vector operator (Equation (3)) for the actual position is unknown.
Submitting the vector operator Equation (3) (4) In this article, the aim of Luenberger observer is to obtain extended back-EMF. Since the sampling current in rotational frame satisfies the Equation (4), hence the Equation (4) Based on reference model Equation (5), the current observer equation is designed as: For the convenience of analysis, using complex vector to transfer Equation (1) into: where u dq s = u d + ju q and i dq s = i d + ji q are the stator voltage and current vectors. Considering the position error, the sampling current can be transferred into estimated rotational frame by using estimated position vector operator (Equation (3)) for the actual position is unknown.
Submitting the vector operator Equation (3) into Equation (2), the PMSM model in the (γ − δ) estimated rotor reference frame is obtained, where the stator voltage equation is shown as follows: In this article, the aim of Luenberger observer is to obtain extended back-EMF. Since the sampling current in rotational frame satisfies the Equation (4), hence the Equation (4) is set as reference model of current observer. Ignoring the parameters fluctuation in PMSM, using command voltage value u γδ * s to replace actual voltage value u γδ s in estimated frame. Considering the speed error keep in steady-state operation is zero (d∆θ/dt = ∆ω = 0), the reference PMSM model Equation (4) can be rewritten as: Based on reference model Equation (5), the current observer equation is designed as: s . Since the current components in rotational frame is direct current (DC), the result is deduced as: Apparently, estimated extended back-EMF is related to current error, hence the amplitude of estimated extended back-EMF can be adjusted by current error, the adjustment strategy is given by: In order to analysis the frequency characteristic of observer system, submitting Equation (8) into Equation (6), and using Laplace transform to figure out the relationship between E γδ andÊ γδ , the relationship is described as:Ê Noting that the observer is a second-order system, which produces hysteresis on high-frequency signal above cut-off frequency. Nevertheless, the estimated variate in rotational frame is a direct signalÊ γδ , hence the hysteresis of the second-order system makes no difference on E γδ . In conclusion, the structure of extended back-EMF Luenberger observer is shown in Figure 2.
Electronics 2020, 9, Apparently, estimated extended back-EMF is related to current error, hence the amplitude of estimated extended back-EMF can be adjusted by current error, the adjustment strategy is given by: In order to analysis the frequency characteristic of observer system, submitting Equation (8) into Equation (6), and using Laplace transform to figure out the relationship between γδ E and ˆγ δ E , the relationship is described as: Noting that the observer is a second-order system, which produces hysteresis on high-frequency signal above cut-off frequency. Nevertheless, the estimated variate in rotational frame is a direct signal ˆγ δ E , hence the hysteresis of the second-order system makes no difference on γδ E . In conclusion, the structure of extended back-EMF Luenberger observer is shown in Figure 2.

Second-Order PLL Rotor Position Tracker
In this article, PLL is applied to obtain the rotor position information from estimated extended back-EMF. Writing the vector ˆγ

Second-Order PLL Rotor Position Tracker
In this article, PLL is applied to obtain the rotor position information from estimated extended back-EMF. Writing the vectorÊ γδ into matrix form as follows: In proposed method, the input of PLL is position error ∆θ, which can be calculated by following equation under the assumption of sin ∆θ ≈ ∆θ: After determining system input, a position tracker structure is designed in Figure 3.

of 17
Electronics 2020, 9, x FOR PEER REVIEW 5 of 18 After determining system input, a position tracker structure is designed in Figure 3. In order to determine the gains 1 K and K 2 , the system is equivalent simplified to a second-order system, which is shown in Figure 4. Hence the transfer function can be described by: e e θ s K s K G s θ s s K s K The parameters 1 K and K 2 of PI (Proportional-Integral) controller are determined by choosing the damping ratio . = ς 0 707 . In the meanwhile, the estimated error of the system is calculated by Equation (13) The rotor speed in a steady-states is a constant, the rotor position is regarded as a ramp signal, whose Laplace form is 2 ( ) / e e θ s ω s = . By using final value theorem, the accuracy of the PLL position tracker can be judged by: Obviously, the result shows high reliability of PLL position tracker. In actual control process, the sampling current is discrete signal. An essential issue should be considered in discrete model, which is that the estimated speed information ˆe ω and operator θˆe j e are unknown during operation. In order to obtain ˆe ω and θˆe j e in operation period, the ˆe ω and ˆe θ in operation PWM period should be calculated by last period. The compensation strategy is designed as: In order to determine the gains K 1 and K 2 , the system is equivalent simplified to a second-order system, which is shown in Figure 4.
Electronics 2020, 9, x FOR PEER REVIEW 5 of 18 After determining system input, a position tracker structure is designed in Figure 3. In order to determine the gains 1 K and K 2 , the system is equivalent simplified to a second-order system, which is shown in Figure 4. Hence the transfer function can be described by: The parameters 1 K and K 2 of PI (Proportional-Integral) controller are determined by choosing the damping ratio . = ς 0 707 . In the meanwhile, the estimated error of the system is calculated by Equation (13) The rotor speed in a steady-states is a constant, the rotor position is regarded as a ramp signal, whose Laplace form is 2 ( ) / e e θ s ω s = . By using final value theorem, the accuracy of the PLL position tracker can be judged by: Obviously, the result shows high reliability of PLL position tracker. In actual control process, the sampling current is discrete signal. An essential issue should be considered in discrete model, which is that the estimated speed information ˆe ω and operator θˆe j e are unknown during operation. In order to obtain ˆe ω and θˆe j e in operation period, the ˆe ω and ˆe θ in operation PWM period should be calculated by last period. The compensation strategy is designed as: Hence the transfer function can be described by: The parameters K 1 and K 2 of PI (Proportional-Integral) controller are determined by choosing the damping rati ς = 0.707. In the meanwhile, the estimated error of the system is calculated by The rotor speed in a steady-states is a constant, the rotor position is regarded as a ramp signal, whose Laplace form is θ e (s) = ω e /s 2 . By using final value theorem, the accuracy of the PLL position tracker can be judged by: Obviously, the result shows high reliability of PLL position tracker.
In actual control process, the sampling current is discrete signal. An essential issue should be considered in discrete model, which is that the estimated speed informationω e and operator e jθ e are unknown during operation. In order to obtainω e and e jθ e in operation period, theω e andθ e in operation PWM period should be calculated by last period. The compensation strategy is designed as: where k is the number of control period and T s is the time of one control period.

Proposed and Conventional Current Prediction Model
In order to realize optimum control on PMSM, the concept of current prediction was referred to track real-time command current value for several years. On the other hand, the input of observer Electronics 2020, 9, 1325 6 of 17 is command voltage value and current value. In order to guarantee the accuracy of observer input, the deadbeat-current predictive control is used. The fundament of current predictive control is discrete model of PMSM. In this article, Forward Euler method is applied to discretize model Equation (1), the discrete mathematical model of PMSM can be deduced as: where, i(k + 1) and i(k) are discrete current vectors, i(k

Conventional Current Predictive Model
In actual computer control system, the command voltage in current time cannot be applied in operating control period. Therefore, the command value in next control period should be predicted in running control period. Equation (17) can be extended to next control period: An assumption is set to be a precondition, which means that the command current value in operating control period can track the command current in control period after next with deadbeat. The relationship can be written as: Submitting Equation (19) into Equation (17) and Equation (18), the conventional current predictive model is obtained as: The conventional current predictive method neglects the fluctuation of current value within predictive period, which influences the predictive control performance. Hence, a new predictive control method is proposed in this article.

Deadbeat Linearized Current Predictive Model
Considering the influence of PWM time-sequence on predictive period, the logic of PWM sampling time point and value updating time point are necessary to be analyzed. The relationship between sampling time point and command updating time point is shown in Figure 5.
Noticing that the sampling time point is the midpoint of per PWM control period, the voltage updating time point is the beginning of every period. According to the closed-loop control strategy of PMSM, the voltage command value is calculated by a controller whose input is current command value. The sampling current i(k − 1) and are used to predict voltage command value u * (k + 1), u * (k + 2), respectively, hence the predictive period T = 1.5T s The predictive period is longer than the control period, which means that the actual current cannot reach to desired command value within one PWM control period based on conventional predictive method. Noticing that the sampling time point is the midpoint of per PWM control period, the voltage updating time point is the beginning of every period. According to the closed-loop control strategy of PMSM, the voltage command value is calculated by a controller whose input is current command value. The sampling current ( ) k− i 1 and are used to predict voltage command value , respectively, hence the predictive period T The predictive period is longer than the control period, which means that the actual current cannot reach to desired command value within one PWM control period based on conventional predictive method. Since the sampling time is usually short, the fluctuation of rotor speed tends to be zero and the change of current within one PWM period is approximated to linear. Moreover, the average voltage value within half of a PWM control period is voltage value in running PWM control period. For instance, the average voltage value in half of period [ , ] Figure 5.
Considering this condition, the voltage can be integrated within one predictive period as following equations: The current command value in next period cannot be obtained in running PWM control period, therefore the command value in running PWM control period is used to replace the unknown value. Submitting Equation (21) into Equation (1), the current predictive model is obtained as: Compared with conventional predictive control method, the method presented in this paper considers the PWM time-sequence with linear current prediction, which improves the performance on tracking current and realizes deadbeat control. Since the sampling time is usually short, the fluctuation of rotor speed tends to be zero and the change of current within one PWM period is approximated to linear. Moreover, the average voltage value within half of a PWM control period is voltage value in running PWM control period. For instance, the average voltage value in half of period [t 1 , Figure 5. Considering this condition, the voltage can be integrated within one predictive period as following equations: The current command value in next period cannot be obtained in running PWM control period, therefore the command value in running PWM control period is used to replace the unknown value. Submitting Equation (21) into Equation (1), the current predictive model is obtained as: Compared with conventional predictive control method, the method presented in this paper considers the PWM time-sequence with linear current prediction, which improves the performance on tracking current and realizes deadbeat control.

Current Compensation
The aim of current compensation is to improve the disturbance rejection behavior. In this article, the vector control strategy is based on i * d = 0, therefore only the i q current predictive model in Equation (22) should be considered. Equation (22) can be rewritten as: Equations (22) and (23) are the most important equations in this article. Ideally, the current in next PWM control period can catch the current command value i * q . However, in actual case the parameters fluctuation of PMSM are not sufficiently small to be ignored, hence the disturbance factors should be analyzed. Equation (23) is established as the ideal model of voltage control, the actual voltage control model is established as Equation (24), which includes actual parameters of PMSM in running time: where R s , L q , λ f are actual stator resistance, inductance in q-axis and PM flux, respectively. In order to explore the relationship between parameters fluctuation and current error, subtract Equation (23) from Equation (24) and the result is: includes the disturbance produced by ∆R s and ∆L q , i q2 (k + 1) includes the disturbance produced by ∆λ f . Noticing that T s is sufficient small, hence one inequation relationship can be established: 3R s T s 4L q . Based on this relationship, Equation (25) is simplified as: The transfer function in discrete form of Equation (26) is described by: Equation (27) indicates that the current can follow the command value within one PWM control period when L q /L q = 1.
Considering more interference factors, the sakes for current error are not limited by parameters fluctuation. Low-order harmonics produced by SVPWM, dead-time and voltage drop of inverter, bus voltage ripple and measurement noise of system can also cause current error. In order to reduce current error, all of the above interference factors in operation control period are regarded as one disturbing term. Considering f q (k), the predictive current in q-axis consists of 3 parts, which is described as: According to Equation (28), this article presents a compensation method by using integrator to obtain the compensation value i * * q (k). Adding the compensation value to current command value, the compensated value can be applied to predict voltage command.
Using new current command value i * * q (k) to replace i * q (k) in Equation (22), the predicted voltage command value is deduced as: Hence, the complete current predictive method is shown in Figure 6. The real current in PMSM will increase with the compensated current command value increases until the real current catch the command value.
Hence, the complete current predictive method is shown in Figure 6. The real current in PMSM will increase with the compensated current command value increases until the real current catch the command value.

System Stability Analysis
In order to determine the gain in integrator in Figure 6, the prediction and compensation system stability should be analyzed. Taking q i component for example and defining that 3 K is the gain in integrator, the discrete mathematical model of current prediction block in Figure 6 can be described by: Using ** ( ) q i k to replace * ( ) q i k in Equation (26) and submitting Equation (30) into Equation (26), the predictive current in next 2 periods control period can be obtained as:

System Stability Analysis
In order to determine the gain in integrator in Figure 6, the prediction and compensation system stability should be analyzed. Taking i q component for example and defining that K 3 is the gain in integrator, the discrete mathematical model of current prediction block in Figure 6 can be described by: Using i * * q (k) to replace i * q (k) in Equation (26) and submitting Equation (30) into Equation (26), the predictive current in next 2 periods control period can be obtained as: Hence the transfer function in discrete form between i * q (z) and i q (z) can be solved by simultaneous Equations (29) and (30), the result is: The system stability is determined by system characteristic Equation (34), which is based on transfer function Equation (33).
According to Juli criterion, the relationship between each coefficient of Equation (34) can be obtained as follows: Electronics 2020, 9,1325 10 of 17 Hence the condition of system stability can be proved from Equations (35) and (34), which is deduced as: The Equation (36) can be regarded as a condition to limit the gain in integrator. Based on all control methods proposed in this article, the whole control strategy can be described by Figure 7.
According to Juli criterion, the relationship between each coefficient of Equation (34) can be obtained as follows: Hence the condition of system stability can be proved from Equations (35) and (34), which is deduced as: The Equation (36) can be regarded as a condition to limit the gain in integrator. Based on all control methods proposed in this article, the whole control strategy can be described by Figure 7.

Experimental Results
The experimental drive setup consists of a 0.4 kW PMSM connected to a dynamometer machine and a load motor used for loading. The PMSM has a rotor with surface-mounted magnets, the parameters of PMSM are shown in Table 1. Both of load motor and tested motor are supplied by a inverter which is connected with a common DC link and controlled by Infineon TC1782, the sampling-time and PWM frequency is f s = 10 kHz and the code is written in C Language. Figure 8 shows the photograph of PMSM experimental platform.
This drive uses the proposed current prediction and compensation controllers in a synchronous frame, with constant d-current reference i * d = 0. Both of load motor and tested motor are supplied by a inverter which is connected with a common DC link and controlled by Infineon TC1782, the sampling-time and PWM frequency is fs = 10 kHz and the code is written in C Language. Figure 8 shows the photograph of PMSM experimental platform.

Experimental Results of Deadbeat-current predictive Control Method
In this article, two current predictive control methods are referred, one is conventional predictive method, the other is deadbeat-current predictive method, which is the new method designed in this article. In order to show the advantage of deadbeat control, the experiments based on same operating condition are conducted by 2 different predictive methods. Figure 9a shows the speed response of 2 current predictive methods in dynamic-state operation at 1000 r/min, 2000 r/min and 3000 r/min with step signal at 5 s, 30 s and 60 s, respectively and a rated load of 0.2 N-m. Figure 9b

Experimental Results of Deadbeat-Current Predictive Control Method
In this article, two current predictive control methods are referred, one is conventional predictive method, the other is deadbeat-current predictive method, which is the new method designed in this article. In order to show the advantage of deadbeat control, the experiments based on same operating condition are conducted by 2 different predictive methods. Figure 9a shows the speed response of 2 current predictive methods in dynamic-state operation at 1000 r/min, 2000 r/min and 3000 r/min with step signal at 5 s, 30 s and 60 s, respectively and a rated load of 0.2 N·m. Figure 9b,c show speed error of 2 current predictive methods in steady-state operation at 3000 r/min with a rated load of 0.2 N·m. Although the rotor speed response of PMSM in Figure 9a based on both 2 methods can follow the speed command value, Figure 9b,c show that the speed error by deadbeat predictive control is smaller than conventional methods and can be controlled within 20 rpm and the error fluctuation is also smaller than conventional methods. As depicted in Figure 9, the fluctuation of current error based on deadbeat-current predictive method is smaller, especially in q-axis and the current error based on new control strategy is within 0.4 A, which can reduce motor system vibration effectively. Figure 10 shows the current error based on two different current predictive control methods with 0.2 N-m load in 3000 rpm steady-state. Figure 10a,b show the current error in d-axis based on two different current predictive control methods and Figure 10c,d show the current error in q-axis based on two different current predictive control methods. As depicted in Figure 9, the fluctuation of current error based on deadbeat-current predictive method is smaller, especially in q-axis and the current error based on new control strategy is within 0.4 A, which can reduce motor system vibration effectively. Figure 10 shows the current error based on two different current predictive control methods with 0.2 N·m load in 3000 rpm steady-state. Figure 10a,b show the current error in d-axis based on two different current predictive control methods and Figure 10c,d show the current error in q-axis based on two different current predictive control methods. As depicted in Figure 9, the fluctuation of current error based on deadbeat-current predictive method is smaller, especially in q-axis and the current error based on new control strategy is within 0.4 A, which can reduce motor system vibration effectively. Figure 10 shows the current error based on two different current predictive control methods with 0.2 N-m load in 3000 rpm steady-state. Figure 10a,b show the current error in d-axis based on two different current predictive control methods and Figure 10c,d show the current error in q-axis based on two different current predictive control methods.  Figure 11 shows the current response in q-axis based on 2 different methods when the speed command value is changed from 2000 rpm to 3000 rpm suddenly. Obviously, the current response based on new methods shows better performance for its smaller fluctuation.

Experimental Results of the Proposed Current Compensation with the System Parameter Fluctuation
In order to prove the current compensation effect to reduce the current error caused by parameter fluctuation, the experiment in this part is conducted by setting the motor parameter in controller as 1.5 times higher than actual PMSM: ' 1.5 Figure 12a, c show experimental current before adding the compensation value, Figure 12b,d show compensated current. All of these experiments are based on inductance change in controller.  Figure 11 shows the current response in q-axis based on 2 different methods when the speed command value is changed from 2000 rpm to 3000 rpm suddenly. Obviously, the current response based on new methods shows better performance for its smaller fluctuation.  Figure 11 shows the current response in q-axis based on 2 different methods when the speed command value is changed from 2000 rpm to 3000 rpm suddenly. Obviously, the current response based on new methods shows better performance for its smaller fluctuation.

Experimental Results of the Proposed Current Compensation with the System Parameter Fluctuation
In order to prove the current compensation effect to reduce the current error caused by parameter fluctuation, the experiment in this part is conducted by setting the motor parameter in controller as 1.5 times higher than actual PMSM: ' 1.5 Figure 12a, c show experimental current before adding the compensation value, Figure 12b,d

Experimental Results of the Proposed Current Compensation with the System Parameter Fluctuation
In order to prove the current compensation effect to reduce the current error caused by parameter fluctuation, the experiment in this part is conducted by setting the motor parameter in controller as 1.5 times higher than actual PMSM: Figure 12a,c show experimental current before adding the compensation value, Figure 12b,d show compensated current. All of these experiments are based on inductance change in controller. Apparently, compared with no compensated current, compensated current is closer to current command value and the steady current error is smaller than no compensated value. Hence, the compensation algorithm proposed in this article shows satisfied performance on suppression of inductances fluctuation. . Figure 13 shows the current response in d-and q-axis, respectively.
Apparently, compared with no compensated current, compensated current is closer to current command value, the steady current error is smaller than no compensated value. Hence, the compensation algorithm proposed in this article shows satisfied performance on suppression of PM flux fluctuation. The experiment for PM flux change of PMSM is also conducted by changing the motor PM flux value in controller as λ f = 1.5λ f . Figure 13 shows the current response in dand q-axis, respectively. Apparently, compared with no compensated current, compensated current is closer to current command value, the steady current error is smaller than no compensated value. Hence, the compensation algorithm proposed in this article shows satisfied performance on suppression of PM flux fluctuation.
Noticing that the current response of the deadbeat-current predictive control in all figures shows a quick overshoot, the overshoot can be decreased by choosing appropriate gain in current compensation integrator, but the system settling time will be extended.
Based on all experimental result, the proposed current compensation method in this article indicates the ability to reduce the current error caused by parameter fluctuation.
The experiment for PM flux change of PMSM is also conducted by changing the motor PM flux value in controller as ' 1.5 . Figure 13 shows the current response in d-and q-axis, respectively.
Apparently, compared with no compensated current, compensated current is closer to current command value, the steady current error is smaller than no compensated value. Hence, the compensation algorithm proposed in this article shows satisfied performance on suppression of PM flux fluctuation. Noticing that the current response of the deadbeat-current predictive control in all figures shows a quick overshoot, the overshoot can be decreased by choosing appropriate gain in current compensation integrator, but the system settling time will be extended.
Based on all experimental result, the proposed current compensation method in this article indicates the ability to reduce the current error caused by parameter fluctuation.

Experimental Results of Rotor Position and Speed Sensorless Observer.
In this article, Luenberger observer and PLL position tracker is applied to estimate rotor position. In order to prove the accuracy of estimated position information, the position error is set as an important index to evaluate the accuracy. The experiment is conducted in 3 different steady-state operations with no load and 0.4 N ⋅ m load, respectively.

Experimental Results of Rotor Position and Speed Sensorless Observer
In this article, Luenberger observer and PLL position tracker is applied to estimate rotor position. In order to prove the accuracy of estimated position information, the position error is set as an important index to evaluate the accuracy. The experiment is conducted in 3 different steady-state operations with no load and 0.4 N·m load, respectively. Figure 14a,c,e shows the position error in 1000 rpm, 2000 rpm and 3000 rpm without load, respectively, and Figure 14b,d,f shows the position error in 1000 rpm, 2000 rpm and 3000 rpm with 0.4 N·m load, respectively. The result indicates that position error can be controlled within ± 4 • in different steady-states, no matter there is load or not. Noticing that the current response of the deadbeat-current predictive control in all figures shows a quick overshoot, the overshoot can be decreased by choosing appropriate gain in current compensation integrator, but the system settling time will be extended.
Based on all experimental result, the proposed current compensation method in this article indicates the ability to reduce the current error caused by parameter fluctuation.

Experimental Results of Rotor Position and Speed Sensorless Observer.
In this article, Luenberger observer and PLL position tracker is applied to estimate rotor position. In order to prove the accuracy of estimated position information, the position error is set as an important index to evaluate the accuracy. The experiment is conducted in 3 different steady-state operations with no load and 0.4 N ⋅ m load, respectively.

Conclusions
In this article, a Luenberger rotor sensorless observer based on a proposed current control method is built for the wide-speed-range operation. The conclusions of this paper are listed as follows: • To solve the phase-lag of estimated extended back-EMF, Luenberger observer and second-order PLL position tracker, in rotational frame, are applied to observe the rotor position and speed, because extended back-EMF signal in rotational frame is a DC signal. • A deadbeat-current predictive method, by predicting current change is linear within one predictive control period, is proposed to obtain real-time voltage command value.

•
Due to the unavoidable parameter fluctuation of PMSM during operation, a current compensation method is proposed, which can reduce current error caused by parameter fluctuation. The method of determining the gain in integrator is designed by system stability.

•
Finally, the effectiveness and the robustness of the proposed control system are verified. The experimental results show the static and dynamic performance of the new current controller. In the meanwhile, the accuracy of estimated position is proved as well, and the minimum position error can be controlled within ±4° whether the load is given or not.

Conclusions
In this article, a Luenberger rotor sensorless observer based on a proposed current control method is built for the wide-speed-range operation. The conclusions of this paper are listed as follows: • To solve the phase-lag of estimated extended back-EMF, Luenberger observer and second-order PLL position tracker, in rotational frame, are applied to observe the rotor position and speed, because extended back-EMF signal in rotational frame is a DC signal. • A deadbeat-current predictive method, by predicting current change is linear within one predictive control period, is proposed to obtain real-time voltage command value.

•
Due to the unavoidable parameter fluctuation of PMSM during operation, a current compensation method is proposed, which can reduce current error caused by parameter fluctuation. The method of determining the gain in integrator is designed by system stability.

•
Finally, the effectiveness and the robustness of the proposed control system are verified. The experimental results show the static and dynamic performance of the new current controller. In the meanwhile, the accuracy of estimated position is proved as well, and the minimum position error can be controlled within ±4 • whether the load is given or not.