<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing with OASIS Tables v3.0 20080202//EN" "https://jats.nlm.nih.gov/nlm-dtd/publishing/3.0/journalpub-oasis3.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:oasis="http://docs.oasis-open.org/ns/oasis-exchange/table" xml:lang="en" dtd-version="3.0" article-type="research-article">
  <front>
    <journal-meta><journal-id journal-id-type="publisher">MS</journal-id><journal-title-group>
    <journal-title>Mechanical Sciences</journal-title>
    <abbrev-journal-title abbrev-type="publisher">MS</abbrev-journal-title><abbrev-journal-title abbrev-type="nlm-ta">Mech. Sci.</abbrev-journal-title>
  </journal-title-group><issn pub-type="epub">2191-916X</issn><publisher>
    <publisher-name>Copernicus Publications</publisher-name>
    <publisher-loc>Göttingen, Germany</publisher-loc>
  </publisher></journal-meta>
    <article-meta>
      <article-id pub-id-type="doi">10.5194/ms-15-99-2024</article-id><title-group><article-title>Improved flux linkage observer for position estimation of permanent magnet synchronous linear motor</article-title><alt-title>Improved flux linkage observer for position estimation of PMSLMs​​​​​​​</alt-title>
      </title-group><?xmltex \runningtitle{Improved flux linkage observer for position estimation of PMSLMs​​​​​​​}?><?xmltex \runningauthor{W. Yu et al.}?>
      <contrib-group>
        <contrib contrib-type="author" corresp="no" rid="aff1">
          <name><surname>Yu</surname><given-names>Wenbin</given-names></name>
          
        </contrib>
        <contrib contrib-type="author" corresp="yes" rid="aff1">
          <name><surname>Yang</surname><given-names>Guolai</given-names></name>
          <email>yangglnjust@126.com</email>
        </contrib>
        <contrib contrib-type="author" corresp="no" rid="aff1">
          <name><surname>Wang</surname><given-names>Liqun</given-names></name>
          
        </contrib>
        <contrib contrib-type="author" corresp="no" rid="aff1">
          <name><surname>Lin</surname><given-names>Darui</given-names></name>
          
        </contrib>
        <contrib contrib-type="author" corresp="no" rid="aff2">
          <name><surname>Al-Zahrani</surname><given-names>Ahmed</given-names></name>
          
        </contrib>
        <aff id="aff1"><label>1</label><institution>School of Mechanical Engineering, Nanjing University of Science &amp; Technology, Nanjing 210014, China</institution>
        </aff>
        <aff id="aff2"><label>2</label><institution>School of Mechanical and Materials Engineering, University of Jeddah, Jeddah 23881, Saudi Arabia​​​​​​​</institution>
        </aff>
      </contrib-group>
      <author-notes><corresp id="corr1">Guolai Yang (yangglnjust@126.com)</corresp></author-notes><pub-date><day>19</day><month>February</month><year>2024</year></pub-date>
      
      <volume>15</volume>
      <issue>1</issue>
      <fpage>99</fpage><lpage>109</lpage>
      <history>
        <date date-type="received"><day>16</day><month>April</month><year>2023</year></date>
           <date date-type="rev-recd"><day>15</day><month>December</month><year>2023</year></date>
           <date date-type="accepted"><day>25</day><month>December</month><year>2023</year></date>
      </history>
      <permissions>
        <copyright-statement>Copyright: © 2024 Wenbin Yu et al.</copyright-statement>
        <copyright-year>2024</copyright-year>
      <license license-type="open-access"><license-p>This work is licensed under the Creative Commons Attribution 4.0 International License. To view a copy of this licence, visit <ext-link ext-link-type="uri" xlink:href="https://creativecommons.org/licenses/by/4.0/">https://creativecommons.org/licenses/by/4.0/</ext-link></license-p></license></permissions><self-uri xlink:href="https://ms.copernicus.org/articles/15/99/2024/ms-15-99-2024.html">This article is available from https://ms.copernicus.org/articles/15/99/2024/ms-15-99-2024.html</self-uri><self-uri xlink:href="https://ms.copernicus.org/articles/15/99/2024/ms-15-99-2024.pdf">The full text article is available as a PDF file from https://ms.copernicus.org/articles/15/99/2024/ms-15-99-2024.pdf</self-uri>
      <abstract><title>Abstract</title>

      <p id="d1e120">With the development of motor technology, sensorless control attracts more and more attention. In this paper, an improved flux linkage observer is proposed to overcome issues including inaccurate initial positions and sampling noise. The voltage and current models are combined, and a sliding-mode observer is designed based on the hybrid model to obtain the compensation voltage. Therefore, the estimated flux linkage after compensation can greatly resist the influence caused by inaccurate initial positions or sampling noise. Phase-locked loop technology is used to process the estimated flux linkage to get the stable estimated speed and position. The proposed scheme has a simple structure and only one parameter. It is easy to use and adjust in practice. The simulation and experimental results verify that the proposed algorithm is effective, and the estimated flux linkage and position is accurate with an inaccurate initial position.</p>
  </abstract>
    
<funding-group>
<award-group id="gs1">
<funding-source>National Postdoctoral Program for Innovative Talents</funding-source>
<award-id>BX2021126</award-id>
</award-group>
<award-group id="gs2">
<funding-source>National Natural Science Foundation of China</funding-source>
<award-id>52105106</award-id>
</award-group>
<award-group id="gs3">
<funding-source>China Postdoctoral Science Foundation</funding-source>
<award-id>2021M701711</award-id>
</award-group>
</funding-group>
</article-meta>
  </front>
<body>
      

<sec id="Ch1.S1" sec-type="intro">
  <label>1</label><title>Introduction</title>
      <p id="d1e132">Permanent magnet synchronous linear motors (PMSLMs) are widely used in various fields such as precision machining, oilfield equipment, rail transportation, etc. due to their special motion, high thrust density and high motion accuracy (Cao et al., 2015; Pellegrino et al., 2012). Because of the unique structure, a long-enough position detection system is always needed in PMSLMs for a good control performance (Lin et al., 2021). On the other hand, linear position encoders are susceptible to damage and need regular maintenance (Cao et al., 2019; Lu et al., 2011). Especially in harsh environments, external impacts or collisions very easily damage the linear position encoder and make the system unable to work properly. These all greatly limit the application and development of PMSLMs. Therefore, a sensorless control algorithm is of great significance to the PMSLM.</p>
      <p id="d1e135">With the development, many sensorless control algorithms have been proposed to estimate the speed and position of rotors accurately. Some scholars achieved real-time estimation by injecting high-frequency pulse signals at zero or low speed (Wang et al., 2019; Ofori et al., 2015). This method makes use of the salient effect of the motor and analyzes the current response, which carries position information. However, this method is powerless for the motor whose salient effect is not obvious, and the high-frequency signal also increases the loss. Model reference adaptation (Cirrincione et al., 2013; Cheng et al., 2015), Kalman filtering (Pasqualotto et al., 2023; Lu et al., 2011), state observer (Accetta et al., 2014; Zhang et al., 2020) and other methods based on intelligent-control theory establish the rotor state observer according to the motor mathematical model and adjust the model in real time according to the operating state to accurately estimate the speed and position. These methods have strong robustness and anti-interference ability and are suitable for observing the motor state under long-term operation. On the other hand, the algorithms are also complex and require high computing capacity. The estimation effect at low speeds is also not ideal. The position estimation based on back electromotive force (EMF) has the advantages of simple structure, less calculation and fast response, but it also has the problems of low detection accuracy and being easily affected by a change in motor parameters (Hussain and Toliyat, 2015; Zhang et al., 2021; Wang and Ge, 2010). Especially at low speeds, the back EMF is too small, which makes the speed and position estimations more sensitive to noise. The<?pagebreak page100?> position estimation based on flux linkage is simple in structure and is easy to implement. The integration process greatly reduces the influence of current noise, but it also increases the influence of direct current (DC) bias. The open-loop detection method also makes the algorithm less anti-interference oriented (Zhao et al., 2019; Xiao et al., 2021). Each method has its advantages and disadvantages. To establish a suitable sensorless position estimation algorithm and to increase the reliability and anti-interference ability on the basis of simple structure and less computation is also one of the purposes of this paper.</p>
      <p id="d1e138">It is a key step to obtain the stator flux linkage by integrating the back EMF of the motor. Theoretically, the flux linkage obtained by integrating the back EMF can still work properly even at low speeds. For the DC bias and high-frequency noise of the current sensor (Pavlitov et al., 2014; Du et al., 2020), some scholars design various types of high-pass or low-pass filters to eliminate the impact (Cupertino et al., 2010; Xiao et al., 2017). Although it brings amplitude attenuation and phase offset, it also achieves good practical results. In fact, the integral process for the flux linkage always requires an accurate initial value of the flux linkage, i.e., the exact initial position of the motor. The inaccurate initial value of integration will lead to errors in the flux linkage estimation and have a negligible effect on the speed and position estimation, which has not been considered in much of the literature (Wang et al., 2012; Yang et al., 2022; Ting et al., 2019). In practical applications, due to the external interference, the direct drive of PMSLMs is more likely to lead to inaccurate initial positions. When the inaccurate initial position is used to calculate the flux linkage, error occurs, which will affect the speed and position estimation.</p>
      <p id="d1e141">To solve this problem and establish a position estimation scheme with good performance, an improved flux linkage observer is proposed in this paper. Sliding-mode control is widely known and used for its good robustness and anti-interference capacity (Afifa et al., 2023). The same sliding-mode observer is also famous because of its strong nonlinearity. Sliding-mode theory can always have a good performance in nonlinearity issues, and it can also be used in this paper. Different from the traditional open-loop calculation method, the flux linkage is calculated through the voltage model firstly and then uses the inaccurate flux linkage to estimate the current through the current model. Finally, a sliding-mode observer is established based on the estimated current and the actual current, and the output is applied to the voltage model as the compensation voltage to form a closed-loop structure. The output of the sliding-mode observer, as the compensation voltage, is used to correct the inaccurate flux linkage. Then, the inaccurate flux linkage will tend towards the actual value under the action of the sliding-mode observer. Compared to the other observers, the proposed observer has the advantages of a simple structure, only one parameter, reliable anti-interference and being easy to implement. The widely known and annoying chattering problem caused by a nonlinear sliding-mode function (Ahmad et al., 2023) can also be effectively suppressed because of the integral item in the modified flux linkage.</p>
      <p id="d1e145">Section 1 of this paper is an overview of position estimation and puts forward the existing problems. Section 2 describes the voltage, current model and common position calculation method of the motor. In Sect. 3, the improved flux linkage observer algorithm proposed in this paper is analyzed, and its stability is proved. Sections 4 and 5 simulate and experiment with the proposed flux linkage observer to further illustrate its performance. Section 6 summarizes the full paper.</p>
</sec>
<sec id="Ch1.S2">
  <label>2</label><title>Mathematical model</title>
      <p id="d1e156">The motor used in this paper is a surface-mounted PMSLM shown in Fig. 1. The windings are internal in the moving part, and permanent magnets are also internal in the fixed shaft. Injecting current into the windings, the moving part can move along the shaft. For good control performance, the position sensor is always needed and should cover the full length. As mentioned before, position information will be replaced by a sensorless scheme in this paper. To simplify the analysis, there are some assumptions that we have outlined below: <list list-type="custom"><list-item><label>a.</label>
      <p id="d1e161">all three phases of the motor have the same resistance and inductance;</p></list-item><list-item><label>b.</label>
      <p id="d1e165">the mutual inductance of the windings are ignored;</p></list-item><list-item><label>c.</label>
      <p id="d1e169">electromagnetic losses and eddy current effects are ignored.</p></list-item></list></p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F1"><?xmltex \currentcnt{1}?><?xmltex \def\figurename{Figure}?><label>Figure 1</label><caption><p id="d1e174">The permanent magnet synchronous tube linear motor used in this paper.</p></caption>
        <?xmltex \igopts{width=241.848425pt}?><graphic xlink:href="https://ms.copernicus.org/articles/15/99/2024/ms-15-99-2024-f01.png"/>

      </fig>

      <p id="d1e183">These simplifications are acceptable and reasonable for the topic in this paper because, under normal operating conditions, the influences of these factors are very small and can be ignored (Zhang et al., 2016). The composite flux linkage produced by permanent magnets and windings is called stator flux linkage, and that produced by permanent magnets is called rotor flux linkage. The two-phase voltage model in a static coordinate system can be expressed as
          <disp-formula id="Ch1.E1" content-type="numbered"><label>1</label><mml:math id="M1" display="block"><mml:mrow><mml:mfenced close="" open="{"><mml:mtable class="array" rowspacing="5.690551pt" columnalign="left"><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mi>R</mml:mi><mml:msub><mml:mi>i</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:mstyle displaystyle="false"><mml:mstyle displaystyle="false"><mml:mfrac style="text"><mml:mrow><mml:mi mathvariant="normal">d</mml:mi><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:mi mathvariant="normal">d</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac></mml:mstyle></mml:mstyle><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mi>R</mml:mi><mml:msub><mml:mi>i</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:mi mathvariant="normal">d</mml:mi><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:mi mathvariant="normal">d</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>.</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mfenced></mml:mrow></mml:math></disp-formula>
       <?pagebreak page101?> The stator current model in static coordinate system can be expressed as
          <disp-formula id="Ch1.E2" content-type="numbered"><label>2</label><mml:math id="M2" display="block"><mml:mrow><mml:mfenced open="{" close=""><mml:mtable rowspacing="5.690551pt" class="array" columnalign="left"><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>L</mml:mi><mml:mi mathvariant="normal">s</mml:mi></mml:msub><mml:msub><mml:mi>i</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="normal">r</mml:mi></mml:msub><mml:mi>cos⁡</mml:mi><mml:mi mathvariant="italic">θ</mml:mi><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>L</mml:mi><mml:mi mathvariant="normal">s</mml:mi></mml:msub><mml:msub><mml:mi>i</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="normal">r</mml:mi></mml:msub><mml:mi>sin⁡</mml:mi><mml:mi mathvariant="italic">θ</mml:mi><mml:mo>.</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mfenced></mml:mrow></mml:math></disp-formula>
        According to Eq. (1),
          <disp-formula id="Ch1.E3" content-type="numbered"><label>3</label><mml:math id="M3" display="block"><mml:mrow><mml:mfenced open="{" close=""><mml:mtable class="array" rowspacing="5.690551pt" columnalign="left"><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mo movablelimits="false">∫</mml:mo><mml:mfenced open="(" close=")"><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:mi>R</mml:mi><mml:msub><mml:mi>i</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:mfenced><mml:mi mathvariant="normal">d</mml:mi><mml:mi>t</mml:mi><mml:mo>+</mml:mo><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mrow><mml:mi mathvariant="italic">α</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:msub><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mo movablelimits="false">∫</mml:mo><mml:mfenced close=")" open="("><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:mi>R</mml:mi><mml:msub><mml:mi>i</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub></mml:mrow></mml:mfenced><mml:mi mathvariant="normal">d</mml:mi><mml:mi>t</mml:mi><mml:mo>+</mml:mo><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mrow><mml:mi mathvariant="italic">β</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:msub><mml:mo>.</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mfenced></mml:mrow></mml:math></disp-formula>
        Calculate the stator flux linkage from Eq. (3) and substitute the result into Eq. (2) to obtain the position angle:
          <disp-formula id="Ch1.E4" content-type="numbered"><label>4</label><mml:math id="M4" display="block"><mml:mrow><mml:mi mathvariant="italic">θ</mml:mi><mml:mo>=</mml:mo><mml:mi>arctan⁡</mml:mi><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mi>L</mml:mi><mml:mi mathvariant="normal">s</mml:mi></mml:msub><mml:msub><mml:mi>i</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mi>L</mml:mi><mml:mi mathvariant="normal">s</mml:mi></mml:msub><mml:msub><mml:mi>i</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula>
        The electrical angle can be obtained from Eq. (4), and the actual rotor angle is <inline-formula><mml:math id="M5" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">r</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mi mathvariant="italic">θ</mml:mi><mml:mo>/</mml:mo><mml:msub><mml:mi>p</mml:mi><mml:mi>n</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>. For linear motors, the actual movement position should be expressed as
          <disp-formula id="Ch1.E5" content-type="numbered"><label>5</label><mml:math id="M6" display="block"><mml:mrow><mml:mi>x</mml:mi><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mi mathvariant="italic">θ</mml:mi><mml:mrow><mml:mn mathvariant="normal">2</mml:mn><mml:mi mathvariant="italic">π</mml:mi><mml:msub><mml:mi>p</mml:mi><mml:mi>n</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mstyle><mml:mi mathvariant="italic">τ</mml:mi><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula>
        Actually, the noise in the current and voltage signal makes the estimated position and speed directly calculated by arctangent very rough. Therefore, the phase-locked loop technology is used in this paper instead of Eq. (4) to obtain a smooth estimated position and speed. The traditional position estimation method based on stator flux linkage and phase-locked loop (PLL) is shown in Fig. 2. To eliminate the influence of DC bias, current and voltage measurement noise, other forms of transfer functions are used to replace the pure integration (Pavlitov et al., 2014; Xiao et al., 2017), and this has achieved good results. Little attention has been paid to the problem of an inaccurate initial position of the motor in the literature.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F2"><?xmltex \currentcnt{2}?><?xmltex \def\figurename{Figure}?><label>Figure 2</label><caption><p id="d1e543">Position estimation based on phase-locked loop (PLL).</p></caption>
        <?xmltex \igopts{width=199.169291pt}?><graphic xlink:href="https://ms.copernicus.org/articles/15/99/2024/ms-15-99-2024-f02.png"/>

      </fig>

</sec>
<sec id="Ch1.S3">
  <label>3</label><title>Improved flux linkage observer</title>
<sec id="Ch1.S3.SS1">
  <label>3.1</label><title>Flux linkage error analysis and observer</title>
      <p id="d1e567">To solve the problem of the initial estimated position being inaccurate due to uncertain interference, the traditional method of stator flux linkage estimation is improved upon in this section. For voltage error caused by inverter nonlinearity, noise in the measuring current, inaccurate initial-phase angle and other interference factors, Eq. (3) can be rewritten as
            <disp-formula id="Ch1.E6" content-type="numbered"><label>6</label><mml:math id="M7" display="block"><mml:mrow><mml:mfenced close="" open="{"><mml:mtable rowspacing="0.2ex 5.690551pt 0.2ex" class="array" columnalign="left left"><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>=</mml:mo><mml:mo movablelimits="false">∫</mml:mo><mml:mfenced open="(" close=")"><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:mi>R</mml:mi><mml:msub><mml:mi>i</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:mfenced><mml:mi mathvariant="normal">d</mml:mi><mml:mi>t</mml:mi><mml:mo>+</mml:mo><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mrow><mml:mi mathvariant="italic">α</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mrow><mml:mi mathvariant="italic">α</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd/><mml:mtd><mml:mrow><mml:mo>-</mml:mo><mml:mo movablelimits="false">∫</mml:mo><mml:mfenced open="[" close="]"><mml:mrow><mml:mfenced close=")" open="("><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi>u</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:mi>R</mml:mi><mml:msub><mml:mi>i</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:mfenced><mml:mo>+</mml:mo><mml:mfenced open="(" close=")"><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:mi>R</mml:mi><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi>i</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:mfenced></mml:mrow></mml:mfenced><mml:mi mathvariant="normal">d</mml:mi><mml:mi>t</mml:mi><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>=</mml:mo><mml:mo movablelimits="false">∫</mml:mo><mml:mfenced close=")" open="("><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:mi>R</mml:mi><mml:msub><mml:mi>i</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub></mml:mrow></mml:mfenced><mml:mi mathvariant="normal">d</mml:mi><mml:mi>t</mml:mi><mml:mo>+</mml:mo><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mrow><mml:mi mathvariant="italic">β</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mrow><mml:mi mathvariant="italic">β</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd/><mml:mtd><mml:mrow><mml:mo>-</mml:mo><mml:mo movablelimits="false">∫</mml:mo><mml:mfenced close="]" open="["><mml:mrow><mml:mfenced open="(" close=")"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi>u</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:mi>R</mml:mi><mml:msub><mml:mi>i</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub></mml:mrow></mml:mfenced><mml:mo>+</mml:mo><mml:mfenced close=")" open="("><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:mi>R</mml:mi><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi>i</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub></mml:mrow></mml:mfenced></mml:mrow></mml:mfenced><mml:mi mathvariant="normal">d</mml:mi><mml:mi>t</mml:mi><mml:mo>.</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mfenced></mml:mrow></mml:math></disp-formula>
          Flux linkages <inline-formula><mml:math id="M8" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M9" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> are affected by errors of voltage, current and initial value of flux linkage. The integral initial-value error of the stator flux linkage of the <inline-formula><mml:math id="M10" display="inline"><mml:mrow><mml:mi mathvariant="italic">α</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:mrow></mml:math></inline-formula> axis may be caused by the inaccurate initial-position estimation, or it may be caused by the inconsistency between the actual value and the motor flux linkage mark. This results in a fixed offset of the estimated flux linkage value. To calculate the stator flux linkage more accurately, it is necessary to compensate for these errors.</p>
      <p id="d1e830">Although <inline-formula><mml:math id="M11" display="inline"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mrow><mml:mi mathvariant="italic">α</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M12" display="inline"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mrow><mml:mi mathvariant="italic">β</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> are fixed, it still can be expressed as a form of the integral of nonlinear time-varying functions <inline-formula><mml:math id="M13" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">δ</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M14" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">δ</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>.
            <disp-formula id="Ch1.E7" content-type="numbered"><label>7</label><mml:math id="M15" display="block"><mml:mrow><mml:mfenced open="{" close=""><mml:mtable rowspacing="5.690551pt" class="array" columnalign="left"><mml:mtr><mml:mtd><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mrow><mml:mi mathvariant="italic">α</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mo movablelimits="false">∫</mml:mo><mml:msub><mml:mi mathvariant="italic">δ</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mi mathvariant="normal">d</mml:mi><mml:mi>t</mml:mi><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mrow><mml:mi mathvariant="italic">β</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mo movablelimits="false">∫</mml:mo><mml:msub><mml:mi mathvariant="italic">δ</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub><mml:mi mathvariant="normal">d</mml:mi><mml:mi>t</mml:mi><mml:mo>.</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mfenced></mml:mrow></mml:math></disp-formula>
          Let
            <disp-formula id="Ch1.E8" content-type="numbered"><label>8</label><mml:math id="M16" display="block"><mml:mrow><mml:mfenced close="" open="{"><mml:mtable rowspacing="5.690551pt" class="array" columnalign="left"><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>e</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mfenced close=")" open="("><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi>u</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:mi>R</mml:mi><mml:msub><mml:mi>i</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:mfenced><mml:mo>+</mml:mo><mml:mfenced open="(" close=")"><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:mi>R</mml:mi><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi>i</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:mfenced><mml:mo>+</mml:mo><mml:msub><mml:mi mathvariant="italic">δ</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>e</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mfenced open="(" close=")"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi>u</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:mi>R</mml:mi><mml:msub><mml:mi>i</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub></mml:mrow></mml:mfenced><mml:mo>+</mml:mo><mml:mfenced open="(" close=")"><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:mi>R</mml:mi><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi>i</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub></mml:mrow></mml:mfenced><mml:mo>+</mml:mo><mml:msub><mml:mi mathvariant="italic">δ</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub><mml:mo>.</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mfenced></mml:mrow></mml:math></disp-formula>
          Thus, Eq. (6) can be further expressed as
            <disp-formula id="Ch1.E9" content-type="numbered"><label>9</label><mml:math id="M17" display="block"><mml:mrow><mml:mfenced open="{" close=""><mml:mtable rowspacing="5.690551pt" class="array" columnalign="left"><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mo movablelimits="false">∫</mml:mo><mml:mfenced close=")" open="("><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:mi>R</mml:mi><mml:msub><mml:mi>i</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mi>e</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:mfenced><mml:mi mathvariant="normal">d</mml:mi><mml:mi>t</mml:mi><mml:mo>+</mml:mo><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mrow><mml:mi mathvariant="italic">α</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:msub><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mo movablelimits="false">∫</mml:mo><mml:mfenced close=")" open="("><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:mi>R</mml:mi><mml:msub><mml:mi>i</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mi>e</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub></mml:mrow></mml:mfenced><mml:mi mathvariant="normal">d</mml:mi><mml:mi>t</mml:mi><mml:mo>+</mml:mo><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mrow><mml:mi mathvariant="italic">β</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:msub><mml:mo>.</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mfenced></mml:mrow></mml:math></disp-formula>
          For various possible errors, there is both <inline-formula><mml:math id="M18" display="inline"><mml:mrow><mml:msub><mml:mi>e</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M19" display="inline"><mml:mrow><mml:msub><mml:mi>e</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> in the flux linkage calculation, and they are always unknown. Actually, the estimated flux linkages <inline-formula><mml:math id="M20" display="inline"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi mathvariant="italic">λ</mml:mi><mml:mo stretchy="false" mathvariant="normal">^</mml:mo></mml:mover><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M21" display="inline"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi mathvariant="italic">λ</mml:mi><mml:mo stretchy="false" mathvariant="normal">^</mml:mo></mml:mover><mml:mi mathvariant="italic">β</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> are calculated by Eq. (9) without <inline-formula><mml:math id="M22" display="inline"><mml:mrow><mml:msub><mml:mi>e</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M23" display="inline"><mml:mrow><mml:msub><mml:mi>e</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> so that, if the compensation voltage <inline-formula><mml:math id="M24" display="inline"><mml:mrow><mml:msub><mml:mi>e</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M25" display="inline"><mml:mrow><mml:msub><mml:mi>e</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> can be observed, the estimated flux linkage will be more accurate.</p>
      <p id="d1e1279">By combining Eqs. (2) and (9), the two-phase current in the static coordinate system calculated by estimated flux linkage can be expressed as
            <disp-formula id="Ch1.E10" content-type="numbered"><label>10</label><mml:math id="M26" display="block"><mml:mrow><mml:mfenced open="{" close=""><mml:mtable rowspacing="5.690551pt" class="array" columnalign="left"><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>i</mml:mi><mml:mo stretchy="false" mathvariant="normal">^</mml:mo></mml:mover><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="false"><mml:mstyle displaystyle="false"><mml:mfrac style="text"><mml:mn mathvariant="normal">1</mml:mn><mml:mrow><mml:msub><mml:mi>L</mml:mi><mml:mi mathvariant="normal">s</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mstyle></mml:mstyle><mml:mfenced open="(" close=")"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi mathvariant="italic">λ</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="normal">r</mml:mi></mml:msub><mml:mi>cos⁡</mml:mi><mml:mi mathvariant="italic">θ</mml:mi></mml:mrow></mml:mfenced><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>i</mml:mi><mml:mo stretchy="false" mathvariant="normal">^</mml:mo></mml:mover><mml:mi mathvariant="italic">β</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mn mathvariant="normal">1</mml:mn><mml:mrow><mml:msub><mml:mi>L</mml:mi><mml:mi mathvariant="normal">s</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mstyle><mml:mfenced open="(" close=")"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi mathvariant="italic">λ</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover><mml:mi mathvariant="italic">β</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="normal">r</mml:mi></mml:msub><mml:mi>sin⁡</mml:mi><mml:mi mathvariant="italic">θ</mml:mi></mml:mrow></mml:mfenced><mml:mo>.</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mfenced></mml:mrow></mml:math></disp-formula>
          For various errors caused by nonlinear factors in Eq. (6), a sliding-mode observer is established using the estimated current and actual current to calculate the compensation voltage.<?pagebreak page102?> The sliding-mode function and sliding-mode observer are expressed as

                <disp-formula specific-use="gather" content-type="numbered"><mml:math id="M27" display="block"><mml:mtable displaystyle="true"><mml:mlabeledtr id="Ch1.E11"><mml:mtd><mml:mtext>11</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:mfenced close="" open="{"><mml:mtable class="array" rowspacing="2.845276pt" columnalign="left"><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>s</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mover accent="true"><mml:mi>i</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mi>i</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>s</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mover accent="true"><mml:mi>i</mml:mi><mml:mo stretchy="false" mathvariant="normal">^</mml:mo></mml:mover><mml:mi mathvariant="italic">β</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mi>i</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub><mml:mo>.</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mfenced></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E12"><mml:mtd><mml:mtext>12</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:mfenced close="" open="{"><mml:mtable class="array" rowspacing="2.845276pt" columnalign="left"><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>e</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mi>k</mml:mi><mml:mrow class="chem"><mml:mi mathvariant="normal">sign</mml:mi></mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>s</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>)</mml:mo><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>e</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mi>k</mml:mi><mml:mrow class="chem"><mml:mi mathvariant="normal">sign</mml:mi></mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>s</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub><mml:mo>)</mml:mo><mml:mo>.</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mfenced></mml:mrow></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math></disp-formula>

            When the error in the system is suppressed by the compensation voltage, the estimated flux linkage will be consistent with the actual value, and the estimated current will also be consistent with the actual value. The improved flux linkage estimation scheme is shown in Fig. 3. The proposed observer is simple in structure, is easy to implement, and has only one parameter. It is very suitable for practical applications and has significant advantages.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F3"><?xmltex \currentcnt{3}?><?xmltex \def\figurename{Figure}?><label>Figure 3</label><caption><p id="d1e1516">Position estimation based on improved flux linkage.</p></caption>
          <?xmltex \igopts{width=241.848425pt}?><graphic xlink:href="https://ms.copernicus.org/articles/15/99/2024/ms-15-99-2024-f03.png"/>

        </fig>

</sec>
<sec id="Ch1.S3.SS2">
  <label>3.2</label><title>Proof of stability</title>
      <p id="d1e1533">Taking the <inline-formula><mml:math id="M28" display="inline"><mml:mi mathvariant="italic">α</mml:mi></mml:math></inline-formula> axis as an example, the stability of the flux linkage observer is proved by the Lyapunov function. We construct the function and its first derivative as

                <disp-formula specific-use="gather" content-type="numbered"><mml:math id="M29" display="block"><mml:mtable displaystyle="true"><mml:mlabeledtr id="Ch1.E13"><mml:mtd><mml:mtext>13</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:msub><mml:mi>V</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mn mathvariant="normal">1</mml:mn><mml:mn mathvariant="normal">2</mml:mn></mml:mfrac></mml:mstyle><mml:msubsup><mml:mi>s</mml:mi><mml:mi mathvariant="italic">α</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msubsup><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E14"><mml:mtd><mml:mtext>14</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:msub><mml:mover accent="true"><mml:mi>V</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>s</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:msub><mml:mover accent="true"><mml:mi>s</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>.</mml:mo></mml:mrow></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math></disp-formula>

            Obviously, the <inline-formula><mml:math id="M30" display="inline"><mml:mrow><mml:msub><mml:mi>V</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is greater than zero.</p>
      <p id="d1e1621">According to Eqs. (9) and (10), the estimated current is
            <disp-formula id="Ch1.E15" content-type="numbered"><label>15</label><mml:math id="M31" display="block"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>i</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mn mathvariant="normal">1</mml:mn><mml:mrow><mml:msub><mml:mi>L</mml:mi><mml:mi mathvariant="normal">s</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mstyle><mml:mfenced open="[" close="]"><mml:mrow><mml:mo movablelimits="false">∫</mml:mo><mml:mfenced open="(" close=")"><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:mi>R</mml:mi><mml:msub><mml:mi>i</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mi>e</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:mfenced><mml:mi mathvariant="normal">d</mml:mi><mml:mi>t</mml:mi><mml:mo>+</mml:mo><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mrow><mml:mi mathvariant="italic">α</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="normal">r</mml:mi></mml:msub><mml:mi>cos⁡</mml:mi><mml:mi mathvariant="italic">θ</mml:mi></mml:mrow></mml:mfenced><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula>
          We substitute Eqs. (11) and (12) into Eq. (15):
            <disp-formula id="Ch1.E16" content-type="numbered"><label>16</label><mml:math id="M32" display="block"><mml:mtable rowspacing="0.2ex" class="split" displaystyle="true" columnalign="right left"><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>s</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>=</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mn mathvariant="normal">1</mml:mn><mml:mrow><mml:msub><mml:mi>L</mml:mi><mml:mi mathvariant="normal">s</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mstyle><mml:mo movablelimits="false">∫</mml:mo><mml:mfenced close="]" open="["><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:mi>R</mml:mi><mml:msub><mml:mi>i</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:mi>k</mml:mi><mml:mrow class="chem"><mml:mi mathvariant="normal">sign</mml:mi></mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>s</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:mfenced><mml:mi mathvariant="normal">d</mml:mi><mml:mi>t</mml:mi><mml:mo>+</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mrow><mml:mi mathvariant="italic">α</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>L</mml:mi><mml:mi mathvariant="normal">s</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mstyle></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd/><mml:mtd><mml:mrow><mml:mo>-</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="normal">r</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>L</mml:mi><mml:mi mathvariant="normal">s</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mstyle><mml:mi>cos⁡</mml:mi><mml:mi mathvariant="italic">θ</mml:mi><mml:mo>-</mml:mo><mml:msub><mml:mi>i</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>.</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
          We differentiate Eq. (16):
            <disp-formula id="Ch1.E17" content-type="numbered"><label>17</label><mml:math id="M33" display="block"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>s</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mn mathvariant="normal">1</mml:mn><mml:mrow><mml:msub><mml:mi>L</mml:mi><mml:mi mathvariant="normal">s</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mstyle><mml:msub><mml:mi>u</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mn mathvariant="normal">1</mml:mn><mml:mrow><mml:msub><mml:mi>L</mml:mi><mml:mi mathvariant="normal">s</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mstyle><mml:mi>R</mml:mi><mml:msub><mml:mi>i</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mn mathvariant="normal">1</mml:mn><mml:mrow><mml:msub><mml:mi>L</mml:mi><mml:mi mathvariant="normal">s</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mstyle><mml:mi>k</mml:mi><mml:mrow class="chem"><mml:mi mathvariant="normal">sign</mml:mi></mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>s</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>)</mml:mo><mml:mo>+</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:mi mathvariant="italic">ω</mml:mi><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="normal">r</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>L</mml:mi><mml:mi mathvariant="normal">s</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mstyle><mml:mi>sin⁡</mml:mi><mml:mi mathvariant="italic">θ</mml:mi><mml:mo>-</mml:mo><mml:msub><mml:mover accent="true"><mml:mi>i</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula>
          We substitute Eq. (17) into Eq. (14):
            <disp-formula id="Ch1.E18" content-type="numbered"><label>18</label><mml:math id="M34" display="block"><mml:mtable class="split" rowspacing="0.2ex" displaystyle="true" columnalign="right left"><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>V</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>=</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msub><mml:mi>s</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:msub><mml:mover accent="true"><mml:mi>s</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mo>=</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mfenced close=")" open="("><mml:mrow><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mn mathvariant="normal">1</mml:mn><mml:mrow><mml:msub><mml:mi>L</mml:mi><mml:mi mathvariant="normal">s</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mstyle><mml:msub><mml:mi>u</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mn mathvariant="normal">1</mml:mn><mml:mrow><mml:msub><mml:mi>L</mml:mi><mml:mi mathvariant="normal">s</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mstyle><mml:mi>R</mml:mi><mml:msub><mml:mi>i</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:mi mathvariant="italic">ω</mml:mi><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="normal">r</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>L</mml:mi><mml:mi mathvariant="normal">s</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mstyle><mml:mi>sin⁡</mml:mi><mml:mi mathvariant="italic">θ</mml:mi><mml:mo>-</mml:mo><mml:msub><mml:mover accent="true"><mml:mi>i</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:mfenced><mml:msub><mml:mi>s</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd/><mml:mtd><mml:mrow><mml:mo>-</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mn mathvariant="normal">1</mml:mn><mml:mrow><mml:msub><mml:mi>L</mml:mi><mml:mi mathvariant="normal">s</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mstyle><mml:mi>k</mml:mi><mml:mfenced close="|" open="|"><mml:mrow><mml:msub><mml:mi>s</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:mfenced><mml:mo>.</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
          From Eq. (18), if <inline-formula><mml:math id="M35" display="inline"><mml:mrow><mml:mover accent="true"><mml:mi>V</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mo>&lt;</mml:mo><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:math></inline-formula>, there is
            <disp-formula id="Ch1.E19" content-type="numbered"><label>19</label><mml:math id="M36" display="block"><mml:mrow><mml:mi>k</mml:mi><mml:mo>&gt;</mml:mo><mml:mfenced open="(" close=")"><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:mi>R</mml:mi><mml:msub><mml:mi>i</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:mi mathvariant="italic">ω</mml:mi><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="normal">r</mml:mi></mml:msub><mml:mi>sin⁡</mml:mi><mml:mi mathvariant="italic">θ</mml:mi><mml:mo>-</mml:mo><mml:msub><mml:mi>L</mml:mi><mml:mi mathvariant="normal">s</mml:mi></mml:msub><mml:msub><mml:mover accent="true"><mml:mi>i</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:mfenced><mml:mrow class="chem"><mml:mi mathvariant="normal">sign</mml:mi></mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>s</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>)</mml:mo><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula>
          When the sliding-mode coefficient <inline-formula><mml:math id="M37" display="inline"><mml:mi>k</mml:mi></mml:math></inline-formula> satisfies Eq. (19), the compensation voltage obtained by the constructed sliding-mode observer will make the estimated current converge with the actual current, and the estimated flux linkage will converge with the actual flux linkage.</p>
      <p id="d1e2167">According to Eqs. (1) and (2),
            <disp-formula id="Ch1.E20" content-type="numbered"><label>20</label><mml:math id="M38" display="block"><mml:mtable rowspacing="0.2ex" class="split" displaystyle="true" columnalign="right left"><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi mathvariant="italic">λ</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>=</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:mi>R</mml:mi><mml:msub><mml:mi>i</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mo>=</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msub><mml:mi>L</mml:mi><mml:mi mathvariant="normal">s</mml:mi></mml:msub><mml:msub><mml:mover accent="true"><mml:mi>i</mml:mi><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:mi mathvariant="italic">ω</mml:mi><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="normal">r</mml:mi></mml:msub><mml:mi>sin⁡</mml:mi><mml:mi mathvariant="italic">θ</mml:mi><mml:mo>.</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
          Substituting Eq. (20) into Eq. (19), we can get
            <disp-formula id="Ch1.E21" content-type="numbered"><label>21</label><mml:math id="M39" display="block"><mml:mrow><mml:mi>k</mml:mi><mml:mo>&gt;</mml:mo><mml:mn mathvariant="normal">0</mml:mn><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula>
          In addition, according to Eqs. (9) and (12), the estimated flux linkage <inline-formula><mml:math id="M40" display="inline"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi mathvariant="italic">λ</mml:mi><mml:mo stretchy="true" mathvariant="normal">^</mml:mo></mml:mover><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> and the compensation voltage <inline-formula><mml:math id="M41" display="inline"><mml:mrow><mml:msub><mml:mi>e</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> have
            <disp-formula id="Ch1.E22" content-type="numbered"><label>22</label><mml:math id="M42" display="block"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mover accent="true"><mml:mi mathvariant="italic">λ</mml:mi><mml:mo stretchy="true" mathvariant="normal">^</mml:mo></mml:mover><mml:mo mathvariant="normal">˙</mml:mo></mml:mover><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>u</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:mi>R</mml:mi><mml:msub><mml:mi>i</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:mi>k</mml:mi><mml:mrow class="chem"><mml:mi mathvariant="normal">sign</mml:mi></mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>s</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>)</mml:mo><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula>
          During the operation of the motor, the stator flux linkage always changes in a sine or cosine manner. The <inline-formula><mml:math id="M43" display="inline"><mml:mi mathvariant="italic">α</mml:mi></mml:math></inline-formula>-axis estimated flux linkage can be expressed as
            <disp-formula id="Ch1.E23" content-type="numbered"><label>23</label><mml:math id="M44" display="block"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi mathvariant="italic">λ</mml:mi><mml:mo stretchy="false" mathvariant="normal">^</mml:mo></mml:mover><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mover accent="true"><mml:mi mathvariant="italic">λ</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover><mml:mi mathvariant="normal">r</mml:mi></mml:msub><mml:mi>cos⁡</mml:mi><mml:mfenced close=")" open="("><mml:mrow><mml:mi mathvariant="italic">ω</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfenced><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula>
          From Eqs. (22) and (23), there is
            <disp-formula id="Ch1.E24" content-type="numbered"><label>24</label><mml:math id="M45" display="block"><mml:mrow><mml:mi>k</mml:mi><mml:mrow class="chem"><mml:mi mathvariant="normal">sign</mml:mi></mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>s</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>)</mml:mo><mml:mo>=</mml:mo><mml:msub><mml:mi>u</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:mi>R</mml:mi><mml:msub><mml:mi>i</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:mi mathvariant="italic">ω</mml:mi><mml:msub><mml:mover accent="true"><mml:mi mathvariant="italic">λ</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover><mml:mi mathvariant="normal">r</mml:mi></mml:msub><mml:mi>sin⁡</mml:mi><mml:mfenced open="(" close=")"><mml:mrow><mml:mi mathvariant="italic">ω</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfenced><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula>
          As the sliding-mode coefficient approaches or exceeds the maximum amplitude of the back EMF, the compensation voltage will become the main signal, and its effect will be greater than that of the back EMF. To ensure the main function of the back EMF, the sliding-mode coefficient <inline-formula><mml:math id="M46" display="inline"><mml:mi>k</mml:mi></mml:math></inline-formula> should also satisfy <inline-formula><mml:math id="M47" display="inline"><mml:mrow><mml:mi>k</mml:mi><mml:mo>&lt;</mml:mo><mml:mi mathvariant="italic">ω</mml:mi><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="normal">r</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>. Therefore, if <inline-formula><mml:math id="M48" display="inline"><mml:mrow><mml:mn mathvariant="normal">0</mml:mn><mml:mo>&lt;</mml:mo><mml:mi>k</mml:mi><mml:mo>&lt;</mml:mo><mml:mi mathvariant="italic">ω</mml:mi><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="normal">r</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, the compensation voltage obtained by the sliding-mode observer can make the estimated flux linkage <inline-formula><mml:math id="M49" display="inline"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi mathvariant="italic">λ</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> converge to the actual value <inline-formula><mml:math id="M50" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, and the proof of the <inline-formula><mml:math id="M51" display="inline"><mml:mi mathvariant="italic">β</mml:mi></mml:math></inline-formula> axis is same.</p>
      <p id="d1e2510">The improved flux linkage observer proposed in this paper has a good suppression effect on the inaccurate position estimation caused by the high-frequency noise of voltage and current, uncertain interference, and inaccurate initial position and has a good anti-interference capacity and stability. Next, the proposed observer will be verified by simulation and experiments.</p>
</sec>
</sec>
<?pagebreak page103?><sec id="Ch1.S4">
  <label>4</label><title>Simulation results</title>
      <p id="d1e2522">To further examine the actual performance of the improved observer, some simulations are carried out according to the mathematical model of PMSLMs. The MATLAB/Simulink is used in this paper, and its version is 2016b. The common space vector pulse width modulation (SVPWM) drive scheme is used to control the motor, and the frequency of the inverter is 10 kHz. The controller adopts a simple PI controller, including speed loop control and current loop control. The control system is shown in Fig. 4, and the motor parameters are shown in Table 1. Actually, the proposed flux linkage observer only needs the voltage and current of the <inline-formula><mml:math id="M52" display="inline"><mml:mrow><mml:mi mathvariant="italic">α</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:mrow></mml:math></inline-formula> axis to estimate the position. The current data are sampled by the current sensor. The voltages <inline-formula><mml:math id="M53" display="inline"><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M54" display="inline"><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> are replaced by the control voltages <inline-formula><mml:math id="M55" display="inline"><mml:mrow><mml:msubsup><mml:mi>u</mml:mi><mml:mi mathvariant="italic">α</mml:mi><mml:mo>∗</mml:mo></mml:msubsup></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M56" display="inline"><mml:mrow><mml:msubsup><mml:mi>u</mml:mi><mml:mi mathvariant="italic">β</mml:mi><mml:mo>∗</mml:mo></mml:msubsup></mml:mrow></mml:math></inline-formula>. Due to the nonlinear effect of the inverter, there is a delay and error between the expected voltage and the actual voltage. Much of the literature has studied the characteristics of inverter and proposed voltage compensation schemes. The flux linkage observer proposed in this paper classifies it as voltage error disturbance, which can also be well compensated for by the nonlinear sliding-mode observer and can be easily implemented.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F4"><?xmltex \currentcnt{4}?><?xmltex \def\figurename{Figure}?><label>Figure 4</label><caption><p id="d1e2585">Control system of PMSLMs.</p></caption>
        <?xmltex \igopts{width=241.848425pt}?><graphic xlink:href="https://ms.copernicus.org/articles/15/99/2024/ms-15-99-2024-f04.png"/>

      </fig>

<?xmltex \floatpos{t}?><table-wrap id="Ch1.T1"><?xmltex \currentcnt{1}?><label>Table 1</label><caption><p id="d1e2597">Parameters of PMSLMs.</p></caption><oasis:table frame="topbot"><oasis:tgroup cols="2">
     <oasis:colspec colnum="1" colname="col1" align="left"/>
     <oasis:colspec colnum="2" colname="col2" align="left"/>
     <oasis:thead>
       <oasis:row rowsep="1">
         <oasis:entry colname="col1">Parameter</oasis:entry>
         <oasis:entry colname="col2">Value</oasis:entry>
       </oasis:row>
     </oasis:thead>
     <oasis:tbody>
       <oasis:row>
         <oasis:entry colname="col1">Pole pair (<inline-formula><mml:math id="M57" display="inline"><mml:mrow><mml:msub><mml:mi>p</mml:mi><mml:mi>n</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>)</oasis:entry>
         <oasis:entry colname="col2">2</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">Pole pitch (<inline-formula><mml:math id="M58" display="inline"><mml:mi mathvariant="italic">τ</mml:mi></mml:math></inline-formula>)</oasis:entry>
         <oasis:entry colname="col2">0.04 mm</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">Permanent magnet flux linkage (<inline-formula><mml:math id="M59" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="normal">r</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>)</oasis:entry>
         <oasis:entry colname="col2">0.3 Wb</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">Phase inductance (<inline-formula><mml:math id="M60" display="inline"><mml:mi>L</mml:mi></mml:math></inline-formula>)</oasis:entry>
         <oasis:entry colname="col2">15 mH</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">Phase resistance (<inline-formula><mml:math id="M61" display="inline"><mml:mi>R</mml:mi></mml:math></inline-formula>)</oasis:entry>
         <oasis:entry colname="col2">9.3 <inline-formula><mml:math id="M62" display="inline"><mml:mi mathvariant="normal">Ω</mml:mi></mml:math></inline-formula></oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">Mover mass (<inline-formula><mml:math id="M63" display="inline"><mml:mi>M</mml:mi></mml:math></inline-formula>)</oasis:entry>
         <oasis:entry colname="col2">1.9 kg</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">Rate current (<inline-formula><mml:math id="M64" display="inline"><mml:mi>I</mml:mi></mml:math></inline-formula>)</oasis:entry>
         <oasis:entry colname="col2">1.5 A</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">Shaft length</oasis:entry>
         <oasis:entry colname="col2">1200 mm</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">Moving-part length</oasis:entry>
         <oasis:entry colname="col2">220 mm</oasis:entry>
       </oasis:row>
     </oasis:tbody>
   </oasis:tgroup></oasis:table><?xmltex \gdef\@currentlabel{1}?></table-wrap>

      <p id="d1e2770">To simulate the possible current noise in actuality, a random disturbance with amplitude 1 A is added to the current<inline-formula><mml:math id="M65" display="inline"><mml:mrow><mml:msub><mml:mi>i</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>. The actual current and the current added random noise are shown in Fig. 5. It can be found that random noise accounts for a large proportion of the output current. Actually, the current noise will not be so large; this is to simulate a sufficiently bad situation. The initial position of the <inline-formula><mml:math id="M66" display="inline"><mml:mi>d</mml:mi></mml:math></inline-formula> axis is zero at the beginning – i.e., <inline-formula><mml:math id="M67" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="normal">r</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M68" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:math></inline-formula> in Eq. (3). This can also be expressed as <inline-formula><mml:math id="M69" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="normal">r</mml:mi></mml:msub><mml:mi>cos⁡</mml:mi><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M70" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="normal">r</mml:mi></mml:msub><mml:mi>sin⁡</mml:mi><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M71" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:math></inline-formula>. In fact, the initial position of the motor cannot always be obtained accurately, which results in the difference between the integral initial value of the stator flux linkage set by the controller and the actual value. It makes the estimated flux linkage difficult to calculate accurately. For this problem, the estimated position of the <inline-formula><mml:math id="M72" display="inline"><mml:mi>d</mml:mi></mml:math></inline-formula> axis is set as 60<inline-formula><mml:math id="M73" display="inline"><mml:msup><mml:mi/><mml:mo>∘</mml:mo></mml:msup></mml:math></inline-formula> in the simulation. Then, the initial flux linkage errors <inline-formula><mml:math id="M74" display="inline"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mrow><mml:mi mathvariant="italic">α</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M75" display="inline"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mrow><mml:mi mathvariant="italic">β</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> in Eq. (6) can be obtained. The convention voltage or current models demonstrate difficulty in dealing with these adverse conditions effectively. The estimated flux linkages <inline-formula><mml:math id="M76" display="inline"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi mathvariant="italic">λ</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover><mml:mrow><mml:mi mathvariant="normal">r</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M77" display="inline"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi mathvariant="italic">λ</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover><mml:mrow><mml:mi mathvariant="normal">r</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> under the convention way are shown in Fig. 6, which shows that the error is very large. The estimated speed and position are completely unusable.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F5"><?xmltex \currentcnt{5}?><?xmltex \def\figurename{Figure}?><label>Figure 5</label><caption><p id="d1e2975">The <inline-formula><mml:math id="M78" display="inline"><mml:mi mathvariant="italic">α</mml:mi></mml:math></inline-formula>-axis actual current and the current with noise.</p></caption>
        <?xmltex \igopts{width=241.848425pt}?><graphic xlink:href="https://ms.copernicus.org/articles/15/99/2024/ms-15-99-2024-f05.png"/>

      </fig>

      <?xmltex \floatpos{t}?><fig id="Ch1.F6"><?xmltex \currentcnt{6}?><?xmltex \def\figurename{Figure}?><label>Figure 6</label><caption><p id="d1e2993">Estimated flux linkage with convention voltage model.</p></caption>
        <?xmltex \igopts{width=241.848425pt}?><graphic xlink:href="https://ms.copernicus.org/articles/15/99/2024/ms-15-99-2024-f06.png"/>

      </fig>

      <?pagebreak page104?><p id="d1e3002"><?xmltex \hack{\newpage}?>Using the improved flux linkage observer proposed in this paper, simulations are performed under constant load and uniform speed. The estimation of the flux linkage is shown in Fig. 7. Under the influence of the current noise and the error of the integral initial value, the estimated value of the flux linkage can still converge to the true value and stabilize quickly. The estimated speed of the PMSLM is shown in Fig. 8. In the initial stage of motor start-up, the estimated speed can quickly and stably converge to the actual value after a short adjustment. The estimated position and position errors are shown in Fig. 9. Similarly, in the initial stage of start-up, the position error is large. After a short adjustment, the estimated position can keep up with the actual position steadily, and the error can converge to near zero.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F7"><?xmltex \currentcnt{7}?><?xmltex \def\figurename{Figure}?><label>Figure 7</label><caption><p id="d1e3009">Estimated flux linkage with the improved flux linkage observer.</p></caption>
        <?xmltex \igopts{width=241.848425pt}?><graphic xlink:href="https://ms.copernicus.org/articles/15/99/2024/ms-15-99-2024-f07.png"/>

      </fig>

      <?xmltex \floatpos{t}?><fig id="Ch1.F8"><?xmltex \currentcnt{8}?><?xmltex \def\figurename{Figure}?><label>Figure 8</label><caption><p id="d1e3020">Estimated speed with the improved flux linkage observer.</p></caption>
        <?xmltex \igopts{width=227.622047pt}?><graphic xlink:href="https://ms.copernicus.org/articles/15/99/2024/ms-15-99-2024-f08.png"/>

      </fig>

      <?xmltex \floatpos{t}?><fig id="Ch1.F9"><?xmltex \currentcnt{9}?><?xmltex \def\figurename{Figure}?><label>Figure 9</label><caption><p id="d1e3031">Estimated position and error with the improved flux linkage observer.</p></caption>
        <?xmltex \igopts{width=241.848425pt}?><graphic xlink:href="https://ms.copernicus.org/articles/15/99/2024/ms-15-99-2024-f09.png"/>

      </fig>

      <p id="d1e3040">Because of the fixed coefficient <inline-formula><mml:math id="M79" display="inline"><mml:mi>k</mml:mi></mml:math></inline-formula> in the sliding-mode observer, the current error is difficult to maintain at zero all the time when approaching the steady state. Therefore, the compensation voltage has high-frequency chattering. This is also the common issue with the sliding-mode theory. However, the integration process in flux linkage calculation has largely eliminated this high-frequency chattering. We can find that the sliding-mode observer combined with the flux linkage estimation is very effective. There is no need to do more processing for the chattering. This is also the advantage of the improved scheme in this paper. As the coefficient <inline-formula><mml:math id="M80" display="inline"><mml:mi>k</mml:mi></mml:math></inline-formula> is small, it means that more fine adjustment can be obtained, but it will also affect the convergence speed.</p>
      <p id="d1e3057">Under the condition of variable speeds, the proposed observer is simulated and tested. The speed of the PMSLM changes from 0.5 to 0.7 m s<inline-formula><mml:math id="M81" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula> at 1 s. As shown in Fig. 10, the estimated flux linkage converges to the actual value after a short adjustment at the beginning. At 1 s, the sudden change in speed changes the estimated flux linkage, but it quickly stabilizes again. The estimated speed is shown in Fig. 11, and it can track the actual speed well. When the actual speed changes, the estimated speed also has a good dynamic respond. Figure 12 shows the estimated position and position error. The sudden change in the speed leads to a small deviation of the position estimation, and it can converge with the actual position quickly. The estimated speed curve under load change is depicted in Fig. 13. The load force has a step change at 1 s. When the load suddenly changes, the estimated speed tracks the actual speed well, and the system returns to stability quickly. Figure 14 shows the estimated position and error. It can be found that the load change has no significant impact on the position estimation. Due to the strong nonlinearity characteristic of the sign function in the sliding-mode observer, the strong anti-disturbance capability is retained. Although the sign function that seems simple and rough will cause chattering in most applications, the integral operation in flux linkage calculation effectively makes up for the shortcoming. These all show that the sensorless scheme has a strong anti-interference ability.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F10"><?xmltex \currentcnt{10}?><?xmltex \def\figurename{Figure}?><label>Figure 10</label><caption><p id="d1e3075">Estimated flux linkage with the improved flux linkage observer under variable speeds.</p></caption>
        <?xmltex \igopts{width=241.848425pt}?><graphic xlink:href="https://ms.copernicus.org/articles/15/99/2024/ms-15-99-2024-f10.png"/>

      </fig>

      <?xmltex \floatpos{t}?><fig id="Ch1.F11"><?xmltex \currentcnt{11}?><?xmltex \def\figurename{Figure}?><label>Figure 11</label><caption><p id="d1e3086">Estimated speed with the improved flux linkage observer under variable speeds.</p></caption>
        <?xmltex \igopts{width=227.622047pt}?><graphic xlink:href="https://ms.copernicus.org/articles/15/99/2024/ms-15-99-2024-f11.png"/>

      </fig>

      <?xmltex \floatpos{t}?><fig id="Ch1.F12"><?xmltex \currentcnt{12}?><?xmltex \def\figurename{Figure}?><label>Figure 12</label><caption><p id="d1e3097">Estimated position and error with the improved flux linkage observer under variable speeds.</p></caption>
        <?xmltex \igopts{width=241.848425pt}?><graphic xlink:href="https://ms.copernicus.org/articles/15/99/2024/ms-15-99-2024-f12.png"/>

      </fig>

      <?xmltex \floatpos{t}?><fig id="Ch1.F13"><?xmltex \currentcnt{13}?><?xmltex \def\figurename{Figure}?><label>Figure 13</label><caption><p id="d1e3108">Estimated speed with the improved flux linkage observer under variable loads.</p></caption>
        <?xmltex \igopts{width=227.622047pt}?><graphic xlink:href="https://ms.copernicus.org/articles/15/99/2024/ms-15-99-2024-f13.png"/>

      </fig>

      <?xmltex \floatpos{t}?><fig id="Ch1.F14"><?xmltex \currentcnt{14}?><?xmltex \def\figurename{Figure}?><label>Figure 14</label><caption><p id="d1e3120">Estimated position and error with the improved flux linkage observer under variable loads.</p></caption>
        <?xmltex \igopts{width=241.848425pt}?><graphic xlink:href="https://ms.copernicus.org/articles/15/99/2024/ms-15-99-2024-f14.png"/>

      </fig>

<?xmltex \hack{\newpage}?>
</sec>
<?pagebreak page105?><sec id="Ch1.S5">
  <label>5</label><title>Experiment results</title>
      <p id="d1e3140">The experiment system used in this paper is shown in Fig. 15, including the tubular PMSLM, position sensor, oscilloscope, computer, controller and driver. The controller and driver are integrated together and communicate with the computer in real time. The control algorithm is set and adjusted through the computer software. The whole control system also adopts the PI control algorithm that is commonly used in practice to establish the current control loop and speed control loop. The position sensor is used to detect the position of the moving part. Commonly, the information from position sensor will feedback to the controller. However, in this paper, the information on position and speed is from the improved sensorless scheme.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F15"><?xmltex \currentcnt{15}?><?xmltex \def\figurename{Figure}?><label>Figure 15</label><caption><p id="d1e3145">Experiment system of the PMSLM in this paper.</p></caption>
        <?xmltex \igopts{width=241.848425pt}?><graphic xlink:href="https://ms.copernicus.org/articles/15/99/2024/ms-15-99-2024-f15.jpg"/>

      </fig>

      <p id="d1e3154">The improved flux linkage observer proposed in this paper is established with the voltage and current data collected by the built-in sensor to estimate the position and speed of the PMSLM. The estimated value is compared with the actual position and speed obtained by the position sensor. Oscilloscope is used to display relevant signals for observation and debugging. Due to the limit of the length of the PMSLM and the experiment equipment, the performance of the proposed flux linkage observer can only be verified under different speeds. The overall movement time is 1 s, and the motor movement speed is changed at 0.5 s. There is an artificial deviation between the motor's actual initial position and the estimated initial position. The flux linkage, speed and position estimations of the motor are tested under the action of the improved flux linkage observer when the permanent magnet synchronous linear motor (PMSLM) is in operation.</p>
      <p id="d1e3158">The inaccuracy of the initial position will cause error in the calculation of the flux linkage during operation. The estimated flux linkage using the convention voltage model without a flux linkage observer is shown in Fig. 16. It can be found that the flux linkage always has a deviation and is<?pagebreak page106?> asymmetrical. This greatly affects the performance of the estimated position and speed. Figure 17 shows the variation of the estimated flux linkage under the action of the proposed observer. Obviously, the initial error can be greatly compensated for by the improved flux linkage observer, and the estimated flux linkage converges with the actual value quickly. At 0.5 s, the controller makes the motor speed change. At the same time, the estimated flux linkage is also adjusted by the observer and soon returns to stability again. It can be found that the improved flux linkage observer has a good compensation effect on the rotor flux linkage estimation error caused by the initial-position error.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F16"><?xmltex \currentcnt{16}?><?xmltex \def\figurename{Figure}?><label>Figure 16</label><caption><p id="d1e3163">The estimated flux linkage with the inaccurate initial-phase angle under the convention voltage model.</p></caption>
        <?xmltex \igopts{width=241.848425pt}?><graphic xlink:href="https://ms.copernicus.org/articles/15/99/2024/ms-15-99-2024-f16.png"/>

      </fig>

      <?xmltex \floatpos{t}?><fig id="Ch1.F17"><?xmltex \currentcnt{17}?><?xmltex \def\figurename{Figure}?><label>Figure 17</label><caption><p id="d1e3174">The estimated flux linkage with the inaccurate initial-phase angle under the improved flux linkage observer.</p></caption>
        <?xmltex \igopts{width=241.848425pt}?><graphic xlink:href="https://ms.copernicus.org/articles/15/99/2024/ms-15-99-2024-f17.png"/>

      </fig>

      <p id="d1e3183">Figures 18 and 19 show the speed and position estimation of the PMSLM based on the corrected rotor flux linkage. In Fig. 18, at the beginning, the estimated speed fluctuates slightly around the actual speed. This is because the estimated flux linkage is inaccurate, and the estimated speed is stable until the initial flux linkage error is well compensated for. At 0.5 s, as the actual speed changes, the estimated speed can still track the actual speed quickly and converge to the actual speed after a short fluctuation. It can be clearly found that the speed fluctuation at 0.5 s is smaller and shorter. This shows that the proposed observer has a good performance in terms of dynamic response after the initial error is compensated for and can respond to system changes in time. The position estimation shown in Fig. 19 also has the same change. At the beginning, the estimated position is inaccurate. As the initial flux linkage error is compensated for, the estimated position constantly approaches the real position. At 0.5 s, the position error also has a small sudden change with the motor speed change, and it returns to stability again. The steady error of the estimated position is about 2 mm, and it gradually approaches the real position.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F18"><?xmltex \currentcnt{18}?><?xmltex \def\figurename{Figure}?><label>Figure 18</label><caption><p id="d1e3189">The actual speed and estimated speed under variable speed conditions.</p></caption>
        <?xmltex \igopts{width=227.622047pt}?><graphic xlink:href="https://ms.copernicus.org/articles/15/99/2024/ms-15-99-2024-f18.png"/>

      </fig>

      <?xmltex \floatpos{t}?><fig id="Ch1.F19"><?xmltex \currentcnt{19}?><?xmltex \def\figurename{Figure}?><label>Figure 19</label><caption><p id="d1e3200">The actual position and estimated position under variable speed conditions.</p></caption>
        <?xmltex \igopts{width=241.848425pt}?><graphic xlink:href="https://ms.copernicus.org/articles/15/99/2024/ms-15-99-2024-f19.png"/>

      </fig>

</sec>
<sec id="Ch1.S6" sec-type="conclusions">
  <label>6</label><title>Conclusions</title>
      <p id="d1e3217">In this paper, an improved flux linkage observer of PMSLM is proposed, which can effectively suppress the speed and<?pagebreak page107?> position estimation errors caused by the inaccurate initial position, high-frequency current noise and voltage error that may exist in practice. Firstly, the voltage model is used to calculate the stator flux linkage preliminarily, and the estimated current is calculated through the preliminary stator flux linkage using the current model. Then, a sliding-mode observer is established according to the estimated current and the actual current. The output of the observer is used as a compensation part of the voltage model to further correct the flux linkage. The theoretical model of the improved flux linkage observer is analyzed. The range of the coefficient in the observer is analyzed, and the proof of system stability is also shown in this paper. At last, the simulation and experiments show that the improved observer can greatly compensate for the estimated error and have a good dynamic response ability, especially for the initial-position error in practice. The proposed improved flux linkage observer has a simple structure and only one parameter. The common issues with regard chattering in sliding mode also have been effectively suppressed. It is suitable for practical applications.</p>
</sec>

      
      </body>
    <back><app-group>

<app id="App1.Ch1.S1">
  <?xmltex \currentcnt{A}?><label>Appendix A</label><title/>
      <p id="d1e3230"><table-wrap id="Taba" position="anchor"><oasis:table><oasis:tgroup cols="2">
     <oasis:colspec colnum="1" colname="col1" align="left"/>
     <oasis:colspec colnum="2" colname="col2" align="justify" colwidth="5.8cm"/>
     <oasis:tbody>
       <oasis:row>
         <oasis:entry colname="col1"><inline-formula><mml:math id="M82" display="inline"><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M83" display="inline"><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M84" display="inline"><mml:mrow><mml:msub><mml:mi>i</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M85" display="inline"><mml:mrow><mml:msub><mml:mi>i</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col2"><inline-formula><mml:math id="M86" display="inline"><mml:mrow><mml:mi mathvariant="italic">α</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:mrow></mml:math></inline-formula>-axis voltage and current</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1"><inline-formula><mml:math id="M87" display="inline"><mml:mrow><mml:msub><mml:mi>e</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M88" display="inline"><mml:mrow><mml:msub><mml:mi>e</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col2">Compensation voltage</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1"><inline-formula><mml:math id="M89" display="inline"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi>u</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M90" display="inline"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi>u</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col2">Voltage error</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1"><inline-formula><mml:math id="M91" display="inline"><mml:mi>R</mml:mi></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col2">Phase resistance</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1"><inline-formula><mml:math id="M92" display="inline"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi>i</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M93" display="inline"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi>i</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col2">Current error</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1"><inline-formula><mml:math id="M94" display="inline"><mml:mrow><mml:msub><mml:mi>L</mml:mi><mml:mi mathvariant="normal">s</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col2"><inline-formula><mml:math id="M95" display="inline"><mml:mrow><mml:mi mathvariant="italic">α</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:mrow></mml:math></inline-formula>-axis inductance</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1"><inline-formula><mml:math id="M96" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M97" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col2">Stator flux linkage component</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1"><inline-formula><mml:math id="M98" display="inline"><mml:mi mathvariant="italic">θ</mml:mi></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col2">Estimated angle</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1"><inline-formula><mml:math id="M99" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mrow><mml:mi mathvariant="italic">α</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M100" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mrow><mml:mi mathvariant="italic">β</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col2">Initial value of flux linkage</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1"><inline-formula><mml:math id="M101" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col2">Initial-phase angle</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1"><inline-formula><mml:math id="M102" display="inline"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mrow><mml:mi mathvariant="italic">α</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M103" display="inline"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mrow><mml:mi mathvariant="italic">β</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col2">Initial-value error of flux linkage</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1"><inline-formula><mml:math id="M104" display="inline"><mml:mi mathvariant="italic">ω</mml:mi></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col2">Estimated speed</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1"><inline-formula><mml:math id="M105" display="inline"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi mathvariant="italic">λ</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M106" display="inline"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi mathvariant="italic">λ</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover><mml:mi mathvariant="italic">β</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col2">Estimated stator flux linkage component</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1"><inline-formula><mml:math id="M107" display="inline"><mml:mrow><mml:msub><mml:mi>s</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M108" display="inline"><mml:mrow><mml:msub><mml:mi>s</mml:mi><mml:mi mathvariant="italic">β</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col2">Sliding-mode function</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1"><inline-formula><mml:math id="M109" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">λ</mml:mi><mml:mi mathvariant="normal">r</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col2">Flux linkage of permanent magnets</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1"><inline-formula><mml:math id="M110" display="inline"><mml:mi>k</mml:mi></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col2">Sliding-mode gain</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1"><inline-formula><mml:math id="M111" display="inline"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi mathvariant="italic">λ</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover><mml:mi mathvariant="normal">r</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col2">Estimated permanent magnet flux linkage</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1"><inline-formula><mml:math id="M112" display="inline"><mml:mi mathvariant="italic">τ</mml:mi></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col2">Pole pitch</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1"><inline-formula><mml:math id="M113" display="inline"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>i</mml:mi><mml:mo stretchy="false" mathvariant="normal">^</mml:mo></mml:mover><mml:mi mathvariant="italic">α</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M114" display="inline"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>i</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover><mml:mi mathvariant="italic">β</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col2">Estimated current</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1"><inline-formula><mml:math id="M115" display="inline"><mml:mrow><mml:msub><mml:mi>p</mml:mi><mml:mi>n</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col2">Pole pair</oasis:entry>
       </oasis:row>
     </oasis:tbody>
   </oasis:tgroup></oasis:table></table-wrap></p>
</app>
  </app-group><notes notes-type="dataavailability"><title>Data availability</title>

      <p id="d1e3779">No data sets were used in this article.</p>
  </notes><notes notes-type="authorcontribution"><title>Author contributions</title>

      <p id="d1e3785">WY conducted theoretical analysis and verification and wrote the manuscript. GY and LW provided experimental assistance. DL and AAZ collected the experimental results and created the figures.</p>
  </notes><notes notes-type="competinginterests"><title>Competing interests</title>

      <p id="d1e3791">The contact author has declared that none of the authors has any competing interests.</p>
  </notes><notes notes-type="disclaimer"><title>Disclaimer</title>

      <p id="d1e3797">Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.</p>
  </notes><notes notes-type="financialsupport"><title>Financial support</title>

      <p id="d1e3803">This research has been supported by the National Postdoctoral Program for Innovative Talents (grant no. BX2021126), the National Natural Science Foundation of China (grant no. 52105106), and the China Postdoctoral Science Foundation (grant no. 2021M701711).</p>
  </notes><?xmltex \hack{\newpage}?><notes notes-type="reviewstatement"><title>Review statement</title>

      <p id="d1e3811">This paper was edited by Zi Bin and reviewed by two anonymous referees.</p>
  </notes><ref-list>
    <title>References</title>

      <ref id="bib1.bib1"><label>1</label><?label 1?><mixed-citation>Accetta, A., Cirrincione, M., Pucci, M., and Vitale, G.: Neural Sensorless Control of Linear Induction Motors by a Full-Order Luenberger Observer Considering the End Effects, IEEE Trans. Ind. Appl., 50, 1891–1904, <ext-link xlink:href="https://doi.org/10.1109/TIA.2013.2288429" ext-link-type="DOI">10.1109/TIA.2013.2288429</ext-link>, 2014.</mixed-citation></ref>
      <ref id="bib1.bib2"><label>2</label><?label 1?><mixed-citation>Afifa, R., Ali, S., Pervaiz, M., and Lqbal, J.: Adaptive backstepping integral sliding mode control of a mimo separately excited DC motor,  Robotics, 12, 105, <ext-link xlink:href="https://doi.org/10.3390/robotics12040105" ext-link-type="DOI">10.3390/robotics12040105</ext-link>, 2023.</mixed-citation></ref>
      <ref id="bib1.bib3"><label>3</label><?label 1?><mixed-citation>Ahmad, S., Uppal, A., Azam, M., and Lqbal, J.: Chattering free sliding mode control and state dependent Kalman filter design for underground gasification energy conversion process, Electronics, 12, 876, <ext-link xlink:href="https://doi.org/10.3390/electronics12040876" ext-link-type="DOI">10.3390/electronics12040876</ext-link>, 2023.</mixed-citation></ref>
      <ref id="bib1.bib4"><label>4</label><?label 1?><mixed-citation>Cao, R., Cheng, M., and Zhang, B.: Speed control of complementary and modular linear flux-switching permanent-magnet motor, IEEE Trans. Ind. Electron., 62, 4056–4064, <ext-link xlink:href="https://doi.org/10.1109/TIE.2015.2390194" ext-link-type="DOI">10.1109/TIE.2015.2390194</ext-link>, 2015.</mixed-citation></ref>
      <ref id="bib1.bib5"><label>5</label><?label 1?><mixed-citation>Cao, R., Jiang, N., Lu, M., Zhang, Y., and Cheng, M.: Sliding mode observer based sensorless vector control of LFSPM motor for long distance drive system, IET Electr. Power App., 13, 643–651, <ext-link xlink:href="https://doi.org/10.1049/iet-epa.2018.5505" ext-link-type="DOI">10.1049/iet-epa.2018.5505</ext-link>, 2019.</mixed-citation></ref>
      <ref id="bib1.bib6"><label>6</label><?label 1?><mixed-citation>Cheng, M., Kong, L., Zhang, B., and Wang, Z.: Sensorless vector control of complementary and modular linear flux-switching permanent magnet motor based on MRAS and SVPWM, 2015 18th International Conference on Electrical Machines and Systems (ICEMS), 25–28 October 2015, Pattaya, Thailand, IEEE, 1597–1602, <ext-link xlink:href="https://doi.org/10.1109/ICEMS.2015.7385296" ext-link-type="DOI">10.1109/ICEMS.2015.7385296</ext-link>, 2015.​​​​​​​</mixed-citation></ref>
      <ref id="bib1.bib7"><label>7</label><?label 1?><mixed-citation>Cirrincione, M., Accetta, A., Pucci, M., and Vitale, G.: MRAS Speed Observer for High-Performance Linear Induction Motor Drives Based on Linear Neural Networks, IEEE Trans. Power Electron., 28, 123–134, <ext-link xlink:href="https://doi.org/10.1109/TPEL.2012.2200506" ext-link-type="DOI">10.1109/TPEL.2012.2200506</ext-link>, 2013.</mixed-citation></ref>
      <ref id="bib1.bib8"><label>8</label><?label 1?><mixed-citation>Cupertino, F., Giangrande, P., Salvatore, L., and Pellegrino, G.: Model based design of a sensorless control scheme for permanent magnet motors using signal injection, 2010 IEEE Energy Conversion Congress and Exposition, 12–16 September 2010, Atlanta, GA, USA, IEEE, 3139–3146, <ext-link xlink:href="https://doi.org/10.1109/ECCE.2010.5618469" ext-link-type="DOI">10.1109/ECCE.2010.5618469</ext-link>, 2010.</mixed-citation></ref>
      <ref id="bib1.bib9"><label>9</label><?label 1?><mixed-citation>Du, B., Zhao, T., Han, S., Song, L., and Cui, S.: Sensorless Control Strategy for IPMSM to Reduce Audible Noise by Variable Frequency Current Injection, IEEE Trans. Ind. Electron., 67, 1149–1159, <ext-link xlink:href="https://doi.org/10.1109/TIE.2019.2898621" ext-link-type="DOI">10.1109/TIE.2019.2898621</ext-link>, 2020.</mixed-citation></ref>
      <ref id="bib1.bib10"><label>10</label><?label 1?><mixed-citation>Hussain, H. A. and Toliyat, H. A.: Back-EMF based sensorless vector control of tubular PM linear motors, 2015 IEEE International Electric Machines &amp; Drives Conference (IEMDC), 10–13 May 2015, Coeur d'Alene, ID, USA, IEEE, 878–883, <ext-link xlink:href="https://doi.org/10.1109/IEMDC.2015.7409164" ext-link-type="DOI">10.1109/IEMDC.2015.7409164</ext-link>, 2015.</mixed-citation></ref>
      <ref id="bib1.bib11"><label>11</label><?label 1?><mixed-citation>Lin, Z., Guo, H., and Xu, J.: Research on a New Sensorless Control of High Power Segmented Permanent Magnet Linear Synchronous Motor by The Improved Sliding Mode Observer, 24th International Conference on Electrical Machines and Systems (ICEMS), 31 October–3 November 2021, Gyeongju, Republic of Korea, IEEE, 862–866, <ext-link xlink:href="https://doi.org/10.23919/ICEMS52562.2021.9634460" ext-link-type="DOI">10.23919/ICEMS52562.2021.9634460</ext-link>, 2021.</mixed-citation></ref>
      <ref id="bib1.bib12"><label>12</label><?label 1?><mixed-citation>Lu, H. G., Jiang, M., Li, S. W., and Liu, B. Y.: Position Sensorless Control Based on Augmented Extended Kalman Filter for Permanent Magnet Linear Synchronous Motor, Appl. Mech. Mater., 88–89, 350–354, <ext-link xlink:href="https://doi.org/10.4028/www.scientific.net/AMM.88-89.350" ext-link-type="DOI">10.4028/www.scientific.net/AMM.88-89.350</ext-link>, 2011.</mixed-citation></ref>
      <ref id="bib1.bib13"><label>13</label><?label 1?><mixed-citation>Ofori, E., Husain, T., Sozer, Y., and Husain, I.: A Pulse-Injection-Based Sensorless Position Estimation Method for a Switched Reluctance Machine Over a Wide Speed Range, IEEE Trans. Ind. Appl., 51, 3867–3876, <ext-link xlink:href="https://doi.org/10.1109/TIA.2015.2420618" ext-link-type="DOI">10.1109/TIA.2015.2420618</ext-link>, 2015.</mixed-citation></ref>
      <ref id="bib1.bib14"><label>14</label><?label 1?><mixed-citation>Pasqualotto, D., Rigon, S., and Zigliotto, M.: Sensorless Speed Control of Synchronous Reluctance Motor Drives Based on Extended Kalman Filter and Neural Magnetic Model, IEEE Trans. Ind. Electron., 70, 1321–1330, <ext-link xlink:href="https://doi.org/10.1109/TIE.2022.3159962" ext-link-type="DOI">10.1109/TIE.2022.3159962</ext-link>, 2023.</mixed-citation></ref>
      <ref id="bib1.bib15"><label>15</label><?label 1?><mixed-citation>Pavlitov, C., Tashev, T., Dimitrov, A., Gorbounov, Y., and Rusinov, R.: Rotor position estimation for a sensorless switched reluctance motor drive with the aid of real time adaptive filter, 2014 16th International Power Electronics and Motion Control Conference and Exposition, 21–24 September 2014, Antalya, Turkey, IEEE, 77–81, <ext-link xlink:href="https://doi.org/10.1109/EPEPEMC.2014.6980576" ext-link-type="DOI">10.1109/EPEPEMC.2014.6980576</ext-link>, 2014.</mixed-citation></ref>
      <ref id="bib1.bib16"><label>16</label><?label 1?><mixed-citation>Pellegrino, G., Vagati, A., Boazzo, B., and Guglielmi, P.: Comparison of induction and PM synchronous motor drives for EV application including design examples, IEEE Trans. Ind. Appl., 48, 2322–2332, <ext-link xlink:href="https://doi.org/10.1109/TIA.2012.2227092" ext-link-type="DOI">10.1109/TIA.2012.2227092</ext-link>, 2012.</mixed-citation></ref>
      <ref id="bib1.bib17"><label>17</label><?label 1?><mixed-citation>Ting, C., Chang, Y., and Chen, Y.: Backstepping direct thrust force control for sensorless PMLSM drive, IET Electr. Power App., 13, 322–331, <ext-link xlink:href="https://doi.org/10.1049/iet-epa.2018.5269" ext-link-type="DOI">10.1049/iet-epa.2018.5269</ext-link>, 2019.</mixed-citation></ref>
      <ref id="bib1.bib18"><label>18</label><?label 1?><mixed-citation> Wang, F. and Ge, Q.: Long Stator Linear Synchronous Motor Sensorless Control Based on Hypothetical Reference Frame, Trans. of China Electrotechnical Society, 25, 37–40, 2010.</mixed-citation></ref>
      <ref id="bib1.bib19"><label>19</label><?label 1?><mixed-citation>Wang, G., Liu, R., Zhao, N., Ding, D., and Xu, D.: Enhanced Linear ADRC Strategy for HF Pulse Voltage Signal Injection-Based Sensorless IPMSM Drives, IEEE Trans. Power Electron., 34, 514–525, <ext-link xlink:href="https://doi.org/10.1109/TPEL.2018.2814056" ext-link-type="DOI">10.1109/TPEL.2018.2814056</ext-link>, 2019.</mixed-citation></ref>
      <ref id="bib1.bib20"><label>20</label><?label 1?><mixed-citation> Wang, L., Lu, Q., and Ye, Y.: Sensorless mover position estimation of permanent magnet linear synchronous motor, Electric Machines and Control, 16, 7–12, 2012.</mixed-citation></ref>
      <ref id="bib1.bib21"><label>21</label><?label 1?><mixed-citation>Xiao, D., Ye, J., Fang, G., Xia, Z., Wang, X., and Emadi, A.: Improved Feature-Position-Based Sensorless Control Scheme for SRM Drives Based on Nonlinear State Observer at Medium and High Speeds, IEEE Trans. Power Electron., 36, 5711–5723, <ext-link xlink:href="https://doi.org/10.1109/TPEL.2020.3030007" ext-link-type="DOI">10.1109/TPEL.2020.3030007</ext-link>, 2021.</mixed-citation></ref>
      <ref id="bib1.bib22"><label>22</label><?label 1?><mixed-citation>Xiao, G., Tu, W., Suo, C., Tang, L., and Yang, K: Research and design of speed control for high speed sensorless brushless DC motor with commutation compensation, 2017 20th International Conference on Electrical Machines and Systems (ICEMS), 11–14 August 2017, Sydney, NSW, Australia, IEEE, 1–5, <ext-link xlink:href="https://doi.org/10.1109/ICEMS.2017.8056026" ext-link-type="DOI">10.1109/ICEMS.2017.8056026</ext-link>, 2017.</mixed-citation></ref>
      <ref id="bib1.bib23"><label>23</label><?label 1?><mixed-citation>Yang, C., Song, B., Xie, Y., Lu, S., and Tang, X.: Stable Simultaneous Inertia and Disturbance Torque Identification for SPMSM Drive Systems Subject to Mismatched Rotor Flux Linkage, IEEE J. Em. Sel. Top. P., 10, 2445–2462, <ext-link xlink:href="https://doi.org/10.1109/JESTPE.2021.3111122" ext-link-type="DOI">10.1109/JESTPE.2021.3111122</ext-link>, 2022.</mixed-citation></ref>
      <ref id="bib1.bib24"><label>24</label><?label 1?><mixed-citation>Zhang, T., Xu, Z., and Gerada, C.: A Nonlinear Extended State Observer for Sensorless IPMSM Drives With Optimized Gains, IEEE Trans. Ind. Appl., 56, 1485–1494, <ext-link xlink:href="https://doi.org/10.1109/TIA.2019.2959537" ext-link-type="DOI">10.1109/TIA.2019.2959537</ext-link>, 2020.</mixed-citation></ref>
      <?pagebreak page109?><ref id="bib1.bib25"><label>25</label><?label 1?><mixed-citation>Zhang, X., Hou, B., and Yang M.: Deadbeat predictive current control of permanent-magnet synchronous motors with stator current and disturbance observer, IEEE Trans. Power Electron., 32, 3818–3834, <ext-link xlink:href="https://doi.org/10.1109/TPEL.2016.2592534" ext-link-type="DOI">10.1109/TPEL.2016.2592534</ext-link>, 2016.</mixed-citation></ref>
      <ref id="bib1.bib26"><label>26</label><?label 1?><mixed-citation>Zhang, Y., Yin, Z., Bai, C., Wang, G., and Liu, J.: A Rotor Position and Speed Estimation Method Using an Improved Linear Extended State Observer for IPMSM Sensorless Drives, IEEE Trans. Power Electron., 36, 14062–14073, <ext-link xlink:href="https://doi.org/10.1109/TPEL.2021.3085126" ext-link-type="DOI">10.1109/TPEL.2021.3085126</ext-link>, 2021. </mixed-citation></ref><?xmltex \hack{\newpage}?>
      <ref id="bib1.bib27"><label>27</label><?label 1?><mixed-citation>Zhao, W., Yang, A., Ji, J., Chen, Q., and Zhu, J.: Modified Flux Linkage Observer for Sensorless Direct Thrust Force Control of Linear Vernier Permanent Magnet Motor, IEEE Trans. Power Electron., 34, 7800–7811, <ext-link xlink:href="https://doi.org/10.1109/TPEL.2018.2879411" ext-link-type="DOI">10.1109/TPEL.2018.2879411</ext-link>, 2019.</mixed-citation></ref>

  </ref-list></back>
    <!--<article-title-html>Improved flux linkage observer for position estimation of permanent magnet synchronous linear motor</article-title-html>
<abstract-html/>
<ref-html id="bib1.bib1"><label>1</label><mixed-citation>
      
Accetta, A., Cirrincione, M., Pucci, M., and Vitale, G.: Neural Sensorless Control of Linear Induction Motors by a Full-Order Luenberger Observer Considering the End Effects, IEEE Trans. Ind. Appl., 50, 1891–1904, <a href="https://doi.org/10.1109/TIA.2013.2288429" target="_blank">https://doi.org/10.1109/TIA.2013.2288429</a>, 2014.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib2"><label>2</label><mixed-citation>
      
Afifa, R., Ali, S., Pervaiz, M., and Lqbal, J.: Adaptive backstepping
integral sliding mode control of a mimo separately excited DC motor,  Robotics, 12, 105, <a href="https://doi.org/10.3390/robotics12040105" target="_blank">https://doi.org/10.3390/robotics12040105</a>, 2023.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib3"><label>3</label><mixed-citation>
      
Ahmad, S., Uppal, A., Azam, M., and Lqbal, J.: Chattering free sliding mode
control and state dependent Kalman filter design for underground
gasification energy conversion process, Electronics, 12, 876, <a href="https://doi.org/10.3390/electronics12040876" target="_blank">https://doi.org/10.3390/electronics12040876</a>, 2023.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib4"><label>4</label><mixed-citation>
      
Cao, R., Cheng, M., and Zhang, B.: Speed control of complementary and
modular linear flux-switching permanent-magnet motor, IEEE Trans. Ind.
Electron., 62, 4056–4064, <a href="https://doi.org/10.1109/TIE.2015.2390194" target="_blank">https://doi.org/10.1109/TIE.2015.2390194</a>, 2015.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib5"><label>5</label><mixed-citation>
      
Cao, R., Jiang, N., Lu, M., Zhang, Y., and Cheng, M.: Sliding mode observer
based sensorless vector control of LFSPM motor for long distance drive
system, IET Electr. Power App., 13, 643–651,
<a href="https://doi.org/10.1049/iet-epa.2018.5505" target="_blank">https://doi.org/10.1049/iet-epa.2018.5505</a>, 2019.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib6"><label>6</label><mixed-citation>
      
Cheng, M., Kong, L., Zhang, B., and Wang, Z.: Sensorless vector control of
complementary and modular linear flux-switching permanent magnet motor based
on MRAS and SVPWM, 2015 18th International Conference on Electrical Machines
and Systems (ICEMS), 25–28 October 2015, Pattaya, Thailand, IEEE, 1597–1602, <a href="https://doi.org/10.1109/ICEMS.2015.7385296" target="_blank">https://doi.org/10.1109/ICEMS.2015.7385296</a>, 2015.​​​​​​​

    </mixed-citation></ref-html>
<ref-html id="bib1.bib7"><label>7</label><mixed-citation>
      
Cirrincione, M., Accetta, A., Pucci, M., and Vitale, G.: MRAS Speed Observer
for High-Performance Linear Induction Motor Drives Based on Linear Neural
Networks, IEEE Trans. Power Electron., 28, 123–134,
<a href="https://doi.org/10.1109/TPEL.2012.2200506" target="_blank">https://doi.org/10.1109/TPEL.2012.2200506</a>, 2013.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib8"><label>8</label><mixed-citation>
      
Cupertino, F., Giangrande, P., Salvatore, L., and Pellegrino, G.: Model based
design of a sensorless control scheme for permanent magnet motors using
signal injection, 2010 IEEE Energy Conversion Congress and Exposition, 12–16 September 2010, Atlanta, GA, USA, IEEE, 3139–3146, <a href="https://doi.org/10.1109/ECCE.2010.5618469" target="_blank">https://doi.org/10.1109/ECCE.2010.5618469</a>, 2010.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib9"><label>9</label><mixed-citation>
      
Du, B., Zhao, T., Han, S., Song, L., and Cui, S.: Sensorless Control Strategy
for IPMSM to Reduce Audible Noise by Variable Frequency Current Injection,
IEEE Trans. Ind. Electron., 67, 1149–1159,
<a href="https://doi.org/10.1109/TIE.2019.2898621" target="_blank">https://doi.org/10.1109/TIE.2019.2898621</a>, 2020.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib10"><label>10</label><mixed-citation>
      
Hussain, H. A. and Toliyat, H. A.: Back-EMF based sensorless vector control
of tubular PM linear motors, 2015 IEEE International Electric Machines &amp;
Drives Conference (IEMDC), 10–13 May 2015, Coeur d'Alene, ID, USA, IEEE, 878–883, <a href="https://doi.org/10.1109/IEMDC.2015.7409164" target="_blank">https://doi.org/10.1109/IEMDC.2015.7409164</a>, 2015.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib11"><label>11</label><mixed-citation>
      
Lin, Z., Guo, H., and Xu, J.: Research on a New Sensorless Control of High Power Segmented Permanent Magnet Linear Synchronous Motor by The Improved Sliding Mode Observer, 24th International Conference on Electrical Machines and Systems (ICEMS), 31 October–3 November 2021, Gyeongju, Republic of Korea, IEEE, 862–866, <a href="https://doi.org/10.23919/ICEMS52562.2021.9634460" target="_blank">https://doi.org/10.23919/ICEMS52562.2021.9634460</a>, 2021.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib12"><label>12</label><mixed-citation>
      
Lu, H. G., Jiang, M., Li, S. W., and Liu, B. Y.: Position Sensorless Control Based on Augmented Extended Kalman Filter for Permanent Magnet Linear Synchronous Motor, Appl. Mech. Mater., 88–89, 350–354, <a href="https://doi.org/10.4028/www.scientific.net/AMM.88-89.350" target="_blank">https://doi.org/10.4028/www.scientific.net/AMM.88-89.350</a>, 2011.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib13"><label>13</label><mixed-citation>
      
Ofori, E., Husain, T., Sozer, Y., and Husain, I.: A Pulse-Injection-Based
Sensorless Position Estimation Method for a Switched Reluctance Machine Over
a Wide Speed Range, IEEE Trans. Ind. Appl., 51, 3867–3876,
<a href="https://doi.org/10.1109/TIA.2015.2420618" target="_blank">https://doi.org/10.1109/TIA.2015.2420618</a>, 2015.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib14"><label>14</label><mixed-citation>
      
Pasqualotto, D., Rigon, S., and Zigliotto, M.: Sensorless Speed Control of
Synchronous Reluctance Motor Drives Based on Extended Kalman Filter and
Neural Magnetic Model, IEEE Trans. Ind. Electron., 70, 1321–1330,
<a href="https://doi.org/10.1109/TIE.2022.3159962" target="_blank">https://doi.org/10.1109/TIE.2022.3159962</a>, 2023.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib15"><label>15</label><mixed-citation>
      
Pavlitov, C., Tashev, T., Dimitrov, A., Gorbounov, Y., and Rusinov, R.: Rotor
position estimation for a sensorless switched reluctance motor drive with
the aid of real time adaptive filter, 2014 16th International Power
Electronics and Motion Control Conference and Exposition, 21–24 September
2014, Antalya, Turkey, IEEE, 77–81, <a href="https://doi.org/10.1109/EPEPEMC.2014.6980576" target="_blank">https://doi.org/10.1109/EPEPEMC.2014.6980576</a>, 2014.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib16"><label>16</label><mixed-citation>
      
Pellegrino, G., Vagati, A., Boazzo, B., and Guglielmi, P.: Comparison of
induction and PM synchronous motor drives for EV application including
design examples, IEEE Trans. Ind. Appl., 48, 2322–2332,
<a href="https://doi.org/10.1109/TIA.2012.2227092" target="_blank">https://doi.org/10.1109/TIA.2012.2227092</a>, 2012.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib17"><label>17</label><mixed-citation>
      
Ting, C., Chang, Y., and Chen, Y.: Backstepping direct thrust force control for sensorless PMLSM drive, IET Electr. Power App., 13, 322–331,
<a href="https://doi.org/10.1049/iet-epa.2018.5269" target="_blank">https://doi.org/10.1049/iet-epa.2018.5269</a>, 2019.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib18"><label>18</label><mixed-citation>
      
Wang, F. and Ge, Q.: Long Stator Linear Synchronous Motor Sensorless Control
Based on Hypothetical Reference Frame, Trans. of China Electrotechnical
Society, 25, 37–40, 2010.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib19"><label>19</label><mixed-citation>
      
Wang, G., Liu, R., Zhao, N., Ding, D., and Xu, D.: Enhanced Linear ADRC
Strategy for HF Pulse Voltage Signal Injection-Based Sensorless IPMSM
Drives, IEEE Trans. Power Electron., 34, 514–525,
<a href="https://doi.org/10.1109/TPEL.2018.2814056" target="_blank">https://doi.org/10.1109/TPEL.2018.2814056</a>, 2019.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib20"><label>20</label><mixed-citation>
      
Wang, L., Lu, Q., and Ye, Y.: Sensorless mover position estimation of permanent magnet linear synchronous motor, Electric Machines and Control, 16, 7–12, 2012.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib21"><label>21</label><mixed-citation>
      
Xiao, D., Ye, J., Fang, G., Xia, Z., Wang, X., and Emadi, A.: Improved
Feature-Position-Based Sensorless Control Scheme for SRM Drives Based on
Nonlinear State Observer at Medium and High Speeds, IEEE Trans. Power
Electron., 36, 5711–5723, <a href="https://doi.org/10.1109/TPEL.2020.3030007" target="_blank">https://doi.org/10.1109/TPEL.2020.3030007</a>, 2021.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib22"><label>22</label><mixed-citation>
      
Xiao, G., Tu, W., Suo, C., Tang, L., and Yang, K: Research and design of
speed control for high speed sensorless brushless DC motor with commutation
compensation, 2017 20th International Conference on Electrical Machines and
Systems (ICEMS), 11–14 August 2017, Sydney, NSW, Australia, IEEE, 1–5,
<a href="https://doi.org/10.1109/ICEMS.2017.8056026" target="_blank">https://doi.org/10.1109/ICEMS.2017.8056026</a>, 2017.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib23"><label>23</label><mixed-citation>
      
Yang, C., Song, B., Xie, Y., Lu, S., and Tang, X.: Stable Simultaneous
Inertia and Disturbance Torque Identification for SPMSM Drive Systems
Subject to Mismatched Rotor Flux Linkage, IEEE J. Em. Sel. Top. P., 10, 2445–2462, <a href="https://doi.org/10.1109/JESTPE.2021.3111122" target="_blank">https://doi.org/10.1109/JESTPE.2021.3111122</a>, 2022.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib24"><label>24</label><mixed-citation>
      
Zhang, T., Xu, Z., and Gerada, C.: A Nonlinear Extended State Observer for
Sensorless IPMSM Drives With Optimized Gains, IEEE Trans. Ind. Appl., 56,
1485–1494, <a href="https://doi.org/10.1109/TIA.2019.2959537" target="_blank">https://doi.org/10.1109/TIA.2019.2959537</a>, 2020.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib25"><label>25</label><mixed-citation>
      
Zhang, X., Hou, B., and Yang M.: Deadbeat predictive current control of
permanent-magnet synchronous motors with stator current and disturbance
observer, IEEE Trans. Power Electron., 32, 3818–3834,
<a href="https://doi.org/10.1109/TPEL.2016.2592534" target="_blank">https://doi.org/10.1109/TPEL.2016.2592534</a>, 2016.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib26"><label>26</label><mixed-citation>
      
Zhang, Y., Yin, Z., Bai, C., Wang, G., and Liu, J.: A Rotor Position and
Speed Estimation Method Using an Improved Linear Extended State Observer for
IPMSM Sensorless Drives, IEEE Trans. Power Electron., 36, 14062–14073,
<a href="https://doi.org/10.1109/TPEL.2021.3085126" target="_blank">https://doi.org/10.1109/TPEL.2021.3085126</a>, 2021.


    </mixed-citation></ref-html>
<ref-html id="bib1.bib27"><label>27</label><mixed-citation>
      
Zhao, W., Yang, A., Ji, J., Chen, Q., and Zhu, J.: Modified Flux Linkage
Observer for Sensorless Direct Thrust Force Control of Linear Vernier
Permanent Magnet Motor, IEEE Trans. Power Electron., 34, 7800–7811,
<a href="https://doi.org/10.1109/TPEL.2018.2879411" target="_blank">https://doi.org/10.1109/TPEL.2018.2879411</a>, 2019.

    </mixed-citation></ref-html>--></article>
