US 6,081,093 AGrant
Sensorless Control Method and Apparatus of Permanent Magnet Synchronous Motor
Issue Date:2000-06-27
•4 Claims
•3 Drawing Sheets
Abstract
The invention provides a sensorless control method and apparatus of a permanent magnet synchronous motor where magnetic axis can be assigned in all speed regions and speed can be controlled continuously irrespective of the speed command. In the control method, the d-q axis being the magnetic axis rotating in the real rotation speed .omega..sub.R of the motor is set and also the .gamma.-.delta. axis being the assigned magnetic axis of the motor is set, and when the rotation speed .omega..sub.R.sbsb..gamma. of the .gamma.-.delta. axis is determined, a distribution gain K1 is set so as to be decreased as the absolute value of the rotation speed command .omega..sub.RREF becomes large and a distribution gain K2 is set so as to be increased as the absolute value of the rotation speed command .omega..sub.RREF becomes large, and the rotation speed command .omega..sub.RREF is multiplied by K1 and the speed estimated value .omega..sub.RP determined from the induced voltage of the synchronous motor or the induced voltage estimated value is multiplied by K2 respectively, and both multiplied values are summed thereby the rotation speed .omega..sub.R.sbsb..gamma. of the .gamma.-.delta. axis being the assigned magnetic axis is determined.
Metadata
Assignee
- Kabushiki Kaisha Yaskawa Denki
Inventors
- Ryuichi Oguro
- Shuichi Fujii
- Masanobu Inazumi
Application Information
Application Number:US 3089193
Filing Date:1999-07-14
Priority Date:1996-12-05
Art Unit:287
Classifications
IPC:
H02P 528H02K 1742
Field of Search:
318700-839;139;245;254;439
Patent Drawings (3 sheets)
Description
Field of Art
The present invention relates to a sensorless control method and apparatus of a permanent magnet synchronous motor which does not have a position detector and a speed detector.
Background Art
When a brushless DC motor using a permanent magnet as a rotor is driven as a synchronous motor, it is necessary that the absolute position of the rotor is obtained and the correct current control is performed. In order to obtain the absolute value of the rotor, a rotor position detector such as an encoder or a resolver is generally used. However, since there are problems regarding complicateness of wiring or structure, price, use environment and the like, methods of estimating the position of the magnetic pole not using a rotor position detector have been proposed.
In the prior art, methods of estimating the position of the magnetic pole of the permanent magnet synchronous motor are known as follows.
[1] The Transaction of the Institute of Electrical Engineers of Japan, Vol. 113-D, No. 5, May 1993, p579 to 586.
[2] The Transaction of the Institute of Electrical Engineers of Japan, Vol. 114-D, No. 5, May, 1994, p591 to 592.
[3] The Transaction of the Institute of Electrical Engineers of Japan, Vol. 115-D, No. 4, April, 1995, p420 to 427.
[1] is a method that when stator currents I.sub..alpha., I.sub..beta. converted to .alpha.-.beta. axis coordinate system set on the stator are made observed values and stator voltages v.sub..alpha., v.sub..beta. are made inputs, magnetic flux .lambda..sub..alpha., .lambda..sub..beta. and rotor speed are estimated using adaptive rule.
[2] is a method that when stator currents i.sub..alpha., i.sub..beta. converted to .alpha.-.beta. axis coordinate system are made observed values and stator voltages v.sub..alpha., v.sub..beta. are made inputs, induced voltage .epsilon..sub..alpha. generated in the .alpha. axis direction in the .alpha.-.beta. axis coordinate system and induced voltage .epsilon..sub..beta. generated in the .beta. axis direction are estimated as disturbance.
[3] is a method that deviation angle .theta.e between the .gamma.-.delta. axis and the d-q axis is estimated from difference between stator currents i.sub..gamma., i.sub..delta. converted to the .gamma.-.delta. axis coordinate system set on the stator and rotating in the synchronous speed and current calculated values i.sub..gamma..spsb.0, i.sub..delta..spsb.0 calculated from model.
However, problems in the above-mentioned conventional methods are as follows.
Regarding [1], when the method is adopted in a permanent magnet synchronous motor having polarity, on the .alpha.-.beta. axis coordinates, the inductance becomes function of the rotor angle .theta.r of the motor and the state equation becomes complicated, and when an observer is to be constituted, the calculation amount is increased and the utilization is difficult. Also since the flux .lambda..sub..alpha., .lambda..sub..beta. is made unknown quantity, the state equation becomes non-observable at the rotor speed being zero and the estimator itself becomes unstable.
Regarding [2], since the induced voltage converted to the .alpha.-.beta. axis becomes alternating quantity, if the pole of the observer is not set large, the phase difference is produced between the actual quantity and the estimated quantity and the method becomes not to be used.
[3] is a simple method in comparison with [1], [2], moreover since the .gamma.-.delta. axis rotating at the angular velocity substantially synchronized with the d-q axis is considered as the reference, when the deviation .theta.e between the d-q axis and the .gamma.-.delta. axis is small, the state equation is not complicated and the method is excellent regarding the utilization. However, since the actual value is compared with the calculated value simply introduced from the model when the .gamma.-.delta. axis is coincident with the d-q axis, the deviation .theta.e can not be always estimated correctly due to the modeling error.
Further in the reference of [3], a method is taken in that the speed estimated value and the error between the d-axis and the .gamma. axis are estimated from the induced voltage or the estimated value of the induced voltage in the high speed region, and the speed control is performed while the .gamma. axis is coincident with the d-axis. However, this method can obtain the precision only in the region where the induced voltage is sufficiently high, and the estimation is impossible at the zero speed.
Therefore in the case of a synchronous motor of salient pole type, utilizing the property that the d-axis direction and the q-axis direction are different in inductance, if the inductance is measured while the motor is stopped, the d-axis can be known from the variation of the value of the inductance.
On the other hand, in the case of a synchronous motor of non-salient pole type, since the inductance is same in any pole, the above-mentioned method can not be used. When the synchronous motor has no load, if DC current flows in certain direction, the magnetic axis of the synchronous machine has property of moving so as to be coincident with the current flowing direction. Consequently if the current flows in the assigned magnetic axis, after lapse of sufficient time, the magnetic axis coincides with the assigned magnetic axis. As a result, the magnetic axis can be known.
However, these methods can be used only at the low speed state, and at the speed higher than certain speed, the estimation method must be changed to that described in the reference of [3]. This changing speed is different depending on sort of machines, and the rapid changing causes variation of the torque, and further different algorithms must be prepared at the low speed state basically. Therefore there is a problem that the design and the control are troublesome.
Disclosure of Invention
Accordingly problems to be solved by the invention are to provide a control method and apparatus where the magnetic axis can be assigned in all speed regions and the speed can be controlled continuously irrespective of the speed command.
In order to solve the above-mentioned problems, in a sensorless control method of a permanent magnet synchronous motor of the present invention, to control a region from zero speed to high speed continuously in a permanent magnet synchronous motor having a permanent magnet as a rotor, in the .alpha.-.beta. space coordinate system where one phase of the stator of the synchronous motor is made .alpha. axis and the forward rotating direction in the electric angle 90.degree. from the .alpha. axis is made .beta. axis, the coordinate d-q axis rotating in the real motor rotation speed .omega..sub.R is set where the magnetic axis of the motor is made d and the axis of leading in 90.degree. from the d-axis is made q, and the .gamma.-.delta. axis is set where the assigned magnetic axis of the motor is made .gamma. and the axis of leading in 90.degree. from the .gamma. axis is made .delta., and when the rotation speed .omega..sub.R.sbsb..gamma. of the .gamma.-.delta. axis is determined, the distribution gain K1 set so as to be decreased as the absolute value of the rotation speed command .omega..sub.RREF becomes large and the distribution gain K2 set so as to be increased as the absolute value of the rotation speed command .omega..sub.RREF becomes large are prepared, and the rotation speed command .omega..sub.RREF is multiplied by K1 and the speed estimated value .omega..sub.RP determined from the induced voltage of the synchronous motor or the induced voltage estimated value is multiplied by K2 respectively, and both multiplied values are summed thereby the rotation speed .omega..sub.R.sbsb..gamma. of the .delta.-.delta. axis of the assigned magnetic axis is determined.
Further when the .gamma. axis current command value i.sub..gamma..spsb.REF is to be determined, the low speed region command value i.sub..gamma..spsb.REFL is multiplied by the distribution gain K1 and the normal speed region command value i.sub..gamma..spsb.REFH is multiplied by the distribution gain K2 respectively, and both multiplied values are summed thereby the .gamma. axis current command value i.sub..gamma..spsb.REF can be determined, and also the .delta. axis current command value i.sub..delta..spsb.REF can be constituted by the proportion control term comprising the speed command values .omega..sub.REF and .omega..sub.RP and the integral control term comprising the speed command values .omega..sub.RREF and .omega..sub.R.sbsb..gamma..
Also a sensorless control apparatus of a permanent magnet synchronous motor of the present invention is provided with a distribution gain generator setting a distribution gain K1 to be decreased as the absolute value of a rotation speed command .omega..sub.RREF becomes large and a distribution gain K2 to be increased as the absolute value of a rotation speed command .omega..sub.RREF becomes large, .gamma. axis current command generator generating command i.sub..gamma..spsb.REF based on .gamma. axis current command i.sub..gamma..spsb.REFL at low speed state, .gamma. axis current command i.sub..gamma..spsb.REFH at high speed state and the distribution gains K1, K2, a speed controller outputting .delta. axis current i.sub..delta..spsb.REF based on rotation speed command .omega..sub.RREF, .delta. axis current controller outputting axis voltage command v.sub..gamma..spsb.REF based on the .delta. axis current i.sub..delta..spsb.REF, .gamma. axis current controller outputting .gamma. axis voltage command v.sub..gamma..spsb.REF based on the .gamma. axis current command i.sub..gamma..spsb.REF, a vector control circuit outputting amount of voltage and position angle to an inverter circuit of the synchronous motor based on the .delta. axis voltage command V.sub..delta..spsb.REF and the .delta. axis voltage command V.sub..gamma..spsb.REF, a three phases/two phases converter generating .delta. axis current i.sub..delta. and .gamma. axis current i.sub..gamma. based on currents i.sub.u, i.sub.w of two phases of the synchronous motor, .gamma.-.delta. axis, current, induced voltage estimator outputting induced voltage estimated values .epsilon..sub..delta..spsb.est (k+1) and .epsilon..gamma..spsb.est(k+1) based on the voltage commands v.sub..delta..spsb.REF, v.sub..gamma..spsb.REF, the .delta. axis current i.sub..delta. and the .gamma. axis current i.sub..gamma. and also outputting .gamma. axis current estimated value i.sub..gamma..spsb.est (k+1) to the .gamma. axis current controller and .delta. axis current estimated value i.sub..delta..spsb.est (k+1) to the .delta. axis current controller respectively, a proportion control speed operator operating command control speed .omega..sub.RPest based on the induced voltage estimated values .epsilon..sub..delta..spsb.est (k+1), .epsilon..sub..gamma..spsb.est (k+1), an integral magnetic axis rotating speed operator outputting the rotating speed .omega..sub.R.spsb..gamma.est of the magnetic axis based on the distribution gains K1, K2 and the command control speed .omega..sub.RPest, and a magnetic axis operator inputting amount of voltage and position angle to the inverter circuit based on the rotating speed .omega..sub.R.sbsb..gamma..sub.est of the magnetic axis.
According to the present invention, following effects are obtained.
(1) In the low speed region where the ratio of K1 is designed larger than that of K2, since the rotation speed .omega..sub.R.sbsb..gamma. of the .gamma. axis is near the speed command value .omega..sub.RREF, the d-axis being the real magnetic axis is rotated at the same speed as that of the .gamma. axis in the mean value and at the speed substantially equal to the speed command value. On the other hand, in the high speed region, since the ratio of K2 is designed sufficiently larger than that of K1 and also .omega..sub.RP and .omega..sub.R.sbsb..gamma. are coincident, the proportion control term and the speed estimated value constituting the integral control term are coincident. As a result, the vector control is executed well.
(2) Since the distribution gain is set continuously, at the intermediate speed region between the low speed region and the high speed region, the error between the .gamma. axis and the d-axis is gradually corrected and the region from the low speed region to the high speed region is executed continuously by the same algorithm.
(3) When the .delta. axis current command value i.sub..delta..spsb.REF is constituted by the proportion control term comprising the speed command values .omega..sub.REF and .omega..sub.RP and the integral control term comprising the speed command values .omega..sub.RREF and .omega..sub.R.sbsb..gamma., the speed estimated value is fed back to the proportion command term comprising .omega..sub.RP and the transient oscillation of the d-axis is suppressed.
Brief Description of Drawings
FIG. 1 is a flow chart showing steps of the invention;
FIG. 2 is a block diagram of control system executing the invention; and
FIG. 3 is an explanation diagram showing relation of distribution gains K1, K2 according to the invention.
Best Mode for Carrying Out the Invention
Embodiments of the present invention will be described according to a flow chart in FIG. 1.
Distribution gains K1, K2 (K1+K2=1) are function of speed command .omega..sub.RREF, and as shown in FIG. 3, setting is performed so that K1 becomes gain 0 at high speed, and varies linearly at the intermediate region between the low speed and the high speed.
Currents i.sub.u, i.sub.w of the U phase and the W phase are inputted from a current sensor (step 1).
Magnetic axis position .theta..sub.est (k) of the .gamma. axis assigned at the previous operation loop is used, and the three phases/two phases conversion is executed and the actual current values i.sub..gamma. (k), i.sub..delta. (k) in the .gamma.-.delta. coordinate system are led (step 2).
Voltage commands v.sub..delta..spsb.REF (k), v.sub..delta..spsb.REF (k) converted to the .gamma.-.delta. coordinate system and outputted at the time of (k+1) seconds are inputted (step 3).
According to formula (1), .gamma.-.delta. axis current estimated values i.sub..gamma..spsb.est (k+1), i .sub..delta..spsb.est (k+1) and induced voltage estimated values .epsilon..sub..gamma..spsb.est (k+1), .epsilon..sub..delta..spsb.est (k+1) at the time of (k+1) Ts seconds are estimated (step 4). ##EQU1##
Where Rs: stator resistance, Lq: q-axis inductance, Ld: d-axis inductance, Ts: sampling time, K.sub.1 to K.sub.8 : observer feedback gain.
As known well, since .epsilon..sub..delta..spsb.est (k+1)=K.epsilon..multidot..omega..sub.R .multidot.sin.theta.e, .epsilon..sub..delta..spsb.est (k+1)=K.epsilon..multidot..omega..sub.R .multidot.cos .theta.e, according to formula (2), temporary speed estimated value .omega..sub.R .o slashed..sub.est (k+1) is estimated (step 5). ##EQU2##
From sign of the induced voltage of the .delta.axis and the induced voltage of the .gamma. axis, according to formula (3), .omega..sub.RP (k+1) to be used in the proportion control is led (step 6).
According to the speed command value, using the distribution gain assigned as shown in FIG. 3, according to formula (4), the rotation speed of the .gamma. axis is assigned and the speed estimated value .omega..sub.R.spsb..gamma.est (k+1) to be used in the integral control is led (step 7).
From these speed estimated values, based on formula (5), position .theta..sub.est (k+1) of the .gamma. axis at the time of (k+1) Ts seconds is assigned (step 8).
Using the speed and the speed command estimated in step 7, according to formula (6), current command i.sub..delta..spsb.REF (k+1) of the .gamma.-.delta. axis is led. According to formula (7), current command i.sub..gamma..spsb.REF of the T axis is led (step 9). ##EQU3##
From the current command and the current, voltage commands v.sub..gamma..spsb.REF (k+1), v.sub..delta..spsb.REF (k+1) to be outputted at the time of (k+1) Ts seconds are estimated (step 10).
In the present invention, if the DC current i.sub..delta..spsb.REFL in the forward direction is generated to the .gamma. axis being arbitrary assigned axis, when the d-axis being the real magnetic axis exists in the phase delayed by the load angle .theta.e from the .gamma. axis, torque being proportional to i.sub..gamma..spsb.REFL sin .theta.e and directed to the .gamma. axis direction is generated in the magnetic axis. Therefore if the load torque is zero, the real magnetic axis is always subjected to the torque directed to the .gamma. axis being the assigned magnetic axis. In a synchronous machine usually having no damper winding, since the damping factor is substantially zero, the d-axis produces simple harmonic oscillation around the .gamma. axis. In the present invention, the speed estimated value is fed back to the proportional control term comprising .omega..sub.RP thereby the transient oscillation of the d-axis is suppressed. The speed used for the integral is coincident with the assigned speed of the .gamma. axis, and in the low speed region where the ratio of K1 is designed larger than that of K2, since the rotation speed .omega..sub.R.sbsb..gamma. of the .gamma. axis is near the speed command value .omega..sub.RREF, the d-axis being the real magnetic axis is rotated at the same speed as that of the .gamma. axis in the mean value and at the speed substantially equal to the speed command value.
On the other hand, in the high speed region, since the induced voltage has sufficient value, the rotation speed .omega..sub.R.sbsb..gamma. of the .gamma. axis is coincident with the rotation speed of the d-axis in the mean value and the setting so as to correct the angle error between the .gamma. axis and the d-axis becomes possible. In the high speed region, since the ratio of K2 is designed sufficiently larger than that of K1 and .omega..sub.RP and .omega..sub.R.sbsb..gamma. are coincident, the proportional control is coincident with the speed estimated value constituting the integral control term. As a result, the .gamma. axis is coincident with the d-axis moreover the speed command is coincident with the speed. In the sensorless vector control, the .gamma. axis is deemed as the magnetic axis but is coincident with the d-axis, thereby the vector control is executed well.
Moreover in the present invention, since the distribution gain is set continuously, in the intermediate region between the low speed region and the high speed region, the error between the .gamma. axis and the d-axis is corrected gradually and the region from the low speed region to the high speed region is executed continuously in the same algorithm.
FIG. 2 is a block diagram of the present invention. In FIG. 2, numeral 1 designates a distribution gain generator, numeral 2 designates a speed controller, numeral 3 designates .gamma. axis current command generator, numeral 4 designates .delta. axis current controller, numeral 5 designates .gamma. axis current controller, numeral 6 designates a vector control circuit, numeral 7 designates an inverter circuit, numeral 8 designates a synchronous motor, numeral 9 designates a three phases/two phases converter, numeral 10 designates .gamma.-.delta. axis, current, induced voltage estimator, numeral 11 designates a proportion control speed operator, numeral 12 designates a magnetic axis operator and numeral 13 designates an integral, magnetic axis rotation speed operator.
In FIG. 2, by the speed command .omega..sub.RREF, the distribution gain generator 1 produces distribution gains K1, K2 in FIG. 3. The .gamma. axis command i.sub..gamma..spsb.REF is produced in that the i.sub..gamma..spsb.REFL : low speed command, the i.sub..gamma..spsb.REFH : high speed command and K1, K2 are inputted to the .gamma. axis current command generator 3. The speed controller 2 outputs the .delta. axis current i.sub.67 .spsb.REF (k+1), and the .delta., .gamma. axes voltage commands v.sub..delta..spsb.REF, v.sub..gamma..spsb.REF are outputted from the .delta. axis current controller 4 and the .gamma. axis current controller 5 and are inputted to the vector control circuit 6. The output .theta..sub.est (k+1) from the magnetic axis operator 12 is inputted to the vector control circuit 6, and the amount of voltage and the position angle are inputted to the inverter circuit 7. The inverter circuit 7 supplies a current to the synchronous motor 8. Currents i.sub.u, i.sub.w are inputted to the three phases/two phases converter 9, and the rotor coordinate system .gamma.-.delta. axis current is produced. This current is inputted to the .gamma.-.delta. axis, current, induced voltage estimator 10, and also the voltage commands v.sub..delta..spsb.REF, v.sub..gamma..spsb.REF are inputted. By the .gamma.-.delta. axis, current, induced voltage estimator 10, the induced voltage estimated values .epsilon..sub..delta..spsb.est (k+1), .epsilon..sub..gamma..spsb.est (k+1) are inputted to the proportion control speed operator 11 and the command control speed .omega..sub.RPest is outputted. On the other hand, i.sub..gamma..spsb.est (k+1), i.sub..delta..spsb.est (k+1) are inputted to the .delta. axis current controller 4 and the .gamma. axis current controller 5, and the voltage command is produced. The .omega..sub.RPest and .omega..sub.RREF, and the distribution gains K1, K2 are inputted to the magnetic axis rotation speed operator 13, and the rotation speed .omega..sub.R.sbsb..gamma..sub.est of the magnetic axis is outputted.
Industrial Applicability
The present invention can be utilized in the field of sensorless speed control of a permanent magnet synchronous motor.
Claims
What is claimed is:
1. A sensorless control method of a permanent magnet synchronous motor wherein a permanent magnet synchronous motor having a permanent magnet as a rotor is controlled in a region from zero speed to high speed continuously, characterized in that in the .alpha.-.beta. space coordinate system where one phase is made .alpha. axis and the forward rotating direction in the electric angle 90.degree. from the .alpha. axis is made .beta. axis, the coordinate d-q axis rotating in the real motor rotation speed .omega..sub.R is set where the magnetic axis of the motor is made d and the axis of leading in 90.degree. from the d-axis is made q, and the .gamma.-.delta. axis is set where the assigned magnetic axis of the motor is made .gamma. and the axis leading in 90.degree. from the .gamma. axis is made .delta., and that when the rotation speed .omega..sub.R.sbsb..gamma. of the .gamma.-.delta. axis is determined, the distribution gain K1 is set so as to be decreased as the absolute value of the rotation speed command .omega..sub.RREF becomes large and the distribution gain K2 is set so as to be increased as the absolute value of the rotation speed command .omega..sub.RREF becomes large, and the rotation speed command .omega..sub.RREF is multiplied by K1 and the speed estimated value .omega..sub.RP determined from the induced voltage of the synchronous motor or the induced voltage estimated value is multiplied by K2 respectively, and both multiplied values are summed thereby the rotation speed .omega..sub.R.sbsb..gamma. of the .gamma.-.delta. axis of the assigned magnetic axis is determined.
2. A sensorless control method of a permanent magnet synchronous motor as set forth in claim 1, wherein when the .gamma. axis current command value i.sub.REFL is multiplied by the gain K1 and the normal speed region command value i.sub..gamma..spsb.REFH is multiplied by the distribution gain K2 respectively, and both multiplied values are summed thereby the .gamma. axis current command value i.sub..gamma..spsb.REF is determined.
3. A sensorless control method of a permanent magnet synchronous motor as set forth in claim 1, wherein the .delta. axis current command value i.sub..gamma..spsb.REF is constituted by the proportion control term comprising the speed command values .omega..sub.REF and .omega..sub.RP and the integral control term comprising the speed command values .omega..sub.RREF and .omega..sub.R.sbsb..gamma..
4. A sensorless control apparatus of a permanent magnet synchronous machine, comprising: a distribution gain generator setting a distribution gain K1 to be decreased as the absolute value of a rotation speed command .omega..sub.RREF becomes large and a distribution gain K2 to be increased as the absolute value of a rotation speed command .omega..sub.RREF becomes large; .gamma. axis current command generator generating the .gamma. axis command i.sub..gamma..spsb.REF based on the .gamma. axis current command i.sub..gamma..spsb.REFL at low speed state, the .gamma. axis current command i.sub..gamma..spsb.REFH at high speed state and the distribution gains K1, K2; a speed controller outputting the .delta. axis current i.sub..delta..spsb.REF based on the rotation speed command .omega..sub.RREF ; .delta. axis current controller outputting the .delta. axis voltage command v.sub..delta..spsb.REF based on the .delta. axis current i.sub..gamma..spsb.REF ; .gamma. axis current controller outputting the .gamma. axis current command v.sub..gamma..spsb.REF based on the .gamma. axis current command i.sub..gamma..spsb.REF ; a vector control circuit outputting amount of voltage and position angle to an inverter circuit of the synchronous motor based on the .delta. axis voltage command V.sub..delta..spsb.REF and the .gamma. axis voltage command v.sub..gamma..spsb.REF ; a three phases/two phases converter generating the .delta. axis current i.sub..delta. and the .gamma. axis current i.gamma. based on currents i.sub.u, i.sub.w of two phases of the synchronous motor; .gamma.-.delta. axis, current, induced voltage estimator outputting the induced voltage estimated values .epsilon..sub..delta..spsb.est (k+1) and .epsilon..sub..gamma..spsb.est (k+1) based on the voltage commands v.sub..delta..spsb.REF, v.sub..gamma..spsb.REF, the .delta. axis current i.sub..delta. and the .gamma. axis current i.sub..gamma. and also outputting the .gamma. axis current estimated value i.sub..gamma..spsb.est (k+1) to the .gamma. axis current controller and the .delta. axis current estimated value i.sub.67 .spsb.est (k+1) to the .delta. axis current controller respectively; a proportion control speed operator operating the command control speed .omega..sub.RPest based on the induced voltage estimated values .epsilon..sub..delta..spsb.est (k+1), .epsilon..sub..gamma..spsb.est (k+1); an integral magnetic axis rotation speed operator outputting the rotation speed .omega..sub.R.sbsb..gamma..sub.est of the magnetic axis based on the distribution gains K1, K2 and the command control speed .omega..sub.RPest ; and a magnetic axis operator inputting amount of voltage and position angle to the inverter circuit based on the rotation speed .omega..sub.R.sbsb..gamma..sub.est of the magnetic axis.