(57)【要約】【課題】無駄な舵取り量の発生を抑えて省エネルギーの効果を高めることができる自動操舵装置及び方法を提供する。【解決手段】予め設定されたバッチ区間Ｂ１，Ｂ２毎に船体運動モデル及び制御ゲインを更新しつつ、バッチ区間よりも短い制御周期毎に針路偏差を解消し得る命令舵角を求めて航走体の針路を自動制御する自動操舵装置及び方法であって、バッチ区間Ｂ１からバッチ区間Ｂ２に切り替わる直前のバッチ終了期間Ｔ１内における針路偏差及び命令舵角を記憶し、バッチ区間Ｂ１からバッチ区間Ｂ２に切り替わったときに、船体運動モデル及び制御ゲインを更新するとともに、記憶した針路偏差及び命令舵角と更新した船体運動モデルとを用いて状態変数を更新し、更新した制御ゲインと更新した状態変数とを用いて新たなバッチ区間Ｂ２における命令舵角を求める。【選択図】図３
2011042284.tif 000002
【特許請求の範囲】
【請求項１】
予め設定されたバッチ区間毎に船体運動モデル及び制御ゲインを更新しつつ、前記バッチ区間よりも短い制御周期毎に針路偏差を解消し得る命令舵角を求めて航走体の針路を自動制御する自動操舵装置において、
前記バッチ区間が切り替わった場合に、前記船体運動モデル及び制御ゲインを更新するとともに、前記バッチ区間が切り替わる直前の所定期間内における針路偏差及び命令舵角と更新後の船体運動モデルとを用いて状態変数を更新し、更新された当該状態変数と更新された前記制御ゲインとを用いて命令舵角を求める演算を行う演算装置を備えることを特徴とする自動操舵装置。
【請求項２】
前記演算装置は、前記バッチ区間が切り替わってから最初の前記制御周期が経過する前に前記状態変数を更新することを特徴とする請求項１記載の自動操舵装置。
【請求項３】
前記所定期間は、少なくとも前記状態変数を更新し得る前記制御周期の数十周期分の期間であることを特徴とする請求項１又は請求項２記載の自動操舵装置。
【請求項４】
予め設定されたバッチ区間毎に船体運動モデル及び制御ゲインを更新しつつ、前記バッチ区間よりも短い制御周期毎に針路偏差を解消し得る命令舵角を求めて航走体の針路を自動制御する自動操舵方法において、
前記バッチ区間が切り替わる直前の所定期間内における針路偏差及び命令舵角を記憶する第１ステップと、
前記バッチ区間が切り替わった場合に、前記船体運動モデル及び制御ゲインを更新する第２ステップと、
前記第１ステップで記憶した前記針路偏差及び命令舵角と前記第２ステップで更新した船体運動モデルとを用いて状態変数を更新する第３ステップと、
前記第２ステップで更新した前記制御ゲインと前記第３ステップで更新した前記状態変数とを用いて命令舵角を求める第４ステップと
を有することを特徴とする自動操舵方法。
【請求項５】
前記第３ステップは、前記バッチ区間が切り替わってから最初の前記制御周期が経過する前に前記状態変数を更新するステップであることを特徴とする請求項４記載の自動操舵方法。
【請求項６】
前記第１ステップにおける前記所定期間は、少なくとも前記第３ステップで前記状態変数を更新し得る前記制御周期の数十周期分の期間であることを特徴とする請求項４又は請求項５記載の自動操舵方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、船舶等の航走体の針路を自動制御する自動操舵装置及び方法に関する。
【背景技術】
【０００２】
現在、大洋を運航する船舶等のように数日間以上に亘って運航されることのある船舶の殆どは、操縦者の睡眠時間の確保、或いは操縦者の負担軽減を図ることを目的として、オートパイロットといわれる自動操舵装置を備えている。この自動操舵装置は、一定周期で針路偏差（予め設定された針路と実際の針路との差）を求め、この針路偏差を解消し得る命令舵角を計算し、得られた命令舵角に基づいた制御を行うことにより船舶の自動操舵を実現している。また、自動操舵装置は、波や風等の外乱を考慮するために、バッチ区間と呼ばれる一定時間（例えば、４００秒）毎に外乱を同定した船体運動モデルを更新し、更には制御ゲインを更新しつつ上記の制御を行っている。
【０００３】
以下の特許文献１には、船舶の舵を最適制御することにより、目標方位角に対する方位角の針路偏差のみならず船体の横揺れをも制御するように船舶を自動操舵する技術が開示されている。具体的には、時々刻々と取得される船体の横揺れ角速度、方位角、舵角に関する各測定値並びに目標方位角及び運動特性行列を制御型多次元自己回帰モデルに適用することにより各時刻における外乱の推定値を算出し、複数の時刻からなる各バッチ区間の各推定値に対して局所定常過程に対するOzaki-Tong法を適用することにより各バッチ区間に対して最適な外乱自己回帰モデルを順次当てはめ、外乱自己回帰モデルに基づく外乱適応制御型多次元自己回帰モデルによって船体の横揺れ及び針路偏差を抑制する最適操舵量を算出して舵を制御している。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特許第３７５１２３９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
ところで、従来は、バッチ区間が切り替わると船体運動モデル及び制御ゲインが即座に更新されるが、命令舵角を求めるために用いられる状態変数の大部分は更新されない。このため、従来はバッチ区間が切り替わっても一定期間は更新される前の船体運動モデルに基づく状態変数を用いて命令舵角が求められているため、バッチ区間の切り替わりの前後において計算される命令舵角に不連続性が生ずるという問題があった。
【０００６】
図５は、従来の自動操舵装置における舵角及び針路偏差の一例を示す図である。図５を参照すると、バッチ区間が切り替わる時刻ｔ１００において、符号Ｒ１００を付した曲線で示す舵角が２度前後から急激に−８度前後まで変化し、その直後に急激に３度前後まで変化しているのが分かる。この舵角の変化は数秒程度の短時間の間に生じるものであるため、図示の通り符号Ｄ１００を付した曲線で示す針路偏差は殆ど変化しておらず、船舶の自動操舵に悪影響を及ぼすものではないと考えられる。しかしながら、図５に示す舵角の変化が生ずると無駄な舵取り量が発生することになるため、省エネルギー等の観点からはこのような舵角の変化は極力避ける必要がある。
【０００７】
本発明は上記事情に鑑みてなされたものであり、無駄な舵取り量の発生を抑えて省エネルギーの効果を高めることができる自動操舵装置及び方法を提供することを目的とする。
【課題を解決するための手段】
【０００８】
上記課題を解決するために、本発明の自動操舵装置は、予め設定されたバッチ区間（Ｂ１、Ｂ２）毎に船体運動モデル（Ａ、Ｂ）及び制御ゲイン（Ｇｎ）を更新しつつ、前記バッチ区間よりも短い制御周期毎に針路偏差（Ｘｎ）を解消し得る命令舵角（Ｙｎ）を求めて航走体の針路を自動制御する自動操舵装置（１）において、前記バッチ区間が切り替わった場合に、前記船体運動モデル及び制御ゲインを更新するとともに、前記バッチ区間が切り替わる直前の所定期間（Ｔ１）内における針路偏差及び命令舵角と更新後の船体運動モデルとを用いて状態変数を更新し、更新された当該状態変数と更新された前記制御ゲインとを用いて命令舵角を求める演算を行う演算装置（１１）を備えることを特徴としている。
また、本発明の自動操舵装置は、前記演算装置が、前記バッチ区間が切り替わってから最初の前記制御周期（ｔ２１）が経過する前に前記状態変数を更新することを特徴としている。
また、本発明の自動操舵装置は、前記所定期間が、少なくとも前記状態変数を更新し得る前記制御周期の数十周期分の期間であることを特徴としている。
上記課題を解決するために、本発明の自動操舵方法は、予め設定されたバッチ区間（Ｂ１、Ｂ２）毎に船体運動モデル（Ａ、Ｂ）及び制御ゲイン（Ｇｎ）を更新しつつ、前記バッチ区間よりも短い制御周期毎に針路偏差（Ｘｎ）を解消し得る命令舵角（Ｙｎ）を求めて航走体の針路を自動制御する自動操舵方法において、前記バッチ区間が切り替わる直前の所定期間（Ｔ１）内における針路偏差及び命令舵角を記憶する第１ステップ（Ｓ１６）と、前記バッチ区間が切り替わった場合に、前記船体運動モデル及び制御ゲインを更新する第２ステップ（Ｓ１７）と、前記第１ステップで記憶した前記針路偏差及び命令舵角と前記第２ステップで更新した船体運動モデルとを用いて状態変数を更新する第３ステップ（Ｓ１８）と、前記第２ステップで更新した前記制御ゲインと前記第３ステップで更新した前記状態変数とを用いて命令舵角を求める第４ステップ（Ｓ１５）とを有することを特徴としている。
また、本発明の自動操舵方法は、前記第３ステップが、前記バッチ区間が切り替わってから最初の前記制御周期（ｔ２１）が経過する前に前記状態変数を更新するステップであることを特徴としている。
また、本発明の自動操舵方法は、前記第１ステップにおける前記所定期間が、少なくとも前記第３ステップで前記状態変数を更新し得る前記制御周期の数十周期分の期間であることを特徴としている。
【発明の効果】
【０００９】
本発明によれば、バッチ区間が切り替わったときに、船体運動モデル及び制御ゲインが更新されるとともに、状態変数が新たな船体運動モデルに基づくものに更新されるため、バッチ区間の切り替わりの前後において計算される命令舵角の不連続性が解消される。これにより、無駄な舵取り量の発生が抑えられ、省エネルギーの効果を高めることができるという効果がある。
【図面の簡単な説明】
【００１０】
【図１】本発明の一実施形態による自動操舵装置の要部構成を示すブロック図である。
【図２】本発明の一実施形態による自動操舵方法を示すフローチャートである。
【図３】本発明の一実施形態による自動操舵方法を説明するための図である。
【図４】本実施形態による自動操舵装置における舵角及び針路偏差の一例を示す図である。
【図５】従来の自動操舵装置における舵角及び針路偏差の一例を示す図である。
【発明を実施するための形態】
【００１１】
以下、図面を参照して本発明の一実施形態による自動操舵装置及び方法について詳細に説明する。図１は、本発明の一実施形態による自動操舵装置の要部構成を示すブロック図である。図１に示す通り、本実施形態の自動操舵装置１は、演算装置１１、針路設定器１２、操作パネル１３、及び表示装置１４を備えており、方位検出器２１、船速検出器２２、及び舵角検出器２３の検出結果に基づいて、切替器３１及びアクチュエータ３２を介して舵３３の舵角を制御することにより、船舶の針路を自動制御する。
【００１２】
方位検出器２１は、回転型或いは振動式等の機械式ジャイロスコープ、光ファイバジャイロ或いはリングレーザジャイロ等の光学式ジャイロスコープ、その他のジャイロスコープであり、船舶の進行に関する方位角を検出して演算装置１１に出力する。船速検出器２２は、電磁式ログ、音響式ログ、回転翼式ログ、若しくは圧力式ログ等であり、船舶の進行速度（船速）を検出して演算装置１１に出力する。舵角検出器２３は、舵３３の舵角を検出して演算装置１１に出力する。
【００１３】
切替器３１は、自動操舵装置１による船舶の自動操舵とマニュアル操舵装置３４によるマニュアル操舵（航海士等の操縦者による手動操舵）とを切り替える。マニュアル操舵装置３４は、舵輪等によって船舶を手動操舵するためのものである。即ち、切替器３１が自動操舵を選択するように操作された場合は、自動操舵装置１から出力される舵角制御量（命令舵角）に基づいてアクチュエータ３２が駆動され、切替器３１が手動操舵を選択するように操作された場合には、マニュアル操舵装置３４から出力される命令舵角に基づいてアクチュエータ３２が駆動される。アクチュエータ３２は、切替器３１を介して出力される命令舵角に基づいて舵３３の舵角を設定する。舵３３は、周知のように、この舵角によって船舶の方位角を規定するものである。
【００１４】
演算装置１１は、針路設定器１２から出力される目標方位と、方位検出器２１、船速検出器２２、及び舵角検出器２３の各検出値とに基づいて、目標方位に対する進行方位の偏差（針路偏差）を求め、この針路偏差を解消し得る命令舵角を算出して切替器３１に出力する。ここで、演算装置１１は、上記の命令舵角を、予め設定された制御周期（例えば、１秒）毎に求める。
【００１５】
また、演算装置１１は、波や風等の外乱を考慮するために、上記の制御周期よりも長い予め設定されたバッチ区間（例えば、４００秒）毎に、外乱を同定した船体運動モデル及び制御ゲインを更新しつつ上記の命令舵角を求める。具体的に、演算装置１１は、以下の（１）式に基づいて命令舵角Ｙｎを求める。
【数１】
2011042284.tif 000003
【００１６】
ここで、上記（１）式中のＺｎは状態変数であり、Ｇｎは最適制御ゲインである。いま、状態変数Ｚｎを以下の（２）式の通り定義する。
【数２】
2011042284.tif 000004
【００１７】
すると、状態変数Ｚｎ及び針路偏差Ｘｎは、状態空間表現によって以下の（３）式の通り表すことができる。
【数３】
2011042284.tif 000005
【００１８】
但し、上記（３）式中のΦは状態遷移マトリクス、Γは制御ベクトル、Ｗｎは外乱、Ｈは観測ベクトルであり、それぞれ以下の（４）〜（７）式で表される。
【数４】
2011042284.tif 000006
【数５】
2011042284.tif 000007
【数６】
2011042284.tif 000008
【数７】
2011042284.tif 000009
【００１９】
尚、上記（４），（５）式、中の変数Ａ，Ｂは船体運動モデルであり、変数Ｍは船体運動モデルの次数である。この変数Ｍは可変であって最大値が「８」である。また、前述した（１）式中の最適制御ゲインＧｎは、以下の（８）式で表される。
【数８】
2011042284.tif 000010
【００２０】
演算装置１１は、前述した（１）式に基づいて制御周期毎に命令舵角Ｙｎを求める。また、バッチ区間が切り替わった場合には、上記（４），（５）式に示された状態遷移マトリクスΦ及び制御ベクトルΓの要素である船体運動モデルＡ，Ｂ及び次数Ｍと、最適制御ゲインＧｎを新たなものに更新した上で命令舵角Ｙｎを求める。
【００２１】
ここで、本実施形態では、バッチ区間の切り替わりの前後において求められる命令舵角Ｙｎの不連続性（図５参照）を解消するために、演算装置１１が、バッチ区間が切り替わる直前の所定期間（以下、バッチ終了期間という）内における針路偏差及び命令舵角を記憶しておき、バッチ区間が切り替わったときに、その針路偏差及び命令舵角と更新した船体運動モデルとを用いて状態変数Ｚｎを更新している。そして、演算装置１１が、更新された状態変数Ｚｎと、パッチ区間の切り替わりによって更新した最適制御ゲインＧｎとを用いて、前述した（１）式により命令舵角Ｙｎを求めている。尚、この処理の詳細については後述する。
【００２２】
針路設定器１２は、船舶の進行に応じた各所（緯度・経度）での目標方位を演算装置１１に出力する。操作パネル１３は、船舶が航行すべき航路である計画航路等の各種情報を演算装置１１に入力するためのものである。表示装置１４は、例えば液晶表示装置等を備えており、上記の計画航路やこれに対する船舶の航海状態（針路、方位、舵角、船速等）を表示する。
【００２３】
次に、上記構成における自動操舵装置１の動作、つまり本実施形態による自動操舵方法について詳細に説明する。図２は本発明の一実施形態による自動操舵方法を示すフローチャートであり、図３は本発明の一実施形態による自動操舵方法を説明するための図である。尚、図２に示す処理は、例えば切替器３１が自動操舵を選択するように操作されたときに開始され、切替器３１が手動操舵を選択するように操作されたときに終了する。
【００２４】
図２に示す処理が開始されると、まず自動操舵を行うために必要な各種の初期処理が演算装置１１で行われる（ステップＳ１１）。例えば、バッチ区間を計時するタイマ（以下、バッチタイマという）を初期化する処理、最初の船体運動モデルＡ，Ｂ及び最適制御ゲインＧｎを決定する処理、外乱Ｗｎを決定する処理等が行われる。尚、船体運動モデルＡ，Ｂの次数Ｍは、自動操舵装置１が搭載される船舶に応じて適宜決定される。
【００２５】
以上の初期処理が終了すると、バッチタイマによる計時が開始されて最初のバッチ区間が開始される。ここでは、図３に示す通り、時刻ｔ１において最初のバッチ区間Ｂ１が開始されたとする。最初のバッチ区間Ｂ１が開始されると、まず演算装置１１によってバッチタイマが参照され、バッチ終了期間であるか否かが判断される（ステップＳ１２）。
【００２６】
バッチ終了期間は、前述の通り、バッチ区間が切り替わる直前の所定期間をいい、図３に示す例では、バッチ区間Ｂ１についてのバッチ終了期間は、時刻ｔ１１〜時刻ｔ２までの符号Ｔ１が付された期間である。このバッチ終了期間は、後述するステップＳ１８で状態変数を更新し得る制御周期の数十周期分の期間に設定される。尚、本実施形態では、説明を簡単にするために、制御周期の２０周期分の期間がバッチ終了期間Ｔ１に設定されているとする。上記のステップＳ１２において、バッチ区間Ｂ１が開始された直後はバッチ終了期間Ｔ１には至っていないため、その判断結果は「ＮＯ」になる。
【００２７】
次に、演算装置１１によってバッチタイマが参照され、バッチ区間が終了したか否かが判断される（ステップＳ１３）。具体的には、バッチタイマによって１つのバッチ区間分の時間（例えば、４００秒）が計時されたか否かが判断される。バッチ区間が終了していない場合には、ステップＳ１３の判断結果は「ＮＯ」になる。すると、針路設定器１２からの目標方位と、方位検出器２１、船速検出器２２、及び舵角検出器２３の各検出値とに基づいて、目標方位に対する進行方位の偏差である針路偏差が演算装置１１で求められるとともに、前述した（３）式を用いて新たな状態変数Ｚｎが算出される（ステップＳ１４）。
【００２８】
新たな状態変数Ｚｎが算出されると、演算装置１１において、前述した（１）式を用いて針路偏差を解消し得る命令舵角Ｙｎが算出される（ステップＳ１５）。尚、演算装置１１によって算出された命令舵角Ｙｎは、切替器３１を介してアクチュエータ３２に出力され、アクチュエータ３２によって命令舵角Ｙｎに基づいた舵３３の舵角が設定される。これにより、船舶の航海状態が変更される。
【００２９】
以上の処理が終了すると、再び演算装置１１によってバッチタイマが参照されて、バッチ終了期間Ｔ１であるか否か判断され（ステップＳ１２）、次いでバッチ区間が終了したか否かが判断され（ステップＳ１３）、何れの判断結果も「ＮＯ」である場合に、針路偏差及び状態変数の算出（ステップＳ１４）及び命令舵角の計算（ステップＳ１５）が行われる。以後同様に、ステップＳ１２，Ｓ１３の判断結果が共に「ＮＯ」である場合には、ステップＳ１２〜Ｓ１５の処理ループが前述した制御周期（例えば、１秒）毎に繰り返される。
【００３０】
ここで、演算装置１１によってバッチタイマが参照され、バッチ終了期間Ｔ１であると判断された場合（ステップＳ１２の判断結果が「ＹＥＳ」の場合）には、算出された針路偏差及び命令舵角が演算装置１１に記憶される（ステップＳ１６：第１ステップ）。具体的には、図３に示す時刻ｔ１１（バッチ終了期間Ｔ１の開始時刻）において算出された針路偏差及び命令舵角を記憶する処理が行われる。
【００３１】
以上の記憶処理が終了すると、演算装置１１によってバッチタイマが参照されて、バッチ区間が終了したか否かが判断され（ステップＳ１３）、この判断結果が「ＮＯ」である場合には、針路偏差及び状態変数の算出（ステップＳ１４）及び命令舵角の計算（ステップＳ１５）が行われる。以後同様に、バッチ終了期間Ｔ１においては、ステップＳ１２の判断結果が「ＹＥＳ」になり、ステップＳ１３の判断結果が「ＮＯ」になるため、ステップＳ１２，Ｓ１６，Ｓ１３〜Ｓ１５の処理ループが前述した制御周期（例えば、１秒）毎に繰り返される。
【００３２】
図３中の記号「×」は制御周期（命令舵角が求めれる時点）を示している。バッチ終了期間Ｔ１では、図３中の記号「×」で示す各制御周期毎に針路偏差及び命令舵角が順次算出され、算出された命令舵角に命令舵角Ｙｎに基づいた舵３３の舵角が順次設定される。加えて、バッチ終了期間Ｔ１では、記号「×」で示す各制御周期毎に算出された針路偏差及び命令舵角が演算装置１１に順次記憶される。
【００３３】
ここで、演算装置１１によってバッチタイマが参照され、バッチ区間が終了したと判断された場合（ステップＳ１３の判断結果が「ＹＥＳ」の場合）には、演算装置１１によって、船体運動モデルＡ，Ｂ及び最適制御ゲインＧｎを更新し、バッチタイマを初期化する処理が行われる（ステップＳ１７：第２ステップ）。かかる処理により、前述した（３）式中の状態遷移マトリクスΦ及び制御ベクトルΓが更新される。
【００３４】
次に、ステップＳ１６で記憶した針路偏差及び命令舵角を用いて状態変数Ｚｎを更新する処理が行われる（ステップＳ１８：第３ステップ）。具体的には、バッチ区間Ｂ１のバッチ終了期間Ｔ１に記憶された制御周期の２０周期分の針路偏差及び命令舵角を、前述した（３）式に順次代入する（合計２０回代入する）ことにより状態変数Ｚｎを更新する処理が行われる。ここで、ステップＳ１７の処理によって、船体運動モデルＡ，Ｂが更新されているため、状態変数Ｚｎは、前のバッチ区間Ｂ１内のバッチ終了期間Ｔ１における針路偏差及び命令舵角と、更新された船体運動モデルＡ，Ｂとを用いて更新されることになる。
【００３５】
尚、以上のステップＳ１８の処理によって、状態変数Ｚｎはバッチ区間が切り替わってから最初の制御周期が経過する前までに更新される。つまり、図３に示す時刻ｔ２が経過してからバッチ区間Ｂ２の最初の制御周期である時刻ｔ２１が経過する前までに、状態変数Ｚｎの更新が行われる。以上の処理が終了すると、ステップＳ１７で更新された最適制御ゲインＧｎと、ステップＳ１８で更新された状態変数Ｚｎとを用いて針路偏差を解消し得る命令舵角Ｙｎが算出される（ステップＳ１５：第４ステップ）。そして、バッチ区間Ｂ２の最初の制御周期である時刻ｔ２１では、その命令舵角Ｙｎに基づいた舵３３の舵角が設定される。
【００３６】
図４は、本実施形態による自動操舵装置における舵角及び針路偏差の一例を示す図である。図４を参照すると、符号Ｒ１０を付した曲線で示す舵角は、全体的に細かな変動が生じているものの、バッチ区間が切り替わる時刻ｔ１０において、図５に示したような急激な変動は生じてないことが分かる。また、図４を参照すると、符号Ｄ１０を付した曲線で示す針路偏差も、全体的に細かな変動が生じているものの、殆ど変化していないことが分かる。尚、図４において、図５で生じていない舵角及び針路偏差の細かな変動が生じているのは、図４の結果を得るために用いた船体運動モデルと、図５の結果を得るために用いた船体運動モデルとが異なるためである。
【００３７】
以上の通り、本実施形態では、バッチ区間Ｂ１のバッチ終了期間Ｔ１における針路偏差及び命令舵角を記憶し、バッチ区間Ｂ１がバッチ区間Ｂ２に切り替わった場合に、船体モデルＡ，Ｂ及び最適制御ゲインＧｎを更新するとともに、記憶した針路偏差及び命令舵角と更新した船体モデルＡ，Ｂとを用いて状態変数Ｚｎを更新している。そして、更新した最適制御ゲインＧｎと更新した状態変数Ｚｎとを用いて命令舵角Ｙｎを求めている。
【００３８】
このように、本実施形態では、バッチ区間が切り替わったときに、船体運動モデル及び制御ゲインが更新されるとともに、状態変数Ｚｎが新たな船体運動モデルに基づくものに更新されるため、バッチ区間の切り替わりの前後において計算される命令舵角に生ずる不連続性を解消することができる。これにより、無駄な舵取り量の発生を抑えることができ、省エネルギーの効果を高めることができる。
【００３９】
以上、本発明の一実施形態による自動操舵装置及び方法について説明したが、本発明は上記実施形態に制限される訳ではなく、本発明の範囲内で自由に変更が可能である。例えば、上記実施形態ではバッチ終了期間が制御周期の２０周期分の期間に設定されている場合を例に挙げたが、必ずしも制御周期の２０周期分の期間に限られる訳ではく、図２中のステップＳ１８で状態変数を更新し得る期間に設定されていれば良い。
【符号の説明】
【００４０】
１            自動舵角装置
１１          演算装置
Ａ，Ｂ        船体運動モデル
Ｂ１，Ｂ２    バッチ区間
Ｇｎ           制御ゲイン
Ｔ１          バッチ終了期間
Ｘｎ           針路偏差
Ｙｎ           命令舵角
【図１】
2011042284.tif 000011
【図２】
2011042284.tif 000012
【図３】
2011042284.tif 000013
【図４】
2011042284.tif 000014
【図５】
2011042284.tif 000015