JP7276110B2 - パワートレーンシステム - Google Patents

Description

この発明は、パワートレーンシステムに関し、より詳細には、外部充電及び外部放電のうちの少なくとも一方を実行可能なハイブリッド車両に搭載され、当該ハイブリッド車両を駆動可能な内燃機関及び電動機と、当該電動機に電力を供給するバッテリとを備えるパワートレーンシステムに関する。

例えば、特許文献１には、ハイブリッド車両が開示されている。このハイブリッド車両は、内燃機関と、電動機と、電動機に電力を供給するバッテリとを備えている。そして、ハイブリッド車両は、内燃機関及び電動機の双方によって駆動される。このように構成されたハイブリッド車両では、出発地から目的地までの全体の経路上で燃料消費量を最小とするためのバッテリの充放電計画が作成される。そして、この充放電計画に基づく最適なハイブリッド車両の走行モードが、エンジン走行モード、電動機走行モード、及び、内燃機関と電動機の双方を使用するアシスト走行モードの中から選択される。また、充放電計画に基づくバッテリの充放電制御も実行される。

特開２０１５－１５５２６１号公報

ところで、ハイブリッド車両を駆動可能な内燃機関及び電動機と、当該電動機に電力を供給するバッテリとを備えるパワートレーンシステムにおいて、パワートレーンシステムの動作点（以下、「システム動作点」と略する）及びバッテリの充放電量を次のように決定することが考えられる。

具体的には、任意の制御期間（例えば、現時刻からトリップ終了までの期間）内で用いられる車速と車両駆動トルクに対して内燃機関の燃料消費量を最小にするために、この制御期間内の将来の各時点でのシステム動作点及び充放電量を特定する最適な制御入力値を先行的に探索して算出することが考えられる。このように将来の各時点の最適な制御入力値を先行的に算出するためには、将来の車速及び車両駆動トルクの予測が必要になる。しかしながら、上記制御期間中の車速及び車両駆動トルクは必ずしも予測通りになるとは限らず、また、何らかの理由によって予測を行えない場合には制御入力値の探索を適切に行えなくなってしまう。

また、外部機器との間で充放電に構成されたハイブリッド車両に対してパワートレーンシステムが搭載される場合には、前回の制御期間と今回の制御期間との間（例えば、前回のトリップと今回のトリップとの間）の制御停止期間中に、外部機器との間でバッテリの充電又は放電が実行され得る。したがって、このようなハイブリッド車両への搭載が想定されるパワートレーンシステムでは、制御停止期間中にバッテリの充電又は放電が実行され得ることをも想定して制御入力値を適切に決定できることが望まれる。

本発明は、上述のような課題に鑑みてなされたものであり、将来の車速及び車両駆動トルクの予測を必要とすることなく、パワートレーン制御の制御期間内に用いられる車速と車両駆動トルクに対して内燃機関の燃料消費量を最小にする制御入力値を決定でき、かつ、制御停止期間中に外部機器との間で充電又は放電がなされた場合であっても当該制御入力値を適切に決定できるようにしたパワートレーンシステムを提供することを目的とする。

本発明に係るパワートレーンシステムは、
ハイブリッド車両を駆動する内燃機関と、
前記ハイブリッド車両を駆動する電動機と、
前記電動機に電力を供給するバッテリと、
前記内燃機関及び前記電動機を制御するパワートレーン制御を行う制御装置と、
を備えるパワートレーンシステムであって、
前記ハイブリッド車両は、１又は複数の外部機器との間で前記バッテリの充電及び放電のうちの少なくとも一方を実行可能に構成されており、
前記制御装置は、前記パワートレーン制御が行われる制御期間中に
前記バッテリの充電率ＳＯＣの動的挙動を制約条件として、前記ハイブリッド車両の速度及び駆動トルクに対して前記制御期間中の前記内燃機関の燃料消費量ｍ_ｆを最小にする最適化問題を解くことにより、前記バッテリの充放電量及び前記パワートレーンシステムのシステム動作点を特定する制御入力値を求める制御入力決定処理と、
前記制御入力決定処理により求めた前記制御入力値を前記パワートレーンシステムに与えることにより、前記内燃機関と前記電動機とを制御するシステム制御処理と、
を実行し、
前記制御入力決定処理は、
前記最適化問題の随伴変数ｐを時間ステップ毎に更新する随伴変数決定処理と、
前記随伴変数決定処理により決定される前記随伴変数ｐを用いて、以下の式により定義されるハミルトニアンＨを最小とする前記制御入力値を前記時間ステップ毎に探索して算出する制御入力算出処理と、を含み、
前記随伴変数決定処理は、前記随伴変数ｐの初期値の基本値と外部充放電補正値との和を前記初期値とする初期値決定処理を含み、
前記基本値は、前回の前記制御期間中の前記随伴変数ｐの最終値又は平均値であり、
前記外部充放電補正値は、今回の前記制御期間の開始時の前記充電率ＳＯＣから前回の前記制御期間の終了時の前記充電率ＳＯＣを引いて得られる外部充放電量に基づいて決定される。

ただし、上記の式の右辺第１項は、前記燃料消費量ｍ_ｆの時間変化率ｄｍ_ｆ／ｄｔであり、右辺第２項は、前記随伴変数ｐと前記充電率ＳＯＣの時間変化率ｄＳＯＣ／ｄｔとの積である。

前記初期値決定処理は、前記外部充放電量が正の場合には、前記外部充放電量がゼロの場合と比べて、前記初期値を大きくしてもよい。

前記初期値決定処理は、前記外部充放電量が負の場合には、前記外部充放電量がゼロの場合と比べて、前記初期値を小さくしてもよい。

前記初期値決定処理による前記外部充放電補正値は、前記外部充放電量と第１ゲインとの積であってもよい。そして、前記第１ゲインは、前記ハイブリッド車両の走行履歴に基づいて設定されてもよい。

前記走行履歴は、前回の前記制御期間中の前記速度及び前記駆動トルクのそれぞれの平均値である第１平均車速及び第１平均駆動トルクのうちの少なくとも一方を含んでもよい。

前記走行履歴が前記第１平均車速を含む場合において、前記第１ゲインは、前記第１平均車速が高いほど大きくてもよい。

前記走行履歴が前記第１平均駆動トルクを含む場合において、前記第１ゲインは、前記第１平均駆動トルクが高いほど大きくてもよい。

前記随伴変数決定処理は、前記充電率ＳＯＣの目標充電率と実充電率との差に基づいて前記随伴変数ｐを修正する随伴変数修正処理を含んでもよい。

前記随伴変数修正処理は、前記目標充電率よりも前記実充電率が低い場合には、前記随伴変数ｐを小さくしてもよい。

前記随伴変数修正処理は、前記目標充電率よりも前記実充電率が高い場合には、前記随伴変数ｐを大きくしてもよい。

前記随伴変数修正処理による前記随伴変数ｐの補正値は、前記差と第２ゲインとの積であってもよい。そして、前記第２ゲインは、前記ハイブリッド車両の走行履歴に基づいて設定されてもよい。

前記走行履歴は、第１所定期間における前記速度の平均値である第２平均車速、及び第２所定期間における前記駆動トルクの平均値である第２平均駆動トルクのうちの少なくとも一方を含んでもよい。

前記走行履歴が前記第２平均車速を含む場合において、前記第２ゲインは、前記第２平均車速が高いほど大きくてもよい。

前記走行履歴が前記第２平均駆動トルクを含む場合において、前記第２ゲインは、前記第２平均駆動トルクが高いほど大きくてもよい。

本発明に係るパワートレーンシステムによれば、上記の式により定義されるハミルトニアンＨが最小となる制御入力値が算出される。したがって、燃料消費量が最小となる最適な制御入力値を時間ステップ毎に決定できるようになる。そして、当該制御入力値は、バッテリの充放電量及びパワートレーンシステムのシステム動作点を特定するものである。このため、決定された最適な制御入力値に応じて、充放電量とシステム動作点とを適切に決定できるようになる。

そして、本発明に係る制御入力決定処理によれば、対象とする最適化問題の随伴変数の初期値の基本値として、前回の制御期間中の随伴変数の最終値又は平均値が用いられる。これにより、将来（制御入力値を最適化したい期間）の車速及び駆動トルクが未知であっても（換言すると、将来の車速及び駆動トルクを別途予測する必要なしに）、充放電量及びシステム動作点を時間ステップ毎に適切に決定できるようになる。そのうえで、初期値決定処理によれば、制御停止期間中の外部充放電量に基づく外部充放電補正値によって、上記基本値が補正される。これにより、制御停止期間中に外部機器との間で充電又は放電がなされた場合であっても、初期値決定処理を備えない例と比べて、充放電量及びシステム動作点を特定する制御入力値をより適切に決定できるようになる。

本発明の実施の形態１に係るパワートレーンシステムの構成例を説明するための模式図である。 本発明の実施の形態１に係るパワートレーン制御に関連して図２に示す制御装置が実行する処理を示すブロック図である。 ハミルトニアンＨを最小とする制御入力値ｕ_ｋ（最適制御入力値ｕ_ｋ ^＊）の探索の概念図である。 最適制御入力値ｕ^＊及びその算出に用いる随伴変数ｐが時間ステップΔｔ毎に繰り返し算出される様子を概念的に表した図である。 本発明の実施の形態１に係るパワートレーン制御に関する処理のルーチンを示すフローチャートである。 制御停止期間ｔｓ中の外部充電又は外部放電の実行に伴う課題について説明するためのタイムチャートである。 本発明の実施の形態１に係る初期値決定処理による補正の効果を表したタイムチャートである。 ＳＯＣ（ｋ＋１）と随伴変数ｐ（ｋ）との関係を表したグラフである。 随伴変数ｐとＳＯＣとの関係を車両走行モード毎に分けて表したグラフである。 ＳＯＣに対する随伴変数ｐの感度、すなわち、比ΔＰ／ΔＳＯＣをモードＡ～モードＧの間で比較して表したグラフである。 本発明の実施の形態２に係るパワートレーン制御に関連して実施の形態２の制御装置が実行する処理を示すブロック図である。 第１平均車速ａｖｅＶ１及び第１平均駆動トルクａｖｅＴｐ１と、比ΔＰ／ΔＳＯＣとの関係を車両走行モード毎に表したグラフである。 本発明の実施の形態２に係る初期値決定処理のサブルーチンを示すフローチャートである。 本発明の実施の形態２に係る初期値決定処理による補正の効果を表したタイムチャートである。 実施の形態１、２に係るパワートレーン制御の課題と、その対策を有する実施の形態３に係るパワートレーン制御とを比較しながら説明するためのタイムチャートである。 本発明の実施の形態３に係るパワートレーン制御に関連して実施の形態３の制御装置が実行する処理を示すブロック図である。 本発明の実施の形態３に係るパワートレーン制御に関する処理のルーチンを示すフローチャートである。 本発明の実施の形態４に係るパワートレーン制御に関連して実施の形態４の制御装置が実行する処理を示すブロック図である。 本発明の実施の形態４に係るパワートレーン制御において、偏差ΔＳＯＣ（ｋ）に応じて随伴変数ｐ（ｋ）を修正するために実行される処理のサブルーチンを示すフローチャートである。 本発明の実施の形態４の随伴変数修正処理との対比のために参照する随伴変数修正処理（比較例１、２）を利用する場合の動作を説明するためのタイムチャートである。 本発明の実施の形態４の随伴変数修正処理を利用する場合の動作を説明するためのタイムチャートである。

以下、添付図面を参照して、本発明の実施の形態について説明する。ただし、各図において共通する要素には、同一の符号を付して重複する説明を省略又は簡略する。以下に示す実施の形態において各要素の個数、数量、量、範囲等の数に言及した場合、特に明示した場合や原理的に明らかにその数に特定される場合を除いて、その言及した数に、この発明が限定されるものではない。また、以下に示す実施の形態において説明する構造やステップ等は、特に明示した場合や明らかに原理的にそれに特定される場合を除いて、この発明に必ずしも必須のものではない。

１．実施の形態１
まず、図１～図８を参照して、本発明の実施の形態１について説明する。

１－１．パワートレーンシステムの構成例
図１は、本発明の実施の形態１に係るパワートレーンシステム１０の構成例を説明するための模式図である。図１に示すパワートレーンシステム１０は、車両の動力源として、内燃機関２０とともに第１モータジェネレータ３０（以下、「ＭＧ１」と略する）及び第２モータジェネレータ３２（以下、「ＭＧ２」と略する）を備えている。すなわち、パワートレーンシステム１０は、一例としてハイブリッド車両に適用されている。

内燃機関２０は、一例として、火花点火式エンジンである。しかしながら、本発明の対象となる内燃機関は、圧縮着火式エンジンであってもよく、また、その気筒数及び気筒配置は特に限定されない。内燃機関２０は、エンジントルクＴｅを制御するためのアクチュエータとして、スロットルバルブ２２と燃料噴射弁２４と点火装置２６とを備えている。スロットルバルブ２２は、吸気通路（図示省略）に配置され、吸入空気流量を制御する。燃料噴射弁２４は、各気筒に配置され、例えば気筒内に直接燃料を噴射する。点火装置２６は、各気筒に配置された点火プラグを用いて、気筒内の混合気に点火する。また、内燃機関２０は、各種エンジン制御に用いられる各種センサを備えている。ここでいう各種センサは、クランク角に応じた信号を出力するクランク角センサ２８を含む。

ＭＧ１及びＭＧ２は、共に発電可能な電動機である。すなわち、ＭＧ１及びＭＧ２は、供給された電力によりトルクを出力する電動機としての機能と、入力された機械的動力を電力に変換する発電機としての機能とを兼ね備えており、例えば交流同期型のモータジェネレータである。図１に示すパワートレーンシステム１０では、ＭＧ１は、主に発電機として用いられ、ＭＧ２は、主に車両を駆動する電動機として用いられる。このため、パワートレーンシステム１０の例では、ＭＧ１及びＭＧ２のうち、ＭＧ２が主に本発明に係る「電動機」の一例に相当する。また、ＭＧ１も、運転条件によっては負トルクを発生させてエンジントルクＴｅを支えることにより車両の駆動トルクを発生させる場合があるので、ＭＧ２とともに本発明に係る「電動機」の一例に相当する。

内燃機関２０、ＭＧ１及びＭＧ２は、動力分割機構３４及び減速機構３６を介して車輪３８と連結されている。動力分割機構３４は、例えばプラネタリギヤユニットであり、内燃機関２０から出力されるエンジントルクＴｅをＭＧ１と車輪３８とに分割する。より詳細には、動力分割機構３４において、サンギヤはＭＧ１の出力軸３０ａに連結され、プラネタリキャリアは内燃機関２０のクランクシャフト２０ａに連結され、リングギヤはＭＧ２の出力軸３２ａに連結されている。リングギヤは、プラネタリキャリアに回転自在に支持された複数のピニオンギヤを介してサンギヤと噛み合っている。内燃機関２０から出力されるエンジントルクＴｅ又はＭＧ２から出力されるＭＧ２トルクＴｍは、減速機構３６を介して車輪３８に伝達される。すなわち、内燃機関２０及びＭＧ２は、車両の駆動力を生じさせ、かつ、車両の駆動力を制御する。ＭＧ１は、動力分割機構３４を介して内燃機関２０から供給されたエンジントルクＴｅを用いて回生発電可能である。このため、ＭＧ１も車両の駆動力を制御するために用いられる。また、ＭＧ２は、車両減速時には発電機として機能し、車両運動エネルギを回収して電力に変換する。

ＭＧ１及びＭＧ２は、インバータ４０及び昇圧コンバータ４２を介してバッテリ４４と電力の授受を行う。インバータ４０は、バッテリ４４に蓄えられた電力を直流から交流に変換してＭＧ２に供給するとともに、ＭＧ１及びＭＧ２によって生成される電力を交流から直流に変換してバッテリ４４に蓄える。このため、バッテリ４４は、ＭＧ１及びＭＧ２で生じた電力によって充電され、ＭＧ２で消費される電力により放電される。昇圧コンバータ４２は、バッテリ４４の電圧を必要に応じて昇圧させる。

パワートレーンシステム１０を搭載するハイブリッド車両は、本発明に係る「外部機器」の一例に相当する外部電源４５及び電気製品４６との間で、バッテリ４４の外部充電及び外部放電をそれぞれ実行可能に構成されている。換言すると、本ハイブリッド車両は、外部給電機能を有するプラグインハイブリッド車両である。具体的には、ハイブリッド車両は、充電器４７とインレット４８とを備えている。バッテリ４４には、外部電源（充電設備）４５からインレット４８及び充電器４７を介して電力が供給される。その結果、バッテリ４４が外部電源４５によって充電される。一方、電気製品４６は、例えば家電製品である。バッテリ４４の電力は、電気製品４６を作動させることによって消費される。すなわち、バッテリ４４が電気製品４６によって放電される。また、外部電源４５からバッテリ４４に供給された電力は、本ハイブリッド車両の駆動にも使用される。なお、このような例に限らず、本発明に係る「パワートレーンシステム」を搭載するハイブリッド車両は、１又は複数の外部機器との間でバッテリの充電及び放電のうちの何れか一方を実行可能に構成されていてもよい。

本実施形態のパワートレーンシステム１０は、さらに、パワートレーン（内燃機関２０、ＭＧ１及びＭＧ２）を制御するための制御装置５０を備えている。制御装置５０は、プロセッサ５０ａとメモリ５０ｂとを有する電子制御ユニット（ＥＣＵ）である。メモリ５０ｂは、パワートレーンシステム１０を制御するためのプログラムを記憶している。プロセッサ５０ａは、メモリ５０ｂからプログラムを読み出して実行する。制御装置５０は、パワートレーンを制御するための各種センサからセンサ信号を取り込む。また、プロセッサ５０ａは、取り込まれたセンサ信号を用いて各種プログラムを実行し、パワートレーンの各種アクチュエータを操作するための操作信号を出力する。

制御装置５０には、上述したクランク角センサ２８等のエンジン制御のための各種センサに加え、アクセルポジションセンサ、ブレーキポジションセンサ及び車速センサを含むセンサ類５２等、パワートレーンの制御に用いられる各種センサが電気的に接続されている。制御装置５０は、クランク角センサ２８からの信号を用いてエンジン回転数Ｎｅを算出できる。

また、制御装置５０には、上述した内燃機関２０（スロットルバルブ２２、燃料噴射弁２４及び点火装置２６）、ＭＧ１及びＭＧ２を含むパワートレーンを制御するための各種アクチュエータが電気的に接続されている。さらに、制御装置５０は、バッテリ４４の充電率（ＳＯＣ）を算出可能に構成されている。より詳細には、ＳＯＣは、バッテリ４４の満充電量に対する現在の充電量の割合（％）に相当する。ＳＯＣの取得手法は、特に限定されないが、ＳＯＣは、例えば、バッテリ４４に出入りする電流を計測する電流センサを備え、計測された電流を時間で積分することにより取得（算出）できる。なお、制御装置５０は、複数のＥＣＵから構成されていてもよい。

上述の動力分割機構３４を利用する本パワートレーンシステム１０では、車両の要求駆動トルクをＴｐと称し、ＭＧ１のトルク及び回転数をＴｇ及びＮｇと称し、ＭＧ２のトルク及び回転数をＴｍ及びＮｍと称すると、パワートレーンシステム１０による車両の制御のための主要なパラメータの関係は、次の（１）～（４）式のように表すことができる。なお、（１）～（４）式中のｇ１～ｇ６は何れも正の値である。より詳細には、ｇ１～ｇ３及びｇ６は動力分割機構３４のギヤ比であり、ｇ４及びｇ５は動力分割機構３４のギヤ比だけでなく、車速Ｖをリングギヤの回転数に変換するためのデフ比及びタイヤ径をも含めて決定される値である。

上記の（１）～（４）式の関係によれば、ある車速Ｖ及び要求駆動トルク（車両駆動トルクの要求値）Ｔｐの下では、エンジントルクＴｅとエンジン回転数Ｎｅとが定まれば、現在のバッテリ４４の充放電量Ｐｃｈｇと現在のパワートレーンシステム１０の動作点（以下、「システム動作点」と略する）とが特定される。より詳細には、システム動作点が特定されることは、内燃機関２０のトルクＴｅ及び回転数Ｎｅ、ＭＧ１のトルクＴｇ及び回転数Ｎｇ、並びにＭＧ２のトルクＴｍ及び回転数Ｎｍが特定されることと同じである。また、バッテリ４４の充放電量Ｐｃｈｇは、ＭＧ１のトルクＴｇ及び回転数Ｎｇ、並びにＭＧ２のトルクＴｍ及び回転数Ｎｍが特定されることによって特定される。これは、充放電量Ｐｃｈｇは、ＭＧ１のトルクＴｇと回転数Ｎｇとの積と、ＭＧ２のトルクＴｍと回転数Ｎｍとの積に応じた値となるためである。そして、バッテリ４４の充放電量Ｐｃｈｇが特定されると、ＳＯＣを特定できる。

付け加えると、（４）式の関係により、エンジン回転数Ｎｅに代えてＭＧ１回転数ＮｇがエンジントルクＴｅとともに定められた場合にも、充放電量とシステム動作点とを同様に特定できる。また、（２）式の関係により、充放電量Ｐｃｈｇとシステム動作点の特定のためにエンジントルクＴｅに代えてＭＧ１トルクＴｇが用いられてもよい。さらに、（３）式の関係により、車速Ｖに代えてＭＧ２回転数Ｎｍが要求駆動トルクＴｐとともに与えられた場合にも、充放電量Ｐｃｈｇとシステム動作点とを上記と同様に特定できる。なお、車両駆動トルクは、車両の走行負荷に相当する。

１－２．パワートレーン制御
制御装置５０は、車両走行のために内燃機関２０、ＭＧ１及びＭＧ２を制御するパワートレーン制御を実行する。図２は、本発明の実施の形態１に係るパワートレーン制御に関連して制御装置５０が実行する処理を示すブロック図である。制御装置５０がパワートレーン制御に関連して実行する処理は、「制御入力決定処理５４」と「システム制御処理５６」とを含む。これらの処理は、パワートレーン制御が行われる制御期間τ中に行われる。この制御期間τの一例は、車両の１トリップ（システム起動開始時からシステム停止時までの期間）である。ただし、パワートレーン制御が行われる制御期間τは、必ずしも１トリップに限られず、１トリップ中の任意の一部の期間であってもよい。

１－２－１．制御入力決定処理５４
バッテリ４４のＳＯＣは、パワートレーンシステム１０の作動中（車両走行中）に、バッテリ４４の充放電量Ｐｃｈｇが変化することによって変動する。制御入力決定処理５４は、ＳＯＣの動的挙動（ダイナミクス）を制約条件として用いて、車速Ｖ及び要求駆動トルクＴｐに対して制御期間τ中の内燃機関２０の燃料消費量ｍ_ｆを最小にするという最適化問題（最小化問題）を解くことにより、制御入力値ｕを求めるために実行される。バッテリ４４の充放電量Ｐｃｈｇ及びシステム動作点を特定可能なパラメータであれば、当該制御入力値ｕとして用いることができる。図１に示す構成を有するパワートレーンシステム１０の例では、エンジントルクＴｅとエンジン回転数Ｎｅが制御入力値ｕとなる。換言すると、制御入力決定処理５４は、燃料消費量ｍ_ｆを最小とするための制御入力値（以下、「最適制御入力値ｕ^*」と称する）を算出するために実行される。最適制御入力値ｕ^*の算出は、所定の時間ステップΔｔ毎に実行される。その結果、燃料消費量ｍ_ｆを最小にできるバッテリ４４の充放電量Ｐｃｈｇとシステム動作点とが時間ステップΔｔ毎に取得される。

この最適化問題では、以下の（５）式により定義されるハミルトニアンＨが用いられる。（５）式に示されるように、ハミルトニアンＨは、燃料消費量ｍ_ｆの時間変化率である燃料消費率ｄｍ_ｆ／ｄｔと、本最適化問題の随伴変数ｐとＳＯＣの時間変化率ｄＳＯＣ／ｄｔの積との和である。燃料消費率ｄｍ_ｆ／ｄｔは、後述の（９）式中に表されているように、エンジントルクＴｅとエンジン回転数Ｎｅの関数として表すことができる。ＳＯＣの時間変化率ｄＳＯＣ／ｄｔは、ＳＯＣの動的挙動（ダイナミクス）を表す以下の状態方程式（（６）式）のように、ＳＯＣ、エンジン回転数Ｎｅ、エンジントルクＴｅ、車速Ｖ及び要求駆動トルクＴｐの関数ｆによって表される。付け加えると、このように定義されたハミルトニアンＨを最小にすることは、随伴変数ｐによって決定される燃料消費とバッテリ４４の電力消費とのバランスの下で、これらの燃料消費と電力消費についての合計のエネルギ消費を最小にすることを意味する。

以下の（７）式は、（５）式中の随伴変数ｐが満たすべき関係を表している。（７）式中のＮｅ^＊及びＴｅ^＊は、最適制御入力値ｕ^＊に該当する時のエンジン回転数Ｎｅ及びエンジントルクＴｅである。また、（８）式は、本最適化問題の対象となる制御入力値ｕに該当するエンジン回転数Ｎｅ及びエンジントルクＴｅの値の範囲（集合Ｕ）を示している。（８）式に示されるように、内燃機関２０を停止してＭＧ２のみを用いて車両を走行するＥＶモードでは、エンジン回転数Ｎｅ及びエンジントルクＴｅは共にゼロとなる。また、集合Ｕに含まれるエンジン回転数Ｎｅ及びエンジントルクＴｅの値には、内燃機関２０が作動している車両走行モードであるＨＶモード（engine on range）において取り得る値も含まれる。具体的には、集合Ｕには、内燃機関２０のハードウェア上の制約、並びにＭＧ１回転数Ｎｇ及び動力分割機構３４のピニオンギヤの回転数の制約の下で取り得る値も含まれる。したがって、（８）式は、上記（６）式とともに、本最適化問題の制約条件に該当する。

（随伴変数決定処理５８と制御入力算出処理６０）
制御入力決定処理５４は、図２に示すように、「随伴変数決定処理５８」と「制御入力算出処理６０」とを含む。制御入力算出処理６０は、（５）式により定義されるハミルトニアンＨを最小とする制御入力値ｕ（ｔ）を時間ステップΔｔ毎に探索して算出する。本最適化問題は、ハミルトニアンＨを最小とする制御入力値ｕ（ｔ）を求めるという問題に置き換えて考えることができ、この問題は、以下の（９）式のように表される。すなわち、（９）式の右辺に表されているハミルトニアンＨは、この問題の目的関数に相当する。各サンプル時刻ｔにおいてハミルトニアンＨを最小とする制御入力値ｕ（ｔ）は、例えば準ニュートン法を利用して（９）式により表される問題を解くことにより算出できる。

付け加えると、パワートレーンシステム１０の例では、制御入力算出処理６０は、ＨＶモードにおけるハミルトニアンＨ_ＨＶ（以下の（１０）式参照）を最小とする最適制御入力値ｕ^＊ＨＶ（ｔ）を時間ステップΔｔ毎に探索して算出するために用いられる。そのうえで、パワートレーンシステム１０は、ＨＶモードとともにＥＶモードを利用可能に構成されている。このように構成されたパワートレーンシステム１０において、制御入力決定処理５４は、次のような手法を用いて、ハミルトニアンＨを最小とする（すなわち、燃料消費量ｍ_ｆを最小とする）最適制御入力値ｕ^＊を決定する。すなわち、制御入力決定処理５４によれば、最適制御入力値ｕ^＊は、制御入力算出処理６０によって算出されるＨＶモードにおける最適制御入力値ｕ^＊ＨＶと、ＥＶモードにおける制御入力値ｕ^ＥＶ（Ｎｅ及びＴｅが共にゼロ）の中から決定される。

より詳細には、ハミルトニアンＨ_ＨＶがＥＶモード時のハミルトニアンＨ_ＥＶ（以下の（１１）式参照）よりも小さい場合には、最適制御入力値ｕ^＊ＨＶが最適制御入力値ｕ^＊として扱われる。したがって、このような最適制御入力値ｕ^＊がＥＶモードの使用中にパワートレーンシステム１０に与えられた場合には、内燃機関２０が始動される。一方、ハミルトニアンＨ_ＥＶの方がハミルトニアンＨ_ＨＶよりも小さい場合には、制御入力値ｕ^ＥＶが最適制御入力値ｕ^＊として扱われる。この場合には、最適制御入力値ｕ^＊に相当する最適エンジン回転数Ｎｅ^＊及び最適エンジントルクＴｅ^＊は共にゼロとなる。したがって、このような最適制御入力値ｕ^＊がＨＶモードの使用中にパワートレーンシステム１０に与えられた場合には、内燃機関２０が停止される。なお、ＥＶモードを利用しないハイブリッド車両のパワートレーンシステムであれば、上述のＨＶモードにおける最適制御入力値ｕ^＊ＨＶをそのまま最適制御入力値ｕ^＊として利用すればよい。

サンプル時刻ｔにおけるＨＶモード時のハミルトニアンＨ_ＨＶ及びＥＶモード時のハミルトニアンＨ_ＥＶは、それぞれ、以下の（１０）及び（１１）式のように表される。なお、（１０）式中の集合Ｕｅは（９）式中のＮｅ及びＴｅの範囲（Ｎｅ、Ｔｅ in engine on range）と同じである。なお、ＨＶモード時のハミルトニアンＨ_ＨＶの算出の際、（１０）式の右辺の燃料消費率ｄｍｆ／ｄｔの関数及びＳＯＣの関数ｆに代入されるエンジン回転数Ｎｅ及びエンジントルクＴｅの値は、ＨＶモードにおける最適制御入力値ｕ^＊ＨＶに相当するＮｅ^＊ＨＶ及びＴｅ^＊ＨＶである。また、ハミルトニアンＨ_ＨＶ及びＨ_ＥＶの算出の際、ＳＯＣの関数ｆに代入されるＳＯＣ、現在車速Ｖ及び要求駆動トルクＴｐとして、後述のステップＳ１００の処理により取得される値を用いることができる。そして、関数ｆに掛け合わされる随伴変数ｐ（ｔ）として、後述のステップＳ１０２の処理により取得される値を用いることができる。

以下の説明では、連続的な時刻を示す場合は「ｔ」が用いられ、離散的な時刻（すなわち、時間ステップΔｔ毎の各時刻）を示す場合は「ｋ」が用いられる。制御入力算出処理６０による制御入力値ｕ（ｔ）の算出のために、現在のＳＯＣ、現在の車速Ｖ及び要求駆動トルクＴｐとともに、随伴変数決定処理５８によって決定される随伴変数ｐ（ｔ）が与えられる。具体的には、随伴変数決定処理５８は、後述の初期値決定処理５９を含む。随伴変数ｐの初期値ｐ_０（すなわち、サンプル時刻ｋ＝０の時の随伴変数ｐ（０））は、初期値決定処理５９によって決定される。

制御期間τの開始後の２回目以降の随伴変数ｐの演算に関しては、随伴変数決定処理５８によれば、時間ステップΔｔ毎に随伴変数ｐ（ｔ）が逐次算出されて更新される。より詳細には、この随伴変数ｐ（ｔ）の更新（すなわち、次回のサンプル時刻（ｋ＋１）の随伴変数ｐ（ｋ＋１）の算出）は、（７）式に対して現在のサンプル時刻ｋのＳＯＣ（ｋ）、最適制御入力値ｕ^＊（ｋ）（＝エンジン回転数Ｎｅ及びエンジントルクＴｅの最適値Ｎｅ^＊（ｋ）及びＴｅ^＊（ｋ））、及び随伴変数ｐ（ｋ）を代入しつつ、（７）式を数学的に解くことにより行われる。なお、最適制御入力値ｕ^＊（ｋ）（＝最適値Ｎｅ^＊（ｋ）及びＴｅ^＊（ｋ））の算出は、ＨＶモードだけでなくＥＶモードをも考慮して行われる。

図３は、ハミルトニアンＨ（より詳細には、Ｈ_ＨＶ）を最小とする制御入力値ｕ_ｋ（最適制御入力値ｕ_ｋ ^＊ＨＶ）の探索の概念図である。図３には、制御入力値ｕに相当するエンジン回転数ＮｅとエンジントルクＴｅとを座標軸とする平面において、ハミルトニアンＨ_ＨＶの等高線が表されている。そして、図３には、サンプル時刻ｋにおける制御入力値ｕ_ｋ（＝ｕ（ｋ））が、図２に示す制御入力算出処理６０による探索によって、最適制御入力値ｕ_ｋ ^＊ＨＶに近づいていく様子が表されている。より詳細には、図３には、ｉ番目の探索値ｕ_ｋ ^ｉに対して所定量Δｕだけ制御入力値ｕ_ｋを移動させた際に、その次の（ｉ＋１）番目の探索値ｕ_ｋ ^ｉ＋１の方がｉ番目の探索値ｕ_ｋ ^ｉと比べてハミルトニアンＨ_ＨＶが小さくなっている様子が表されている。そして、図３に示す一例では、その次の（ｉ＋２）回目の探索において、ハミルトニアンＨ_ＨＶを最小とする最適制御入力値ｕ_ｋ ^＊ＨＶが得られている。このような探索による最適制御入力値ｕ_ｋ ^＊ＨＶの算出は、例えば、準ニュートン法を利用して行うことができる。

次に、図４は、最適制御入力値ｕ^＊及びその算出に用いる随伴変数ｐが時間ステップΔｔ毎に繰り返し算出される様子を概念的に表した図である。図４に示すように、制御期間τの開始時（例えば、パワートレーンシステム１０の起動時）に相当する初回のサンプル時刻（ｋ＝０）では、制御入力算出処理６０は、最適制御入力値ｕ（０）^＊（＝ｕ_０ ^＊）を算出する。この最適制御入力値ｕ_０ ^＊（すなわち、Ｎｅ_０ ^＊及びＴｅ_０ ^＊）は、随伴変数ｐの初期値ｐ_０及び現在の（実）ＳＯＣ_０と、図示省略された現在車速Ｖ及び要求駆動トルクＴｐとを用いて、（９）式で表される最適化問題を解くことにより算出される。なお、図４を参照して説明する最適制御入力値ｕ（ｔ）^＊の算出についても、ＨＶモードだけでなくＥＶモードをも考慮して行われている。

また、初回のサンプル時刻（ｋ＝０）では、随伴変数決定処理５８は、最適制御入力値ｕ_０ ^＊（Ｎｅ_０ ^＊及びＴｅ_０ ^＊）とＳＯＣ_０と初期値ｐ_０とが代入された（７）式を解く。これにより、次回（すなわち、時間ステップΔｔが１つ後）のサンプル時刻（ｋ＝１）の演算に用いられる随伴変数ｐ_１が算出される。

次のサンプル時刻（ｋ＝１）では、制御入力算出処理６０は、上記のように算出された随伴変数ｐ_１及び現在のＳＯＣ_１と、現在車速Ｖ及び要求駆動トルクＴｐとを用いて、最適制御入力値ｕ_１ ^＊を算出する。以後、同様の演算が時間ステップΔｔ毎に繰り返し実行される。

（初期値決定処理５９）
初期値決定処理５９は、初期値ｐ_０の基本値として、車両の前回の制御期間τ中の随伴変数ｐの最終値ｐ_ｆを用いる。そして、初期値決定処理５９は、以下の（１２）式に示すように、この基本値（最終値ｐ_ｆ）と外部充放電補正値ｐ_ｇとの和を初期値ｐ_０として算出（決定）する。なお、当該基本値として、最終値ｐ_ｆに代え、前回の制御期間τ中の各サンプル時刻ｋにおいて算出された随伴変数ｐの平均値ｐ_ＡＶＥが用いられてもよい。

外部充放電補正値ｐ_ｇは、外部充放電量ΔＳＯＣ_ｇに基づいて決定される。ここで、説明の便宜上、今回の制御期間τの開始時の充電率ＳＯＣを、単に「今回の初期ＳＯＣ_０ＮＥＷ」とも称し、前回の制御期間τの終了時の前記充電率ＳＯＣを、単に「前回の終端ＳＯＣ_ｆＯＬＤ」とも称する。外部充放電量ΔＳＯＣ_ｇは、以下の（１３）式に示すように、今回の初期ＳＯＣ_０ＮＥＷから前回の終端ＳＯＣ_ｆＯＬＤを引くことにより得られる。そして、外部充放電補正値ｐ_ｇは、より詳細には、（１４）式に示すように、外部充放電量ΔＳＯＣ_ｇと補正ゲインＧ１との積によって表される。補正ゲインＧ１は正の値である。また、本実施形態では、補正ゲインＧ１の一例として、所定の固定値が用いられている。なお、補正ゲインＧ１は、本発明に係る「第１ゲイン」の一例に相当する。

前回の制御期間τの終了時から今回の制御期間τの開始時までの期間を、便宜上、「制御停止期間ｔｓ」と称する（後述の図６参照）。外部充放電量ΔＳＯＣ_ｇは、この制御停止期間ｔｓ中に外部充電及び外部放電の何れもなされなかった場合はゼロとなる。このため、制御停止期間ｔｓ中に外部充放電がなされなかった場合には、外部充放電補正値ｐ_ｇはゼロとなり、初期値ｐ_０は補正されずに基本値（最終値ｐ_ｆ）のままとなる。

一方、制御停止期間ｔｓ中に外部充電がなされた場合には、外部充放電量ΔＳＯＣ _ｇ は正となり、外部充放電補正値ｐ_ｇ も正となる。その結果、（１２）式より、初期値ｐ_０は、正の外部充放電補正値ｐ_ｇが基本値に加算されるので、外部充放電補正値ｐ_ｇがゼロの場合と比べて大きくなる。付け加えると、随伴変数ｐは、基本的には負の値をとる。このため、初期値ｐ_０が負であるとした場合、初期値ｐ_０は、当該補正により、絶対値としては小さくなる。より詳細には、（１４）式によれば、外部充放電補正値ｐ_ｇは、補正ゲインＧ１が同じであれば、外部充放電量ΔＳＯＣ_ｇ（正の値）が大きいほど大きくなる。その結果、初期値ｐ_０は、（１２）式より、外部充放電量ΔＳＯＣ_ｇが大きいほど大きくなるように補正される。

また、制御停止期間ｔｓ中に外部放電がなされた場合には、外部充放電量ΔＳＯＣ _ｇ は負となり、外部充放電補正値ｐ_ｇ も負となる。その結果、（１２）式より、初期値ｐ_０は、負の外部充放電補正値ｐ_ｇが基本値に加算されるので、外部充放電補正値ｐ_ｇがゼロの場合と比べて小さくなる。付け加えると、初期値ｐ_０が負であるとした場合、初期値ｐ_０は、当該補正により、絶対値としては大きくなる。より詳細には、（１４）式によれば、外部充放電補正値ｐ_ｇは、補正ゲインＧ１が同じであれば、負の値である外部充放電量ΔＳＯＣ_ｇの絶対値が大きいほど小さくなる（負側で大きくなる）。その結果、初期値ｐ_０は、（１２）式より、負の値である外部充放電量ΔＳＯＣ_ｇの絶対値が大きいほど小さくなるように補正される。

１－２－２．システム制御処理５６
図２に示すシステム制御処理５６は、制御入力決定処理５４によって上述のように求められた最適制御入力値ｕ（ｔ）^＊（すなわち、Ｎｅ^＊及びＴｅ^＊）をパワートレーンシステム１０に与えることにより、内燃機関２０、ＭＧ１及びＭＧ２を制御する。より詳細には、決定された最適制御入力値ｕ（ｔ）^＊（すなわち、Ｎｅ^＊及びＴｅ^＊）に基づいて、サンプル時刻ｋにおけるエンジントルクＴｅ、ＭＧ１トルクＴｇ及びＭＧ２トルクＴｍのそれぞれの目標値が、一例として、次のような手法を利用して決定される。

以下の（１５）式は、パワートレーンシステム１０の駆動系の運動方程式を示している。換言すると、（１５）式は、パワートレーンシステム１０によって制御される状態量（制御量）に相当する要求駆動トルクＴｐ（Ｎｍ）、充放電量Ｐｃｈｇ（Ｗ）及びエンジン回転数Ｎｅの変化率ｄＮｅ／ｄｔ（ｒａｄ／ｓ^２）と、トルクデバイス（内燃機関２０、ＭＧ１、ＭＧ２）の操作量である各トルクＴｅ、Ｔｇ、Ｔｍ（Ｎｍ）との関係を規定する線形の状態方程式である。

（１５）式において、ｃ（ｃ_１１、ｃ_１２、…）は、充放電量ＰｃｈｇとトルクＴｇ、Ｔｍとに関するｃ_２２及びｃ_２３を除き、パワートレーンシステム１０のハードウェア諸元（例えば、各部のイナーシャ及びギヤ比）に応じて定まる定数である。ｃ_２２及びｃ_２３については、運転中の回転数Ｎｇ、Ｎｍの変化に応じて変化する。

制御入力決定処理５４によって最適エンジン回転数Ｎｅ^＊及び最適エンジントルクＴｅ^＊が決定されると、これらに応じたシステム動作点（エンジン回転数Ｎｅを含む）と充放電量Ｐｃｈｇとが特定される。そして、エンジン回転数の変化率ｄＮｅ／ｄｔは、最適エンジン回転数Ｎｅ^＊から現在のエンジン回転数Ｎｅ（ｋ）を引いて得られる差（Ｎｅ^＊－Ｎｅ（ｋ））を時間ステップΔｔで除することにより得られる。（１５）式では、このように特定されたエンジン回転数変化率ｄＮｅ／ｄｔ（最適エンジン回転数Ｎｅ^＊に基づく値）及び充放電量Ｐｃｈｇが、ドライバーからのトルク要求（アクセルペダルの踏み込み量）に基づく要求駆動トルクＴｐとともに、目標状態量として利用されている。したがって、（１５）式によって表された連立方程式を解くことにより、これらの目標状態量を満足する操作量（Ｔｅ、Ｔｇ、Ｔｍ）を得ることができる。得られた操作量は、トルクデバイス（内燃機関２０、ＭＧ１、ＭＧ２）に出力される。その結果、内燃機関２０では、与えられた操作量（Ｔｅ）を実現するために、エンジントルク制御のための各アクチュエータ（スロットルバルブ２２、燃料噴射弁２４及び点火装置２６）が制御される。また、ＭＧ１とＭＧ２の制御に関しては、与えられた操作量（Ｔｇ、Ｔｍ）を実現するように、ＭＧ１及びＭＧ２がインバータ４０によって制御される。

１－２－３．制御ルーチン
図５は、本発明の実施の形態１に係るパワートレーン制御に関する処理のルーチンを示すフローチャートである。本ルーチンは、一例としてパワートレーンシステム１０の起動中（すなわち、１トリップ中）に所定の時間ステップΔｔ毎に繰り返し実行される。すなわち、本実施形態の制御期間τの一例は、車両の１トリップである。

なお、本実施形態の「制御入力決定処理５４」は以下のステップＳ１００～Ｓ１０６の処理と対応しており、「システム制御処理５６」はステップＳ１０８の処理と対応している。より詳細には、「初期値決定処理５９」を含む「随伴変数決定処理５８」はステップＳ１０２及びＳ１０６の処理に対応し、また、ステップＳ１０４の処理は、ＨＶモード用の演算を行う「制御入力算出処理６０」を含む。

図５に示すルーチンでは、制御装置５０は、まず、ステップＳ１００において、現在の（実）ＳＯＣを取得（算出）するとともに、センサ類５２を用いて現在の車速Ｖ及び要求駆動トルクＴｐを取得する。より詳細には、現在の車速Ｖは、例えば車速センサを用いて取得され、要求駆動トルクＴｐは、アクセルポジションセンサにより検出されるアクセルペダルの踏み込み量に応じた値として取得される。その後、処理はステップＳ１０２に進む。

ステップＳ１０２では、制御装置５０は、現在のサンプル時刻ｋの随伴変数ｐ（ｋ）を取得する。制御装置５０のメモリ５０ｂには、前回の制御期間τ（本実施形態では、前回トリップ）中の随伴変数ｐの最終値ｐ_ｆ及び前回の終端ＳＯＣ_ｆＯＬＤが記憶されている。制御装置５０は、本パワートレーン制御の開始後（例えば、パワートレーンシステム１０の起動後）の初回の演算では、初期値ｐ_０の決定のために上述の初期値決定処理５９を実行する。初期値決定処理５９では、前回トリップの最終値ｐ_ｆ（記憶値）が初期値ｐ_０の基本値とされる。ステップＳ１００において取得した現在のＳＯＣ（すなわち、今回の初期ＳＯＣ_０ＮＥＷ）と前回の終端ＳＯＣ_ｆＯＬＤ（記憶値）とを用いて外部充放電量ΔＳＯＣ_ｇが算出されるとともに、外部充放電量ΔＳＯＣ_ｇと補正ゲインＧ１とを用いて外部充放電補正値ｐ_ｇが算出される。そして、基本値（初期値ｐ_０）と外部充放電補正値ｐ_ｇとの和が初期値ｐ_０として算出される。

なお、既述したように、初期値ｐ_０は、前回の制御期間τ中の随伴変数ｐの平均値ｐ_ＡＶＥであってもよい。平均値ｐ_ＡＶＥを用いる例では、個々の制御期間τ中に各サンプル時刻ｋの随伴変数ｐの平均値ｐ_ＡＶＥが算出され、かつ、個々の制御期間τの終了時に最新の平均値ｐ_ＡＶＥがメモリ５０ｂに記憶される。

また、ステップＳ１０２では、本パワートレーン制御の開始後の２回目以降の演算においては、後述のステップＳ１０６の処理によって前回のサンプル時刻（ｋ－１）において算出された随伴変数ｐ（ｋ）が、今回の随伴変数ｐ（ｋ）として用いられる。その後、処理はステップＳ１０４に進む。

ステップＳ１０４では、制御装置５０は、上記（９）式により表される最適化問題を解くことにより、ハミルトニアンＨを最小とする最適制御入力値ｕ^＊（ｋ）（すなわち、Ｎｅ^＊（ｋ）及びＴｅ^＊（ｋ））を探索して算出する。具体的には、ステップＳ１０４の処理は、既述したようにＨＶモード時のハミルトニアンＨ_ＨＶを最小とする最適制御入力値ｕ^＊ＨＶ（ｋ）を探索して算出する制御入力算出処理６０を含む。そして、最適制御入力値ｕ^＊（ｋ）は、制御入力算出処理６０によって算出される最適制御入力値ｕ^＊ＨＶ（ｋ）と、ＥＶモードにおける制御入力値ｕ^ＥＶ（ｋ）（Ｎｅ及びＴｅが共にゼロ）の中から、ハミルトニアンＨを最小とするように決定される。

ステップＳ１０４の後に、処理はステップＳ１０６に進む。ステップＳ１０６では、制御装置５０は、次のサンプル時刻（ｋ＋１）の演算に用いられる随伴変数ｐ（ｋ＋１）を算出する。随伴変数ｐ（ｋ＋１）の算出手法は、（７）式を参照して既述した通りである。算出された随伴変数ｐ（ｋ＋１）は、メモリ５０ｂに記憶される。その後、処理はステップＳ１０８に進む。

ステップＳ１０８では、制御装置５０は、ステップＳ１０４の処理により算出した最適制御入力値ｕ^＊（ｋ）（すなわち、Ｎｅ^＊（ｋ）及びＴｅ^＊（ｋ））をパワートレーンシステム１０に適用するための処理（システム制御処理５６）を実行する。その結果、最適制御入力値ｕ^＊（ｋ）に対応するシステム動作点及びバッテリ４４の充放電量Ｐｃｈｇが特定される。そして、既述したように、最適制御入力値ｕ^＊（ｋ）に基づいて内燃機関２０、ＭＧ１及びＭＧ２が制御される。

１－３．効果
以上説明した本実施形態のパワートレーンシステム１０によれば、（５）式により定義されるハミルトニアンＨが最小となる最適制御入力値ｕ^＊（ｋ）がサンプル時刻ｋ毎に取得される。したがって、燃料消費量ｍ_ｆが最小となるように各サンプル時刻ｋの最適制御入力値ｕ^＊（ｋ）を決定できるようになる。その結果、決定された最適制御入力値ｕ^＊（ｋ）に応じて、バッテリ４４の充放電量Ｐｃｈｇとシステム動作点とを適切に決定できるようになる。このため、本実施形態のパワートレーン制御によれば、各サンプル時刻ｋにおいて燃料消費量ｍ_ｆを低減させられるようにパワートレーンシステム１０を制御することができる。

そして、本実施形態の制御入力決定処理５４によれば、対象とする最適化問題の随伴変数ｐの初期値ｐ_０の基本値として、前回の制御期間τ中の随伴変数ｐの最終値ｐ_ｆ（又は平均値ｐ_ＡＶＥ）が用いられる。これにより、将来（制御入力値ｕを最適化したい期間）の車速Ｖ及び要求駆動トルクＴｐが未知であっても、各サンプル時刻ｋの充放電量Ｐｃｈｇ及びシステム動作点を適切に決定できるようになる。付け加えると、将来の車速Ｖ及び要求駆動トルクＴｐを別途予測する必要がないため、プロセッサ５０ａの演算負荷の低減を図りつつ、各サンプル時刻ｋの充放電量Ｐｃｈｇ及びシステム動作点を適切に決定できる。

付け加えると、随伴変数ｐの初期値ｐ_０の基本値として前回の制御期間τにおける最終値ｐ_ｆ（又は平均値ｐ_ＡＶＥ）を用いることは、次のような意義を有している。まず、制御期間τの途中で実ＳＯＣがＳＯＣ制御中心値から大きく外れると、強制発電又は強制放電が必要になって燃費が悪化する可能性がある。このような事態を回避するためには、終端ＳＯＣ_ｆ（制御期間τの終了時の実ＳＯＣ）がＳＯＣ制御中心値に近づくように制御期間τ中の実ＳＯＣを制御することが望ましい。その理由は、ＳＯＣ制御中心値に近づくように終端ＳＯＣ_ｆを制御できていれば、制御期間τ中の実ＳＯＣをＳＯＣ制御中心値を中心として推移させられるといえるためである。ここで、初期値ｐ_０は終端ＳＯＣ_ｆに影響を及ぼす。すなわち、初期値ｐ_０が変わると、終端ＳＯＣ_ｆが変化する。また、サンプル時刻ｋ毎に算出される随伴変数ｐ（ｋ）は、車両走行履歴（制御期間τの車速Ｖ及び要求駆動トルクＴｐの履歴）に応じて変化する。そして、ハミルトニアンＨに関する（５）式によれば、随伴変数ｐは、燃料消費に対する電力消費の重みに相当する。したがって、前回の制御期間τの初期の実ＳＯＣ（初期ＳＯＣ_０）及び終端ＳＯＣ_ｆがともにＳＯＣ制御中心値（例えば、６０％）である例において、今回の制御期間τの初期の実ＳＯＣ_０が６０％であって今回の制御期間τ中の車両走行履歴（Ｖ、Ｔｐの履歴）が前回の制御期間τのものと同様であれば、初期値ｐ_０の基本値として上記の最終値ｐ_ｆ又は平均値ｐ_ＡＶＥを用いることにより、今回の制御期間τにおいても前回と同様に６０％の終端ＳＯＣ_ｆを狙える（すなわち、今回の制御期間τ中の実ＳＯＣがＳＯＣ制御中心値を中心として推移する）ことを期待できる。以上のことから、初期値ｐ_０の基本値として前回の制御期間τ中の最終値ｐ_ｆ（又は平均値ｐ_ＡＶＥ）を使用するということは、過去（前回）の走行履歴に基づいて燃料とバッテリ４４の電力に関する将来の適切な使用比率を予測していることに相当するといえる。換言すると、初期値ｐ_０の基本値として前回の制御期間τ中の最終値ｐ_ｆ又は平均値ｐ_ＡＶＥを用いることは、前回の制御期間τにおいて得られた情報から、今回の制御期間τにおいて狙いとする（すなわち、ＳＯＣ制御中心値に近づく）終端ＳＯＣ_ｆを実現可能な初期値ｐ_０の基本値を予測しているといえる。

次に、外部充放電量ΔＳＯＣ_ｇに基づく外部充放電補正値ｐ_ｇを用いて行われる初期値ｐ_０の基本値（前回の制御期間τの最終値ｐ_ｆ又は平均値ｐ_ＡＶＥ）の補正の効果について説明する。

図６は、制御停止期間ｔｓ中の外部充電又は外部放電の実行に伴う課題について説明するためのタイムチャートである。図６には、前回のトリップ（制御期間τ）において、初期ＳＯＣ_０ＯＬＤ及び終端ＳＯＣ_ｆＯＬＤがともにＳＯＣ制御中心値と一致するように実ＳＯＣが推移した例（比較例）が表されている。今回トリップの車両走行履歴（すなわち、車両の走り方）が前回トリップのものと同様であれば、このような前回トリップを経験した後の次のトリップ（今回トリップ）における実ＳＯＣ波形は次のようになる。すなわち、初期値ｐ_０の基本値として前回の制御期間τ中の最終値ｐ_ｆ（又は平均値ｐ_ＡＶＥ）を使用することにより、今回トリップの実ＳＯＣ波形は、図６中に破線で示すように前回のものと同様となる。このため、今回トリップにおいて、前回トリップと同様の終端ＳＯＣ_ｆを狙えることを期待できる。

ただし、上記の効果は、図６中の破線の例のように、今回の初期ＳＯＣ_０ＮＥＷが前回の終端ＳＯＣ_ｆＯＬＤと等しいこと、すなわち、制御停止期間ｔｓ中に外部充電又は外部放電がなされていないことを条件として得られる。例えば、制御停止期間ｔｓ中に外部充電がなされた場合には、今回の初期ＳＯＣ_０ＮＥＷが前回の終端ＳＯＣ_ｆＯＬＤに対して増加する。その結果、今回トリップの走り方が前回トリップのものと同じであっても、図６中に実線で示すように今回トリップの実ＳＯＣ波形は、外部充電の有無による初期ＳＯＣ_０ＮＥＷの差分だけ破線の実ＳＯＣ波形からＳＯＣ増加側に平行移動したような波形となる。このことは、制御停止期間ｔｓ中に外部放電がなされたために今回の初期ＳＯＣ_０ＮＥＷが前回の終端ＳＯＣ_ｆＯＬＤに対して減少した場合についても同様である。すなわち、走り方が同じであっても、今回トリップの実ＳＯＣ波形は、外部放電の有無による初期ＳＯＣ_０ＮＥＷの差分だけ破線の実ＳＯＣ波形からＳＯＣ減少側に平行移動したような波形となる。

上述のように、初期値ｐ_０として前回の制御期間τ中の最終値ｐ_ｆ（又は平均値ｐ_ＡＶＥ）をそのまま使用した場合には、制御停止期間ｔｓ中に外部充電又は外部放電がなされると、走り方が同じであっても、今回トリップの終端ＳＯＣ_ｆをＳＯＣ制御中心値に合わせることができなくなる。この理由について説明を補足する。初期値ｐ_０の基本値（最終値ｐ_ｆ又は平均値ｐ_ＡＶＥ）は、前回トリップにおける車両の走り方（代表的には車速Ｖ及び車両駆動トルク（要求駆動トルクＴｐ）の変化の仕方）及び充放電収支（終端ＳＯＣ_ｆと初期ＳＯＣ_０との差）を前提に決定されたものであり、今回トリップにおいて前回トリップと同じ走り方で同じ充放電収支を実現したとしても、外部充電又は外部放電があれば初期ＳＯＣ_０ＮＥＷが変化するので狙いの終端ＳＯＣ_ｆが得られないためである。

このような課題に鑑み、本実施形態の初期値決定処理５９によれば、制御停止期間ｔｓ中に外部充電又は外部放電がなされた場合には、外部充放電量ΔＳＯＣ_ｇに基づく外部充放電補正値ｐ_ｇによって初期値ｐ_０の基本値（前回トリップの最終値ｐ_ｆ又は平均値ｐ_ＡＶＥ）が補正される。

図７は、本発明の実施の形態１に係る初期値決定処理５９による補正の効果を表したタイムチャートである。図８は、ＳＯＣ（ｋ＋１）と随伴変数ｐ（ｋ）との関係を表したグラフである。図８に示すように、サンプル時刻ｋの随伴変数ｐ（ｋ）と次回のサンプル時刻（ｋ＋１）のＳＯＣ（ｋ＋１）とは、互いに反比例の関係にある。より詳細には、現在のサンプル時刻ｋの随伴変数ｐ（ｋ）及び（実）ＳＯＣ（ｋ）に応じて最適制御入力値ｕ^＊（ｋ）が決定されると、決定した最適制御入力値ｕ^＊（ｋ）に応じて次回のサンプル時刻（ｋ＋１）の（実）ＳＯＣ（ｋ＋１）が変化する。具体的には、随伴変数ｐ（ｋ）が大きくなると、ＳＯＣ（ｋ＋１）はＳＯＣ（ｋ）に対して低くなる。逆に、随伴変数ｐ（ｋ）が小さくなると、ＳＯＣ（ｋ＋１）はＳＯＣ（ｋ）に対して高くなる。

初期値決定処理５９によれば、（１２）～（１４）式によって表されるように、制御停止期間ｔｓ中に外部充電がなされた場合には、初期値ｐ_０は、外部充放電補正値ｐ_ｇがゼロの場合と比べて大きくなる。このように初期値ｐ_０（初回のサンプル時刻（ｋ＝０）の随伴変数ｐの値）が大きくなると、（７）式より、次のサンプル時刻（ｋ＝１）の随伴変数ｐ_１も大きくなる。このことは、その次のサンプル時刻（ｋ＝２）以降の随伴変数ｐについても同様である。随伴変数ｐ（ｋ）が大きくなることは、図８に示すようにＳＯＣ（ｋ＋１）が低くなること（すなわち、バッテリ４４の放電が促進されること）を意味する。したがって、外部充放電補正値ｐ_ｇを用いて初期値ｐ_０を大きくする補正が行われると、図７に示すように、今回トリップ中の各時点の実ＳＯＣが、当該補正が行われない場合（破線の波形）と比べて放電側（ＳＯＣを小さくする側）の値をとるようになる。

同様に、初期値決定処理５９によれば、制御停止期間ｔｓ中に外部放電がなされた場合には、初期値ｐ_０は、外部充放電補正値ｐ_ｇがゼロの場合と比べて小さくなる。その結果、次のサンプル時刻（ｋ＝１）以降の随伴変数ｐ（ｋ）も小さくなるので、ＳＯＣ（ｋ＋１）が高くなる（すなわち、バッテリ４４の充電が促進される）。したがって、外部充放電補正値ｐ_ｇを用いて初期値ｐ_０を小さくする補正が行われると、図７に示すように、今回トリップ中の各時点の実ＳＯＣが、当該補正が行われない場合（破線の波形）と比べて充電側（ＳＯＣを大きくする側）の値をとるようになる。

以上説明した理由により、初期値決定処理５９によれば、制御停止期間ｔｓ中に外部充電又は外部放電がなされた場合に、図７に示すように、当該補正がなされない場合と比べて、今回トリップの終端ＳＯＣ_ｆをＳＯＣ制御中心値に近づけることが可能となる。そして、このことは、当該補正がなされない場合と比べて、今回トリップ中の実ＳＯＣをＳＯＣ制御中心値により近づけて推移させられることを意味する。その結果、制御停止期間ｔｓ中に外部電源４５又は電気製品４６との充電又は放電がなされた場合であっても、上記補正がなされない場合と比べて、制御入力値ｕ（ｋ）をより適切に決定できるようになる。

２．実施の形態２
次に、図９～図１４を参照して、本発明の実施の形態２について説明する。

２－１．パワートレーン制御
実施の形態２は、以下の点において、上述した実施の形態１と相違している。すなわち、実施の形態１では、補正ゲインＧ１として固定値が用いられている。これに対し、本実施形態では、外部充放電の影響を考慮して初期値ｐ_０をより適切に設定するために、補正ゲインＧ１はハイブリッド車両の走行履歴に基づいて設定される。

図９は、随伴変数ｐとＳＯＣとの関係を車両走行モード毎に分けて表したグラフである。図９中には、モードＡ～Ｇが例示されている。図９中の各直線は、各車両走行モードの使用中に所定時間毎にサンプリングされた随伴変数ｐ及び（実）ＳＯＣのデータのプロット値の近似直線に相当する。

図９に示す例では、モードＡ～Ｇを用いて多様な車両走行パターンの情報が得られている。より詳細には、モードＡは、市街地モード、郊外モード及び高速道路モードを含むＷＬＴＣ（World harmonized Light vehicles Test Cycle）モードに相当する。モードＢは、米国の高速道路モードの１つであるＨＷＹモードに相当する。モードＣは、日本の市街地の渋滞走行を模擬したモードに相当する。モードＤは、大都市の渋滞走行を模擬した他の渋滞走行モードに相当する。モードＥは、米国の他の高速道路モードであるＵＳ０６モードに相当する。モードＦは、ロサンゼルス郊外の走行パターンを模擬したＬＡ＃４モードに相当する。モードＧは、日本で過去に用いられていたＪＣ０８モードに相当する。

車両走行モードが異なると、車両の走り方（代表的には車速Ｖ及び車両駆動トルク（要求駆動トルクＴｐ）の変化の仕方）が異なるものとなる。その結果、図９に示すように、随伴変数ｐとＳＯＣとの関係が車両走行モードに応じて異なるものとなる。より詳細には、図９中に示す各近似直線の傾き（比ΔＰ／ΔＳＯＣ）の相違から分かるように、車両走行モードの相違に起因して、ＳＯＣに対する随伴変数ｐの感度が変化する。図１０は、ＳＯＣに対する随伴変数ｐの感度、すなわち、比ΔＰ／ΔＳＯＣをモードＡ～モードＧの間で比較して表したグラフである。

上述した実施の形態１の初期値決定処理５９によれば、外部充放電量ΔＳＯＣ_ｇに応じて初期値ｐ_０が補正される。その一方で、図９及び図１０を参照して上述したように、ＳＯＣに対する随伴変数ｐの感度は車両走行モード（車両走行パターン）に応じて変化する。このことは、外部充放電量ΔＳＯＣ_ｇから外部充放電補正値ｐ_ｇに変換するために用いられる補正ゲインＧ１の最適値が車両走行パターン（すなわち、車両の走り方）に応じて変化することを意味する。

したがって、上記のような初期値決定処理５９において用いられる補正ゲインＧ１が固定値であると、前回トリップの車両の走り方が補正ゲインＧ１の適合時に想定した車両の走り方と異なる場合には補正ゲインＧ１が適切なものでなくなる可能性がある。このことは、ハイブリッド車両の燃費悪化に繋がり得る。付け加えると、用いられる補正ゲインＧ１が固定値でなかったとしても、ＳＯＣに対する随伴変数ｐの感度が前回トリップの車両の走り方に応じて変化する点を考慮して補正ゲインＧ１が変更されるようになっていなければ、上記と同様の課題が生じ得る。

図１１は、本発明の実施の形態２に係るパワートレーン制御に関連して制御装置６２が実行する処理を示すブロック図である。制御装置６２の処理は、制御入力決定処理５４’に含まれる随伴変数決定処理５８’が初期値決定処理５９に代えて初期値決定処理６４を含む点において、制御装置５０の処理と相違している。

初期値決定処理６４は、補正ゲインＧ１をハイブリッド車両の走行履歴に基づいて設定する。走行履歴の具体例としては、前回の制御期間τ（前回トリップ）中の車速Ｖ及び要求駆動トルクＴｐのそれぞれの平均値である第１平均車速ａｖｅＶ１及び第１平均駆動トルクａｖｅＴｐ１が用いられる。すなわち、本実施形態では、補正ゲインＧ１の設定のために、前回トリップの走行履歴が用いられる。

図１２は、第１平均車速ａｖｅＶ１及び第１平均駆動トルクａｖｅＴｐ１と、比ΔＰ／ΔＳＯＣとの関係を車両走行モード毎に表したグラフである。より詳細には、各車両走行モード（モードＡ～Ｇ）に関するプロット点は、それぞれのモードの使用中の全体のデータ（ａｖｅＶ１，ａｖｅＴｐ１，ΔＰ／ΔＳＯＣ）の平均値に対応している。したがって、モードＡ～Ｇは、それぞれのモード中の平均的な車両走行負荷（車両駆動トルク）及び平均的な車速に応じて、図１２中に示されるように分布する。

図１２中の近似直線Ｌは、モードＡ～Ｇの各プロット点に関する近似直線に相当する。この近似直線Ｌによって表されるように、比ΔＰ／ΔＳＯＣは、第１平均車速ａｖｅＶ１が高いほど高くなる。また、比ΔＰ／ΔＳＯＣは、第１平均駆動トルクａｖｅＴｐ１が高いほど高くなる。このように、図１２に示す３次元的な関係から、比ΔＰ／ΔＳＯＣは、第１平均車速ａｖｅＶ１及び第１平均駆動トルクａｖｅＴｐ１と高い相関を有することが分かる。

上述の知見より、第１平均車速ａｖｅＶ１及び第１平均駆動トルクａｖｅＴｐ１を利用することで、前回トリップの走行履歴に基づく比ΔＰ／ΔＳＯＣ（すなわち、ＳＯＣに対する随伴変数ｐの感度）を適切に取得できるようになる。

以下の（１６）式は、本実施形態の初期値決定処理６４において用いられる補正ゲインＧ１の算出式の一例を示している。（１６）式に示すように、本実施形態では、第１平均車速ａｖｅＶ１と第１平均駆動トルクａｖｅＴｐ１の関数ｆ（ａｖｅＶ１，ａｖｅＴｐ１）である比ΔＰ／ΔＳＯＣが補正ゲインＧ１として用いられる。

（１６）式の例のように、補正ゲインＧ１、すなわち、比ΔＰ／ΔＳＯＣの近似式としては、第１平均車速ａｖｅＶ１及び第１平均駆動トルクａｖｅＴｐ１を独立変数とする一次式を用いることができる。この近似式中の係数Ａ、Ｂ及びＣは、例えば、最小二乗法等の重回帰分析を用いて予め設定される。このような近似式の利用により、前回トリップの走行履歴に基づく比ΔＰ／ΔＳＯＣ（すなわち、ＳＯＣに対する随伴変数ｐの感度）に応じた適切な補正ゲインＧ１を取得（推定）できるようになる。

付け加えると、図１２に示すような関係を有する比ΔＰ／ΔＳＯＣの利用により、補正ゲインＧ１は、第１平均車速ａｖｅＶ１が高いほど大きくなるように設定されることになる。また、補正ゲインＧ１は、第１平均駆動トルクａｖｅＴｐ１が高いほど大きくなるように設定されることになる。

なお、補正ゲインＧ１の設定のために用いられる近似式は、上記の一次式の例に限られず、二次以上の高次の式であってもよい。また、補正ゲインＧ１の設定のために、このような近似式に代え、第１平均車速ａｖｅＶ１及び第１平均駆動トルクａｖｅＴｐ１に対する比ΔＰ／ΔＳＯＣの関係を定めたマップ（図示省略）が用いられてもよい。

２－１－１．制御ルーチン
図１３は、本発明の実施の形態２に係る初期値決定処理６４のサブルーチンを示すフローチャートである。制御装置６２は、このサブルーチンの処理を、上述の図５に示すルーチン（メインルーチン）と並行して実行する。そして、このサブルーチンの処理は、上記メインルーチンのステップＳ１０２において初期値ｐ_０を算出する場合に実行される。

図１３に示すサブルーチンでは、制御装置６２は、まず、ステップＳ２００において、前回の終端ＳＯＣ_ｆＯＬＤ（記憶値）とステップＳ１００において取得した現在のＳＯＣ（すなわち、今回の初期ＳＯＣ_０ＮＥＷ）とを用いて外部充放電量ΔＳＯＣ_ｇを算出する。その後、処理は、ステップＳ２０２に進む。

制御装置６２は、各トリップの終了時に、トリップ中の車両走行履歴（第１平均車速ａｖｅＶ１及び第１平均駆動トルクａｖｅＴｐ１）を算出してメモリに記憶している。ステップＳ２０２では、制御装置６２は、そのような前回トリップの第１平均車速ａｖｅＶ１及び第１平均駆動トルクａｖｅＴｐ１の記憶値を取得する。

第１平均車速ａｖｅＶ１の算出に用いられる車速Ｖは、例えば車速センサを用いて取得可能である。第１平均駆動トルクａｖｅＴｐ１の算出に用いられる車両駆動トルクは、例えば要求駆動トルクＴｐである。要求駆動トルクＴｐは、既述したように、例えば、アクセルポジションセンサにより検出されるアクセルペダルの踏み込み量に応じた値として取得可能である。また、第１平均駆動トルクａｖｅＴｐ１の算出のために用いられる車両駆動トルクは、実駆動トルクであってもよい。

ステップＳ２０２の後に、処理はステップＳ２０４に進む。ステップＳ２０４では、制御装置６２は、補正ゲインＧ１を決定する。具体的には、制御装置６２は、上記（１６）式を用いて、ステップＳ２０２の処理により取得した第１平均車速ａｖｅＶ１及び第１平均駆動トルクａｖｅＴｐ１に応じた補正ゲインＧ１を算出する。その後、処理はステップＳ２０６に進む。

ステップＳ２０６では、制御装置６２は、ステップＳ２００において取得した外部充放電量ΔＳＯＣ_ｇとステップＳ２０４において決定された補正ゲインＧ１との積を外部充放電補正値ｐ_ｇとして算出する。その後、処理は、ステップＳ２０８に進む。

ステップＳ２０８では、制御装置６２は、上記（１２）式に従い、基本値（前回トリップの最終値ｐ_ｆ）とステップＳ２０６の処理により算出した外部充放電補正値ｐ_ｇとの和を算出し、その算出値を初期値ｐ_０として決定する。

２－２．効果
図１４は、本発明の実施の形態２に係る初期値決定処理６４による補正の効果を表したタイムチャートである。以上説明した本実施形態の初期値決定処理６４によれば、補正ゲインＧ１は、前回トリップの車両走行履歴に基づいて設定される。既述したように、ＳＯＣに対する随伴変数ｐの感度（すなわち、比ΔＰ／ΔＳＯＣ）は、車両の走り方に応じて変化する。このため、前回トリップの車両走行履歴（車両の走り方）を考慮して補正ゲインＧ１を決定することにより、補正ゲインＧ１を固定値とする例と比べて、前回トリップの走り方に応じた適切な補正ゲインＧ１を用いつつ外部充放電量ΔＳＯＣ_ｇに基づく外部充放電補正値ｐ_ｇを適切に算出できるようになる。その結果、図１４に示すように、補正ゲインＧ１を固定値とする例（図７参照）と比べて、今回トリップの終端ＳＯＣ_ｆをＳＯＣ制御中心値により近づけることが可能となる。

より詳細には、初期値決定処理６４によれば、補正ゲインＧ１は、第１平均車速ａｖｅＶ１が高いほど大きくなるように設定される。これにより、ＳＯＣに対する随伴変数ｐの感度を示す比ΔＰ／ΔＳＯＣが、第１平均車速ａｖｅＶ１が高いほど高くなるという特性を考慮して、補正ゲインＧ１を適切に設定できるようになる。さらに、補正ゲインＧ１は、第１平均駆動トルクａｖｅＴｐ１が高いほど大きくなるように設定される。これにより、上記感度を示す比ΔＰ／ΔＳＯＣが、第１平均駆動トルクａｖｅＴｐ１が高いほど高くなるという特性を考慮して、補正ゲインＧ１を適切に設定できるようになる。

３．実施の形態３
次に、図１５～図１７を参照して、本発明の実施の形態３について説明する。

３－１．パワートレーンシステムの構成例
実施の形態３に係るパワートレーンシステムは、制御装置５０に代え、後述の図１６に示す制御装置７０を備えている点を除き、実施の形態１に係るパワートレーンシステム１０と同様である。

３－２．パワートレーン制御
３－２－１．実施の形態１、２に係る制御入力決定処理５４、５４’の課題
上述した制御入力決定処理５４、５４’によれば、随伴変数ｐの初期値ｐ_０の基本値として、前回の制御期間τ中の随伴変数ｐの最終値ｐ_ｆ（又は平均値ｐ_ＡＶＥ）が用いられる。既述したように、初期値ｐ_０の基本値として前回の制御期間τ中の最終値ｐ_ｆ（又は平均値ｐ_ＡＶＥ）を使用するということは、過去（前回）の走行履歴に基づいて燃料とバッテリ４４の電力に関する将来の適切な使用比率を予測していることに相当するといえる。換言すると、前回の制御期間τの最終値ｐ_ｆ又は平均値Ｐ_ＡＶＥの使用は、将来の走行パターン（Ｖ、Ｔｐの推移）を把握せずに、将来の走行パターンが過去の走行パターンと同様になると想定して随伴変数ｐを決定することに相当するといえる。

しかしながら、過去の走行パターンと将来の走行パターンとは、必ずしも同じになるとは限られない。具体的には、制御入力決定処理５４が使用される例では、仮に将来の走行パターンが過去の走行パターンに対して異なるものになると、実ＳＯＣ軌道（時間に対して変動する実ＳＯＣの変動中心線）がＳＯＣ制御中心値（後述の図１５参照）からずれてしまう可能性がある。

図１５は、上述した実施の形態１、２に係るパワートレーン制御の課題と、その対策を有する実施の形態３に係るパワートレーン制御とを比較しながら説明するためのタイムチャートである。図１５中に破線で示す実ＳＯＣの波形は、実施の形態１、２に係るパワートレーン制御の実行中にＳＯＣ制御中心値に対する実ＳＯＣの軌道のずれが生じた例を示している。より詳細には、図１５は、時間の経過とともに、ＳＯＣ制御中心値（例えば、６０％）に対して実ＳＯＣが低くなる側（すなわち、放電側）に実ＳＯＣ軌道がずれた例を示している。実ＳＯＣ軌道は、図１５に示す例とは逆に、ＳＯＣ制御中心値に対して実ＳＯＣが高くなる側（すなわち、充電側）にずれることもある。ここでいうＳＯＣ制御中心値は、本発明に係る「目標充電率」の一例に相当する。

３－２－２．過去の走行パターンに対する将来の走行パターンの変化を考慮したパワートレーン制御の概要
図１６は、本発明の実施の形態３に係るパワートレーン制御に関連して制御装置７０が実行する処理を示すブロック図である。制御装置７０は、以下に説明する点において、実施の形態１の制御装置５０と相違している。すなわち、制御装置７０がパワートレーン制御に関連して実行する処理は、「制御入力決定処理５４」に代え、「制御入力決定処理７２」を含む。そして、この制御入力決定処理７２に含まれる「随伴変数決定処理７４」は、図１５を参照して説明した課題に鑑みて「随伴変数修正処理７４ａ」を含み、かつ初期値決定処理５９を含むという点において、実施の形態１の随伴変数決定処理５８と相違している。このように、本実施形態においても、制御停止期間ｔｓ中に外部充電又は外部放電がなされた場合には、随伴変数ｐの初期値ｐ_０の基本値は、外部充放電量ΔＳＯＣ_ｇに基づく外部充放電補正値ｐ_ｇによって補正される。なお、制御装置７０の処理は、初期値決定処理５９に代え、初期値決定処理６４を含んでいてもよい。このことは、後述の実施の形態４についても同様である。

随伴変数修正処理７４ａは、ＳＯＣの偏差ΔＳＯＣ_fbに基づいて随伴変数ｐを修正するために実行される。偏差ΔＳＯＣ_ｆｂは、実ＳＯＣ軌道からＳＯＣ制御中心値を引いて得られる差に相当する。より詳細には、随伴変数修正処理７４ａによれば、図１５に示すように、一例として時間ステップΔｔ毎に（すなわち、本実施形態のパワートレーン制御の制御周期で）、偏差ΔＳＯＣ_ｆｂが算出される。実ＳＯＣ軌道は、例えば、実ＳＯＣの移動平均線を算出することによって取得できる。なお、偏差ΔＳＯＣ_ｆｂは、本発明に係る「目標充電率と実充電率との差」の一例に相当する。

随伴変数修正処理７４ａは、上記のように算出される偏差ΔＳＯＣ_ｆｂに基づいて、一例として時間ステップΔｔ毎に（制御周期で）随伴変数ｐを修正する。

次に、偏差ΔＳＯＣ_ｆｂに基づく随伴変数ｐの修正手法の一例を示す。随伴変数ｐは、以下の（１７）式に従って修正される。すなわち、修正後の随伴変数ｐは、修正前の随伴変数ｐに対して、偏差ΔＳＯＣ_ｆｂに所定のＦＢゲインＧ２を乗じて得られる補正値を加えることにより算出される。ＦＢゲインＧ２は正の値である。したがって、この補正値の絶対値は、偏差ΔＳＯＣ_ｆｂの絶対値が大きいほど大きくなる。なお、ＦＢゲインＧ２は、本発明に係る「第２ゲイン」の一例に相当する。

偏差ΔＳＯＣ_ｆｂ（＝実ＳＯＣ軌道－ＳＯＣ制御中心値）は、ＳＯＣ制御中心値に対して実ＳＯＣ軌道が放電側にずれている場合（実ＳＯＣ＜ＳＯＣ制御中心値）には負となる。したがって、（１７）式によれば、実ＳＯＣ軌道が放電側にずれている場合には、負の補正値が修正前の随伴変数ｐに加算されるので、随伴変数ｐが小さくなるように修正される。より詳細には、随伴変数ｐは、偏差ΔＳＯＣ_ｆｂ（の絶対値）が大きいほど小さくなるように修正される。修正後の随伴変数ｐ（ｋ）は、同一時間ステップΔｔにおける最適制御入力値ｕ^＊（ｋ）の算出に直ちに反映される。これにより、修正前と比べて実ＳＯＣの低下が抑制されるように（すなわち、バッテリ４４の電力消費を抑えるために内燃機関２０の作動時間が増えるように）パワートレーンシステム１０が制御されることになる。その結果、偏差ΔＳＯＣ_ｆｂが減少する。

一方、実ＳＯＣ軌道が充電側にずれている場合（実ＳＯＣ＞ＳＯＣ制御中心値）には、偏差ΔＳＯＣ_ｆｂは正となる。したがって、（１７）式によれば、実ＳＯＣ軌道が充電側にずれている場合には、正の補正値が修正前の随伴変数ｐに加算されるので、修正前の随伴変数ｐが大きくなるように修正される。より詳細には、随伴変数ｐは、偏差ΔＳＯＣ_ｆｂ（の絶対値）が大きいほど大きくなるように修正される。これにより、修正前と比べて実ＳＯＣの低下（バッテリ４４の電力消費）が促進されるようにパワートレーンシステム１０が制御されることになる。その結果、偏差ΔＳＯＣ_ｆｂが減少する。

既述したように、初期値ｐ_０の基本値として前回の制御期間τ中の最終値ｐ_ｆ（又は平均値ｐ_ＡＶＥ）を使用するということは、過去（前回）の走行履歴に基づいて燃料とバッテリ４４の電力に関する将来の適切な使用比率を予測していることに相当するといえる。この点に関し、随伴変数修正処理７４ａによれば、そのような使用比率の予測値（前回の最終値ｐ_ｆ又は平均値ｐ_ＡＶＥ）を修正するフィードバック処理が偏差ΔＳＯＣ_ｆｂに基づいて時間ステップΔｔ毎に行われることになるといえる。

３－２－３．制御ルーチン
図１７は、本発明の実施の形態３に係るパワートレーン制御に関する処理のルーチンを示すフローチャートである。なお、図１７に示すルーチン中のステップＳ１００～Ｓ１０８の処理については、実施の形態１において既述した通りである。本ルーチンのステップＳ３００及びＳ３０２の処理は、「随伴変数修正処理７４ａ」に対応している。

図１７に示すルーチンでは、処理はステップＳ１０２の後にステップＳ３００に進む。ステップＳ３００では、制御装置７０は、現在のサンプル時刻ｋにおける偏差ΔＳＯＣ_ｆｂ（ｋ）を算出する。偏差ΔＳＯＣ_ｆｂの算出手法の一例は、図１５、１６を参照して説明した通りである。その後、処理はステップＳ３０２に進む。

ステップＳ３０２では、制御装置７０は、ステップＳ３００において算出した偏差ΔＳＯＣ_ｆｂ（ｋ）に応じて、ステップＳ１０２において取得された随伴変数ｐ（ｋ）を修正する。具体的には、一例として、上述の（１７）式を利用して補正値が算出され、算出された補正値によって随伴変数ｐ（ｋ）が修正される。その結果、ＳＯＣ制御中心値に対して実ＳＯＣ軌道が放電側にずれている場合には、随伴変数ｐ（ｋ）が小さくされ、逆に、実ＳＯＣ軌道が充電側にずれている場合には、随伴変数ｐ（ｋ）が大きくされる。なお、偏差ΔＳＯＣ_ｆｂ（ｋ）がゼロである場合には、随伴変数ｐ（ｋ）は修正されない。

その後、処理はステップＳ１０４に進む。したがって、ステップＳ３０２の処理による修正後の随伴変数ｐ（ｋ）が、ステップＳ１０４の処理による最適制御入力値ｕ^＊（ｋ）の演算に反映される。

３－３．効果
以上説明したように、本実施形態の制御装置７０によって実行される制御入力決定処理７２の随伴変数決定処理７４は、随伴変数修正処理７４ａを含む。随伴変数修正処理７４ａによれば、パワートレーン制御の実行中に、偏差ΔＳＯＣ_ｆｂに応じて随伴変数ｐが時間ステップΔｔ毎に修正される。これにより、図１５中の２つの実ＳＯＣ波形を比較すると分かるように、実線の例（ＳＯＣのフィードバック（随伴変数修正処理７４ａ）を伴う実施の形態３）によれば、破線の例（ＳＯＣのフィードバックを伴わない実施の形態１、２）と比べて、パワートレーン制御の実行中の実ＳＯＣ軌道の中心値をＳＯＣ制御中心値に近づけることが可能となる。換言すると、今回の制御期間τの終端ＳＯＣ_ｆをＳＯＣ制御中心値に近づけることが可能となる。このため、制御期間τの途中で強制充電又は強制放電が必要となることを抑制できる。以上のように、本実施形態のパワートレーン制御によれば、実施の形態１、２のパワートレーン制御と同様に車速Ｖ及び要求駆動トルクＴｐの予測を不要としつつ、過去（前回の制御期間τ）の走行パターンに対して今回の制御期間τの走行パターンが異なる場合であっても燃料消費量ｍ_ｆを最小とするシステム動作点及びバッテリ４４の充放電量Ｐｃｈｇをより適切に決定できるようになる。

３－４．随伴変数修正処理の実行時期の他の例
上述した実施の形態３においては、時間ステップΔｔ毎に随伴変数ｐを修正する随伴変数修正処理７４ａを例に挙げた。しかしながら、本発明に係る「随伴変数修正処理」は、上記の例に代え、例えば、時間ステップΔｔ（パワートレーン制御の制御周期）よりも長い任意の周期で、パワートレーン制御が行われる制御期間中に繰り返し実行されてもよい。

４．実施の形態４
次に、図１８～図２１を参照して、本発明の実施の形態４について説明する。

４－１．パワートレーンシステムの構成例
実施の形態４に係るパワートレーンシステムは、制御装置７０に代え、後述の図１８に示す制御装置８０を備えている点を除き、実施の形態３に係るパワートレーンシステムと同様である。

４－２．パワートレーン制御
４－２－１．実施の形態３に係る随伴変数修正処理７４ａの課題
上述した実施の形態３の随伴変数修正処理７４ａによれば、ＳＯＣの偏差ΔＳＯＣ_ｆｂ（より詳細には、狙いとするＳＯＣ軌道（ＳＯＣ制御中心値）に対する実ＳＯＣ軌道のずれ量）に応じて随伴変数ｐが修正される（フィードバック処理）。その一方で、図９及び図１０を参照して実施の形態２で既に説明したように、ＳＯＣに対する随伴変数ｐの感度は車両走行モード（車両走行パターン）に応じて変化する。このことは、偏差ΔＳＯＣ_ｆｂから補正値（上述の（１７）式の右辺第２項）に変換するために用いられるＦＢゲインＧ２の最適値が車両走行パターンに応じて変化することを意味する。

したがって、上記のようなフィードバック処理において用いられるＦＢゲインＧ２が固定値であると、ＦＢゲインＧ２の適合時に想定した車両走行パターン以外の車両走行パターンの下では、ＦＢゲインＧ２が適切なものでなくなる可能性がある。このことは、ハイブリッド車両の燃費悪化に繋がり得る。付け加えると、用いられるＦＢゲインＧ２が固定値でなかったとしても、ＳＯＣに対する随伴変数ｐの感度が車両走行モードに応じて変化する点を考慮してＦＢゲインＧ２が変更されるようになっていなければ、上記と同様の課題が生じ得る。

４－２－２．実施の形態４に係る随伴変数修正処理８４ａの概要
図１８は、本発明の実施の形態４に係るパワートレーン制御に関連して制御装置８０が実行する処理を示すブロック図である。制御装置８０は、以下に説明する点において、実施の形態３の制御装置７０と相違している。すなわち、制御装置８０がパワートレーン制御に関連して実行する処理は、「制御入力決定処理８２」と「システム制御処理５６」とを含む。そして、この制御入力決定処理８２に含まれる「随伴変数決定処理８４」は、図９及び図１０を参照して説明した課題に鑑みて「随伴変数修正処理７４ａ」に代えて「随伴変数修正処理８４ａ」を含み、かつ初期値決定処理５９を含むという点において、実施の形態３の随伴変数決定処理７４と相違している。

随伴変数修正処理８４ａは、以下のように設定されるＦＢゲインＧ２を用いるという点において、随伴変数修正処理７４ａと相違している。具体的には、本実施形態で用いられるＦＢゲインＧ２は、ハイブリッド車両の走行履歴に基づいて設定される。走行履歴の具体例として、第２平均車速ａｖｅＶ２と第２平均駆動トルクａｖｅＴｐ２とが用いられる。

第２平均車速ａｖｅＶ２は、所定期間における車速Ｖの平均値に相当する。第２平均駆動トルクａｖｅＴｐ２は、所定期間における車両駆動トルク（要求駆動トルクＴｐ）の平均値に相当する。これらの所定期間の一例は、現時点から遡って３０秒である。このように、本実施形態では、ＦＢゲインＧ２の設定のために、現時点に対して直前の走行履歴が用いられる。付け加えると、ＦＢゲインＧ２の設定手法は、実施の形態２の補正ゲインＧ１の設定手法と類似している。ただし、補正ゲインＧ１の設定に用いられる第１平均車速ａｖｅＶ１及び第１平均駆動トルクａｖｅＴｐ１は、それぞれ「前回トリップ中の車速ｖ及び要求駆動トルクＴｐの平均値である」点において、本実施形態で用いられる第２平均車速ａｖｅＶ２及び第２平均駆動トルクａｖｅＴｐ２と相違している。

なお、第２平均車速ａｖｅＶ２の算出に関する所定期間は本発明に係る「第１所定期間」の一例に相当し、第２平均駆動トルクａｖｅＴｐ２の算出に関する所定期間は本発明に係る「第２所定期間」の一例に相当する。上記の例では、第１所定期間と第２所定期間とは同じであるが、これらは必ずしも同じでなくてもよい。

ＦＢゲインＧ２の設定に用いられる第２平均車速ａｖｅＶ２及び第２平均駆動トルクａｖｅＴｐ２は、既述したように、車速ｖ及び要求駆動トルクＴｐの各平均値の算出対象となる期間において、第１平均車速ａｖｅＶ１及び第１平均駆動トルクａｖｅＴｐ１と異なっている。しかしながら、第２平均車速ａｖｅＶ２及び第２平均駆動トルクａｖｅＴｐ２についても、図１２に表されるような第２平均車速ａｖｅＶ２及び第２平均駆動トルクａｖｅＴｐ２と比ΔＰ／ΔＳＯＣとの関係が同様に得られる。そして、図１２から得られる上述の知見（実施の形態２参照）より、第２平均車速ａｖｅＶ２及び第２平均駆動トルクａｖｅＴｐ２を利用することで、現時点に対する直前の走行履歴に基づく比ΔＰ／ΔＳＯＣ（すなわち、ＳＯＣに対する随伴変数ｐの感度）を適切に取得できるようになる。

以下の（１８）式は、本実施形態の随伴変数修正処理８４ａにおいて用いられるＦＢゲインＧ２の算出式の一例を示している。（１８）式に示すように、本実施形態では、第２平均車速ａｖｅＶ２と第２平均駆動トルクａｖｅＴｐ２の関数である比ΔＰ／ΔＳＯＣがＦＢゲインＧ２として用いられる。

（１８）式の例のように、ＦＢゲインＧ２、すなわち、比ΔＰ／ΔＳＯＣの近似式としては、第２平均車速ａｖｅＶ２及び第２平均駆動トルクａｖｅＴｐ２を独立変数とする一次式を用いることができる。この近似式中の係数Ａ、Ｂ及びＣは、例えば、最小二乗法等の重回帰分析を用いて予め設定される。このような近似式の利用により、現時点に対する直前の走行履歴に基づく比ΔＰ／ΔＳＯＣ（すなわち、ＳＯＣに対する随伴変数ｐの感度）に応じた適切なＦＢゲインＧ２を取得（推定）できるようになる。

付け加えると、図１２に示すような関係を有する比ΔＰ／ΔＳＯＣの利用により、ＦＢゲインＧ２は、第２平均車速ａｖｅＶ２が高いほど大きくなるように設定されることになる。また、ＦＢゲインＧ２は、第２平均駆動トルクａｖｅＴｐ２が高いほど大きくなるように設定されることになる。

なお、ＦＢゲインＧ２の設定のために用いられる近似式は、上記の一次式の例に限られず、二次以上の高次の式であってもよい。また、ＦＢゲインＧ２の設定のために、このような近似式に代え、第２平均車速ａｖｅＶ２及び第２平均駆動トルクａｖｅＴｐ２に対する比ΔＰ／ΔＳＯＣの関係を定めたマップ（図示省略）が用いられてもよい。

４－２－３．制御ルーチン
本実施形態に係るパワートレーン制御に関する処理のメインルーチン（図示省略）は、ステップＳ３０２の処理に代えて次の図１９に示すサブルーチンの処理が実行される点において、実施の形態３の図１７に示すルーチンと相違している。

図１９は、本発明の実施の形態４に係るパワートレーン制御において、偏差ΔＳＯＣ_ｆｂ（ｋ）に応じて随伴変数ｐ（ｋ）を修正するために実行される処理のサブルーチンを示すフローチャートである。なお、本実施形態では、ステップＳ３００の処理（図１７参照）及び図１９に示すサブルーチンの処理は、「随伴変数修正処理８４ａ」に対応している。

図１９に示すサブルーチンでは、制御装置８０は、まず、ステップＳ４００において、偏差ΔＳＯＣ_ｆｂ（ｋ）を取得する。取得される偏差ΔＳＯＣ_ｆｂ（ｋ）は、ステップＳ３００の処理による算出値である。その後、処理は、ステップＳ４０２に進む。

ステップＳ４０２では、制御装置８０は、現時点に対する直前の車両走行履歴を取得する。具体的には、制御装置８０は、上述の所定期間（例えば、現時点を遡って３０秒）における第２平均車速ａｖｅＶ２及び第２平均駆動トルクａｖｅＴｐ２を取得する。これらは、より詳細には移動平均値である。すなわち、本サブルーチンの処理は、時間ステップΔｔ毎に繰り返し実行されるため、本ステップＳ４００による第２平均車速ａｖｅＶ２の算出に用いられる車速ｖの時系列データが時間ステップΔｔ毎に更新されていく。このことは、第２平均駆動トルクａｖｅＴｐ２についても同様である。

時間ステップΔｔ毎の第２平均車速ａｖｅＶ２及び第２平均駆動トルクａｖｅＴｐ２の値の取得手法の例は、実施の形態２において説明した第１平均車速ａｖｅＶ１及び第１平均駆動トルクａｖｅＴｐ１のものと同じである。付け加えると、ハイブリッド車両の新たなトリップ（今回のトリップ）の開始から上記所定期間の終了までの期間では、今回のトリップの開始後に取得された車速Ｖと要求駆動トルクＴｐのデータのみで上記所定期間中の第２平均車速ａｖｅＶ２及び第２平均駆動トルクａｖｅＴｐ２を満足に算出できない。このような期間中には、前回のトリップ中の終了直前の車速Ｖと要求駆動トルクＴｐのデータをも使用して第２平均車速ａｖｅＶ２及び第２平均駆動トルクａｖｅＴｐ２が算出されてもよい。しかしながら、前回のトリップ中のデータは使用せずに、今回のトリップの開始後に取得されたデータのみを使用して第２平均車速ａｖｅＶ２及び第２平均駆動トルクａｖｅＴｐ２が算出されてもよい。また、第２平均車速ａｖｅＶ２及び第２平均駆動トルクａｖｅＴｐ２の算出に用いられる上記所定期間は、現時点に対する直前の期間（すなわち、現時点から遡るように特定される期間）に必ずしも限られず、現時点に対して所定時間前の時点から遡るように特定される期間であってもよい。

ステップＳ４０２の後に、処理はステップＳ４０４に進む。ステップＳ４０４では、制御装置８０は、ＦＢゲインＧ２を決定する。具体的には、制御装置８０は、上記（１８）式を用いて、ステップＳ４０２の処理により取得した第２平均車速ａｖｅＶ２及び第２平均駆動トルクａｖｅＴｐ２に応じたＦＢゲインＧ２を算出する。その後、処理はステップＳ４０６に進む。

ステップＳ４０６では、制御装置８０は、随伴変数ｐ（ｋ）の補正値Ｃ（ｋ）を算出する。ここでいう補正値Ｃ（ｋ）は、上記（１７）式の右辺第２項に相当する。したがって、制御装置８０は、ステップＳ４００において取得した偏差ΔＳＯＣ_ｆｂ（ｋ）とステップＳ４０４において決定されたＦＢゲインＧ２との積を補正値Ｃ（ｋ）として算出する。その後、処理は、ステップＳ４０８に進む。

ステップＳ４０８では、制御装置８０は、上記（１７）式に従い、ステップＳ４０６の処理により算出した補正値Ｃ（ｋ）を用いて、ステップＳ１０２において取得された随伴変数ｐ（ｋ）を修正する。

４－３．効果
図２０（Ａ）及び図２０（Ｂ）は、本発明の実施の形態４の随伴変数修正処理８４ａとの対比のために参照する随伴変数修正処理（比較例１、２）を利用する場合の動作を説明するためのタイムチャートである。これらの比較例１及び２における随伴変数修正処理では、ＦＢゲインＧ２として固定値が用いられている。

まず、図２０（Ａ）に示す比較例１において、時点ｔ１～時点ｔ２までの期間は、この期間中の車両の走り方に対してＦＢゲインＧ２（固定値）が適している期間に相当する。このため、この期間では、ＳＯＣ制御中心値（狙いのＳＯＣ軌道）に対して実ＳＯＣ軌道が良好に追従している。

次に、時点ｔ２は、車両の走り方が変化した時点に相当する。時点ｔ２よりも後の期間では、変更後の走り方が継続されているものとする。そして、比較例１における時点ｔ２よりも後の期間では、走り方の変化に起因して、ＦＢゲインＧ２（固定値）は、理想的なＦＢゲインＧ２に対して小さ過ぎるものになっているものとする。その結果、比較例１では、随伴変数修正処理による随伴変数ｐの修正が満足に働かずに、図２０（Ａ）に示すように時間の経過とともにＳＯＣ制御中心値に対して実ＳＯＣが低くなる側（すなわち、放電側）に実ＳＯＣ軌道が徐々に変化してしまっている。そして、時点ｔ３では、強制的な充電の実施が必要なレベルにまで実ＳＯＣが低下している。また、強制的な充電が実際に行われた場合には、内燃機関２０の燃費が悪化する。

なお、比較例１は、ＦＢゲインＧ２が小さ過ぎるために実ＳＯＣが放電側に変化した状況を表している。しかしながら、状況によっては、ＦＢゲインＧ２が小さ過ぎるために実ＳＯＣが充電側に変化する場合もあり、そして、その結果として、強制的な放電の実施が必要とされる場合もある。また、強制的な放電が実際に行われた場合には、バッテリ４４からの不必要な放電がなされることになるため、内燃機関２０の燃費悪化に繋がる。

一方、図２０（Ｂ）に示す比較例２における時点ｔ２よりも後の期間では、走り方の変化に起因して、ＦＢゲインＧ２（固定値）は、理想的なＦＢゲインＧ２に対して大き過ぎるものになっているものとする。その結果、比較例２では、随伴変数修正処理による随伴変数ｐの修正が過剰に働き、図２０（Ｂ）に示すように実ＳＯＣ軌道がＳＯＣ制御中心値に対して大きくハンチングしている。その結果、バッテリ４４の不必要な放電／充電が生じてしまう。このことは、内燃機関２０の燃費悪化及びバッテリ４４の劣化の少なくとも一方に繋がり得る。

図２１は、本発明の実施の形態４の随伴変数修正処理８４ａを利用する場合の動作を説明するためのタイムチャートである。本実施形態の随伴変数修正処理８４ａによれば、上述したように、ＦＢゲインＧ２は、ハイブリッド車両の走行履歴（第２平均車速ａｖｅＶ２及び第２平均駆動トルクａｖｅＴｐ２）に基づいて設定（変更）される。より詳細には、図２１中に示す時点ｔ２において車両の走り方が変更された後には、第２平均車速ａｖｅＶ２及び第２平均駆動トルクａｖｅＴｐ２が時間ステップΔｔ毎に更新されていく。つまり、第２平均車速ａｖｅＶ２及び第２平均駆動トルクａｖｅＴｐ２のそれぞれの値に対して、走り方の変更後のデータが徐々に反映されていく。

そして、第２平均車速ａｖｅＶ２及び第２平均駆動トルクａｖｅＴｐ２の更新に伴い、ＦＢゲインＧ２である比ΔＰ／ΔＳＯＣも徐々に更新されていく。つまり、変更後の走り方が反映されるように比ΔＰ／ΔＳＯＣ（ＦＢゲインＧ２）の学習が進んでいく。その結果、ＦＢゲインＧ２が、変更後の走り方に適した値に徐々に近づいていく。なお、図２１に示す例では、変更後の走り方に適した値となるようにＦＢゲインＧ２が徐々に大きくなっている。

本実施形態の随伴変数修正処理８４ａによれば、上述のように、走り方が変更された場合には、変更後の走り方に適した値となるようにＦＢゲインＧ２が変更される。その結果、図２１に例示されるように、走り方が変化した場合であっても、ＳＯＣ制御中心値に対する実ＳＯＣ軌道のずれが過剰に大きくなることが抑制される。そして、ＦＢゲインＧ２が変更後の走り方に適した値に収束した後には、図２１に示されるように、実ＳＯＣ軌道をＳＯＣ制御中心値に適切に収束させられるようになる。

また、本実施形態の随伴変数修正処理８４ａによれば、ＦＢゲインＧ２は、第２平均車速ａｖｅＶ２が高いほど大きくなるように設定される。これにより、ＳＯＣに対する随伴変数ｐの感度を示す比ΔＰ／ΔＳＯＣが、第２平均車速ａｖｅＶ２が高いほど高くなるという特性を考慮して、ＦＢゲインＧ２を適切に設定できるようになる。さらに、ＦＢゲインＧ２は、第２平均駆動トルクａｖｅＴｐ２が高いほど大きくなるように設定される。これにより、上記感度を示す比ΔＰ／ΔＳＯＣが、第２平均駆動トルクａｖｅＴｐ２が高いほど高くなるという特性を考慮して、ＦＢゲインＧ２を適切に設定できるようになる。

５．他の実施の形態
５－１．ゲインＧ１、Ｇ２の設定に用いられる走行履歴の他の例
ここでは、ＦＢゲインＧ２を例に挙げて説明するが、以下の例は、補正ゲインＧ１の設定に用いられる走行履歴についても同様に適用可能である。

上述した実施の形態４においては、第２平均車速ａｖｅＶ２及び第２平均駆動トルクａｖｅＴｐ２の双方がＦＢゲインＧ２の設定のために用いられている。しかしながら、このような例に代え、第２平均車速ａｖｅＶ２及び第２平均駆動トルクａｖｅＴｐ２の何れか一方がＦＢゲインＧ２の設定のために用いられてもよい。また、第２平均駆動トルクａｖｅＴｐ２（Ｎｍ）に代え、所定期間における車両駆動出力の平均値である平均駆動出力（ｋＷ）が用いられてもよい。

さらに、ＦＢゲインＧ２の設定に用いられる走行履歴は、例えば、第２平均車速及び第２平均駆動トルクのうちの少なくとも一方に加え、平均車両加速度を含んでもよい。ここでいう平均車両加速度は、所定期間におけるハイブリッド車両の加速度の平均値であり、加速度は、正の加速度だけでなく、負の加速度（減速度）をも含む。平均車両加速度の算出に用いられる上記所定期間は、第２平均車速及び第２平均駆動トルクのうちの少なくとも一方の算出に用いられる所定期間と同じでもよいし、異なっていてもよい。そして、一例として、第２平均車速ａｖｅＶ２及び第２平均駆動トルクａｖｅＴｐ２に加えて平均車両加速度ａｖｅＡを使用する場合、ゲイン（ＦＢゲインＧ２）は、上述の（１８）式と同様の考えに基づき、例えば次の（１９）式に従って決定されてもよい。

また、ＦＢゲインＧ２の設定に用いられる走行履歴は、例えば、第２平均車速及び第２平均駆動トルクのうちの少なくとも一方に加え、平均車両加加速度を含んでもよく、或いは、第２平均車速及び第２平均駆動トルクのうちの少なくとも一方に加え、平均車両加速度と平均車両加加速度とを含んでもよい。ここでいう平均車両加加速度は、所定期間におけるハイブリッド車両の加加速度の平均値である。加加速度は、正の加加速度だけでなく、負の加加速度をも含む。平均車両加加速度の算出に用いられる上記所定期間は、第２平均車速及び第２平均駆動トルクのうちの少なくとも一方の算出に用いられる所定期間、並びに平均車両加速の算出に用いられる所定期間と同じでもよいし、異なっていてもよい。そして、一例として、第２平均車速ａｖｅＶ２及び第２平均駆動トルクａｖｅＴｐ２に加えて平均車両加速度ａｖｅＡ及び平均車両加加速度ａｖｅＪを使用する場合、ＦＢゲインＧ２は、上述の（１８）式と同様の考えに基づき、例えば次の（２０）式に従って決定されてもよい。

付け加えると、上述の平均車両加速度の算出に用いられる加速度は、加速度センサ（Ｇセンサ）を用いて取得されてもよい。また、平均車両加加速度の算出に用いられる加加速度は、加速度センサの出力を時間で微分することによって取得されてもよい、或いは、加加速度センサを用いて取得されてもよい。また、これらの加速度及び加加速度は、例えば以下のように、アクセルペダルの操作情報に基づく推定値として取得されてもよい。

具体的には、加速度は、例えば、アクセルペダルの踏み込み量（又はその踏力）と加速度との関係を定める関係情報（例えば、関係式又はマップ）を用いて、踏み込み量（又は踏力）に応じた値として推定されてもよい。また、加加速度は、例えば、アクセルペダルの踏み込み速度と加加速度との関係を定める関係情報（例えば、関係式又はマップ）を用いて、踏み込み速度に応じた値として推定されてもよい。

５－２．パワートレーンシステムの他の例
上述した実施の形態１～４においては、図１に示す動力分割機構３４を利用する動力分割方式のハイブリッドシステムに相当するパワートレーンシステム１０について説明した。しかしながら、本発明の対象となる「パワートレーンシステム」は、ハイブリッド車両を駆動可能な内燃機関及び電動機と、当該電動機に電力を供給するバッテリを備えていれば、上記の例に限られない。すなわち、本発明に係るパワートレーンシステムは、例えば、車両を駆動する内燃機関と電動機とが並列に連結されたパラレル方式のハイブリッドシステムとして構成されてもよい。

５－３．パワートレーンシステムの制御入力値の他の例
上述した実施の形態１～４においては、制御入力値ｕとして、エンジン回転数Ｎｅ及びエンジントルクＴｅが用いられた。しかしながら、本発明に係る最適化対象の「制御入力値」は、パワートレーンシステムの構成次第では、上記の例に限られない。すなわち、例えば有段式の自動変速機と組み合わされた内燃機関を備えるパワートレーンシステムでは、エンジン回転数Ｎｅに代えて自動変速機のギヤ段がエンジントルクＴｅとともに制御入力値として用いられてもよい。これは、選択されているギヤ段が分かると、車速Ｖに基づいてエンジン回転数Ｎｅが分かるためである。なお、有段式の自動変速機を備える例では、ギヤ段が変更されると、エンジン回転数Ｎｅが非連続に（ステップ的に）変化する。このため、最適な制御入力値の探索の際に、ギヤ段毎にハミルトニアンを算出し、算出された複数のハミルトニアンのうちで最小となるハミルトニアンに対応するギヤ段が選択されるように自動変速機が制御されてもよい。

以上説明した各実施の形態に記載の例及び他の各変形例は、明示した組み合わせ以外にも可能な範囲内で適宜組み合わせてもよいし、また、本発明の趣旨を逸脱しない範囲で種々変形してもよい。

１０ パワートレーンシステム
２０ 内燃機関
２８ クランク角センサ
３０ モータジェネレータ（ＭＧ１）
３２ モータジェネレータ（ＭＧ２）
３４ 動力分割機構
３８ 車輪
４０ インバータ
４４ バッテリ
４５ 外部電源
４６ 電気製品
５０、６２、７０、８０ 制御装置
５２ センサ類

Claims (14)

1. ハイブリッド車両を駆動する内燃機関と、
   前記ハイブリッド車両を駆動する電動機と、
   前記電動機に電力を供給するバッテリと、
   前記内燃機関及び前記電動機を制御するパワートレーン制御を行う制御装置と、
   を備えるパワートレーンシステムであって、
   前記ハイブリッド車両は、１又は複数の外部機器との間で前記バッテリの充電及び放電のうちの少なくとも一方を実行可能に構成されており、
   前記制御装置は、前記パワートレーン制御が行われる制御期間中に
   前記バッテリの充電率ＳＯＣの動的挙動を制約条件として、前記ハイブリッド車両の速度及び駆動トルクに対して前記制御期間中の前記内燃機関の燃料消費量ｍｆを最小にする最適化問題を解くことにより、前記バッテリの充放電量及び前記パワートレーンシステムのシステム動作点を特定する制御入力値を求める制御入力決定処理と、
   前記制御入力決定処理により求めた前記制御入力値を前記パワートレーンシステムに与えることにより、前記内燃機関と前記電動機とを制御するシステム制御処理と、
   を実行し、
   前記制御入力決定処理は、
   前記最適化問題の随伴変数ｐを時間ステップ毎に更新する随伴変数決定処理と、
   前記随伴変数決定処理により決定される前記随伴変数ｐを用いて、以下の式により定義されるハミルトニアンＨを最小とする前記制御入力値を前記時間ステップ毎に探索して算出する制御入力算出処理と、を含み、
   前記随伴変数決定処理は、前記随伴変数ｐの初期値の基本値と外部充放電補正値との和を前記初期値とする初期値決定処理を含み、
   前記基本値は、前回の前記制御期間中の前記随伴変数ｐの最終値又は平均値であり、
   前記外部充放電補正値は、今回の前記制御期間の開始時の前記充電率ＳＯＣから前回の前記制御期間の終了時の前記充電率ＳＯＣを引いて得られる外部充放電量に基づいて決定され、
   前記初期値決定処理による前記外部充放電補正値は、前記外部充放電量と第１ゲインとの積であり、
   前記第１ゲインは、前記ハイブリッド車両の走行履歴に基づいて設定される
   ことを特徴とするパワートレーンシステム。
   

   

   ただし、上記の式の右辺第１項は、前記燃料消費量ｍｆの時間変化率ｄｍｆ／ｄｔであり、右辺第２項は、前記随伴変数ｐと前記充電率ＳＯＣの時間変化率ｄＳＯＣ／ｄｔとの積である。
2. 前記初期値決定処理は、前記外部充放電量が正の場合には、前記外部充放電量がゼロの場合と比べて、前記初期値を大きくする
   ことを特徴とする請求項１に記載のパワートレーンシステム。
3. 前記初期値決定処理は、前記外部充放電量が負の場合には、前記外部充放電量がゼロの場合と比べて、前記初期値を小さくする
   ことを特徴とする請求項１又は２に記載のパワートレーンシステム。
4. 前記走行履歴は、前回の前記制御期間中の前記速度及び前記駆動トルクのそれぞれの平均値である第１平均車速及び第１平均駆動トルクのうちの少なくとも一方を含む
   ことを特徴とする請求項１～３の何れか１つに記載のパワートレーンシステム。
5. 前記走行履歴が前記第１平均車速を含む場合において、前記第１ゲインは、前記第１平均車速が高いほど大きい
   ことを特徴とする請求項４に記載のパワートレーンシステム。
6. 前記走行履歴が前記第１平均駆動トルクを含む場合において、前記第１ゲインは、前記第１平均駆動トルクが高いほど大きい
   ことを特徴とする請求項４又は５に記載のパワートレーンシステム。
7. 前記随伴変数決定処理は、前記充電率ＳＯＣの目標充電率と実充電率との差に基づいて前記随伴変数ｐを修正する随伴変数修正処理を含む
   ことを特徴とする請求項１～６の何れか１つに記載のパワートレーンシステム。
8. ハイブリッド車両を駆動する内燃機関と、
   前記ハイブリッド車両を駆動する電動機と、
   前記電動機に電力を供給するバッテリと、
   前記内燃機関及び前記電動機を制御するパワートレーン制御を行う制御装置と、
   を備えるパワートレーンシステムであって、
   前記ハイブリッド車両は、１又は複数の外部機器との間で前記バッテリの充電及び放電のうちの少なくとも一方を実行可能に構成されており、
   前記制御装置は、前記パワートレーン制御が行われる制御期間中に
   前記バッテリの充電率ＳＯＣの動的挙動を制約条件として、前記ハイブリッド車両の速度及び駆動トルクに対して前記制御期間中の前記内燃機関の燃料消費量ｍｆを最小にする最適化問題を解くことにより、前記バッテリの充放電量及び前記パワートレーンシステムのシステム動作点を特定する制御入力値を求める制御入力決定処理と、
   前記制御入力決定処理により求めた前記制御入力値を前記パワートレーンシステムに与えることにより、前記内燃機関と前記電動機とを制御するシステム制御処理と、
   を実行し、
   前記制御入力決定処理は、
   前記最適化問題の随伴変数ｐを時間ステップ毎に更新する随伴変数決定処理と、
   前記随伴変数決定処理により決定される前記随伴変数ｐを用いて、以下の式により定義されるハミルトニアンＨを最小とする前記制御入力値を前記時間ステップ毎に探索して算出する制御入力算出処理と、を含み、
   前記随伴変数決定処理は、前記随伴変数ｐの初期値の基本値と外部充放電補正値との和を前記初期値とする初期値決定処理を含み、
   前記基本値は、前回の前記制御期間中の前記随伴変数ｐの最終値又は平均値であり、
   前記外部充放電補正値は、今回の前記制御期間の開始時の前記充電率ＳＯＣから前回の前記制御期間の終了時の前記充電率ＳＯＣを引いて得られる外部充放電量に基づいて決定され、
   前記随伴変数決定処理は、前記充電率ＳＯＣの目標充電率と実充電率との差に基づいて前記随伴変数ｐを修正する随伴変数修正処理を含む
   ことを特徴とするパワートレーンシステム。
   

   

   ただし、上記の式の右辺第１項は、前記燃料消費量ｍｆの時間変化率ｄｍｆ／ｄｔであり、右辺第２項は、前記随伴変数ｐと前記充電率ＳＯＣの時間変化率ｄＳＯＣ／ｄｔとの積である。
9. 前記随伴変数修正処理は、前記目標充電率よりも前記実充電率が低い場合には、前記随伴変数ｐを小さくする
   ことを特徴とする請求項８に記載のパワートレーンシステム。
10. 前記随伴変数修正処理は、前記目標充電率よりも前記実充電率が高い場合には、前記随伴変数ｐを大きくする
    ことを特徴とする請求項８又は９に記載のパワートレーンシステム。
11. 前記随伴変数修正処理による前記随伴変数ｐの補正値は、前記差と第２ゲインとの積であり、
    前記第２ゲインは、前記ハイブリッド車両の走行履歴に基づいて設定される
    ことを特徴とする請求項８～１０の何れか１つに記載のパワートレーンシステム。
12. 前記走行履歴は、第１所定期間における前記速度の平均値である第２平均車速、及び第２所定期間における前記駆動トルクの平均値である第２平均駆動トルクのうちの少なくとも一方を含む
    ことを特徴とする請求項１１に記載のパワートレーンシステム。
13. 前記走行履歴が前記第２平均車速を含む場合において、前記第２ゲインは、前記第２平均車速が高いほど大きい
    ことを特徴とする請求項１２に記載のパワートレーンシステム。
14. 前記走行履歴が前記第２平均駆動トルクを含む場合において、前記第２ゲインは、前記第２平均駆動トルクが高いほど大きい
    ことを特徴とする請求項１２又は１３に記載のパワートレーンシステム。

Images