Next Article in Journal
Optimal Planning of Sustainable Buildings: Integration of Life Cycle Assessment and Optimization in a Decision Support System (DSS)
Previous Article in Journal
Insights on Energy Transitions in Mexico from the Analysis of Useful Exergy 1971–2009
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

State Estimation of Permanent Magnet Synchronous Motor Using Improved Square Root UKF

1
The School of Electrical and Information Engineering, Jiangsu University, Zhenjiang 212013, China
2
Key Laboratory of Facility Agriculture Measurement and Control Technology and Equipment of Machinery Industry, Jiangsu University, Zhenjiang 212013, China
*
Author to whom correspondence should be addressed.
Energies 2016, 9(7), 489; https://doi.org/10.3390/en9070489
Submission received: 9 April 2016 / Revised: 18 May 2016 / Accepted: 20 June 2016 / Published: 24 June 2016

Abstract

:
This paper focuses on an improved square root unscented Kalman filter (SRUKF) and its application for rotor speed and position estimation of permanent magnet synchronous motor (PMSM). The approach, which combines the SRUKF and strong tracking filter, uses the minimal skew simplex transformation to reduce the number of the sigma points, and utilizes the square root filtering to reduce computational errors. The time-varying fading factor and softening factor are introduced to self-adjust the gain matrices and the state forecast covariance square root matrix, which can realize the residuals orthogonality and force the SRUKF to track the real state rapidly. The theoretical analysis of the improved SRUKF and implementation details for PMSM state estimation are examined. The simulation results show that the improved SRUKF has higher nonlinear approximation accuracy, stronger numerical stability and computational efficiency, and it is an effective and powerful tool for PMSM state estimation under the conditions of step response or load disturbance.

1. Introduction

The extended Kalman filter (EKF) has been successfully implemented as a state observer for induction motor (IM) drives in various areas [1,2,3,4,5,6,7,8]. However, EKF needs to calculate the nonlinear equation of the Jacobian matrix, which is sub-optimal and can easily lead to divergence. In Ref. [9], Unscented Kalman filter yielded performance equivalent to the Kalman filter for linear systems, yet generalized elegantly to nonlinear systems without the linearization steps required by the EKF. Also in this paper, a symmetric sigma point solution, which used 2n + 1 points to match the mean and covariance of an n-dimensional random variable, was presented. With this set of points, the unscented transform guaranteed the same performance as the truncated second order filter, with the same order of calculations as an EKF but without the need to calculate any approximations or derivatives. In real-time applications, it is critical that both the computational costs and storage requirements are minimized. In [10], the minimal skew simplex points were introduced, which used the minimum number of sigma points, and reduced the computation time. In [11], UKF was used to estimate the synchronous generator state. Compared with EKF, this method can effectively improve the precision of state estimation. However, covariance asymmetry and non-positive definite defects all exist in both EKF and UKF. A square root unscented Kalman filter (SRUKF) algorithm solves the problem of filtering divergence caused by non-positive of error covariance matrix in general EKF and UKF, and the stability of the algorithm is improved. In [12], SRUKF was used to estimate the speed and rotor space of permanent magnet synchronous motor, and effectively overcame the numerical calculation defects. The problems of bad robustness for the model parameter change, slow convergence, and lower tracking ability to abrupt state still exist. Strong tracking filter (STF) can adjust the state prediction error covariance matrix and filter gain matrix, which also has a strong robustness of the model mismatch and a unique strong tracking ability. It can effectively improve the tracking performance at a steady and abrupt state [13].
Combining the ordinary SRUKF and strong tracking filtering, this paper presents an improved SRUKF for rotor speed estimation of permanent magnet synchronous motor (PMSM) sensorless drives. On the one hand, the time-varying fading factor and softening factor based on STF theory are introduced to self-adjust the gain matrices and the state forecast covariance square root matrix, which can realize the residuals orthogonality and force the SRUKF track the real state rapidly. On the other hand, Cholesky and QR decomposition are introduced in SRUKF, which can effectively avoid filter divergence caused by the error covariance matrix negative, and improve the convergence speed and stability. The vector control system for PMSM without a speed sensor is set up based on this approach. Rotor position and speed estimators are designed by this method. Simulation results confirm the validity of the proposed approach.
This paper is organized as follows. Section 2 presents the SRUKF. Section 3 discusses the design of improved SRUKF. The PMSM model and observer are described in Section 4, and Section 5 shows the simulation results. Conclusions and future work are given in Section 6.

2. Square-Root UKF

Consider the following a general discrete-time nonlinear dynamic system:
{ x k + 1 = f ( x k , u k ) + w k y k = h ( x k , u k ) + v k
where f(,) and h(,) denote the nonlinear function with one-order continuous partial derivative. xk represents the unobserved state of the system. uk is a known exogenous input. yk is the observed measurement signal. wk and vk are uncorrelated from each other with zero-mean and covariance, respectively. Qk and Rk are the noise covariance of the process and measurement, respectively.

2.1. The Minimal Skew Simplex Points

The minimal skew sigma points are a simplex set, which is chosen to match the mean, covariance and minimize their skew. These points have the important property that they can be constructed recursively. The point selection algorithm for the simplex unscented transform is as follows:
(1)
Choose the initial weight value:
0 W 0 1
(2)
Choose weight sequence:
W i = { 1 W 0 2 n for   i = 1 W 1 for   i = 2 2 i 1 W 1 for   i = 3 , , n + 1
(3)
Initialize the vector sequences of sigma points as:
χ 0 1 = [ 0 ] ,     χ 1 1 = [ 1 2 W 1 ]   and   χ 2 1 = [ 1 2 W 1 ]
(4)
Expand vector sequences for j = 2,…,n, according to:
χ i j + 1 = { [ χ 0 j 0 ] for   i = 0 [ χ i j 1 2 W j ] for   i = 1 , , j [ 0 j 1 2 W j ] for   i = j + 1

2.2. Square-Root UKF for State-Estimation

The square-root UKF makes use of three powerful linear algebra techniques, QR decomposition, Cholesky factor updating and efficient least squares. The algorithmic description of SRUKF is as follows [12]:
(1)
Initialize with
x ^ 0 = E [ x 0 ]
P 0 = E [ ( x 0 x ^ 0 ) ( x 0 x ^ 0 ) T ]
S 0 = c h o l   ( P 0 )
where x ^ 0 denotes an initial guess for the initial state x 0 , and P 0 is the initial estimation variance. Chol (·) represents the matrix square-root using lower triangle Cholesky decomposition.
(2)
Sigma point calculation and time update
χ k 1 = x ^ k 1 + S k 1 χ k 1
χ ^ k / k 1 = f ( χ i , k 1 , u k 1 )
x ^ k = i = 0 2 n W i ( m ) χ ^ i , k / k 1
S k = q r { [ W 1 ( c ) ( χ ^ 1 : 2 n , k / k 1 x ^ k )   Q k ] }
S k = c h o l u p d a t e   { S k , χ ^ 0 , k x ^ k , W 0 ( c ) ] }
where qr{·} and cholupdate{·} represents QR decomposition and Cholesky factor updating.
y ^ k / k 1 = h ( χ ^ k / k 1 )
y k = i = 0 2 n W i ( m ) y ^ k / k 1
(3)
Measurement update
S y k = q r { [ W 1 ( c ) [ y ^ 1 : 2 n , k y k ]   R k ] }
S y k = c h o l u p d a t e { [ S y k , y ^ 0 , k y k , W 0 ( c ) ] }
P x k y k = i = 0 2 n W i ( c ) { χ ^ i , k / k 1 x ^ k } { y ^ i , k / k 1 y k } T
K k / k = P x k y k P y k 1 = ( P x k y k / S y k T ) / S y k
x ^ k / k = x ^ k + K k [ y k y k ]
S k = c h o l u p d a t e { S k , K k / k S y k , 1 }

3. Improved SRUKF

SRUKF has the benefit of numerical stability, but this method still has bad robustness, slow tracking capacity, and low convergence to the model parameters variation and abrupt state. So, combined with strong tracking filter, improved SRUKF filtering algorithm is presented. The orthogonality conditions, which can make the SRUKF filter become a strong tracking filter [14], are given as follows:
E [ x k x ^ k / k ] [ x k x ^ k / k ] T = M I N
E [ γ k + j γ k T ] = 0
where MIN in Equation (22) means the minimum mean square error and γ k = y k y k .
SRUKF can satisfy the requirements of Equation (22). In order to satisfy Equation (23), selection principle of the gain matrix Kk/k is deduced.
Lemma 1. 
Define ε ( k ) = x k x ^ k , x ^ k is the estimated state with the improved SRUKF algorithm, if O [ | ε ( k ) | 2 ] < < O [ | ε ( k ) | ] , then the following equation exist.
E ( γ k + j γ k T ) H ( x ^ k + j / k + j 1 ) F ( u k + j 1 , x ^ k + j 1 / k + j 1 ) × ( I K k + j 1 H k + j 1 ) F ( u k + j 2 , x ^ k + j 1 / k + j 2 ) × × ( I K k + 1 H k + 1 ) F ( u k , x ^ k / k ) ( P x k y k K k / k C k )
where Ck = Ek γkT). H and F are Jacobian matrixs of h(k, xk) and f(k, xk, uk), respectively.
Proof of Lemma 1. 
E ( γ k + j γ k T ) = E [ ( y k + j y k + j ) γ k T ] = E { [ y k + j i = 0 2 n W i ( m ) h ( χ ^ k + j / k + j 1 ) ] γ k T } E { [ H ( x ^ k + j / k + j 1 ) ( x k + j i = 0 2 n W i ( m ) χ ^ k + j / k + j 1 ] γ k T } E { [ H ( x ^ k + j / k + j 1 ) ( x k + j x ^ k + j / k + j 1 ) ] γ k T } E { [ H ( x ^ k + j / k + j 1 ) ( F ( u k + j 1 , x ^ k + j 1 / k + j 1 ) × ( x k + j 1 x ^ k + j 1 / k + j 1 ) + w k + j 1 ) ] γ k T } E { [ H ( x ^ k + j / k + j 1 ) ( F ( u k + j 1 , x ^ k + j 1 / k + j 1 ) × ( x k + j 1 x ^ k + j 1 / k + j 2 K k + j 1 H k + j 1 ) × ( x k + j 1 x ^ k + j 1 / k + j 2 ) ] γ k T } E { [ H ( x ^ k + j / k + j 1 ) ( F ( u k + j 1 , x ^ k + j 1 / k + j 1 ) × ( I K k + j 1 H k + j 1 ) ( x k + j 1 x ^ k + j 1 / k + j 2 ) γ k T } E { [ H ( x ^ k + j / k + j 1 ) F ( u k + j 1 , x ^ k + j 1 / k + j 1 ) × ( I K k + j 1 H k + j 1 ) F ( u k + j 2 , x ^ k + j 1 / k + j 2 ) × × ( I K k + 1 H k + 1 ) F ( u k , x ^ k / k ) ( x k x ^ k / k ) ] γ k T } [ H ( x ^ k + j / k + j 1 ) F ( u k + j 1 , x ^ k + j 1 / k + j 1 ) × ( I K k + j 1 H k + j 1 ) F ( u k + j 2 , x ^ k + j 1 / k + j 2 ) × × ( I K k + 1 H k + 1 ) F ( u k , x ^ k / k ) ] E [ ( x k x ^ k / k ) γ k T ]
where
E [ ( x k x ^ k / k ) γ k T ] = E [ ( x k x ^ k / k 1 K k / k γ k ) ( y k y k ) T ] = E [ ( x k x ^ k / k 1 ) ( y k y k ) T K k / k γ k γ k T ] = P x k y k K k / k C k
So Lemma 1 is satisfied. In the proof, the noise and signal is statistical independence, and the three times and above error are ignored.
In order to weaken the effect of old data, and achieve fast tracking to the abrupt status, time-varying fading factor λk is introduced as follows
K k / k = P x k y k ( λ k P y k ) 1
where λ k = { λ k , λ k > 1 1 , λ k 1 .
From Lemma 1, if P x k y k K k / k C k = 0 , i.e., I ( λ k P y k ) 1 C k = 0 , then the following equation can be deduced.
λ k S y k S y k T = C k
Solving the trace of Equation (27),
t r ( λ k S y k S y k T ) = t r C k
In order to avoid excessive regulation of λk, and achieve more smooth of the state estimation, softening factor η is introduced, then Equation (28) becomes
λ k = t r ( C k η R k ) t r ( S y k S y k T )
where C k = { γ 0 γ 0 T k = 0 ρ C k 1 + γ k γ k T 1 + ρ k 1 , ρ is forgetting factor, 0 < ρ ≤ 0.95. Generally, ρ = 0.95. When the model is more accurate, λk = 1, the improved SRUKF degrades into SRUKF.

4. Mathematical Model of PMSM and System Implementation

4.1. Mathematical Model of PMSM

The mathematical model of PMSM in α–β coordinate system can be represented as follows
[ d i α d t d i β d t d ω d t d θ d t ] = [ R L d i α + ψ f L d P 1 ω sin   θ R L q i β ψ L q P 1 ω cos   θ 1 J ( T e F ω T m ) P 1 ω ] + [ 1 L d 0 0 1 L d 0 0 0 0 ] [ u α u β ]
y = [ i α i β ] T = [ 1 0 0 0 0 1 0 0 ] [ i α i β ω θ ]
where uα, uβ and iα, iβ are the stator voltage and current respectively. Ψf is magnet flux linkage. Ld and Lq are the stator inductance in dq coordinate system, respectively. F is the friction coefficient of rotor and load. R is the stator resistance. Te and Tm are the electromagnetic torque and mechanical torque, respectively. P1 is the number of pole pairs in torque winding. J is the total moment of inertia of the mechanical system. θ is the angle of rotor position. ω is the rotor speed. Assuming the sampling period Ts, the Equations (30) and (31) are discretized. We set the state variable Xk = [i,k i,k ωk θk]T, control variables Uk = [u,k u,k]T, output variables Yk = [i,k i,k]T. From the discretization Formula (32),
X ˙ 1 T s ( X k + 1 X k )
Equations (30) and (31) can be transformed to
X k + 1 = [ i α , k + 1 i β , k + 1 ω k + 1 θ k + 1 ] = [ ( 1 T s R L d ) i α , k + T s ψ f L d P 1 ω k sin θ k ( 1 T s R L q ) i β , k T s ψ f L q P 1 ω k cos θ k + [ T s L d 0 0 T s L d 0 0 0 0 ] [ u α , k u β , k ] ( 1 T s J F ) ω k + T s J ( T e T m ) T s P 1 ω k + θ k ]
Y k = [ i α , k i β , k ] T = [ 1 0 0 0 0 1 0 0 ] [ i α , k i β , k ω θ ]
H k = [ 1 0 0 0 0 1 0 0 ]

4.2. System Implementation

The basic configuration of speed sensorless vector control system for PMSM based on the improved SRUKF is shown in Figure 1. The space vector pulse width modulation (SVPWM) inverter is used. Parameters of PMSM used in simulation experiment are shown in Table 1.
The simulation solver mode is ode3 with fixed step. Simulation time is 0.1s. The following values are chosen in simulation: x0 = 0, P0 = diag(0.2, 0.2, 180, 20), W0 = 0.25, Qk = diag(10−6, 10−6, 10−2, 10−5), Rk = diag(0.1, 0.1), η = diag(3.2, 3.2).

5. Simulation Results and Error Analysis

5.1. Simulation Results

The performance of the improved SRUKF observer is validated subject to a speed command. Experimental results are shown in Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11, Figure 12, Figure 13, Figure 14, Figure 15 and Figure 16. In which, “Real” indicates the measurement speed, “SRUKF”, and “improved SRUKF” indicate the speed or space position estimation and error with the SRUKF and improved SRUKF, respectively. Δθ=θ− θ ^ , θ is the measurement space position of the rotor, θ ^ is the estimation space position, Δθ is the space position estimation error. In Figure 2, the rotor speed ω reference is 50 rad/s, during this test, the motor is loaded at 3 N·m. The estimated and measured speeds match well, which indicates that the observer provides a good estimate at low speeds. The corresponding space position estimation and error are shown in Figure 3 and Figure 4. The error of improved SRUKF is smaller than SRUKF.
In Figure 5, the drive runs at a reference speed of 800 rad/s. In Figure 8, the reference speed is 500 rad/s. The reference speed is abruptly changed to 100 rad/s at 0.04 s, and then, the reference speed is abruptly returned to 500 rad/s. In Figure 11, the drive runs at a reference speed of 700 rad/s, the disturbance is loaded at 0.04 s–0.06 s. The corresponding space position estimation and error are shown in Figure 6, Figure 7, Figure 9, Figure 10, Figure 12 and Figure 13, respectively.
Considering the change of motor physical parameters during PMSM operation, approximate parameter values are used in the observer. Waveforms for operation with detuned stator and rotor resistances are shown in Figure 14, Figure 15 and Figure 16, respectively. The stator resistance and d, q axis inductance are detuned by increasing 25%, with respect to the values in the Table 1. The results show that the observer remains stable and can handle well.
The experimental results show that the estimated speed follows the actual speed closely with rapid dynamics, and small estimation error of rotor space is achieved. It indicates that the improved SRUKF is an accurate and fast state estimator. Also, compared with the ordinary SRUKF, the improved SRUKF exhibits a strong robustness and estimation accuracy under the conditions of step response or load disturbance.

5.2. Error Analysis

Set the root mean square error (RMSE) as the system estimation quality evaluation criteria
R M S E ( x ) = E [ | x x ^ | 2 ] = 1 N i = 1 N | x x ^ | 2
where x is the actual measurement. x ^ is an estimation. N is the sampling number.
Table 2 and Table 3 show the RMSE of speed estimation and rotor space position estimation, respectively. From the tables, it can be found that system performance is little difference between SRUKF and improved SRUKF under the conditions of steady state. However, under the conditions of step response or load disturbance, speed estimation error of improved SRUKF reduced about 55% more than SRUKF, and space position estimation error reduced about 65%.
Figure 17 shows the values of the time-varying fading factor at speed accelerated from 500 rad/s to 100 rad/s, then to 500 rad/s. Figure 18 shows the values of the time-varying fading factor at parameters disturbance. Time-varying fading factor is adaptive, which can adjust the filter gain and error covariance square root matrix. During steady-state operation at a constant speed, time-varying fading factor equals 1, and the improved SRUKF filter degrades to ordinary SRUKF.

6. Conclusions

This paper has proposed and investigated an improved SRUKF filter for state estimation in sensorless PMSM drives. The main contribution is the combination of the SRUKF and a strong tracking filter. The SRUKF reduces the computational errors by propagating the SR of matrices instead of the matrices themselves. In order to realize the residuals orthogonality and force the SRUKF filter to track the real state rapidly, the time-varying fading factor and softening factor are introduced to self-adjust gain matrices and the state forecast covariance square root matrix.
A speed sensorless vector control system for PMSM based on the improved SRUKF is implemented. The simulation results illustrate that the proposed method has higher nonlinear approximation accuracy, stronger numerical stability, and computational efficiency. Strong robustness is achieved under the conditions of low speed, high speed, step response, and load disturbance. It can achieve a precise estimate of the speed, space position, and ensure the rotor suspends stably. Going forward, learning how to apply this algorithm on an actual system based on DSP will be the next project.

Acknowledgments

This work is supported by Natural Science Foundation of Jiangsu Province (grant number BK20150530), Jiangsu Province Colleges and Universities Natural Science Research Project (grant number 12KJB470002), A Project Funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD), and the Professional Research Foundation for Advanced Talents of Jiangsu University (grant number 14JDG077).

Author Contributions

All authors have worked on this manuscript together and all authors have read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflicts of interest.

Nomenclature

P1Number of poles pairs in torque winding
RStator resistance
Ld (Lq)Stator inductance in dq coordinate system
ψfMagnet flux linkage
JRotor moment of inertia
MRotor mass
VRated voltage of torque winding
IRated current
TeElectromagnetic torque
TMechanical time constant
TmRated mechanical torque
TsSampling period
θRotor angle position
ωRotor speed
FFriction coefficient
uα (uβ)Stator voltage in α–β coordinate system
iα (iβ)Stator current in α–β coordinate system

References

  1. Kalman, R.E. A new approach to linear filtering and prediction problems. J. Basic Eng. 1960, 82, 35–45. [Google Scholar] [CrossRef]
  2. Alonge, F.; D’Ippolito, F.; Sferlazza, A. Sensorless control of induction-motor drive based on robust Kalman filter and adaptive speed estimation. IEEE Trans. Ind. Electron. 2014, 61, 1444–1453. [Google Scholar] [CrossRef]
  3. Salvatore, N.; Caponio, A.; Neri, F.; Stasi, S.; Cascella, G.L. Optimization of delayed-state Kalman-filter-based algorithm via differential evolution for sensorless control of induction motors. IEEE Trans. Ind. Electron. 2010, 57, 385–394. [Google Scholar] [CrossRef]
  4. Masi, A.; Butcher, M.; Martino, M.; Picatos, R. An application of the extended Kalman filter for a sensorless stepper motor drive working with long cables. IEEE Trans. Ind. Electron. 2012, 59, 4217–4225. [Google Scholar] [CrossRef]
  5. Laamari, Y.; Chafaa, K.; Athamena, B. Particle swarm optimization of an extended Kalman filter for speed and rotor flux estimation of an induction motor drive. Electr. Eng. 2015, 97, 129–138. [Google Scholar] [CrossRef]
  6. Bolognani, S.; Oboe, R.; Zigliotto, M. Sensorless full-digital PMSM drive with EKF estimation of speed and rotor position. IEEE Trans. Ind. Electron. 1999, 46, 184–191. [Google Scholar] [CrossRef]
  7. Lu, K.Y.; Lei, X.; Blaabjerg, F. Artificial inductance concept to compensate nonlinear inductance effects in the back EMF-based sensorless control method for PMSM. IEEE Trans. Energy Convers. 2013, 28, 593–600. [Google Scholar] [CrossRef]
  8. Kumar, S.; Prakash, J.; Kanagasabapathy, P. A critical evaluation and experimental verification of Extended Kalman Filter, Unscented Kalman Filter and Neural State Filter for state estimation of three phase induction motor. Appl. Soft Comput. 2011, 11, 3199–3208. [Google Scholar] [CrossRef]
  9. Julier, J.S.; Uhlmann, J.; Whyte, H.D. A new method for the nonlinear transformation of means and covariance in filters and estimators. IEEE Trans. Autom. Control 2000, 45, 477–482. [Google Scholar] [CrossRef]
  10. Julier, J.S.; Uhlmann, J. Reduced sigma point filters for the propagation of means and covariances through nonlinear transformations. In Proceedings of the American Control Conference, Anchorage, AK, USA, 8–10 May 2002; pp. 887–892.
  11. Ghahremani, E.; Kamwa, I. Online state estimation of a synchronous generator using unscented Kalman filter from phasor measurements units. IEEE Trans. Energy Convers. 2011, 26, 1099–1108. [Google Scholar] [CrossRef]
  12. Jafarzadeh, S.; Lascu, C.; Fadali, M.S. Square root unscented Kalman filters for state estimation of induction motor drives. IEEE Trans. Ind. Appl. 2013, 49, 92–99. [Google Scholar] [CrossRef]
  13. He, X.; Wang, Z.D.; Wang, X.F. Networked strong tracking filtering with multiple packet dropouts: Algorithms and applications. IEEE Trans. Ind. Electron. 2014, 61, 1454–1463. [Google Scholar] [CrossRef]
  14. Zhou, D.H.; Ye, Y.Z. Modern Fault Diagnosis and Fault-Tolerant Control; Tsinghua University Press: Beijing, China, 2000. [Google Scholar]
Figure 1. Configuration of the vector control system based on improved SRUKF.
Figure 1. Configuration of the vector control system based on improved SRUKF.
Energies 09 00489 g001
Figure 2. Speed estimation at 50 rad/s.
Figure 2. Speed estimation at 50 rad/s.
Energies 09 00489 g002
Figure 3. Space position estimation at 50 rad/s.
Figure 3. Space position estimation at 50 rad/s.
Energies 09 00489 g003
Figure 4. Space position estimation error at 50 rad/s.
Figure 4. Space position estimation error at 50 rad/s.
Energies 09 00489 g004
Figure 5. Speed estimation at 800 rad/s.
Figure 5. Speed estimation at 800 rad/s.
Energies 09 00489 g005
Figure 6. Space position estimation at 800 rad/s.
Figure 6. Space position estimation at 800 rad/s.
Energies 09 00489 g006
Figure 7. Space position estimation error at 800 rad/s.
Figure 7. Space position estimation error at 800 rad/s.
Energies 09 00489 g007
Figure 8. Speed estimation at speed step response.
Figure 8. Speed estimation at speed step response.
Energies 09 00489 g008
Figure 9. Space position estimation at speed step response.
Figure 9. Space position estimation at speed step response.
Energies 09 00489 g009
Figure 10. Space position estimation error at speed step response.
Figure 10. Space position estimation error at speed step response.
Energies 09 00489 g010
Figure 11. Speed estimation at load disturbance.
Figure 11. Speed estimation at load disturbance.
Energies 09 00489 g011
Figure 12. Space position estimation at load disturbance.
Figure 12. Space position estimation at load disturbance.
Energies 09 00489 g012
Figure 13. Space position estimation error at load disturbance.
Figure 13. Space position estimation error at load disturbance.
Energies 09 00489 g013
Figure 14. Speed estimation at parameters disturbance.
Figure 14. Speed estimation at parameters disturbance.
Energies 09 00489 g014
Figure 15. Space position estimation error at parameters disturbance.
Figure 15. Space position estimation error at parameters disturbance.
Energies 09 00489 g015
Figure 16. Space position estimation error at parameters disturbance.
Figure 16. Space position estimation error at parameters disturbance.
Energies 09 00489 g016
Figure 17. Fading factors of speed step response.
Figure 17. Fading factors of speed step response.
Energies 09 00489 g017
Figure 18. Fading factors of parameters disturbance.
Figure 18. Fading factors of parameters disturbance.
Energies 09 00489 g018
Table 1. Parameter of PMSM.
Table 1. Parameter of PMSM.
ParameterSymbolValue
Number of poles pairs in torque windingP14
Stator resistanceR4.025 Ω
d-axis inductanceLd11.9 mH
q-axis inductanceLq11.9 mH
Magnet flux linkageψf0.245 Wb
Rotor inertiaJ1.0 × 10−4 kg·m2
Rotor massm1.2 kg
Rated voltage of torque windingV240 V
Rated currentI5.86 A
The maximum electromagnetic torqueTe18.82 N·m
Mechanical time constantT0.25 s
Rated load torqueTm9.41 N·m
Speed rangeω0–1000 rad/s
Table 2. RMSE of speed estimation.
Table 2. RMSE of speed estimation.
MethodSRUKFImproved SRUKF
Speed
50 rad/s6.78656.1343
800 rad/s61.453225.1143
700 rad/s load disturbance55.727025.3146
500 rad/s to 100 rad/s to 500 rad/s71.621931.5823
parameters disturbance55.204243.0836
Table 3. RMSE of space position estimation.
Table 3. RMSE of space position estimation.
MethodSRUKFImproved SRUKF
Speed
50 rad/s0.00660.0053
800 rad/s0.05110.0133
700 rad/s load disturbance0.04580.0156
500 rad/s to 100 rad/s to 500 rad/s0.15410.0187
parameters disturbance0.10840.0451

Share and Cite

MDPI and ACS Style

Xu, B.; Mu, F.; Shi, G.; Ji, W.; Zhu, H. State Estimation of Permanent Magnet Synchronous Motor Using Improved Square Root UKF. Energies 2016, 9, 489. https://doi.org/10.3390/en9070489

AMA Style

Xu B, Mu F, Shi G, Ji W, Zhu H. State Estimation of Permanent Magnet Synchronous Motor Using Improved Square Root UKF. Energies. 2016; 9(7):489. https://doi.org/10.3390/en9070489

Chicago/Turabian Style

Xu, Bo, Fangqiang Mu, Guoding Shi, Wei Ji, and Huangqiu Zhu. 2016. "State Estimation of Permanent Magnet Synchronous Motor Using Improved Square Root UKF" Energies 9, no. 7: 489. https://doi.org/10.3390/en9070489

APA Style

Xu, B., Mu, F., Shi, G., Ji, W., & Zhu, H. (2016). State Estimation of Permanent Magnet Synchronous Motor Using Improved Square Root UKF. Energies, 9(7), 489. https://doi.org/10.3390/en9070489

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop