# Sensorless SPMSM Position Estimation Using Position Estimation Error Suppression Control and EKF in Wide Speed Range

Authors: Zhanshan Wang; Longhu Quan; Xiuchong Liu
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

The control of a high performance alternative current (AC) motor drive under sensorless operation needs the accurate estimation of rotor position. In this paper, one method of accurately estimating rotor position by using both motor complex number model based position estimation and position estimation error suppression proportion integral (PI) controller is proposed for the sensorless control of the surface permanent magnet synchronous motor (SPMSM). In order to guarantee the accuracy of rotor position estimation in the flux-weakening region, one scheme of identifying the permanent magnet flux of SPMSM by extended Kalman filter (EKF) is also proposed, which formed the effective combination method to realize the sensorless control of SPMSM with high accuracy. The simulation results demonstrated the validity and feasibility of the proposed position/speed estimation system.


---

## Body

## 1. Introduction

Recently, permanent magnet synchronous motor (PMSM) has been widely used in many automation control fields due to its advantages of superior power density, maintenance-free operation, and high controllability [1–3]. The conventional motor servo control needs a position sensor such as incremental encoder to measure the rotor position and feed it back to the controller for ensuring the accurate control. This sensor increases the cost and size while decreasing the reliability of the drive. Therefore, a position sensorless control for PMSM drive has become popular research direction [4–10].The primary methods for position estimation under the sensorless condition can be divided into several main categories: adaptive approaches, back-electromagnetic force (back-EMF) based methods, reduced order observer methods, EKF based methods, and signal injection methods for low speed range [3–10]. In these studies, the estimation result is very dependent on the accuracy of system parameters that include the motor parameters like inductance, resistance, and magnetic flux and the unmodeled parameters like inverter dead-time effect and inductance cross-saturation. However, the system parameters vary continuously during the motor dynamic, so the estimation result may be biased or cause the vibration.In order to realize the accurate position estimation, estimation of rotor position and identifying of system parameters must be employed simultaneously. Reference [11] presents a sensorless position estimation method using the combination of identifying the stator resistance by EKF and estimating the rotor position by back-EMF. References [12–15] present a sensorless position estimation method using the combination of identifying the motor parameters by EKF and estimating the rotor angular speed by model reference adaptive system (MRAS), in which EKF has become one of the most effective ways to identify the parameters of nonlinear system like AC motor [16, 17]. These studies are concentrated to improve the estimation accuracy of rotor position by identifying one or two of the key parameters like inductance, resistance, and flux. However, because of unmodeled nonlinearities, for example, inverter dead-time effect, inductance cross-saturation, and other motor parameters drift and variation, there will be the uncertain phenomena in a close-loop sensorless control with estimated position and speed.Now, in order to eliminate the influence of system uncertainty that affects the control result, the artificial intelligent techniques such as ANN were applied to estimate the position and speed of sensorless AC motor [18, 19]. These methods can increase the robustness of the position and speed estimation through making a nonlinear function between the input and output of the control system. But they require the relative long operation time for the real-time process in an actual control system. Many researchers tried to reduce the control error that resulted from the system uncertainty through using the complex number analysis of a control system [20–22]. By the use of the instantaneous power that separated into the real part and imaginary part by the complex number analysis, a rotor position estimation method unrelated to the actual speed and flux linkage in a rated speed area was studied in [22]. In [22], the system parameters used in the position estimation were reduced to the phase inductance and resistance in comparison to the conventional methods. This approach not only is simple in the system analysis, but also can easily separate the system variables into the real part and the error (imaginary) part caused by the system uncertainty. In this paper, a novel method that eliminates the influence of entire system uncertainty in the position estimation of sensorless SPMSM is realized by the motor complex number model. However, severe demagnetization in the flux-weakening region can break estimating the position accurately and can also cause the instability problem in the whole control system [10]. The flux of motor is following a Gaussian distribution in the flux-weakening region [1]. References [5, 23] revealed that EKF is useful in identifying the accidental variables that are following the Gaussian distribution. In the previous studies, identifying the permanent magnet flux by EKF has focused on the interior permanent magnet synchronous motor (IPMSM) and asynchronous motor, but SPMSM has not been studied.This paper presents a novel robust position estimation method in which its dynamic position estimation performance is not affected by system parameter variations. The rotor position angle is obtained by a novel position estimator over the wide speed range including low speed region and flux-weakening region. First, the mathematical model of position estimator composed of the real value (real part) and error value (imaginary part) of the estimated rotor position angle is presented. Second, in order to suppress the position estimation error to zero, a position type PI controller is applied to imaginary part of the proposed position estimator. The proposed method can enhance the accuracy and robustness of position estimation by suppressing to zero the estimation error that may be caused by the uncertainty of system parameters in a rated speed range including the low speed region. Third, in order to expand the robustness to the flux-weakening region, the permanent magnet flux identification model of SPMSM by EKF is proposed; the affection that is caused by uncertainty of permanent magnet flux during the motor at the flux-weakening region is eliminated.The paper is organized as follows. The mathematical model of SPMSM and the principle of the proposed position estimation scheme are shown in Section2. The combination of a complex number model based position estimator and an error suppression PI control is given in Section 3. An online identification scheme of permanent magnet flux based on EKF is presented in Section 4. The SPMSM vector control system and detailed control scheme are given in Section 5. The simulation block and simulation results are presented in Section 6, which verify the effectiveness of the proposed control method.


## 2. Basic Control Equations of SPMSM and Novel Sensorless Control Scheme

### 2.1. Basic Control Equations of SPMSM

In the paper, the vector control of SPMSM is considered in a rotor-fixed frame (d-q frame) [8, 10] as follows:
(1)ud=RSid+LSddtid-LSiqωk,uq=RSiq+LSddtiq+LSidωk+ψmωk,(2)ψd=LSid+ψm,ψq=LSiq,
where ud, uq, id, and iq are d-q components of the voltage and current in stator and ωk is an electrical angle speed of rotor in d-q frame. RS and LS are the resistance and inductance of stator, respectively. ψm is a rotor permanent magnet flux; ψd and ψq are the d-q components of the actual air gap flux of motor.The electrical angle speed of rotor is(3)ωk=dθkdt,
where θk is the field angle between the d axis of rotor-fixed frame and the α axis of stator-fixed frame (α-β frame) [8, 10].Considering the mechanical load, the dynamic equation of SPMSM drive system is obtained by(4)Jddtωk=Te-TL-Bωk,
where Te is an electromagnetic torque, TL is a load torque, J is an inertia moment, and B is a viscous damping coefficient.The electromagnetic torque of motor is given as follows:(5)Te=pψmiq,
where p is the number of pole pairs.In the SPMSM vector control,dTL/dωk>0 is supposed for all practical speeds [8]. The mentioned system is usually used in simulations as well as in design of the current controller and speed controller. In SPMSM vector control of this paper, the d axis current reference is controlled to be zero and the motor torque is controlled by the q axis current.


### 2.2. Novel Sensorless Estimation Scheme

The proposed novel position estimation scheme is shown in Figure1. Uabc, i^, and Δi^ refer to a three-phase stator voltage of SPMSM, an estimated current, and an error between command current and estimated current; U^ and ψ^ refer to an output of current controller and an identified flux obtained from EKF; θ^ and ω^ refer to an estimated rotor position angle and an estimated rotor angular speed.Figure 1
Block diagram of proposed estimation scheme.In the previous studies [24–28], the “rotor position estimation block” calculated the position estimation value directly using i^,U^, and motor parameters without considering the estimation error. The estimation result was too dependent on the accuracy of system parameters. The system parameters contain a certain degree of uncertainty, so the estimation result could not be well. In Figure 1, the “rotor position estimation block” is composed of the real position estimation value and the position estimation error that can be caused by the uncertainty of system parameters. In order to realize this separation, we have derived the SPMSM complex number model based position estimator that can easily separate the real value and error value of position estimation. The estimation result is composed of the sum of real position valve (“real part”) and estimation error value (“imaginary part”). The position estimation error is controlled to be zero by the error suppression PI controller. As a result, the accurate estimation result may be obtained. By identifying the permanent magnet flux by EKF, the affection accepted by the magnetic flux change in the flux-weakening region is eliminated.


## 2.1. Basic Control Equations of SPMSM

In the paper, the vector control of SPMSM is considered in a rotor-fixed frame (d-q frame) [8, 10] as follows:
(1)ud=RSid+LSddtid-LSiqωk,uq=RSiq+LSddtiq+LSidωk+ψmωk,(2)ψd=LSid+ψm,ψq=LSiq,
where ud, uq, id, and iq are d-q components of the voltage and current in stator and ωk is an electrical angle speed of rotor in d-q frame. RS and LS are the resistance and inductance of stator, respectively. ψm is a rotor permanent magnet flux; ψd and ψq are the d-q components of the actual air gap flux of motor.The electrical angle speed of rotor is(3)ωk=dθkdt,
where θk is the field angle between the d axis of rotor-fixed frame and the α axis of stator-fixed frame (α-β frame) [8, 10].Considering the mechanical load, the dynamic equation of SPMSM drive system is obtained by(4)Jddtωk=Te-TL-Bωk,
where Te is an electromagnetic torque, TL is a load torque, J is an inertia moment, and B is a viscous damping coefficient.The electromagnetic torque of motor is given as follows:(5)Te=pψmiq,
where p is the number of pole pairs.In the SPMSM vector control,dTL/dωk>0 is supposed for all practical speeds [8]. The mentioned system is usually used in simulations as well as in design of the current controller and speed controller. In SPMSM vector control of this paper, the d axis current reference is controlled to be zero and the motor torque is controlled by the q axis current.


## 2.2. Novel Sensorless Estimation Scheme

The proposed novel position estimation scheme is shown in Figure1. Uabc, i^, and Δi^ refer to a three-phase stator voltage of SPMSM, an estimated current, and an error between command current and estimated current; U^ and ψ^ refer to an output of current controller and an identified flux obtained from EKF; θ^ and ω^ refer to an estimated rotor position angle and an estimated rotor angular speed.Figure 1
Block diagram of proposed estimation scheme.In the previous studies [24–28], the “rotor position estimation block” calculated the position estimation value directly using i^,U^, and motor parameters without considering the estimation error. The estimation result was too dependent on the accuracy of system parameters. The system parameters contain a certain degree of uncertainty, so the estimation result could not be well. In Figure 1, the “rotor position estimation block” is composed of the real position estimation value and the position estimation error that can be caused by the uncertainty of system parameters. In order to realize this separation, we have derived the SPMSM complex number model based position estimator that can easily separate the real value and error value of position estimation. The estimation result is composed of the sum of real position valve (“real part”) and estimation error value (“imaginary part”). The position estimation error is controlled to be zero by the error suppression PI controller. As a result, the accurate estimation result may be obtained. By identifying the permanent magnet flux by EKF, the affection accepted by the magnetic flux change in the flux-weakening region is eliminated.


## 3. Complex Number Model of SPMSM and Proposed Position Estimation Scheme

Stator voltage in the stator-fixed frame (α-β frame) [1] is defined as
(6)US=RSiS+ψSdt,
where
(7)US=Uα+jUβ,iS=iα+jiβ,ψS=ψα+jψβ,
where RS is a stator resistance. ψS, US, and iS are rotor flux vector, stator voltage vector, and stator current vector in the α-β frame, respectively.Now, theα-β frame rotates with an angular speed ωk, and all quantities can be transformed from the stator-fixed frame to the rotor-fixed frame (d-q frame) [1]; the results are given as follows:
(8)US=UkejθK,iS=ikejθK,ψS=ψkejθK,dψSdt=dψkdtejθK+jωkψkejθK,
where ψk, Uk, and ik are rotor flux vector, stator voltage vector, and stator current vector in the d-q frame, respectively. By substituting (8) into (6), the equation of stator voltage in the d-q frame can be expressed as
(9)Uk=RSik+dψkdt+jωkψk.Equation (9) is usually used for a high performance SPMSM drive. However, this d-q frame cannot be utilized directly in the position sensorless system because the model does not contain any information about the rotor position error. It is necessary to use another reference frame, noted as a δ-γ frame, which lags by position estimation error Δθk from the d-q frame as shown in Figure 2. Δθk is not a constant and varies continuously during the motor dynamic. The SPMSM model in the δ-γ frame can be derived by transforming the voltage equation (9) of d-q frame into the δ-γ frame as follows:
(10)Um=RSim+dψmdt+jωSψm,
where ωS is a rotor electrical angle speed in the δ-γ frame. ωS is a sum of the actual electrical angle speed ωk and speed estimation error Δωk. ψm, Um, and im are rotor flux vector, stator voltage vector, and stator current vector in the δ-γ frame, respectively. If the sensorless control is correct, θS of δ-γ frame is convergent to θk of d-q frame, the position error Δθk is also convergent to zero, and all electromagnetic parameters of δ-γ frame can be processed on the d-q frame. In order to simplify consideration of the position estimation, the position error Δθk is ignored in δ-γ frame and the SPMSM model in d-q frame is employed to estimate the rotor position in position sensorless SPMSM drive system.Figure 2
Reference frame diagram of SPMSM.Equation (10) can be written as
(11)Uk=RSik+LSdikdt+jωkψm(exp⁡j(Δθk)).
Using the Euler formula, we can eliminate j from (11):
(12)Uk=RSik+LSdikdt+dθkdtψm(exp⁡j(Δθk+π2)).We can compute the differential of positionθk from (12):
(13)dθkdt=1ψm·(Uk-RSik-LSdikdt)·exp⁡(-j(Δθk+π2)).The position estimation equation can be obtained from Figure1 and (13) as follows:
(14)dθ^kdt=1ψm·(U^k-RSi^k-LSdi^kdt)·exp⁡(-j(Δθ^k+π2)),
where the U^k,i^k,andΔθ^k of (14) correspond to the U^,i^,andΔθ^ of Figure 1. Index k means that the parameters in estimation equation are values of d-q frame.In order to simplify the consideration, let each term of (14) be as follows:
(15)U^k-RSi^k-LSdi^kdt=v,exp⁡(-j(Δθ^k+π2))=b.Rotor position estimation equation (14) can be written by the modulus (v·b) as follows:
(16)dθ^kdt=1ψm(v·b).For practical implementation, it is necessary to obtain the algorithm in discrete form. Equation (16) can be solved as follows:
(17)θ^k(n)=θ^k(n-1)+Tψm(v(n)·b(n))=θ^k(n-1)+Tψm{[i^k(n)-i^k(n-1)T]U^k(n)-RSi^k(n)-LS[i^k(n)-i^k(n-1)T]}·exp⁡{-j[Δθ^k(n)+π2]}=θ^k(n-1)+Tψm{[jcos⁡(Δθ^k(n))vd(n)+jvq(n)]×[-sin(Δθ^k(n))-jcos⁡(Δθ^k(n))]}=θ^k(n-1)+Tψm{Re^[v(n)b(n)]+Im⁡^[v(n)b(n)]}=θ^k(n-1)+Δθ^k(n).
Now, the estimation θ^k(n) of the actual rotor position θk(n) is determined by adding the estimated angle increment Δθ^k(n) in one sampling period to the previous estimation θ^k(n-1). There are the real part and imaginary part in the angle increment. Real part is the actual increment of estimation, and the imaginary part is the estimation error. If there are no system parameter errors and the rotor position angle is correctly obtained, the imaginary part of (v·b) is zero. Because the system parameters vary continuously during the motor operation, the position estimation error exists in a certain degree. In order to suppress the position estimation error to zero, the proposed estimation method uses a position type PI controller to the imaginary part of position estimation equation.The estimated angle increment is obtained as follows:(18)Δθ^k(n)=Tψm{(1-KiKp)Re[v(n)·b(n)]+Kp[(1-KiKp)Im⁡(v(n)b(n))-Im⁡(v(n-1)b(n-1))(1-KiKp)]+Δθcorr(n-1)(1-KiKp)},
where T, Kp, and Ki are sample time of estimation algorithm, proportional gain of PI controller, and T/Ti, respectively. Ti is an integral time constant. The second term of (18) carries out the function of PI control so that the imaginary part of (17) converges to 0. The position estimation error that may be caused by the uncertainty of system parameters is eliminated through this control.The estimated angular speed incrementΔω^k(n) of the rotation speed is obtained by differentiating estimated position angle increment Δθ^k(n) during the sample time T. The estimated speed ω^k(n) of the rotor is obtained by
(19)ω^k(n)=ω^k(n-1)+Δω^k(n)=ω^k(n-1)+dΔθ^k(n)dt.The imaginary part of (18) must be multiplied by the sign of the estimated speed because the angle between v and b is π when speed is negative direction.In this part, one method of improving the position estimation accuracy and robustness is studied about the sensorless SPMSM vector control. The position estimation error caused by the uncertainty of system parameters could be regulated to zero by this method. The real value and error value of rotor position estimation are correctly separated from the position estimator based on the complex number model of SPMSM. The position estimation error is removed by error suppression PI operation. Using this estimation method, not only is the dynamic estimation of the position and speed fast but also convergence is well, and then the proper estimation accuracy can be obtained.The position estimation method derived in this part is robust on the uncertainty of system parameters in the rated speed range including low speed. However, severe demagnetization in the flux-weakening area can break estimating the position accurately and can cause the instability problem in the whole control system. The permanent magnet fluxψm was processed as a constant in position estimation equation (17). In the next part, we present one method that obtains the identified flux ψ^m of ψm by using EKF, in order to guarantee robustness and estimation accuracy in the flux-weakening area.


## 4. Online Identification of Permanent Magnet Flux Based on EKF

### 4.1. Introduction to EKF

EKF presents the general problem of trying to estimate the statesx∈Rn that is dominated by the nonlinear stochastic difference equation as follows:
(20)xk=f(xk-1,uk,σk-1),
and the measurement output y∈Rn is as follows:
(21)yk=h(xk,ρk),
where uk is the control vector, the random variables σk-1 and ρk represent the process and measurement noise, and the covariance matrices of these noises are defined as
(22)Q=E[σ·σT],R=E[ρ·ρT].The solutions to EKF can be divided into two groups.(1)
Time update equation (predict), which uses the current state and error covariance to obtain the a priori estimations for the next time step.
(i)
Project the state ahead:(23)x^k-=f(x^k-1,uk,0).(ii)
Project the error covariance ahead:(24)Pk-=AkPk-1AkT+WkQk-1WkT.(2)
Measurement update equation (correct), which incorporates new measurements into the a priori estimations to obtain improved a posteriori estimations.
(i)
Compute EKF gain:(25)Kk=Pk-HkT(HkPk-HkT+VkRVkT)-1.(ii)
Update estimation with measurement:(26)x^k=x^k-+Kk(yk-h(x^k-,0)).(iii)
Update the error covariance:(27)Pk=(1-KkHk)Pk-,where x^k- is a priori state estimation at step k given knowledge of the process prior to step k, x^k is a posteriori state estimation at step k given measurement yk, and Pk- and Pk are the covariance of the a priori and a posteriori estimation errors.Ak is a Jacobian matrix of partial derivatives of f with respect to x; Ak[i,j]=(∂f[i]/∂x[j])(x^k-1,uk,0).Wk is a Jacobian matrix of partial derivatives of f with respect to σ; Wk[i,j]=(∂f[i]/∂σ[j])(x^k-1,uk,0).Hk is a Jacobian matrix of partial derivatives of h with respect to x; Hk[i,j]=(∂h[i]/∂x[j])(x^k,0).Vk is a Jacobian matrix of partial derivatives of h with respect to ρ; Vk[i,j]=(∂h[i]/∂ρ[j])(x^k,0).


### 4.2.  Permanent Magnet Flux Identification of SPMSM Based on EKF

Here, we present one method that identifies the magnetic flux change that can be generated during the motor dynamic. In Section3 of this paper, the complex number model of SPMSM is used for estimating the accurate position. In (17), the imaginary part of position estimator caused by the uncertainty of motor parameters is controlled to be zero by the error suppression PI controller. But, the demagnetization amount must not be controlled to be zero in the flux-weakening area. Calculating the demagnetization amount accurately is one of the most important problems to realize the accurate position estimation and vector control in the flux-weakening area. The complex number model composed of real part and imaginary part may cause the complexity in the flux identification. Therefore, the rotor flux is identified through the combination of SPMSM model (1) and EKF in this paper. Equation (1) with considering the magnetic flux change can be rewritten as follows:
(28)ddti^d=-RSi^dLS+i^qω^k+ω^kψ^mqLS+u^dLS,ddti^q=-RSi^qLS-i^dω^k-ω^kψ^mdLS+u^qLS.In order to estimate the magnetic flux, the magnetic flux is chosen as one state variable. Because the magnetic flux linkage cannot change sharply during the motor dynamic [1], its differential is assumed to be zero:
(29)dψ^mddt=0,dψ^mqdt=0.By choosing set of the electric current variables and magnetic flux of theq axis as the state variables, a fourth-order EKF is constructed for identifying the magnetic flux. The state, input, and output vectors are described as follows:
(30)x=[i^di^qψ^mdψ^mq]T,u=[u^dLSu^qLS]T,y=[i^di^q]T.The system and output measurements are(31)ddtx=(-RSLSω^k0ω^kLS-ω^k-RSLS-ω^kLS000000000)x+(10010000)u,(32)y=(10000100)x.In order to simplify the description, (31) and (32) can be rewritten, respectively, as follows:
(33)ddtx=Lx+Mu,(34)y=hx.Considering the noise in process and measurement, approximate discrete equations (33) and (34) can be expanded as
(35)xk=(I+LT)xk-1+MTuk-1+σk-1,(36)yk=hxk+ρk,
where I is a unity matrix and T is a sample time. Equations (35) and (36) are equivalent to (20) and (21). Therefore, using (35) and (36) can realize (23)–(27) of the EKF recursive processing. The identified flux obtained through EKF is ψ^md and ψ^mq. The identified real flux is given as follows:
(37)ψ^m=(ψ^md)2+(ψ^mq)2.The accurate estimation of rotor position and speed in the flux-weakening area is obtained from (17) and (19) by using identified flux ψ^m. From (30)–(32), EKF for magnetic flux identification is a fourth-order linear system. After choosing [idiqψmdψmq]T as state variable and [idiq]T as output, the proposed magnetic flux identification system is linear to the input and output.In this part, SPMSM mathematical model based on EKF is induced, and the permanent magnet flux of SPMSM is identified by recurrent operation of EKF. The influence on uncertainties in magnetic magnet flux is eliminated from the position estimation.


## 4.1. Introduction to EKF

EKF presents the general problem of trying to estimate the statesx∈Rn that is dominated by the nonlinear stochastic difference equation as follows:
(20)xk=f(xk-1,uk,σk-1),
and the measurement output y∈Rn is as follows:
(21)yk=h(xk,ρk),
where uk is the control vector, the random variables σk-1 and ρk represent the process and measurement noise, and the covariance matrices of these noises are defined as
(22)Q=E[σ·σT],R=E[ρ·ρT].The solutions to EKF can be divided into two groups.(1)
Time update equation (predict), which uses the current state and error covariance to obtain the a priori estimations for the next time step.
(i)
Project the state ahead:(23)x^k-=f(x^k-1,uk,0).(ii)
Project the error covariance ahead:(24)Pk-=AkPk-1AkT+WkQk-1WkT.(2)
Measurement update equation (correct), which incorporates new measurements into the a priori estimations to obtain improved a posteriori estimations.
(i)
Compute EKF gain:(25)Kk=Pk-HkT(HkPk-HkT+VkRVkT)-1.(ii)
Update estimation with measurement:(26)x^k=x^k-+Kk(yk-h(x^k-,0)).(iii)
Update the error covariance:(27)Pk=(1-KkHk)Pk-,where x^k- is a priori state estimation at step k given knowledge of the process prior to step k, x^k is a posteriori state estimation at step k given measurement yk, and Pk- and Pk are the covariance of the a priori and a posteriori estimation errors.Ak is a Jacobian matrix of partial derivatives of f with respect to x; Ak[i,j]=(∂f[i]/∂x[j])(x^k-1,uk,0).Wk is a Jacobian matrix of partial derivatives of f with respect to σ; Wk[i,j]=(∂f[i]/∂σ[j])(x^k-1,uk,0).Hk is a Jacobian matrix of partial derivatives of h with respect to x; Hk[i,j]=(∂h[i]/∂x[j])(x^k,0).Vk is a Jacobian matrix of partial derivatives of h with respect to ρ; Vk[i,j]=(∂h[i]/∂ρ[j])(x^k,0).


## 4.2.  Permanent Magnet Flux Identification of SPMSM Based on EKF

Here, we present one method that identifies the magnetic flux change that can be generated during the motor dynamic. In Section3 of this paper, the complex number model of SPMSM is used for estimating the accurate position. In (17), the imaginary part of position estimator caused by the uncertainty of motor parameters is controlled to be zero by the error suppression PI controller. But, the demagnetization amount must not be controlled to be zero in the flux-weakening area. Calculating the demagnetization amount accurately is one of the most important problems to realize the accurate position estimation and vector control in the flux-weakening area. The complex number model composed of real part and imaginary part may cause the complexity in the flux identification. Therefore, the rotor flux is identified through the combination of SPMSM model (1) and EKF in this paper. Equation (1) with considering the magnetic flux change can be rewritten as follows:
(28)ddti^d=-RSi^dLS+i^qω^k+ω^kψ^mqLS+u^dLS,ddti^q=-RSi^qLS-i^dω^k-ω^kψ^mdLS+u^qLS.In order to estimate the magnetic flux, the magnetic flux is chosen as one state variable. Because the magnetic flux linkage cannot change sharply during the motor dynamic [1], its differential is assumed to be zero:
(29)dψ^mddt=0,dψ^mqdt=0.By choosing set of the electric current variables and magnetic flux of theq axis as the state variables, a fourth-order EKF is constructed for identifying the magnetic flux. The state, input, and output vectors are described as follows:
(30)x=[i^di^qψ^mdψ^mq]T,u=[u^dLSu^qLS]T,y=[i^di^q]T.The system and output measurements are(31)ddtx=(-RSLSω^k0ω^kLS-ω^k-RSLS-ω^kLS000000000)x+(10010000)u,(32)y=(10000100)x.In order to simplify the description, (31) and (32) can be rewritten, respectively, as follows:
(33)ddtx=Lx+Mu,(34)y=hx.Considering the noise in process and measurement, approximate discrete equations (33) and (34) can be expanded as
(35)xk=(I+LT)xk-1+MTuk-1+σk-1,(36)yk=hxk+ρk,
where I is a unity matrix and T is a sample time. Equations (35) and (36) are equivalent to (20) and (21). Therefore, using (35) and (36) can realize (23)–(27) of the EKF recursive processing. The identified flux obtained through EKF is ψ^md and ψ^mq. The identified real flux is given as follows:
(37)ψ^m=(ψ^md)2+(ψ^mq)2.The accurate estimation of rotor position and speed in the flux-weakening area is obtained from (17) and (19) by using identified flux ψ^m. From (30)–(32), EKF for magnetic flux identification is a fourth-order linear system. After choosing [idiqψmdψmq]T as state variable and [idiq]T as output, the proposed magnetic flux identification system is linear to the input and output.In this part, SPMSM mathematical model based on EKF is induced, and the permanent magnet flux of SPMSM is identified by recurrent operation of EKF. The influence on uncertainties in magnetic magnet flux is eliminated from the position estimation.


## 5. System Control Structure

The simulation block of sensorless SPMSM drive system by using the proposed position/speed estimation and EKF based permanent magnet flux identification is shown in Figure3. The rotor-flux-oriented vector control is employed as the basic control strategy for SPMSM drive. The reference voltage vectors are obtained by current proportional integral (PI) regulators. As shown in Figure 3, “space vector pulsewidth modulation (SVPWM)” module realizes SVPWM or overmodulation algorithm according to the amplitude of the reference voltage vector. “Sampling” module detects only 3-phase motor voltage for estimating the position and speed of sensorless control. “Position estimation, flux identification” module identifies the permanent magnet flux by EKF and estimates the position and speed by the proposed estimation method. “Flux identification using EKF” module identifies the permanent magnet flux, whose result (ψ^m) is sent to “position estimation” module (17) for rotor position estimation. In this system, the rotor position estimation by using the combination of motor complex number model based position estimation and position estimation error suppression PI controller and the permanent magnet flux identification by EKF are very interrelated to each other. The d-q axis reference voltages, currents, and estimated speed are used in flux identification. The detailed contents of “flux identification using EKF” block are shown in Figure 4.Figure 3
Simulation block of sensorless SPMSM drive system by using proposed position/speed estimation and flux identification.Figure 4
Simulation block of identifying PM flux by EKF.


## 6. Simulation Study

In order to verify the proposed control strategy, a model of pulse width modulation (PWM) inverter-fed SPMSM drive system is built in MATLAB9.0/SIMULINK platform. Simulation studies are carried out with a 10 kHz control frequency. The parameters of SPMSM used in the simulation are listed in Table1.Table 1
Parameters for the motor and inverter systems in simulations.
Parameters
Value
d axis inductance LS
8.5
e
-
3 H
Stator resistanceRS
2.875Ω
Back-EMF coefficientKE
0.2256 V/(rad/s)
Pole pairsp
4
DC link voltage Vdc
540 V
Load torqueTL
10 Nm
Maximum allowable current of  motor and inverter Iqmax⁡
10 A
Inertia momentl
0.8
e
-
3 Nm2
Magnet flux linkageψm
0.175 wb
Motor rated speedω
1000 rpmFigure5 shows the error between the estimated position by proposed estimation method and the real position from position sensor under no load condition. It is clear that the real position is consistent with the estimated position which is obtained by the proposed estimation method when the motor speed comes to steady state. Figure 6 shows that the real flux and the flux that is identified by EKF are well matched. Figure 7 shows the d axis and q axis values of the identified flux ψ^m. We know that ψmd converges to ψ^m and ψmq converges to 0. It is clear that the proposed permanent magnet flux identification scheme is correct.Figure 5
Error between actual position and estimated position.Figure 6
Wave of assumed flux and identified flux.Figure 7
Rotor flux ind-q frame.As mentioned above, the SPMSM drive system that includes permanent magnet flux identification by EKF exhibits a very good steady-state performance. Simulations are conducted to examine the dynamic performance of the proposed identification scheme in the case of load torque change. Figure8 shows the estimation error of rotor speed (the error between the estimated speed and the real speed from the rotor position sensor) and Figure 9 shows the estimated speed under the condition of load torque change. Figure 10 shows the real flux and the identified flux when the load torque is changed. Figure 11 shows the d-q axis current curves.Figure 8
Speed estimation error with load torque from 2 Nm to 4 Nm.Figure 9
Speed estimation with load torque from 2 Nm to 4 Nm.Figure 10
Flux identification with load torque from 2 Nm to 4 Nm.Figure 11
d
-
q frame currents with load torque from 2 Nm to 4 Nm.Figure12 shows the speed estimation curves obtained from the back-EMF based estimator [6, 11] and our proposed estimator. The motor is accelerated from speed 0 to speed 1200 r/min⁡. Estimation time is about 0.28 s in the case of using back-EMF based method and is about 0.25 s in the case of applying our proposed method. In the steady state, vibrations are about 2% when back-EMF based estimator is used and are about 0.4% when our method is used. The performance of our proposed method is quite similar to the reference curve compared to back-EMF based method. Figure 13 shows the assumed flux and identified flux in the demagnetization state. Figure 14 shows the position estimation error curves obtained from the back-EMF based estimator [6, 11] and our proposed estimator. The position estimation error generated from the proposed position estimation method is reduced to 25% with respect to when the back-EMF based method is used. The simulation results show that both estimating the rotor position and speed by our proposed method and identifying the flux by EKF are done effectively.Figure 12
Estimation response of step speed command.Figure 13
Identified flux in the demagnetization state.Position estimation error. (a) Back-EMF based estimator, (b) proposed estimator.
(a)
(b)


## 7. Conclusion

In this paper, the rotor position angle is estimated by a novel position estimator based on complex number model of SPMSM. By using an error suppression PI controller that is combined with imaginary part of estimator, the estimation error of rotor position angle caused by the uncertainty of system parameters is converged to zero quickly. By identifying the permanent magnet flux by EKF, the sensitivity to uncertainty of permanent magnet flux that may be generated in the flux-weakening region is eliminated. Simulation results demonstrated that both the estimation of the rotor position and the accurate identification of permanent magnet flux are realized with satisfying performance. The estimation error of rotor position is limited within 0.4%, and the identification error of permanent magnet flux is realized with 0.8%. The step response rate has decreased about 10% and the position estimation error is reduced to 25% compared to using the back-EMF based position estimation method under the same control condition.



---

*Source: 480640-2014-08-03.xml*