(57)【要約】【課題】操作者の感性をフィードバック制御できる制御装置を提供する。【解決手段】感性メータ２０は、対象機器１０の出力ｕ（ｔ）に対応する操作者の生体情報ｘ（ｔ）を検出し、当該生体情報ｘ（ｔ）に基づき、操作者の快適度ｙ（ｔ）を決定する。第１制御部３０は、快適度ｙ（ｔ）に関する第１目標値ｒ（ｔ）と、快適度ｙ（ｔ）との差分ｅ（ｔ）に基づき、出力ｕ（ｔ）に関する第２目標値ｗ（ｔ）を決定する。第２制御部４０は、第２目標値ｗ（ｔ）と出力ｕ（ｔ）との差分ε（ｔ）に基づき、対象機器１０に対する制御入力ｖc （ｔ）を決定する。δ設定部５０は、操作者による対象機器１０に対する操作入力ｖh （ｔ）、及び、制御入力ｖc （ｔ）のそれぞれについて、操作者の操作レベルに応じた重み付けを行う。重み付けされた操作入力ｖh （ｔ）及び制御入力ｖc （ｔ）が合算されて対象機器１０に入力される。【選択図】図１
2020071766.tif 000002
【特許請求の範囲】
【請求項１】
操作者によって操作される対象機器と、
前記対象機器の出力ｕ（ｔ）に対応する前記操作者の生体情報ｘ（ｔ）を検出し、当該生体情報ｘ（ｔ）に基づき、前記操作者の快適度ｙ（ｔ）を決定する感性メータと、
前記快適度ｙ（ｔ）に関する第１目標値ｒ（ｔ）と、前記快適度ｙ（ｔ）との差分に基づき、前記出力ｕ（ｔ）に関する第２目標値ｗ（ｔ）を決定する第１制御部と、
前記第２目標値ｗ（ｔ）と前記出力ｕ（ｔ）との差分に基づき、前記対象機器に対する制御入力ｖc （ｔ）を決定する第２制御部と、
前記操作者による前記対象機器に対する操作入力ｖh （ｔ）、及び、前記制御入力ｖc （ｔ）のそれぞれについて、前記操作者の操作レベルに応じた重み付けを行う重み付け設定部と
を備え、
前記重み付け設定部によりそれぞれ重み付けされた前記操作入力ｖh （ｔ）及び前記制御入力ｖc （ｔ）が合算されて前記対象機器に入力される、
感性フィードバック制御装置。
【請求項２】
請求項１に記載の感性フィードバック制御装置において、
前記感性メータは、快・不快、活性・非活性、期待感の３つの情報を前記生体情報ｘ（ｔ）として検出し、当該３つの情報の相関性から得られる感性値を前記快適度ｙ（ｔ）として決定する、
感性フィードバック制御装置。
【請求項３】
請求項１又は２に記載の感性フィードバック制御装置において、
前記第１目標値ｒ（ｔ）、前記第２目標値ｗ（ｔ）及び前記快適度ｙ（ｔ）が逐次蓄積されるデータベースをさらに備え、
前記第１制御部は、前記データベースに蓄積されているデータを用いて、制御パラメータの調整を行いながら、前記第２目標値ｗ（ｔ）を決定する、
感性フィードバック制御装置。
【請求項４】
請求項１〜３のいずれか１項に記載の感性フィードバック制御装置において、
前記対象機器は、建設機械である、
感性フィードバック制御装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、感性フィードバック制御装置に関するものである。
【背景技術】
【０００２】
内閣府の調査によれば、我が国は、国内総生産（ＧＤＰ）が高いにも関わらず、幸福度が低いとされている。すなわち、ＧＤＰに関する「物の豊かさ」と、幸福度に関する「心の豊かさ」との間には大きなギャップが存在している。
【０００３】
このようなギャップを埋めるための方策の１つとして、既に高度化されている「物」（福祉支援機器など）が、人の感性を考慮した動作を行うことによって、心の豊かさを向上させることが考えられている（例えば特許文献１）。また、「感性」の可視化技術についても、社会実装を想定した研究が行われている（例えば特許文献２）。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１８−３６７７３号公報
【特許文献２】特開２０１７−７４３５６号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、感性に関するほとんどの研究は、製品のデザイン評価・設計などの静的な分野であり、特に操作者によって操作される自動車や建設機械等の機器を対象として操作者の感性を動的に制御する研究はほとんど行われていない。
【０００６】
そこで、本発明は、操作者によって操作される機器において操作者の感性をフィードバック制御できる感性フィードバック制御装置を提供することを目的とする。
【課題を解決するための手段】
【０００７】
前記の目的を達成するために、本願発明者らは、人の感性は、時変系や非線形系であると考えられるため、そのモデル化は困難であることに着目し、以下のような発明を想到した。すなわち、操作者によって操作される対象機器の出力（例えば油圧ショベルの応答速度）に対応する人の感性を計測し、制御対象値として「快適度」を設定し、「快適度」が目標値を満足するように対象機器の制御を行うという発明である。ここで、「快適度」は、操作者が、脳内の対象機器の目標出力（例えば油圧ショベルの目標応答速度）と同様な出力が得られるように現実の対象機器の操作ができたときに、最も高い値を持つ。
【０００８】
以上のように、本発明は、(i) 操作者のモデルを必要とせず、(ii)感性情報がフィードバックされることから、本発明を「感性フィードバック制御」と称するものとする。また、本発明の基本構成は、対象機器に関する制御系である内側ループに、感性フィードバック制御系である外側ループを付加したカスケード制御系である。
【０００９】
具体的には、本発明に係る感性フィードバック制御装置は、操作者によって操作される対象機器と、対象機器の出力ｕ（ｔ）に対応する操作者の生体情報ｘ（ｔ）を検出し、当該生体情報ｘ（ｔ）に基づき、操作者の快適度ｙ（ｔ）を決定する感性メータと、快適度ｙ（ｔ）に関する第１目標値ｒ（ｔ）と、快適度ｙ（ｔ）との差分に基づき、出力ｕ（ｔ）に関する第２目標値ｗ（ｔ）を決定する第１制御部と、第２目標値ｗ（ｔ）と出力ｕ（ｔ）との差分に基づき、対象機器に対する制御入力ｖc （ｔ）を決定する第２制御部と、操作者による対象機器に対する操作入力ｖh （ｔ）、及び、制御入力ｖc （ｔ）のそれぞれについて、操作者の操作レベルに応じた重み付けを行う重み付け設定部とを備え、重み付け設定部によりそれぞれ重み付けされた操作入力ｖh （ｔ）及び制御入力ｖc （ｔ）が合算されて対象機器に入力される。
【００１０】
本発明に係る感性フィードバック制御装置によると、感性メータが、対象機器の出力ｕ（ｔ）に対応する操作者の生体情報ｘ（ｔ）を検出し、当該生体情報ｘ（ｔ）に基づき、操作者の快適度ｙ（ｔ）を決定する。また、第１制御部は、快適度ｙ（ｔ）に関する第１目標値ｒ（ｔ）と、快適度ｙ（ｔ）との差分に基づき、出力ｕ（ｔ）に関する第２目標値ｗ（ｔ）を決定する。また、第２制御部は、第２目標値ｗ（ｔ）と出力ｕ（ｔ）との差分に基づき、対象機器に対する制御入力ｖc （ｔ）を決定する。このため、感性メータ、第１制御部及び第２制御部を用いて、対象機器において操作者の感性（具体的には快適度）をフィードバックした制御を行うことができる。
【００１１】
また、本発明に係る感性フィードバック制御装置によると、重み付け設定部は、操作者による対象機器に対する操作入力ｖh （ｔ）、及び、制御入力ｖc （ｔ）のそれぞれについて、操作者の操作レベルに応じた重み付けを行い、重み付けされた操作入力ｖh （ｔ）及び制御入力ｖc （ｔ）が合算されて対象機器に入力される。このため、重み付け設定部を用いて、感性フィードバックの程度を操作者の操作レベルに応じて設定することができる。これにより、例えば、未熟な操作者に対しては感性フィードバックの程度を高くして、対象機器における当該操作者の操作入力の反映度を低くしたり、逆に、熟練操作者に対しては感性フィードバックの程度を低くして、対象機器における当該操作者の操作入力の反映度を高くしたりすることができる。
【００１２】
また、本発明に係る感性フィードバック制御装置において、感性メータは、快・不快、活性・非活性、期待感の３つの情報を生体情報ｘ（ｔ）として検出し、当該３つの情報の相関性から得られる感性値を快適度ｙ（ｔ）として決定してもよい。
【００１３】
このようにすると、感性である「快適度」を定量的に評価できるので、感性フィードバック制御の精度を向上させることができる。
【００１４】
また、本発明に係る感性フィードバック制御装置において、第１目標値ｒ（ｔ）、第２目標値ｗ（ｔ）及び快適度ｙ（ｔ）が逐次蓄積されるデータベースをさらに備え、第１制御部は、データベースに蓄積されているデータを用いて、制御パラメータの調整を行いながら、第２目標値ｗ（ｔ）を決定してもよい。
【００１５】
このようにすると、非線形系である操作者の感性をモデル化することなく、データベース駆動型制御を用いて、感性フィードバック制御を行うことができる。
【００１６】
また、本発明に係る感性フィードバック制御装置において、対象機器は、建設機械であってもよい。
【００１７】
このようにすると、操作者が、脳内の建設機械の目標出力（例えば油圧ショベルの目標応答速度）と同様な出力が得られるように現実の建設機械を操作することが可能になる。
【発明の効果】
【００１８】
本発明によると、操作者によって操作される機器において操作者の感性をフィードバック制御できる感性フィードバック制御装置を提供することができる。
【図面の簡単な説明】
【００１９】
【図１】実施形態に係る感性フィードバック制御装置の構成図である。
【図２】実施形態に係る感性フィードバック制御装置において用いられる感性多軸モデルの模式図である。
【図３】実施形態に係る感性フィードバック制御装置において用いられるデータベース駆動型制御を説明する図である。
【図４】実施形態に係る感性フィードバック制御装置によって得られる効果を示す図である。
【図５】図１に示す感性フィードバック制御装置においてδ値が１である場合の構成を簡単化して示す図である。
【図６】実施形態に係る感性フィードバック制御装置において用いられるＦＲＩＴのブロック線図である。
【図７】実施形態に係る感性フィードバック制御装置における快適度と出力誤差との関係を示す図である。
【図８】実施形態に係る感性フィードバック制御装置における快適度、出力及び入力の時間変化をシミュレーションした結果を示す図である。
【図９】実施形態に係る感性フィードバック制御装置における制御パラメータの時間変化をシミュレーションした結果を示す図である。
【図１０】比較例に係る感性フィードバック制御装置の構成図である。
【図１１】比較例に係る感性フィードバック制御装置における快適度、出力及び入力の時間変化をシミュレーションした結果を示す図である。
【発明を実施するための形態】
【００２０】
以下、図面を参照しながら、本発明の実施形態に係る感性フィードバック制御装置について説明する。尚、本発明の範囲は、以下の実施の形態に限定されず、本発明の技術的思想の範囲内で任意に変更可能である。
【００２１】
図１は、一実施形態に係る感性フィードバック制御装置の構成図である。図１に示すように、感性フィードバック制御装置１００は、対象機器１０と、感性メータ２０と、第１制御部３０と、第２制御部４０と、重み付け（δ）設定部５０とを備えている。対象機器１０は、操作者によって操作される機器、例えば、自動車や建設機械等である。感性メータ２０は、対象機器１０の出力ｕ（ｔ）に対応する操作者の生体情報ｘ（ｔ）を検出し、当該生体情報ｘ（ｔ）に基づき、操作者の快適度ｙ（ｔ）を決定する。対象機器１０の出力ｕ（ｔ）とは、例えば、油圧ショベルの応答速度である。第１制御部３０は、快適度ｙ（ｔ）に関する第１目標値（目標快適度）ｒ（ｔ）と、快適度ｙ（ｔ）との差分ｅ（ｔ）に基づき、出力ｕ（ｔ）に関する第２目標値（目標出力）ｗ（ｔ）を決定する。第２制御部４０は、目標出力ｗ（ｔ）と出力ｕ（ｔ）との差分ε（ｔ）に基づき、対象機器１０に対する制御入力ｖc （ｔ）を決定する。δ設定部５０は、操作者による対象機器１０に対する操作入力ｖh （ｔ）、及び、制御入力ｖc （ｔ）のそれぞれについて、操作者の操作レベルに応じた重み付けを行う。
【００２２】
（感性メータ）
本実施形態では、対象機器の出力ｕ（ｔ）に対応する操作者の生体情報ｘ（ｔ）に基づく快適度ｙ（ｔ）が制御パラメータに用いられる。感性メータ２０は、操作者の生体情報ｘ（ｔ）を検出する手段として、例えば、ＣＣＤカメラによる顔表情の検出、心拍センサによる心拍の検出、ＣＣＤカメラによる呼吸（回数や深さ）の検出、皮膚インピーダンスセンサによる皮膚抵抗の検出、筋電センサによる上肢筋電、下肢筋電の検出、マイクによる音声情報の検出等から適宜選択される一つ或いは二つ以上の手段を用いてもよい。
【００２３】
ところで、人は何かを見たり、聞いたり、或いは何かに触れたり、触れられたりしたときに、わくわくしたり、うきうきしたり、はらはらしたり、どきどきしたりといった感情或いは情動を抱くことがある。そのような感情或いは情動は人の複雑で高次の脳活動によってもたらされると考えられ、感情或いは情動の形成には、運動神経及び感覚神経を含む体性神経系、交感神経及び副交感神経を含む自律神経系、さらには記憶や経験などが深く関与していると考えられる。そこで、感性を、外受容感覚情報（体性神経系）と内受容感覚情報（自律神経系）を統合し、過去の経験、記憶と照らし合わせて生じる情動反応を、より上位のレベルで俯瞰する高次脳機能と定義する。換言すると、感性は、予測（イメージ）と結果（感覚情報）とのギャップを経験・知識と比較することによって“はっ”と気付く高次脳機能であると言える。
【００２４】
このような高次脳機能である感性は、種々の観点から総合的に捉える必要がある。例えば、人が快い、快適、或いは心地よいと感じているか、反対に人が気持ち悪い、不快、或いは心地よくないと感じているかといった「快／不快」の観点から感性を捉えることができる。また、例えば、人が覚醒、興奮、或いは活性状態にあるか、反対に人がぼんやり、沈静、或いは非活性状態にあるかといった「活性／非活性」の観点から感性を捉えることができる。また、例えば、人が何かを期待或いは予期してわくわくしているか、反対にわくわく感が阻害されているかといった「期待感」の観点から感性を捉えることができる。
【００２５】
快／不快及び活性／非活性を２軸に表したラッセル（Russell）の円環モデルが知られており、この円環モデルで感情を表すことができる。しかし、感性は予測（イメージ）と結果（感覚情報）とのギャップを経験・知識と比較する高次脳機能であるので、快／不快及び活性／非活性の２軸からなる既存の円環モデルでは十分に表し得ない。そこで、本実施形態では、例えば、ラッセルの円環モデルに時間軸（例えば期待感）を第３軸として加えた感性多軸モデルとして、感性をとらえる。
【００２６】
図２は、本実施形態において一例として採用する感性多軸モデルの模式図である。図２に示すように、感性多軸モデルは、「快／不快」を第１軸、「活性／非活性」を第２軸、「時間（期待感）」を第３軸として表すことができる。感性を多軸モデル化することのメリットは、各軸に評価値を算出し、それらを総合することで、漠然と広い概念の感性を定量的に評価する、つまり、可視化することができる点にある。具体的には、感性多軸モデルの各軸の脳生理情報から、各軸の脳生理指標値（EEG快，EEG活性，EEG期待感）を求め、被験者の主観的な統計データから得られた、感性多軸モデルの各軸の重み付け係数（a，b，c）を示す主観心理軸を用いて、次のような計算式で、感性値を評価できることが実証されている。
【００２７】
感性値＝[主観心理軸]＊[脳生理指標]＝a＊EEG快＋b＊EEG活性＋c＊EEG期待感
感性メータ２０は、例えば、このような計算式に従って得られた感性値を、快適度ｙ（ｔ）として用いてもよい。
【００２８】
（第１制御部）
第１制御部３０の入出力関係については、例えば、ＫP 、ＫI 、ＫD をそれぞれＰＩＤパラメータとして、下記の式（１）に示すＰＩＤ制御がなされる。
【００２９】
【数１】
2020071766.tif 000003
【００３０】
本実施形態では、快適度ｙ（ｔ）を向上させることが目的となるが、フィードバック制御において個々に適した目標出力ｗ（ｔ）を設定することは困難である。そこで、図１に示すようなカスケード制御系を採用することにより、目標快適度ｒ（ｔ）を与えることによって、個々に適した目標出力ｗ（ｔ）を自動生成する。このとき、人の感性は時変系、非線形系であると考えられるため、第１制御部３０は、図１に示すように、データベース６０を利用したフィードバック制御（データベース駆動型制御：例えば「脇谷伸 他、FRIT法を用いた非線形PID制御系の設計、計測と制御、Vol.52、No.10、pp.885-891 (2013)」参照）を行う。データベース６０には、第１制御部３０のＰＩＤパラメータの調整に必要なデータ、例えば、目標快適度ｒ（ｔ）、快適度ｙ（ｔ）及び目標出力ｗ（ｔ）等が逐次蓄積される。
【００３１】
図３は、第１制御部３０のデータベース駆動型制御を説明する図である。第１制御部３０のデータベース駆動型制御については、後記の実施例で詳述するが、例えば、勾配法を用いた学習機能を利用してもよいし、さらには、ＦＲＩＴ（Fictitious Reference Iterative Tuning）法を組み合わせてオフライン学習を行ってもよい。
【００３２】
（第２制御部）
第２制御部４０の入出力関係については、例えば、ＫP ’、ＫI ’、ＫD ’をそれぞれＰＩＤパラメータとして、下記の式（２）に示すＰＩＤ制御がなされる。
【００３３】
【数２】
2020071766.tif 000004
【００３４】
第２制御部４０のＰＩＤパラメータには、対象機器１０に応じた既存のパラメータを用いることが可能である。既存のパラメータは、固定した定数であってもよいし、或いは、可変パラメータであってもよい。
【００３５】
（δ設定部）
δ設定部５０は、基本的には、スキルの高い操作者には相対的に低い重み（δ値）を設定し、スキルの低い操作者には相対的に高いδ値を設定する。例えば、操作者のスキルレベルを５段階で予め評価しておき、δ設定部５０は、各スキルレベルに応じたδ値を定数として保持しておき、操作者のスキルレベルが入力されると、δ値を一意に決定してもよい。或いは、δ設定部５０は、操作者にテスト操作を行わせた結果に基づき、当該操作者の操作レベルを評価し、その評価結果に応じたδ値を設定してもよい。
【００３６】
（効果）
以上に説明したように、本実施形態の感性フィードバック制御装置１００によると、感性メータ２０が、対象機器１０の出力ｕ（ｔ）に対応する操作者の生体情報ｘ（ｔ）を検出し、当該生体情報ｘ（ｔ）に基づき、操作者の快適度ｙ（ｔ）を決定する。また、第１制御部３０は、快適度ｙ（ｔ）に関する第１目標値（目標快適度）ｒ（ｔ）と、快適度ｙ（ｔ）との差分に基づき、出力ｕ（ｔ）に関する第２目標値（目標出力）ｗ（ｔ）を決定する。また、第２制御部４０は、目標出力ｗ（ｔ）と出力ｕ（ｔ）との差分に基づき、対象機器１０に対する制御入力ｖc （ｔ）を決定する。このように、感性メータ２０、第１制御部３０及び第２制御部４０を用いて、対象機器１０において操作者の感性（具体的には快適度ｙ（ｔ））をフィードバックした制御を行うことができる。
【００３７】
図４は、本実施形態の感性フィードバック制御装置１００によって得られる効果を示す図である。図４に示すように、対象機器１０、例えば油圧ショベルにおいて操作者の感性をフィードバックした制御を行うことができるので、元来操作者のスキルレベルと無関係に画一的に設計されている対象機器１０に、多種多様なスキルレベルを持つ操作者に応じて出力応答させることが可能となる。従って、多種多様なスキルレベルを持つ操作者の快適度を向上させることができる。
【００３８】
また、本実施形態の感性フィードバック制御装置１００によると、δ設定部５０は、操作者による対象機器１０に対する操作入力ｖh （ｔ）、及び、制御入力ｖc （ｔ）のそれぞれについて、操作者の操作レベルに応じた重み付けを行い、重み付けされた操作入力ｖh （ｔ）及び制御入力ｖc （ｔ）が合算されて対象機器１０に入力される。このため、δ設定部５０を用いて、感性フィードバックの程度を操作者の操作レベルに応じて設定することができる。これにより、例えば、未熟な操作者に対しては感性フィードバックの程度を高くして、対象機器１０における当該操作者の操作入力の反映度を低くしたり、逆に、熟練操作者に対しては感性フィードバックの程度を低くして、対象機器１０における当該操作者の操作入力の反映度を高くしたりすることができる。
【００３９】
また、本実施形態の感性フィードバック制御装置１００において、感性メータ２０は、快・不快、活性・非活性、期待感の３つの情報を生体情報ｘ（ｔ）として検出し、当該３つの情報の相関性から得られる感性値を快適度ｙ（ｔ）として決定してもよい。このようにすると、感性である「快適度」を定量的に評価できるので、感性フィードバック制御の精度を向上させることができる。
【００４０】
また、本実施形態の感性フィードバック制御装置１００において、目標快適度ｒ（ｔ）、目標出力ｗ（ｔ）及び快適度ｙ（ｔ）が逐次蓄積されるデータベース６０をさらに備え、第１制御部３０は、データベース６０に蓄積されているデータを用いて、制御パラメータの調整を行いながら、目標出力ｗ（ｔ）を決定してもよい。このようにすると、非線形系である操作者の感性をモデル化することなく、データベース駆動型制御を用いて、感性フィードバック制御を行うことができる。
【００４１】
また、本実施形態の感性フィードバック制御装置１００において、対象機器１０が建設機械であると、操作者が、脳内の建設機械の目標出力（例えば油圧ショベルの目標応答速度）と同様な出力が得られるように現実の建設機械を操作することが可能になる。
【００４２】
また、本実施形態の感性フィードバック制御装置１００において、第１制御部３０のデータベース駆動型制御が学習機能を利用するものであると、操作者による対象機器１０の操作時間が長くなるに従って、対象機器１０は、当該操作者にとってより適した出力応答を行うように徐々に変化する。これにより、操作者の快適度をより一層向上させることができる。
【００４３】
（実施例）
以下、感性フィードバック制御装置１００においてδ値を「１」に設定し、対象機器１０に制御入力ｖc （ｔ）のみが入力される場合を例として、主として第１制御部３０のデータベース駆動型制御について説明する。
【００４４】
図５は、図１に示す感性フィードバック制御装置１００においてδ値が１である場合の構成を簡単化して示す図である。尚、図５において、図１と同じ構成要素には同じ符号を付し、δ設定部５０等の図示を省略している。
【００４５】
本実施例では、対象機器１０は油圧ショベル、対象機器１０の出力ｕ（ｔ）は油圧ショベルのバケットの応答速度、対象機器１０に対する制御入力ｖc （ｔ）はトルクである。
【００４６】
また、本実施例では、制御対象である操作者の快適度ｙ（ｔ）を非線形系と捉え、データベース駆動型制御を適用して第１制御部３０を構成する。尚、内側ループの第２制御部４０は、ＫP ’＝１．５、ＫI ’＝０．１、ＫD ’＝０．１として、式（２）に示すＰＩＤ制御を行うものとし、本実施例では調整対象とはしない。
【００４７】
具体的には、制御対象である快適度ｙ（ｔ）として，下記の式（３）で表される離散時間非線形システムを想定する。
【００４８】
【数３】
2020071766.tif 000005
【００４９】
式（３）において、ｙ（ｔ）はシステム出力、ｆ（）は非線形関数、φ（ｔ−１）は、システムの時刻ｔより前の状態を表しており、情報ベクトルと呼ぶ。また、情報ベクトルφ（ｔ−１）は、下記の式（４）で定義される。
【００５０】
【数４】
2020071766.tif 000006
【００５１】
ここで，ｗ（ｔ）（対象機器１０の目標出力ｗ（ｔ））は、快適度ｙ（ｔ）の制御におけるシステム入力であり、ｎy 、ｎw はそれぞれシステム出力及びシステム入力の次数である。データベース駆動型制御では、各操業データが式（４）の形式でデータベース６０（図３参照）に蓄積される。また、現在のシステムの状態を表す情報ベクトルφ（ｔ）を要求点(クエリ)と呼ぶ。
【００５２】
また、第１制御部３０の制御則として、下記の式（５）で表される速度型Ｉ−ＰＤ制御則を用いる。
【００５３】
【数５】
2020071766.tif 000007
【００５４】
但し、式（５）において、ｅ（ｔ）は制御誤差信号であり、ｒ（ｔ）を目標快適度とすると、下記の式（６）で定義することができる。
【００５５】
【数６】
2020071766.tif 000008
【００５６】
また、式（５）において、ＫP （ｔ）、ＫI （ｔ）、ＫD （ｔ）はそれぞれ、各ステップ（時刻ｔ）における比例ゲイン、積分ゲイン、微分ゲインを表している。さらに、Δ（：＝１−ｚ-1）は差分演算子を、ｚ-1は遅延演算子を示している。
【００５７】
本実施例では、ＦＲＩＴ法を用いて、１回の動作実験によって得られた入出力データｗ0 （ｔ）、ｙ0 （ｔ）、及び、これらの実験データから生成される擬似参照入力〜ｒ（ｔ）によって、第１制御部３０の制御パラメータを直接的に算出する。図６は、本実施例で用いるＦＲＩＴのブロック線図を示す。図６において、Systemは、第２制御部４０と対象機器１０と感性メータ２０とからなるシステムであり、Ｃ（ｚ-1）は制御器（第１制御部３０）であり、下記の式（７）で表すことができる。
【００５８】
【数７】
2020071766.tif 000009
【００５９】
式（７）において、ｎは制御則の次数を示し、ＰＩＤ制御則の場合はｎ＝２となる。また、図６に示すように、Ｃ（ｚ-1）の入出力関係を下記の式（８）で表すことができる。
【００６０】
【数８】
2020071766.tif 000010
【００６１】
式（８）より、擬似参照入力〜ｒは、下記の式（９）のように算出される。
【００６２】
【数９】
2020071766.tif 000011
【００６３】
尚、ＦＲＩＴ法では、図６に示すように、予め所望の特性を有する参照モデルＧｍ（ｚ-1）を設計しておき、擬似参照入力〜ｒに対する参照モデルＧｍ（ｚ-1）の出力を〜ｙm （ｔ）とし、〜ｙm （ｔ）とｙ0 （ｔ）との誤差が小さくなるように、制御パラメータを算出する。参照モデルＧｍ（ｚ-1）は、下記の式（１０）、式（１１）の特性多項式で表すことができる。
【００６４】
【数１０】
2020071766.tif 000012
【００６５】
【数１１】
2020071766.tif 000013
【００６６】
式（１１）において、α、σはそれぞれ制御系の立ち上がり特性、減衰特性に関係するパラメータを示しており、任意に設定可能である。また、Ｔs はサンプリング時間である。具体的には、σは制御系の出力がステップ状の目標値の約６０％に達するまでの時間を示している。また、αは０以上２．０以下の範囲で設定されることが望ましく、α＝０ではBinomialモデルの応答を示し、α＝１ではButterworthモデルの応答を示す。
【００６７】
以下、データベース駆動型ＰＩＤコントローラの設計手順について説明する。
【００６８】
＜初期データベースの作成＞
データベース駆動型制御では、過去の蓄積データが存在しない場合、原理的に局所コントローラの設計を行うことができない。従って、本実施例では、ある平衡点周りで得られた入出力データから、Zieglar & Nichols（ＺＮ）法（「J.G.Zieglar 他、Optimum settings for automatic controllers、Trans.ASME、Vol.64、No.8、pp.759-768 (1942)」参照）や、Chien，Hrones & Reswick（ＣＨＲ）法（「K.L. Chien 他、On the Automatic Control of Generalized Passive Systems、Trans.ASME、Vol.74、pp.175-185 (1972)」参照）などを用いてＰＩＤゲインを算出し、これらのＰＩＤゲインと前述の入出力データとからなる情報ベクトル（下記の式（１２）で表される）によって初期データベースを作成する。
【００６９】
【数１２】
2020071766.tif 000014
【００７０】
式（１２）において、ｊ＝１、２、・・・、Ｎ（０）であり、ｉ＝１、２、・・・Ｍであり、−φ（ｊ）、Ｋ（ｊ）はそれぞれ下記の式（１３）、式（１４）で与えられる。
【００７１】
【数１３】
2020071766.tif 000015
【００７２】
【数１４】
2020071766.tif 000016
【００７３】
また、Ｎ（０）は初期データ数（初期データベースにおける情報ベクトルの数)を表し、Ｍは要素数を表す。初期データベースにおけるＰＩＤゲインは固定であるので、Ｋ（１）＝Ｋ（２）＝・・・＝Ｋ（Ｎ（０））である。
【００７４】
＜距離の算出、近傍の選択＞
要求点−φ（ｔ）と、データベースに蓄えられている情報ベクトル−φ（ｊ）との距離を、下記の式（１５）で表される重みつきＬ1 ノルムにより求める。
【００７５】
【数１５】
2020071766.tif 000017
【００７６】
ここで、Ｎ（ｔ）は時刻ｔにおいてデータベースに蓄えられているデータ数（情報ベクトル）を表している。また、−φl （ｊ）は第ｊ番目の情報ベクトルの第ｌ番目の要素を表している。同じく、−φl （ｔ）は時刻ｔにおける要求点の第ｌ番目の要素を表している。さらに、ｍａｘ−φl （ｍ）は、データベースに存在する全ての情報ベクトル（−φ（ｊ）：ｊ＝１、２、・・・、Ｎ（ｔ））の第ｌ番目の要素の中で最も大きな要素を示しており、ｍｉｎ−φl （ｍ）は、同第ｌ番目の要素のなかで最も小さな要素を示している。
【００７７】
本実施例では、式（１５）により求められた、距離ｄが小さい方からｋ個の情報ベクトルを選択し、当該選択されたデータ集合を近傍として定義する。
【００７８】
＜局所コントローラの構成＞
次に、前述のように選択された近傍に対して、下記の式（１６）で示される、重みつき局所線形平均法（Linearly Weighted Average：ＬＷＡ）により局所コントローラを構成する。
【００７９】
【数１６】
2020071766.tif 000018
【００８０】
ここで、ｗi は選択された第ｉ番目の情報ベクトルに含まれるＫ（ｉ）に対する重みであり、下記の式（１７）で与えられる。
【００８１】
【数１７】
2020071766.tif 000019
【００８２】
以上の手順により各時刻ｔにおけるＰＩＤゲインを算出することができる。さらに、データベース駆動型ＰＩＤ制御系が各平衡点において適切にＰＩＤゲインを調整できようにするためには、データベースの学習（制御パラメータの更新）を行う必要がある。そこで、本実施例では、ＦＲＩＴを適用し、データベースの構築に用いた初期データから、学習によってデータベース内の各データセットにおけるＰＩＤゲインをオフラインで更新する。
【００８３】
＜ＦＲＩＴを用いたデータベース駆動型ＰＩＤ制御のオフライン学習＞
以下、ＦＲＩＴを用いたデータベース駆動型ＰＩＤ制御のオフライン学習について，具体的に説明する。まず、閉ループデータでの要求点−φ0 （ｔ）におけるＰＩＤゲインを算出するために、式（１５）によって、要求点とデータベース内の情報ベクトルとの間の距離を計算し、ｋ個の近傍データを選択する。続いて、式（１６）によってＰＩＤゲインを算出し、算出したＰＩＤゲインを用い、下記の式（１８）、式（１９）で表される最急降下法によってＰＩＤゲインＫold （ｔ）の学習を行い、新たにＫnew を導出する。
【００８４】
【数１８】
2020071766.tif 000020
【００８５】
【数１９】
2020071766.tif 000021
【００８６】
式（１８）、式（１９）において、ηは学習係数、Ｊ（ｔ＋１）は、下記の式（２０）、式（２１）で定義される評価規範を表している。
【００８７】
【数２０】
2020071766.tif 000022
【００８８】
【数２１】
2020071766.tif 000023
【００８９】
但し、〜ｙm （ｔ）は、下記の式（２２）のように設計される。
【００９０】
【数２２】
2020071766.tif 000024
【００９１】
式（２２）において、Ｇm （１）＝１＋ｐ1 ＋ｐ2 である（式（１０）参照）。
【００９２】
また、式（１８）の右辺第２項のそれぞれの偏微分は、式（２３）のように展開される。
【００９３】
【数２３】
2020071766.tif 000025
【００９４】
式（２３）において、Γ（ｔ）は、下記の式（２４）のように表される。
【００９５】
【数２４】
2020071766.tif 000026
【００９６】
式（２３）に擬似参照入力〜ｒ（ｔ）が含まれていることから、式（１８）は、ＦＲＩＴによるオフライン学習となっている。式（１８）により得られたＫnew （ｔ）を用いて、データベース上の各近傍データの更新を行う。この手順を、式（２０）で表される評価規範が十分小さくなるまで繰り返すことによって、最適なデータベースを取得することができる。システムに対してデータ駆動型制御を適用する際は、各ステップ（時刻）ごとに、前述の「初期データベースの作成」、「距離の算出、近傍の選択」、「局所コントローラの構成」の手順に従って局所コントローラを構成することにより、非線形システムに対して、より有効な制御性能を得ることが可能となる。
【００９７】
＜数値例＞
以下、本実施例による感性フィードバック制御の数値例について説明する。
【００９８】
ここで、図５に示す対象機器（操作者が操作する機器）１０は、油圧ショベルであるとし、下記の式（２５）の一次遅れ系で表す。
【００９９】
【数２５】
2020071766.tif 000027
【０１００】
一方、快適度ｙ（ｔ）は、ウェーバー・フェヒナーの法則（「I.P.Herman、Physics of the Human Body:Biological and Medical Physics、Biomedical Engineering、Springer-Verlag GmbH & CO.KG (2007)」参照）を用いて、快適度ｙ（ｔ）の最大値が１となるように、下記の式（２６）、式（２７）で表す。
【０１０１】
【数２６】
2020071766.tif 000028
【０１０２】
【数２７】
2020071766.tif 000029
【０１０３】
式（２６）、式（２７）において、ｗh （ｔ）は、操作者が脳内で持つ油圧ショベルの目標速度であるが、これは制御系にとって未知である。また、ｅh （ｔ）は、操作者の脳内で感じる油圧ショベルの速度誤差である。
【０１０４】
式（２６）から、操作者の脳内での速度誤差ｅh （ｔ）が完全にゼロであれば、快適度ｙ（ｔ）は最大の１となる。尚、Ｅ（ｔ）は、快適度ｙ（ｔ）に関する変数であり、操作者によって異なる値を持つ。図７は、操作者の脳内での速度誤差ｅh （ｔ）と、快適度ｙ（ｔ）との関係を示す。図７に示すように、速度誤差ｅh （ｔ）が大きくなればなるほど、快適度ｙ（ｔ）が低下している。また、Ｅ（ｔ）が大きくなればなるほど、快適度ｙ（ｔ）の低下率が大きくなる。
【０１０５】
尚、以下に説明する数値例における各設定パラメータは、ｒ＝０．８、ｗh ＝４０、σ＝１０、α＝０、η＝[８０、６０、８０]とした。
【０１０６】
また、内側ループの第２制御器４０には、既存の制御パラメータを用い、外側ループの第１制御部３０の制御パラメータの調整には、非線形システムに対して有効なデータベース駆動型制御（「脇谷伸  他、ＦＲＩＴ法を用いた非線形ＰＩＤ制御系の設計、計測と制御、Vol.52、No.10、pp.885-891 (2013)」参照）を用いた。
【０１０７】
また、初期データ｛ｕ0 、ｙ0 ｝を取得するための外側ループ（第１制御部３０）のＰＩＤゲインは、ＫP ＝３．５、ＫI ＝０．５、ＫD ＝３．５とし、内側ループ（第２制御器４０）のＰＩＤゲインは、ＫP ’＝１．５、ＫI ’＝０．１、ＫD ’＝０．１とした。
【０１０８】
以上の設定においてデータベース駆動型制御により第１制御部３０の制御パラメータの調整を行いながら感性フィードバック制御を行った結果（各ステップ（時刻）における快適度ｙ（ｔ）、速度ｕ（ｔ）、トルクｖ（ｔ））を図８に、また、各ステップ（時刻）における調整された制御パラメータ（ＰＩＤゲイン）を図９にそれぞれ示す。
【０１０９】
図８に示すように、ＰＩＤゲインを固定した初期データｙ0 は、式（２６）で表される快適度ｙ（ｔ）が非線形系であるため、目標快適度ｒ（＝０．８）に追従していない。一方、データベース駆動型制御によりＰＩＤゲインの調整を行った場合（本実施例）の快適度ｙは、目標快適度ｒに追従している。すなわち、本実施例によって感性フィードバック制御が実現されている。これは、図９に示すように、ＰＩＤゲインが適応的に調整されているためであると考えられる。
【０１１０】
また、操作者の脳内の目標速度ｗh （＝４０）は制御系にとっては未知であるが、本実施例によって得られた目標速度ｗは、最終的には４０になっている。このことから、データベースを解析すれば、操作者の脳内の目標速度を推定することが可能となることが示唆される。
【０１１１】
また、図８に示すように、負荷が大きいときには、それに伴って大きいトルクｖ（対象機器１０に対する入力）が自動算出されていることが分かる。
【０１１２】
（比較例）
以下、感性フィードバック制御を外側ループのみにより実現した場合、つまり、対象機器に関する内側の制御ループが無い場合を比較例として、本発明の感性フィードバック制御におけるカスケード制御系（外側ループ＋内側ループ）の有効性について説明する。
【０１１３】
図１０は、比較例に係る感性フィードバック制御装置の構成図である。図１０に示すように、比較例に係る感性フィードバック制御装置２００は、対象機器１０と、感性メータ２０と、制御部３０Ａとを備えている。対象機器１０は、前述の実施例と同様に、操作者によって操作される油圧ショベルである。感性メータ２０は、図１に示す本実施形態の感性フィードバック制御装置２００の感性メータ２０と同様に、対象機器１０の出力ｕ（ｔ）に対応する操作者の生体情報ｘ（ｔ）を検出し、当該生体情報ｘ（ｔ）に基づき、操作者の快適度ｙ（ｔ）を決定する。制御部３０Ａは、快適度ｙ（ｔ）に関する目標値（目標快適度）ｒ（ｔ）と、快適度ｙ（ｔ）との差分ｅ（ｔ）に基づき、対象機器１０に対する制御入力ｖ（ｔ）を決定する。
【０１１４】
制御部３０Ａの入出力関係については、例えば、ＫP ’’、ＫI ’’、ＫD ’’をそれぞれＰＩＤパラメータとして、下記の式（２８）に示すＰＩＤ制御がなされる。
【０１１５】
【数２８】
2020071766.tif 000030
【０１１６】
また、制御部３０Ａは、図１０に示すように、データベース６０Ａを利用したデータベース駆動型制御を行う。データベース６０Ａには、制御部３０ＡのＰＩＤパラメータの調整に必要なデータ、例えば、目標快適度ｒ（ｔ）、快適度ｙ（ｔ）及び制御入力ｖ（ｔ）等が逐次蓄積される。
【０１１７】
尚、本比較例でも、対象機器１０の出力ｕ（ｔ）は油圧ショベルのバケットの応答速度、対象機器１０に対する制御入力ｖc （ｔ）はトルクである。
【０１１８】
データベース駆動型制御により制御部３０Ａの制御パラメータの調整を行いながら感性フィードバック制御を行った結果（各ステップ（時刻）における快適度ｙ（ｔ）、速度ｕ（ｔ）、トルクｖ（ｔ））を図１１に示す。尚、初期データ｛ｕ0 、ｙ0 ｝を取得するための制御部３０ＡのＰＩＤゲインは、前述の実施例と同様に、ＫP ＝３．５、ＫI ＝０．５、ＫD ＝３．５とし、その他のパラメータも、前述の実施例と同様に、ｒ＝０．８、ｗh ＝４０、σ＝１０、α＝０、η＝[８０、６０、８０]とした。
【０１１９】
図１１に示すように、比較例においては、データベース駆動型制御によりＰＩＤゲインの調整を行っても、快適度ｙは目標快適度ｒに追従していない。すなわち、比較例に係る感性フィードバック制御は、十分な制御性能を実現できていない。
【０１２０】
この結果から、本発明の感性フィードバック制御においてはカスケード制御系により内側ループ（第２制御部４０）の即応性を向上させていることによって、制御系の設計が容易になり、感性フィードバック制御の精度が向上したものと推認される。
【０１２１】
以上、本発明についての実施形態（実施例を含む）を説明したが、本発明は前述の実施形態のみに限定されず、発明の範囲内で種々の変更が可能である。すなわち、前述の実施形態の説明は、本質的に例示に過ぎず、本発明、その適用物或いはその用途を制限することを意図するものではない。
【０１２２】
例えば、本実施形態においては、油圧ショベル（建設機械）を例として、感性フィードバック制御装置を説明したが、操作者により操作される他の機器についても同様に適用可能である。
【符号の説明】
【０１２３】
１０    対象機器
２０    感性メータ
３０    第１制御部
４０    第２制御部
５０    重み付け（δ）設定部
６０    データベース
１００    感性フィードバック制御装置
【図１】
2020071766.tif 000031
【図２】
2020071766.tif 000032
【図３】
2020071766.tif 000033
【図４】
2020071766.tif 000034
【図５】
2020071766.tif 000035
【図６】
2020071766.tif 000036
【図７】
2020071766.tif 000037
【図８】
2020071766.tif 000038
【図９】
2020071766.tif 000039
【図１０】
2020071766.tif 000040
【図１１】
2020071766.tif 000041