(57)【特許請求の範囲】
【請求項１】
負荷装置を駆動するモータを、任意波形のトルク指令値に出力制御するトルク制御回路を設けたモータのトルク制御装置において、
モータのトルク出力からトルク周期性外乱をオブザーバで直接に推定し、この推定値でトルク指令値を補正してトルク周期性外乱を抑制する周期性外乱抑制制御手段と、
前記トルク出力の周波数成分による電気・機械系の共振現象を、前記トルク指令値からトルク検出値までの周波数成分別のシステム伝達関数として同定し、この周波数成分別の同定結果の逆数を用いた共振抑制テーブルまたは共振抑制演算でトルク指令値を生成する機械共振抑制制御手段を備えたことを特徴とするモータのトルク制御装置。
【請求項２】
前記機械共振抑制制御手段は、前記システム同定結果を、システム特性を近似的な数式で表現することを特徴とする請求項１に記載のモータのトルク制御装置。
【請求項３】
前記トルク制御回路は、前記トルク指令値をベクトル制御における回転座標系のｄ軸ｑ軸電流指令値に変換してモータ電流を制御する構成とし、
前記周期性外乱抑制制御手段は、前記ｑ軸電流指令値から卜ルク検出値までのシステム同定結果に基づいて、その逆数を用いて周波数成分ごとに前記ｑ軸電流指令値のみのトルク指令値を生成することを特徴とする請求項１または２に記載のモータのトルク制御装置。
【請求項４】
前記機械共振抑制制御手段は、直流成分を含む複数の周波数成分を有する任意波形のトルク指令値を直流分のトルク成分と直流以外のトルク成分に分離し、直流分のトルク指令値はｄ軸ｑ軸電流指令値の直流分として生成し、直流以外のトルク指令値は前記共振抑制テーブルまたは共振抑制演算で共振抑制を考慮したｑ軸電流指令値に変換することを特徴とする請求項１〜３のいずれか１項に記載のモータのトルク制御装置。
【請求項５】
前記機械共振抑制制御手段は、直流成分を含む複数の周波数成分を有する任意波形のトルク指令値を直流分のトルク成分と直流以外のトルク成分に分離し、直流分のトルク成分はｄ軸ｑ軸電流指令値に変換し、直流以外のトルク成分は前記共振抑制テーブルまたは共振抑制演算により機械共振特性を考慮したトルク指令値に変換することを特徴とする請求項１〜４のいずれか１項に記載のモータのトルク制御装置。
【請求項６】
前記機械共振抑制制御手段は、直流成分を含む複数の周波数成分を有する任意波形のトルク指令値を直流分のトルク成分と直流以外のトルク成分に分離し、直流分のトルク成分はｄ軸ｑ軸電流指令値に変換し、直流以外のトルク成分は前記オブザーバが周期性外乱を複素ベクトルで表現して周期性外乱を求める場合の実部成分ｄＩAn*（余弦係数成分）と虚部成分ｄＩBn*（正弦係数成分）に分解して該オブザーバに周期性外乱電流指令値として入力することを特徴とする請求項１〜５のいずれか１項に記載のモータのトルク制御装置。
【請求項７】
前記機械共振抑制制御手段は、直流成分を含む複数の周波数成分を有する任意波形のトルク指令値を直流分のトルク成分と直流以外のトルク成分に分離し、直流分のトルク成分はｄ軸ｑ軸電流指令値に変換し、直流以外のトルク成分から前記トルクリプル周期性外乱と同期する周波数成分のみを抽出して前記オブザーバが周期性外乱を複素ベクトルで表現して周期性外乱を求める場合の実部成分ｄＩAn*（余弦係数成分）と虚部成分ｄＩBn*（正弦係数成分）に分解して該オブザーバに周期性外乱電流指令値として入力することを特徴とする請求項１〜５のいずれか１項に記載のモータのトルク制御装置。
【請求項８】
前記機械共振抑制制御手段および前記周期性外乱抑制制御手段は、直流成分を含む複数の周波数成分を有する任意波形のトルク指令値を重ね合わせて任意波形のトルク出力を発生させる場合、並列かつ独立に任意周波数成分別のトルク指令値、複素ベクトル表現の電流指令値を生成することを特徴とする請求項１〜７のいずれか１項に記載のモータのトルク制御装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、モータのトルク制御装置に係り、特にモータに任意波形のトルク出力を発生させるモータの制御装置におけるトルク出力の周期性外乱抑制および電気・機械共振抑制に関する。
【背景技術】
【０００２】
モータは原理的にトルクリプル（トルク脈動）を発生し、振動・騒音・乗り心地への悪影響、電気・機械共振等の種々の問題を引き起こす。特に、近年普及が進んでいる埋込磁石ＰＭモータ（ＩＰＭＳＭ）は、コギングトルクリプルとリラクタンストルクリプルが複合的に発生する。その対策として、トルクリプルを打ち消す補償電流をトルク指令に重畳する種々の方式が検討されている。
【０００３】
しかしながら、例えば数式解析モデルを用いてフィードフォワード補償を行なう方式では、解析誤差の影響が懸念される。また、定常動作点でのフィードバック学習制御結果を記憶してフィードフォワード補償する方式では動作点毎に制御パラメータを適切に調整するための時間を要し、オンライン補償が難しくなる。また、電流リプルを低減する方式は、トルクリプルの観点からは最適に抑制されているとは限らない。そのほか、トルクリプルオブザーバ補償方式も検討されているが、可変速運転時の特性やオンラインフィードバック抑制の検証が不充分なものであった。
【０００４】
上記の課題に対し、電気・機械共振の元凶であるトルクリプルを高精度に抑制する目的で、軸トルクメータによるフィードバック抑制制御法を本願発明者等は既に提案している（例えば、非特許文献１参照）。この制御法は、トルクリプルの周期性に着目して脈動周波数成分毎に補償する制御系を構築すると共に、システム同定結果を用いて動作状態変化に即応するようパラメータを自動調整する機能を設けたものであり、この制御法の詳細を以下に説明する。
【０００５】
（１）トルクリプル抑制制御装置の基本構成
図８は、従来のトルクリプル抑制制御装置の制御ブロック図である。同図は、モータによって負荷をトルク加振するシステムに適用したものである。
【０００６】
トルクリプルの発生源となるモータ１と、何らかの負荷装置２をシャフト３で結合し、その軸トルクをトルクメータ４で計測してトルクリプル抑制装置５に入力する。また、ロータリエンコーダ等の回転位置センサ６を用いてモータの回転子位置（位相）情報を入力する。トルクリプル抑制装置５は、トルク脈動抑圧手段を搭載し、トルク指令値（あるいは速度指令値）に基づいて生成された電流指令値に、トルク脈動補償電流を上乗せした指令値をインバータ７に与える。図８の例では、インバータ７で電流ベクトル制御することを考慮して、モータの回転に同期した回転座標（直交ｄｑ軸）上のｄ軸、ｑ軸電流指令値ｉｄ*、ｉｑ*を与えている。
【０００７】
トルクリプル（トルク脈動）はモータの構造上、回転子位置に応じて周期的に発生することが知られている。そこで、モータ回転に同期してトルクリプル周波数成分を抽出する手段を用い、任意次数ｎのトルクリプルを余弦、正弦係数ＴAn，ＴBn［Ｎｍ］に変換する。トルクリプル周波数成分の厳密な計測手段にはフーリエ変換などあがるが、演算容易性を重視すれば回転位相θ［ｒａｄ］を基準とした単相の高調波回転座標系に低域通過フィルタを通すことでトルクリプル周波数成分を抽出することができる。
【０００８】
トルクリプル抑制装置５では、上記の余弦、正弦係数ＴAn，ＴBnを用いてトルクリプル抑制制御を行い、任意周波数成分の補償電流ｉｑｃ*［アンペア］の余弦／正弦係数ＩAn，ＩBn［アンペア］を生成する。補償電流ｉｑｃ*への変換は、変換時と同じ回転位相θを用いて以下の（２）式で求め、この補償電流はｑ軸電流指令値に重畳して通常のベクトル制御を行う。
【０００９】
【数１】
0005488043.tif 000002
【００１０】
図９は、前記の軸トルクＴdetと回転位相θからトルクリプルの余弦、正弦係数ＴAn，ＴBnを求め、これらと回転速度ωから補償電流余弦係数ＩAnおよび正弦係数ＩBnを求めてトルク周期性外来抑制の補償電流を求めるトルクリプル抑制装置の制御ブロック図である。図中の記号は、以下の意味である。
【００１１】
Ｔ*：トルク指令値、Ｔdet：トルク検出値、ＴAn：ｎ次トルク脈動抽出成分（余弦係数）、ＴBn：ｎ次トルク脈動抽出成分（正弦係数）、ω：回転数検出値、θ：回転位相検出値、ｉｑｃ*：トルクリプル補償電流、ｉｄ：ｄ軸電流検出値、ｉｄ*：ｄ軸電流指令値、ｉｑ：ｑ軸電流検出値、ｉｑ*：ｑ軸電流指令値、ｉｕ．ｉｖ．ｉｗ：ｕ、ｖ、ｗ相電流、ｉｑｏ*：ｑ軸電流指令値（補償電流重畳前）、ＩAn：ｎ次補償電流余弦係数、ｌBn：ｎ次補償電流正弦係数、ａｂｚ：回転センサ信号。なお、添え字のｎはｎ次成分トルクリプルである。
【００１２】
図９において、指令値変換部１１は、トルク指令値Ｔ*から、ベクトル制御における回転ｄｑ座標系のｄ軸およびｑ軸電流指令値Ｉｄ*、Ｉｑｏ*に変換する。一般には、最大トルク／電流制御を実現するような変換数式やテーブルなどが用いられる。電流ベクトル制御部１２は、ｑ軸電流指令値ｉｑｏ*にトルク脈動補償電流ｉｑｃ*を重畳したものをｑ軸電流指令値ｉｑ*としてトルク脈動を抑制する。図９の例ではｑ軸電流指令値に補償電流ｉqcを重畳しているが、ｄ軸電流、あるいはｄ軸とｑ軸の両方に与えても良い。あるいは、ｄｑ軸電流の干渉が問題にならないシステムであれば、トルク指令値に対して直接的にトルク脈動補償信号を重畳しても良い。
【００１３】
電流ベクトル制御部１２は、一般的な直交回転座標系ｄ軸ｑ軸において、電流ベクトル制御の動作を行い、モータ（ＩＰＭＳＭ）１３をベクトル制御で駆動することで、負荷装置１４を駆動する。座標変換部１５は、電流センサ１６で検出する３相交流電流ｉｕ、ｉｖ、ｉｗとモータ回転位相θから、モータ回転座標に同期したｄｑ軸直交回転座標系の電流ｉｄ，ｉｑに変換する。回転位相／速度検出部１７は、エンコーダ等の回転位置センサ１８の回転センサ信号ａｂｚから回転速度ωおよび回転位相θの情報に変換する。
【００１４】
トルク脈動周波数成分抽出部１９では、軸トルクメータ２０で検出する軸トルク検出値Ｔdetと回転位相θからトルク周期性外乱を脈動周波数成分ごとに抽出する。代表的な手段としてフーリエ変換がある。脈動成分の抽出方法は任意であるが、演算簡素化の例として（３）〜（５）式のような手段を用いることにする。
【００１５】
【数２】
0005488043.tif 000003
【００１６】
£：ラプラス変換、ωf：脈動抽出ローパスフィルタ遮断周波数［ｒａｄ／ｓ］、ｓ：ラプラス演算子
トルクリプル抑制制御部２１では、トルクリプル周波数成分に同期した座標系で、その余弦係数成分ＩAnと正弦係数成分ＩBnのそれぞれを求める。このトルクリプル抑制制御部２１の代表形態として周期性外乱オブザーバ補償法または補償電流フーリエ係数学習制御法を用いることができる。
【００１７】
補償電流生成部２２では、前記の（２）式で補償電流指令値ｉｑｃ*を生成し、ｑ軸電流指令値に重畳する。なお、同図ではｑ軸電流指令値に補償電流指令値を重畳するが、ｄ軸電流指令値、ｄ軸とｑ軸の双方の電流指令値、トルク指令値などに置き換えることも可能である。
【００１８】
（２）システムの同定
図８に示すようなシステム構成は、ＰＭモータ１，負荷装置２，卜ルクメータ４，カップリング３等の慣性モーメントにより多慣性軸ねじれ共振系となる。軸トルク検出値をフィードバックする場合は複数の共振・反共振周波数があるため、動作状態に応じて適切に抑制制御パラメータを決定しなければならない。制御パラメータの学習時間か長いと電気・機械共振現象を増大させる危険があるため、速やかな自動調整機能が必要である。
【００１９】
そこで、非特許文献１では、回転速度変化に適応する可変ノミナル制御パラメータを導出するために、図８のトルクリプル抑制装置５の出力値から入力値までのシステム伝達特性、すなわち図９における補償電流指令値ｉｑｃ*から軸トルク検出器４の検出値Ｔdetまでの周波数伝達特性を同定する。システム同定手法は任意であるが、開ループでｉｑｃ*にガウス性ノイズ信号を与えた時の軸トルク検出値Ｔdetを演算周期１００μｓで２０秒間計測し、入出力のパワースペクトル密度の比から周波数伝達関数をノンパラメトリックに推定した結果を図１０に示す（機械系，インバータ電流応答，トルクメータ応答，無駄時間などを含んだ実機の特性）。また、周波数伝達特性の傾向から４慣性系に近似した場合のパラメトリック同定結果も併せて示す。近似のための最適化手法にも種々の方式があるが、周波数領域でｌｋＨｚまでの振幅特性の誤差を評価し、制約付き非線形最小化（逐次２次計画法）を行った。
【００２０】
図９ではトルクリプル周波数に同期した座標で制御系を構築するため、図１０のシステム同定結果から任意の周波数伝達特性のみを抽出する。定常状態において、トルクリプル周波数に同期したシステムの振幅・位相伝達特性は１次元複素ベクトルで表現できるので、図９の制御系でのシステム特性Ｐsysを下記の（６）式のように定義する。
【００２１】
【数３】
0005488043.tif 000004
【００２２】
ＰAm：システム特性の実部、ＰBm：システム特性の虚部、ｍ：システム同定テーブルの周波数要素番号
例えば、１〜１０００［Ｈｚ］までのシステム特性をＩＨｚ毎に（３）式で表現した場合、１０００個の複素ベクトルの要素からシステム同定テーブルを構築することができる。制御系で使用されるのは常に１つの複素ベクトルに限られ、回転速度変化（トルクリプル周波数変化）に応じて同定テーブルから瞬時にＰAmとＰBmを読み出し、線形補間を施して複素ベクトル化された同定結果を抑制制御に適用する。なお、回転位相を基準とした実部と虚部の軸を定義するため、（５）式における余弦係数は実部成分、正弦係数は虚部成分に対応する。
【００２３】
（３）補償電流フーリエ係数学習制御法
前記の非特許文献１に方式１として記載されるトルクリプル抑制制御法であり、トルクリプル周波数成分のフーリエ係数として求め、これから前記（２）式の演算で補償電流ｉｑｃ*を求める。この制御法では、トルクリプル周波数に同期した周波数成分のシステム伝達関数を１次元複素ベクトルで表現し、任意周波数成分のトルクリプルの実部・虚部をフーリエ変換等で抽出している。その余弦・正弦フーリエ係数を複素ベクトルの実部・虚部に当てはめて、フィードバック抑制制御系を構築する。
【００２４】
補償電流係数はＩ−Ｐ（比例・積分）学習制御方式で求める。Ｉ−Ｐ抑制制御系の閉ループ特性が、モデルマッチング法によって任意の標準系規範モデルの極配置と一致するように比例・積分ゲインを決定する。また、それらは前記のシステム同定結果と回転速度情報を用いて自動的にパラメータを適応させるので、多慣性共振系システムヘの実装を容易にする。
【００２５】
任意の定常動作点（定常トルク・定常回転数）において、抑制完了したときの補償電流信号の振幅・位相を保存し、それを複数の動作点で実施して、トルク・回転数の２次元テーブルとして実装する。この際、トルク・回転数情報をテーブルに入力し、読み出した補償電流振幅・位相データから補償電流を生成してフィードフォワード抑制することが可能となる。
【００２６】
（４）周期性外乱オブザーバ補償法
前記の非特許文献１に方式２として記載されるトルクリプル抑制制御法である。前記補償電流フーリエ係数学習制御法での制御パラメータ自動調整手法は、Ｉ−Ｐ制御ゲインの調整を介して外乱を抑制するため、可変速運転に対する即応性が低下する。そのため、学習した結果を予めテーブル化しておき、フィードフォワード抑制する使い方を推奨している。
【００２７】
本方式の周期性外乱オブザーバ補償法は、周期性外乱オブザーバという考え方を用いて直接的にトルクリプル外乱を推定するため、上記の即応性の問題を改善できる。したがって、可変速・負荷変動を有するシステムに対しても常時オンラインフィードバックでトルクリプルを抑制できる。また、補償電流フーリエ係数学習制御法と同様に１次元複素ベクトルで表現したシステム同定結果を用いて、周期性外乱オブザーバの逆モデルを自動調整する機能を持たせることで、多慣性共振系システムヘの実装も容易になる特徴を持つ。
【００２８】
図１１は周期性外乱オブザーバ補償法の演算ブロック図である。同図中、ｄＩAn*：ｎ次周期性外乱電流実部（余弦係数）指令値、ｄＩBn*：ｎ次周期性外乱電流虚部（正弦係数）指令値、ｉｑｃｎ：ｎ次補償電流、ｄＩAn：ｎ次周期性外乱の実部（余弦）成分推定値、ｄＩBn：ｎ次周期性外乱の虚部（正弦）成分推定値、ＩAn：ｎ次補償電流実部（余弦係数）、ＩBn：ｎ次補償電流虚部（正弦係数）。
【００２９】
同図は、トルクリプル周波数に同期した制御系のみを示しているので、システムの特性は１次元複素ベクトルで示される。すなわち、以下の（７）式で実システムを表現し、そのシステム同定結果を（８）式とする。なお、（８）式中の推定値には「＾」記号をＰの頂部に付して示すが、明細書中では「Ｐ＾」のように記す。
【００３０】
【数４】
0005488043.tif 000005
【００３１】
ＰAn：ｎ次トルクリプル周波数成分のシステム実部、ＰBn：ｎ次トルクリプル周波数成分のシステム虚部、Ｐ＾An：ｎ次トルクリプル周波数成分のシステム実部推定値、Ｐ＾Bn：ｎ次トルクリプル周波数成分のシステム虚部推定値
図中のトルクリプル抽出部のローパスフィルタ伝達関数は前記（５）式より、以下の（９）式となる。
【００３２】
【数５】
0005488043.tif 000006
【００３３】
図１１は、一般的な従来の外乱オブザーバと同じ構成であるが、周期性外乱のみに着目しているため、システム特性は（７）式のように１次元複素ベクトルで表され、オブザーハ部の逆システムＰ＾ｓｙｓ-1はシステム同定結果を用いて（８）式の逆数で簡単に示すことができる。
【００３４】
【数６】
0005488043.tif 000007
【００３５】
周期性外乱オブザーバは、上記の複素ベクトル演算した後に、周期性外乱（トルクリプル）の実部成分と虚部成分をそれぞれ推定する。周期性外乱推定値実部ｄＩAnおよび虚部ｄＩBnは、図１１のごとく、その外乱成分を打ち消すように補償電流として与えればよい。通常は、周期性外乱電流指令値ＩAn*、ＩBn*をゼロとすることで、その周波数成分の外乱電流を抑制できる。
【先行技術文献】
【非特許文献】
【００３６】
【非特許文献１】只野  他、「ＰＭモータの周期性外乱に着目したトルクリプル抑制制御法の検討」、平成２１年電気学会産業応用部門大会、Ｉ−６１５〜６１８、会期：平成２１年８月３１日〜９月２日、会場：三重大学
【発明の概要】
【発明が解決しようとする課題】
【００３７】
複数の周波数成分の正弦波トルク指令で加振・重畳して任意トルク波形を生成するモータの制御装置において、トルク加振周波数が電気機械系システムの共振点に一致すると、トルク波形が共振して振動・騒音・機器破損・計測障害などの種々の問題を引き起こす。また、モータのトルク周期性外乱も機械共振によって同様の問題を引き起こす。
【００３８】
非特許文献１に記載される補償電流フーリエ係数学習制御法および周期性外乱オブザーバ補償法では、基本的にトルクリプル外乱を抑制できるが、上記のトルク加振などでシステムにトルク共振が発生してしまう。すなわち、周期性外乱の抑制には対応しているが、トルク加振などで加えられるトルク指令値による電気・機械系の共振現象の抑制には対応していない。
【００３９】
さらに、補償電流フーリエ係数学習制御法では、オンラインでトルクリプルを抑制可能とするには補償電流をテーブル化する必要があるが、この補償電流テーブル方式ではシステムの制御パラメータの変動には対応できない。
【００４０】
本発明の目的は、トルク加振など、モータに任意波形のトルク出力を発生させるモータ制御装置におけるモータのトルク出力の周期性外乱および電気・機械共振をオンラインで抑制できるモータのトルク制御装置を提供することにある。
【課題を解決するための手段】
【００４１】
本発明は、前記の課題を解決するため、任意波形のトルク指令値でモータの出力を制御するモータのトルク制御装置において、モータのトルク出力からトルク周期性外乱をオブザーバで直接に推定してトルク指令値補正してトルク周期性外乱を抑制し、トルク出力の周波数成分による電気・機械系の共振現象を、トルク指令値からトルク検出値までの周波数成分別のシステム伝達関数として同定し、この周波数成分別の同定結果の逆数を用いた共振抑制テーブルまたは共振抑制演算でトルク指令値を生成する。
【００４２】
また、任意波形のトルク指令値を直流分のトルク成分と直流以外のトルク成分に分離し、直流分のトルク成分はｄ軸ｑ軸電流指令値に変換し、直流以外のトルク成分はオブザーバが周期性外乱を複素ベクトルで表現して周期性外乱を求める場合の実部成分と虚部成分に分解して該オブザーバに周期性外乱電流指令値として入力する。
【００４３】
以上のことから、本発明は、以下の構成を特徴とする。
【００４４】
（１）負荷装置を駆動するモータを、任意波形のトルク指令値に出力制御するトルク制御回路を設けたモータのトルク制御装置において、
モータのトルク出力からトルク周期性外乱をオブザーバで直接に推定し、この推定値でトルク指令値を補正してトルク周期性外乱を抑制する周期性外乱抑制制御手段と、
前記トルク出力の周波数成分による電気・機械系の共振現象を、前記トルク指令値からトルク検出値までの周波数成分別のシステム伝達関数として同定し、この周波数成分別の同定結果の逆数を用いた共振抑制テーブルまたは共振抑制演算でトルク指令値を生成する機械共振抑制制御手段を備えたことを特徴とする。
【００４５】
（２）前記機械共振抑制制御手段は、前記システム同定結果を、システム特性を近似的な数式で表現することを特徴とする。
【００４６】
（３）前記トルク制御回路は、前記トルク指令値をベクトル制御における回転座標系のｄ軸ｑ軸電流指令値に変換してモータ電流を制御する構成とし、
前記周期性外乱抑制制御手段は、前記ｑ軸電流指令値から卜ルク検出値までのシステム同定結果に基づいて、その逆数を用いて周波数成分ごとに前記ｑ軸電流指令値のみのトルク指令値を生成することを特徴とする。
【００４７】
（４）前記機械共振抑制制御手段は、直流成分を含む複数の周波数成分を有する任意波形のトルク指令値を直流分のトルク成分と直流以外のトルク成分に分離し、直流分のトルク指令値はｄ軸ｑ軸電流指令値の直流分として生成し、直流以外のトルク指令値は前記共振抑制テーブルまたは共振抑制演算で共振抑制を考慮したｑ軸電流指令値に変換することを特徴とする。
【００４８】
（５）前記機械共振抑制制御手段は、直流成分を含む複数の周波数成分を有する任意波形のトルク指令値を直流分のトルク成分と直流以外のトルク成分に分離し、直流分のトルク成分はｄ軸ｑ軸電流指令値に変換し、直流以外のトルク成分は前記共振抑制テーブルまたは共振抑制演算により機械共振特性を考慮したトルク指令値に変換することを特徴とする。
【００４９】
（６）前記機械共振抑制制御手段は、直流成分を含む複数の周波数成分を有する任意波形のトルク指令値を直流分のトルク成分と直流以外のトルク成分に分離し、直流分のトルク成分はｄ軸ｑ軸電流指令値に変換し、直流以外のトルク成分は前記オブザーバが周期性外乱を複素ベクトルで表現して周期性外乱を求める場合の実部成分ｄＩAn*（余弦係数成分）と虚部成分ｄＩBn*（正弦係数成分）に分解して該オブザーバに周期性外乱電流指令値として入力することを特徴とする。
【００５０】
（７）前記機械共振抑制制御手段は、直流成分を含む複数の周波数成分を有する任意波形のトルク指令値を直流分のトルク成分と直流以外のトルク成分に分離し、直流分のトルク成分はｄ軸ｑ軸電流指令値に変換し、直流以外のトルク成分から前記トルクリプル周期性外乱と同期する周波数成分のみを抽出して前記オブザーバが周期性外乱を複素ベクトルで表現して周期性外乱を求める場合の実部成分ｄＩAn*（余弦係数成分）と虚部成分ｄＩBn*（正弦係数成分）に分解して該オブザーバに周期性外乱電流指令値として入力することを特徴とする。
【００５１】
（８）前記機械共振抑制制御手段および前記周期性外乱抑制制御手段は、直流成分を含む複数の周波数成分を有する任意波形のトルク指令値を重ね合わせて任意波形のトルク出力を発生させる場合、並列かつ独立に任意周波数成分別のトルク指令値、複素ベクトル表現の電流指令値を生成することを特徴とする。
【発明の効果】
【００５２】
以上のとおり、本発明によれば、モータのトルク出力からトルク周期性外乱をオブザーバで直接に推定してトルク指令値補正してトルク周期性外乱を抑制し、トルク出力の周波数成分による電気・機械系の共振現象を、トルク指令値からトルク検出値までの周波数成分別のシステム伝達関数として同定し、この周波数成分別の同定結果の逆数を用いた共振抑制テーブルまたは共振抑制演算でトルク指令値を生成するようにしたため、トルク加振装置など、モータに任意波形のトルク出力を発生させるモータ制御装置におけるモータのトルク出力の周期性外乱抑制および機械共振抑制をオンラインでできる。
【００５３】
また、任意波形のトルク指令値を直流分のトルク成分と直流以外のトルク成分に分離し、直流分のトルク成分はｄ軸ｑ軸電流指令値に変換し、直流以外のトルク成分はオブザーバが周期性外乱を複素ベクトルで表現して周期性外乱を求める場合の実部成分と虚部成分に分解して該オブザーバに周期性外乱電流指令値として入力するようにしたため、トルクリプル周波数とトルク指令周波数が一致した場合の干渉による不安定動作を解消できる。
【図面の簡単な説明】
【００５４】
【図１】トルク制御装置の制御ブロック図（実施形態１）。
【図２】周期性外乱オブザーバのブロック図。
【図３】トルク制御装置の制御ブロック図（実施形態２）。
【図４】トルク制御装置の制御ブロック図（実施形態３）。
【図５】周期性外乱オブザーバのブロック図。
【図６】トルク制御装置の制御ブロック図（実施形態４）。
【図７】トルク制御装置の制御ブロック図（実施形態５）。
【図８】従来のトルク制御装置の構成図。
【図９】従来のトルク制御装置の制御ブロック図。
【図１０】トルク制御装置の周波数伝達関数の例。
【図１１】周期性外乱オブザーバによるシステム同定のブロック図。
【発明を実施するための形態】
【００５５】
（実施形態１）
本実施形態の基本構成図を図１に示す。図２は、図１におけるトルクリプル抑制制御部の代表形態である周期性外乱オブザーバの演算ブロック構成を示す。
【００５６】
本実施形態は、図１が図９と異なる部分は、図２に示す周期性外乱オブザーバによるトルクの周期性外乱抑制を行うと同時に、電気・機械系の共振を抑制する機械共振抑制制御手段として軸トルク指令Ｔｓｈ*から機械共振特性を考慮したトルク指令値Ｔ*に変換する共振抑制テーブル２３を設けた点にある。
【００５７】
ここで、共振抑制テーブル２３の実現方法を説明する。共振抑制テーブル２３は、システム同定結果に基づいて構成する（同定結果でなく、システム特性を近似式等で表現してもよいが、本実施形態では精密に得られた同定結果を用いる例を考える）、同定入力「トルク指令値Ｔ*」から同定出力「トルク検出値Ｔdet」までの周波数成分別のシステム伝達関数として同定し、この周波数成分別の同定結果の逆数を用いて構成する。システム同定は一般的な技術であるため、手法の詳細は特に限定しない。例えば、トルク指令値Ｔ*にホワイトノイズを同定入力信号して与え、そのときの出力信号（Ｔdet）を観測して記録する。入力・出力の記録データから周波数伝達関数を求めれば、周波数成分ごとのシステム伝達特性（同定結果）が入手できる。
【００５８】
同定結果は下記の（１２）式のとおり、周波数成分ごとに１次元複素ベクトルで表現してテーブル化する。ｍは周波数成分を表すテーブル要素番号、Ｎはテーブルの要素数を示している。実際に使用するときは所望の加振周波数成分に対応する１次元複素ベクトルのみをテーブルから抽出して使用する。
【００５９】
【数７】
0005488043.tif 000008
【００６０】
（１２）式中、Ｎはテーブルの周波数分解能を決定づける任意自然数で、例えばＮ＝１０００個でＩＨｚごとにテーブルを構成した場合は、１〜１０００Ｈｚまでを１Ｈｚ刻みで表現できる。ＰAm：システム同定結果の実部ｍ番目要素、ＰBm：システム同定結果の虚部ｍ番目要素。
【００６１】
上記の（１２）式で表現したＰTmはｍ番目要素に対応する周波数成分におけるトルク指令Ｔ*からトルク検出値Ｔdetまでの伝達特性を示しているので、その逆算値、すなわち逆特性ＰTm-1を介して所望の加振トルク指令値Ｔｓｈ*を入力すれば、システム共振特性を予め考慮したトルク指令値Ｔ*を生成することができる。（１２）式の逆数であるＰTm-1は、１次元複素ベクトルの逆数の集合データとなるので容易に計算できる。
【００６２】
したがって、図１の共振抑制テーブル２３には、予め計算したシステム逆特性ＰTm-1のテーブルを与える。共振抑制テーブル２３から抽出する周波数成分は、加振トルク指令値の周波数成分に一致させる。なお、上記の共振抑制テーブル２３は回転速度ωをパラメータとした関数式で表現し、この関数式に従った演算でトルク指令Ｔ*をオンラインで求めることも可能である。
【００６３】
トルクリプル抑制制御部２１については、前記の周期性外乱オブザーバ補償法と同等の構成であり、外乱電流指令値ｄＩAn＝０、ｄＩBn＝０としておけば、周期性外乱をオンラインで学習して抑制できる。
【００６４】
本実施形態によれば、加振トルク指令が電気−機械系の共振点に一致しても、危険なトルク共振現象を発生することなく、安全かつ良好に軸トルク波形を所望値に制御することができる。さらに、トルクの周期性外乱についてもオブザーバによってオンラインで抑制制御を同時に行う構成とすることで、周期性のある外乱抑制・目標値応答の双方を同時に実現できる。
【００６５】
一般に多慣性共振系システムは高次の複雑なシステムとなり、制御器も複雑な演算を要求される場合が多い。本実施形態を用いれば、逆モデル特性も１次元複素ベクトルで表されるため、比較的容易に抑制制御系を設計できる利点がある。
【００６６】
（実施形態２）
モータが埋込磁石同期電動機（ＩＰＭＳＭ：Ｉｎｔｅｒｉｏｒ  Ｐｅｒｍａｎｅｎｔ  Ｍａｇｎｅｔ  Ｓｙｎｃｈｒｏｎｏｕｓ  Ｍｏｔｏｒ）の場合は、ｄ軸電流を用いてリラクタンストルクを有効活用するのが通例である。例えば最大トルクが得られるようにｄ軸電流とｑ軸電流の割合を制御する。一般には理論的な変換数式や、調整可能な変換テーブルを用いてトルク指令からｉｄ、ｉｑ電流指令値に変換している。
【００６７】
上記の実施形態では、加振周波数成分を含むトルク指令をｉｄ、ｉｑ電流指令値に変換してベクトル制御を行っている。つまり、ｉｄ、ｉｑ電流指令の双方が加振周波数成分を有することになる。ＩＰＭＳＭは構造上、負荷や回転数に依存して複雑なトルク発生メカニズムを持っており、トルク誤差を発生しやすい。そのようなモータに対してｄ軸、ｑ軸電流を加振した場合、上述の誤差やトルクリプル特性などの影響を受けて、より複雑な周波数成分のトルクを発生してしまう可能性がある。
【００６８】
ＩＰＭＳＭの最も単純なトルク式は（１３）式で示されるが、右辺ではｉｄとｉｑの乗算が行われるため、双方に加振周波数成分があると２倍の周波数を発生する。また、実際にはトルクリプルがあるので．（１３）式のような単純な近似式ではトルクに含まれる周波数成分を正確に表現できない。
【００６９】
【数８】
0005488043.tif 000009
【００７０】
Ｔ：トルク、Φ：永久磁石による鎖交磁束数、Ｌｄ：ｄ軸インダクタンス、Ｌｑ：ｑ軸インダクタンス、ｐ：極対数
このように、上述する種々の原因から、モータが発生するトルクには複雑な周波数成分を含まれている。そこで、本実施形態では、ｄ軸とｑ軸電流の加振周波数の誤差と干渉による悪影響を低減するため、ｑ軸電流のみを用いて加振する。図３は、本実施形態の基本構成図であり、図１と異なる部分は、軸トルク指令分離部２４を追加した点にある。
【００７１】
軸トルク指令分離部２４は、任意波形（すなわち複数の周波数成分を含む）の軸トルク指令値Ｔｓｈ*から、直流成分Ｔ0*と加振成分Ｔｖｉｂ*（直流以外の複数の周波数成分）を分離する。分離方法は任意である。例えば、任意波形の指令値を生成する際、予め直流分とそれ以外の周波数成分の指令値を分けて指令してもよいし、混合された指令値であっても高域通過フィルタ（ハイパスフィルタ）を用いて直流以外の成分を分離することができる。
【００７２】
上記で分離したトルク指令の直流成分Ｔ0*は、通常のトルク−ｉｄ、ｉｑ電流変換部１１を介してｄ軸電流指令値ｉｄ*、ｑ軸電流指令値（直流分）ｉｑｏ*に変換する。一方、加振トルク指令値Ｔｖｉｂ*は、共振抑制テーブル２３Ａを用いて、加振トルク指令値による共振現象を抑制する。ただし、ここで用いる共振抑制テーブル２３Ａは前述の実施形態とは異なる。
【００７３】
本実施形態では、加振トルク指令からｑ軸電流指令値に変換するため、下記の（１４）式に示す加振電流指令値ｉｑｖｉｂ*からトルク検出値Ｔdetまでのシステム同定結果Ｐiqを用いて、その逆数Ｐiq-1（逆モデル）を共振抑制テーブルとして与える。（基本原理説明は実施形態１に示した通り。）
【００７４】
【数９】
0005488043.tif 000010
【００７５】
ＰAiqm：システム同定結果の実部ｍ番目要素、ＰBiqm：システム同定結果の虚部ｍ番目要素
なお、トルクリプル周期性外乱抑制については、実施形態１と同様に、外乱電流指令値ｄｌAn、ｄＩBnをゼロとしておけば周期性外乱をオンラインで学習して抑制できる。そして、図３に示すとおり、共振抑制を考慮した加振電流指令値ｉｑｖｉｂ*とトルクリプル補償電流ｉｑｃ*は、共にｑ軸電流指令値に重畳してｉｑ*を生成する。
【００７６】
本実施形態によれば、直流トルク成分と分離しつつ、ｑ軸電流指令によるトルク加振を用いることで、トルクリプルやｄ軸・ｑ軸電流加振の干渉によるトルク波形ひずみやトルク誤差の悪影響を低減することができる。
【００７７】
（実施形態３）
上記の実施形態１，２ではオンラインでトルクリプル周期性外乱を抑制すると同時に、加振トルク指令値に対しても電気・機械共振を起こさないシステムを実現できる。しかしながら、トルクリプル周波数と加振トルク指令周波数が一致した場合、加振したい周波数成分まで周期性外乱と見なしてトルクリプル抑制制御系が機能してしまう。そのため、上記の条件では所望の卜ルクが得られず、不安定な動作となる可能性がある。
【００７８】
本実施形態以降の実施形態では、上記のトルクリプル抑制制御系と加振トルク共振抑制系の干渉問題を解決できるトルク制御装置を提案するものである。
【００７９】
本実施形態の基本構成図を図４に示し、図５には図４におけるトルクリプル抑制制御部２１の代表形態である周期性外乱オブザーバの演算ブロックを示す。図４において、加振周波数成分抽出部２５は、任意波形（すなわち複数の周波数成分を含む）の軸トルク指令値Ｔｓｈ*から、直流成分Ｔ0*と加振成分Ｔｖｉｂ*（直流以外の複数の周波数成分）を分離する。分離方法は任意である。例えば、任意波形の指令値を生成する際、予め直流分とそれ以外の周波数成分の指令値を分けて指令してもよいし、混合された指令値であっても高域通過フィルタ（ハイパスフィルタ）を用いる、あるいはフーリエ変換や（３）〜（５）式で示したような周波数成分抽出方法を用いて直流以外の成分を分離することができる。
【００８０】
上記で分離した軸トルク指令値の直流成分Ｔ0*は、トルク−ｉｄ、ｉｑ電流変換部１１を介してｄ軸電流指令値ｉｄ*、ｑ軸電流指令値（直流分）ｉｑｏ*に変換する。一方、加振トルク指令値Ｔｖｉｂ*は、共振抑制テーブル２３Ｂを用いて、加振トルク指令値による共振現象を抑制する。ただし、ここで用いる共振抑制テーブル２３Ｂは前述の実施形態とは異なり、軸トルク指令値Ｔｓｈ*から機械共振特性を考慮したトルク指令値に変換する。
【００８１】
共振抑制テーブル２３Ｂは、システム同定結果に基づいて構成し（同定結果でなく、システム特性を示す近似式等でもよいが、本実施形態では精密に得られた同定結果を用いる例を考える）、同定入力「補償電流指令値ｉｑｃ*」から同定出力「トルク検出値Ｔdet」までの周波数成分別のシステム伝達関数として同定し、この周波数成分別の同定結果の逆数を用いて構成する。システム同定は一般的な技術であるため、手法の詳細は特に限定しない。例えば、補償電流指令値ｉｑｃ*には同定入力信号してホワイトノイズを与え、そのときの出力信号（Ｔdet）を観測して記録する。入力・出力の記録データから周波数伝達関数を求めれば、周波数成分ごとのシステム伝達特性（同定結果）が入手できる。
【００８２】
同定結果は下記の（１５）式のとおり、周波数成分ごとに１次元複素ベクトルで表現してテーブル化する。ｍは周波数成分を表すテーブル要素番号、Ｎはテーブルの要素数を示している。実際に使用するときは所望の加振周波数成分に対応する１次元複素ベクトルのみをテーブルから抽出して使用する。
【００８３】
【数１０】
0005488043.tif 000011
【００８４】
（１５）式中、Ｎはテーブルの周波数分解能を決定づける任意自然数で、例えばＮ＝１０００個で１Ｈｚごとにテーブルを構成した場合は、１〜１０００Ｈｚまでを１Ｈｚ刻みで表現できる。ＰAiqm：システム同定結果の実部ｍ番目要素、ＰBiqm：システム同定結果の虚部ｍ番目要素
上記の（１５）式で表現したＰiqは、ｍ番目要素に対応する周波数成分における補償電流指令ｉｑｃ*からトルク検出値Ｔdetまでの伝達特性を示しているので、その逆算値、すなわち逆特性Ｐiq-1を介して所望の加振トルク指令値Ｔｓｈ*を入力すれば、システム共振特性を予め考慮したトルク指令値を生成することができる。（１５）式の逆数であるＰiq-1は、１次元複素ベクトルの逆数の集合データとなるので容易に計算できる。
【００８５】
したがって、図４の共振抑制テーブル２３Ｂには、予め計算したシステム逆特性Ｐiq-1のテーブルを与える。共振抑制テーブル２３Ｂから抽出する周波数成分は、加振トルク指令値の周波数成分に一致させる。なお、上記共振抑制テーブルは回転速度ωをパラメータとした関数式で表現することも可能である。
【００８６】
上記で分離・抽出したトルク指令の直流成分Ｔ0*は、通常のトルク−ｉｄ、ｉｑ電流変換部１１を介してｄ軸電流指令値ｉｄ*、ｑ軸電流指令値（直流分）ｉｑ0*に変換する。一方、加振成分Ｔｖｉｂ*は、上述の共振抑制テーブル２３Ｂを用いて機械共振を考慮した加振電流指令値ｉｑｖｉｂ*に変換する。図３に示す実施形態２の手法では、加振電流指令値ｉｑｖｉｂ*をそのままｑ軸電流指令値に重畳していたが、本実施形態では実部・虚部分離部２６により、任意周波数成分の加振電流指令値を複素ベクトルして表現した場合の実部成分ｄＩAn*（余弦係数成分）と虚部成分ｄＩBn*（正弦係数成分）に分解し、トルクリプル抑制制御部２１に周期性外乱電流指令値として入力する。つまり、任意ｎ次周波数成分の加振電流指令値が下記の（１６）式で示されるとき、ｉｑAn＝ｄＩAn*、ｉｑBn＝ｄｌBn*として周期性外乱電流指令値を与える。トルクリプル抑制制御部２１については、実施形態１，２で提案した周期性外乱オブザーバと同等である。
【００８７】
【数１１】
0005488043.tif 000012
【００８８】
上記で説明したように、加振電流指令値ｉｑｖｉｂ*を周期性外乱電流指令値ｉｑAn＝ｄＩAn*、ｉｑBn＝ｄｌBn*とすることで、周期性外乱オブザーバによるトルクリプル補償は所望の加振電流を生成しつつ、加振成分以外の外乱を抑制するように機能する。外乱抑制された加振電流指令値の実部成分ＩAnおよび虚部成分ＩBnは補償電流生成部２２で下記の（１６）式に従った演算で補償電流指令値ｉｑｃ＊に変換し、ｑ軸電流指令値ｉｑ０*に重畳する。
【００８９】
本実施形態によれば、加振トルク周波数成分とトルクリプル周波数成分が一致する条件であっても、加振トルク共振抑制制御系とトルクリプル抑制制御系の干渉を回避して共振現象を抑制し、軸トルクを所望の値に制御することができる。
【００９０】
（実施形態４）
実施形態３は、軸トルク加振周波数成分が予め分かっている状態で、その加振電流指令値の実部と虚部を分離してトルクリプル抑制制御系の周期性外乱電流指令値としていた。本実施形態では、軸トルク指令値の周波数成分が分からない場合にもトルクリプル抑制制御系との干渉問題を回避するシステムを提案する。
【００９１】
図６は、本実施形態の基本構成図である。基本的には実施形態３と同様であるが、共振抑制テーブル２３Ｂの出力である加振電流指令値ｉｑｖｉｂ*をそのままｑ軸電流指令値に重畳するとともに、トルク脈動同期周波数成分抽出部２７では加振電流指令値ｉｑｖｉｂ*からトルクリプル周期性外乱と同期する周波数成分のみを抽出してトルクリプル抑制制御の周期性外乱電流指令値ｄＩAn*、ｄｌBn*とする。
【００９２】
トルク脈動同期周波数成分抽出部２７は、基本的にトルク脈動抽出部と同じ処理を行えばよい。すなわち、回転位相θを用いて以下のように抽出する。
【００９３】
【数１２】
0005488043.tif 000013
【００９４】
£：ラプラス変換、ωf：脈動抽出カットオフ周波数、ｉｑAn：加振電流指令値の実部成分（余弦係数）、ｉｑBn：加振電流指令値の虚部成分（正弦係数）
上記の（１７）〜（１９）式で求めた加振電流指令値の実部・虚部成分は、そのままトルクリプル抑制制御系の周期性外乱電流指令値の実部・虚部成分として与えればよい。すなわち、ｄＩAn*＝ｉｑAn、ｄｌBn*＝ｉｑBnとすればよい。
【００９５】
本実施形態によれば、周波数成分が分からない、あるいは非周期的な周波数成分を含む加振トルク指令値が与えられても、基本的に加振トルク指令値による共振抑制を考慮した電流指令値を生成できる。また、トルクリプル周波数と同期する瞬間があっても、そのときの加振電流の周波数成分を抽出して周期性外乱電流指令値を与えるため、抑制制御系が干渉する問題も同時に回避することができる。
【００９６】
（実施形態５）
本実施形態では、実施形態３あるいは実施形態４の周波数成分を並列かつ同時に実現する方法を提案する。基本構成図を図７に示し、加振周波数成分抽出部２５Ａ、共振抑制テーブル２３Ｃ、トルク脈動同期周波数成分抽出部２７Ａ、トルクリプル抑制制御部２１Ａ、トルク脈動周波数成分抽出部１９Ａ、補償電流生成部２２Ａでは周波数成分別の指令値やトルク脈動等を求める。また、加算器２８では周波数成分の加振電流指令値ｉｑｖｉｂｌ*〜ｉｑｖｉｂｎ*の加算によって加振電流指令値ｉｑｖｉｂ*を求める。
【００９７】
図７中の各記号は、Ｔｖｉｂｌ*〜Ｔｖｉｂｎ*は異なる任意周波数成分のトルク加振指令値、ｉｑｖｉｂｌ*〜ｉｑｖｉｂｎ*は異なる任意周波数成分の加振電流指令値、ｄｉA1*〜ｄｉAn*は異なる任意周波数成分の周期性外乱電流指令値の実部成分（余弦係数）、ｄｉB1*〜ｄｉBn*は異なる任意周波数成分の周期性外乱電流指令値の虚部成分（正弦係数）、ＴA1*〜ＴAn*は異なる任意周波数成分の周期性外乱トルクリプルの実部成分（余弦係数）、ＴB1*〜ＴBn*は異なる任意周波数成分の周期性外乱トルクリプルの虚部成分（正弦係数）、ＩA1*〜ＩAn*は周期性外乱抑制された加振電流指令値の実部成分（余弦係数）、ＩB1*〜ＩBn*は周期性外乱抑制された加振電流指令値の虚部成分（正弦係数）、ｉｑｃ*は周期性外乱補償された加振電流指令値である。
【００９８】
任意波形は、複数の周波数成分の重ね合わせで実現できる。そこで、本実施形態では、実施形態３または実施形態４の構成を複数の周波数成分のトルク加振指令値に対応させる。ここでは代表形態として実施形態４の構成を並列化する。トルクリプル周期性外乱抑制と加振信号による共振現象の抑制を両立する方法は、実施形態４に記載の手法と同様である。
【００９９】
本実施形態では、図７に示すとおり、軸トルク指令値が複数の周波数成分を持つ場合を想定しているので、その周波数成分Ｔｖｉｂｌ*〜Ｔｖｉｂｎ*を抽出して共振抑制テーブル２３Ｃに入力し、このテーブルには共振特性を考慮した加振電流指令値ｉｑｖｉｂｌ*〜ｉｑｖｉｂｎ*が得られるので、それらを加算器２８で加算した加振電流指令値ｉｑｖｉｂ*をｑ軸電流指令値に重畳する。次に、トルクリプル抑制制御系との干渉を防止するために、トルク脈動同期周波数成分抽出部１９Ａでｉｑｖｉｂ*に含まれるトルク脈動同期成分を抽出する。抽出方法は実施形態４と同様で、例えば前記の（１７）〜（１９）式を用いればよい。抽出された実部・虚部電流信号はそのまま周期性外乱電流指令値ｄｉA1*〜ｄｉAn*、ｄｉA1*〜ｄｉAn*としてトルクリプル抑制制御系に入力する。周期性外乱補償電流ＩA1*〜ＩAn*、ＩB1*〜ＩBn*は、前記の（２）式を用いてそれぞれの次数で生成・加算したものを補償電流ｉｑｃ*とし、ｑ軸電流指令値に重畳する。そのほかの構成部については、基本的に実施形態４と同様である。
【０１００】
以上、本実施形態を用いれば、複数の周波数成分の加振トルクを重ね合わせて、共振現象を発生することなく任意トルク波形を生成することが可能となる。併せて、トルクリプル周期性外乱を同時に抑制することができる。
【０１０１】
なお、各実施形態においては、モータ出力にトルク加振を得る制御装置の場合を示すが、モータに任意波形のトルク出力を発生させるモータの制御装置に適用して同等の作用効果を得ることができる。
【符号の説明】
【０１０２】
１１  指令値変換部
１２  電流ベクトル制御部
１３  モータ
１４  負荷装置
１５  座標変換部
１６  電流センサ
１７  回転位相／速度検出部
１８  回転位置センサ
１９，１９Ａ  トルク脈動周波数成分抽出部
２０  軸トルク検出器
２１，２１Ａ  トルクリプル抑制制御部
２２，２２Ａ  補償電流生成部
２３，２３Ａ，２３Ｂ、２３Ｃ  共振抑制テーブル
２４  軸トルク指令分離部
２５，２５Ａ  加振周波数成分抽出部
２６  実部・虚部分離部
２７，２７Ａ  トルク脈動同期周波数成分抽出部
２８  加算器
【図１】
0005488043.tif 000014
【図２】
0005488043.tif 000015
【図３】
0005488043.tif 000016
【図４】
0005488043.tif 000017
【図５】
0005488043.tif 000018
【図６】
0005488043.tif 000019
【図７】
0005488043.tif 000020
【図８】
0005488043.tif 000021
【図９】
0005488043.tif 000022
【図１０】
0005488043.tif 000023
【図１１】
0005488043.tif 000024