(57)【要約】  対象物を位置誤差に従い、目標位置に位置制御するとともに、外乱に応じた外乱適応制御を行う位置制御装置において、外乱の角周波数ωとし、サンプル周期をＴとした場合に、適応則に従い、ｓｉｎ（ωＴ／２）を計算し、外乱周波数を推定し、ｓｉｎ（ωＴ／２）の値により、制御変数を変更して、デジタルの外乱適応制御するので、デジタル制御において、外乱の周波数の推定精度を向上でき、且つ広い周波数範囲の外乱適応制御を精度よく実現できる。
2008142727.tif 000002
【特許請求の範囲】
【請求項１】
対象物の位置に応じて、アクチュエータにより、前記対象物を位置制御する位置制御装置において、
目標位置と前記対象物の現在位置との位置誤差に応じて、内部制御定数を使用して、外乱適応制御を含む位置制御を、所定のサンプル周期で実行するデジタル制御ブロックと、
前記外乱の角周波数ωとし、前記サンプル周期をＴとした場合に、ｓｉｎ（ωＴ／２）の値に応じた前記内部制御定数を格納するテーブルとを有し、
前記デジタル制御ブロックは、前記位置誤差に従い、前記アクチュエータの制御値を求めるとともに、前記位置誤差を基準にした信号から、適応則に従い、前記ｓｉｎ（ωＴ／２）の値を計算し、前記計算したｓｉｎ（ωＴ／２）の値で、前記テーブルから対応する前記内部制御定数を読み出し、前記内部制御定数を更新して、前記外乱適応制御を含む位置制御を行う
ことを特徴とする位置制御装置。
【請求項２】
前記デジタル制御ブロックは、オブザーバ制御により、前記外乱適応制御を含む位置制御を行う
ことを特徴とする請求項１の位置制御装置。
【請求項３】
前記デジタル制御ブロックは、前記オブザーバの状態変数と、推定ゲインと、推定位置誤差とから誤差を計算し、前記ｓｉｎ（ωＴ／２）に加算して、前記ｓｉｎ（ωＴ／２）を逐次更新する
ことを特徴とする請求項２の位置制御装置。
【請求項４】
前記デジタル制御ブロックは、前記オブザーバの現サンプルの状態変数と、次のサンプルの状態変数とから誤差を計算し、前記ｓｉｎ（ωＴ／２）を逐次更新する
ことを特徴とする請求項２の位置制御装置。
【請求項５】
前記デジタル制御ブロックは、前記計算した前記ｓｉｎ（ωＴ／２）の値を、所定の範囲に収めるためにクリップする
ことを特徴とする請求項１の位置制御装置。
【請求項６】
前記デジタル制御ブロックは、
前記位置誤差に応じて、前記アクチュエータの制御値を演算するコントロールブロックと、
前記適応則に従い、前記ｓｉｎ（ωＴ／２）を計算し、前記計算されたｓｉｎ（ωＴ／２）の値により、前記テーブルから読み出された内部制御定数を使用して、外乱適応制御値を計算する外乱適応制御ブロックとを有し、
前記コントロールブロックの制御値と前記外乱適応制御値とで、前記アクチュエータを駆動する
ことを特徴とする請求項１の位置制御装置。
【請求項７】
前記テーブルは、前記推定した前記ｓｉｎ（ωＴ／２）の値に対応した前記オブザーバの状態推定ゲインを格納し、
前記デジタル制御ブロックは、前記内部制御定数として、前記オブザーバの状態推定ゲインを、前記テーブルから読み出された状態推定ゲインで更新する
ことを特徴とする請求項２の位置制御装置。
【請求項８】
前記デジタル制御ブロックは、前記オブザーバの（ｋ＋１）サンプル時の状態変数ｚ１［ｋ＋１］，ｚ２［ｋ＋１］から、０．５サンプル戻った（ｋ＋０．５）サンプル時の状態変数ｚ１［ｋ＋０．５］，ｚ２［ｋ＋０．５］を計算し、前記オブザーバの推定位置誤差とから誤差を計算し、前記ｓｉｎ（ωＴ／２）に加算して、前記ｓｉｎ（ωＴ／２）を逐次更新する
ことを特徴とする請求項３乃至４のいずれかの位置制御装置。
【請求項９】
前記デジタル制御ブロックは、下記式（３８）及び（３９）を計算して、前記ｓｉｎ（ωＴ／２）を求める
ことを特徴とする請求項３の位置制御装置。
【数３８】
2008142727.tif 000071
【数３９】
2008142727.tif 000072
【請求項１０】
前記デジタル制御ブロックは、下記式（４０）乃至（４２）を計算して、前記ｓｉｎ（ωＴ／２）を求める
ことを特徴とする請求項４の位置制御装置。
【数４０】
2008142727.tif 000073
【数４１】
2008142727.tif 000074
【数４２】
2008142727.tif 000075
【請求項１１】
ディスクから情報を読み取るヘッドと、
前記ヘッドを前記ディスクのトラック横断方向に移動するアクチュエータと、
目標位置と前記ヘッドの現在位置との位置誤差に応じて、内部制御定数を使用して、外乱適応制御を含む位置制御を、所定のサンプル周期で実行するデジタル制御ブロックと、
前記外乱の角周波数ωとし、前記サンプル周期をＴとした場合に、ｓｉｎ（ωＴ／２）の値に応じた前記内部制御定数を格納するテーブルとを有し、
前記デジタル制御ブロックは、前記位置誤差に従い、前記アクチュエータの制御値を求めるとともに、前記位置誤差を基準にした信号から、適応則に従い、前記ｓｉｎ（ωＴ／２）の値を計算し、前記計算したｓｉｎ（ωＴ／２）の値で、前記テーブルから対応する前記内部制御定数を読み出し、前記内部制御定数を更新して、前記外乱適応制御を含む位置制御を行う
ことを特徴とするディスク装置。
【請求項１２】
前記デジタル制御ブロックは、オブザーバ制御により、前記外乱適応制御を含む位置制御を行う
ことを特徴とする請求項１１のディスク装置。
【請求項１３】
前記デジタル制御ブロックは、前記オブザーバの状態変数と、推定ゲインと、推定位置誤差とから誤差を計算し、前記ｓｉｎ（ωＴ／２）に加算して、前記ｓｉｎ（ωＴ／２）を逐次更新する
ことを特徴とする請求項１２のディスク装置。
【請求項１４】
前記デジタル制御ブロックは、前記オブザーバの現サンプルの状態変数と、次のサンプルの状態変数とから誤差を計算し、前記ｓｉｎ（ωＴ／２）を逐次更新する
ことを特徴とする請求項１２のディスク装置。
【請求項１５】
前記デジタル制御ブロックは、前記計算した前記ｓｉｎ（ωＴ／２）の値を、所定の範囲に収めるためにクリップする
ことを特徴とする請求項１１のディスク装置。
【請求項１６】
前記デジタル制御ブロックは、
前記位置誤差に応じて、前記アクチュエータの制御値を演算するコントロールブロックと、
前記適応則に従い、前記ｓｉｎ（ωＴ／２）を計算し、前記計算されたｓｉｎ（ωＴ／２）の値により、前記テーブルから読み出された内部制御定数を使用して、外乱適応制御値を計算する外乱適応制御ブロックとを有し、
前記コントロールブロックの制御値と前記外乱適応制御値とで、前記アクチュエータを駆動する
ことを特徴とする請求項１１のディスク装置。
【請求項１７】
前記テーブルは、前記推定した前記ｓｉｎ（ωＴ／２）の値に対応した前記オブザーバの状態推定ゲインを格納し、
前記デジタル制御ブロックは、前記内部制御定数として、前記オブザーバの状態推定ゲインを、前記テーブルから読み出された状態推定ゲインで更新する
ことを特徴とする請求項１２のディスク装置。
【請求項１８】
前記デジタル制御ブロックは、前記オブザーバの（ｋ＋１）サンプル時の状態変数ｚ１［ｋ＋１］，ｚ２［ｋ＋１］から、０．５サンプル戻った（ｋ＋０．５）サンプル時の状態変数ｚ１［ｋ＋０．５］，ｚ２［ｋ＋０．５］を計算し、前記オブザーバの推定位置誤差とから誤差を計算し、前記ｓｉｎ（ωＴ／２）に加算して、前記ｓｉｎ（ωＴ／２）を逐次更新する
ことを特徴とする請求項１３乃至１４のいずれかのディスク装置。
【請求項１９】
前記デジタル制御ブロックは、下記式（３８）及び（３９）を計算して、前記ｓｉｎ（ωＴ／２）を求める
ことを特徴とする請求項１３のディスク装置。
【数３８】
2008142727.tif 000076
【数３９】
2008142727.tif 000077
【請求項２０】
前記デジタル制御ブロックは、下記式（４０）乃至（４２）を計算して、前記ｓｉｎ（ωＴ／２）を求める
ことを特徴とする請求項１４のディスク装置。
【数４０】
2008142727.tif 000078
【数４１】
2008142727.tif 000079
【数４２】
2008142727.tif 000080
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、外乱による対象物の位置ずれを抑制するための位置制御装置及びディスク装置に関し、特に、未知の外乱周波数による位置ずれを抑制するための位置制御装置及びディスク装置に関する。
【背景技術】
【０００２】
近年、対象物を高精度で位置制御する装置は、広く利用されている。例えば、回転する記憶媒体（例えば、ディスク媒体）に、ヘッドで、データを記録、再生するディスク装置では、ヘッドを記憶媒体の所望のトラックに位置決めし、そのトラックに、データをリード／ライトする。このようなディスク装置、例えば、磁気ディスク装置や光ディスク装置において、ヘッドを目標トラックに正確に位置決めすることが、記録密度向上のために極めて重要である。【０００３】
この位置決めを阻害する要因の１つとして、ディスク上の同心円状のサーボ信号の中心が、モータの回転中心と異なることにより生じる偏心がある。この偏心を補正する技術として、オブザーバを利用した制御方法がある（例えば、特許文献１参照）。【０００４】
偏心は、回転周波数の整数倍に同期した正弦波状の位置ゆれを発生させる。前記のオブザーバ制御方法は、このような正弦波状の位置揺れを抑圧し、目標トラックに正確に位置決めすることができる。しかしながら、このような偏心補正は、あらかじめ補正する周波数を事前に把握していなければならない。たとえば、回転数の整数倍、１倍や２倍でなければならない。【０００５】
位置決めを阻害する２つ目の要因として、ディスク装置の外部から加わる振動がある。この振動はさまざまな波形があるが、ここでは、正弦波状の振動に対処する方法について検討する。前記の偏心補正制御を応用することで、回転数の整数倍以外の周波数にも対処することは可能になる。【０００６】
前記の従来の構成は、外乱の周波数が既知であることを前提としている。ところが、外部振動は、どのような振動が加わるのかは、制御系の設計時点では未知であり、事前にその周波数を把握することはできない。それゆえ、なんらかの未知の周波数の検出手段が必要になり、周波数さえ検出できれば、前述の特許文献１に示したような制御方法を用いて、外部振動による位置ゆれを抑圧する（外部振動に適応追従する）ことが可能になる。【０００７】
図１８は、従来の外乱周波数を検出し、所定の周波数の正弦波状の外乱を抑圧し、位置追従するための制御系の構成図である。図１８に示すように、目標位置ｒと現在位置ｙとの位置誤差ｅを演算ブロック１２０で演算し、制御ブロック（Ｃｎ）１２１で制御演算し、制御量Ｕｎを計算し、プラント１２２（例えば、ディスク装置のＶＣＭ）を駆動する。ディスク装置では、磁気ヘッドからのサーボ信号を、復調し、プラント１２２の現在位置ｙを計算し、演算ブロック１２０にフィードバックする。【０００８】
周波数（ω）推定器１２４は、位置誤差ｅと外乱抑圧補償器１２３（Ｃｄ）の内部変数を用いて、外部振動の角周波数ωを推定する。補償用テーブル１２５は、各周波数ωに対応した外部振動抑圧用の補償器１２３（Ｃｄ）の定数を格納する。外乱抑圧補償器１２３（Ｃｄ）は、周波数推定器１２４の角周波数ωにより、補償用テーブル１２５から読み出された定数により内部定数を修正し、位置誤差ｅから外乱抑圧制御量Ｕｄを計算する。加算ブロック１２６は、制御量Ｕｎと外乱抑圧制御量Ｕｄとを加算して、プラント１２２に出力する。【０００９】
このように、従来は、角周波数（外乱周波数）ωを推定し、角周波数ωの値に応じて、補償器Ｃｄの内部定数を修正し、補償器Ｃｄの動作を、広い周波数範囲にわたり、常に最適に保っていた（例えば、特許文献２参照）。【００１０】
この従来の角周波数ωの推定のために、未知の外乱周波数に適応するため、推定した角周波数ωを逐次修正する。このための適応則は、アナログ制御系で、以下のように導出していた。即ち、周期性外乱を補償する補償器の回転ベクトルを、式（１）の連続系で表現する。【００１１】
【数１】
2008142727.tif 000003
【００１２】
式（１）において、ｚ１、ｚ２は、回転ベクトルのＸ軸成分、Ｙ軸成分である。【００１３】
式（１）のｚ１、ｚ２の時間微分を求めると、下記（２）式が得られる。但し、ωは、時間に対して可変である時間関数である。【００１４】
【数２】
2008142727.tif 000004
【００１５】
この回転ベクトルの角度は、ｓｉｎ成分ｚ１，ｃｏｓ成分ｚ２のｔａｎｇｅｎｔであり、角周波数ωは、その微分値であるから、この回転ベクトルの角周波数ωは、次式（３）で求められる。【００１６】
【数３】
2008142727.tif 000005
【００１７】
ここで、補償器は、次の（４）式を利用する。尚、式（４）は、ωが定常値（ｃｏｎｓｔａｎｔ）であることを仮定している。【００１８】
【数４】
2008142727.tif 000006
【００１９】
式（４）において、ｚ１は、外乱のｓｉｎ成分、ｚ２は、外乱のｃｏｓ成分、Ｌ４，Ｌ５は、入力ゲイン、ｅは、位置誤差、ｓは、ラプラス演算子である。【００２０】
式（４）のアナログの補償器の状態方程式を展開し、ｚ１、ｚ２の微分式を求め、ｚ１，ｚ２の微分式を，式（３）に代入すると、次の式（５）が得られる。【００２１】
【数５】
2008142727.tif 000007
【００２２】
推定した未知外乱の角周波数が、正しければ、補償器１２３（Ｃｄ）は、適切に外乱を抑圧できる。結果として、位置誤差ｅ又はオブザーバの推定位置誤差が、ゼロになる。即ち、式（５）において、補償器２３（Ｃｄ）の対象とする角周波数ωと、推定した外乱の角周波数ωとが一致する条件は、式（５）のω‘、即ち、右辺の位置誤差ｅの項が、ゼロであれば、良い。このため、（５）式を、角周波数の時間微分式で表現し、次式（６）の適応則（積分補償則）を得る。【００２３】
【数６】
2008142727.tif 000008
【００２４】
この式を用いて、ωの値を逐次修正する。（６）式を、積分形式に直し、デジタル制御の式で表現して、次式（７）が得られる。【００２５】
【数７】
2008142727.tif 000009
【００２６】
尚、Ｋは、適応ゲインである。上式（７）は、加算形式で適応則を用いている。これを乗算形式で表現することもできる。この適応則を用いると、図１９に示すように、外乱の振動周波数に対し、推定周波数ωは、リニアーな関係が得られる。
【特許文献１】日本特許第３，４６０，７９５号公報
【特許文献２】日本特許公開２００７−００４９４６号公報
【発明の開示】
【発明が解決しようとする課題】
【００２７】
一方、近年、このような制御系をデジタル制御、即ち、プロセッサのデータ処理で実行する方式が採用されている。前述の従来技術では、アナログの式を元に、アナログの適応則を導出し、これをデジタルの適応則に変換していた。このため、式（７）をデジタル制御に用いると、デジタルの外乱適応制御で、外乱を良好に抑圧できない（外乱に位置追従できない）という問題が生じる。即ち、推定角周波数ωを直接導出しているが、これは、アナログの式に基づくからである。【００２８】
従って、本発明の目的は、デジタル処理に適した外乱適応制御を行うための位置制御装置及びディスク装置を提供することにある。【００２９】
又、本発明の他の目的は、デジタル処理に適した外乱適応則を利用して、広い周波数範囲の外乱適応制御するための位置制御装置及びディスク装置を提供することにある。【００３０】
本発明の更に他の目的は、オブザーバ制御に適した外乱適応制御を行い、より正確に外乱による位置ずれを抑圧するための位置制御装置及びディスク装置を提供することにある。
【課題を解決するための手段】
【００３１】
この目的の達成のため、本発明は、対象物の位置に応じて、アクチュエータにより、前記対象物を位置制御する位置制御装置において、目標位置と前記対象物の現在位置との位置誤差に応じて、内部制御定数を使用して、外乱適応制御を含む位置制御を、所定のサンプル周期で実行するデジタル制御ブロックと、前記外乱の角周波数ωとし、前記サンプル周期をＴとした場合に、ｓｉｎ（ωＴ／２）の値に応じた前記内部制御定数を格納するテーブルとを有し、前記デジタル制御ブロックは、前記位置誤差に従い、前記アクチュエータの制御値を求めるとともに、前記位置誤差を基準にした信号から、適応則に従い、前記ｓｉｎ（ωＴ／２）ωＴの値を計算し、前記計算したｓｉｎ（ωＴ／２）の値で、前記テーブルから対応する前記内部制御定数を読み出し、前記内部制御定数を更新して、前記外乱適応制御を含む位置制御を行う。【００３２】
又、本発明のディスク装置は、ディスクから情報を読み取るヘッドと、前記ヘッドを前記ディスクのトラック横断方向に移動するアクチュエータと、目標位置と前記ヘッドの現在位置との位置誤差に応じて、内部制御定数を使用して、外乱適応制御を含む位置制御を、所定のサンプル周期で実行するデジタル制御ブロックと、前記外乱の角周波数ωとし、前記サンプル周期をＴとした場合に、ｓｉｎ（ωＴ／２）の値に応じた前記内部制御定数を格納するテーブルとを有し、前記デジタル制御ブロックは、前記位置誤差に従い、前記アクチュエータの制御値を求めるとともに、前記位置誤差を基準にした信号から、適応則に従い、前記ｓｉｎ（ωＴ／２）の値を計算し、前記計算したｓｉｎ（ωＴ／２）の値で、前記テーブルから対応する前記内部制御定数を読み出し、前記内部制御定数を更新して、前記外乱適応制御を含む位置制御を行う。【００３３】
更に、本発明では、好ましくは、前記デジタル制御ブロックは、オブザーバ制御により、前記外乱適応制御を含む位置制御を行う。【００３４】
更に、本発明では、好ましくは、前記デジタル制御ブロックは、前記オブザーバの状態変数と、推定ゲインと、推定位置誤差とから誤差を計算し、前記ｓｉｎ（ωＴ／２）に加算して、前記ｓｉｎ（ωＴ／２）を逐次更新する。【００３５】
更に、本発明では、好ましくは、前記デジタル制御ブロックは、前記オブザーバの現サンプルの状態変数と、次のサンプルの状態変数とから誤差を計算し、前記ｓｉｎ（ωＴ／２）を逐次更新する。【００３６】
更に、本発明では、好ましくは、前記デジタル制御ブロックは、前記計算した前記ｓｉｎ（ωＴ／２）の値を、所定の範囲に収めるためにクリップする。【００３７】
更に、本発明では、好ましくは、前記デジタル制御ブロックは、前記位置誤差に応じて、前記アクチュエータの制御値を演算するコントロールブロックと、前記適応則に従い、前記ｓｉｎ（ωＴ／２）を計算し、前記計算されたｓｉｎ（ωＴ／２）の値により、前記テーブルから読み出された内部制御定数を使用して、外乱適応制御値を計算する外乱適応制御ブロックとを有し、前記コントロールブロックの制御値と前記外乱適応制御値とで、前記アクチュエータを駆動する。【００３８】
更に、本発明では、好ましくは、前記テーブルは、前記推定した前記ｓｉｎ（ωＴ／２）の値に対応した前記オブザーバの状態推定ゲインを格納し、前記デジタル制御ブロックは、前記内部制御定数として、前記オブザーバの状態推定ゲインを、前記テーブルから読み出された状態推定ゲインで更新する。【００３９】
更に、本発明では、好ましくは、前記デジタル制御ブロックは、前記オブザーバの（ｋ＋１）サンプル時の状態変数ｚ１［ｋ＋１］，ｚ２［ｋ＋１］から、０．５サンプル戻った（ｋ＋０．５）サンプル時の状態変数ｚ１［ｋ＋０．５］，ｚ２［ｋ＋０．５］を計算し、前記オブザーバの推定位置誤差とから誤差を計算し、前記ｓｉｎ（ωＴ／２）に加算して、前記ｓｉｎ（ωＴ／２）を逐次更新する。【００４０】
更に、本発明では、好ましくは、前記デジタル制御ブロックは、下記式（３８）及び（３９）を計算して、前記ｓｉｎ（ωＴ／２）を求める。【００４１】
【数３８】
2008142727.tif 000010
【００４２】
【数３９】
2008142727.tif 000011
【００４３】
更に、本発明では、好ましくは、前記デジタル制御ブロックは、下記式（４０）乃至（４２）を計算して、前記ｓｉｎ（ωＴ／２）を求める。【００４４】
【数４０】
2008142727.tif 000012
【００４５】
【数４１】
2008142727.tif 000013
【００４６】
【数４２】
2008142727.tif 000014
【発明の効果】
【００４７】
外乱の角周波数ωとし、前記サンプル周期をＴとした場合に、ｓｉｎ（ωＴ／２）を計算し、外乱適応制御するので、デジタル制御において、外乱の周波数の推定精度を向上でき、外乱適応制御を精度よく実現できる。
【図面の簡単な説明】
【００４８】
【図１】本発明の一実施の形態のディスク記憶装置の構成図である。
【図２】図１の構成のディスクの位置信号の説明図である。
【図３】図２の位置信号の詳細説明図である。
【図４】図１のヘッドのシーク動作の説明図である。
【図５】本発明の第１の実施の形態の外乱抑圧補償器を付加した位置制御系のブロック図である。
【図６】図５のｓｉｎ（ωＴ／２）推定の説明図である。
【図７】図５のｓｉｎ（ωＴ）と振動周波数との関係図である。
【図８】図６のｓｉｎ（ωＴ／２）と振動周波数との関係図である。
【図９】本発明の第２の実施の形態の外乱抑圧補償器を付加した位置制御系のブロック図である。
【図１０】図９の構成におけるテーブルの構成図である。
【図１１】図９の構成における他のｓｉｎ（ωＴ／２）推定器の構成図である。
【図１２】本発明の第３の実施の形態の外乱抑圧補償器を付加した位置制御系のブロック図である。
【図１３】図１２の構成における他のｓｉｎ（ωＴ／２）推定器の構成図である。
【図１４】本発明の第４の実施の形態の外乱抑圧補償器を付加した位置制御系のブロック図である。
【図１５】本発明の第５の実施の形態の外乱抑圧補償器を付加した位置制御系のブロック図である。
【図１６】本発明の第６の実施の形態の外乱抑圧補償器を付加した位置制御系のブロック図である。
【図１７】図１６のマルチレート制御の説明図である。
【図１８】従来の外乱適応制御を付加した位置制御系のブロック図である。
【図１９】図１８の従来の外乱周波数の推定の説明図である。
【符号の説明】
【００４９】
１  アクチュエータ
２  スピンドルモータの回転軸
３  ヘッド
４  ディスク
５  スピンドルモータ
６  アクチュエータのＶＣＭ駆動回路
７  位置復調回路
８  スピンドルモータの駆動回路
９  バス
１０  データの記録再生回路
１１  ハードディスクコントローラ
１２  ＭＣＵのＲＡＭ
１３  ＭＣＵのＲＯＭ
１４  マイクロコントローラユニット（ＭＣＵ）
１５  ハードディスクコントローラのＲＡＭ
１６  位置信号
２０  位置誤差計算ブロック
２１  コントローラ
２２  プラント（ヘッド、アクチュエータ）
２３  外乱抑圧補償器
２７  ｓｉｎ（ωＴ／２）推定ブロック
２８  テーブル（定数補正ブロック）
２６  加算ブロック
【発明を実施するための最良の形態】
【００５０】
以下、本発明を図面に従い、ディスク装置、位置制御系の第１の実施の形態、位置制御系の第２の実施の形態、位置制御系の第３の実施の形態、位置制御系の第４の実施の形態、位置制御系の第５の実施の形態、位置制御系の第６の実施の形態、他の実施の形態の順で説明するが、本発明は、下記実施の形態に限らず、種々の変形が可能である。【００５１】
［ディスク装置］
図１は、本発明の一実施の形態のディスク装置の構成図、図２は、図１の磁気ディスクの位置信号の配置図、図３は、図１及び図２の磁気ディスクの位置信号の構成図、図４は、図１のヘッド位置制御の説明図である。【００５２】
図１は、ディスク装置として、磁気ディスク装置を示す。図１に示すように、磁気記憶媒体である磁気ディスク４が、スピンドルモータ５の回転軸２に設けられている。スピンドルモータ５は、磁気ディスク４を回転する。アクチュエータ（ＶＣＭ）１は、そのアーム先端に磁気ヘッド３を備え、磁気ヘッド３を磁気ディスク４の半径方向に移動する。【００５３】
アクチュエータ１は、回転軸を中心に回転するボイスコイルモータ（ＶＣＭ）で構成される。図では、磁気ディスク装置に、２枚の磁気ディスク４が搭載され、４つの磁気ヘッド３が、同一のアクチュエータ１で同時に駆動される。【００５４】
磁気ヘッド３は、リード素子と、ライト素子とからなる。磁気ヘッド３は、スライダに、磁気抵抗（ＭＲ）素子を含むリード素子を積層し、その上にライトコイルを含むライト素子を積層して、構成される。【００５５】
位置検出回路７は、磁気ヘッド３が読み取った位置信号（アナログ信号）をデジタル信号に変換する。リード／ライト（Ｒ／Ｗ）回路１０は、磁気ヘッド３の読み取り及び書込みを制御する。スピンドルモータ（ＳＰＭ）駆動回路８は、スピンドルモータ５を駆動する。ボイスコイルモータ（ＶＣＭ）駆動回路６は、ボイスコイルモータ（ＶＣＭ）１に駆動電流を供給し、ＶＣＭ１を駆動する。【００５６】
マイクロコントローラ（ＭＣＵ）１４は、位置検出回路７からのデジタル位置信号から現在位置を検出（復調）し、検出した現在位置と目標位置との誤差に従い、ＶＣＭ駆動指令値を演算する。即ち、位置復調とサーボ制御を行う。リードオンリーメモリ（ＲＯＭ）１３は、ＭＣＵ１４の制御プログラム等を格納する。ランダムアクセスメモリ（ＲＡＭ）１２は、ＭＣＵ１４の処理のためのデータ等を格納する。【００５７】
ハードディスクコントローラ（ＨＤＣ）１１は、サーボ信号のセクタ番号を基準にして，１周内の位置を判断し，データを記録・再生する。バッファ用ランダムアクセスメモリ（ＲＡＭ）１５は、リードデータやライトデータを一時格納する。ＨＤＣ１１は、ＵＳＢ，ＡＴＡやＳＣＳＩ等のインターフェイスＩＦで、ホストと通信する。バス９は、これらを接続する。【００５８】
図２に示すように、磁気ディスク４には、外周から内周に渡り、各トラックにサーボ信号（位置信号）１６が、円周方向に等間隔に配置される。尚、各トラックは、複数のセクタで構成され、図２の実線は、サーボ信号１６の記録位置を示す。図３に示すように、位置信号は，サーボマークＳｅｒｖｏ  Ｍａｒｋと、トラック番号Ｇｒａｙ  Ｃｏｄｅと、インデックスＩｎｄｅｘと、オフセット情報（サーボバースト）ＰｏｓＡ，ＰｏｓＢ，ＰｏｓＣ，ＰｏｓＤとからなる。尚、図３の点線は、トラックセンターを示す。【００５９】
図３に示すトラック番号Ｇｒａｙ  Ｃｏｄｅとオフセット情報ＰｏｓＡ，ＰｏｓＢ，ＰｏｓＣ，ＰｏｓＤを使い，磁気ヘッドの半径方向の位置を検出する。さらに、インデックス信号Ｉｎｄｅｘを元にして，磁気ヘッドの円周方向の位置を把握する。【００６０】
例えば，インデックス信号を検出したときのセクタ番号を０番に設定し、サーボ信号を検出する毎に、カウントアップして、トラックの各セクタのセクタ番号を得る。このサーボ信号のセクタ番号は，データの記録再生を行うときの基準となる。尚、インデックス信号は、１周に１つである、又、インデックス信号の代わりに、セクタ番号を設けることもできる。【００６１】
図４は、図１のＭＣＵ１４が行うアクチュエータのシーク制御例である。図１の位置検出回路７を通じて、ＭＣＵ１４が、アクチュエータの位置を確認して，サーボ演算し、適切な電流をＶＣＭ１に供給する。図４では、あるトラック位置から目標トラック位置へヘッド３を移動するシーク開始時からの制御の遷移と、アクチュエータ１の電流、アクチュエータ(ヘッド）の速度、アクチュエータ(ヘッド）の位置を示す。【００６２】
即ち、シーク制御は、コアース制御、整定制御及びフォローイング制御と遷移することで，目標位置まで移動させることができる。コアース制御は、基本的に速度制御であり、整定制御、フォローイング制御は、基本的に位置制御であり、いずれも、ヘッドの現在位置を検出する必要がある。【００６３】
このような，位置を確認するためには，前述の図２のように、磁気ディスク上にサーボ信号を事前に記録しておく。即ち、図３に示したように、サーボ信号の開始位置を示すサーボマーク，トラック番号を表すグレイコード，インデックス信号，オフセットを示すＰｏｓＡ〜ＰｏｓＤといった信号が記録されている。この信号を磁気ヘッドで読み出し、このサーボ信号を、位置検出回路７が、デジタル値に変換する。【００６４】
［位置制御系の第１の実施の形態］
図５は、本発明の位置制御系（サーボ制御系）の第１の実施の形態のブロック図、図６は、図５のｓｉｎ（ωＴ／２）推定の説明図、図７は、振動周波数と推定値ｓｉｎ（ωＴ）との関係図、図８は、図５の振動周波数と推定値ｓｉｎ（ωＴ／２）との関係図である。【００６５】
図５は、ＭＣＵ１４が実行するサーボ制御系の演算ブロック図である。即ち、目標位置ｒと現在位置ｙとの位置誤差ｅを演算ブロック２０で演算し、制御ブロック（Ｃｎ）２１に入力する。制御ブロック（Ｃｎ）２１は、位置誤差ｅに従い、制御演算し、制御量Ｕｎを計算し、プラント２２であるＶＣＭ１、３を駆動する。プラントの位置は、磁気ヘッド３からのサーボ信号を、復調し、現在位置ｙを計算し、演算ブロック２０にフィードバックする。【００６６】
周波数推定器２７は、位置誤差ｅと外乱抑圧補償器２３（Ｃｄ）の内部変数を用いて、外部振動の角周波数ωのｓｉｎ（ωＴ／２）を推定する。補償用テーブル２８は、角周波数ωのｓｉｎ（ωＴ／２）に対応した外部振動抑圧用の補償器２３（Ｃｄ）の定数を格納する。【００６７】
外乱抑圧補償器２３（Ｃｄ）は、周波数推定器２７の角周波数ωのｓｉｎ（ωＴ／２）により、補償用テーブル２８から読み出された定数により内部定数を修正し、位置誤差ｅから外乱抑圧制御量Ｕｄを計算する。加算ブロック２６は、制御量Ｕｎと外乱抑圧制御量Ｕｄとを加算して、プラント２２に出力する。外乱抑圧制御量とは、外乱に対し、ディスクとヘッドとの位置関係がずれないように制御するための制御量であり、外乱が印加されると、通常ヘッドの位置ずれが生じるため、ヘッドを駆動するＶＣＭの駆動値を適応制御する。【００６８】
ここで、本発明では、周波数推定器２７は、ｓｉｎ（ωＴ／２）を導出している。この理由を説明する。制御ブロック２１、補償器２３、周波数推定器２４を、オブザーバで構成した制御系で説明する。先ず、予測オブザーバの適応則を説明する。【００６９】
予測オブザーバでの補償器２３は、下記式（８）で表される。尚、式（８）において、ｚ１は、外乱のｓｉｎ成分、ｚ２は、外乱のｃｏｓ成分、Ｌ４，Ｌ５は、入力ゲイン、ｅは、位置誤差である。【００７０】
【数８】
2008142727.tif 000015
【００７１】
ここで、外乱モデルが、正弦波の場合には、式（８）のＡ行列は、下記式（９）で表される。【００７２】
【数９】
2008142727.tif 000016
【００７３】
式（８）を展開すると、式（１０）となる。【００７４】
【数１０】
2008142727.tif 000017
【００７５】
式（２）のアナログの時間微分は、デジタルでは、引き算であるので、式（１０）のＺ２［ｋ＋１］の式から、Ｚ２［ｋ］を引き、Ｚ１［ｋ］を掛けると、式（１１）が得られる。【００７６】
【数１１】
2008142727.tif 000018
【００７７】
同様に、式（１０）のＺ１［ｋ＋１］の式から、Ｚ１［ｋ］を引き、Ｚ２［ｋ］を掛けると、式（１）が得られる。【００７８】
【数１２】
2008142727.tif 000019
【００７９】
式（３）と同様に、アナログの微分は、デジタルの引き算であるから、式（１１）と式（１２）の差分を取ると、式（１３）が得られる。【００８０】
【数１３】
2008142727.tif 000020
【００８１】
又、外乱モデルが正弦波である時は、式（９）のように、ａ１１＝ａ２２（＝ｃｏｓ（ωＴ））、ａ１２＝−ａ２２（＝ｓｉｎ（ωＴ））であるから、式（１３）は、下記式（１４）に変形できる。【００８２】
【数１４】
2008142727.tif 000021
【００８３】
式（１４）において、ｚ１とｚ２の二乗和で除算して、式（１５）を得る。【００８４】
【数１５】
2008142727.tif 000022
【００８５】
式（５）と比較すると、式（５）では、外乱として、ωを推定するのに対し、式（１５）では、ｓｉｎ（ωＴ）を推定することを示す。即ち、デジタル制御においては、ｓｉｎωＴを推定する必要があり、適応則は、式（１５）の右辺の第２項を「０」にするように、ｓｉｎωＴを推定することになる。【００８６】
次に、現在オブザーバの適応則を説明する。現在オブザーバは、デジタル制御において、予測オブザーバより広く利用されている。【００８７】
現在オブザーバでの補償器２３は、下記式（１６）で表される。尚、式（１６）において、ｚ１は、外乱のｓｉｎ成分、ｚ２は、外乱のｃｏｓ成分、Ｌ４，Ｌ５は、入力ゲイン、ｅは、位置誤差である。【００８８】
【数１６】
2008142727.tif 000023
【００８９】
式（１６）の２つの式を纏めると、式（１７）が得られる。【００９０】
【数１７】
2008142727.tif 000024
【００９１】
ここで、外乱モデルが、正弦波の場合には、式（１７）のＡ行列は、下記式（１８）で表される。尚、式の簡略化のため、ｃｏｓ（ωＴ），ｓｉｎ（ωＴ）を、「ｃ」、「ｓ」で表すこととする。【００９２】
【数１８】
2008142727.tif 000025
【００９３】
式（１８）のＡ行列の逆行列は、次式（１９）となる。【００９４】
【数１９】
2008142727.tif 000026
【００９５】
式（１７）を、式（１８）のＡ行列で割り、式（１９）のＡ行列の逆行列に置き換えると、式（２０）が得られる。【００９６】
【数２０】
2008142727.tif 000027
【００９７】
式（２０）から、ｋサンプル時のｚ１、ｚ２を、求めると、式（２１）が得られる。【００９８】
【数２１】
2008142727.tif 000028
【００９９】
式（２１）のｚ２［ｋ＋１］を掛け、ｚ２［ｋ］に、ｚ１［ｋ＋１］を掛けると、式（２２）が得られる。この時、ｚ２［ｋ］に、左辺のマイナスを付与している。【０１００】
【数２２】
2008142727.tif 000029
【０１０１】
式（２２）の上下の式の和をとると、式（２３）が得られる。【０１０２】
【数２３】
2008142727.tif 000030
【０１０３】
式（２３）から、現在オブザーバの構成でも、ｓｉｎ（ωＴ）を推定すべきである。即ち、アナログの式のように、外乱周波数又は外乱角周波数ωを推定するのではなく、デジタル制御においては、直接推定すべきは、ｓｉｎ（ωＴ）である。式（２３）から、ｓｉｎ（ωＴ）は、下記式（２４）で得られる。【０１０４】
【数２４】
2008142727.tif 000031
【０１０５】
図７は、振動周波数と推定値ｓｉｎ（ωＴ）との関係図である。推定値ｓｉｎ（ωＴ）は実線で、従来の推定値ωは、一点鎖線で示す。横軸の振動周波数に対し、従来の推定値ωは、直線状の特性を描くが、推定値ｓｉｎ（ωＴ）と比べると、振動周波数に対し、推定周波数のずれがあり、推定精度が低い。このため、デジタル制御では、推定値ωを推定しても、正確に振動周波数を推定できず、振動周波数に適応する位置制御が困難であった。【０１０６】
一方、本願は、基本的に、振動周波数を、ｓｉｎ（ωＴ）で推定しているため、デジタル制御において、正確に振動周波数を推定でき、振動周波数成分に適応する位置制御が実現できる。【０１０７】
又、デジタル制御において、ｓｉｎ（ωＴ）を、振動周波数ｆの推定値に使用することは、デジタル制御のサンプリング周波数をＦｓとすると、次式（２５）の周波数を推定することになる。【０１０８】
【数２５】
2008142727.tif 000032
【０１０９】
このことは、図７に示したように、外乱周波数ｆとｓｉｎ（ωＴ）とが、単調な関係を示す範囲は、０からＦｓ／４の領域であり、Ｆｓ／４以上では、適応追従ができない。例えば、サンプリング周波数Ｆｓを１２ＫＨｚとすると、３ＫＨｚまで推定できる。【０１１０】
一方、使用用途によっては、更に推定周波数範囲を拡大したいとの要望がある。例えば、ナイキスト周波数であるＦｓ／２まで拡張したいとの要望がある。この周波数範囲の拡張方法を説明する。【０１１１】
現在オブザーバの適応則式（２４）の導出で示した、ｚ１［ｋ＋１］，ｚ２［ｋ＋１］は、あくまでも、推定値である。従って、この推定値ｚ１［ｋ＋１］，ｚ２［ｋ＋１］を使って、０．５サンプル分戻した推定値ｚ１［ｋ＋０．５］，ｚ２［ｋ＋０．５］を作成する。【０１１２】
ｚ１［ｋ＋１］，ｚ２［ｋ＋１］は、振動周波数の回転ベクトルのｓｉｎ，ｃｏｓ成分であるから、図６に示すように、ｋ＋１サンプル時点のｚ１［ｋ＋１］，ｚ２［ｋ＋１］を、０．５サンプル戻す、即ち、ｋ＋０．５サンプル時点の回転ベクトルｚ１［ｋ＋０．５］，ｚ２［ｋ＋０．５］を作成することは、回転ベクトルをωＴ／２戻せば良い。従って、ｚ１［ｋ＋１］，ｚ２［ｋ＋１］から、下記式（２６）により、０．５サンプル戻したｚ１［ｋ＋０．５］，ｚ２［ｋ＋０．５］を計算できる。【０１１３】
【数２６】
2008142727.tif 000033
【０１１４】
この式（２６）の左辺と右辺を入れ替えて、更に、ｚ１［ｋ＋１］，ｚ２［ｋ＋１］が、単独で表されるように、変換すると、次式（２７）が得られる。【０１１５】
【数２７】
2008142727.tif 000034
【０１１６】
前述のｚ１［ｋ＋１」、ｚ２［ｋ＋１］の式（１７）に、式（２７）を代入すると、式（２８）が得られる。【０１１７】
【数２８】
2008142727.tif 000035
【０１１８】
更に、式を整理するため、式（２８）の左辺、右辺ともに、下記式（２９）を掛ける。【０１１９】
【数２９】
2008142727.tif 000036
【０１２０】
この時、ｃｏｓ，ｓｉｎの行列間の掛け算では、下記式（３０）が成立する。【０１２１】
【数３０】
2008142727.tif 000037
【０１２２】
式（３０）を利用して、式（２８）の左辺、右辺に、式（２９）を掛けると、式（３１）が得られる。【０１２３】
【数３１】
2008142727.tif 000038
【０１２４】
即ち、ｓｉｎωＴの式（１７）と同様に、ｓｉｎ（ωＴ／２）の関係式が得られた。この式（３１）を、ｓｉｎωＴの導出と同様に展開する。即ち、ｓｉｎ（ωＴ／２），ｃｏｓ（ωＴ／２）を、それぞれ、「ｓ」、「ｃ」と簡略化すると、式（３１）は、式（３２）に変換される。【０１２５】
【数３２】
2008142727.tif 000039
【０１２６】
式（３２）から、ｋサンプル時のｚ１［ｋ］、ｚ２［ｋ］を、求めると、式（３３）が得られる。【０１２７】
【数３３】
2008142727.tif 000040
【０１２８】
式（３３）の上側に、ｚ２［ｋ＋０．５］を掛け、下側に、−ｚ１［ｋ＋０．５］を掛けると、式（３４）が得られる。【０１２９】
【数３４】
2008142727.tif 000041
【０１３０】
式（３４）の上下の式の和をとると、式（３５）が得られる。【０１３１】
【数３５】
2008142727.tif 000042
【０１３２】
式（３５）から、ｓｉｎ（ωＴ／２）は、下記式（３６）で表される。これを適応則に利用する。【０１３３】
【数３６】
2008142727.tif 000043
【０１３４】
式（３６）から、デジタル制御においては、外乱周波数を直接推定するため、ｓｉｎ（ωＴ／２）を推定すればよい。【０１３５】
図８は、振動周波数と推定値ｓｉｎ（ωＴ／２）との関係図である。デジタル制御において、ｓｉｎ（ωＴ／２）を、振動周波数ｆの推定値に使用することは、デジタル制御のサンプリング周波数をＦｓとすると、次式（３７）の周波数を推定することになる。【０１３６】
【数３７】
2008142727.tif 000044
【０１３７】
このことは、図８に示したように、外乱周波数ｆとｓｉｎ（ωＴ／２）とが、単調な関係を示す範囲は、０からＦｓ／２の領域であり、ナイキスト周波数まで推定できる。又、Ｆｓ／２以上では、適応追従ができない。例えば、サンプリング周波数Ｆｓを１２ＫＨｚとすると、６ＫＨｚまで推定できる。【０１３８】
従って、低い周波数側の推定精度を高くして、且つ高い周波数まで、推定できる。又、ディスク装置の使用形態では、外部振動が、通常、低い周波数側で問題になることが多い。例えば、数１００Ｈｚ程度の外部振動が、外部から加わることが多い。そのため、ディスク装置では、低い周波数の振動を精度よく推定できるｓｉｎ（ωＴ／２）推定式を利用することが、都合が良い。【０１３９】
次に、図５のｓｉｎ（ωＴ／２）推定器２７は、次のように、ｓｉｎ（ωＴ／２）を推定して、適応制御する。【０１４０】
第１の適応制御は、正しく、推定できて、制御系が外乱周波数を抑圧できれば、式（３６）のｅ［ｋ］は、「０」になるので、式（３６）の右辺第２項は、「０」となるはずである。【０１４１】
そこで、ｓｉｎ（ωＴ／２）推定の場合には、式（３６）の右辺第２項Ｅ［ｋ］を，次式（３８）で計算する。【０１４２】
【数３８】
2008142727.tif 000045
【０１４３】
次に、式（３８）のＥ［ｋ］に適応ゲインＫを掛け、次式（３９）のように、ｓｉｎ（ωＴ／２）を更新する。【０１４４】
【数３９】
2008142727.tif 000046
【０１４５】
上記のようにして、ｚ１［ｋ＋１］，ｚ２［ｋ＋１］，ｅ［ｋ］，Ｌ４，Ｌ５の５つの値を利用して、ｓｉｎ（ωＴ／２）を、適応的に推定できる。【０１４６】
第２の適応制御の計算方法を説明する。固定小数点方式のプロセッサを使用する場合には、桁数が限られる。上記の、ｚ１［ｋ＋１］，ｚ２［ｋ＋１］，ｅ［ｋ］，Ｌ４，Ｌ５の５つの値を利用する場合に、ｚ１［ｋ＋１］，ｚ２［ｋ＋１］の値と，ｅ［ｋ］，Ｌ４，Ｌ５の値との桁が相違する。そこで、限られた桁数を有効に利用するため、以下の手順で、計算する。【０１４７】
ｓｉｎ（ωＴ）推定の場合には、式（３６）の右辺第１項Ｇ［ｋ］を，次式（４０）で計算する。【０１４８】
【数４０】
2008142727.tif 000047
【０１４９】
次に、現在の設定値であるｓｉｎ（ωＴ／２）と、式（４０）のずれＥ［ｋ］を、次式（４１）で計算する。【０１５０】
【数４１】
2008142727.tif 000048
【０１５１】
この式（４１）のＥ［ｋ］に適応ゲインＫを掛け、次式（４２）のように、ｓｉｎ（ωＴ／２）を更新する。【０１５２】
【数４２】
2008142727.tif 000049
【０１５３】
上記のようにして、ｚ１［ｋ＋０．５］，ｚ１［ｋ］，ｚ２［ｋ＋０．５］，ｚ［ｋ］の４つの値を利用するだけで、ｓｉｎ（ωＴ／２）を、適応的に推定できる。【０１５４】
次に、図５で示した補償器２３は、オブザーバで構成され、下記式（４３）により、外乱抑圧制御量（電流）Ｕｄｉｓｔ（又は、Ｕｄ）を計算する。【０１５５】
【数４３】
2008142727.tif 000050
【０１５６】
図５の制御定数テーブル２８は、図１０で説明するように、式（４３）で使用するｓｉｎ（ωＴ／２）の値と、対応するＬ４、Ｌ５、ｃｏｓωＴ又はｓｉｎωＴの値のテーブルで構成される。ｃｏｓωＴ又はｓｉｎωＴの値は、式（４３）のＨの計算に使用する。【０１５７】
従って、サンプル毎に、ｓｉｎ（ωＴ／２）演算器２７で計算されたｓｉｎωＴ／２に対応するＬ４、Ｌ５、ｃｏｓωＴ，ｓｉｎωＴの値が、テーブル２８から読み出され、補償器２３にセットされる。補償器２３は、セットされた制御定数と、位置誤差から、式（４３）を計算し、外部抑圧制御量Ｕｄｉｓｔを計算し、加算ブロック２６に出力する。【０１５８】
このようにして、デイタル位置制御系において、外乱周波数に適応した位置制御を行う場合に、外乱周波数を精度良く推定し、外乱周波数に適応した位置制御が可能となる。【０１５９】
［位置制御系の第２の実施の形態］
図９は、本発明の第２の実施の形態の位置制御系のブロック図、図１０は、図９の制御定数テーブルの説明図である。図９は、コントローラ２１、補償器２３を、一体の現在オブザーバで構成した例である。又、ｓｉｎωＴ／２演算部２７を、第１の適応制御の計算方法を適用した例である。【０１６０】
図９の位置制御系は、外乱周波数を検出し、外乱を適応制御により、抑圧するため、下記式（４４）、（４５）、（４６）で示される外乱抑圧を行う現在オブザーバである。【０１６１】
【数４４】
2008142727.tif 000051
【０１６２】
【数４５】
2008142727.tif 000052
【０１６３】
【数４６】
2008142727.tif 000053
【０１６４】
尚、式（４６）において、ａ１１〜ａ２２は、式（９）、式（１８）で示したｓｉｎ，ｃｏｓの値である。【０１６５】
即ち、この実施の形態は、コントローラのモデルと外乱モデルを含むオブザーバによる位置制御系の例である。図９において、図５で示したものと同一のものは、同一の記号で示してある。【０１６６】
図９において、第１の演算ブロック２０は、ヘッド３（図１参照）が読み取った前述のサーボ情報を復調して得た観測位置ｙ［ｋ］から目標位置ｒを差し引き、実位置誤差ｅｒ［ｋ］を演算する。【０１６７】
第２の演算ブロック３２は、実位置誤差ｅｒ［ｋ］からオブザーバの推定位置ｘ［ｋ］を差し引き、推定位置誤差ｅ［ｋ］を演算する。状態推定ブロック３４は、この推定位置誤差ｅ［ｋ］と、内部変数である式（４４）の推定ゲインＬ（Ｌ１，Ｌ２、Ｌ３，Ｌ４，Ｌ５）を用いて、推定修正値（式（４４）の右辺第２項）を、演算する。【０１６８】
そして、加算ブロック３６は、この推定修正値と、遅延ブロック４６から状態量（式（４４）の右辺第１項）ｘ［ｋ］，ｖ［ｋ］、ｂ［ｋ］，ｚ１［ｋ］，ｚ２［ｋ］とを加算し、式（４４）のように、推定位置ｘ［ｋ］，推定速度ｖ［ｋ］、推定バイアス値ｂ［ｋ］，推定外乱抑圧値ｚ１［ｋ］，ｚ２［ｋ］を得る。尚、式（４４）では、推定位置誤差ｅ［ｋ］を、（ｙ［ｋ］−ｘ［ｋ］）で表示する。【０１６９】
この推定値は、第４の演算ブロック３８で、状態フィードバックゲイン（一Ｆａ＝Ｆ１，Ｆ２，Ｆ３，Ｆ４，Ｆ５）を乗算され、式（４５）のように、アクチュエータ１の駆動値ｕ［ｋ］を得る。一方、加算ブロック３６からの式（４４）の推定値は、第５の演算ブロック４２で、推定ゲインＡａ（式（４６）の右辺第１項の行列）を乗じられる。【０１７０】
第４の演算ブロック３８の駆動値ｕ［ｋ］は、第６の演算ブロック４０で、推定ゲインＢ（式（４６）の右辺第２項のｕ［ｋ］に乗じる値）を乗じられる。両乗算結果は、加算ブロック４４で、加算され、式（４６）の次のサンプルの推定状態量ｘ［ｋ＋１］，ｖ［ｋ＋１］、ｂ［ｋ＋１］，ｚ１［ｋ＋１］，ｚ２［ｋ＋１］を得る。【０１７１】
この次のサンプルの推定状態量は、前述のように、遅延ブロック４６に入力し、状態推定ブロック３４で、推定修正値で、修正される。そして、加算ブロック３６からの式（４４）の推定値は、第７の演算ブロック４８で、推定位置ｘ［ｋ］が取り出され、前述の第２の演算ブロック３２に入力する。【０１７２】
この現在オブザーバに、周波数推定器２７と、制御定数テーブル２８とを付与する。周波数推定器２７は、ｓｉｎ（ωＴ／２）の第１の適応則である式（３８）、（３９）を演算し、ｓｉｎ（ωＴ／２）を計算する。【０１７３】
即ち、周波数推定器２７では、状態時点変換ブロック２７−１は、式（２６）に従い、ｚ１［ｋ＋１］，ｚ２［ｋ＋１］を、ｚ１［ｋ＋０．５］，ｚ２［ｋ＋０．５］に変換する。推定誤差（Ｅ［ｋ］）演算器２７−２が、推定位置誤差ｅ［ｋ］、推定状態量ｚ１［ｋ＋０．５］，ｚ２［ｋ＋０．５］、推定ゲインＬ４，Ｌ５を受け、式（３８）を計算し、推定誤差Ｅ［ｋ］を出力する。【０１７４】
ゲイン乗算器２７−３は、推定誤差Ｅ［ｋ］にゲインＫを乗算し、式（３９）の右辺第２項を計算する。加算器２７−５は、遅延ブロック２７−４からの１サンプル前のｓｉｎ（ω［ｋ］Ｔ／２）と、ゲイン乗算器２７−３の出力とを加算し、式（３９）のｓｉｎ（ω［ｋ＋１］Ｔ／２）を出力する。このｓｉｎ（ω［ｋ＋１］Ｔ／２）は、遅延ブロック２７−４で１サンプル遅延され、加算器２７−５に入力する。【０１７５】
周波数推定器２７の出力ｓｉｎ（ω［ｋ＋１］Ｔ／２）は、制御定数テーブル２８に入力する。図１０に示すように、制御定数テーブル２８は、推定値ｓｉｎωＴ／２に対応する、推定ゲインＬ１、Ｌ２,Ｌ３,Ｌ４,Ｌ５、ｃｏｓωＴ，ｓｉｎωＴの値を格納する。ｓｉｎωＴ／２は、図８で説明したように、０〜Ｆｓ／２までを使用するため、図１０のテーブルは、３３個のｓｉｎωＴ／２の値を、０〜Ｆｓ／２の範囲に割り振る。【０１７６】
そして、制御定数テーブル２８は、推定値ｓｉｎωＴ／２に対応する推定ゲインＬ、ｃｏｓωＴ、ｓｉｎωＴを、引き出し、且つブロック３４，４２に設定する。これにより、ブロック３４の式（４４）の推定ゲインＬ１〜Ｌ５、ブロック４２の式（４６）の制御定数ａ１１〜ａ２２（＝ｃｏｓωＴ又はｓｉｎωＴ）が更新され、オブザーバの外乱抑圧特性を、外乱周波数に応じて、変更する。【０１７７】
このように、現在オブザーバで構成すると、設計が容易であり、且つ外乱抑圧機能を容易に設定できる。【０１７８】
図１１は、図９の周波数推定器２７の他の実施の形態のブロック図である。図１１において、図９で示したものと同一のものは、同一の記号で示してある。【０１７９】
前述のように、０〜Ｆｓ／２の周波数の範囲に対応するｓｉｎ（ωＴ／２）の値は、０〜１である。従って、推定したｓｉｎ（ωＴ／２）の値は、常に「０」から「１」の範囲に制限する必要がある。この実施の形態は、加算ブロック２７−５の出力を、「０」から「１」の範囲に制限するクリップブロック２７−６を設けたものである。【０１８０】
クリップブロック（飽和処理ブロック）２７−６は、入力されたｓｉｎ（ωＴ／２）が、「０」より小さい場合は、「０」に、「１」より大きい場合は、「１」に制限する。このため、状態変数からｓｉｎ（ωＴ／２）を推定しても、図８の周波数範囲の外乱周波数成分を得ることができる。【０１８１】
［位置制御系の第３の実施の形態］
図１２は、本発明の第３の実施の形態の位置制御系のブロック図である。図１２は、コントローラ２１、補償器２３を、一体の現在オブザーバで構成した例である。又、ｓｉｎ（ωＴ／２）演算部２７Ａに、第２の適応制御の計算方法を適用した例である。【０１８２】
図１２において、図９で示したものと同一のものは、同一の記号で示してある。図１２の位置制御系は、図９と同様に、外乱周波数を検出し、外乱を適応制御により、抑圧するため、前述の式（４４）、（４５）、（４６）で示されるバイアス補償を含む現在オブザーバである。従って、現在オブザーバ２１，２３，２６の説明は、省略する。【０１８３】
図１２に示すように、この現在オブザーバに、周波数推定器２７Ａと、制御定数テーブル２８とを付与する。周波数推定器２７Ａは、ｓｉｎ（ωＴ／２）の第２の適応則である式（４０）乃至式（４２）を演算し、ｓｉｎ（ωＴ／２）を計算する。【０１８４】
即ち、周波数推定器２７Ａでは、状態時点変換ブロック２７−１は、式（２６）に従い、ｚ１［ｋ＋１］，ｚ２［ｋ＋１］を、ｚ１［ｋ＋０．５］，ｚ２［ｋ＋０．５］に変換する。推定誤差（Ｇ［ｋ］）演算器２７−７が、推定状態量ｚ１［ｋ］,ｚ１［ｋ＋０．５］，ｚ２［ｋ］，ｚ２［ｋ＋０．５］を受け、式（４０）を計算し、第１の推定誤差Ｇ［ｋ］を出力する。【０１８５】
第１の演算ブロック２７−８は、第１の推定誤差Ｇ［ｋ］から、遅延ブロック２７−４からの１サンプル前のｓｉｎ（ω［ｋ］Ｔ／２）を引き、式（４１）の第２の推定誤差Ｅ［ｋ］を出力する。【０１８６】
ゲイン乗算器２７−８は、第２の推定誤差Ｅ［ｋ］にゲインＫを乗算し、式（４２）の右辺第２項を計算する。加算器２７−５は、遅延ブロック２７−４からの１サンプル前のｓｉｎ（ω［ｋ］Ｔ／２）と、ゲイン乗算器２７−３の出力とを加算し、式（４２）のｓｉｎ（ω［ｋ＋１］Ｔ／２）を出力する。このｓｉｎ（ω［ｋ＋１］Ｔ／２）は、遅延ブロック２７−４で１サンプル遅延され、演算ブロック２７−８、加算器２７−５に入力する。【０１８７】
この周波数推定器２７Ａの出力ｓｉｎ（ω［ｋ＋１］Ｔ／２）は、制御定数テーブル２８に入力する。制御定数テーブル２８は、図１０で示したように、推定値ｓｉｎ（ωＴ／２）に対応する、推定ゲインＬ１、Ｌ２,Ｌ３,Ｌ４,Ｌ５、ｓｉｎωＴ，ｃｏｓωＴの値を格納する。ｓｉｎ（ωＴ／２）は、図８で説明したように、０〜Ｆｓ／２までを使用するため、図１０のテーブルは、３３個のｓｉｎ（ωＴ／２）の値を、０〜Ｆｓ／２の範囲に割り振る。【０１８８】
そして、制御定数テーブル２８は、推定値ｓｉｎ（ωＴ／２）に対応する推定ゲインＬ、ｓｉｎωＴ，ｃｏｓωＴを、引き出し、且つブロック３４，４２に設定する。これにより、ブロック３４の式（４４）の推定ゲインＬ１〜Ｌ５、ブロック４２の式（４６）の制御定数ａ１１〜ａ２２（＝ｃｏｓωＴ又はｓｉｎωＴ）が更新され、オブザーバの外乱抑圧特性を、外乱周波数に応じて、変更する。【０１８９】
このように、現在オブザーバで構成すると、設計が容易であり、且つ外乱抑圧機能を容易に設定できる。又、第２の適応式を用いているので、入力の桁数に差異がなく、固定小数点方式のプロセッサでも、精度良く、計算できる。【０１９０】
図１３は、図１２の周波数推定器２７Ａの他の実施の形態のブロック図である。図１３において、図９、図１１、図１２で示したものと同一のものは、同一の記号で示してある。【０１９１】
前述のように、０〜Ｆｓ／２の周波数の範囲に対応するｓｉｎ（ωＴ／２）の値は、「０」〜「１」である。従って、推定したｓｉｎ（ωＴ／２）の値は、常に「０」から「１」の範囲に制限する必要がある。この実施の形態は、図１２の構成に対し、加算ブロック２４−５の出力を、「０」から「１」の範囲に制限するクリップブロック２４−６を設けたものである。【０１９２】
クリップブロック（飽和処理ブロック）２４−６は、入力されたｓｉｎ（ωＴ／２）が、「０」より小さい場合は、「０」に、「１」より大きい場合は、「１」に制限する。このため、状態変数からｓｉｎ（ωＴ／２）を推定しても、図８の周波数範囲の外乱周波数成分を得ることができる。【０１９３】
［位置制御系の第４の実施の形態］
図１４は、本発明の第４の実施の形態の位置制御系のブロック図である。図１４は、コントローラ２１と、補償器２３を、分離した現在オブザーバで構成した例である。又、ｓｉｎ（ωＴ／２）演算部２７に、第１の適応制御の計算方法を適用した例である。【０１９４】
図１４の位置制御系は、外乱周波数を検出し、外乱を適応制御により、抑圧するため、下記式（４７）、（４８）、（４９）、（５０）、（５１）で示される外乱抑圧を行う現在オブザーバである。【０１９５】
【数４７】
2008142727.tif 000054
【０１９６】
【数４８】
2008142727.tif 000055
【０１９７】
【数４９】
2008142727.tif 000056
【０１９８】
【数５０】
2008142727.tif 000057
【０１９９】
【数５１】
2008142727.tif 000058
【０２００】
尚、式（５１）において、ａ１１〜ａ２２は、式（９）、式（１８）で示したｓｉｎ，ｃｏｓの値である。【０２０１】
即ち、この実施の形態は、コントローラのモデルと外乱モデルを含むオブザーバによる位置制御系の例である。図１４において、図５、図９、図１２で示したものと同一のものは、同一の記号で示してある。【０２０２】
図１４において、第１の演算ブロック２０は、ヘッド３（図１参照）が読み取った前述のサーボ情報を復調して得た観測位置ｙ［ｋ］から目標位置ｒを差し引き、実位置誤差ｅｒ［ｋ］を演算する。【０２０３】
第２の演算ブロック３２は、実位置誤差ｅｒ［ｋ］からオブザーバの推定位置ｘ［ｋ］を差し引き、推定位置誤差ｅ［ｋ］を演算する。状態推定ブロック３４−１は、この推定位置誤差ｅ［ｋ］と、内部変数である式（４３）の推定ゲインＬ（Ｌ１，Ｌ２）を用いて、推定修正値（式（４７）の右辺第２項）を、演算する。【０２０４】
そして、加算ブロック３６−１は、この推定修正値と、遅延ブロック４６−１から状態量（式（４７）の右辺第１項）ｘ［ｋ］，ｖ［ｋ］とを加算し、式（４７）のように、推定位置ｘ［ｋ］，推定速度ｖ［ｋ］を得る。尚、式（４７）では、推定位置誤差ｅ［ｋ］を、（ｙ［ｋ］−ｘ［ｋ］）で表示する。【０２０５】
この推定値は、第４の演算ブロック３８−１で、状態フィードバックゲイン（一Ｆａ＝Ｆ１，Ｆ２）を乗算され、式（４８）のように、アクチュエータ１の駆動値ｕ［ｋ］を得る。一方、加算ブロック３６−１からの式（４７）の推定値は、第５の演算ブロック４２−１で、推定ゲインＡａ（式（５０）の右辺第１項の行列）を乗じられる。【０２０６】
第４の演算ブロック３８−１の駆動値ｕ［ｋ］は、第６の演算ブロック４０で、推定ゲインＢ（式（５０）の右辺第２項のｕ［ｋ］に乗じる値）を乗じられる。両乗算結果は、加算ブロック４４で、加算され、式（５０）の次のサンプルの推定状態量ｘ［ｋ＋１］，ｖ［ｋ＋１］を得る。【０２０７】
この次のサンプルの推定状態量は、前述のように、遅延ブロック４６−１に入力し、状態推定ブロック３４−１の推定修正値で、修正される。そして、加算ブロック３６−１からの式（４７）の推定値は、第５の演算ブロック４８で、推定位置ｘ［ｋ］が取り出され、前述の第２の演算ブロック３２に入力する。【０２０８】
同様に、外乱モデルにおいても、第２の状態推定ブロック３４−２は、この推定位置誤差ｅ［ｋ］と、内部変数である式（４７）の推定ゲインＬｄ１（Ｌ３、Ｌ４，Ｌ５）を用いて、推定修正値（式（４７）の右辺第２項）を、演算する。【０２０９】
そして、第２の加算ブロック３６−２は、この推定修正値と、第２の遅延ブロック４６−２から状態量（式（４７）の右辺第１項）ｂ［ｋ］，ｚ１［ｋ］，ｚ２［ｋ］とを加算し、式（４７）のように、推定バイアスｂ［ｋ］，推定外乱量ｚ１［ｋ］，ｚ２［ｋ］を得る。【０２１０】
この推定値は、第６の演算ブロック３８−２で、状態フィードバックゲイン（一Ｆｄ１＝Ｆ３、Ｆ４，Ｆ５）を乗算され、式（４９）のように、アクチュエータ１の外乱抑圧駆動値ｕｄ［ｋ］を得る。一方、第２の加算ブロック３６−２からの式（４７）の推定値は、第７の演算ブロック４２−２で、推定ゲインＡｄ１（式（５１）の右辺の行列）を乗じられる。【０２１１】
第８の演算ブロック５０は、第５の演算ブロック３８−２の駆動値ｕｄ［ｋ］と、第４の演算ブロック３８−１の駆動値ｕ［ｋ］とを加算し、アクチュエータの駆動値ｕ［ｋ］を出力する。【０２１２】
この例でも、この現在オブザーバに、周波数推定器２７と、制御定数テーブル２８とを付与する。周波数推定器２７は、図９、図１１で説明したものである。又、制御定数テーブル２８は、図１０で説明したものである。【０２１３】
そして、制御定数テーブル２８は、推定値ｓｉｎωＴに対応する推定ゲインＬ１、Ｌ２を、状態推定ブロック３４−１に、推定ゲインＬ３，Ｌ４，Ｌ５，ｓｉｎωＴ，ｃｏｓωＴを、ブロック３４−２，４２−２に設定する。これにより、ブロック３４−１，３４−２の式（４７）の推定ゲインＬ１〜Ｌ５、ブロック４２−２の式（５１）の制御定数ａ１１〜ａ２２（＝ｃｏｓωＴ又はｓｉｎωＴ）が更新され、オブザーバの外乱抑圧特性を、外乱周波数に応じて、変更する。【０２１４】
このように、外乱モデルを分離した現在オブザーバで構成すると、外乱モデルの設計が容易であり、且つ外乱抑圧機能を容易に設定できる。【０２１５】
［位置制御系の第５の実施の形態］
図１５は、本発明の第５の実施の形態の位置制御系のブロック図である。図１５は、コントローラ２１と、補償器２３を、分離した現在オブザーバで構成した例である。又、ｓｉｎ（ωＴ／２）演算部２７Ａに、第２の適応制御の計算方法を適用した例である。【０２１６】
図１５において、図１２、図１４で示したものと同一のものは、同一の記号で示してある。図１５の位置制御系は、外乱周波数を検出し、外乱を適応制御により、抑圧するため、前述の式（４７）、（４８）、（４９）、（５０）、（５１）で示される外乱抑圧を行う現在オブザーバである。【０２１７】
即ち、この実施の形態は、図１４と同様に、コントローラのモデルと外乱モデルを含むオブザーバによる位置制御系の例である。従って、オブザーバの説明は、省略する。【０２１８】
この例でも、この現在オブザーバに、周波数推定器２７Ａと、制御定数テーブル２８とを付与する。周波数推定器２７Ａは、図１２、図１３で説明したものである。又、制御定数テーブル２８は、図１０で説明したものである。【０２１９】
そして、制御定数テーブル２８は、推定値ｓｉｎ（ωＴ／２）に対応する推定ゲインＬ１、Ｌ２を、状態推定ブロック３４−１に、推定ゲインＬ３，Ｌ４，Ｌ５，ｓｉｎωＴ，ｃｏｓωＴを、ブロック３４−２，４２−２に設定する。これにより、ブロック３４−１，３４−２の式（４７）の推定ゲインＬ１〜Ｌ５、ブロック４２−２の式（５１）の制御定数ａ１１〜ａ２２（＝ｃｏｓωＴ又はｓｉｎωＴ）が更新され、オブザーバの外乱抑圧特性を、外乱周波数に応じて、変更する。【０２２０】
このように、外乱モデルを分離した現在オブザーバで構成すると、外乱モデルの設計が容易であり、且つ外乱抑圧機能を容易に設定できる。【０２２１】
［位置制御系の第６の実施の形態］
図１６は、本発明の位置制御系（サーボ制御系）の第６の実施の形態のブロック図、図１７は、図１６のタイムチャート図である。【０２２２】
一般に、サンプル周波数を高くすれば、位置決め精度の改善が期待できる。しかし、磁気ディスク装置においては、サーボ信号の記録領域を増やすことにつながり、その分、データを記録する面積が減少してしまう。【０２２３】
その改善策として、磁気ディスク装置においてマルチレート制御が用いられることがある。マルチレート制御は、フィルタを用いて実現することもできる。磁気ディスク装置においては、オブザーバ制御を用いてマルチレート制御を実現する方法が一般的である。【０２２４】
この実施の形態は、マルチレート制御を用いたオブザーバ制御のコントローラの構成を使用した実施の形態である。即ち、マルチレート制御は、入力のサンプル周波数以上の周波数の出力を実現するための方法である。この例では、２倍のマルチレート制御を構成している。すなわち、入力である観測位置ｙの検出サンプル周波数の２倍の周波数で電流Ｕを出力する。【０２２５】
図１６のマルチレート制御を用いたオブザーバ制御のコントローラ２１Ａは、式（４４）乃至（４６）から、以下の式（５２）で、構成される。【０２２６】
【数５２】
2008142727.tif 000059
【０２２７】
式（５２）で示されるように、基本的には、観測位置ｙの検出サンプルで、式（４４）乃至（４６）を２回実行し、ｕ［ｋ］，ｕ［ｋ＋０．５］を計算し、２倍の周波数で、電流ｕ［ｋ］，ｕ［ｋ＋０．５］を出力する。【０２２８】
即ち、図１７にも示すように、観測位置の検出サンプル（サーボゲート）で、観測位置ｙを取得した後、オブザーバの推定位置誤差、すなわち観測位置ｙ［ｋ］と、１サンプル前に計算した予測位置ｘ＿ｂａｒ［ｋ］との差分値により、状態変数ｘ＿ｈａｔ［ｋ］、ｖ＿ｈａｔ［ｋ］、ｂ＿ｈａｔ［ｋ］、ｚ１＿ｈａｔ［ｋ］，ｚ２＿ｈａｔ［ｋ］を修正する。【０２２９】
その上で、５つの状態変数ｘ＿ｈａｔ［ｋ］、ｖ＿ｈａｔ［ｋ］、ｂ＿ｈａｔ［ｋ］、ｚ１＿ｈａｔ［ｋ］，ｚ２＿ｈａｔ［ｋ］に、フィードバックゲインＦｘ，Ｆｖ，Ｆｂ，Ｆｚ１，Ｆｚ２をかけて、電流ｕ［ｋ］を出力する。【０２３０】
電流出力後に、０．５サンプル後のサンプルの状態変数の値ｘ＿ｈａｔ［ｋ＋０．５］，ｖ＿ｈａｔ［ｋ＋０．５］，ｂ＿ｈａｔ［ｋ＋０．５］、ｚ１＿ｈａｔ［ｋ＋０．５］，ｚ２＿ｈａｔ［ｋ＋０．５］を，ｕ［ｋ］により、推定する。【０２３１】
その上で、５の状態変数ｘ＿ｈａｔ［ｋ＋０．５］、ｖ＿ｈａｔ［ｋ＋０．５］、ｂ＿ｈａｔ［ｋ＋０．５］、ｚ１＿ｈａｔ［ｋ＋０．５］，ｚ２＿ｈａｔ［ｋ＋０．５］に、フィードバックゲインＦｘ，Ｆｖ，Ｆｂ、Ｆｚ１,Ｆｚ２をかけて、電流ｕ［ｋ＋０．５］を出力する。【０２３２】
電流出力後、次のサンプルの状態変数の値ｘ＿ｂａｒ［ｋ＋１］，ｖ＿ｂａｒ［ｋ＋１］、ｂ＿［ｋ＋１］、ｚ１＿ｈａｔ［ｋ＋１］，ｚ２＿ｈａｔ［ｋ＋１］を推定する。【０２３３】
具体的には、観測位置ｙの検出サンプル周波数の２倍の周波数のクロック（マルチレートサンプルクロック）３０を、コントローラ２１Ａ，２１Ｂに与え、Ｕｎ１（Ｕ［ｋ］），Ｕｎ２（Ｕ［ｋ＋０．５］）を，このクロックに同期して、Ｕｎとして、出力する。【０２３４】
一方、適応制御のブロック２７，２８は、オブザーバ２１Ａが、０．５サンプル後の状態変数ｚ１＿ｈａｔ［ｋ＋０．５］，ｚ２＿ｈａｔ［ｋ＋０．５］を計算しているので、これを式（３６）のｚ１［ｋ＋０．５］，ｚ２［ｋ＋０．５］に使用し、式（３６）を計算して、ｓｉｎ（ωＴ／２）の値を計算する。【０２３５】
そして、制御定数テーブル２８（図１０参照）から、ｓｉｎ（ωＴ／２）の値に対応するので、オブザーバ２１Ａに、推定ゲインＬ１，Ｌ２、Ｌ３，Ｌ４，Ｌ５，ｓｉｎωＴ，ｃｏｓωＴを、設定する。これにより、同様に、オブザーバの外乱抑圧特性を、外乱周波数に応じて、変更する。【０２３６】
［他の実施の形態］
以上の説明では、補償器２３に、オブザーバの構成で説明した。しかしながら、図５の補償器２３にフィルタを使用する場合にも、適用できる。以下に、その適応則を説明する。先ず、フィルタを使用する場合には、入力ゲインＬを補正するフィルタを使用し、元となる式は、式（８）、（９）の予測オブザーバの式と同じである。又、０．５サンプル戻すｚ１［ｋ＋０．５］、ｚ２［ｋ＋０．５］，ｚ１［ｋ＋１］，ｚ２［ｋ＋１］の関係式は、式（２７）と同一である。【０２３７】
式（２７）を、式（９）に代入すると、以下の式（５３）が得られる。【０２３８】
【数５３】
2008142727.tif 000060
【０２３９】
更に、式を整理するため、式（２８）の左辺、右辺ともに、下記式（５４）を掛ける。【０２４０】
【数５４】
2008142727.tif 000061
【０２４１】
この時、ｃｏｓ，ｓｉｎの行列間の掛け算では、下記式（５５）が成立する。【０２４２】
【数５５】
2008142727.tif 000062
【０２４３】
式（５５）を利用して、式（５３）の左辺、右辺に、式（５４）を掛けると、式（５６）が得られる。【０２４４】
【数５６】
2008142727.tif 000063
【０２４５】
即ち、ｓｉｎωＴの式（１７）と同様に、ｓｉｎ（ωＴ／２）の関係式が得られた。この式（５６）を、ｓｉｎωＴの導出と同様に展開する。即ち、ｓｉｎ（ωＴ／２），ｃｏｓ（ωＴ／２）を、それぞれ、「ｓ」、「ｃ」と簡略化すると、式（５６）は、式（５７）に変換される。【０２４６】
【数５７】
2008142727.tif 000064
【０２４７】
式（５７）を展開して、整理すると、式（５８）が得られる。【０２４８】
【数５８】
2008142727.tif 000065
【０２４９】
式（５８）において、ｓｉｎのみを消去するため、式（５８）の１行目に、（ｚ２［ｋ］＋Ｌ５・ｅ［ｋ］）を掛け、２行目に、（ｚ１［ｋ］＋Ｌ４・ｅ［ｋ］）を掛け、且つ１行目から２行目を引くと、式（５９）が得られる。【０２５０】
【数５９】
2008142727.tif 000066
【０２５１】
式（５９）のｓｉｎのみを右辺に残すように、変形すると、式（６０）が得られる。【０２５２】
【数６０】
2008142727.tif 000067
【０２５３】
式（６０）を整理すると、式（６１）が得られる。【０２５４】
【数６１】
2008142727.tif 000068
【０２５５】
式（６１）で、分母の式は、以下の式（６２）のように、ｅ［ｋ］への項の影響が充分小さい。【０２５６】
【数６２】
2008142727.tif 000069
【０２５７】
式（６２）の関係を利用して、式（６１）を変形すると、式（６３）が得られる。これを適応則に利用する。【０２５８】
【数６３】
2008142727.tif 000070
【０２５９】
以上のように、フィルタで構成した場合には、ｓｉｎ，ｃｏｓとも、分母のｚ１［ｋ］，ｚ２［ｋ］が、現在オブザーバの式（３６）とは、異なる。【０２６０】
又、適応式は、式（３８）乃至式（４２）の現在オブザーバの場合と同様である。【０２６１】
更に、ヘッドの位置制御装置で、説明したが、他の位置制御装置にも適用できる。同様に、ディスク装置を磁気ディスク装置で説明したが、他のデータ記憶装置にも適用できる。
【産業上の利用可能性】
【０２６２】
外乱の角周波数ωとし、前記サンプル周期をＴとした場合に、ｓｉｎ（ωＴ／２）を計算し、外乱適応制御するので、デジタル制御において、外乱の周波数の推定精度を向上でき、且つ外乱周波数の推定範囲を拡大でき、広い周波数範囲で、外乱適応制御を精度よく実現できる。このため、高精度の位置制御が要求される装置の、外乱による位置ずれを抑制することができる。
【図１】
2008142727.tif 000081
【図２】
2008142727.tif 000082
【図３】
2008142727.tif 000083
【図４】
2008142727.tif 000084
【図５】
2008142727.tif 000085
【図６】
2008142727.tif 000086
【図７】
2008142727.tif 000087
【図８】
2008142727.tif 000088
【図９】
2008142727.tif 000089
【図１０】
2008142727.tif 000090
【図１１】
2008142727.tif 000091
【図１２】
2008142727.tif 000092
【図１３】
2008142727.tif 000093
【図１４】
2008142727.tif 000094
【図１５】
2008142727.tif 000095
【図１６】
2008142727.tif 000096
【図１７】
2008142727.tif 000097
【図１８】
2008142727.tif 000098
【図１９】
2008142727.tif 000099
【国際調査報告】
2008142727.tif 000100
2008142727.tif 000101