JP2004303085A - 船舶等の制御装置及び方法 - Google Patents

船舶等の制御装置及び方法 Download PDF

Info

Publication number
JP2004303085A
JP2004303085A JP2003097319A JP2003097319A JP2004303085A JP 2004303085 A JP2004303085 A JP 2004303085A JP 2003097319 A JP2003097319 A JP 2003097319A JP 2003097319 A JP2003097319 A JP 2003097319A JP 2004303085 A JP2004303085 A JP 2004303085A
Authority
JP
Japan
Prior art keywords
control
behavior
sequentially
time range
control amount
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2003097319A
Other languages
English (en)
Other versions
JP3683890B2 (ja
Inventor
Retsu Yamakawa
烈 山川
Hitoshi Maeno
仁 前野
Eiji Uchino
英治 内野
Hirohiko Morita
博彦 森田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
FUZZY LOGIC SYSTEMS
Original Assignee
FUZZY LOGIC SYSTEMS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by FUZZY LOGIC SYSTEMS filed Critical FUZZY LOGIC SYSTEMS
Priority to JP2003097319A priority Critical patent/JP3683890B2/ja
Priority to GB0406640A priority patent/GB2400189B/en
Priority to US10/809,399 priority patent/US7403841B2/en
Publication of JP2004303085A publication Critical patent/JP2004303085A/ja
Application granted granted Critical
Publication of JP3683890B2 publication Critical patent/JP3683890B2/ja
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B63SHIPS OR OTHER WATERBORNE VESSELS; RELATED EQUIPMENT
    • B63HMARINE PROPULSION OR STEERING
    • B63H25/00Steering; Slowing-down otherwise than by use of propulsive elements; Dynamic anchoring, i.e. positioning vessels by means of main or auxiliary propulsive elements
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/02Control of position or course in two dimensions
    • G05D1/0206Control of position or course in two dimensions specially adapted to water vehicles
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B11/00Automatic controllers
    • G05B11/01Automatic controllers electric
    • G05B11/36Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/02Control of position or course in two dimensions
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Theoretical Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Mathematical Physics (AREA)
  • Databases & Information Systems (AREA)
  • Data Mining & Analysis (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Mechanical Engineering (AREA)
  • Ocean & Marine Engineering (AREA)
  • Feedback Control In General (AREA)

Abstract

【課題】船舶等の制御対象の挙動を判別して、その判別結果を用いて効率的な制御を行うことができる制御装置及び方法を提供すること。
【解決手段】方位センサにより取得される方位データに基づいて船舶により行われる所定種類の挙動の時間範囲を順次判断する挙動検出部52と、挙動検出部52により順次判断される時間範囲において方位センサにより取得される方位データに基づいて、各時間範囲において船舶により行われた前記所定種類の挙動が当てはまる挙動パターンを複数の挙動パターンの中から順次判断するパターンマッチング部58と、パターンマッチング部58により順次判断される、前記所定種類の挙動が当てはまる挙動パターンに対応づけて記憶手段に記憶される制御パラメータを順次読み出し、該制御パラメータに基づいて船舶の制御を行わせる制御パラメータ設定部62と、を含むことを特徴とする。
【選択図】 図4

Description

【0001】
【発明の属する技術分野】
本発明は制御装置及び方法に関し、特に複数の挙動パターンのいずれかに当てはまる挙動を繰り返す制御対象を好適に制御する技術に関する。
【0002】
【背景技術】
従来一般のリアルタイム制御技術では、ある時刻の制御量を、例えば該制御量と目標値との偏差によって評価し、それにより次の時刻で用いる操作量を決定している。しかしながら、周波数帯域が制御量出力信号と重なっていたり、周波数帯域が不定(未知又は時変)であったりして、従来一般のフィルタリングでは除去困難な外乱が制御量に重畳されている場合には、制御量を適切に評価することが困難であり、それ故、制御対象を適切に制御することが困難であるという問題がある。
【0003】
本発明は上記背景のもとでなされたものであって、その目的は、複数の挙動パターンのいずれかに当てはまる挙動を繰り返す制御対象を、従来一般のフィルタリングでは除去困難な外乱の影響を可及的に排して、好適に制御することができる制御装置及び方法を提供することにある。
【0004】
【課題を解決するための手段】
上記課題を解決するために、本発明に係る制御装置は、制御対象に関する制御量を制御する制御装置であって、制御量を順次取得する制御量取得手段と、前記制御量取得手段により取得される制御量に基づいて、前記制御対象により行われる所定種類の挙動の時間範囲を順次判断する時間範囲判断手段と、前記時間範囲判断手段により順次判断される時間範囲において前記制御量取得手段により取得される制御量に基づいて、各時間範囲において前記制御対象により行われた前記所定種類の挙動が当てはまる挙動パターンを複数の挙動パターンの中から順次判断する挙動パターン判断手段と、前記複数の挙動パターンの各々に対応づけて制御パラメータを記憶する制御パラメータ記憶手段と、前記挙動パターン判断手段により順次判断される、前記所定種類の挙動が当てはまる挙動パターンに対応づけて前記制御パラメータ記憶手段に記憶される制御パラメータを順次読み出す制御パラメータ読み出し手段と、前記制御パラメータ読み出し手段により順次読み出される制御パラメータに基づいて、前記制御対象を制御する制御手段と、を含むことを特徴とする。
【0005】
また、本発明に係る制御方法は、制御対象に関する制御量を制御する制御方法であって、制御量を順次取得する制御量取得ステップと、前記制御量取得ステップで取得される制御量に基づいて、前記制御対象により行われる所定種類の挙動の時間範囲を順次判断する時間範囲判断ステップと、前記時間範囲判断ステップで順次判断される時間範囲において前記制御量取得ステップで取得される制御量に基づいて、各時間範囲において前記制御対象により行われた前記所定種類の挙動が当てはまる挙動パターンを複数の挙動パターンの中から順次判断する挙動パターン判断ステップと、前記複数の挙動パターンの各々に対応づけて制御パラメータを記憶する制御パラメータ記憶手段から、前記挙動パターン判断ステップで順次判断される、前記所定種類の挙動が当てはまる挙動パターンに対応づけて記憶される制御パラメータを順次読み出す制御パラメータ読み出しステップと、前記制御パラメータ読み出しステップで順次読み出される制御パラメータに基づいて、前記制御対象を制御する制御ステップと、を含むことを特徴とする。
【0006】
本発明によれば、制御量が順次取得され、該制御量に基づいて制御対象により行われる所定種類の挙動の時間範囲が順次判断される。制御対象は船舶、自動車、二輪車、航空機、ロボット、建造物等、何であってもよく、制御量は位置、姿勢、温度、濃度等、何であってもよい。また、所定種類の挙動の時間範囲は、例えば制御量が極大値をとるタイミングを始期とし、次に極大値をとるタイミングを終期とする時間範囲である。その他、制御量が極小値をとるタイミングから次に極小値をとるタイミングまでの時間範囲であってもよいし、制御量が極小値をとるタイミングから次に極大値をとるタイミングまでの時間範囲、制御量が変曲するタイミングから次に変曲するタイミングまでの時間範囲、制御量が目標値に一致するタイミングから次に一致するタイミングまでの時間範囲等であってもよい。いずれにしても、所定種類の挙動の時間範囲は、順次取得される制御量に基づいて順次判断されるようになっている。
【0007】
そして、各時間範囲において取得される制御量に基づき、その時間範囲において行われる挙動が当てはまる挙動パターンが複数の挙動パターンの中から順次判断される。そして、判断される挙動パターンに対応づけて記憶されている制御パラメータに基づいて制御対象が制御されるようになっている。こうして、本発明によれば、制御対象の挙動を判別して、その判別結果を用いて効率的な制御を行うことができるようになる。
【0008】
また、本発明の一態様では、前記時間範囲判断手段により順次判断される時間範囲において前記制御量取得手段により取得される制御量に基づいて、各時間範囲において前記制御対象により行われた前記所定種類の挙動に対する評価値を順次算出する挙動評価値算出手段と、前記挙動評価値算出手段により順次算出される評価値に基づいて、前記制御パラメータ記憶手段に記憶される制御パラメータを更新する制御パラメータ更新手段と、をさらに含む。こうすれば、適応制御を行うことができるようになる。なお、評価値は、例えば挙動の大きさを表すもの等であり、挙動の時間範囲における制御量の変化量等を採用することができる。また、この態様では、前記制御パラメータ更新手段は、前記挙動評価値算出手段により順次算出される評価値に対応する前記所定種類の挙動が行われた時間範囲において用いられた制御パラメータを更新するようにしてもよい。
【0009】
また、本発明の一態様では、前記時間範囲判断手段は、前記制御量取得手段により取得される制御量に基づいて、前記制御量取得手段により取得される制御量が極値をとるタイミングを前記所定種類の挙動の時間範囲の始期及び終期と判断する。こうすれば、比較的簡単に所定種類の挙動の時間範囲を判断することができる。
【0010】
また、本発明の一態様では、前記制御対象は船舶であり、前記制御量はその針路であり、前記制御手段は前記船舶に備えられた舵機を制御する。こうすれば、船舶の進行方向を好適に制御することができる。
また、本発明に係る制御装置は、制御対象に関する制御量を制御する制御装置であって、制御量を順次取得する制御量取得手段と、前記制御量取得手段により取得される制御量に基づいて、前記制御対象により行われる所定種類の挙動の時間範囲を順次判断する時間範囲判断手段と、前記時間範囲判断手段により順次判断される時間範囲において前記制御量取得手段により取得される制御量に基づいて、各時間範囲において前記制御対象により行われた前記所定種類の挙動が当てはまる挙動パターンを複数の挙動パターンの中から順次判断する挙動パターン判断手段と、前記挙動パターン判断手段により順次判断される、前記所定種類の挙動が当てはまる挙動パターンに対応する制御パラメータを順次取得する制御パラメータ取得手段と、前記制御パラメータ取得手段により順次取得される制御パラメータに基づいて、前記制御対象を制御する制御手段と、を含むことを特徴とする。
また、本発明に係る制御方法は、制御対象に関する制御量を制御する制御方法であって、制御量を順次取得する制御量取得ステップと、前記制御量取得ステップで取得される制御量に基づいて、前記制御対象により行われる所定種類の挙動の時間範囲を順次判断する時間範囲判断ステップと、前記時間範囲判断ステップで順次判断される時間範囲において前記制御量取得ステップで取得される制御量に基づいて、各時間範囲において前記制御対象により行われた前記所定種類の挙動が当てはまる挙動パターンを複数の挙動パターンの中から順次判断する挙動パターン判断ステップと、前記挙動パターン判断ステップで順次判断される、前記所定種類の挙動が当てはまる挙動パターンに対応する制御パラメータを順次取得する制御パラメータ取得ステップと、前記制御パラメータ取得ステップで順次取得される制御パラメータに基づいて、前記制御対象を制御する制御ステップと、を含むことを特徴とする。
本発明によれば、制御対象の挙動を判別して、その判別結果を用いて効率的な制御を行うことができるようになる。なお、制御パラメータ取得手段及び制御パラメータ取得ステップでは、ファジィ推論、ニューラルネットワーク、ウェーブレットネットワーク、RBF(Radial Basis Function)ネットワーク等を用いて制御パラメータを取得してもよい。
【0011】
【発明の実施の形態】
以下、本発明の実施形態について図面に基づき詳細に説明する。
【0012】
図1は、本発明の一実施形態に係る自動操舵制御装置の構成を示す図である。同図に示すように、この自動操舵制御装置10は、船舶(制御対象)に備えられるものであり、針路設定部12、方位センサ14、舵角センサ15、舵機16、加算器18,24、制御パラメータ決定部20、操舵量決定部22及びデッドバンド処理部26を含んで構成されており、舵機16を自動操作して、船舶の船首方位を制御するようになっている。
【0013】
針路設定部12は、本船舶の進むべき方位θ(目標値)を出力する。この方位θは、方位設定つまみにより手動設定されてもよいし、例えば衛星測位システムを含んで構成される公知の船舶機器により演算生成されてもよい。この方位θは、真北を基準とし、そこからのずれ角を右(東)回りに0°〜360°の範囲で表したものである。針路設定部12の出力は負符号を与えられて加算器18に供給されている。
【0014】
方位センサ14は、本船舶の舳先が現在向いている方位θ(制御量)を所定時間毎に出力する。この船首方位θも、真北を基準とし、そこからのずれ角を右(東)回りに0°〜360°の範囲で表したものであり、加算器18及び制御パラメータ決定部20に供給されている。なお、船首方位θは、例えばローパスフィルタを施してから出力してもよい(加算器18の出力に対してローパスフィルタを施してもよい)。
【0015】
舵機16は、本船舶に備えられた舵を例えば油圧ポンプやシリンダ等により駆動する舵機駆動部と、実舵角を指令舵角に一致させる舵機制御部と、を含んで構成された公知の船舶機器である。舵機16には、舵角センサ15が接続されており、該舵角センサ15から現在の舵角、すなわち実舵角δrが出力されるようになっている。実舵角δrは加算器24及びデッドバンド処理部26に供給されている。また、舵機16は、デッドバンド処理部26から指令舵角δr+δDを入力することにより、実舵角をδrからδr+δDに変化させるようになっている。
【0016】
加算器18は、方位センサ14から出力される船首方位θと針路設定部12から出力される基準方位θとの偏差を生成し、それを操舵量決定部22に供給している。なお、加算器18から出力される偏差θ−θは、±180°の範囲に正規化される。
【0017】
操舵量決定部22はPID制御(比例制御+積分制御+微分制御)を行う公知の制御手段であって、加算器18及び制御パラメータ決定部20の出力が供給されており、制御パラメータ決定部20から供給される制御パラメータ(比例係数KP、積分係数KI及び微分係数KD)に基づいて、加算器18から供給される方位偏差θ−θ(±180°の範囲に正規化したもの)から操舵量δPIDを算出する。この操舵量決定部22はハードウェアのみによって構成されてもよいし、コンピュータと該コンピュータが実行するプログラムによって構成されてもよく、具体的には、図2に示すように、操舵量決定部22は、積分器30、微分器36、係数乗算器32,34,38、加算器40を含んで構成されている。そして、加算器18から出力される偏差θ−θに対して、積分器30により積分するとともに、積分結果に対して係数乗算器32により積分係数KIを乗算する。また、並行して係数乗算器34により比例係数KPを乗算する。さらに、並行して微分器36により微分するとともに、微分結果に対して係数乗算器38により微分係数KDを乗算する。加算器40には、それらの演算結果が入力されており、総和を操舵量δPIDとして出力するようになっている。
【0018】
この操舵量δPIDは加算器24に供給されている。加算器24には舵角センサ15から出力される実舵角δrも負符号を与えられて供給されており、操舵量δPIDと実舵角δrとの差分が演算されるようになっている。演算結果はデッドバンド処理部26に供給されている。
【0019】
デッドバンド処理部26は、入力値である操舵量δPIDと実舵角δrとの差分の絶対値が所定値DB未満である場合には内部値δDを零とし、それ以上の場合には入力値をそのまま内部値δDとする処理を行う。デッドバンド処理部26には舵角センサ15から実舵角δrも入力されている。そして、デッドバンド処理部26は、この実舵角δrと内部値δDとの和を出力し、舵機16に供給している。こうして、操舵量決定部22の出力に対する舵角δr+δDの不感帯を設けている。このデッドバンド処理部26も、ハードウェアのみによって構成されてもよいし、コンピュータと該コンピュータが実行するプログラムによって構成されてもよい。デッドバンド処理部26における入力値と内部値との関係は、図3に示される通りである。
【0020】
制御パラメータ決定部20は、操舵量決定部22において操舵量δPIDを決定するために用いる制御パラメータ(比例係数KP、積分係数KI及び微分係数KD)を決定する処理を行うものであり、方位センサ14から出力される船首方位θと図示しない衛星測位装置又は船速センサから供給されている現在船速Vとが入力されている。制御パラメータ決定部20は、これらのデータに基づいて制御パラメータを決定する。この制御パラメータ決定部20も、ハードウェアのみによって構成されてもよいし、コンピュータと該コンピュータが実行するプログラムによって構成されてもよい。
【0021】
具体的には、制御パラメータ決定部20は、図4に示されるように、挙動検出部52、挙動特徴量算出部56、挙動評価値算出部54、パターンマッチング部58,制御パラメータ更新部60及び制御パラメータ設定部62を含んで構成されている。挙動検出部52には、加算器18から出力される偏差θ−θが入力されており、船舶の所定挙動(偏差θ−θが極大値をとるタイミングから再び極大値をとるタイミングまでの挙動(水平方向の船首揺,ヨーイング))の時間範囲を順次判断する。
【0022】
例えば、挙動検出部52に順次入力される偏差θ−θに基づいて、最新の偏差θ−θから直前の偏差θ−θの差分を順次計算し、その値が正から負に変化するタイミングを制御量である船首方位θが極大値をとるタイミングであると判断する。そして、このタイミングを直前の挙動の終了タイミング、且つ次の挙動の開始タイミングと判断する。すなわち、図5に示されるように、加算器18から出力される偏差θ−θは一般には増減を繰り返しており、図中波線で示される、偏差θ−θが極大値をとるタイミング、すなわち船首方位θが極大値をとるタイミングを検知し、それをある挙動の開始タイミング、且つ次の挙動の終了タイミングとして挙動特徴量算出部56に供給するようになっている。なお、挙動検出部52は、その他、船首方位θが極小値をとるタイミングや、船首方位θの二階微分が符号反転するタイミング等に基づいて、挙動の時間範囲の開始タイミングや終了タイミングを決定してもよい。
【0023】
挙動特徴量算出部56には、挙動検出部52から各挙動の開始タイミング及び終了タイミングが順次供給されるとともに、加算器18から偏差θ−θが順次供給されている。挙動特徴量算出部56は記憶手段を備えており、少なくとも1挙動分の偏差θ−θが記憶されるようになっている。そして、各挙動の時間範囲において加算器18から供給された偏差θ−θに基づいて、該挙動の特徴量を算出するようになっている。なお、挙動特徴量算出部56に記憶手段を設けることなく、逐次各種挙動特徴量を算出するようにしてもよい。
【0024】
ここでは、挙動特徴量算出部56は、各挙動の特徴量として、1つの挙動中に取得される偏差θ−θの平均値θc、1つの挙動中に取得される偏差θ−θの一階微分値ωの平均値ωc、1つの挙動中に取得される偏差θ−θの最大値と最小値との差Δθ、1つの挙動中に取得される偏差θ−θの一階微分ωの最大値と最小値との差Δωを算出するようにしている。図6は、ある挙動に対応する偏差θ−θを、横軸が偏差θ−θであり、縦軸がその一階微分であるωである位相面に表したものであり、挙動特徴量算出部56は、同図に示されるθc、ωc、Δθ及びΔωを特徴量として算出するようにしている。
【0025】
こうして算出される各特徴量は挙動評価値算出部54とパターンマッチング部58に供給されている。具体的には、平均値θc及びωcはパターンマッチング部58に供給され、平均値θc及びωcと偏差θ−θ及び一階微分ωの各最大値と最小値の差Δθ及びΔωは挙動評価値算出部54に供給されている。パターンマッチング部58には、図示しない衛星測位装置や船速センサから本船舶の速度、すなわち船速Vも供給されている。なお、ここでは挙動特徴量として上記のものを採用したが、もちろん他の量を挙動特徴量としてもよい。
【0026】
パターンマッチング部58は、挙動検出部52により順次判断される時間範囲において方位センサ14により取得される船首方位θに基づいて、各時間範囲において本船舶により行われた挙動が当てはまる挙動パターンを、複数の挙動パターンの中から順次判断する。このパターンマッチング部58によるパターンマッチング処理は、ここでは後述するメンバーシップ関数群及びファジィルールテーブル群を用いたファジィ推論により行うこととするが、その他、ニューラルネットワーク、RBF(Radial Basis Function)ネットワーク、ウェーブレットネットワーク等により行うこととしてもよい。また、より単純に、各挙動パターンに対応する各挙動特徴量の数値範囲を予め決定しておき、実際の挙動に対して算出される挙動特徴量がどの数値範囲に属するかにより、その挙動に対応する挙動パターンを判断してもよい。
【0027】
また、パターンマッチング部58は、複数の挙動パターンの各々に対応づけて制御パラメータの組(比例係数KP、積分係数KI及び微分係数KD)を記憶する、制御パラメータ記憶部59を含んで構成されている。より具体的には、制御パラメータ記憶部59は、各挙動パターン(前件部)に対応づけて制御パラメータの組(後件部シングルトン)を複数記憶しており、各組は船速Vに対応づけられている。すなわち、制御パターン記憶部59は、挙動パターンと船速Vの組合せに対応づけて制御パラメータの組を記憶している。この制御パラメータ記憶部59は、例えばRAMやROM、或いはハードディスク記憶装置等の公知の記憶手段を含んで構成されてよい。そして、パターンマッチング部58は、こうして順次判断される、各挙動が当てはまる挙動パターンに対応づけて、制御パラメータ記憶部59に記憶される制御パラメータの組を順次読み出し、それを制御パラメータ設定部62に供給している。なお、各挙動が当てはまる挙動パターンが複数ある場合には、それらに対応する制御パラメータの重心等を制御パラメータ設定部62に供給すればよい。
【0028】
ここで、パターンマッチング部58による処理をファジィ推論により行う例について説明する。図7は、パターンマッチング部58に記憶されている、メンバーシップ関数群の一例を示す図である。パターンマッチング部58に入力される特徴量θc及びωcは、同図(a)及び(b)に示されるメンバーシップ関数に照査され、それぞれの特徴量θc及びωcと各ファジィ集合との適合度が算出される。同様に、パターンマッチング部58に入力される船速Vも、同図(c)に示されるメンバーシップ関数に照査され、それと各ファジィ集合との適合度が算出される。その後、こうして算出される適合度に基づいて1又は複数のIF−THENルールが選出される。
【0029】
図8は、制御パラメータ記憶部59に記憶されている、ファジィルールテーブル群の一例を示す図である。同図に示すように、制御パラメータ記憶部59では、特徴量θcと特徴量ωcと船速Vとを前件部とし、制御パラメータの組Pmn(シングルトン)を後件部とする3入力1出力ファジィルールテーブルを複数記憶している(3入力1出力ファジィルールテーブルが制御パラメータ(KP,KI,KD)毎に記憶されている)。なお、ここでは船速Vが負の場合、手動による操舵が行われる。同テーブル群から、特徴量θc及びωc並びに船速Vに対応する制御パラメータの組Pmnを取得できるようになっている。こうして取得される制御パラメータの組は、制御パラメータ設定部62によって、操舵量決定部22に含まれる係数乗算部32,34,38にそれぞれ設定されるようになっている(後述するように、適応制御とする場合には、制御パラメータにそれぞれ所定値を加算したものを操舵量決定部22に供給することになる)。
【0030】
挙動評価値算出部54は、上述のように挙動特徴量算出部56から各挙動の特徴量を取得する。そして、それら特徴量に基づいて各挙動に対する評価値fを算出する。この評価値fは、例えば次式に従って算出され、制御パラメータ更新部60に供給される。
【0031】
【数1】
f=θc+ωc+Δθ+Δω
【0032】
ここで、特徴量θc及びωcは各挙動の原点からの全体的なずれを表しており、これらの値が小さいほど望ましい挙動であると判断できる。また、特徴量Δθ及びΔωは各挙動の大きさを表しており、これらの値が小さいほど望ましい挙動であると判断できる。そこで、本自動操舵制御装置10では、これらの値の増加関数である上記評価値fを最小化するよう、制御パラメータ記憶部59に記憶されている制御パラメータを更新するようにしている。
【0033】
すなわち、制御パラメータ設定部62は、制御パラメータ記憶部59から読み出された制御パラメータ(比例係数KP、積分係数KI及び微分係数KD)のうち、最も実際の挙動に対する適合度が高かった挙動パターンに対応するものを更新量ΔKP,ΔKI及びΔKDだけそれぞれ増加させ、制御パラメータ記憶部59の記憶内容を更新する。そして、パターンマッチング部58により、再度挙動パターンに対応する制御パラメータの組を読み出し、その重心等を制御パラメータ設定部62に供給する。そして、この制御パラメータは、制御パラメータ設定部62により操舵量決定部22に供給される。なお、更新量ΔKP、ΔKI及びΔKDは、固定値であってもよいが、船首方位θやそれに基づく挙動特徴量等に基づいて算出されるようにすれば好適である。例えば、連続する挙動の各評価値fn+1及びfnの差分に比例する量としてもよい。なお、この場合は、ΔKP、ΔKI及びΔKDが負になることもあり得る。更新量ΔKP等が負値を取りうる場合には、以下の説明において、比例係数KP等の増加は、比例係数KP等に更新量ΔKP等を加算することを意味し、比例係数KP等の減少は、比例係数KP等から更新量ΔKP等を減算することを意味する。
【0034】
制御パラメータ更新部60は、上記のようにして供給される制御パラメータに基づいて操舵量決定部22により制御を行った結果、次の挙動の評価が向上していれば(評価値が減少していれば)、ΔKP,ΔKI及びΔKDだけ各制御パラメータを増加させたことが制御の改善に貢献したと判断し、前回の挙動に対する適合度が最も高かった挙動パターンに対応する各制御パラメータを、さらにΔKP,ΔKI及びΔKDだけ増加させる。また、次の挙動の評価が低下していれば(評価値が増加していれば)、ΔKP,ΔKI及びΔKDだけ各制御パラメータを増加させたことが制御の悪化に繋がったと判断し、2ΔKP,2ΔKI及び2ΔKDだけ各制御パラメータを減少させる。
【0035】
こうして、上記評価値fを最小化するよう、制御パラメータ記憶部59に記憶されている各制御パラメータを順次更新するようにしている。これにより、事前のパラメータ設定を詳細にしなくても、各種船舶において好適に挙動を制御できる制御パラメータを次第に得ることができる。なお、以上の説明では、実際の挙動に対して最も適合度が高かった挙動パターン、いわゆる勝者パターンに対応する制御パラメータのみを更新する方式を説明したが、本発明はこれに限定されるものではない。
【0036】
ここで、本実施形態に係る自動操舵制御装置10の動作について説明する。
【0037】
図9及び図10は、本実施形態に係る自動操舵制御装置10の動作を示すフロー図である。同図に示すように、この自動操舵制御装置10では、まず各挙動を特定する数値を格納する変数nをメモリに確保し、そこに初期値1を設定する(図9;S101)。次に、予め用意された制御パラメータ(比例係数KP、積分係数KI及び微分係数KD)を操舵量決定部22に供給し、それにより挙動Y1(Y:Yawing)の間、舵機16を制御する(S102)。そして、挙動検出部52は、偏差θ−θにより挙動Y1の終了タイミングを監視する(S103)。そして、挙動検出部52が挙動Y1の終了タイミングを検知すると、それを挙動特徴量算出部56に通知し、それにより挙動特徴量算出部56が挙動Y1の特徴量(θc,ωc,Δθ,Δω)を算出する(S104)。
【0038】
こうして、算出される特徴量はパターンマッチング部58及び挙動評価値算出部54に供給され、パターンマッチング部58は挙動Y1が最も当てはまる挙動パターンR1を決定する(S105)。また、挙動評価値算出部54は挙動Y1の評価値f1を算出する(S106)。そして、パターンマッチング部58は、S105において決定された挙動パターンR1に対応する制御パラメータを制御パラメータ記憶部59から読み出し、制御パラメータ設定部62は、それを上述した更新量ΔKP,ΔKI及びΔKDだけ増加させる。そして、こうして更新した制御パラメータを挙動パターンR1に対応づけて制御パラメータ記憶部59に記憶させる。すなわち、挙動パターンR1に対応する制御パラメータを更新する。その後、パターンマッチング部58により更新後の制御パラメータに基づいて再度ファジィ推論された制御パラメータを取得し、それを操舵量決定部22に供給する(S107)。
【0039】
その後、操舵量決定部22は、S107において制御パラメータ設定部62から供給された制御パラメータを用いて、挙動Y2の間、舵機16を制御する(図10;S108)。そして、挙動検出部52は、偏差θ−θにより挙動Y2の終了タイミングを監視する(S109)。そして、挙動検出部52が挙動Y2の終了タイミングを検知すると、それを挙動特徴量算出部56に通知し、それにより挙動特徴量算出部56が挙動Y2の特徴量(θc,ωc,Δθ,Δω)を算出する(S110)。こうして、算出される特徴量はパターンマッチング部58及び挙動評価値算出部54に供給され、パターンマッチング部58は挙動Y2が最も当てはまる挙動パターンR2を決定する(S111)。また、挙動評価値算出部54は挙動Y2の評価値f2を算出する(S112)。
【0040】
その後、制御パラメータ更新部60は、挙動Y1の評価値f1よりも挙動Y2の評価値f2の方が小さくなっているか否かを判断する(S113)。そして、小さくなっていれば、挙動Y1における制御が良好に進んだとして、挙動パターンR1に対応する制御パラメータをさらに更新量ΔKP,ΔKI及びΔKDだけ増加させるよう、制御パラメータ記憶部59の記憶内容を更新する(S114)。一方、S113において、挙動Y2の評価値f2が挙動Y1の評価値f1以上であると判断すると、挙動パターンR1に対応する制御パラメータを更新量2ΔKP,2ΔKI及び2ΔKDだけ減少させるよう、制御パラメータ記憶部59の記憶内容を更新する(S115)。
【0041】
その後、S111での決定に従い、挙動パターンR2に対応する制御パラメータを制御パラメータ記憶部59から読み出し、制御パラメータ設定部62は、それを上述した更新量ΔKP,ΔKI及びΔKDだけ増加させる。また、制御パラメータ設定部62は、こうして更新した制御パラメータを、挙動パターンR2に対応づけて制御パラメータ記憶部59に記憶させる。すなわち、挙動パターンR2に対応する制御パラメータを更新する。そして、パターンマッチング部58により更新後の制御パラメータに基づいて再度ファジィ推論された制御パラメータを取得し、それを操舵量決定部22に供給する(S116)。その後、変数nを1だけ増加し(S117)、S108の処理に戻る。こうして、以降はS108乃至S117の処理を繰り返し、後続の挙動Ynにおいて舵機16を制御していく。
【0042】
以上説明した自動操舵装置10によれば、目標値である基準方位θと制御量である船首方位θとの偏差が増減を繰り返す場合に、その増減パターンに対応する制御パラメータを用いて、船舶の制御を好適に行うことができる。また、各挙動に対する評価を行い、それにより制御パラメータを自動更新して適応制御を行うようにしているので、複雑なパラメータ調整作業を経ることなく、様々な船舶の制御を行うことができるようになる。
【0043】
なお、本発明は上記実施形態に限定されるものではない。例えば、ここでは船舶の操舵制御に本発明を適用する例を示したが、他の移動体にも本発明を同様に適用することができる。また、移動体制御以外にも本発明を同様に適用することができる。さらに、移動体の進行方向のみならず、姿勢制御や速度制御等にも適用することができ、移動体以外であっても、温度や濃度等の物理量の制御に適用することができる。また、以上の説明では制御パラメータKP,KI,KDを同時に更新するようにしたが、所定挙動数毎に順に、1種類ずつ順に制御パラメータを更新するようにすれば、制御パラメータ毎の制御に与える影響の違いを反映させた、さらに適切な学習を行うことができる。
【0044】
【発明の効果】
以上説明したように、本発明では、所定種類の挙動の時間範囲を順次判断し、それら時間範囲において取得される制御量に基づいて、各時間範囲において前記制御対象により行われた前記所定種類の挙動が当てはまる挙動パターンを複数の挙動パターンの中から順次判断している。そして、前記所定種類の挙動が当てはまる挙動パターンに対応づけて記憶手段に記憶される制御パラメータを順次読み出し、それにより前記制御対象を制御している。こうして、制御対象の挙動を判別し、それを考慮して効率的な制御を行うことができる。
【図面の簡単な説明】
【図1】本発明の実施形態に係る自動操舵制御装置の構成を示す図である。
【図2】操舵量決定部の構成を示す図である。
【図3】デッドバンド処理部の処理内容を説明する図である。
【図4】制御パラメータ決定部の構成を示す図である。
【図5】方位差(実方位と目標方位との差)の経時変化を示す図である。
【図6】挙動あたりの方位差とその一階微分の推移を位相面で示す図である。
【図7】メンバーシップ関数群の一例を示す図である。
【図8】ファジィルールテーブル群の一例を示す図である。
【図9】本発明の実施形態に係る自動操舵制御装置の動作を示すフロー図である。
【図10】本発明の実施形態に係る自動操舵制御装置の動作を示すフロー図である。
【符号の説明】
10 自動操舵制御装置、12 針路設定部、14 方位センサ、15 舵角センサ、16 舵機、18,24,40 加算器、20 制御パラメータ決定部、22 操舵量決定部、26 デッドバンド処理部、30 積分器、32 積分パラメータ乗算部、34 比例パラメータ乗算部、36 微分器、38 微分パラメータ乗算部、52 挙動検出部、54 挙動評価値算出部、56 挙動特徴量算出部、58 パターンマッチング部、60 制御パラメータ更新部、62 制御パラメータ設定部。

Claims (8)

  1. 制御対象に関する制御量を制御する制御装置であって、
    制御量を順次取得する制御量取得手段と、
    前記制御量取得手段により取得される制御量に基づいて、前記制御対象により行われる所定種類の挙動の時間範囲を順次判断する時間範囲判断手段と、
    前記時間範囲判断手段により順次判断される時間範囲において前記制御量取得手段により取得される制御量に基づいて、各時間範囲において前記制御対象により行われた前記所定種類の挙動が当てはまる挙動パターンを複数の挙動パターンの中から順次判断する挙動パターン判断手段と、
    前記複数の挙動パターンの各々に対応づけて制御パラメータを記憶する制御パラメータ記憶手段と、
    前記挙動パターン判断手段により順次判断される、前記所定種類の挙動が当てはまる挙動パターンに対応づけて前記制御パラメータ記憶手段に記憶される制御パラメータを順次読み出す制御パラメータ読み出し手段と、
    前記制御パラメータ読み出し手段により順次読み出される制御パラメータに基づいて、前記制御対象を制御する制御手段と、
    を含むことを特徴とする制御装置。
  2. 請求項1に記載の制御装置において、
    前記時間範囲判断手段により順次判断される時間範囲において前記制御量取得手段により取得される制御量に基づいて、各時間範囲において前記制御対象により行われた前記所定種類の挙動に対する評価値を順次算出する挙動評価値算出手段と、
    前記挙動評価値算出手段により順次算出される評価値に基づいて、前記制御パラメータ記憶手段に記憶される制御パラメータを更新する制御パラメータ更新手段と、
    をさらに含むことを特徴とする制御装置。
  3. 請求項2に記載の制御装置において、
    前記制御パラメータ更新手段は、前記挙動評価値算出手段により順次算出される評価値に対応する前記所定種類の挙動が行われた時間範囲において用いられた制御パラメータを更新する
    ことを特徴とする制御装置。
  4. 請求項1乃至3のいずれかに記載の制御装置において、
    前記時間範囲判断手段は、前記制御量取得手段により取得される制御量に基づいて、前記制御量取得手段により取得される制御量が極値をとるタイミングを前記所定種類の挙動の時間範囲の始期及び終期と判断する
    ことを特徴とする制御装置。
  5. 請求項1乃至4のいずれかに記載の制御装置において、
    前記制御対象は船舶であり、
    前記制御量はその針路であり、
    前記制御手段は前記船舶に備えられた舵機を制御する
    ことを特徴とする制御装置。
  6. 制御対象に関する制御量を制御する制御方法であって、
    制御量を順次取得する制御量取得ステップと、
    前記制御量取得ステップで取得される制御量に基づいて、前記制御対象により行われる所定種類の挙動の時間範囲を順次判断する時間範囲判断ステップと、
    前記時間範囲判断ステップで順次判断される時間範囲において前記制御量取得ステップで取得される制御量に基づいて、各時間範囲において前記制御対象により行われた前記所定種類の挙動が当てはまる挙動パターンを複数の挙動パターンの中から順次判断する挙動パターン判断ステップと、
    前記複数の挙動パターンの各々に対応づけて制御パラメータを記憶する制御パラメータ記憶手段から、前記挙動パターン判断ステップで順次判断される、前記所定種類の挙動が当てはまる挙動パターンに対応づけて記憶される制御パラメータを順次読み出す制御パラメータ読み出しステップと、
    前記制御パラメータ読み出しステップで順次読み出される制御パラメータに基づいて、前記制御対象を制御する制御ステップと、
    を含むことを特徴とする制御方法。
  7. 制御対象に関する制御量を制御する制御装置であって、
    制御量を順次取得する制御量取得手段と、
    前記制御量取得手段により取得される制御量に基づいて、前記制御対象により行われる所定種類の挙動の時間範囲を順次判断する時間範囲判断手段と、
    前記時間範囲判断手段により順次判断される時間範囲において前記制御量取得手段により取得される制御量に基づいて、各時間範囲において前記制御対象により行われた前記所定種類の挙動が当てはまる挙動パターンを複数の挙動パターンの中から順次判断する挙動パターン判断手段と、
    前記挙動パターン判断手段により順次判断される、前記所定種類の挙動が当てはまる挙動パターンに対応する制御パラメータを順次取得する制御パラメータ取得手段と、
    前記制御パラメータ取得手段により順次取得される制御パラメータに基づいて、前記制御対象を制御する制御手段と、
    を含むことを特徴とする制御装置。
  8. 制御対象に関する制御量を制御する制御方法であって、
    制御量を順次取得する制御量取得ステップと、
    前記制御量取得ステップで取得される制御量に基づいて、前記制御対象により行われる所定種類の挙動の時間範囲を順次判断する時間範囲判断ステップと、
    前記時間範囲判断ステップで順次判断される時間範囲において前記制御量取得ステップで取得される制御量に基づいて、各時間範囲において前記制御対象により行われた前記所定種類の挙動が当てはまる挙動パターンを複数の挙動パターンの中から順次判断する挙動パターン判断ステップと、
    前記挙動パターン判断ステップで順次判断される、前記所定種類の挙動が当てはまる挙動パターンに対応する制御パラメータを順次取得する制御パラメータ取得ステップと、
    前記制御パラメータ取得ステップで順次取得される制御パラメータに基づいて、前記制御対象を制御する制御ステップと、
    を含むことを特徴とする制御方法。
JP2003097319A 2003-03-31 2003-03-31 船舶等の制御装置及び方法 Expired - Lifetime JP3683890B2 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2003097319A JP3683890B2 (ja) 2003-03-31 2003-03-31 船舶等の制御装置及び方法
GB0406640A GB2400189B (en) 2003-03-31 2004-03-24 Control system and method for controlling a moving body
US10/809,399 US7403841B2 (en) 2003-03-31 2004-03-26 Control system and method for controlling a moving body

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2003097319A JP3683890B2 (ja) 2003-03-31 2003-03-31 船舶等の制御装置及び方法

Publications (2)

Publication Number Publication Date
JP2004303085A true JP2004303085A (ja) 2004-10-28
JP3683890B2 JP3683890B2 (ja) 2005-08-17

Family

ID=32212186

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2003097319A Expired - Lifetime JP3683890B2 (ja) 2003-03-31 2003-03-31 船舶等の制御装置及び方法

Country Status (3)

Country Link
US (1) US7403841B2 (ja)
JP (1) JP3683890B2 (ja)
GB (1) GB2400189B (ja)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009179263A (ja) * 2008-01-31 2009-08-13 Yokogawa Denshikiki Co Ltd 自動操舵装置
JP2011235771A (ja) * 2010-05-11 2011-11-24 Furuno Electric Co Ltd 操舵制御装置、船舶、操舵制御方法、および操舵制御プログラム
JP2012071621A (ja) * 2010-09-27 2012-04-12 Tokyo Keiki Inc 船舶用自動操舵装置
JPWO2020246568A1 (ja) * 2019-06-06 2020-12-10

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3677274B2 (ja) * 2003-03-31 2005-07-27 財団法人ファジィシステム研究所 制御装置及び方法
GB2442020A (en) * 2006-09-20 2008-03-26 Northrop Gruman Sperry Marine Control system for a marine vessel
CN105675195B (zh) * 2016-02-03 2018-04-24 中国海洋石油总公司 基于船舶姿态测量的工程船作业实时分析系统
JP7531377B2 (ja) * 2020-12-03 2024-08-09 日本発條株式会社 船舶、船首方位制御方法、船首方位制御装置およびプログラム
CN112578801B (zh) * 2021-03-01 2021-06-01 北京海兰信数据科技股份有限公司 一种船舶航向抗干扰控制方法
CN113238488B (zh) * 2021-07-12 2021-10-12 北京海兰信数据科技股份有限公司 一种船舶模型参数的获得方法及装置
JP2023135995A (ja) * 2022-03-16 2023-09-29 古野電気株式会社 自動操舵装置、自動操舵システム、自動操舵方法、および自動操舵プログラム

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3517285A (en) * 1968-09-04 1970-06-23 Anschuetz & Co Gmbh Electrical automatic pilot for rudder-controlled vehicles,particularly for ships
US3656043A (en) * 1969-02-08 1972-04-11 Tokyo Keika Seizosho Kk Steering system with automatic increase of rate time
JPS5030915B1 (ja) * 1969-08-30 1975-10-04
US3604907A (en) * 1970-03-02 1971-09-14 Sperry Rand Corp Steering command computer for navigable craft
IT973336B (it) * 1972-02-01 1974-06-10 Licentia Gmbh Regolatore di rotta elettrico per navi
DE3486373T2 (de) * 1983-11-04 1995-06-29 Nissan Motor Fahrzeugmotorsteuersystem mit der Fähigkeit den Betriebszustand des Motors zu vermitteln und das passende Betriebsschema zu wählen.
GB2173019B (en) * 1983-12-12 1988-08-24 Foxboro Co Pattern-recognizing self-tuning controller
JPH01153664A (ja) 1987-12-11 1989-06-15 Ajinomoto Co Inc 新規アミノ酸誘導体及びその製造、使用
US5014206A (en) * 1988-08-22 1991-05-07 Facilitech International Incorporated Tracking system
US6546363B1 (en) * 1994-02-15 2003-04-08 Leroy G. Hagenbuch Apparatus for tracking and recording vital signs and task-related information of a vehicle to identify operating patterns
JP3911056B2 (ja) 1996-12-20 2007-05-09 株式会社トキメック 操船支援装置
JP4660870B2 (ja) 1999-11-08 2011-03-30 ソニー株式会社 脚式移動ロボット及びその制御方法
JP2002298150A (ja) 2001-03-30 2002-10-11 Dainippon Printing Co Ltd 画像生成装置および画像生成方法

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009179263A (ja) * 2008-01-31 2009-08-13 Yokogawa Denshikiki Co Ltd 自動操舵装置
JP2011235771A (ja) * 2010-05-11 2011-11-24 Furuno Electric Co Ltd 操舵制御装置、船舶、操舵制御方法、および操舵制御プログラム
JP2012071621A (ja) * 2010-09-27 2012-04-12 Tokyo Keiki Inc 船舶用自動操舵装置
JPWO2020246568A1 (ja) * 2019-06-06 2020-12-10
WO2020246568A1 (ja) * 2019-06-06 2020-12-10 日本発條株式会社 自動設定装置、自動設定方法およびプログラム
JP7157945B2 (ja) 2019-06-06 2022-10-21 日本発條株式会社 自動設定装置、自動設定方法およびプログラム
US11958583B2 (en) 2019-06-06 2024-04-16 Nhk Spring Co., Ltd. Automatic setting device, automatic setting method, and program

Also Published As

Publication number Publication date
US20040193332A1 (en) 2004-09-30
US7403841B2 (en) 2008-07-22
JP3683890B2 (ja) 2005-08-17
GB0406640D0 (en) 2004-04-28
GB2400189A (en) 2004-10-06
GB2400189B (en) 2006-07-05

Similar Documents

Publication Publication Date Title
CN108549208B (zh) 一种基于因子自适应模糊pid的四旋翼飞行器姿态控制方法
CN105785999B (zh) 无人艇航向运动控制方法
JP2669031B2 (ja) 自律走行車両
JP2004303085A (ja) 船舶等の制御装置及び方法
CN111766777B (zh) 一种pid控制器及pid控制方法
JP2004303086A (ja) 制御装置及び方法、並びに制御状態判定装置及び方法
CN108983774B (zh) 一种基于模糊状态观测器的单喷泵推进无人水面艇自适应航向控制方法
CN113031614A (zh) 一种远洋船舶航向控制复合优化节油方法
Sutton et al. A design study of a self-organizing fuzzy autopilot for ship control
CN111061262B (zh) 一种减少无效操舵的航向控制方法
JP2006321455A (ja) 船舶用自動操舵装置
JP4543179B2 (ja) 移動体等の制御装置及び方法
JP3346006B2 (ja) Rovの自動制御方法
JP3629760B2 (ja) Rovの自動制御方法
JP4897450B2 (ja) 船舶用自動操舵装置
CN113741183A (zh) 基于阻尼比模型的工业机器人自适应导纳控制方法
WO2021149435A1 (ja) 自動操縦ロボットの制御装置及び制御方法
JP2002258905A (ja) 自動制御装置の制御則決定方法
JP2800582B2 (ja) 水中航走体の姿勢制御方式
JPH07210202A (ja) 水中ロボット姿勢制御装置
JPH05284771A (ja) Pi制御器のオートチューニング方法
JPH07165180A (ja) 水中航走体の深度制御方式
KR102605858B1 (ko) Pid 제어 방법 및 장치
JPH057238B2 (ja)
JPH045598B2 (ja)

Legal Events

Date Code Title Description
A871 Explanation of circumstances concerning accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A871

Effective date: 20041203

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20050118

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20050318

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20050517

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20050526

R150 Certificate of patent or registration of utility model

Ref document number: 3683890

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313113

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080603

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090603

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090603

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100603

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100603

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110603

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120603

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130603

Year of fee payment: 8

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140603

Year of fee payment: 9

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

EXPY Cancellation because of completion of term