# Patent

## Title
Rotary electric machine controller, rotary electric machine control method, and method of creating control map

## Abstract

A rotary electric machine controller according to the present application includes: a superimposer for superimposing a high-frequency signal on one of a voltage and a current of a motor having saliency; an estimator for calculating an estimate value of a magnetic-pole position of a rotor included in the rotary electric machine based on a high-frequency component appearing in the other one of the voltage and the current of the rotary electric machine; a retainer for retaining information relating to a correction amount corresponding to the estimate value of the magnetic-pole position of the rotor; and a corrector for correcting the estimate value of the magnetic-pole position of the rotor based on the information relating to the correction amount.

## Background
INCORPORATION BY REFERENCE
The present disclosure contains subject matter related to that disclosed in Japanese Priority Patent Application JP2014-006548 filed in the Japan Patent Office on Jan. 17, 2014, the entire contents of which are hereby incorporated by reference.
BACKGROUND OF THE INVENTION
1. Field of the Invention
The present application relates to a rotary electric machine controller, a rotary electric machine control method, and a method of creating a control map.
2. Description of the Related Art
In Japanese Patent Application Laid-open No. 2013-090552, the following rotary electric machine controller is disclosed. The rotary electric machine controller applies a high-frequency current to a rotary electric machine including a rotor having saliency and estimates a magnetic-pole direction of the rotor based on a high-frequency component contained in a voltage command as a component in response to the high-frequency current. The rotary electric machine controller calculates an error of an estimate value of the magnetic-pole direction, which is generated due to a d-q axis magnetic flux interference.
The inventors of the present application have found that the error of the estimate value of the magnetic-pole direction, which is generated due to the d-q axis magnetic flux interference, is not only a steady error as described in Japanese Patent Application Laid-open No. 2013-090552 but may also be generated depending on a magnetic-pole position of the rotor.
The present application has been made to solve the problem described above, and has an object to provide a rotary electric machine controller, a rotary electric machine control method, and a method of creating a control map, which are capable of suppressing an error of an estimate value corresponding to a magnetic-pole position of a rotor.
SUMMARY OF THE INVENTION
In order to solve the above-mentioned problem, according to one embodiment of the present application, there is provided a rotary electric machine controller including: a superimposer for superimposing a high-frequency signal on one of a voltage and a current of a rotary electric machine having saliency; an estimator for calculating an estimate value of a magnetic-pole position of a rotor included in the rotary electric machine based on a high-frequency component appearing in the other one of the voltage and the current of the rotary electric machine; a retainer for retaining information relating to a correction amount corresponding to the magnetic-pole position of the rotor; and a corrector for correcting the estimate value of the magnetic-pole position of the rotor based on the information relating to the correction amount.
Further, according to one embodiment of the present application, there is provided a rotary electric machine control method including: superimposing a high-frequency signal on one of a voltage and a current of a rotary electric machine having saliency; calculating an estimate value of a magnetic-pole position of a rotor included in the rotary electric machine based on a high-frequency component appearing in the other one of the voltage and the current of the rotary electric machine; retaining information relating to a correction amount corresponding to the magnetic-pole position of the rotor; and correcting the estimate value of the magnetic-pole position of the rotor based on the information relating to the correction amount.
Further, according to one embodiment of the present application, there is provided a method of creating a control map to be retained in a rotary electric machine controller, the method including: generating, for each of at least two points of a magnetic-pole position of a rotor included in a rotary electric machine having saliency, information relating to a correction amount for correcting an estimate value of the magnetic-pole position of the rotor, the correction amount being calculated based on a high-frequency component appearing in one of a voltage and a current of the rotary electric machine when a high-frequency signal is superimposed on the other one of the voltage and the current, the correction amount corresponding to the magnetic-pole position of the rotor; and associating, on the control map, the points of the magnetic-pole position of the rotor and the information relating to the correction amount generated for each of the points with each other.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a diagram illustrating a configuration example of a motor controller according to an embodiment of the present application.
FIG. 2 is a graph showing a coordinate system.
FIG. 3 is a diagram illustrating a configuration example of a magnetic-pole position corrector.
FIG. 4 is a diagram illustrating a content example of a control map.
FIG. 5 is a diagram illustrating a relationship between high-frequency inductances.
FIG. 6 is a flowchart illustrating an operation example of the motor controller.
FIG. 7A is a graph showing a relationship between Δθ and λ corresponding to a magnetic-pole position when Id=0 [A] and Iq=0 [A].
FIG. 7B is a graph showing a relationship between Δθ and λ corresponding to the magnetic-pole position when Id=2 [A] and Iq=0 [A].
FIG. 7C is a graph showing a relationship between Δθ and λ corresponding to the magnetic-pole position when Id=0 [A] and Iq=7 [A].
FIG. 7D is a graph showing a relationship between Δθ and λ corresponding to the magnetic-pole position when Id=2 [A] and Iq=7 [A].
FIG. 8 is a diagram illustrating a method of creating a control map with an actual measurement.



DESCRIPTION OF THE EMBODIMENTS
An embodiment of the present application is described referring to the drawings.
FIG. 1 is a diagram illustrating a configuration example of a motor controller 1 according to an embodiment of the present application. The motor controller 1 is an example of a rotary electric machine controller. Each of blocks included in the motor controller 1 is realized by, for example, a microprocessor executing information processing in accordance with a program stored in a memory.
A motor 9 to be controlled by the motor controller 1 is an example of a rotary electric machine, and is, for example, a synchronous motor. The motor 9 is a sensorless motor without a detector for detecting a position of a rotor. The motor controller 1 uses saliency of the motor 9 to estimate a magnetic-pole position of the rotor.
The motor controller 1 includes a subtractor 11, a position controller 13, a subtractor 15, a speed controller 17, subtractors 19d and 19q, a current controller 21, an adder 23, a PWM controller 25, a three-to-two phase converter 31, a d-q converter 33, and a notch filter 35 as blocks for controlling the rotation of the motor 9.
Further, the motor controller 1 includes a high-frequency superimposer 41, an envelope extractor 43, a magnetic-pole position estimator 45, a magnetic-pole position corrector 51, a bandpass filter 53, a map retainer 55, and a position/speed calculator 57 as blocks for estimating a magnetic-pole position of the rotor of the motor 9.
The subtractor 11 calculates a differential value between a position command value P* output from, for example, a high-order system (not shown) and a position estimate value P^ output from the position/speed calculator 57, and then outputs the calculated differential value to the position controller 13.
The position controller 13 calculates a speed command value ω* for suppressing the differential value output from the subtractor 11, and then outputs the calculated speed command value ω* to the subtractor 15.
The subtractor 15 calculates a differential value between the speed command value ω* output from the position controller 13 and a speed estimate value ω^ output from the position/speed calculator 57, and then outputs the calculated differential value to the speed controller 17.
The speed controller 17 calculates a q-axis current command value Iq* for suppressing the differential value output from the subtractor 15, and then outputs the calculated q-axis current command value Iq* to the subtractor 19q. 
The subtractor 19d calculates a differential value between a d-axis current command value Id* output from, for example, a high-order system and a d-axis current detection value Id′ output from the notch filter 35, and then outputs the calculated differential value to the current controller 21.
The subtractor 19q calculates a differential value between the q-axis current command value Iq* output from the speed controller 17 and a q-axis current detection value Iq′ output from the notch filter 35, and then outputs the calculated differential value to the current controller 21.
The current controller 21 calculates a d-axis voltage command value Vd* and a q-axis voltage command value Vq* for suppressing the differential values respectively output from the subtractors 19d and 19q. Then, the current controller 21 outputs the d-axis voltage command value Vd* to the adder 23 and the q-axis voltage command value Vq* to the PWM controller 25.
The adder 23 adds the d-axis voltage command value Vd* output from the current controller 21 to a high-frequency signal Vsig output from the high-frequency superimposer 41, and then outputs a d-axis voltage total value obtained by the addition to the PWM controller 25.
The PWM controller 25 generates a PWM signal for driving the motor 9 based on the d-axis voltage total value output from the adder 23 and the q-axis voltage command value Vq* output from the current controller 21, and then outputs the thus generated PWM signal to the motor 9. When generating the PWM signal, the PWM controller 25 uses a correction value θa^ of the magnetic-pole position (magnetic-pole phase) output from the magnetic-pole position corrector 51.
The three-to-two phase converter 31 converts current detection values Iu, Iv, and Iw of three-phase currents flowing through the motor 9 into current detection values Iα and Iβ of a fixed α-β coordinate system, and then outputs the thus obtained current detection values Iα and Iβ to the d-q converter 33 and the envelope extractor 43. In the fixed α-β coordinate system, the α-axis and the β-axis are orthogonal to each other, and the α-axis is fixed in a direction of a U-phase, for example.
The d-q converter 33 converts the current detection values Iα and Iβ of the fixed α-β coordinate system, which are output from the three-to-two phase converter 31, into current detection values Id and Iq of a rotating d-q coordinate system, and then outputs the thus obtained current detection values Id and Iq to the notch filter 35, the magnetic-pole position corrector 51, and the bandpass filter 53. When performing the d-q conversion, the d-q converter 33 uses the correction value θa^ of the magnetic-pole position output from the magnetic-pole position corrector 51.
In the rotating d-q coordinate system, the d-axis is a magnetic-pole direction of the rotor of the motor 9, and the q-axis is an axis orthogonal to the d-axis. The d-axis current Id is also called “field current”, whereas the q-axis current Iq is also called “torque current”. The details of the rotating d-q coordinate system are described later.
The notch filter (band-stop filter) 35 suppresses high-frequency components of the current detection values Id and Iq output from the d-q converter 33, and then outputs the current detection values Id′ and Iq′ whose high-frequency components are suppressed to the subtractors 19d and 19q. 
Now, the blocks for estimating the magnetic-pole position of the rotor of the motor 9 are described.
The high-frequency superimposer 41 is an example of a superimposer. The high-frequency superimposer 41 outputs the high-frequency signal Vsig to the adder 23 so as to superimpose a high-frequency signal on the d-axis voltage command value Vd* output from the current controller 21. A frequency of the high-frequency signal Vsig is sufficiently higher than a fundamental frequency of the voltage applied to the motor 9.
The mode of superimposition is not limited to superimposing the high-frequency signal only on the d-axis voltage command value Vd*. The high-frequency signal may be superimposed on the q-axis voltage command value Vq*, and may also be superimposed on both the d-axis voltage command value Vd* and the q-axis voltage command value Vq*.
Alternatively, the high-frequency signal may be superimposed on one or both of the d-axis current command value Id* and the q-axis current command value Iq*. In this case, the magnetic-pole position can be estimated by using a response component which appears in the voltage applied to the motor 9.
The envelope extractor 43 extracts envelopes of the current detection values Iα and Iβ output from the three-to-two phase converter 31 based on the high-frequency signal Vsig output from the high-frequency superimposer 41, and then outputs envelope detection values Isin and Icos to the magnetic-pole position estimator 45.
The magnetic-pole position estimator 45 is an example of an estimator. The magnetic-pole position estimator 45 calculates an estimate value θr^ of the magnetic-pole position based on the envelope detection values Isin and Icos output from the envelope extractor 43, and then outputs the thus calculated estimate value θr^ to the magnetic-pole position corrector 51.
A method of calculating the estimate value θr^ of the magnetic-pole position is not limited to the above-mentioned mode in which the high-frequency signal is superimposed on the d-axis voltage command value are the envelopes are extracted. For example, there is a method of estimating the magnetic-pole position based on a parallel component and an orthogonal component of the motor current in response to an alternating voltage (or alternating current) applied to the motor. Also, there is a method of estimating the magnetic-pole position at a low speed including that in a stopped state by using a carrier frequency component. Further, there is a method of estimating the magnetic-pole position by detecting the motor current while changing an application voltage in each predetermined cycle of a carrier wave and obtaining a vector of the amount of change in current for each predetermined cycle. A method of calculating the estimate value θr^ of the magnetic-pole position is publicly known, and therefore the detailed description thereof is herein omitted.
The rotating d-q coordinate system estimated by the magnetic-pole position estimator 45 is now described referring to FIG. 2. The magnetic-pole position estimator 45 calculates the estimate value θr^ of the magnetic-pole position, and thereby an estimated d-axis and an estimated q-axis are defined. The estimate value θr^ of the magnetic-pole position is given as, for example, an angle between the estimated d-axis and the α-axis. A symbol θre shown in FIG. 2 denotes the estimate value θr^ of the magnetic-pole position (both “e” and “^” indicate that it is an estimate value).
However, the estimated d-axis and the estimated q-axis defined by the magnetic-pole position estimator 45 may be shifted from actual d-axis and q-axis of the motor 9 in some cases. A symbol θr shown in FIG. 2 indicates an actual magnetic-pole position given as the angle between the actual d-axis and the α-axis. A symbol Δθ shown in FIG. 2 indicates an error (axial shift) between the estimated d-axis and the actual d-axis. Such an error Δθ is generated by d-q axis magnetic flux interference.
In Japanese Patent Application Laid-open No. 2013-090552, it is described that a steady error generated in accordance with a magnitude of the current flowing through the motor is corrected.
The inventors of the present application have found that the error Δθ of the estimate value θr^ of the magnetic-pole position contains not only the steady component as described in Japanese Patent Application Laid-open No. 2013-090552 but also a component depending on the magnetic-pole position of the rotor.
Therefore, in this embodiment, the error Δθ contained in the estimate value θr^ of the magnetic-pole position is suppressed by the operation of the magnetic-pole position corrector 51 described below.
Returning to the description referring to FIG. 1, the magnetic-pole position corrector 51 is an example of a corrector. Based on a control map retained by the map retainer 55, the magnetic-pole position corrector 51 corrects the estimate value θr^ of the magnetic-pole position, which is output from the magnetic-pole position estimator 45, and then outputs a correction value θa^ of the magnetic-pole position obtained by the correction to the position/speed calculator 57. The details of the magnetic-pole position corrector 51 are described later.
The bandpass filter 53 allows the high-frequency components of the d-axis current Id and the q-axis current Iq output from the d-q converter 33 to pass therethrough and suppresses the other frequency components. Specifically, the bandpass filter 53 outputs a high-frequency component of the d-axis current Id (d-axis current high-frequency component Idh) and a high-frequency component of the q-axis current Iq (q-axis current high-frequency component Iqh) to the magnetic-pole position corrector 51.
The map retainer 55 is an example of a retainer. The map retainer 55 retains the control map for correcting the estimate value θr^ of the magnetic-pole position, and returns a parameter λ relating to a correction amount in response to a request from the magnetic-pole position corrector 51 to the magnetic-pole position corrector 51. The details of the control map are described later.
The position/speed calculator 57 calculates the position estimate value P^ and the speed estimate value ω^ based on the correction value θa^ of the magnetic-pole position, which is output from the magnetic-pole position corrector 51. Then, the position/speed calculator 57 outputs the position estimate value P^ to the subtractor 11 and the speed estimate value ω^ to the subtractor 15.
FIG. 3 is a diagram illustrating a configuration example of the magnetic-pole position corrector 51. The magnetic-pole position corrector 51 includes a λ-acquiring unit 61, a multiplier 62, an adder 63, an integrator 65, a PI controller 67, and an adder 69.
The estimate value θr^ of the magnetic-pole position output from the magnetic-pole position estimator 45, the d-axis current Id and the q-axis current Iq output from the d-q converter 33, and the d-axis current high-frequency component Idh and the q-axis current high-frequency component Iqh output from the bandpass filter 53 are input to the magnetic-pole position corrector 51. The d-axis current and the q-axis current may be input from the notch filter 35 also.
The λ-acquiring unit 61 acquires the parameter λ relating to the correction amount corresponding to the estimate value θr^ of the magnetic-pole position, the d-axis current Id, and the q-axis current Iq, from the control map. The control map is a four-dimensional map in which the magnetic-pole position θ, the d-axis current Id, and the q-axis current Iq, and the parameter λ relating to the correction amount are associated with each other (see FIG. 4).
The acquisition of the parameter λ is not limited to the mode in which the control map is retained. Information of mathematical formulae for calculating the parameter λ relating to the correction amount from the magnetic-pole position θ, the d-axis current Id, and the q-axis current Iq may be retained instead.
The multiplier 62 multiplies the parameter λ output from the λ-acquiring unit 61 by the d-axis current high-frequency component Idh, and then outputs λIdh obtained by the multiplication to the adder 63.
The adder 63 adds λIdh output from the multiplier 62 to the q-axis current high-frequency component Iqh, and then outputs Iqh+λIdh obtained by the addition to the integrator 65.
The integrator 65 integrates Iqh+λIdh output from the adder 63, and then outputs the thus obtained integrated value of Iqh+λIdh to the PI controller 67. The integrator 65 is a block provided to stabilize the control, and is not indispensable.
The PI controller 67 performs PI control so that the integrated value of Iqh+λIdh, which is output from the integrator 65, approaches zero, and then outputs a correction amount obtained by the PI control to the adder 69. The control is not limited to the PI control, and may also be control for increasing or decreasing a predetermined value in accordance with a positive/negative sign of the integrated value.
The adder 69 adds the estimate value θr^ of the magnetic-pole position to the correction amount output from the PI controller 67, and then outputs the correction value θa^ of the magnetic-pole position obtained by the addition.
The correction amount output from the PI controller 67 is not completely identical with the error Δθ (see FIG. 2) contained in the estimate value θr^ of the magnetic-pole position as described later. However, because Iqh+λIdh is made to approach zero, the correction amount is approximated to the error Δθ to the extent that the correction amount and the error Δθ become approximately the same.
The reason why Iqh+λIdh is made to approach zero and the details of the parameter λ relating to the correction amount are described with mathematical formulae.
In the following description of the mathematical formulae, the superscript “e” follows Vd, Vq, Id, Iq, Idh, Iqh, and θr relating to the estimated d- and q-axes, clearly marking them as different from those relating to actual d- and q-axes.
A voltage equation on the d- and q-axes in consideration of the d- and q-axes mutual interference is expressed by Mathematical Formula 2.




  
    
      
        
          [
          
            
              
                
                  V
                  d
                
              
            
            
              
                
                  V
                  q
                
              
            
          
          ]
        
        =
        
          
            
              [
              
                
                  
                    
                      
                        R
                        s
                      
                      +
                      
                        
                          L
                          dh
                        
                        ⁢
                        p
                      
                    
                  
                  
                    
                      
                        
                          -
                          
                            ω
                            r
                          
                        
                        ⁢
                        
                          L
                          q
                        
                      
                      +
                      
                        
                          L
                          dqh
                        
                        ⁢
                        p
                      
                    
                  
                
                
                  
                    
                      
                        
                          ω
                          r
                        
                        ⁢
                        
                          L
                          d
                        
                      
                      +
                      
                        
                          L
                          qdh
                        
                        ⁢
                        p
                      
                    
                  
                  
                    
                      
                        R
                        s
                      
                      +
                      
                        
                          L
                          qh
                        
                        ⁢
                        p
                      
                    
                  
                
              
              ]
            
            ⁡
            
              [
              
                
                  
                    
                      I
                      d
                    
                  
                
                
                  
                    
                      I
                      q
                    
                  
                
              
              ]
            
          
          +
          
            [
            
              
                
                  0
                
              
              
                
                  
                    
                      ω
                      r
                    
                    ⁢
                    
                      Φ
                      m
                    
                  
                
              
            
            ]
          
        
      
    
    
      
        [
        
          Math
          .
          
            
          
          ⁢
          2
        
        ]
      
    
  




In Mathematical Formula 2, Rs is a phase resistance of a winding of a stator, Φm is an interlinkage magnetic flux for one phase, which is generated by a permanent magnet, ωr is a rotation speed, p is a differential operator, Vd is the d-axis voltage, Vq is the q-axis voltage, Id is the d-axis current, Iq is the q-axis current, Ld is a d-axis self-inductance, and Lq is a q-axis self-inductance.
In addition, Ldh is a d-axis high-frequency inductance, Lqh is a q-axis high-frequency inductance, Ldqh is a high-frequency inductance due to interference of the d-axis with the q-axis, and Lqdh is a high-frequency inductance due to interference of the q-axis with the d-axis. The inductances Ldh, Lqh, Ldqh, and Lqdh are expressed by Mathematical Formulae 3.




  
    
      
        
          
            L
            dh
          
          =
          
            
              
                
                  Ψ
                  d
                
                ⁡
                
                  (
                  
                    
                      
                        I
                        d
                      
                      +
                      
                        Δ
                        ⁢
                        
                          
                        
                        ⁢
                        
                          I
                          d
                        
                      
                    
                    ,
                    
                      I
                      q
                    
                    ,
                    
                      Φ
                      m
                    
                  
                  )
                
              
              -
              
                
                  Ψ
                  d
                
                ⁡
                
                  (
                  
                    
                      I
                      d
                    
                    ,
                    
                      I
                      q
                    
                    ,
                    
                      Φ
                      m
                    
                  
                  )
                
              
            
            
              Δ
              ⁢
              
                
              
              ⁢
              
                I
                d
              
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            L
            qh
          
          =
          
            
              
                
                  Ψ
                  q
                
                ⁡
                
                  (
                  
                    
                      I
                      d
                    
                    ,
                    
                      
                        I
                        q
                      
                      +
                      
                        Δ
                        ⁢
                        
                          
                        
                        ⁢
                        
                          I
                          q
                        
                      
                    
                    ,
                    
                      Φ
                      m
                    
                  
                  )
                
              
              -
              
                
                  Ψ
                  q
                
                ⁡
                
                  (
                  
                    
                      I
                      d
                    
                    ,
                    
                      I
                      q
                    
                    ,
                    
                      Φ
                      m
                    
                  
                  )
                
              
            
            
              Δ
              ⁢
              
                
              
              ⁢
              
                I
                q
              
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            L
            dqh
          
          =
          
            
              
                
                  Ψ
                  d
                
                ⁡
                
                  (
                  
                    
                      I
                      d
                    
                    ,
                    
                      
                        I
                        q
                      
                      +
                      
                        Δ
                        ⁢
                        
                          
                        
                        ⁢
                        
                          i
                          q
                        
                      
                    
                    ,
                    
                      Φ
                      m
                    
                  
                  )
                
              
              -
              
                
                  Ψ
                  d
                
                ⁡
                
                  (
                  
                    
                      I
                      d
                    
                    ,
                    
                      I
                      q
                    
                    ,
                    
                      Φ
                      m
                    
                  
                  )
                
              
            
            
              Δ
              ⁢
              
                
              
              ⁢
              
                I
                q
              
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            L
            qdh
          
          =
          
            
              
                
                  Ψ
                  q
                
                ⁡
                
                  (
                  
                    
                      
                        I
                        d
                      
                      +
                      
                        Δ
                        ⁢
                        
                          
                        
                        ⁢
                        
                          I
                          d
                        
                      
                    
                    ,
                    
                      I
                      q
                    
                    ,
                    
                      Φ
                      m
                    
                  
                  )
                
              
              -
              
                
                  Ψ
                  q
                
                ⁡
                
                  (
                  
                    
                      I
                      d
                    
                    ,
                    
                      I
                      q
                    
                    ,
                    
                      Φ
                      m
                    
                  
                  )
                
              
            
            
              Δ
              ⁢
              
                
              
              ⁢
              
                I
                d
              
            
          
        
      
    
    
      
        [
        
          Math
          .
          
            
          
          ⁢
          3
        
        ]
      
    
  




In Mathematical Formulae 3, Ψd is a magnetic flux of the d-axis and Ψq is a magnetic flux of the q-axis (magnetic fluxes obtained by converting the interlinkage magnetic fluxes of the U-, V-, and W-phases into those of the rotating d-q coordinate system). As expressed by Mathematical Formulae 3, Ldh, Lqh, Ldqhf and Lqdh depend on both Id and Iq due to magnetic saturation by the currents Id and Iq.
When only the high-frequency component is extracted to analyze the error Δθ due to the high-frequency signal (error generated by the d- and q-axes mutual interference), Mathematical Formula 2 can be approximated by Mathematical Formula 4.




  
    
      
        
          [
          
            
              
                
                  V
                  dh
                
              
            
            
              
                
                  V
                  qh
                
              
            
          
          ]
        
        =
        
          
            [
            
              
                
                  
                    L
                    dh
                  
                
                
                  
                    L
                    dqh
                  
                
              
              
                
                  
                    L
                    qdh
                  
                
                
                  
                    L
                    qh
                  
                
              
            
            ]
          
          ⁢
          
            p
            ⁡
            
              [
              
                
                  
                    
                      I
                      dh
                    
                  
                
                
                  
                    
                      I
                      qh
                    
                  
                
              
              ]
            
          
        
      
    
    
      
        [
        
          Math
          .
          
            
          
          ⁢
          4
        
        ]
      
    
  




Here, a rotation matrix T of Mathematical Formula 5 is defined to convert the voltage equation of the d-q axis coordinate system of Mathematical Formula 4 into that of the estimated d-q axis coordinate system, and Mathematical Formula 6 is obtained.




  
    
      
        
          T
          ⁡
          
            (
            Δθ
            )
          
        
        =
        
          [
          
            
              
                
                  cos
                  ⁡
                  
                    (
                    Δθ
                    )
                  
                
              
              
                
                  sin
                  ⁡
                  
                    (
                    Δθ
                    )
                  
                
              
            
            
              
                
                  -
                  
                    sin
                    ⁡
                    
                      (
                      Δθ
                      )
                    
                  
                
              
              
                
                  cos
                  ⁡
                  
                    (
                    Δθ
                    )
                  
                
              
            
          
          ]
        
      
    
    
      
        [
        
          Math
          .
          
            
          
          ⁢
          5
        
        ]
      
    
  




Here, θre is an estimated phase (estimate value of the magnetic-pole position), θr is a true phase (actual magnetic-pole position), and the error Δθ is θre-θr (see FIG. 2).




  
    
      
        
          
            
              
                [
                
                  
                    
                      
                        V
                        dh
                        e
                      
                    
                  
                  
                    
                      
                        V
                        qh
                        e
                      
                    
                  
                
                ]
              
              =
              
              ⁢
              
                
                  T
                  ⁡
                  
                    (
                    Δθ
                    )
                  
                
                ⁡
                
                  [
                  
                    
                      
                        
                          V
                          dh
                        
                      
                    
                    
                      
                        
                          V
                          qh
                        
                      
                    
                  
                  ]
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                [
                
                  
                    
                      
                        
                          L
                          avg
                        
                        -
                        
                          
                            L
                            dif
                            *
                          
                          ⁢
                          
                            cos
                            ⁡
                            
                              (
                              
                                
                                  2
                                  ⁢
                                  Δθ
                                
                                +
                                
                                  θ
                                  m
                                
                              
                              )
                            
                          
                        
                      
                    
                    
                      
                        
                          
                            L
                            dif
                            *
                          
                          ⁢
                          
                            sin
                            ⁡
                            
                              (
                              
                                
                                  2
                                  ⁢
                                  Δθ
                                
                                +
                                
                                  θ
                                  m
                                
                              
                              )
                            
                          
                        
                        -
                        
                          L
                          dif
                          *
                        
                      
                    
                  
                  
                    
                      
                        
                          
                            L
                            dif
                            *
                          
                          ⁢
                          
                            sin
                            ⁡
                            
                              (
                              
                                
                                  2
                                  ⁢
                                  Δθ
                                
                                +
                                
                                  θ
                                  m
                                
                              
                              )
                            
                          
                        
                        +
                        
                          L
                          dif
                          *
                        
                      
                    
                    
                      
                        
                          L
                          avg
                        
                        +
                        
                          
                            L
                            dif
                            *
                          
                          ⁢
                          
                            cos
                            ⁡
                            
                              (
                              
                                
                                  2
                                  ⁢
                                  Δθ
                                
                                +
                                
                                  θ
                                  m
                                
                              
                              )
                            
                          
                        
                      
                    
                  
                
                ]
              
            
          
        
        
          
            
              
              ⁢
              
                p
                ⁡
                
                  [
                  
                    
                      
                        
                          I
                          dh
                          e
                        
                      
                    
                    
                      
                        
                          I
                          qh
                          e
                        
                      
                    
                  
                  ]
                
              
            
          
        
      
    
    
      
        [
        
          Math
          .
          
            
          
          ⁢
          6
        
        ]
      
    
  




Here, Lavg, Ldif, θm, and the like contained in Mathematical Formula 6 are expressed by Mathematical Formulae 7. The relationship between Lavg, Ldif, θm, and the like is shown in FIG. 5.




  
    
      
        
          
            L
            avg
          
          =
          
            
              
                L
                qh
              
              +
              
                L
                dh
              
            
            2
          
        
        ⁢
        
          

        
        ⁢
        
          
            L
            dif
          
          =
          
            
              
                L
                qh
              
              -
              
                L
                dh
              
            
            2
          
        
        ⁢
        
          

        
        ⁢
        
          
            θ
            m
          
          =
          
            
              
                tan
                
                  -
                  1
                
              
              ⁢
              
                
                  L
                  avg
                  *
                
                
                  L
                  dif
                
              
            
            =
            
              
                tan
                
                  -
                  1
                
              
              ⁡
              
                (
                
                  
                    
                      L
                      dqh
                    
                    +
                    
                      L
                      qdh
                    
                  
                  
                    
                      L
                      qh
                    
                    -
                    
                      L
                      dh
                    
                  
                
                )
              
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            L
            dif
            *
          
          =
          
            
              
                L
                dif
                2
              
              +
              
                L
                avg
                
                  *
                  2
                
              
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            L
            avg
            *
          
          =
          
            
              
                L
                dqh
              
              +
              
                L
                qdh
              
            
            2
          
        
        ⁢
        
          

        
        ⁢
        
          
            L
            dif
            *
          
          =
          
            
              
                L
                qdh
              
              +
              
                L
                dqh
              
            
            2
          
        
      
    
    
      
        [
        
          Math
          .
          
            
          
          ⁢
          7
        
        ]
      
    
  




In Mathematical Formulae 7, θm is a phase shift due to the presence of the d- and q-axes mutual interference. While Mathematical Formula 6 is as shown above, the voltage equation when the high-frequency voltage Vsig is applied to the estimated d-axis is expressed by Mathematical Formula 8.




  
    
      
        
          [
          
            
              
                
                  V
                  sig
                
              
            
            
              
                0
              
            
          
          ]
        
        =
        
          
            [
            
              
            
            ⁢
            
              
                
                  
                    
                      L
                      avg
                    
                    -
                    
                      
                        L
                        dif
                        *
                      
                      ⁢
                      
                        cos
                        ⁡
                        
                          (
                          
                            
                              2
                              ⁢
                              Δθ
                            
                            +
                            
                              θ
                              m
                            
                          
                          )
                        
                      
                    
                  
                
                
                  
                    
                      
                        L
                        dif
                        *
                      
                      ⁢
                      
                        sin
                        ⁡
                        
                          (
                          
                            
                              2
                              ⁢
                              Δθ
                            
                            +
                            
                              θ
                              m
                            
                          
                          )
                        
                      
                    
                    -
                    
                      L
                      dif
                      *
                    
                  
                
              
              
                
                  
                    
                      
                        L
                        dif
                        *
                      
                      ⁢
                      
                        sin
                        ⁡
                        
                          (
                          
                            
                              2
                              ⁢
                              Δθ
                            
                            +
                            
                              θ
                              m
                            
                          
                          )
                        
                      
                    
                    +
                    
                      L
                      dif
                      *
                    
                  
                
                
                  
                    
                      L
                      avg
                    
                    +
                    
                      
                        L
                        dif
                        *
                      
                      ⁢
                      
                        cos
                        ⁡
                        
                          (
                          
                            
                              2
                              ⁢
                              Δθ
                            
                            +
                            
                              θ
                              m
                            
                          
                          )
                        
                      
                    
                  
                
              
            
            ⁢
            
              
            
            ]
          
          ⁢
          
             
            
              p
              ⁡
              
                [
                
                  
                    
                      
                        I
                        dh
                        e
                      
                    
                  
                  
                    
                      
                        I
                        qh
                        e
                      
                    
                  
                
                ]
              
            
          
        
      
    
    
      
        [
        
          Math
          .
          
            
          
          ⁢
          8
        
        ]
      
    
  




By simplifying Mathematical Formula 8 for the current, Mathematical Formula 9 is obtained.




  
    
      
        
          
            
              
                p
                ⁡
                
                  [
                  
                    
                      
                        
                          I
                          dh
                          e
                        
                      
                    
                    
                      
                        
                          I
                          qh
                          e
                        
                      
                    
                  
                  ]
                
              
              =
              
              ⁢
              
                
                  
                    [
                    
                      
                    
                    ⁢
                    
                      
                        
                          
                            
                              L
                              avg
                            
                            -
                            
                              
                                L
                                dif
                                *
                              
                              ⁢
                              
                                cos
                                ⁡
                                
                                  (
                                  
                                    
                                      2
                                      ⁢
                                      Δθ
                                    
                                    +
                                    
                                      θ
                                      m
                                    
                                  
                                  )
                                
                              
                            
                          
                        
                        
                          
                            
                              
                                L
                                dif
                                *
                              
                              ⁢
                              
                                sin
                                ⁡
                                
                                  (
                                  
                                    
                                      2
                                      ⁢
                                      Δθ
                                    
                                    +
                                    
                                      θ
                                      m
                                    
                                  
                                  )
                                
                              
                            
                            -
                            
                              L
                              dif
                              *
                            
                          
                        
                      
                      
                        
                          
                            
                              
                                L
                                dif
                                *
                              
                              ⁢
                              
                                sin
                                ⁡
                                
                                  (
                                  
                                    
                                      2
                                      ⁢
                                      Δθ
                                    
                                    +
                                    
                                      θ
                                      m
                                    
                                  
                                  )
                                
                              
                            
                            +
                            
                              L
                              dif
                              *
                            
                          
                        
                        
                          
                            
                              L
                              avg
                            
                            +
                            
                              
                                L
                                dif
                                *
                              
                              ⁢
                              
                                cos
                                ⁡
                                
                                  (
                                  
                                    
                                      2
                                      ⁢
                                      Δθ
                                    
                                    +
                                    
                                      θ
                                      m
                                    
                                  
                                  )
                                
                              
                            
                          
                        
                      
                    
                    ⁢
                    
                      
                    
                    ]
                  
                  
                    -
                    1
                  
                
                ⁡
                
                  [
                  
                    
                      
                        
                          V
                          sig
                        
                      
                    
                    
                      
                        0
                      
                    
                  
                  ]
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  
                    V
                    sig
                  
                  
                    
                      L
                      avg
                      2
                    
                    -
                    
                      L
                      dif
                      
                        *
                        2
                      
                    
                    +
                    
                      L
                      dif
                      
                        *
                        2
                      
                    
                  
                
                ⁡
                
                  [
                  
                    
                      
                        
                          
                            L
                            avg
                          
                          +
                          
                            
                              L
                              dif
                              *
                            
                            ⁢
                            
                              cos
                              ⁡
                              
                                (
                                
                                  
                                    2
                                    ⁢
                                    Δθ
                                  
                                  +
                                  
                                    θ
                                    m
                                  
                                
                                )
                              
                            
                          
                        
                      
                    
                    
                      
                        
                          
                            
                              -
                              
                                L
                                dif
                                *
                              
                            
                            ⁢
                            
                              sin
                              ⁡
                              
                                (
                                
                                  
                                    2
                                    ⁢
                                    Δθ
                                  
                                  +
                                  
                                    θ
                                    m
                                  
                                
                                )
                              
                            
                          
                          -
                          
                            L
                            dif
                            *
                          
                        
                      
                    
                  
                  ]
                
              
            
          
        
      
    
    
      
        [
        
          Math
          .
          
            
          
          ⁢
          9
        
        ]
      
    
  




By integrating both sides of Mathematical Formula 9, Mathematical Formulae 10 are obtained as high-frequency currents of the estimated d- and q-axes.




  
    
      
        
          
            I
            dh
            e
          
          =
          
            
              
                V
                sig
              
              
                p
                ⁡
                
                  (
                  
                    
                      L
                      avg
                      2
                    
                    -
                    
                      L
                      dif
                      
                        *
                        2
                      
                    
                    +
                    
                      L
                      dif
                      
                        *
                        2
                      
                    
                  
                  )
                
              
            
            ⁢
            
              {
              
                
                  L
                  avg
                
                +
                
                  
                    L
                    dif
                    *
                  
                  ⁢
                  
                    cos
                    ⁡
                    
                      (
                      
                        
                          2
                          ⁢
                          Δθ
                        
                        +
                        
                          θ
                          m
                        
                      
                      )
                    
                  
                
              
              }
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            I
            qh
            e
          
          =
          
            
              
                -
                
                  V
                  sig
                
              
              
                p
                ⁡
                
                  (
                  
                    
                      L
                      avg
                      2
                    
                    -
                    
                      L
                      dif
                      
                        *
                        2
                      
                    
                    +
                    
                      L
                      dif
                      
                        *
                        2
                      
                    
                  
                  )
                
              
            
            ⁢
            
              {
              
                
                  
                    L
                    dif
                    *
                  
                  ⁢
                  
                    sin
                    ⁡
                    
                      (
                      
                        
                          2
                          ⁢
                          Δθ
                        
                        +
                        
                          θ
                          m
                        
                      
                      )
                    
                  
                
                +
                
                  L
                  dif
                  *
                
              
              }
            
          
        
      
    
    
      
        [
        
          Math
          .
          
            
          
          ⁢
          10
        
        ]
      
    
  




Here, when the control is performed so that Iqhe=0 is satisfied, that is, when the high-frequency voltage is applied only to the estimated d-axis, Mathematical Formula 11 is obtained.




  
    
      
        
          I
          qh
          e
        
        =
        
          
            
              
                -
                
                  V
                  sig
                
              
              
                p
                ⁡
                
                  (
                  
                    
                      L
                      avg
                      2
                    
                    -
                    
                      L
                      dif
                      
                        *
                        2
                      
                    
                    +
                    
                      L
                      dif
                      
                        *
                        2
                      
                    
                  
                  )
                
              
            
            ⁢
            
              {
              
                
                  
                    L
                    dif
                    *
                  
                  ⁢
                  
                    sin
                    ⁡
                    
                      (
                      
                        
                          2
                          ⁢
                          Δθ
                        
                        +
                        
                          θ
                          m
                        
                      
                      )
                    
                  
                
                +
                
                  L
                  dif
                  *
                
              
              }
            
          
          =
          0
        
      
    
    
      
        [
        
          Math
          .
          
            
          
          ⁢
          11
        
        ]
      
    
  




Deduced from Mathematical Formula 11, the error Δθ generated due to the d- and q-axes mutual interference is expressed by Mathematical Formula 12.




  
    
      
        Δθ
        =
        
          
            
              1
              2
            
            ⁢
            
              
                sin
                
                  -
                  1
                
              
              ⁡
              
                (
                
                  -
                  
                    
                      L
                      dif
                      *
                    
                    
                      L
                      dif
                      *
                    
                  
                
                )
              
            
          
          -
          
            
              1
              2
            
            ⁢
            
              
                tan
                
                  -
                  1
                
              
              ⁡
              
                (
                
                  
                    
                      L
                      dqh
                    
                    +
                    
                      L
                      qdh
                    
                  
                  
                    
                      L
                      qh
                    
                    -
                    
                      L
                      dh
                    
                  
                
                )
              
            
          
        
      
    
    
      
        [
        
          Math
          .
          
            
          
          ⁢
          12
        
        ]
      
    
  




Further, for a motor in practical use, it can be considered that Ldqh and Lqdh are approximately equal to each other. Therefore, Mathematical Formula 13 is obtained.




  
    
      
        Δθ
        =
        
          
            -
            
              1
              2
            
          
          ⁢
          
            
              tan
              
                -
                1
              
            
            ⁡
            
              (
              
                
                  2
                  ⁢
                  
                    L
                    dqh
                  
                
                
                  
                    L
                    qh
                  
                  -
                  
                    L
                    dh
                  
                
              
              )
            
          
        
      
    
    
      
        [
        
          Math
          .
          
            
          
          ⁢
          13
        
        ]
      
    
  




According to Mathematical Formula 13, when there is no d- and q-axes mutual interference (Ldqh=Lqdh=0), the error Δθ becomes zero.
Here, it is conceivable to perform control by using the error Δθ in Mathematical Formula 12 as the correction amount. In such a case, however, the control becomes non-linear control, and therefore is difficult to perform. In consideration of performing the correction which is started at a point at which θre starts slightly shifting from θr (at a point at which Δθ is approximately zero), the control can be linearized, supposing that the error Δθ is approximately zero.
Therefore, supposing that Δθ is approximately zero, cos (Δθ) can be approximated as 1 and sin (Δθ) can be approximated as Δθ. Thus, Mathematical Formulae 14 are obtained.

cos(Δθ)≈1,sin(Δθ)≈Δθ

cos(2Δθ+θm)≈ cos(θm)−2Δθ·sin(θm)

sin(2Δθ+θm)≈2Δθ·cos(θm)+sin(θm)  [Math. 14]

As a result, Idhe and Iqhe in Mathematical Formulae 10 are approximated as expressed by Mathematical Formulae 15.




  
    
      
        
          
            
              
                I
                dh
                e
              
              =
              
              ⁢
              
                
                  
                    V
                    sig
                  
                  
                    p
                    ⁡
                    
                      (
                      
                        
                          L
                          avg
                          2
                        
                        -
                        
                          L
                          dif
                          
                            *
                            2
                          
                        
                        +
                        
                          L
                          dif
                          
                            *
                            2
                          
                        
                      
                      )
                    
                  
                
                ⁢
                
                  {
                  
                    
                      L
                      avg
                    
                    +
                    
                      
                        L
                        dif
                        *
                      
                      ⁢
                      
                        cos
                        ⁡
                        
                          (
                          
                            
                              2
                              ⁢
                              Δθ
                            
                            +
                            
                              θ
                              m
                            
                          
                          )
                        
                      
                    
                  
                  }
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  
                    V
                    sig
                  
                  
                    p
                    ⁡
                    
                      (
                      
                        
                          L
                          avg
                          2
                        
                        -
                        
                          L
                          dif
                          
                            *
                            2
                          
                        
                        +
                        
                          L
                          dif
                          
                            *
                            2
                          
                        
                      
                      )
                    
                  
                
                ⁢
                
                  (
                  
                    
                      L
                      qh
                    
                    -
                    
                      2
                      ⁢
                      
                        
                          L
                          avg
                          *
                        
                        ·
                        Δθ
                      
                    
                  
                  )
                
              
            
          
        
      
    
    
      
        [
        
          Math
          .
          
            
          
          ⁢
          15
        
        ]
      
    
  
  
    
      
        
          
            
              
                I
                qh
                e
              
              =
              
              ⁢
              
                
                  
                    -
                    
                      V
                      sig
                    
                  
                  
                    p
                    ⁡
                    
                      (
                      
                        
                          L
                          avg
                          2
                        
                        -
                        
                          L
                          dif
                          
                            *
                            2
                          
                        
                        +
                        
                          L
                          dif
                          
                            *
                            2
                          
                        
                      
                      )
                    
                  
                
                ⁢
                
                  {
                  
                    
                      
                        L
                        dif
                        *
                      
                      ⁢
                      
                        sin
                        ⁡
                        
                          (
                          
                            
                              2
                              ⁢
                              Δθ
                            
                            +
                            
                              θ
                              m
                            
                          
                          )
                        
                      
                    
                    +
                    
                      L
                      dif
                      *
                    
                  
                  }
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  
                    V
                    sig
                  
                  
                    p
                    ⁡
                    
                      (
                      
                        
                          L
                          avg
                          2
                        
                        -
                        
                          L
                          dif
                          
                            *
                            2
                          
                        
                        +
                        
                          L
                          dif
                          
                            *
                            2
                          
                        
                      
                      )
                    
                  
                
                ⁢
                
                  (
                  
                    
                      L
                      dqh
                    
                    -
                    
                      2
                      ⁢
                      
                        
                          L
                          dif
                        
                        ·
                        Δθ
                      
                    
                  
                  )
                
              
            
          
        
      
    
    
      
        
      
    
  




Here, it is defined that λ=Ldqh/Lqh, and then Mathematical Formula 16 is obtained from Mathematical Formulae 15.




  
    
      
        
          
            I
            qh
            e
          
          +
          
            λ
            ⁢
            
              
            
            ⁢
            
              I
              dh
              e
            
          
        
        =
        
          
            
              V
              sig
            
            
              p
              ⁡
              
                (
                
                  
                    L
                    avg
                    2
                  
                  -
                  
                    L
                    dif
                    
                      *
                      2
                    
                  
                  +
                  
                    L
                    dif
                    
                      *
                      2
                    
                  
                
                )
              
            
          
          ⁢
          
            
              (
              
                
                  
                    -
                    2
                  
                  ⁢
                  
                    L
                    dif
                  
                
                -
                
                  2
                  ⁢
                  λ
                  ⁢
                  
                    
                  
                  ⁢
                  
                    L
                    avg
                    *
                  
                
              
              )
            
            ·
            Δθ
          
        
      
    
    
      
        [
        
          Math
          .
          
            
          
          ⁢
          16
        
        ]
      
    
  




Then, deduced from this Mathematical Formula 16, Δθ is given by Mathematical Formulae 17.




  
    
      
        
          
            
              
                Δθ
                =
                
                ⁢
                
                  
                    
                      (
                      
                        
                          L
                          avg
                          2
                        
                        -
                        
                          L
                          dif
                          
                            *
                            2
                          
                        
                        +
                        
                          L
                          dif
                          
                            *
                            2
                          
                        
                      
                      )
                    
                    ⁢
                    
                      p
                      ⁡
                      
                        (
                        
                          
                            I
                            qh
                            e
                          
                          +
                          
                            λ
                            ⁢
                            
                              
                            
                            ⁢
                            
                              I
                              dh
                              e
                            
                          
                        
                        )
                      
                    
                  
                  
                    
                      -
                      2
                    
                    ⁢
                    
                      
                        V
                        sig
                      
                      ⁡
                      
                        (
                        
                          
                            L
                            dif
                          
                          -
                          
                            λ
                            ⁢
                            
                              
                            
                            ⁢
                            
                              L
                              avg
                              *
                            
                          
                        
                        )
                      
                    
                  
                
              
            
          
          
            
              
                =
                
                ⁢
                
                  
                    k
                    
                      H
                      ⁢
                      
                        
                      
                      ⁢
                      F
                    
                    Δθ
                  
                  ⁡
                  
                    (
                    
                      
                        I
                        qh
                        e
                      
                      +
                      
                        λ
                        ⁢
                        
                          
                        
                        ⁢
                        
                          I
                          dh
                          e
                        
                      
                    
                    )
                  
                
              
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            k
            
              H
              ⁢
              
                
              
              ⁢
              F
            
            Δθ
          
          =
          
            
              -
              
                
                  
                    L
                    avg
                    2
                  
                  -
                  
                    L
                    dif
                    
                      *
                      2
                    
                  
                  +
                  
                    L
                    dif
                    
                      *
                      2
                    
                  
                
                
                  2
                  ⁢
                  
                    
                      V
                      sig
                    
                    ⁡
                    
                      (
                      
                        
                          L
                          dif
                        
                        -
                        
                          λ
                          ⁢
                          
                            
                          
                          ⁢
                          
                            L
                            avg
                            *
                          
                        
                      
                      )
                    
                  
                
              
            
            ⁢
            p
          
        
      
    
    
      
        [
        
          Math
          .
          
            
          
          ⁢
          17
        
        ]
      
    
  




As described above, when the control is performed so as to satisfy: Ieqh+λIedh=0, the error Δθ becomes zero. The positions of the actual d- and q-axes are unknown at the time of the control, and therefore it is obvious that Iegh and Iedh are currents on the estimated d- and q-axes. Thus, the superscript “e” is omitted from Iegh and Iedh in the following description.
However, the differential operator p is contained in k^HFΔθ in Mathematical Formulae 17, and hence there is a fear that the value of Iqh+λIdh may oscillate to make the control difficult. Therefore, it is preferred to perform the control so that the integrated value of Iqh+λIdh becomes zero.
As to the inductances contained in the above-mentioned value k^HFΔθ, the relationships: Lavg>L^dif*, Ldif>0 and Lavg*>0 are satisfied. Therefore, the parts of above-mentioned k^HFΔθ except for the differential operator p become always negative. Therefore, when the integrated value of Iqh+λIdh is positive, the error Δθ is negative. On the other hand, when the integrated value of Iqh+λIdh is negative, the error Δθ is positive. Thus, the PI feedback control can be applied, considering the integrated value of Iqh+λIdh as a deviation.
FIG. 6 is a flowchart illustrating an operation example of the motor controller 1. FIG. 6 illustrates an operation for correcting the estimate value of the magnetic-pole position among various operations realized by the motor controller 1.
In Step S1, the motor controller 1 applies the high-frequency voltage Vsig to the estimated d-axis. This operation is realized by the high-frequency superimposer 41 illustrated in FIG. 1.
In Step S2, the motor controller 1 convers the current detection values detected from the motor 9 from those of the three phases into those of the d-q axes. This operation is realized by the three-to-two phase converter 31 and the d-q converter 33 illustrated in FIG. 1.
In Step S3, the motor controller 1 acquires a superimposed component from the current detection values obtained by the conversion through the bandpass filter. This operation is realized by the bandpass filter 53 illustrated in FIG. 1.
In Step S4, the motor controller 1 performs the control so that Iqh+λIdh becomes zero. This operation is realized by the magnetic-pole position corrector 51 illustrated in FIGS. 1 and 3.
FIGS. 7A to 7D are graphs each showing a relationship between Δθ and λ corresponding to the magnetic-pole position. The horizontal axis indicates the magnetic-pole position, whereas the vertical axis indicates the error Δθ (axial shift) and the parameter λ (=Ldqh/Lqh) relating to the correction amount. The error Δθ is indicated by the broken line, whereas the parameter λ is indicated by the solid line.
As shown in FIGS. 7A to 7D, a change of the parameter λ in accordance with the magnetic-pole position corresponds to a change of the error Δθ in accordance with the magnetic-pole position. Therefore, the correction amount for the error Δθ corresponding to the magnetic-pole position can be generated by using the parameter λ.
The fact that the error Δθ changes depending on the magnetic-pole position is caused by the fact that the high-frequency inductances Ldqh and Lqh contained in the parameter λ changes depending on the magnetic-pole position. This is especially true in a concentrated-winding motor among various motors. Therefore, the correction of the estimate value of the magnetic-pole position in this embodiment is suitable for the case where the concentrated-winding motor is to be controlled.
A steady magnitude of the error Δθ depends on magnitudes of the d-axis current Id and the q-axis current Iq. Specifically, the error Δθ contains an oscillating component depending on the magnetic-pole position and a steady component depending on the magnitude of the current. Therefore, by using the parameter λ corresponding not only to the magnetic-pole position but also to the d-axis current Id and the q-axis current Iq, both of the components can be suppressed.
Now, a method of creating the control map is described.
First Example
A first example of the method of creating the control map is a method of calculating the parameter λ relating to the correction amount by using magnetic-field analysis.
The parameter λ(=Ldqh/Lqh) relating to the correction amount is calculated based on, for example, the d-axis interlinkage magnetic flux Ψd due to the d-axis current Id, the q-axis interlinkage magnetic flux Ψq due to the q-axis current Iq, the interlinkage magnetic flux Φm due to the magnet included in the motor 9, the d-axis current high-frequency component Idh, and the q-axis current high-frequency component Iqh. The parameter λ is calculated for each of the case where the magnetic-pole position θ is changed, the case where the d-axis current Id is changed, and the case where the q-axis current Iq is changed.
A specific description is now given. First, the d-axis interlinkage magnetic flux Ψd due to the d-axis current Id and the q-axis interlinkage magnetic flux Ψq due to the q-axis current Iq are calculated. For example, three-phase interlinkage magnetic fluxes Ψu, Ψv, and Ψw are calculated from the three-phase currents Iu, Iv, and Iw, and are further converted into the currents on the d- and q-axes. In this manner, the d-axis interlinkage magnetic flux Ψd and the q-axis interlinkage magnetic flux Ψq are calculated. The interlinkage magnetic flux Φm due to the magnet included in the motor 9 has an approximately constant value.
Next, by substituting ΔId and ΔIq for the d-axis current high-frequency component Idh and the q-axis current high-frequency component Iqh, which are high-frequency components appearing in the current when a high-frequency voltage (for example, 8 kHz and 40 V) is superimposed on the voltage applied to the motor 9, the high-frequency inductances Ldh, Lqh, Ldqh, and Lqdh are calculated from Mathematical Formulae 3. In this manner, the parameter λ(=Ldqh/Lqh) is calculated.
The parameter λ depends on the magnetic-pole position θ, the d-axis current Id, and the q-axis current Iq. The parameter λ is calculated by changing each of the three parameters (on which the parameter λ depends) independently by a certain fixed amount. An interval to change the magnetic-pole position θ can be, for example, 5 degrees of an electric angle. An interval to change the d-axis current Id and the q-axis current Iq can be, for example, 1/10th to 1⁄5th of their respective maximum magnitude. Then, the calculated values are associated with each other on the control map.
Second Example
A second example of the method of creating the control map is a method of calculating the parameter λ relating to the correction amount by using actually measured values.
The parameter λ(=Ldqh/Lqh) relating to the correction amount is calculated based on the d-axis current high-frequency component Idh and the q-axis current high-frequency component Iqh appearing in the current flowing through the motor 9 when the high-frequency voltage Vsig is superimposed on the voltage applied to the motor 9 with an encoder 99 mounted thereon as illustrated in FIG. 8, for example. The parameter λ is calculated for each of the case where the magnetic-pole position θ is changed, the case where the d-axis current Id is changed, and the case where the q-axis current Iq is changed.
A specific description is now given. When the high-frequency voltage Vsig is applied only to the d-axis in Mathematical Formula 4, Mathematical Formula 18 is obtained.




  
    
      
        
          [
          
            
              
                
                  V
                  sig
                
              
            
            
              
                0
              
            
          
          ]
        
        =
        
          
            [
            
              
                
                  
                    L
                    dh
                  
                
                
                  
                    L
                    dqh
                  
                
              
              
                
                  
                    L
                    dqh
                  
                
                
                  
                    L
                    qh
                  
                
              
            
            ]
          
          =
          
            [
            
              
                
                  
                    I
                    dh
                  
                
              
              
                
                  
                    I
                    qh
                  
                
              
            
            ]
          
        
      
    
    
      
        [
        
          Math
          .
          
            
          
          ⁢
          18
        
        ]
      
    
  




The d-axis current high-frequency component Idh and the q-axis current high-frequency component Iqh in this case are expressed by Mathematical Formulae 19 based on Mathematical Formulae 15.




  
    
      
        
          
            I
            dh
          
          =
          
            
              
                V
                sig
              
              
                p
                ⁡
                
                  (
                  
                    
                      L
                      avg
                      2
                    
                    -
                    
                      L
                      dif
                      2
                    
                  
                  )
                
              
            
            ⁢
            
              L
              qh
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            I
            qh
          
          =
          
            
              
                V
                sig
              
              
                p
                ⁡
                
                  (
                  
                    
                      L
                      avg
                      2
                    
                    -
                    
                      L
                      dif
                      2
                    
                  
                  )
                
              
            
            ⁢
            
              L
              dqh
            
          
        
      
    
    
      
        [
        
          Math
          .
          
            
          
          ⁢
          19
        
        ]
      
    
  




Deduced from these Mathematical Formulae 19, the parameter λ(=Ldqh/Lqh) is expressed by the following Mathematical Formula 20.




  
    
      
        λ
        =
        
          
            
              L
              dqh
            
            
              L
              qh
            
          
          =
          
            -
            
              
                I
                qh
              
              
                I
                dh
              
            
          
        
      
    
    
      
        [
        
          Math
          .
          
            
          
          ⁢
          20
        
        ]
      
    
  




Therefore, when the d-axis current high-frequency component Idh and the q-axis current high-frequency component Iqh are measured in a state in which the magnetic-pole position θ can be detected, the parameter λ can be calculated.
The d-axis current high-frequency component Idh and the q-axis current high-frequency component Iqh are measured in a state in which the encoder 99 is mounted to the motor 9 so that the magnetic-pole position θ is detectable, as illustrated in FIG. 8. A shaft of the motor 9 is coupled to a shaft of another motor 102 through an intermediation of a shaft connecting portion 109. The motor 102 and an encoder 101 are provided to fix the shaft of the motor 9 so as to prevent the shaft of the motor 9 from rotating.
Then, the high-frequency voltage Vsig is applied to the motor 9 in a state in which the predetermined magnetic-pole position θ, d-axis current Id, and q-axis current Iq are given. The d-axis current high-frequency component Idh and the q-axis current high-frequency component Iqh are measured. Thereafter, the parameter (=Ldqh/Lqh=−Iqh/Idh) is calculated.
As in the first example described above, the parameter λ is calculated by changing each of the three parameters (on which the parameter λ depends) independently by a certain fixed amount. An interval to change the magnetic-pole position θ can be, for example, 5 degrees of an electric angle. An interval to change the d-axis current Id and the q-axis current Iq can be, for example, 1⁄5th to 1/10th of their respective maximum magnitude. Then, the calculated values are associated with each other on the control map.
It should be understood by those skilled in the art that various modifications, combinations, sub-combinations and alterations may occur depending on design requirements and other factors insofar as they are within the scope of the appended claims or the equivalents thereof.

## Claims

1. A rotary electric machine controller, comprising:
a superimposer for superimposing a high-frequency signal on one of a voltage and a current of a rotary electric machine having saliency;
an estimator for calculating an estimate value of a magnetic-pole position of a rotor included in the rotary electric machine based on a high-frequency component appearing in the other one of the voltage and the current of the rotary electric machine;
a retainer for retaining correction amount relating values which are different from each other for at least two points of the magnetic-pole position of the rotor; and
a corrector for acquiring one of the correction amount relating values retained by the retainer corresponding to the estimate value of the magnetic-pole position of the rotor and correcting the estimate value of the magnetic-pole position of the rotor based on the acquired one of the correction amount relating values.

2. The rotary electric machine controller according to claim 1, wherein:
the retainer retains the correction amount relating values associated with respective points of a torque current (q-axis current), respective points of a field current (d-axis current), and respective points of the magnetic-pole position of the rotor of the rotary electric machine; and
the corrector acquires the one of the correction amount relating values corresponding to a detected value of the torque current, a detected value of the field current, and the estimate value of the magnetic-pole position of the rotor.

3. The rotary electric machine controller according to claim 1, wherein the correction amount relating values contains a high-frequency inductance due to d- and q-axes mutual interference, which changes in accordance with the magnetic-pole position of the rotor.

4. The rotary electric machine controller according to claim 1, wherein the corrector corrects the estimate value of the magnetic-pole position of the rotor so that D expressed by Mathematical Formula (1) approaches zero,




  
    
      
        [
        
          Math
          .
          
            
          
          ⁢
          1
        
        ]
      
    
    
      
        
      
    
  
  
    
      
        D
        =
        
          
            I
            qh
            e
          
          +
          
            λ
            ⁢
            
              
            
            ⁢
            
              I
              dh
              e
            
          
        
      
    
    
      
        (
        1
        )
      
    
  
  
    
      
        λ
        =
        
          
            L
            dqh
          
          
            L
            qh
          
        
      
    
    
      
        (
        2
        )
      
    
  



where λ is the one of the correction amount relating values and is a variable expressed by Mathematical Formula (2), Idhe is a high-frequency component of an estimated d-axis current, Iqhe is a high-frequency component of an estimated q-axis current, Ldqh is a high-frequency inductance due to d- and q-axes mutual interference, and Lqh is a high-frequency inductance of a q-axis.

5. The rotary electric machine controller according to claim 4, wherein the corrector corrects the estimate value of the magnetic-pole position of the rotor so that an integrated value of the D approaches zero.

6. The rotary electric machine controller according to claim 1, wherein the correction amount relating values are generated and retained for each of at least two points of the magnetic-pole position of the rotor based on a d-axis interlinkage magnetic flux due to a d-axis current, a q-axis interlinkage magnetic flux due to a q-axis current, an interlinkage magnetic flux due to a magnet included in the rotary electric machine, a high-frequency component of the d-axis current, and a high-frequency component of the q-axis current.

7. The rotary electric machine controller according to claim 1, wherein the correction amount relating values are generated and retained for each of at least two points of the magnetic-pole position of the rotor based on the high-frequency component appearing in the other one of the voltage and the current of the rotary electric machine provided with a magnetic-pole position detector when the high-frequency signal is superimposed on the one of the voltage and the current.

8. The rotary electric machine controller according to claim 1, wherein the rotary electric machine comprises a concentrated-winding motor.

9. A rotary electric machine control method, which is executed on a microprocessor executing information processing in accordance with a program stored in memory, the method comprising:
superimposing a high-frequency signal on one of a voltage and a current of a rotary electric machine having saliency;
calculating an estimate value of a magnetic-pole position of a rotor included in the rotary electric machine based on a high-frequency component appearing in the other one of the voltage and the current of the rotary electric machine;
retaining correction amount relating values which are different from each other for at least two points of the magnetic-pole position of the rotor; and
acquiring one of the correction amount relating values retained corresponding to the estimate value of the magnetic-pole position of the rotor; and
correcting the estimate value of the magnetic-pole position of the rotor based on the acquired one of the correction amount relating values.

10. A method of creating a control map to be retained in a rotary electric machine controller, which is executed on a microprocessor executing information processing accordance with in memory, the method comprising:
generating, for each of at least two points of a magnetic-pole position of a rotor included in a rotary electric machine having saliency, correction amount relating values for correcting an estimate value of the magnetic-pole position of the rotor, correction amount relating values being calculated based on a high-frequency component appearing in one of a voltage and a current of the rotary electric machine when a high-frequency signal is superimposed on the other one of the voltage and the current; and
associating, on the control map, the points of the magnetic-pole position of the rotor and the information relating to the correction amount relating values generated for each of the points with each other.