Control device for alternating-current motor

Abstract

In one coordinate system of arbitrary rectangular coordinate systems including a coordinate system in which the position of a stator is fixed, a coordinate system in which the position of a rotor is fixed, and a coordinate system which rotates at a rotational frequency which is n times (n is an integer which is not 0 or 1) that of the rotor, a filter on any other coordinate system is defined. Further, driving of a motor is controlled by utilizing this filter. As a result, a coordinate transformation operation need not be individually performed with respect to control variables.

Description

BACKGROUND OF THE INVENTION

1. Field of the Invention

The present invention relates to control over a driving current which is supplied to an alternating-current motor including a rotor and a stator.

2. Description of the Related Art

A typical driving current supplied to a three phase alternating-current motor including a rotor and a stator is a current having three phases of iu, iv and iw. This three-phase driving current is controlled based on an output torque command from the motor. In a heretofore common technique for controlling such a motor, a current having respective phases (the respective phases of u, v and w) is converted into currents of d and q axis coordinate systems of an exciting current axis (a d axis) and a torque current axis (a q axis), and each converted axis current is controlled to match with an axis command value obtained from a torque command of the motor.

Such control assumes that a motor driving current basically conforms to a sine wave, and such control is directed at this fundamental wave component only. However, in reality, a magnetic flux generated in accordance with a motor driving current is distorted, or a higher harmonic component is generated in a motor driving current due to various situations such as characteristics at the time of inverter switching or the like.

Therefore, in order to perform further accurate control, control must be carried out taking a higher harmonic component into consideration.

As a control method for a higher harmonic current as a frequency component n times that of a fundamental wave component, there has been proposed preparation of a coordinate system by which a current can be processed as a direct current with respect to all current components as control targets and execution of control over a current transformed on the coordinate system and a command value (See, for example, Japanese Patent Application Laid-open No. 2002-223600).

However, in the above-described art, control over a higher harmonic current component is executed on a coordinate system for the higher harmonic current component. Therefore, after performing coordinate transformation from three phases of u, v and w (an αβ phase: a coordinate system with a stator fixed) into an axis coordinate system or the like, control or the like is performed, and an output obtained by this control is again subjected to reverse coordinate transformation. Therefore, because many coordinate transformation operations must be carried out, calculations are undesirably complicated.

SUMMARY OF THE INVENTION

According to the present invention, a higher harmonic component can be also controlled. Therefore, a higher harmonic wave contained in a motor driving current can be suppressed, and a copper loss can be thereby reduced. Further, by applying an appropriate higher harmonic current to a higher harmonic component contained a magnet, an increase in an output torque can be expected. Furthermore, because control over any coordinate axis can be executed on a single coordinate axis, a coordinate transformation operation need not be performed each time, thereby effectively executing control. For example, PI control over a higher harmonic current which is n times the rotational frequency of a rotor can be executed in a dq axis coordinate system in which the position of the rotor is fixed, or PI control over a dq current in the dq axis coordinate system with the rotor fixed can be executed on an αβ coordinate system in which the position of the stator is fixed.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a view showing a configuration in a dq axis coordinate system for PI control over an nth-order higher harmonic wave;

FIG. 2 is a view showing another configuration in the dq axis coordinate system for PI control over the nth-order higher harmonic wave;

FIG. 3 is a view showing a structural example of a control system on a dq axis;

FIG. 4 is a view showing a configuration in an αβ coordinate system for PI control over a dq axis current;

FIG. 5 is a view showing a structural example of a control system on αβ;

FIG. 6 is a view showing a change in a motor current on the dq axis caused by a conventional technique;

FIG. 7 is a view showing a change in a motor current caused by a technique with the dq axis prepared for each higher harmonic wave;

FIG. 8 is a view showing a change in a motor current caused by a technique according to the embodiment which realizes on the dq axis the control on an ef axis (n=6); and

FIG. 9 is a view showing a configuration of a higher harmonic component in the coordinate system shown in FIG. 1.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT

A preferred embodiment of the present invention will be described hereinafter.

“Coordinate System”

An example will be considered of a rotor which rotates at a fixed angular velocity ω and a rotation angle θ. An αβ coordinate system in which the position of a stator is fixed, a dq axis coordinate system in which the position of a is rotor fixed, and an ef axis coordinate system rotating at a rotational velocity which is n times that of the rotor will be introduced.

Respective state quantities (column vectors) (x_α, x_β), (x_d, x_q), (x_e, x_f) on the αβ coordinate system, the dq axis coordinate system and the ef axis coordinate system have the following relationship through a transformation matrix T(θ).

$\begin{array}{cc}\left(\begin{array}{c}{x}_d\\ x_q\end{array}\right)=T\left(\theta \right)\left(\begin{array}{c}{x}_{\alpha }\\ x_{\beta }\end{array}\right)& \left(1\right)\\ \left(\begin{array}{c}{x}_e\\ x_f\end{array}\right)=T\left(n\theta \right)\left(\begin{array}{c}{x}_{\alpha }\\ x_{\beta }\end{array}\right)& \left(2\right)\\ T\left(n\theta \right)=\left(\begin{array}{cc}\cos n\theta & \sin n\theta \\ -\sin n\theta & \cos n\theta \end{array}\right)& \left(3\right)\end{array}$

PI control in each coordinate system will now be prepared as follows. It is to be noted that Kp is a constant for proportional control and Ki is a constant for integration control. Furthermore, a suffix r denotes a target value. Moreover, in the description of this specification, same fonts may be used to denote scalars, vectors and matrices.

(α–β Coordinate System)

$\begin{array}{cc}\left(\begin{array}{c}{v}_{\alpha }\\ v_{\beta }\end{array}\right)=\left(\begin{array}{cc}{K}_{p\alpha }& 0\\ 0& K_{p\beta }\end{array}\right)\left\{\left(\begin{array}{c}{i}_{\alpha r}\\ i_{\beta r}\end{array}\right)-\left(\begin{array}{c}{i}_{\alpha }\\ i_{\beta }\end{array}\right)\right\}+\left(\begin{array}{cc}{K}_{i\alpha }& 0\\ 0& K_{i\beta }\end{array}\right)\left(\begin{array}{c}{e}_{\alpha }\\ e_{\beta }\end{array}\right)& \left(4\right)\\ \frac{d}{dt}\left(\begin{array}{c}{e}_{\alpha }\\ e_{\beta }\end{array}\right)=\left(\begin{array}{c}{i}_{\alpha r}\\ i_{\beta r}\end{array}\right)-\left(\begin{array}{c}{i}_{\alpha }\\ i_{\beta }\end{array}\right)& \left(5\right)\end{array}$ (d–q Axis Coordinate System)

$\begin{array}{cc}\left(\begin{array}{c}{v}_d\\ v_q\end{array}\right)=\left(\begin{array}{cc}{K}_{\mathrm{pd}}& 0\\ 0& K_{\mathrm{pq}}\end{array}\right)\left\{\left(\begin{array}{c}{i}_{\mathrm{dr}}\\ i_{\mathrm{qr}}\end{array}\right)-\left(\begin{array}{c}{i}_d\\ i_q\end{array}\right)\right\}+\left(\begin{array}{cc}{K}_{\mathrm{id}}& 0\\ 0& K_{\mathrm{iq}}\end{array}\right)\left(\begin{array}{c}{e}_d\\ e_q\end{array}\right)& \left(6\right)\\ \frac{d}{dt}\left(\begin{array}{c}{e}_d\\ e_q\end{array}\right)=\left(\begin{array}{c}{i}_{\mathrm{dr}}\\ i_{\mathrm{qr}}\end{array}\right)-\left(\begin{array}{c}{i}_d\\ i_q\end{array}\right)& \left(7\right)\end{array}$ (e–f Axis Coordinate System)

$\begin{array}{cc}\left(\begin{array}{c}{v}_e\\ v_f\end{array}\right)=\left(\begin{array}{cc}{K}_{\mathrm{pe}}& 0\\ 0& K_{\mathrm{pf}}\end{array}\right)\left\{\left(\begin{array}{c}{i}_{\mathrm{er}}\\ i_{\mathrm{fr}}\end{array}\right)-\left(\begin{array}{c}{i}_e\\ i_f\end{array}\right)\right\}+\left(\begin{array}{cc}{K}_{\mathrm{ie}}& 0\\ 0& K_{\mathrm{if}}\end{array}\right)\left(\begin{array}{c}{e}_e\\ e_f\end{array}\right)& \left(8\right)\\ \frac{d}{dt}\left(\begin{array}{c}{e}_e\\ e_f\end{array}\right)=\left(\begin{array}{c}{i}_{\mathrm{er}}\\ i_{\mathrm{fr}}\end{array}\right)-\left(\begin{array}{c}{i}_e\\ i_f\end{array}\right)& \left(9\right)\end{array}$ “Higher Harmonic Control in dq Axis Coordinate System (PI Control Over ef Axis Coordinate System in dq Axis Coordinate System)”

The dq axis coordinate system and the ef axis coordinate system have the relationships represented by Expressions (1) and (2). Considering a voltage v, a current i and an error integration value (see Expression (5)) e of the current as state quantities, they can be expressed as follows. It is to be noted that T′((n−1)θ) is a transposed matrix of T((n−1)θ) and the dq axis coordinate system advances by θ with respect to αβ. Therefore, a difference in order between the both coordinate systems is n−1.

$\begin{array}{cc}\left(\begin{array}{c}{v}_d\\ v_q\end{array}\right)=T^{\prime }\left(\left(n-1\right)\theta \right)\left(\begin{array}{c}{v}_e\\ v_f\end{array}\right)& \left(10\right)\\ \left(\begin{array}{c}{i}_d\\ i_q\end{array}\right)=T^{\prime }\left(\left(n-1\right)\theta \right)\left(\begin{array}{c}{i}_e\\ i_f\end{array}\right)& \left(11\right)\\ \left(\begin{array}{c}{e}_d\\ e_q\end{array}\right)=T^{\prime }\left(\left(n-1\right)\theta \right)\left(\begin{array}{c}{e}_e\\ e_f\end{array}\right)& \left(12\right)\end{array}$

The PI control over the ef coordinate system represented by Expressions (8) and (9) is transformed into the dq axis coordinate system by utilizing Expressions (10) to (12).

$\begin{array}{cc}\left(\begin{array}{c}{v}_d\\ v_q\end{array}\right)=T^{\prime }\left(\left(n-1\right)\theta \right)\left(\begin{array}{cc}{K}_{\mathrm{pe}}& 0\\ 0& K_{\mathrm{pf}}\end{array}\right)T\left(\left(n-1\right)\theta \right)\left\{\left(\begin{array}{c}{i}_{\mathrm{dr}}\\ i_{\mathrm{qr}}\end{array}\right)-\left(\begin{array}{c}{i}_d\\ i_q\end{array}\right)\right\}+T^{\prime }\left(\left(n-1\right)\theta \right)\left(\begin{array}{cc}{K}_{\mathrm{ie}}& 0\\ 0& K_{\mathrm{if}}\end{array}\right)T\left(\left(n-1\right)\theta \right)\left(\begin{array}{c}{e}_d\\ e_q\end{array}\right)=\left(\begin{array}{cc}\frac{K_{\mathrm{pe}}+K_{\mathrm{pf}}}{2}+\frac{K_{\mathrm{pe}}-K_{\mathrm{pf}}}{2}\cos \left(2\left(n-1\right)\theta \right)& \frac{K_{\mathrm{pe}}-K_{\mathrm{pf}}}{2}\sin \left(2\left(n-1\right)\theta \right)\\ \frac{K_{\mathrm{pe}}-K_{\mathrm{pf}}}{2}\sin \left(2\left(n-1\right)\theta \right)& \frac{K_{\mathrm{pe}}+K_{\mathrm{pf}}}{2}-\frac{K_{\mathrm{pe}}-K_{\mathrm{pf}}}{2}\cos \left(2\left(n-1\right)\theta \right)\end{array}\right)\left\{\left(\begin{array}{c}{i}_{\mathrm{dr}}\\ i_{\mathrm{qr}}\end{array}\right)-\left(\begin{array}{c}{i}_d\\ i_q\end{array}\right)\right\}+\left(\begin{array}{cc}\frac{K_{\mathrm{ie}}+K_{\mathrm{if}}}{2}+\frac{K_{\mathrm{ie}}-K_{\mathrm{if}}}{2}\cos \left(2\left(n-1\right)\theta \right)& \frac{K_{\mathrm{ie}}-K_{\mathrm{if}}}{2}\sin \left(2\left(n-1\right)\theta \right)\\ \frac{K_{\mathrm{ie}}-K_{\mathrm{if}}}{2}\sin \left(2\left(n-1\right)\theta \right)& \frac{K_{\mathrm{ie}}+K_{\mathrm{if}}}{2}-\frac{K_{\mathrm{ie}}-K_{\mathrm{if}}}{2}\cos \left(2\left(n-1\right)\theta \right)\end{array}\right)\left(\begin{array}{c}{e}_d\\ e_q\end{array}\right)& \left(14\right)\\ \frac{d}{dt}\left(\begin{array}{c}{e}_d\\ e_q\end{array}\right)=\frac{d}{dt}\left\{T^{\prime }\left(\left(n-1\right)\theta \right)\left(\begin{array}{c}{e}_e\\ e_f\end{array}\right)\right\}=\frac{d}{dt}\left\{T^{\prime }\left(\left(n-1\right)\theta \right)\right\}T\left(\left(n-1\right)\theta \right)\left(\begin{array}{c}{e}_d\\ e_q\end{array}\right)+T^{\prime }\left(\left(n-1\right)\theta \right)\frac{d}{dt}\left(\begin{array}{c}{e}_e\\ e_f\end{array}\right)=\left(n-1\right)\omega \left(\begin{array}{cc}0& -1\\ 1& 0\end{array}\right)\left(\begin{array}{c}{e}_d\\ e_q\end{array}\right)+T^{\prime }\left(\left(n-1\right)\theta \right)\left\{\left(\begin{array}{c}{i}_{\mathrm{er}}\\ i_{\mathrm{fr}}\end{array}\right)-\left(\begin{array}{c}{i}_e\\ i_f\end{array}\right)\right\}=\left(n-1\right)\omega \left(\begin{array}{cc}0& -1\\ 1& 0\end{array}\right)\left(\begin{array}{c}{e}_d\\ e_q\end{array}\right)+\left(\begin{array}{c}{i}_{\mathrm{dr}}\\ i_{\mathrm{qr}}\end{array}\right)-\left(\begin{array}{c}{i}_d\\ i_q\end{array}\right)& \left(15\right)\end{array}$

Expression (15) can be converted into Expression (17). Additionally, assuming that K_pe=K_pf, K_ie=K_if, Expression (16) can be obtained by organizing the above expressions.

$\begin{array}{cc}\left(\begin{array}{c}{v}_d\\ v_q\end{array}\right)=K_{\mathrm{pe}}\left\{\left(\begin{array}{c}{i}_{\mathrm{dr}}\\ i_{\mathrm{qr}}\end{array}\right)-\left(\begin{array}{c}{i}_d\\ i_q\end{array}\right)\right\}+K_{\mathrm{ie}}\left(\begin{array}{c}{e}_d\\ e_q\end{array}\right)& \left(16\right)\\ \frac{d}{dt}\left(\begin{array}{c}{e}_d\\ e_q\end{array}\right)=\left(n-1\right)\omega \left(\begin{array}{cc}0& -1\\ 1& 0\end{array}\right)\left(\begin{array}{c}{e}_d\\ e_q\end{array}\right)+\left(\begin{array}{c}{i}_{\mathrm{dr}}\\ i_{\mathrm{qr}}\end{array}\right)-\left(\begin{array}{c}{i}_d\\ i_q\end{array}\right)& \left(17\right)\end{array}$

Based on a current error (i_dr−i_d, i_qr−i_q) in Expression (17), a transfer function of the integration value (e_d, e_q) is represented by Expression (18).

$\begin{array}{cc}F\left(\left(\begin{array}{c}{e}_d\\ e_q\end{array}\right)\right)=\frac{1}{s^2+{\{\left(n-1\right)\omega \}}^2}\left(\begin{array}{cc}s& \left(n-1\right)\omega \\ -\left(n-1\right)\omega & s\end{array}\right)F\left(\left(\begin{array}{c}{i}_{\mathrm{dr}}\\ i_{\mathrm{qr}}\end{array}\right)-\left(\begin{array}{c}{i}_d\\ i_q\end{array}\right)\right)& \left(18\right)\end{array}$

According to this expression, a diagonal component constitutes a band pass filter by which a pass band has a rotational frequency ω, and a non-diagonal component constitutes a low pass filter by which a cutoff frequency has a rotational frequency ω. Therefore, there is the possibility of robust properties with respect to the rotational frequency ω, performing transformation as follows can be considered.

$\begin{array}{cc}\begin{array}{c}F\left(\left(\begin{array}{c}{e}_d\\ e_q\end{array}\right)\right)=\frac{1}{s^2+2\zeta \left(n-1\right)\omega s+{\left\{\left(n-1\right)\omega \right\}}^2}\\ \left(\begin{array}{cc}s+\zeta \left(n-1\right)\omega & \left(n-1\right)\omega \\ -\left(n-1\right)\omega & s+\zeta \left(n-1\right)\omega \end{array}\right)\\ F\left(\left(\begin{array}{c}{i}_{\mathrm{dr}}\\ i_{\mathrm{qr}}\end{array}\right)-\left(\begin{array}{c}{i}_d\\ i_q\end{array}\right)\right)\end{array}& \left(19\right)\end{array}$

In this expression, F( ) is a Laplace transform, s is a Laplace operator, and ζ is a constant corresponding to damping and 0<ζ<0.7 can be considered as an appropriate value.

$\begin{array}{cc}\frac{d}{dt}\left(\begin{array}{c}{e}_d\\ e_q\end{array}\right)=\left(n-1\right)\omega \left(\begin{array}{cc}\frac{-\zeta }{\left(n-1\right)\omega }& -1\\ 1& \frac{-\zeta }{\left(n-1\right)\omega }\end{array}\right)\left(\begin{array}{c}{e}_d\\ e_q\end{array}\right)+\left(\begin{array}{c}{i}_{\mathrm{dr}}\\ i_{\mathrm{qr}}\end{array}\right)-\left(\begin{array}{c}{i}_d\\ i_q\end{array}\right)& \left(20\right)\end{array}$

Assuming that the above-described control system is a control system for an nth-order higher harmonic wave (ef) and Expressions (6) and (7) are control systems for a fundamental wave (dq), a control system for the fundamental wave and the nth-order higher harmonic wave is as shown in FIG. 1 or 2. Further, a configuration which realizes this control system is shown in FIG. 3.

In FIG. 1, processing with respect to a signal which has been transmitted through a low pass filter in order to take out a fundamental wave component corresponds to fundamental wave processing, and processing with respect to a signal which has been transmitted through a high pass filter in order to take out an nth-order higher harmonic component corresponds to nth-order higher harmonic wave processing.

i_dr and i_d, as well as i_qr and i_q of the fundamental wave transmitted through the low pass filter, are input to a subtracter where error values (i_dr−i_d) and (i_qr−i_q) are calculated. The obtained error values are each multiplied by K_pd, thereby calculating proportionals in the PI control over the fundamental wave. Furthermore, the error values (i_dr−i_d) and (i_qr−i_q) are subjected to integration (1/s) and then multiplied by K_id, thereby obtaining integration terms of the PI control. Furthermore, these results are added so that control voltages vd and vq of the fundamental wave can be obtained.

Moreover, i_dr and i_d, as well as i_qr and i_q of the higher harmonic component transmitted through the high pass filter, are input to the subtracter where error values (i_dr−i_d) and (i_qr−i_q) are calculated. The obtained error values are multiplied by K_p, thereby calculating proportionals in the PI control over the fundamental wave. Additionally, the error values (i_dr−i_d) and (i_qr−i_q) are subjected to integration (1/s) and then multiplied by K_ie, thereby obtaining integration terms of the PI control. However, an adder (a subtracter) is provided before each integrator. The integration term of the d axis is multiplied by (n−1)ω and then subtracted from the error of the q axis, and an obtained result is input to the integrator of the q axis. The integration term of the q axis is multiplied by (n−1)ω and the added to the error of the d axis, and an obtained result is input to the integrator of the d axis. Consequently, the control represented by Expression (15) and the like can be executed.

With such a configuration, the proportionals and the integration terms of the nth-order higher harmonic wave can be obtained with respect to the d axis and the q axis, and the proportional and the integration term of the fundamental wave and the proportional and the integration term of the higher harmonic wave are added in the adders in accordance with the d axis and the q axis, thereby obtaining the control voltage vd of the d axis and the control voltage vq of the q axis.

In this manner, the PI control in the ef axis can be collectively performed as the control of dq. It is to be noted that K_pd=K_pq, K_id=K_iq, K_pe=K_pf and K_ie=K_if are determined in this example.

FIG. 9 shows a processing part with respect to the higher harmonic wave in FIG. 1.

Further, in FIG. 2, a coefficient of proportional control is determined as K_pd=K_pe, and proportionals of the fundamental wave and the higher harmonic wave are altogether calculated. That is, i_dr and i_d of the d axis are added without using a filter, and an obtained result is multiplied by K_pd, thereby obtaining a proportional. Furthermore, i_qr and i_q of the q axis are added without using a filter, and an obtained result is multiplied by K_pd, thereby obtaining a proportional.

FIG. 3 shows the overall configuration of the system. The above-described configurations of FIGS. 1 and 2 are employed in a current compensator. This current compensator outputs vd and vq as voltage commands of the d axis and the q axis, and they are input to a dq/uvw converting section. The dq/uvw converting section converts the voltage command values of the dq axes into a switching command for an inverter which outputs each phase voltage driving voltage, and outputs the switching command. The switching command is input to a PWM inverter. Motor driving voltages v_u, v_v and v_w corresponding to vd and vq are supplied to respective phase coils of the three-phase motor 3 in accordance with the PWM inverter switching command.

On the other hand, a rotor rotating position of the motor is detected by a position sensor. The position sensor may be of a type which detects a change in any other motor current of a hall element. A detected value from the position sensor is input to an angle and angular velocity calculator where an angle θ and an angular velocity ω of the rotor are calculated from the rotor position detection result. This rotor angle θ is input to the uvw/dq converter. Motor currents having a v phase and a w phase (which may be any two phases or three phases) detected by the current detector are supplied to this uvw/dq converter where an exciting current id and a torque current iq in the dq axis coordinate system are calculated.

Moreover, id and iq from this uvw/dq converter and the angular velocity ω from the angle and angular velocity calculator are supplied to the current compensator. That is, a target value idr of the exciting current, a target value iqr of the torque current, and the detection results id, iq and ω are input to this current compensator, and hence vd and vq are calculated by such configurations as shown in FIGS. 1 and 2.

With such configurations, the motor driving control taking the higher harmonic wave into consideration can be executed without performing the coordinate transformation operation.

“dq Axis Current Control in αβ Coordinate System (PI Control Over dq Axis Coordinate System in αβ Coordinate System)”

When the PI control with the rotor fixed in the dq axis coordinate system is transformed into the αβ coordinate system with the stator fixed, the following expressions can be achieved.

$\begin{array}{cc}\left(\begin{array}{c}{v}_{\alpha }\\ v_{\beta }\end{array}\right)={\kappa }_{\mathrm{pd}}\left\{\left(\begin{array}{c}{i}_{\alpha r}\\ i_{\beta r}\end{array}\right)-\left(\begin{array}{c}{i}_{\alpha }\\ i_{\beta }\end{array}\right)\right\}+{\kappa }_{\mathrm{id}}\left(\begin{array}{c}{e}_{\alpha }\\ e_{\beta }\end{array}\right)& \left(21\right)\\ \frac{d}{dt}\left(\begin{array}{c}{e}_{\alpha }\\ e_{\beta }\end{array}\right)=\omega \left(\begin{array}{cc}0& -1\\ 1& 0\end{array}\right)\left(\begin{array}{c}{e}_{\alpha }\\ e_{\beta }\end{array}\right)+\left(\begin{array}{c}{i}_{\alpha }\\ i_{\beta r}\end{array}\right)-\left(\begin{array}{c}{i}_{\alpha }\\ i_{\beta }\end{array}\right)& \left(22\right)\end{array}$

They correspond to Expressions (16) and (17), a difference in order between the αβ coordinate system and the dq axis coordinate system is 1, and ω is adopted in place of (n−1)ω in Expression (17).

Therefore, FIG. 4 shows a block diagram for this control. Exciting current and torque current commands (target values) i_dr and i_qr are input to the dq/αβ converter where these values are converted into iα_r and iβ_r. It is to be noted that a rotor angle θ is required for this conversion, and this θ is also input to the dq/αβ converter.

iα_r and iα, along with iβ_r and iβ, are input to the subtracter where error values (iα_r−iα) and (iβ_r−iβ) are respectively calculated. The obtained error values are multiplied by K_pd, thereby calculating proportionals in the PI control of the fundamental wave. Additionally, the error values (iα_r−iα) and (iβ_r−iβ) are subjected to integration (1/s) and then multiplied by Kid, thereby obtaining integration terms of the PI control. However, an adder or a subtracter is provided before each integrator, and the integration term concerning the α axis is multiplied by ω and then added to the error of the β axis, and the obtained result is input to the integrator of the β axis. Further, the integration term concerning the α axis is multiplied by ω and then subtracted from the error of the β axis, and the obtained result is input to the integrator of the β axis. Consequently, the control represented by Expression (15) and the like can be executed.

With such a configuration, the proportional and the integration term of the PI control in the dq axis coordinate system can be obtained with respect to the αβ coordinate system, and the control voltage vα of the α axis and the control voltage vβ of the β axis an be obtained.

FIG. 5 shows the entire control system. This configuration is basically the same as that shown in FIG. 3. iα_r, iβ_r, iα, iβ and ω are input to the current compensator, and vα and vβ are obtained by the configuration shown in FIG. 4. vα and vβ are input to the αβ/uvw converter where commands of respective phases u, v and w are created, and the motor is thereby driven. Furthermore, motor currents iv and i_w are converted into iα and iβ in the uvw/aβ converter.

In this manner, the PI control in the dq axis coordinate system can be realized in the αβ coordinate system. “Higher Harmonic Control Method in dq Coordinate System (General Filter of ef Coordinate System in dq Coordinate System)”

In the above, an example in which the PI control (the control which calculates proportionals and integration terms of errors) in the ef axis coordinate system is carried out in the dq axis coordinate system was described. Next, an example which realizes in the dq axis coordinate system a general filter in the ef coordinate system (the coordinate system of n rotations) will be described.

It is to be noted that although in the following description the order of the filter is the fourth order, the order of the filter can be readily extended to the sixth order, the eighth order, etc. It is to be noted that m=n−1 is determined.

First, transformation into the dq axis coordinate system by the fourth-order filter corresponding to Expressions (10) to (12) can be represented as follows.

$\begin{array}{cc}\left(\begin{array}{c}{v}_d\\ v_q\end{array}\right)=T^{\prime }\left(m\theta \right)\left(\begin{array}{c}{v}_e\\ v_f\end{array}\right)& \left(23\right)\\ \left(\begin{array}{c}{i}_d\\ i_q\end{array}\right)=T^{\prime }\left(m\theta \right)\left(\begin{array}{c}{i}_e\\ i_f\end{array}\right)& \left(24\right)\\ \left(\begin{array}{c}{e}_d\\ e_q\end{array}\right)=T^{\prime }\left(m\theta \right)\left(\begin{array}{c}{e}_e\\ e_f\end{array}\right)& \left(25\right)\\ \left(\begin{array}{c}{i}_r\\ i_s\end{array}\right)=T^{\prime }\left(m\theta \right)\left(\begin{array}{c}{i}_g\\ i_h\end{array}\right)& \left(26\right)\\ \left(\begin{array}{c}{e}_r\\ e_s\end{array}\right)=T^{\prime }\left(m\theta \right)\left(\begin{array}{c}{e}_g\\ e_h\end{array}\right)& \left(27\right)\end{array}$

On the other hand, the fourth-order filter in the ef axis coordinate system can be represented as follows.

$\begin{array}{cc}\left(\begin{array}{c}{v}_e\\ v_f\end{array}\right)={\kappa }_p\left\{\left(\begin{array}{c}{i}_{\mathrm{er}}\\ i_{\mathrm{fr}}\end{array}\right)-\left(\begin{array}{c}{i}_e\\ i_f\end{array}\right)\right\}+{\kappa }_i\left(\begin{array}{c}{e}_e\\ e_f\\ e_g\\ e_h\end{array}\right)& \left(29\right)\\ \frac{d}{dt}\left(\begin{array}{c}{e}_e\\ e_f\\ e_g\\ e_h\end{array}\right)=B\left\{\left(\begin{array}{c}{i}_{\mathrm{er}}\\ i_{\mathrm{fr}}\end{array}\right)-\left(\begin{array}{c}{i}_e\\ i_f\end{array}\right)\right\}+A\left(\begin{array}{c}{e}_e\\ e_f\\ e_g\\ e_h\end{array}\right)& \left(30\right)\end{array}$

The general filter of the ef coordinate system represented by Expressions (29) and (30) is transformed into the dq coordinate system. In this transformation, φ2×2 is a zero matrix of two rows and two columns.

$\begin{array}{cc}\begin{array}{c}\left(\begin{array}{c}{v}_d\\ v_q\end{array}\right)=T^{\prime }\left(m\theta \right){\kappa }_{p}T\left(m\theta \right)\left\{\left(\begin{array}{c}{i}_{\mathrm{dr}}\\ i_{\mathrm{qr}}\end{array}\right)-\left(\begin{array}{c}{i}_d\\ i_q\end{array}\right)\right\}+\\ T^{\prime }\left(m\theta \right){\kappa }_i\left(\begin{array}{cc}T\left(m\theta \right)& {\varphi }_{2\times 2}\\ {\varphi }_{2\times 2}& T\left(m\theta \right)\end{array}\right)\left(\begin{array}{c}{e}_d\\ e_q\\ e_r\\ e_s\end{array}\right)\end{array}& \left(31\right)\\ \begin{array}{c}\frac{d}{dt}\left(\begin{array}{c}{e}_d\\ e_q\\ e_r\\ e_s\end{array}\right)=\frac{d}{dt}\left\{\left(\begin{array}{cc}{T}^{\prime }\left(m\theta \right)& {\varphi }_{2\times 2}\\ {\varphi }_{2\times 2}& T^{\prime }\left(m\theta \right)\end{array}\right)\left(\begin{array}{c}{e}_e\\ e_f\\ e_g\\ e_h\end{array}\right)\right\}\\ =\frac{d}{dt}\left\{\left(\begin{array}{cc}{T}^{\prime }\left(m\theta \right)& {\varphi }_{2\times 2}\\ {\varphi }_{2\times 2}& T^{\prime }\left(m\theta \right)\end{array}\right)\right\}\left(\begin{array}{cc}T\left(m\theta \right)& {\varphi }_{2\times 2}\\ {\varphi }_{2\times 2}& T\left(m\theta \right)\end{array}\right)\left(\begin{array}{c}{e}_d\\ e_q\\ e_r\\ e_s\end{array}\right)+\\ \left(\begin{array}{cc}{T}^{\prime }\left(m\theta \right)& {\varphi }_{2\times 2}\\ {\varphi }_{2\times 2}& T^{\prime }\left(m\theta \right)\end{array}\right)\frac{d}{dt}\left(\begin{array}{c}{e}_e\\ e_f\\ e_g\\ e_h\end{array}\right)\\ =\frac{d}{dt}\left\{\left(\begin{array}{cc}{T}^{\prime }\left(m\theta \right)& {\varphi }_{2\times 2}\\ {\varphi }_{2\times 2}& T^{\prime }\left(m\theta \right)\end{array}\right)\right\}\left(\begin{array}{cc}T\left(m\theta \right)& {\varphi }_{2\times 2}\\ {\varphi }_{2\times 2}& T\left(m\theta \right)\end{array}\right)\left(\begin{array}{c}{e}_d\\ e_q\\ e_r\\ e_s\end{array}\right)+\\ \left(\begin{array}{cc}{T}^{\prime }\left(m\theta \right)& {\varphi }_{2\times 2}\\ {\varphi }_{2\times 2}& T^{\prime }\left(m\theta \right)\end{array}\right)\left[B\left\{\left(\begin{array}{c}{i}_{\mathrm{er}}\\ i_{\mathrm{fr}}\end{array}\right)-\left(\begin{array}{c}{i}_e\\ i_f\end{array}\right)\right\}+A\left(\begin{array}{c}{e}_e\\ e_f\\ e_g\\ e_h\end{array}\right)\right]\\ =\frac{d\theta }{dt}\left(\begin{array}{cc}\left(\begin{array}{cc}0& -1\\ 1& 0\end{array}\right)& {\varphi }_{2\times 2}\\ {\varphi }_{2\times 2}& \left(\begin{array}{cc}0& -1\\ 1& 0\end{array}\right)\end{array}\right)\left(\begin{array}{c}{e}_d\\ e_q\\ e_r\\ e_s\end{array}\right)+\\ \left(\begin{array}{cc}{T}^{\prime }\left(m\theta \right)& {\varphi }_{2\times 2}\\ {\varphi }_{2\times 2}& T^{\prime }\left(m\theta \right)\end{array}\right)\mathrm{BT}\left(m\theta \right)\left\{\left(\begin{array}{c}{i}_{\mathrm{dr}}\\ i_{\mathrm{qr}}\end{array}\right)-\left(\begin{array}{c}{i}_d\\ i_q\end{array}\right)\right\}+\\ \left(\begin{array}{cc}{T}^{\prime }\left(m\theta \right)& {\varphi }_{2\times 2}\\ {\varphi }_{2\times 2}& T^{\prime }\left(m\theta \right)\end{array}\right)A\left(\begin{array}{cc}T\left(m\theta \right)& {\varphi }_{2\times 2}\\ {\varphi }_{2\times 2}& T\left(m\theta \right)\end{array}\right)\left(\begin{array}{c}{e}_d\\ e_q\\ e_r\\ e_s\end{array}\right)\end{array}& \left(32\right)\end{array}$

In Expression (32), the first term in the right side member is a part corresponding to an interference of the integration term, and the second and third terms are a part corresponding to the regular filter. That is, in case of the above-described PI control, this corresponds to a part where (n−1)ω is multiplied and a result is added to or subtracted from an error of the other axis. As described above, in this example, since the first term in the right side member exists, the appropriate control can be achieved in cases where the nth-order filter is realized in the dq axis coordinate system.

Moreover, in this example, the filter on the ef axis coordinate system is realized on the dq axis coordinate system. However, the filter on one coordinate system can be realized with another coordinate by the same technique.

“Simulation Example”

A simulation result when the sixth-order higher harmonic wave is applied to a back electromotive force in a motor voltage equation represented by the dq axes will now be described with reference to FIGS. 6 to 8. In these drawings, the top portion represents an axis current, the lower portion represents an axis current, and the right side is an expanded view of the left side.

FIG. 6 shows a conventional technique by which the PI control is executed with respect to a fundamental wave in the dq axes, and it can be understood that a higher harmonic component cannot be suppressed. FIG. 7 shows a technique by which the dq axes are prepared in accordance with each higher harmonic wave, and FIG. 8 shows a technique according to this embodiment which also performs the control over a higher harmonic wave having a frequency six times that of a rotor as the control of the dq axis coordinate system, in which the control technique itself is transformed into the dq axis coordinate system. As described above, it can be understood that the technique according to this embodiment or the method which prepares the coordinate system in accordance with each higher harmonic wave can suppress a harmonic wave six times greater than possible with the conventional method.

As described above, according to this embodiment, the PI control in the ef coordinate system which realizes a rotational frequency which is n times that of the rotor can be carried out for control in a dq axis coordinate system. Further, the nth-order filter in one coordinate system can be handled as a filter on any other coordinate system which relatively rotates with respect to the former coordinate. Therefore, in the motor driving control, the control taking a higher harmonic wave or the like into consideration can be achieved without performing coordinate transformation of all variables each time.

Claims

1. A method of controlling an alternating-current motor including a rotor and a stator, wherein in one coordinate system of arbitrary rectangular coordinate systems including a coordinate system fixed to the stator, a coordinate system fixed to the rotor, and a coordinate system which rotates at a rotational frequency which is n times (n is an integer which is not 0 or 1) that of the rotor, defining a filter on any other coordinate system, and controlling driving of the motor by utilizing the filter.

2. A method of controlling an alternating-current motor including a rotor and a stator, the method comprising: obtaining via a first device a current control output which controls a fundamental wave component of a motor current in an axis coordinate system which rotates in synchronization with the rotor with respect to a d axis current which is mainly an exciting current component of a phase current of the alternating-current motor and a q axis current which is mainly a torque current component of the same;controlling via a second device an e axis current and an f axis current, the e axis current and the f axis current having a frequency which is an integer multiple of a frequency of the fundamental wave component of the alternating-current motor, and realizing via the second device the e and f axes current controls on the d axis current control and the q axis current control; andadding an output from the first device and an output from the second device, and outputting a motor driving current,wherein PI control in the second device is also executed as current control on the d and q axes.

3. The method according to claim 2, wherein the step of controlling via the second device comprises: performing via a third device subtraction or addition with respect to a constant multiple of an output from an f axis current error integration controller and an e axis current error;performing via a fourth device addition or subtraction with respect to a constant multiple of an output from an e axis current error integration controller and the e axis current error;subjecting an output from the third device to proportional multiplication via an e axis current error proportional controller;subjecting an output from the fourth device to proportional multiplication via an f axis current error proportional controller;adding a constant multiple of the e axis current error and an output from the e axis current error proportional controller;adding a constant multiple of the f axis current error and an output from the f axis current error proportional controller;subjecting an output from the third device to integration via the e axis current error integration controller; andsubjecting an output from the fourth device to integration and proportional multiplication via the f axis current error integration controller.

4. The method according to claim 3, wherein the constants are the same value and the value is a product of n−1 and a rotation angular velocity based on an electrical angle, wherein n is a rotation order of the e and f axes which rotate at a rotation frequency which is n times that of the d and q axes, wherein n is an integer which is not 0 or 1.

5. The method according to claim 2, wherein the step of realizing via the second device further comprises, realizing on the d and q axes PI control on the e and f axes by utilizing expressions obtained by transforming expressions of the PI control on the e and f axes into expressions on the d and q axes.

6. A method of controlling current by separating an α phase current and a β phase current as current components of axes which are fixed to a stator and orthogonal to each other with respect to an alternating-current motor including a rotor and the stator, and controlling the respective phase currents based on an α phase current error and a β phase current error which are errors from command values of motor currents of the α phase current and the β phase current, the method comprising: multiplying the α phase current error by a predetermined proportional coefficient to obtain an α phase proportional control component;integrating via a first device a difference value between the α phase current error and a coefficient multiple based on a rotational frequency of the rotor as an output from a second device;integrating via the second device an added value of the β phase current error and a coefficient multiple based on a rotational frequency as an output from the first device;adding an α phase integration control component obtained by multiplying an output from the first device by a predetermined proportional coefficient and the α phase proportional control component, and outputting an α phase voltage command;adding a β phase integration control component obtained by multiplying an output from the second device by a predetermined proportional coefficient and a β phase proportional control component as an output from a β phase proportional unit, and outputting a β phase voltage command; anddriving the alternating-current motor based on the α phase voltage command and the β phase voltage command.