(57)【特許請求の範囲】
【請求項１】
電動機の固定子上に設定された直交座標系でのベクトルの向きが前回生成した重畳電圧指令に対して９０度異なる重畳電圧指令を所定周期で生成する重畳成分生成部と、
前記重畳電圧指令を重畳した駆動電圧指令に基づいた駆動電圧を前記電動機へ出力するインバータ部と、
前記電動機の各相に流れる電流を前記所定周期で検出して、電流検出値を出力する電流検出部と、
前記電流検出値の変化量に基づき、前記電動機の磁極位置を検出する磁極位置推定部と
を備える電動機の制御装置。
【請求項２】
前記重畳成分生成部は、
前記前回生成した重畳電圧指令に対して前記直交座標系でのベクトルの向きが正方向または逆方向に９０度ずつ異なるように前記重畳電圧指令を繰り返し生成する
ことを特徴とする請求項１に記載の電動機の制御装置。
【請求項３】
前記直交座標系の２軸をα軸およびβ軸とした直交座標軸とし、
前記重畳電圧指令のベクトルの向きは、前記α軸または前記β軸に対して平行である
ことを特徴とする請求項２に記載の電動機の制御装置。
【請求項４】
前記直交座標系の軸をα軸およびβ軸とした直交座標軸とし、
前記重畳電圧指令のベクトルの向きは、前記α軸または前記β軸に対して４５度の位相差をもつ
ことを特徴とする請求項２に記載の電動機の制御装置。
【請求項５】
前記電流検出部は、
前記電動機の各相に流れる電流の値を前記直交座標系の電流検出値に変換して出力し、
前記磁極位置推定部は、
前記電流検出部から出力される電流検出値の１階差分値を演算する電流差分演算部と、
前記１階差分値に基づき、前記電動機の磁極位置を求める磁極位置演算部と
を備えることを特徴とする請求項３に記載の電動機の制御装置。 
【請求項６】
前記電流検出部は、
前記電動機の各相に流れる電流の値を前記直交座標系の電流検出値に変換して出力し、
前記磁極位置推定部は、
前記電流検出部から出力される電流検出値の２階差分値を演算する電流差分演算部と、
前記２階差分値に基づき、前記電動機の磁極位置を求める磁極位置演算部と
を備えることを特徴とする請求項４に記載の電動機の制御装置。 
【請求項７】
前記重畳電圧指令を重畳した駆動電圧指令に基づき、前記駆動電圧指令毎に、空間ベクトル変調法によって時間周期内で異なる複数の駆動信号を生成して前記インバータ部へ出力する駆動信号生成部を備え、
前記電流検出部は、
前記時間周期のｎＴ／２（ｎは自然数）で前記電流検出値を前記磁極位置推定部へ出力する
ことを特徴とする請求項１〜６のいずれか１項に記載の電動機の制御装置。 
【請求項８】
前記重畳電圧指令を重畳した駆動電圧指令と搬送波との比較に基づく搬送波比較変調法によって、駆動信号を生成して前記インバータ部へ出力する駆動信号生成部を備え、
前記電流検出部は、
前記電流検出値を前記三角波状の搬送波の山および／または谷のタイミングで前記磁極位置推定部へ出力する
ことを特徴とする請求項１〜６のいずれか１項に記載の電動機の制御装置。 
【請求項９】
前記電動機は、
埋め込み型永久磁石同期電動機である
ことを特徴とする請求項１〜８のいずれか１項に記載の電動機の制御装置。
【請求項１０】
電動機の固定子上に設定された直交座標系でのベクトルの向きが前回生成した重畳電圧指令に対して９０度異なる重畳電圧指令を所定周期で生成する重畳成分生成部と、
前記重畳電圧指令を重畳した駆動電圧指令に基づいて駆動される電動機の各相に流れる電流を前記所定周期で検出して、電流検出値を出力する電流検出部と、
前記電流検出値の変化量に基づき、前記電動機の磁極位置を検出する磁極位置推定部と
を備える電動機の磁極位置推定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
開示の実施形態は、電動機の磁極位置推定装置およびそれを用いた制御装置に関する。
【背景技術】
【０００２】
従来、突極性を有する電動機の磁極位置を推定する磁極位置推定装置が知られており、かかる磁極位置推定装置として、高調波成分をＰＷＭ信号に重畳させることによって、電動機の磁極位置を推定するものが知られている。【０００３】
例えば、特許文献１に記載の磁極位置推定装置では、ＰＷＭ信号を得る際に、三相巻線の各相に、搬送波の連続する３周期分を１期間とし、１／３の期間（最初の搬送波１周期）で本来の指令値を３倍し、残りの２／３の期間（残り搬送波２周期）で高周波成分を重畳させて電動機の磁極位置を推定する。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特許第４６７００４４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、特許文献１に記載の磁極位置推定装置では、電動機を駆動するための本来の指令値が搬送波の３周期毎でしか更新できない。【０００６】
実施形態の一態様は、上記に鑑みてなされたものであって、電動機を駆動するための電圧指令の更新周期に影響を与えることなく磁極位置の推定を行うこができる磁極位置推定装置および電動機の制御装置を提供することを目的とする。
【課題を解決するための手段】
【０００７】
実施形態の一態様に係る電動機の制御装置は、重畳成分生成部と、インバータ部と、電流検出部と、磁極位置推定部とを備える。前記重畳成分生成部は、電動機の固定子上に設定された座標系でのベクトルの向きが前回生成した重畳電圧指令に対して９０度異なる重畳電圧指令を所定周期で生成する。前記インバータ部は、前記重畳電圧指令を重畳した駆動電圧指令に基づいた駆動電圧を前記電動機へ出力する。前記電流検出部は、前記電動機の各相に流れる電流を前記所定周期で検出して、電流検出値を出力する。前記磁極位置推定部は、前記電流検出値の変化量に基づき、前記電動機の磁極位置を検出する。
【発明の効果】
【０００８】
実施形態の一態様によれば、電動機を駆動するための電圧指令の更新周期に影響を与えることなく磁極位置の推定を行うことができる電動機の磁極位置推定装置およびそれを用いた制御装置を提供することができる。
【図面の簡単な説明】
【０００９】
【図１】図１は、第１の実施形態に係る電動機の制御装置の構成を示す図である。
【図２Ａ】図２Ａは、空間ベクトル変調法を説明するための図である。
【図２Ｂ】図２Ｂは、空間ベクトル変調法を説明するための図である。
【図３】図３は、αβ軸座標系における電圧ベクトルとφとの関係の一例を示す図である。
【図４】図４は、パイロット電圧の時間周期における変化を示す図である。
【図５】図５は、重畳成分生成部のパイロット電圧生成処理を示すフローチャートである。
【図６】図６は、重畳成分生成部から出力されるパイロット電圧のαβ軸成分の変化を示す図である。
【図７】図７は、磁極位置推定部の構成を示す図である。
【図８】図８は、磁極位置推定部の磁極位置推定処理を示すフローチャートである。
【図９】図９は、パイロット電圧と搬送波との関係を示す図である。
【図１０】図１０は、第２の実施形態に係る制御装置の構成を示す図である。
【図１１】図１１は、第２の実施形態に係る電圧ベクトルとφとの関係の一例を示す図である。
【図１２】図１２は、第２の実施形態に係るパイロット電圧の時間周期における変化を示す図である。
【図１３】図１３は、第２の実施形態に係るパイロット電圧のαβ軸成分の変化を示す図である。
【図１４】図１４は、第２の実施形態に係る重畳成分生成部のパイロット電圧生成処理を示すフローチャートである。
【図１５】図１５は、第２の実施形態に係る磁極位置推定部の構成を示す図である。
【図１６】図１６は、第２の実施形態に係る磁極位置推定部の磁極位置推定処理を示すフローチャートである。
【発明を実施するための形態】
【００１０】
以下、図面を参照して、本願の開示する電動機の磁極位置推定装置およびそれを用いた制御装置の実施形態を詳細に説明する。なお、電動機の制御装置は、磁極位置推定装置を含み、以下、電動機の制御装置を単に「制御装置」と記載する。また、以下に示す実施形態によりこの発明が限定されるものではない。
【００１１】
（第１の実施形態）
まず、第１の実施形態に係る制御装置について説明する。図１は、第１の実施形態に係る制御装置の構成を示す図である。
【００１２】
図１に示すように、第１の実施形態に係る制御装置１は、インバータ部１０と、電流検出部２０と、制御部３０とを備える。かかる制御装置１は、直流電源２と電動機３との間に接続され、直流電源２から供給される直流電圧Ｖｄｃをインバータ部１０によって交流電圧へ変換して電動機３へ出力することによって電動機３を動作させる。
【００１３】
電動機３は、突極性を有する電動機であり、電機子巻線を有する固定子３ａと、複数の永久磁石を回転子コアに埋め込んで周方向に配置した回転子３ｂとを含む。かかる電動機３としては、例えば、埋め込み型永久磁石同期電動機（Interior  Permanent  Magnet  Synchronous  Motor）などがある。
【００１４】
インバータ部１０は、スイッチング素子Ｑ１〜Ｑ６と、ダイオードＤ１〜Ｄ６と、平滑コンデンサＣ１とを備える。スイッチング素子Ｑ１〜Ｑ６は三相ブリッジ接続され、スイッチング素子Ｑ１〜Ｑ６毎にそれぞれダイオードＤ１〜Ｄ６が逆並列に接続される。
【００１５】
かかるインバータ部１０は、制御部３０から出力される駆動信号Ｓ１〜Ｓ６に基づいてスイッチング素子Ｑ１〜Ｑ６が制御され、駆動信号Ｓ１〜Ｓ６に応じた電圧（以下、駆動電圧と記載する）を電動機３へ出力する。なお、スイッチング素子Ｑ１〜Ｑ６として、例えば、ＩＧＢＴやＭＯＳＦＥＴなどの半導体素子が用いられる。
【００１６】
電流検出部２０は、インバータ部１０と電動機３との間のＵ相、Ｖ相およびＷ相の各相に流れる電流を検出し、その検出結果として電流検出値ＩＵ、ＩＶ、ＩＷを出力する。電流検出値ＩＵは、Ｕ相電流の瞬時値であり、電流検出値ＩＶは、Ｖ相電流の瞬時値であり、電流検出値ＩＷは、Ｗ相電流の瞬時値である。なお、電流検出部２０として、例えば、磁電変換素子であるホール素子を利用した電流センサを用いることができる。
【００１７】
制御部３０は、３相／２相座標変換部３１と、回転座標変換部３２と、電流指令出力部３３と、電流制御部３４と、カウンタ部３５と、重畳成分生成部３６と、加算部３７、３８と、駆動信号生成部３９と、磁極位置推定部４０とを備える。制御部３０は、磁極位置推定装置としての機能を備え、電動機３における回転子３ｂの磁極位置θ（以下、電動機３の磁極位置θと記載する場合がある）を推定することができる。
【００１８】
３相／２相座標変換部３１は、電流検出値Ｉｕ、ＩＶ、ＩＷを固定座標上の直交した２軸のαβ軸成分へ変換して、α軸方向の電流検出値Ｉαとβ軸方向の電流検出値Ｉβとをベクトル成分とするαβ軸座標系の電流ベクトルを求める。αβ軸座標系は、電動機３の固定子３ａ上に設定される直交座標系であり、固定子座標系とも呼ばれる。なお、電流検出部２０に、３相／２相座標変換部３１を設けるようにしてもよい。
【００１９】
回転座標変換部３２は、電動機３の磁極位置θに基づき、３相／２相座標変換部３１から出力されるαβ軸座標系の成分を座標変換する。これにより、回転座標変換部３２は、回転子３ｂに対応して回転するｄｑ軸回転座標系のｑ軸成分およびｄ軸成分であるｑ軸電流検出値Ｉｑおよびｄ軸電流検出値Ｉｄを求め、電流制御部３４へ出力する。なお、電動機３の磁極位置θは、αβ軸座標系からみた磁極位置である。
【００２０】
電流指令出力部３３は、外部から入力されるトルク指令Ｔ*に基づき、ｑ軸電流指令Ｉｑ*と、ｄ軸電流指令Ｉｄ*とを生成し、電流制御部３４へ出力する。ｑ軸電流指令Ｉｑ*はｑ軸成分の電流指令であり、ｄ軸電流指令Ｉｄ*は、ｄ軸成分の電流指令である。
【００２１】
電流制御部３４は、ｑ軸電流指令Ｉｑ*、ｄ軸電流指令Ｉｄ*、ｑ軸電流検出値Ｉｑ、ｄ軸電流検出値Ｉｄおよび電動機３の磁極位置θに基づき、αβ軸座標系のα軸指令成分Ｖα*およびβ軸指令成分Ｖβ*を駆動電圧として生成して出力する。
【００２２】
具体的には、電流制御部３４は、ｑ軸電流指令Ｉｑ*とｑ軸電流検出値Ｉｑとの偏差を零とするようにｑ軸電圧指令Ｖｑ*を調整し、ｄ軸電流指令Ｉｄ*とｄ軸電流検出値Ｉｄとの偏差を零とするようにｄ軸電圧指令Ｖｄ*を調整する。さらに、電流制御部３４は、電動機３の磁極位置θに基づき、ｑ軸電圧指令Ｖｑ*およびｄ軸電圧指令Ｖｄ*をαβ軸座標系の成分へ変換して、駆動電圧指令であるα軸指令成分Ｖα*およびβ軸指令成分Ｖβ*を生成する。
【００２３】
カウンタ部３５は、後述する時間周期Ｔを４周期分として、内部のカウント値ＣＮＴをカウントアップまたはリセットし、０〜３のカウント値ＣＮＴを生成する。具体的には、カウンタ部３５は、時刻ｔ＝ｎＴ（ｎ＝０、１、２、３、・・・）の場合にＣＮＴ＝０とし、時刻ｔ＝Ｔ／４＋ｎＴの場合にＣＮＴ＝１とし、時刻ｔ＝Ｔ／２＋ｎＴの場合にＣＮＴ＝２とし、時刻ｔ＝３Ｔ／４＋ｎＴの場合にＣＮＴ＝３とする。カウンタ部３５は、生成したカウント値ＣＮＴを重畳成分生成部３６および磁極位置推定部４０へ出力する。
【００２４】
重畳成分生成部３６は、α軸指令成分Ｖα*およびβ軸指令成分Ｖβ*へ重畳するための重畳電圧指令として電圧ベクトルＶｈ（以下、パイロット電圧Ｖｈと記載することもある）を生成する。具体的には、重畳成分生成部３６は、所定周期毎に、αβ軸座標系でのベクトルの向きが前回生成したパイロット電圧Ｖｈに対して９０度異なるパイロット電圧Ｖｈを生成する。かかるパイロット電圧Ｖｈは、α軸成分であるα軸指令成分Ｖαｈと、β軸成分であるβ軸指令成分Ｖβｈとによって構成される。
【００２５】
加算部３７は、電流制御部３４から入力されるα軸指令成分Ｖα*と重畳成分生成部３６から入力されるα軸指令成分Ｖαｈとを加算して生成したα軸指令成分Ｖα１*を駆動信号生成部３９へ出力する。また、加算部３８は、電流制御部３４から入力されるβ軸指令成分Ｖβ*と重畳成分生成部３６から入力されるβ軸指令成分Ｖβｈとを加算して生成したβ軸指令成分Ｖβ１*を駆動信号生成部３９へ出力する。
【００２６】
駆動信号生成部３９は、加算部３７、３８から出力されるα軸指令成分Ｖα１*およびβ軸指令成分Ｖβ１*に基づき、空間ベクトル変調法を用いてインバータ部１０を駆動する駆動信号Ｓ１〜Ｓ６を生成する。
【００２７】
図２Ａおよび図２Ｂは、空間ベクトル変調法を説明するための図であって、これらの図では、任意に設定した時間周期ＴＳが２ＴＶで表されている。図２Ａには、空間ベクトル変調における電圧ベクトルＶ０〜Ｖ７が示されており、この図に示される例では、α軸指令成分Ｖα１*およびβ軸指令成分Ｖβ１*によって規定される電圧指令ベクトルＶｓ*が２つの隣接する電圧ベクトルＶ１、Ｖ２を用いて形成される。
【００２８】
電圧ベクトルＶ１（１００）は、インバータ部１０においてＵ相の上側スイッチング素子Ｑ１をＯＮにし、Ｕ相の下側スイッチング素子Ｑ４をＯＦＦにし、Ｖ相およびＷ相の上側スイッチング素子Ｑ２、Ｑ３をＯＦＦにし、Ｖ相およびＷ相の下側スイッチング素子Ｑ５、Ｑ６をＯＮにする。一方、電圧ベクトルＶ２（１１０）は、Ｕ相およびＶ相の上側スイッチング素子Ｑ１、Ｑ２をＯＮにし、Ｕ相およびＶ相の下側スイッチング素子Ｑ４、Ｑ５をＯＦＦとし、Ｗ相の上側スイッチング素子Ｑ３をＯＦＦにし、Ｗ相の下側スイッチング素子Ｑ６をＯＮにする。
【００２９】
時間周期Ｔにおける各相のスイッチング状態は図２Ｂに示すようになる。図２Ｂに示すＯＮ時間ｔ１、ｔ２は下記式（１）、（２）により算出される。なお、下記式（１）、（２）において、θＶは電圧ベクトルＶ１から電圧指令ベクトルＶｓ*までの位相角、｜Ｖｓ*｜は電圧指令ベクトルＶｓ*の振幅、Ｖmaxはインバータ部１０の出力最大電圧である。
【００３０】
【数１】
0005900600.tif 000002
【００３１】
駆動信号生成部３９において、電圧指令ベクトルＶｓ*の更新は、任意に設定した時間周期ＴＳの半周期ＴＶ毎に行われ、ＯＦＦ時間ｔ０、ｔ７は、半周期ＴＶからＯＮ時間ｔ１、ｔ２を引いた残りの時間を分割して設定される。その他の象限も同様に、電圧指令ベクトルＶｓ*が設定される。
【００３２】
磁極位置推定部４０は、３相／２相座標変換部３１から入力される電流検出値Ｉα、Ｉβの変化量に基づき、電動機３の磁極位置θを推定し、かかる推定結果を回転座標変換部３２および電流制御部３４へ出力する。
【００３３】
制御装置１は、重畳成分生成部３６が所定周期で生成する９０度ずつ位相の異なるパイロット電圧Ｖｈと、電動機３の各相に流れる電流を用いて、電動機３を駆動するための駆動電圧指令の更新周期に影響を与えることなく、電動機３の磁極位置θを容易に推定できるようにしている。以下、磁極位置θの推定についてさらに具体的に説明する。
【００３４】
まず、突極性を有する電動機３に関し、αβ軸座標系における電流・電圧方程式は、下記式（３）のように表すことができる。なお、Ｖα、Ｖβは電動機３の電圧のαβ軸座標系における各成分を示す。Ｉα、Ｉβは電動機３に流れる電流のαβ軸座標系における各成分を示す。Ｒは電動機３の電機子抵抗を、Ｌおよびｌは電動機３の電機子反作用のインダクタンスを示す。また、Ｋeは電動機３の誘起電圧定数を、ωrは磁極回転速度を示す。
【００３５】
【数２】
0005900600.tif 000003
【００３６】
ここで、インバータ部１０の出力電圧として、駆動電圧指令であるｑ軸電圧指令Ｖｑ*およびｄ軸電圧指令Ｖｄ*に基づく駆動電圧に、パイロット電圧Ｖｈが含まれるようにすることを考える。例えば、図３に示すように、αβ軸座標系において、電圧ベクトルＶｈ０を重畳し、一定時間後に電圧ベクトルＶｈ１を重畳する。図３は、αβ軸座標系における電圧ベクトルＶｈ０、Ｖｈ１とφとの関係の一例を示す図である。
【００３７】
電圧ベクトルＶｈ０は、α軸の正の向きに位相差φ［ｄｅｇ］で、かつ、振幅Ｖinj1（＞０）のベクトルであり、下記式（４）に示すように表される。また、電圧ベクトルＶｈ１は、α軸の正の向きに位相差φ＋９０［ｄｅｇ］で、かつ、振幅Ｖinj1（＞０）のベクトルであり、下記式（５）に示すように表される。
【００３８】
【数３】
0005900600.tif 000004
【００３９】
重畳する２つの電圧ベクトル間の時間間隔が、駆動電圧周期に比べて十分に短い場合、上記式（３）において駆動電圧周波数の成分を取り除いた電流・電圧方程式は、下記式（６）のように近似できる。
【００４０】
【数４】
0005900600.tif 000005
【００４１】
さらに、電動機３の回転数が低速または中速の場合、電動機３の誘起電圧を無視することができ、上記式（６）は下記式（７）のように近似できる。
【００４２】
【数５】
0005900600.tif 000006
【００４３】
上記式（７）は、電流に関して下記式（８）に示すように書き直すことができる。
【００４４】
【数６】
0005900600.tif 000007
【００４５】
ここで、上記式（４）で示される電圧ベクトルＶｈ０の重畳時刻を時刻ｔｈ０とし、上記式（５）で示される電圧ベクトルＶｈ１の重畳時刻を時刻ｔｈ１とすると、時刻ｔｈ０と時刻ｔｈ１での電流・電圧方程式はそれぞれ下記式（９）、（１０）に示すように表すことができる。
【００４６】
【数７】
0005900600.tif 000008
【００４７】
さらに、上記式（９）、（１０）から、下記式（１１）の関係式を導くことができる。
【００４８】
【数８】
0005900600.tif 000009
【００４９】
上記式（１１）は、電動機３の磁極位置θに関して下記式（１２）に示すように書き直すことができ、したがって、時刻ｔｈ０と時刻ｔｈ１での各相の電流変化を検出すれば電動機３の磁極位置θを推定することができる。
【００５０】
【数９】
0005900600.tif 000010
【００５１】
ここで、φを０、９０、１８０、２７０［ｄｅｇ］のいずれかに設定すると、電動機３の磁極位置θをより簡単な演算式で推定することができる。すなわち、上記式（１２）は、φ＝０［ｄｅｇ］の場合、下記式（１３）のように、φ＝９０［ｄｅｇ］の場合、下記式（１４）のように、φ＝１８０［ｄｅｇ］の場合、下記式（１５）のように、φ＝２７０［ｄｅｇ］の場合、下記式（１６）のようにそれぞれ表すことができる。
【００５２】
【数１０】
0005900600.tif 000011
【００５３】
したがって、電圧ベクトルＶｈ０を、φ＝０、９０、１８０、２７０［ｄｅｇ］のいずれかの状態になるように重畳し、その後、電圧ベクトルＶｈ１を重畳することによって、電動機３の磁極位置θをより簡単な演算式で推定できる。
【００５４】
また、時刻ｔｈ１でα軸の正の向きに位相差φ＋９０［ｄｅｇ］の電圧ベクトルＶｈ１を重畳する場合を説明したが、時刻ｔｈ１でα軸の正の向きに位相差φ−９０［ｄｅｇ］の電圧ベクトルＶｈ１を重畳しても同様に、電動機３の磁極位置θをより簡単な演算式で推定できる。
【００５５】
このように、時刻ｔｈ０においてφを０、９０、１８０、２７０［ｄｅｇ］のいずれかの位相に電圧ベクトルＶｈ０を重畳した後、時刻ｔｈ１でφ＋９０［ｄｅｇ］またはφ−９０［ｄｅｇ］の位相に電圧ベクトルＶｈ１を重畳する。そして、時刻ｔｈ０と時刻ｔｈ１での各相の電流変化を検出し、例えば、上記式（１３）〜（１６）に基づいた演算を用いることで、電動機３の磁極位置θをより簡単に推定できる。
【００５６】
そこで、重畳成分生成部３６は、下記式（１７）に示すように、時間周期Ｔ＝４ＴＶのパイロット電圧Ｖｈを生成し、α軸指令成分Ｖα*およびβ軸指令成分Ｖβ*で規定される駆動電圧指令に重畳するようにしている。ＴＶは、空間ベクトル変調を行う駆動信号生成部３９の更新周期である。
【００５７】
【数１１】
0005900600.tif 000012
【００５８】
上記式（１７）に示すように重畳成分生成部３６から出力するパイロット電圧Ｖｈのα軸指令成分Ｖαｈおよびβ軸指令成分Ｖβｈを設定することで、図４に示すように、α軸およびβ軸で打ち消し合う方向にパイロット電圧Ｖｈを出力できる。図４は、パイロット電圧Ｖｈの時間周期Ｔにおける変化を示す図である。そのため、時間周期Ｔにおける平均電圧を０とすることができ、これにより、トルクリップルを低減することができる。
【００５９】
ここで、重畳成分生成部３６の構成および動作について図５および図６を参照してさらに具体的に説明する。図５は、重畳成分生成部３６から出力されるパイロット電圧Ｖｈのαβ軸成分であるＶαｈ、Ｖβｈの変化を示す図であり、図６は、重畳成分生成部３６のパイロット電圧生成処理を示すフローチャートである。
【００６０】
図５に示すように、重畳成分生成部３６は、カウント値ＣＮＴとして０、１、２、３の順に繰り返し出力するカウンタ部３５のカウント値ＣＮＴに基づき、パイロット電圧Ｖｈのα軸成分およびβ軸成分であるα軸指令成分Ｖαｈおよびβ軸指令成分Ｖβｈを生成する。
【００６１】
具体的には、図６に示すように、重畳成分生成部３６は、カウンタ部３５のカウント値ＣＮＴを判定し（ステップＳ１０）、カウント値ＣＮＴ＝０の場合、Ｖαｈ＝＋Ｖｉｎｊ１およびＶβｈ＝０としたパイロット電圧Ｖｈを生成して、加算部３７、３８へ出力する（ステップＳ１１）。カウント値ＣＮＴ＝１の場合、重畳成分生成部３６は、Ｖαｈ＝０およびＶβｈ＝＋Ｖｉｎｊ１としたパイロット電圧Ｖｈを生成して、加算部３７、３８へ出力する（ステップＳ１２）。
【００６２】
また、カウント値ＣＮＴ＝２の場合、重畳成分生成部３６は、Ｖαｈ＝−Ｖｉｎｊ１およびＶβｈ＝０としたパイロット電圧Ｖｈを生成して、加算部３７、３８へ出力する（ステップＳ１３）。カウント値ＣＮＴ＝３の場合、重畳成分生成部３６は、Ｖαｈ＝０およびＶβｈ＝−Ｖｉｎｊ１としたパイロット電圧Ｖｈを生成して、加算部３７、３８へ出力する（ステップＳ１４）。
【００６３】
このように、重畳成分生成部３６は、αβ軸座標系で、時間周期Ｔｖ毎に前回生成したパイロット電圧Ｖｈに対して９０度の方向となるパイロット電圧Ｖｈを生成するようにしている。なお、重畳成分生成部３６は、αβ軸座標系で、時間周期Ｔｖ毎に前回生成したパイロット電圧Ｖｈに対して−９０度の方向となるパイロット電圧Ｖｈを生成することもできる。
【００６４】
磁極位置推定部４０は、時刻ｔ＝ｎＴ（ｎ＝１、２、３、・・・）での磁極位置θ、時刻ｔ＝Ｔ／４＋ｎＴでの磁極位置θ、時刻ｔ＝Ｔ／２＋ｎＴでの磁極位置θ、時刻ｔ＝３Ｔ／４＋ｎＴでの磁極位置θを推定する。なお、磁極位置推定部４０は、重畳成分生成部３６によって生成されたパイロット電圧Ｖｈに応じた電流を検出することができるように、カウンタ部３５のカウント値ＣＮＴが変化してから所定時間経過後に、カウント値ＣＮＴの判定処理を行う。
【００６５】
図７は、磁極位置推定部４０の構成を示す図である。図７に示すように、磁極位置推定部４０は、差分演算部５０、５１と、１サンプル遅延部５２、５３と、θ演算部５４とを備え、Ｔ／４をサンプル周期として動作する。
【００６６】
差分演算部５０は、３相／２相座標変換部３１からα軸成分電流検出値Ｉαを取得し、かかる電流検出値Ｉαと１サンプル時間（＝Ｔ／４）前に取得した電流検出値Ｉαとの差分値ΔＩαを演算して出力する。また、差分演算部５１は、３相／２相座標変換部３１からβ軸成分電流検出値Ｉβを取得し、かかる電流検出値Ｉβと１サンプル時間前に取得した電流検出値Ｉβとの差分値ΔＩβを演算して出力する。
【００６７】
１サンプル遅延部５２は、差分演算部５０から出力される差分値ΔＩαを１サンプル時間分遅延させて、θ演算部５４へ出力する。また、１サンプル遅延部５３は、差分演算部５１から出力される差分値ΔＩβを１サンプル時間分遅延させて、θ演算部５４へ出力する。
【００６８】
θ演算部５４は、差分演算部５０、５１から出力される差分値ΔＩα、ΔＩβと、１サンプル遅延部５２、５３から出力される１サンプル時間前の差分値ΔＩα、ΔＩβとに基づいて、磁極位置θを推定する。
【００６９】
図８は、磁極位置推定部４０の磁極位置推定処理を示すフローチャートである。以下において、ΔＩαｈ２＝Ｉαｈ２−Ｉαｈ１、ΔＩβｈ２＝Ｉβｈ２−Ｉβｈ１、ΔＩαｈ１＝Ｉαｈ１−Ｉαｈ０、ΔＩβｈ１＝Ｉβｈ１−Ｉβｈ０である。
【００７０】
図８に示すように、磁極位置推定部４０は、カウンタ部３５のカウント値ＣＮＴを判定する（ステップＳ２０）。
【００７１】
カウント値ＣＮＴ＝０の場合、磁極位置推定部４０は、下記式（１８）に基づいて、磁極位置θを演算する（ステップＳ２１）。下記式（１８）は、上記式（１３）において電流１階微分を電流差分で近似したものである。
【００７２】
【数１２】
0005900600.tif 000013
【００７３】
上記式（１８）において、電流検出値Ｉαｈ２、Ｉαｈ１、Ｉαｈ０は、それぞれ時刻ｔ＝ｎＴ、−Ｔ／４＋ｎＴ、−Ｔ／２＋ｎＴにおけるα軸成分の電流検出値である。また、電流検出値Ｉβｈ２、Ｉβｈ１、Ｉβｈ０は、それぞれ時刻ｔ＝ｎＴ、−Ｔ／４＋ｎＴ、−Ｔ／２＋ｎＴにおけるβ軸成分の電流検出値である。
【００７４】
また、カウント値ＣＮＴ＝１の場合、磁極位置推定部４０は、下記式（１９）に基づいて、磁極位置θを演算する（ステップＳ２２）。下記式（１９）は、上記式（１４）において電流１階微分を電流差分で近似したものである。
【００７５】
【数１３】
0005900600.tif 000014
【００７６】
なお、カウント値ＣＮＴ＝０の場合、電流検出値Ｉαｈ２、Ｉαｈ１、Ｉαｈ０は、それぞれ時刻ｔ＝Ｔ／４＋ｎＴ、ｎＴ、−Ｔ／４＋ｎＴにおけるα軸成分の電流検出値である。また、電流検出値Ｉβｈ２、Ｉβｈ１、Ｉβｈ０は、それぞれ時刻ｔ＝Ｔ／４＋ｎＴ、ｎＴ、−Ｔ／４＋ｎＴにおけるβ軸成分の電流検出値である。
【００７７】
また、カウント値ＣＮＴ＝２の場合、磁極位置推定部４０は、下記式（２０）に基づいて、磁極位置θを演算する（ステップＳ２３）。下記式（２０）は、上記式（１５）において電流１階微分を電流差分で近似したものである。
【００７８】
【数１４】
0005900600.tif 000015
【００７９】
なお、カウント値ＣＮＴ＝２の場合、電流検出値Ｉαｈ２、Ｉαｈ１、Ｉαｈ０は、それぞれ時刻ｔ＝Ｔ／２＋ｎＴ、Ｔ／４＋ｎＴ、ｎＴにおけるα軸成分の電流検出値である。また、電流検出値Ｉβｈ２、Ｉβｈ１、Ｉβｈ０は、それぞれ時刻ｔ＝Ｔ／２＋ｎＴ、Ｔ／４＋ｎＴ、ｎＴにおけるβ軸成分の電流検出値である。
【００８０】
また、カウント値ＣＮＴ＝３の場合、磁極位置推定部４０は、下記式（２１）に基づいて、磁極位置θを演算する（ステップＳ２４）。下記式（２１）は、上記式（１６）において電流１階微分を電流差分で近似したものである。
【００８１】
【数１５】
0005900600.tif 000016
【００８２】
なお、カウント値ＣＮＴ＝３の場合、電流検出値Ｉαｈ２、Ｉαｈ１、Ｉαｈ０は、それぞれ時刻ｔ＝３Ｔ／４＋ｎＴ、Ｔ／２＋ｎＴ、Ｔ／４＋ｎＴにおけるα軸成分の電流検出値である。また、電流検出値Ｉβｈ２、Ｉβｈ１、Ｉβｈ０は、それぞれ時刻ｔ＝３Ｔ／４＋ｎＴ、Ｔ／２＋ｎＴ、Ｔ／４＋ｎＴにおけるβ軸成分の電流検出値である。
【００８３】
以上のように、第１の実施形態に係る制御装置１では、重畳成分生成部３６と、磁極位置推定部４０とを備える。重畳成分生成部３６は、αβ軸座標系のα軸またはβ軸に対して平行な方向のベクトルとしてパイロット電圧Ｖｈを生成し、かつ、ベクトルの向きが前回生成したパイロット電圧Ｖｈに対して９０度異なるパイロット電圧Ｖｈを所定周期で生成する。
【００８４】
重畳成分生成部３６で生成したパイロット電圧Ｖｈは重畳電圧指令として駆動電圧指令に重畳されてインバータ部１０へ入力される。これにより、インバータ部１０は、重畳電圧指令が重畳された駆動電圧指令に基づいて電動機３に対して駆動電圧を出力する。磁極位置推定部４０は、電動機３の各相に流れる電流を所定周期で検出し、各相の電流変化量に基づき、電動機３の磁極位置θを推定する。
【００８５】
このように、制御装置１では、所望のパイロット電圧Ｖｈを所定周期毎に生成し、電動機３の各相の電流変化量に基づき、磁極位置θの推定を行うことができる。
【００８６】
したがって、第１の実施形態に係る制御装置１では、磁極位置θの推定を容易に行うことが可能となり、また、電動機３を駆動するための駆動電圧指令の更新周期に影響を与えることもない。しかも、パイロット電圧Ｖｈの重畳を開始してから時間周期Ｔの３／４の時間電動機３の磁極位置θを推定することができ、また、時間周期Ｔの１／４の時間で繰り返し、電動機３の磁極位置θを推定することができるため、磁極位置推定の応答性を向上させることができる。
【００８７】
上述においては、駆動信号生成部３９において空間ベクトル変調法を用いてインバータ部１０を駆動する駆動信号Ｓ１〜Ｓ６を生成する例を説明したが、駆動信号生成部３９において搬送波比較変調法（ＰＷＭ変調法）を用いて駆動信号Ｓ１〜Ｓ６を生成してもよい。
【００８８】
この場合、カウンタ部３５は、時間周期Ｔを搬送波の時間周期Ｔｃの２倍の長さとし、時間周期Ｔを４周期分として、カウント値ＣＮＴを変更する。カウンタ部３５がカウント値ＣＮＴを変更するタイミングは、図９に示すように、搬送波の山および谷である。図９は、駆動信号生成部３９において搬送波比較変調法を用いて駆動信号Ｓ１〜Ｓ６を生成した場合のパイロット電圧Ｖｈ（α軸指令成分Ｖαｈおよびβ軸指令成分Ｖβｈ）と搬送波との関係を示す図である。
【００８９】
重畳成分生成部３６は、カウンタ部３５からのカウント値ＣＮＴに基づき、α軸指令成分Ｖα*およびβ軸指令成分Ｖβ*へ重畳する時間周期Ｔ＝２Ｔｃを４周期分とするパイロット電圧Ｖｈを生成する。カウント値ＣＮＴは、搬送波の山および谷で変更されることから、パイロット電圧Ｖｈのα軸指令成分Ｖαｈおよびβ軸指令成分Ｖβｈは、図５に示すフローチャートに従って搬送波の山および谷で変更される。
【００９０】
駆動信号生成部３９は、加算部３７、３８から入力されるα軸指令成分Ｖα１*およびβ軸指令成分Ｖβ１*に基づき、出力電圧指令ＶＵ*、ＶＶ*、ＶＷ*を求める。そして、駆動信号生成部３９は、出力電圧指令ＶＵ*、ＶＶ*、ＶＷ*と時間周期Ｔｃの搬送波とを比較してＰＷＭ信号である駆動信号Ｓ１〜Ｓ６を生成し、インバータ部１０へ出力する。
【００９１】
磁極位置推定部４０は、カウンタ部３５から入力されるカウント値ＣＮＴに基づき、図８に示すフローチャートに従って時刻ｔ＝ｎＴ、Ｔ／４＋ｎＴ、Ｔ／２＋ｎＴ、３Ｔ／４＋ｎＴでのそれぞれの電動機３の磁極位置θを推定する。
【００９２】
このように、駆動信号生成部３９において搬送波比較変調法を用いて駆動信号Ｓ１〜Ｓ６を生成した場合であっても、空間ベクトル変調法を用いて駆動信号Ｓ１〜Ｓ６を生成した場合と同様に、電動機３の磁極位置θを容易に推定することができる。
【００９３】
なお、上述においては、空間ベクトル変調法を用いた場合には、時間周期Ｔ＝４ＴＶとし、搬送波比較変調法を用いた場合には、時間周期Ｔ＝２Ｔｃとしたが、時間周期Ｔの長さは、かかる例に限定されるものではない。例えば、空間ベクトル変調法を用いた場合に、時間周期Ｔ＝４ｍＴＶ（ｍは２以上の整数）としたり、搬送波比較変調法を用いた場合には、時間周期Ｔ＝２ｍＴｃ（ｍは２以上の整数）とすることができる。
【００９４】
また、上述においては、重畳成分生成部３６は、正方向に位相差が９０度ずつ異なるようにパイロット電圧Ｖｈを所定周期で生成するようにしたが、重畳成分生成部３６によるパイロット電圧Ｖｈの生成方法はこれに限られない。例えば、重畳成分生成部３６は、αβ軸座標系でのベクトルの向きが負方向に９０度ずつの位相差をもつようにパイロット電圧Ｖｈを所定周期で生成するようにしてもよく、これによっても、トルクリップルを低減でき、また、電動機３の磁極位置θをより簡単な演算式で推定できる。
【００９５】
また、重畳成分生成部３６は、位相差９０度とするパイロット電圧Ｖｈの向きを、周期ｍＴ（ｍは自然数）毎に正方向と負方向とで逆転させるようにしてもよく、これによっても、トルクリップルを低減することができる。
【００９６】
また、重畳成分生成部３６は、位相差を正方向と負方向で交互に異なるようにパイロット電圧Ｖｈを生成することもできる。なお、この場合、トルクリップルの低減効果は少ないが、磁極位置推定部４０における磁極位置θの推定を、上記演算式に基づいて行うことができるため、磁極位置θの推定を容易に行うことができる。
【００９７】
また、第１の実施形態に係る制御装置１では、電動機３の回転数が低速または中速の場合、電動機３の誘起電圧を無視することができることから、上記式（６）に近似する上記式（７）を用いている。しかし、電動機３の回転数が高速な場合であっても、磁極位置推定部４０において、電動機３の誘起電圧成分を除去するフィルタを設けることによって、同様に磁極位置θの推定を行うことができる。
【００９８】
また、上述においては、カウンタ部３５からカウント値ＣＮＴを重畳成分生成部３６および磁極位置推定部４０へ出力するようにしたが、重畳成分生成部３６および磁極位置推定部４０にそれぞれカウント値ＣＮＴを生成するカウンタを設けてもよい。
【００９９】
また、上述においては、インバータ部１０の一例として２レベルインバータの構成（図１）を説明したが、インバータ部１０として３レベルインバータなどのマルチレベルインバータやマトリクスコンバータを用いることもでき、また、その他種々の変更が可能である。
【０１００】
（第２の実施形態）
次に、第２の実施形態に係る制御装置について説明する。第２の実施形態に係る制御装置は、第１の実施形態に係る制御装置１に対し、重畳成分生成部によるパイロット電圧Ｖｈの生成方法と、磁極位置推定部による磁極位置θの推定方法とが異なる。
【０１０１】
図１０は、第２の実施形態に係る制御装置の構成を示す図である。なお、図１０においては、重複説明を避けるため、第１の実施形態に係る制御装置１と異なる制御部３０Ａの構成のみを示しており、その他の構成は省略している。また、第１の実施形態と同様の構成については同一符号を付し、適宜重複説明を省略するものとする。
【０１０２】
図１０に示すように、第２の実施形態に係る制御装置１Ａは、３相／２相座標変換部３１と、回転座標変換部３２と、電流指令出力部３３と、電流制御部３４と、カウンタ部３５と、重畳成分生成部３６Ａと、加算部３７、３８と、駆動信号生成部３９と、磁極位置推定部４０Ａとを備える。制御装置１Ａは、重畳成分生成部３６Ａおよび磁極位置推定部４０Ａを除き、第１の実施形態に係る制御装置１と同様の構成である。
【０１０３】
上述したように、突極性を有する電動機３において、電動機３の固定子３ａ上に設定されるαβ軸座標系における電流・電圧方程式は、上記式（３）のように表すことができる。さらに、重畳する２つの電圧ベクトル間の時間間隔が、駆動電圧周期に比べて十分に短い場合、駆動電圧周波数の成分を取り除いた上記式（３）の電流・電圧方程式は、下記式（２２）のように近似できる。この式（２２）は、上記式（６）と同一である。
【０１０４】
【数１６】
0005900600.tif 000017
【０１０５】
上記式（２２）を時間微分すると、下記式（２３）を得ることができる。
【０１０６】
【数１７】
0005900600.tif 000018
【０１０７】
電動機３における回転加速度が、上記式（２３）の右辺第１項に比べて十分小さいと仮定すると、上記式（２３）は、下記式（２４）のように近似できる。
【０１０８】
【数１８】
0005900600.tif 000019
【０１０９】
上記式（２４）は、電流に関して下記式（２５）に示すように書き直すことができる。
【０１１０】
【数１９】
0005900600.tif 000020
【０１１１】
ここで、インバータ部１０の出力電圧として、駆動電圧指令であるｑ軸電圧指令Ｖｑ*およびｄ軸電圧指令Ｖｄ*に基づく駆動電圧に、パイロット電圧Ｖｈが含まれるようにすることを考える。例えば、図１１に示すように、αβ軸座標系において、電圧ベクトルＶｈ０、Ｖｈ１、Ｖｈ２をこの順に一定周期で重畳する。図１１は、αβ軸座標系における電圧ベクトルＶｈ０、Ｖｈ１、Ｖｈ２とφとの関係の一例を示す図である。
【０１１２】
電圧ベクトルＶｈ０は、α軸の正の向きから位相差φ［ｄｅｇ］のベクトル、電圧ベクトルＶｈ１はα軸の正の向きから位相差φ＋９０［ｄｅｇ］のベクトル、電圧ベクトルＶｈ２は、α軸の正の向きから位相差φ＋１８０［ｄｅｇ］のベクトルである。また、電圧ベクトルＶｈ０、Ｖｈ１、Ｖｈ２は、振幅Ｖｉｎｊ２（＞０）のベクトルである。
【０１１３】
電圧ベクトルＶｈ０と電圧ベクトルＶｈ１の差分ベクトルΔＶｈ０（＝Ｖｈ０―Ｖｈ１）は、下記式（２６）に示すように表され、電圧ベクトルＶｈ１と電圧ベクトルＶｈ２の差分ベクトルΔＶｈ１（＝Ｖｈ１―Ｖｈ２）は、下記式（２７）に示すように表される。
【０１１４】
【数２０】
0005900600.tif 000021
【０１１５】
電圧ベクトルＶｈ０の重畳時刻を時刻ｔｈ０とし、電圧ベクトルＶｈ１の重畳時刻を時刻ｔｈ１＝ｔｈ０＋Δｔとし、電圧ベクトルＶｈ０の重畳時刻を時刻ｔｈ２＝ｔｈ０＋２Δｔとする。また、電圧微分ベクトルを、差分ベクトルで近似する。この場合、上記式（２６）、（２７）に示す差分ベクトルから、時刻ｔｈ１と時刻ｔｈ２での電流・電圧方程式はそれぞれ下記式（２８）、（２９）に示すように表すことができる。
【０１１６】
【数２１】
0005900600.tif 000022
【０１１７】
さらに、上記式（２８）、（２９）から、下記式（３０）の関係式を導くことができる。
【０１１８】
【数２２】
0005900600.tif 000023
【０１１９】
上記式（３０）は、電動機３の磁極位置θに関して下記式（３１）に示すように書き直すことができる。したがって、時刻ｔｈ１と時刻ｔｈ２での各相の電流変化を検出すれば電動機３の磁極位置θを推定することができる。
【０１２０】
【数２３】
0005900600.tif 000024
【０１２１】
ここで、φを４５、１３５、２２５、３１５［ｄｅｇ］のいずれかに設定すると、後述するように、電動機３の磁極位置θをより簡単な演算式で推定することができる。すなわち、上記式（３１）は、φ＝４５［ｄｅｇ］の場合、下記式（３２）のように、φ＝１３５［ｄｅｇ］の場合、は下記式（３３）のように、φ＝２２５［ｄｅｇ］の場合、下記式（３４）のように、φ＝３１５［ｄｅｇ］の場合、下記式（３５）のように、それぞれ表すことができる。
【０１２２】
【数２４】
0005900600.tif 000025
【０１２３】
したがって、電圧ベクトルＶｈ０を、φ＝４５、１３５、２２５、３１５［ｄｅｇ］のいずれかの状態になるように重畳し、一定期間後、電圧ベクトルＶｈ０に対し９０［ｄｅｇ］の位相差とした電圧ベクトルＶｈ１を重畳する。さらに一定期間後、電圧ベクトルＶｈ１に対し９０［ｄｅｇ］の位相差とした電圧ベクトルＶｈ２を重畳する。これにより、電動機３の磁極位置θをより簡単な演算式で推定できる。
【０１２４】
また、時刻ｔｈ１でα軸の正の向きから位相差φ＋９０［ｄｅｇ］の電圧ベクトルＶｈ１を重畳し、時刻ｔｈ２でα軸の正の向きから位相差φ＋１８０［ｄｅｇ］の電圧ベクトルＶｈ１を重畳した場合を説明したが、これに限定されるものではない。例えば、時刻ｔｈ１でα軸の正の向きから位相差φ−９０［ｄｅｇ］の電圧ベクトルＶｈ１を重畳し、時刻ｔｈ２でα軸の正の向きから位相差φ−１８０［ｄｅｇ］の電圧ベクトルＶｈ１を重畳しても同様に、電動機３の磁極位置θを簡単な演算式で推定できる。
【０１２５】
このように、時刻ｔｈ０においてφを０、９０、１８０、２７０［ｄｅｇ］のいずれかに設定して電圧ベクトルＶｈ０を重畳した後、時刻ｔｈ１でφ＋９０［ｄｅｇ］（またはφ−９０［ｄｅｇ］）として電圧ベクトルＶｈ１を重畳する。さらに、時刻ｔｈ２でφ＋１８０［ｄｅｇ］（またはφ−１８０［ｄｅｇ］）として電圧ベクトルＶｈ２を重畳する。そして、時刻ｔｈ０と時刻ｔｈ１での各相の電流変化と、時刻ｔｈ１と時刻ｔｈ２での各相の電流変化とから２階差分値を検出し、上記式（３２）〜（３５）などの簡単な演算式に基づいた演算によって、電動機３の磁極位置θを推定できる。
【０１２６】
そこで、重畳成分生成部３６Ａは、下記式（３６）に示すように、時間周期Ｔ＝４ＴＶのパイロット電圧Ｖｈを生成し、α軸指令成分Ｖα*およびβ軸指令成分Ｖβ*で規定される駆動電圧指令に重畳するようにしている。ＴＶは、空間ベクトル変調を行う駆動信号生成部３９の更新周期である。
【０１２７】
【数２５】
0005900600.tif 000026
【０１２８】
上記式（３６）に示すように重畳成分生成部３６Ａから出力するパイロット電圧Ｖｈを設定することで、図１２に示すように、α軸およびβ軸で打ち消し合う方向にパイロット電圧Ｖｈを出力できる。図１２は、パイロット電圧Ｖｈの時間周期Ｔにおける変化を示す図である。そのため、時間周期Ｔにおける平均電圧を０とすることができ、これにより、トルクリップルを低減することができる。
【０１２９】
重畳成分生成部３６Ａの構成および動作について図１３および図１４を参照してさらに具体的に説明する。図１３は、重畳成分生成部３６Ａから出力されるパイロット電圧Ｖｈのαβ軸成分であるＶαｈ、Ｖβｈの変化を示す図であり、図１４は、重畳成分生成部３６Ａのパイロット電圧生成処理を示すフローチャートである。
【０１３０】
図１３に示すように、重畳成分生成部３６Ａは、カウント値ＣＮＴとして０、１、２、３の順に繰り返し出力するカウンタ部３５のカウント値ＣＮＴに基づき、パイロット電圧Ｖｈのα軸成分およびβ軸成分であるα軸指令成分Ｖαｈおよびβ軸指令成分Ｖβｈを生成する。
【０１３１】
具体的には、図１４に示すように、重畳成分生成部３６Ａは、カウンタ部３５のカウント値ＣＮＴを判定し（ステップＳ３０）、カウント値ＣＮＴ＝０の場合、Ｖαｈ＝＋Ｖｉｎｊ２およびＶβｈ＝＋Ｖｉｎｊ２としたパイロット電圧Ｖｈを生成して、加算部３７、３８へ出力する（ステップＳ３１）。カウント値ＣＮＴ＝１の場合、重畳成分生成部３６Ａは、Ｖαｈ＝−Ｖｉｎｊ２およびＶβｈ＝＋Ｖｉｎｊ２としたパイロット電圧Ｖｈを生成して、加算部３７、３８へ出力する（ステップＳ３２）。
【０１３２】
また、カウント値ＣＮＴ＝２の場合、重畳成分生成部３６Ａは、Ｖαｈ＝−Ｖｉｎｊ２およびＶβｈ＝−Ｖｉｎｊ２としたパイロット電圧Ｖｈを生成して、加算部３７、３８へ出力する（ステップＳ３３）。カウント値ＣＮＴ＝３の場合、重畳成分生成部３６Ａは、Ｖαｈ＝＋Ｖｉｎｊ２およびＶβｈ＝−Ｖｉｎｊ２としたパイロット電圧Ｖｈを生成して、加算部３７、３８へ出力する（ステップＳ３４）。
【０１３３】
このように、重畳成分生成部３６は、αβ軸座標系で、時間周期Ｔｖ毎に前回生成したパイロット電圧Ｖｈに対して９０度の方向となるパイロット電圧Ｖｈを繰り返し生成するようにしている。なお、重畳成分生成部３６Ａは、αβ軸座標系で、時間周期Ｔｖ毎に前回生成したパイロット電圧Ｖｈに対して−９０度の方向となるパイロット電圧Ｖｈを繰り返し生成することもできる。
【０１３４】
磁極位置推定部４０Ａは、時刻ｔ＝ｎＴ（ｎ＝１、２、３、・・・）での磁極位置θ、時刻ｔ＝Ｔ／４＋ｎＴでの磁極位置θ、時刻ｔ＝Ｔ／２＋ｎＴでの磁極位置θ、時刻ｔ＝３Ｔ／４＋ｎＴでの磁極位置θを推定する。
【０１３５】
図１５は、磁極位置推定部４０Ａの構成を示す図である。図１５に示すように、磁極位置推定部４０Ａは、差分演算部６０〜６３と、１サンプル遅延部６４、６５と、θ演算部６６とを備え、Ｔ／４をサンプル周期として動作する。
【０１３６】
差分演算部６０は、３相／２相座標変換部３１からα軸成分の電流検出値Ｉαを取得し、かかる電流検出値Ｉαと１サンプル時間（＝Ｔ／４）前に取得した電流検出値Ｉαとの差分値ΔＩαを演算して出力する。また、差分演算部６１は、３相／２相座標変換部３１からβ軸成分の電流検出値Ｉβを取得し、かかる電流検出値Ｉβと時間Ｔ／４前に取得した電流検出値Ｉβとの差分値ΔＩβを演算して出力する。
【０１３７】
差分演算部６２は、差分演算部６０から差分値ΔＩαを取得し、かかる差分値ΔＩαと１サンプル時間前に取得した差分値ΔＩαとの２階差分値Δ２Ｉαを演算して出力する。また、差分演算部６３は、差分演算部６１から差分値ΔＩβを取得し、かかる差分値ΔＩβと１サンプル時間前に取得した差分値ΔＩβとの２階差分値Δ２Ｉβを演算して出力する。
【０１３８】
１サンプル遅延部６４は、差分演算部６２から出力される２階差分値Δ２Ｉαを１サンプル時間分遅延させて、θ演算部６６へ出力する。また、１サンプル遅延部６５は差分演算部６３から出力される２階差分値Δ２Ｉβを１サンプル時間分遅延させて、θ演算部６６へ出力する。
【０１３９】
θ演算部６６は、差分演算部６２、６３から出力される２階差分値Δ２Ｉα、Δ２Ｉβと、１サンプル遅延部６４、６５から出力される１サンプル時間前の２階差分値Δ２Ｉα、Δ２Ｉβとに基づいて、磁極位置θを推定する。
【０１４０】
図１６は、磁極位置推定部４０Ａの磁極位置推定処理を示すフローチャートである。なお、以下において、Δ２Ｉαｈ２＝ΔＩαｈ２−ΔＩαｈ１、Δ２Ｉβｈ２＝ΔＩβｈ２−ΔＩβｈ１、Δ２Ｉαｈ１＝ΔＩαｈ１−ΔＩαｈ０、Δ２Ｉβｈ１＝ΔＩβｈ１−ΔＩβｈ０である。また、ΔＩαｈ２＝Ｉαｈ２−Ｉαｈ１、ΔＩβｈ２＝Ｉβｈ２−Ｉβｈ１、ΔＩαｈ１＝Ｉαｈ１−Ｉαｈ０、ΔＩβｈ１＝Ｉβｈ１−Ｉβｈ０、ΔＩαｈ０＝Ｉαｈ０−Ｉαｈ−１、ΔＩβｈ０＝Ｉβｈ０−Ｉβｈ−１である。
【０１４１】
図１６に示すように、磁極位置推定部４０Ａは、カウンタ部３５のカウント値ＣＮＴを判定する（ステップＳ４０）。
【０１４２】
カウント値ＣＮＴ＝０の場合、磁極位置推定部４０Ａは、下記式（３７）に基づいて、電動機３の磁極位置θを推定する（ステップＳ４１）。下記式（３７）は、上記式（３２）において電流２階微分を電流差分で近似したものである。
【０１４３】
【数２６】
0005900600.tif 000027
【０１４４】
なお、カウント値ＣＮＴ＝０の場合、電流検出値Ｉαｈ２、Ｉαｈ１、Ｉαｈ０、Ｉαｈ−１は、それぞれ時刻ｔ＝ｎＴ、−Ｔ／４＋ｎＴ、−Ｔ／２＋ｎＴ、−３Ｔ／４＋ｎＴにおけるα軸成分の電流検出値である。また、電流検出値Ｉβｈ２、Ｉβｈ１、Ｉβｈ０、Ｉβｈ−１は、それぞれ時刻ｔ＝ｎＴ、−Ｔ／４＋ｎＴ、−Ｔ／２＋ｎＴ、−３Ｔ／４＋ｎＴにおけるβ軸成分の電流検出値である。
【０１４５】
また、カウント値ＣＮＴ＝１の場合、磁極位置推定部４０Ａは、下記式（３８）に基づいて、電動機３の磁極位置θを推定する（ステップＳ４２）。なお、下記式（３８）は、上記式（３３）において電流２階微分を電流差分で近似したものである。
【０１４６】
【数２７】
0005900600.tif 000028
【０１４７】
なお、カウント値ＣＮＴ＝１の場合、電流検出値Ｉαｈ２、Ｉαｈ１、Ｉαｈ０、Ｉαｈ−１は、それぞれ時刻ｔ＝Ｔ／４＋ｎＴ、ｎＴ、−Ｔ／４＋ｎＴ、−Ｔ／２＋ｎＴにおけるα軸成分の電流検出値である。また、電流検出値Ｉβｈ２、Ｉβｈ１、Ｉβｈ０、Ｉβｈ−１は、それぞれ時刻ｔ＝Ｔ／４＋ｎＴ、ｎＴ、−Ｔ／４＋ｎＴ、−Ｔ／２＋ｎＴにおけるβ軸成分の電流検出値である。
【０１４８】
また、カウント値ＣＮＴ＝２の場合、磁極位置推定部４０Ａは、下記式（３９）に基づいて、電動機３の磁極位置θを推定する（ステップＳ４３）。なお、下記式（３９）は、上記式（３４）において電流２階微分を電流差分で近似したものである。
【０１４９】
【数２８】
0005900600.tif 000029
【０１５０】
なお、カウント値ＣＮＴ＝２の場合、電流検出値Ｉαｈ２、Ｉαｈ１、Ｉαｈ０、Ｉαｈ−１は、それぞれ時刻ｔ＝Ｔ／２＋ｎＴ、Ｔ／４＋ｎＴ、ｎＴ、−Ｔ／４＋ｎＴにおけるα軸成分の電流検出値である。また、電流検出値Ｉβｈ２、Ｉβｈ１、Ｉβｈ０、Ｉβｈ−１は、それぞれ時刻ｔ＝Ｔ／２＋ｎＴ、Ｔ／４＋ｎＴ、ｎＴ、−Ｔ／４＋ｎＴにおけるβ軸成分の電流検出値である。
【０１５１】
また、カウント値ＣＮＴ＝３の場合、磁極位置推定部４０Ａは、下記式（４０）に基づいて、電動機３の磁極位置θを推定する（ステップＳ４４）。なお、下記式（４０）は、上記式（３５）において電流２階微分を電流差分で近似したものである。
【０１５２】
【数２９】
0005900600.tif 000030
【０１５３】
なお、カウント値ＣＮＴ＝３の場合、電流検出値Ｉαｈ２、Ｉαｈ１、Ｉαｈ０、Ｉαｈ−１は、それぞれ時刻ｔ＝３Ｔ／４＋ｎＴ、Ｔ／２＋ｎＴ、Ｔ／４＋ｎＴ、ｎＴにおけるα軸成分の電流検出値である。また、電流検出値Ｉβｈ２、Ｉβｈ１、Ｉβｈ０、Ｉβｈ−１は、それぞれ時刻ｔ＝３Ｔ／４＋ｎＴ、Ｔ／２＋ｎＴ、Ｔ／４＋ｎＴ、ｎＴにおけるβ軸成分の電流検出値である。
【０１５４】
ステップＳ４１〜Ｓ４４において、磁極位置推定部４０Ａは、差分演算部６２、６３から出力される差分値Δ２Ｉα、Δ２ＩβをΔ２ｉαｈ１、Δ２ｉβｈ１とし、１サンプル遅延部６４、６５から出力される１サンプル時間前の差分値Δ２Ｉα、Δ２ＩβをΔ２ｉαｈ２、Δ２ｉβｈ２として、上記式（３７）〜（４０）の演算を行うことで、電動機３の磁極位置θを推定する。
【０１５５】
以上のように、第２の実施形態に係る制御装置１Ａでは、重畳成分生成部３６Ａと、磁極位置推定部４０Ａとを備える。重畳成分生成部３６Ａは、αβ軸座標系のα軸またはβ軸に対して４５度の位相差をもつベクトルとしてパイロット電圧Ｖｈを生成し、かつ、ベクトルの向きが前回生成したパイロット電圧Ｖｈに対して９０度異なるパイロット電圧Ｖｈを所定周期で生成する。
【０１５６】
したがって、磁極位置推定部４０Ａにおける磁極位置θの推定を、簡易な演算式に基づいて行うことができ、電動機３を駆動するための駆動電圧指令の更新周期に影響を与えることもない。しかも、パイロット電圧Ｖｈの重畳を開始してから時間周期Ｔで電動機３の磁極位置θを推定でき、さらに、時間周期Ｔの１／４の時間で繰り返し、電動機３の磁極位置θを推定することができるため、磁極位置推定の応答性を向上させることができる。
【０１５７】
上述においては、駆動信号生成部３９において空間ベクトル変調法を用いてインバータ部１０を駆動する駆動信号Ｓ１〜Ｓ６を生成する例を説明したが、第１の実施形態の欄において説明した例と同様の方法で、搬送波比較変調法を用いることもできる。この場合、カウンタ部３５は、時間周期Ｔを搬送波の時間周期Ｔｃの２倍の長さとし、時間周期Ｔを４周期分として、カウント値ＣＮＴを変更する。カウンタ部３５がカウント値ＣＮＴを変更するタイミングは、搬送波の山および谷である。
【０１５８】
また、第１の実施形態の欄において説明したその他の種々の変形についても、同様に、上述した第２の実施形態に係る制御装置１Ａに適用することができる。
【０１５９】
さらなる効果や変形例は、当業者によって容易に導き出すことができる。このため、本発明のより広範な態様は、以上のように表しかつ記述した特定の詳細および代表的な実施形態に限定されるものではない。したがって、添付の特許請求の範囲およびその均等物によって定義される総括的な発明の概念の精神または範囲から逸脱することなく、様々な変更が可能である。
【符号の説明】
【０１６０】
１、１Ａ  電動機の制御装置
２  直流電源
３  電動機
３ａ  電動機の固定子
１０  インバータ部
２０  電流検出部
３０、３０Ａ  制御部
３１  ３相／２相座標変換部３１（電流検出部）
３６、３６Ａ  重畳成分生成部
３９  駆動信号生成部
４０、４０Ａ  磁極位置推定部
５０、５１、６０〜６３  差分演算部（電流差分演算部）
５４、６６  θ演算部（磁極位置演算部）
【図１】
0005900600.tif 000031
【図２Ａ】
0005900600.tif 000032
【図２Ｂ】
0005900600.tif 000033
【図３】
0005900600.tif 000034
【図４】
0005900600.tif 000035
【図５】
0005900600.tif 000036
【図６】
0005900600.tif 000037
【図７】
0005900600.tif 000038
【図８】
0005900600.tif 000039
【図９】
0005900600.tif 000040
【図１０】
0005900600.tif 000041
【図１１】
0005900600.tif 000042
【図１２】
0005900600.tif 000043
【図１３】
0005900600.tif 000044
【図１４】
0005900600.tif 000045
【図１５】
0005900600.tif 000046
【図１６】
0005900600.tif 000047