Abstract
In this paper, in order to solve the problem of mismatch between the noise covariance matrix of the traditional Kalman filter and the actual system, a Sage-Husa noise estimator is incorporated to adaptively tune the system and process noise matrices, thus improving the robustness and convergence speed. An extended inverse potential observer based on Kalman filtering is used to construct the PMSM full-order state equation, which combines NQPLL and LPF filtering techniques. A modification of the Sage-Husa inertial weight equation introduces stabilization coefficients to ensure the semi-positive character of the noise covariance matrix, thus improving the system stability. In addition, sliding mode control is integrated into the velocity loop to improve immunity and response speed, while the high torque per ampere module optimizes the torque-to-current ratio control to minimize the motor losses. Simulation and experimental results show that the proposed observer achieves high estimation accuracy and adaptability in all velocity domains.
Similar content being viewed by others
Introduction
Permanent magnet synchronous motors (PMSM) are increasingly recognized for their wide range of applications in industrial automation and renewable energy1. This recognition stems from its high power output, low energy consumption, excellent performance, long lifetime and minimal noise level2. To address the challenges faced in optimizing the performance of PMSMs, in the field of position sensorless control3,4 such as Extended Kalman Filter (EKF)5,6, Model-Referenced Adaptive Control (MRAC)7, Sliding Mode Observer (SMO)8, and Nonlinear Magnetic Chain Observer (NMCO)9,10 have been explored by an increasing number of scholars.
Compared to the Lomborg sliding film observer, the EKF exhibits superior convergence performance at low speeds, although its effectiveness relies heavily on the precise tuning of the noise matrix parameters, which often requires extensive simulation and calibration work. Bolognani et al. proposed an all-digital position sensor-less drive technique that is reduced to be limited to current and DC link voltage sensors, which greatly reduces the use of sensors and improves the reliability and economy of the system11. Barut M et al. implemented a two-input extended Kalman filter estimator that simultaneously estimates the stator resistance, rotor resistance, and load torque variations based on two extended induction motor models12. Kim et al. proposed a high performance embedded permanent magnet synchronous motor drive without rotary position sensor using a reduced order EKF with lower order state variables greatly reducing the computational load and computational error13.Bolognani S et al. combined a novel self-tuning process to bring the EKF drive closer to the actual product, with an algorithm based on full normalization14.
Studies15,16 have utilized real-coded genetic algorithms to augment sensorless control systems for PMSMs, resulting in the development of a fully normalized EKF algorithm that improves the adaptability of various control systems. However, these optimization techniques are very complex and require a large amount of training data and long training time, posing challenges for practical engineering applications. Literature17 introduced fuzzy logic system into EKF for estimating the noise covariance matrices Q and R in sensorless DC brushless control systems, but the control rules of the fuzzy logic system will directly affect the estimation accuracy of KF. Literature18 proposed a new interest-based method for estimating the noise covariance matrix of a system, which can adjust Q in real time, and verified its effectiveness in a DC sensorless control system19. Compared with other improved methods, the main advantages of this method are its simple structure and small computational effort20. However, estimating the noise covariance still requires traversing new intervals at each historical moment, which is not practical for the implementation of the algorithm21. Literature22 applied the SageHusa Kalman filtering algorithm to position sensorless control systems for PMSMs and permanent magnet synchronous linear motors, both of which use the SageHusa noise estimator to achieve a recursive estimation of R, which is then further solved to find Q. However, the conventional SageHusa noise estimator leads to the failure of observer convergence.
The aim of this paper is to address the limitations of traditional methods by improving the Sage-Husa noise estimator integrated within the EKF framework of a PMSM sensorless control system. While previous implementations utilize the Sage-Husa estimator for recursive estimation of the noise covariance matrices Q and R, they often encounter the problem of non-convergence. Our proposed improvement retains the advantages of fast convergence and adaptability, while significantly improving the stability and estimation accuracy of the observer. The next sections describe our approach in detail, show experimental results, and discuss the implications of our findings for future PMSM applications.
PMSM mathematical model
In the d-q axis reference coordinate system, if the magnetic saturation effect is not considered, the extended reverse potential Eq. (1) of the PMSM is shown below:
where \(\:{\text{u}}_{\text{d}}\) and \(\:{\text{u}}_{\text{q}}\) represent the d-q axis voltage, while \(\:{\text{i}}_{\text{d}}\) and \(\:{\text{i}}_{\text{q}}\) represent the d-q axis current. Rs denotes the stator winding resistance; \(\:{\text{L}}_{\text{d}}\) and \(\:{\text{L}}_{\text{q}}\) stand for the d-q axis inductance. \(\:{{\upomega\:}}_{\text{e}}\:\)signifies the rotational speed of the magnetic flux, and Eq. (2)\(\:\:{\text{E}}_{\text{a}}\:\)denotes the amplitude of the extended counter electromotive force of the PMSM expressed as:
By performing the inverse Park transformation on Eq. (1), the extended inverse electromotive force model in the \(\:{\upalpha\:}-{\upbeta\:}\) axis coordinate system can be obtained:
The Eq. (4) of motion for PMSM is:
Assumptions
In the theoretical calculation process, motor damping is typically neglected. It is assumed that the external load torque changes slowly and can be compensated for by the PI controller in the integral link. Consequently, the motor load is disregarded to simplify the model, as demonstrated below:
where\(\:{\:\text{T}}_{\text{e}}\:\)is the motor torque,\(\:{\:\text{p}}_{\text{n}}\)is the number of motor pole pairs,\(\:{\:\text{J}}_{\text{m}}\) is the motor inertia momentum, and ψf is the motor magnetic chain.
Observer design
In this paper, Kalman Filtering is employed to estimate EMF of PMSM with high accuracy and robustness. By choosing\(\:{\:\text{i}}_{{\upalpha\:}}\)and\(\:{\:\text{i}}_{{\upbeta\:}}\)as state variables and incorporating them as extended state variables, the full-order state equations of the PMSM are formulated as follows:
In Eq. (6),\(\:{\text{L}}_{\text{d}}^{{\prime\:}}\)presents the dynamic inductance, which maintains the same value as the steady-state inductance when magnetic saturation of the magnetic chain is not taken into account. \(\:{e}_{\alpha\:}\) and \(\:{e}_{\beta\:}\) indicate the back electromotive force of the \(\:{\upalpha\:}-{\upbeta\:}\) axis.
Improved sage-huse for full-order EKF observer
The state space equation of the PMSM nonlinear system can be expressed as:
The input state variables, control variables, and output state variables needed to extend the Kalman observer are selected as:\(\:\text{x}\left(\text{t}\right)={[{\text{i}}_{{\upalpha\:}},{\text{i}}_{{\upbeta\:}},{\text{e}}_{{\upalpha\:}},{\text{e}}_{{\upbeta\:}}]}^{\text{T}}\), \(\:\text{u}\left(\text{t}\right)={[{\text{u}}_{{\upalpha\:}},{\text{u}}_{{\upbeta\:}}]}^{\text{T}}\), \(\:\text{Z}\left(\text{t}\right)={[{\text{i}}_{{\upalpha\:}},{\text{i}}_{{\upbeta\:}}]}^{\text{T}}\).\(\:{\upsigma\:}\left(\text{t}\right)\) is the system noise and \(\:{\upmu\:}\left(\text{t}\right)\) is the measurement noise.
where \({t_k}=k{T_s},\,\,{t_{k - 1}}=(k - 1){T_s}\), Ts is the sampling time, and linearization is performed to obtain the Jacobi matrix F and H, and the state transfer matrix\(\:{\Phi\:}\), respectively:
After formulating a discrete model of the non-linear PMSM system, the Sage-Husa adaptive EKF algorithm is utilized for estimating the state variables using the measured input and output data. The Sage-Husa algorithm demonstrates improved adaptability to the system’s non-linear and non-Gaussian characteristics compared to conventional Kalman Filtering methods, thus enhancing the rate of noise estimation convergence while accommodating the adaptive capabilities of the dynamic system.
Step 1: Forecasting.
In the following Eq. (15),\(\:{\stackrel{\sim}{x}}_{k}\) represents the a priori state value, while\(\:{\:\widehat{q}}_{k}\:\)signifies the system noise’s mean value computed using the Sage-Huse noise estimator.
In the following Eq. (16), \(\:{\stackrel{\sim}{\text{P}}}_{\text{k}}\) represents the a priori error covariance matrix, while\(\:{\widehat{Q}}_{k}\) denotes the system noise matrix at the kth iteration.
The difference between the best predicted value and the measured value is called new interest:
Step 2: Update.
In the following Eq. (18),\(\:{\widehat{R}}_{k}\) represents the measurement noise matrix at the kth iteration, and the Kalman Filter gain is:
The updated optimal estimate is:
The updated optimal covariance matrix is:
The Sage-Husa noise estimator is a commonly used noise mean and covariance estimation method, which is combined with the KF to form the Sage-Husa Kalman Filter. The traditional Sage-Husa noise estimator simultaneously estimates the mean of system noise\(\:{\:\widehat{\text{q}}}_{\text{k}}\:\)and measurement noise\(\:{\:\widehat{\text{r}}}_{\text{k}}\)in the Kalman Filter.The recursive formula for these estimates involves subtraction operations to ensure unbiased covariance estimation.However, this algorithm’s straightforward implementation can compromise the semi-positive definiteness of the noise covariance matrix, resulting in observer dispersion. To address this issue, a stabilization coefficient\(\:{\lambda\:}_{s}({\lambda\:}_{s}>0)\)is introduced to enhance observer stability. Its respective recursive expression is:
The variables\(\:{d}_{k}\)and \(\:{b}_{k}\)and represent the weighting coefficients in the noise covariance matrix, while\(\:\gamma\:\)represents the reserve coefficient for the noise mean. The coefficients\(\:{d}_{k}\)are determined based on the values provided in Eq. 25, primarily influenced by the magnitude of the new interest rate.
In order to prevent observer dispersion, the covariance matching principle is used to determine the:
When\(\gamma =1\), constitutes the most stringent convergence criterion. If Eq. (24) holds it means that the observer has diverged, at which point the system noise covariance\({\widehat {Q}_k}\)computed in Eq. (22) is not usable and the estimate\({\widehat {Q}_{k - 1}}\)from the previous moment can be chosen instead.
NQPLL
The extended inverse potential amplitude is essential in traditional PLL parameter design, leading to the utilization of a NQPLL for extracting rotor speed and position information. The normalized extended inverse electromotive force is subsequently integrated into the QPLL. Block diagram of NQPLL structure is as shown in Fig. 1.
The normalized error signal can be expressed as:
When \({\theta _e} - {\widehat {\theta }_e}<\pi /2\), the normalized error signal can be simplified as:
The simplified block diagram of the structure represents the transfer function as:
In the block diagram of the NQPLL structure,\(\:{k}_{pp\:}\)and\(\:{k}_{pi\:}\)represent the proportional and integral gains of the PI regulator. The PLL parameters are fine-tuned using the pole configuration method. Here,\(\:{\omega\:}_{p}\:\)specifies the bandwidth of the NQPLL, with a typical selection of\(\:{\:k}_{pp}=2{\omega\:}_{p}\) and \(\:{k}_{pi}={\omega\:}_{p}^{2}\).
LPF
Due to the nonlinearity of the inverter with magnetic field space harmonics, which lead to the presence of high-frequency harmonics in the extended reaction potential, the estimated position may incorporate sub-harmonics. To mitigate this concern, a low-pass filter is utilized to dampen resonance in the signal. The transfer function is expressed as follows:
where\({k_F}\)is the coefficient in the low-pass filter.
Speed loop controller design
SMC
The error between the rotor speed estimated by the observer and the reference speed to the speed loop controller, whose controller uses a slip film controller with greater immunity and robustness, and whose slip film controller is designed as follows:
Define the state variables of the SMC controller of the PMSM system as:
where \(\:{x}_{1}\:\)represents the rotor speed error,\(\:{\omega\:}_{ref}\)denotes the set reference speed, and \(\:{\omega\:}_{m}\)signifies the estimated rotor speed. In theoretical calculations, motor damping is often disregarded to derive the equation of state for the velocity error:
By making the above equation\({\dot {i}_q}=u\), the equation of state can be made as follows:
where u is the SMC control law, \({K_t}=3/2 \cdot {p_n}{\psi _f}\), J is the inertial momentum;
Define the synovial surface as:
In this context, \(\:{x}_{2}\) represents the first-order derivative of the rotational speed error, with c being the adjustable parameter. The following equation is derived:
The convergence rate of the slip flux controller is designed as \(\:\dot{s}=(-{k}_{sp}\bullet\:\sqrt{\left|s\right|}-{k}_{si}/S),{k}_{sp}、{k}_{si}>0\). The integral term S, denoted by \(\:{k}_{si}/S\), indicates integration, while all other terms are associated with sliding mode surfaces.
The SMC control law expression is:
MTPA
Utilizing the PMSM electromagnetic torque formula, the Lagrange pole method is applied to determine an optimal set of combinations\(\:{i}_{d}\)and\(\:{i}_{q}\), enabling the implementation of maximum torque-current ratio control.
The maximum current formula is as follows:
Considering the electromagnetic torque equation, the partial derivative with respect to the independent variable \(\:{i}_{q}\) is:
The final derivation\({i_d},{i_q}\)of equation:
Simulation results
In Fig. 2, the speed loop contains two modules: SMC and MTPA. It takes the observed rotor speed estimation error \(\:{\varDelta\:\omega\:}_{m}\) as input, calculated the electromagnetic torque\(\:\:{T}_{e}\:\)through the slip film controller designed in this paper, and then realize the maximum torque current ratio control through MTPA to find the optimum\(\:{i}_{d}\)and\(\:{i}_{q}\). Initially, the sampling of motor parameters \(\:{i}_{a}\), \(\:{i}_{b}\), and \(\:{i}_{c}\)is conducted. Following the Clark transformation and Park transformation to get a comparison, the PI regulator in the current loop outputs axis voltage. Subsequently, after inverse Park transformation into the voltage of the axis, the SVPWM algorithm modulation, written using a MATLAB function, generates the coded values for the state of the three half-bridges at that particular moment in time. These coded values are then utilized to control the switching state of the IGBT inverter, thereby driving the motor. The observer contains three modules: ISH-FEKF, NQPLL & LPF, whose derivation equations are shown above. All of them are simulation models built by math module.
The parameters of the two PMSM control system are shown in Table 1, and since the maximum speed of PMSM1 is 1600 rpm, PMSM2 is used to test the estimation accuracy of the observer at 2000 rpm and 3000 rpm. The PI parameters of the speed and current loops are parameterized using a typical type II system, which suppresses the overshoot and unstable fluctuations in the speed waveform so that the current system’s can perform well.
The convergence and stability of the conventional EKF observer have been rigorously evaluated. The control system parameters are simulated and tested using a parameter identification algorithm to establish the initial values of the process noise matrix Q and the measurement noise matrix R. These initial values are set in a manner that ensures optimal performance across different speed bands. The noise covariance matrix (QR) of the conventional EKF observer for PMSM1 and PMSM2 is presented in Eqs. (41) and (42).
The comparison of observed speeds within each speed domain of 500 rpm, 1000 rpm, and 1500 rpm is examined, as illustrated in Fig. 3. The total sampling time of the motor system is set to 10 s, with a sampling period of 1e-4. fpwm is 10 kHz. In Fig. 3, the black waveform represents the actual RPM of the motor controlled by the Hall sensor, the red waveform illustrates the observed RPM from the conventional EKF, the green waveform depicts the observed RPM from the EKF incorporating Sage-Husa, and the red waveform exhibits the observed RPM from the full-order EKF enhanced with the proposed Sage-Husa modifications in this study.
At 500 rpm, the motor reaches steady-state within 0.4 s after startup and acceleration, with minimal overshooting. The traditional EKF displays evident error fluctuations of approximately 3 rpm, while both SH-EKF and ISH-FEKF exhibit negligible observation errors. However, the SH-EKF observation waveforms demonstrate noticeable distortion, whereas the ISH-FEKF observation waveforms exhibit slight shifts.
As the motor speed increases to 1000 rpm, the traditional EKF shows an observation error of around 5 rpm. The performance of the traditional Sage-Husa variant deteriorates, resulting in increased waveform distortion and observation errors ranging from 3 to 5 rpm. The ISH-FEKF demonstrates optimal observation performance stability in this scenario. Upon reaching 1500 rpm, the observation error of the traditional EKF observer escalates to about 10 rpm. The SH-EKF observer struggles to maintain stability in the high-speed domain, whereas the ISH-FEKF observer showcases robust observation performance and stability at 1500 rpm.
Upon reaching 1500 rpm, the observation error of the traditional EKF observer escalates to about 10 rpm. The SH-EKF observer struggles to maintain stability in the high-speed domain, whereas the ISH-FEKF observer showcases robust observation performance and stability at 1500 rpm.
Figure 4 illustrates the comparison of motor torque among different observers within the PMSM position sensor-less system, wherein the torque reflects the error in the observed speed. At a motor speed of 500 rpm, the motor torque values for the traditional EKF observer range from 0.145 to 0.15 rad, for the SH-EKF observer range from 0.14 to 0.15 rad, and for the ISH-FEKF observer range from 0.128 to 0.132 rad. The ISH-FEKF observer demonstrates improved stability and reduced motor torque by approximately 0.015 rad compared to the SH-EKF observer.
As the motor speed increases to 1000 rpm, the motor torque controlled by the traditional EKF observer ranges from 0.2 to 0.218 rad, by the SH-EKF ranges from 0.182 to 0.212 rad, and by the ISH-FEKF ranges from 0.182 to 0.195 rad. Compared to the traditional EKF observer, the ISH-FEKF observer reduces the motor torque by around 0.023 rad, and compared to the SH-EKF observer, the reduction is about 0.015 rad, with significantly decreased fluctuation range observed in the SH-EKF.
Upon reaching a motor speed of 1500 rpm, the motor torque controlled by the traditional EKF observer ranges from 0.228 to 0.252 rad, the SH-EKF observer exhibits instability, and the ISH-FEKF observer ranges from 0.202 to 0.22 rad. The ISH-FEKF observer achieves a torque reduction of 0.026 rad compared to the traditional EKF observer.
Figure 5 illustrates the ISH-FEKF observation of the rotor position error within each speed domain during sudden speed changes. The estimation accuracy of the rotor position can be maintained at 0.0001 rad across varying speeds. A specific setting incorporated in the ISH-FEKF observer is 0.01, enhancing the observer’s convergence performance. This setting ensures improved performance in achieving observer convergence during sudden speed changes.
Figure 6a displays the observed speeds of the traditional EKF at 2000 rpm and 3000 rpm. In the graph, the blue curve represents the actual motor speed, while the red curve depicts the observed speed. At a reference speed of 2000 rpm, the EKF observation curve fluctuates between 1958 and 1978 rpm, resulting in an error of approximately 20–40 rpm compared to the actual motor speed. When the reference speed is set to 3000 rpm, the EKF observation curve fluctuates within the range of 2940–2990 rpm, with an error ranging from 5 to 55 rpm.
Figure 6b illustrates the observation speeds of the ISH-FEKF using the filtered curve. By employing filtering technology, the FEKF observation speed exhibits reduced fluctuation in the observation speed curve. At 2000 rpm, the error with the actual motor speed is approximately 24 rpm, while at 3000 rpm, this error decreases to 0-18 rpm, leading to a reduction of about 30 rpm compared to the traditional EKF observer.
Figure 7 presents a comparative analysis of motor torque under sudden load fluctuations, with the implementation of a sliding mode controller to address disruptive factors such as load variations and enhance response speed. In the context of simulation testing, a specific load is applied to the motor with a reference speed of 1000 rpm. The conventional PI controller results in motor speed lag, where adjusting the proportional term to suppress overshooting concurrently has a significant impact on response speed, leading to a noticeable decline in speed post-load addition.The utilization of the sliding mode controller effectively alleviates the aforementioned issues, with poor parameter settings potentially causing sliding mode chattering.
The value of c in the SMC module is 26,\(\:{k}_{si}\) is 3.2, and \(\:{k}_{sp}\:\)is 0.225. The duration for the rotor to transition from start-up acceleration to a steady state is reduced from 0.4 s to 0.25 s.
Figure 8 displays the results of the motor forward and reverse tests conducted at 200 rpm and 500 rpm. The speed tracking performance is excellent in the low-speed range at both 200 rpm and 500 rpm, as evidenced by the electrical angle changes depicted in Fig. 9a,b respectively. Notably, the response time for speed increase and decrease at 500 rpm is under 0.04 s, significantly quicker than that observed at 200 rpm. However, it is observed that the observer exhibits fluctuations during forward and reverse movements at medium and high rotational speeds, with a higher likelihood of fluctuations occurring at higher rotational speeds.
Physical results
The experimental platform of this study, as illustrated in Fig. 10, primarily comprises three components: the motor body, the drive control circuit, and the host computer. In this research, an adaptive extended Kalman Filter observer is employed in the sensorless control system for PMSM. Specific configuration operations carried out on the main MCU STM32G4 include enabling the Uart serial port to output waveforms of three-phase currents and motor speed for display on the VOFA host computer. Three-phase current sampling is performed using a combination of op-amps and ADC. TIM1, ADC, COMP analog comparator, and DAC are utilized to generate three-phase complementary PWM waves with dead-time and basic overcurrent protection. Additionally, GPIO controls LED and captures input key signals. The DAC converts digital signals within the MCU to analog outputs for debugging in an oscilloscope. The rated input voltage of the drive control circuit board is 24 V. It is connected to the motor’s three-phase lines and sequentially linked to ST-Link’s TX, RX, DIO, CLK for communication with the PC. Finally, the sensorless algorithm developed in MATLAB Simulink is used to generate embedded code, which is then burned into the drive control circuit board for experimental testing.
The tracking performance of the conventional EKF and ISH-FEKF observers at 500 rpm, 1000 rpm and 1500 rpm is tested under the PMSM1 control system. Figure 11a shows the real-time RPM curve observed by the EKF when the reference speed is set to 500 rpm, and it can be seen that the real-time RPM curve is controlled in the range of 480–518 rpm for 4 s, with about rpm error. Figure 11b shows the real-time speed curve observed by ISH-FEKF when the reference speed is set to 500 rpm, and it can be seen that the real-time speed curve is controlled to float in the range of 482-517 rpm, and there is about rpm error. The experimental results show that at 500 rpm, the observation accuracy of the two observers is similar.
Figure 12a shows the real-time speed curve observed by EKF when the reference speed is set to 1000 rpm, and it can be seen that the real-time speed curve in 2s is controlled to float in the range of 982–1018 rpm, and there is about rpm error. Figure 12b shows the real-time speed curve observed by ISH-FEKF when the reference speed is set to 1000 rpm, and it can be seen that the real-time speed curve is controlled to float in the range of 985–1012 rpm. The experimental results show that at 1000 rpm, the error interval of ISH-FEKF observation is narrowed down by 3–6 rpm compared with that of the conventional EKF observer, and as the rotational speed increases, the fluctuation of the observed rotational speed curves increases, and the real-time rotational speed curves observed by ISH-FEKF are obviously larger than that of the conventional EKF.
Figure 13a depicts the real-time speed curve observed by the EKF when the reference speed is set to 1500 rpm. It is evident that the real-time speed curve fluctuates within the 1478–1522 rpm range within 1 s, with an approximate rpm error. In Fig. 13b, the real-time speed curve observed by the ISH-FEKF at a reference speed of 1500 rpm shows fluctuations within the range of 1488–1512 rpm, with a similar approximate rpm error. The experimental findings demonstrate that at 1500 rpm, the error interval observed by the ISH-FEKF observer is reduced by 10 rpm compared to the traditional EKF observer. While the accuracy of the ISH-FEKF observer improves with higher speeds, the inevitable level of fluctuations also escalates due to the adaptive parameters. The test results closely align with the simulation outcomes.
It is tested whether the observer designed in this paper can converge quickly to ensure the normal operation of the motor at 3000 rpm. As shown in Fig. 14 (a) (b), the experimental results show that at 2000 rpm, the error interval of the traditional EKF observation is about 20–40 rpm, and the error interval of the ISH-FEKF observation is about 10–15 rpm, which is narrowed down by 10–25 rpm compared with the conventional EKF observation, and the optimization of the LPF parameter should alleviate the problem, which is probably due to the adaptive parameter leading to the observation of the waveform appearing high-frequency spike signals several times. Optimization of the LPF parameters should alleviate the problem.
As shown in Fig. 15a, b, the experimental results show that at 3000 rpm, the error interval of traditional EKF observation is about 55–70 rpm, and the error interval of ISH-FEKF observer is about rpm under the ideal condition of stable signal without sharp pulse signal, which is narrowed down by rpm in comparison. In Figs. 14 and 15, the red waveform is the stator current and the green waveform is the electrical angle, and the electrical angle observed by ISH-FEKF is smoother in comparison.
Conclusion
In this research, we introduce an improved Sage-Huse adaptive Kalman observer for sensorless position control of PMSM. The study involves testing the traditional EKF and the adaptive EKF with the addition of the Sage-Huse noise estimator using MATLAB Simulink. The estimation performance of the conventional EKF observer in the system is evaluated, highlighting the issue of increased dispersion after incorporating Sage-Huse. Subsequently, we identify the dispersion problem and enhance the observer accordingly. The sensorless control algorithm is developed by integrating SMC and MTPA, followed by simulation tests under various conditions such as fixed load and variable speed, fixed speed and variable load, as well as forward and reverse rotation. The effectiveness of the ISH-FEKF observer is demonstrated in experiments conducted on the PMSM drive control platform utilizing the STM32G4 master MCU. Through experimental tests ranging from 500 to 3000 rpm, the results demonstrate that the ISH-FEKF observer can sustain a high level of estimation accuracy despite fluctuations in reference speed and the system environment. However, it may also exhibit potential instability, resulting in heightened fluctuations. At 3000 rpm, while ensuring the motor’s normal operation, the ISH-FEKF observer demonstrates a narrower error margin in rpm compared to the traditional EKF observations. The ISH-FEKF observer maintains excellent estimation accuracy across various speed ranges and control systems by incorporating new information to update the noise covariance between Q and R based on the error between predicted and measured values. The observer showcases superior performance in terms of accuracy, adaptability, and rapid convergence compared to the conventional EKF, which struggles with varying speed conditions due to its reliance on QR integration based on the current system state. The proposed ISH-FEKF observer stands out for its high accuracy, adaptability, and efficient convergence properties.
Data availability
All data generated or analyzed during this study are included in this published article.
References
Pillay, P. & Krishnan, R. Application Characteristisc of Permanent Magnet Synchronous and Brushless DC Motor for Servo drives. IEEE Trans. Industrial Appl. 27 (5), 986–996 (1991).
Perera, P. et al. A sensorless, stable v/f control method for permanentmagnet synchronous motor drives[J]. IEEE Trans. Ind. Appl. 39 (3), 783–791 (2003).
Xiufeng Liu, Y. et al. Composite control based on FNTSMC and adaptive neural network for PMSM system,ISA Transactions, (2024).
Pasqualotto, D., Rigon, S. & Zigliotto, M. Sensorless speed control of synchronous reluctance motor drives based on extended kalman filter and neural magnetic model[J]. IEEE Trans. Industr. Electron. 70 (2), 1321–1330 (2023).
Jiang, F. et al. Robustness improvement of modelbased sensorless spmsm drivers based on an adaptive extended state observer and an enhanced quadrature pll[J]. IEEE Trans. Power Electron. 36 (4), 4802–4814 (2021).
Zhang, Y. et al. A rotor position and speed estimation method using an improved linear extended state observer for ipmsm sensorless drives[J]. IEEE Trans. Power Electron. 36 (12), 14062–14073 (2021).
Zhang, G. et al. Adalinenetworkbased pll for position sensorless interior permanent magnet synchronous motor drives[J]. IEEE Trans. Power Electron. 31 (2), 1450–1460 (2016).
Kim, H., Son, J. & Lee, J. A highspeed slidingmode observer for the sensorless speed control of a pmsm[J]. IEEE Trans. Industr. Electron. 58 (9), 4069–4077 (2011).
Wu, C., Sun, X. & Wang, J. A rotor flux observer of permanent magnet synchronous motors with adaptive flux compensation[J]. IEEE Trans. Energy Convers. 34 (4), 2106–2117 (2019).
Xu, W. et al. Improved nonlinear flux observerbased secondorder soifo for pmsm sensorless control[J]. IEEE Trans. Power Electron. 34 (1), 565–579 (2019).
Bolognani, S., Oboe, R. & Zigliotto, M. Sensorless fulldigital pmsm drive with ekf estimation of speed and rotor position[J]. IEEE Trans. Industr. Electron. 46 (1), 184–191 (1999).
Barut, M. et al. Realtime implementation of Bi inputextended kalman filterbased estimator for speedsensorless control of induction motors[J]. IEEE Trans. Industr. Electron. 59 (11), 4197–4206 (2012).
Kim, Y. H. & Kook, Y. S. High performance ipmsm drives without rotational position sensors using reduced order ekf[J]. IEEE Trans. Energy Convers. 14 (4), 868–873 (1999).
Bolognani, S., Tubiana, L. & Zigliotto, M. Extended kalman filter tuning in sensorless pmsm drives[J]. IEEE Trans. Ind. Appl. 39 (6), 1741–1747 (2003).
Salvatore, N. et al. Optimization of delayedstate kalmanfilterbased algorithm via differential evolution for sensorless control of induction motors[J]. IEEE Trans. Industr. Electron. 57 (1), 385–394 (2010).
Zerdali, E. & Barut, M. The comparisons of optimized extended kalman filters for speedsensorless control of induction motors[J]. IEEE Trans. Industr. Electron. 64 (6), 4340–4351 (2017).
Yin, Z. et al. A speed and flux estimation method of induction motor using fuzzy extended kalman filter[C]. In 2014 International Power Electronics and Application Conference and Exposition, pp. 693698. (2014).
Zerdali, E. Adaptive extended kalman filter for speedsensorless control of induction motors[J]. IEEE Trans. Energy Convers. 34 (2), 789–800 (2019).
Jiancheng, F. & Sheng, Y. Study on innovation adaptive ekf for inflight alignment of airborne pos[J]. IEEE Trans. Instrum. Meas. 60 (4), 1378–1388 (2011).
Larik, N. A., Li, M. S. & Wu, Q. H. Enhanced Fault detection and localization strategy for high-speed protection in medium-voltage DC distribution networks using extended Kalman Filtering Algorithm. IEEE Access 12, 30329–30344. https://doi.org/10.1109/ACCESS (2024).
Larik, N., Ali, M., Li & Wu, Q. Triple-indexed Passive Islanding Detection Strategy for grid-connected Distributed Generation Networks Using an Extended Kalman Filter ( IET Generation, Transmission & Distribution, 2024).
Wang, T. et al. Adaptive extended kalman filter based dynamic equivalent method of pmsg wind farm cluster[J]. IEEE Trans. Ind. Appl. 57 (3), 2908–2917 (2021).
Acknowledgements
This work was supported in part by 2024 Ministry of Education Humanities and Social Sciences Research Planning Fund Project (24YJAZH116).
Author information
Authors and Affiliations
Contributions
Sang Yingjun: Part of the writing and editing, Wu Zhenglong: Main writing, Fan Yuanyuan: Data collection and editing.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
About this article
Cite this article
Yingjun, S., Zhenglong, W. & Yuanyuan, F. An adaptive extended Kalman filter observer-for permanent magnet synchronous motor position sensorless control systems. Sci Rep 15, 11605 (2025). https://doi.org/10.1038/s41598-025-85787-5
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-025-85787-5