(57)【要約】【課題】システム同定の精度を向上させうるシステム同定装置を提供する。【解決手段】線形要素２２と非線形要素２４とを含むシステム２を同定するシステム同定装置１であって、システム２に対して入力される、第１入力と第１入力との間に所定の関係を有する第２入力とを生成する生成部１４と、システム２から、第１入力に対応する第１出力と、第２入力に対応する第２出力とを取得する取得部１６と、第１出力と第２出力との差分信号に基づいて、システム２の線形要素２２を同定する制御部１０とを備えるシステム同定装置１。【選択図】図１
2018190288.tif 000002
【特許請求の範囲】
【請求項１】
線形要素と非線形要素とを含むシステムを同定するシステム同定装置であって、
前記システムに対して入力される、第１入力と、前記第１入力との間に所定の関係を有する第２入力とを生成する生成部と、
前記システムから、前記第１入力に対応する第１出力と、前記第２入力に対応する第２出力とを取得する取得部と、
前記第１出力と前記第２出力との差分信号に基づいて、前記システムの前記線形要素を同定する制御部と
を備えるシステム同定装置。
【請求項２】
前記非線形要素は、前記システムに対する入力に関する斉次多項式で表される、請求項１に記載のシステム同定装置。
【請求項３】
前記生成部は、前記差分信号において前記非線形要素に係る出力が消去されるように、前記第１入力と前記第２入力とを生成する、請求項１又は２に記載のシステム同定装置。
【請求項４】
前記制御部は、前記線形要素の同定結果に基づいて、前記非線形要素を同定する、請求項１乃至３のいずれか一項に記載のシステム同定装置。
【請求項５】
線形要素と非線形要素とを含むシステムを同定するシステム同定方法であって、
前記システムに対して入力される、第１入力と前記第１入力との間に所定の関係を有する第２入力とを生成するステップと、
前記システムから、前記第１入力に対応する第１出力と、前記第２入力に対応する第２出力とを取得するステップと、
前記第１出力と前記第２出力との差分信号に基づいて、前記システムの前記線形要素を同定するステップと
を含むシステム同定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、システム同定装置及びシステム同定方法に関する。
【背景技術】
【０００２】
線形要素と非線形要素とを含むシステムを同定する装置が知られている（例えば特許文献１参照）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１４−８６３１３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
システムに含まれる線形要素と非線形要素とが同時に同定される場合、それぞれの要素が個別に同定される場合と比較して、互いの要素が影響を及ぼしうる。線形要素と非線形要素とが互いに影響を及ぼすことによって、システム同定の精度が向上されにくい。
【０００５】
かかる事情に鑑みてなされた本発明の目的は、システム同定の精度を向上させうるシステム同定装置及びシステム同定方法を提供することにある。
【課題を解決するための手段】
【０００６】
上記課題を解決するために、第１の観点に係るシステム同定装置は、
線形要素と非線形要素とを含むシステムを同定するシステム同定装置であって、
前記システムに対して入力される、第１入力と前記第１入力との間に所定の関係を有する第２入力とを生成する生成部と、
前記システムから、前記第１入力に対応する第１出力と、前記第２入力に対応する第２出力とを取得する取得部と、
前記第１出力と前記第２出力との差分信号に基づいて、前記システムの前記線形要素を同定する制御部と
を備える。
【０００７】
上記課題を解決するために、第２の観点に係るシステム同定方法は、
線形要素と非線形要素とを含むシステムを同定するシステム同定方法であって、
前記システムに対して入力される、第１入力と前記第１入力との間に所定の関係を有する第２入力とを生成するステップと、
前記システムから、前記第１入力に対応する第１出力と、前記第２入力に対応する第２出力とを取得するステップと、
前記第１出力と前記第２出力との差分信号に基づいて、前記システムの前記線形要素を同定するステップと
を含む。
【発明の効果】
【０００８】
第１の観点に係るシステム同定装置によれば、システム同定の精度が向上されうる。
【０００９】
第２の観点に係るシステム同定方法によれば、システム同定の精度が向上されうる。
【図面の簡単な説明】
【００１０】
【図１】一実施形態に係るシステム同定装置の構成例を示すブロック図である。
【図２】入力の例を示すグラフである。（Ａ）第１入力（Ｂ）第２入力。
【図３】出力の例を示すグラフである。
【図４】一実施形態に係るシステム同定方法の例を示すフローチャートである。
【図５】バッテリのモデルを近似する等価回路の一例を示す回路図である。
【図６】バッテリのモデルを近似する等価回路の一例を示す回路図である。
【図７】バッテリのＳＯＣ−ＯＣＶ特性の一例を示すグラフである。
【図８】システム同定の例を示すグラフである。（Ａ）入力（Ｂ）線形要素に係る出力（Ｃ）非線形要素に係る出力。
【発明を実施するための形態】
【００１１】
図１に示されるように、一実施形態に係るシステム同定装置１は、線形要素２２と非線形要素２４とを含む対象システム２に入力を与え、対象システム２から出力を取得する。対象システム２の入力及び出力はそれぞれ、ｕ（ｔ）及びｙ（ｔ）で表されるものとする。システム同定装置１は、制御部１０と、記憶部１２と、生成部１４と、取得部１６とを備える。
【００１２】
制御部１０は、システム同定装置１の各構成部と接続され、各構成部を制御する。制御部１０は、例えばプロセッサとして構成されてよい。制御部１０は、記憶部１２に、各構成部を動作させるためのプログラム等を格納してよい。
【００１３】
記憶部１２は、例えば半導体メモリ等で構成されてよい。記憶部１２は、制御部１０のワークメモリとして機能してよい。記憶部１２は、制御部１０に含まれてよい。
【００１４】
生成部１４は、対象システム２に与える入力を生成する。入力は、信号として表されるものとする。入力を表す信号は、入力信号ともいう。生成部１４は、例えばプロセッサとして構成されてよく、制御部１０に含まれてもよい。
【００１５】
取得部１６は、対象システム２からの出力を取得する。出力は、信号として表されるものとする。出力を表す信号は、出力信号ともいう。取得部１６は、例えばプロセッサとして構成されてよく、制御部１０に含まれてもよい。
【００１６】
対象システム２は、入力信号に応じた出力信号を生成する。ｕ（ｔ）で表される入力信号は、システム同定装置１から入力され、線形要素２２及び非線形要素２４のそれぞれに入力される。線形要素２２及び非線形要素２４はそれぞれ、入力信号に応じた信号を出力する。線形要素２２及び非線形要素２４のそれぞれから出力された信号は、加算され、ｙ（ｔ）で表される出力信号となり、システム同定装置１に出力される。
【００１７】
対象システム２において、ｕ（ｔ）に対する線形要素２２の応答を表すｙＬ（ｔ）は、以下の式（１）で表されるものとする。
【数１】
2018190288.tif 000003
Ｇ（ｐ）は、伝達関数である。ｐは、微分オペレータである。
【００１８】
対象システム２において、ｕ（ｔ）に対する非線形要素２４の応答を表すｙｎ（ｔ）は、以下の式（２）及び式（３）で表されるものとする。
【数２】
2018190288.tif 000004
ｘ（ｔ）は、対象システム２の状態変数である。ｆ（・）及びｇ（・）は、所定の関数である。
【００１９】
対象システム２の出力信号を表すｙ（ｔ）は、以下の式（４）で表されるものとする。
【数３】
2018190288.tif 000005
【００２０】
本実施形態に係るシステム同定装置１は、ｕ（ｔ）とｙ（ｔ）とに基づいて、Ｇ（ｐ）で表される伝達関数を推定しうる。
【００２１】
非線形要素２４は、入力に対して斉次なシステムであるものとする。式（２）及び式（３）によって表されるシステムが入力に対して斉次なシステムである場合、以下の式（５）が成り立つものとする。
【数４】
2018190288.tif 000006
α及びｄは、実数である。ｄは、斉次性の次数ともいう。
【００２２】
入力に対して斉次なシステムは、入力に対応する変数の次数が同じ次数となっている多項式で表されてよい。入力に対応する変数の次数が同じ次数となっている多項式は、斉次多項式ともいう。例えば、システムを表す多項式の各項にｕが含まれ、ｕ２等の他の次数の項が含まれない場合、システムは、ｄ＝１の入力斉次なシステムであるといえる。システムを表す多項式の各項にｕ２が含まれ、ｕ又はｕ３等の他の次数の項が含まれない場合、システムは、ｄ＝２の入力に対して斉次なシステムであるといえる。
【００２３】
入力に対して斉次なシステムの応答は、２つの異なる入力に対する応答の差分をとることによって、消去されうる。例えば、ｕ１（ｔ）で表される第１入力と、ｕ２（ｔ）で表される第２入力とが対象システム２に入力されると仮定する。第２入力は、第１入力との間に所定の関係を有するものとする。対象システム２にｕ１（ｔ）が入力されるときの対象システム２からの出力である第１出力は、ｙ１（ｔ）で表されるものとする。対象システム２にｕ２（ｔ）が入力されるときの対象システム２からの出力である第２出力は、ｙ２（ｔ）で表されるものとする。
【００２４】
ｕ１（ｔ）とｕ２（ｔ）との間に成り立つ所定の関係は、以下の式（６）で表される関係であるものとする。
【数５】
2018190288.tif 000007
【００２５】
式（６）が成り立つ場合、ｕ１（ｔ）及びｕ２（ｔ）はそれぞれ、図２（Ａ）及び図２（Ｂ）に例示されるようなグラフとなりうる。図２（Ａ）において、ｕ１（ｔ）は、時刻（ｔ）が０から１までの範囲で定義され、０から１までの値をとる。図２（Ｂ）において、ｕ２（ｔ）は、時刻（ｔ）が０から１／αｄまでの範囲で定義され、０からαｄまでの値をとる。つまり、ｕ２（ｔ）は、ｕ１（ｔ）に対して、１／αｄ倍の時間軸と、αｄ倍の値とを有する。
【００２６】
対象システム２にｕ１（ｔ）が入力されるときの非線形要素２４に含まれる状態変数は、ｘ１（ｔ）で表されるものとする。対象システム２にｕ２（ｔ）が入力されるときの非線形要素２４に含まれる状態変数は、ｘ２（ｔ）で表されるものとする。ｘ１（ｔ）及びｘ２（ｔ）について、以下の式（７）が成り立つものとする。
【数６】
2018190288.tif 000008
【００２７】
非線形要素２４にｕ１（ｔ）が入力されるときの非線形要素２４からの出力は、ｙｎ，１（ｔ）で表されるものとする。非線形要素２４にｕ２（ｔ）が入力されるときの非線形要素２４からの出力は、ｙｎ，２（ｔ）で表されるものとする。式（６）及び式（７）が成り立つ場合、ｙｎ，１（ｔ）とｙｎ，２（ｔ）との間には、以下の式（８）で表される関係が成り立つ。
【数７】
2018190288.tif 000009
【００２８】
線形要素２２にｕ１（ｔ）が入力されるときの線形要素２２からの出力は、ｙＬ，１（ｔ）で表されるものとする。ｙ１（ｔ）は、ｙＬ，１（ｔ）とｙｎ，１（ｔ）との和である。線形要素２２にｕ２（ｔ）が入力されるときの線形要素２２からの出力は、ｙＬ，２（ｔ）で表されるものとする。ｙ２（ｔ）は、ｙＬ，２（ｔ）とｙｎ，２（ｔ）との和である。
【００２９】
以下の式（９）で表される差分信号が計算されうる。
【数８】
2018190288.tif 000010
【００３０】
図３に示されるように、ｙ２（ｔ）の時間軸は、ｙ１（ｔ）の時間軸に対して１／αｄ倍となっている。式（９）は、ｙ２（ｔ）の時間軸がｙ１（ｔ）の時間軸に合わせるように伸縮されることを示している。
【００３１】
式（８）が成り立つ場合、式（９）からｙｎ，１（ｔ）の項とｙｎ，２（ｔ）の項とが消去される。つまり、入力に対して斉次なシステムである非線形要素２４に係る出力が消去される。非線形要素２４に係る出力が消去される結果、以下の式（１０）で表されるように、線形要素２２の出力だけで表される関係が導かれる。
【数９】
2018190288.tif 000011
【００３２】
式（１０）は、以下の式（１１）のようにも表されうる。
【数１０】
2018190288.tif 000012
Ｇ’（ｐ）は、Ｇ（ｐ）のｐとして、αｄｐを代入したものである。つまり、Ｇ’（ｐ）＝Ｇ（αｄｐ）である。
【００３３】
システム同定装置１は、式（６）に基づいて第１入力及び第２入力を生成し、各入力に応じた出力の差分信号を計算しうる。システム同定装置１は、差分信号と第１入力とを式（１１）に適用することによって、Ｇ（ｐ）で表される伝達関数を計算しうる。
【００３４】
システム同定装置１は、例えば図４のフローチャートに示される処理によって、対象システム２を同定しうる。
【００３５】
システム同定装置１は、生成部１４によって第１入力（ｕ１（ｔ））を生成し、第１入力を対象システム２に印加する（ステップＳ１）。システム同定装置１は、制御部１０によって、第１入力を記憶部１２に格納してよい。
【００３６】
システム同定装置１は、取得部１６によって、第１入力に応じた対象システム２からの第１出力（ｙ１（ｔ））を取得する（ステップＳ２）。取得部１６は、第１出力を、所定のサンプリング時刻で、Ｎ回にわけて取得するものとする。Ｎは、適宜決められうる数である。所定のサンプリング時刻は、等間隔で設定されてよいし、異なる間隔で設定されてよい。取得部１６は、第１出力をＮ回にわけて取得することによって、Ｎ個のデータからなるデータセットを取得する。第１出力のデータセット（Ｙ１）は、Ｙ１＝ｙ１（ｔｋ）で表されるものとする。ｋは、Ｎ以下の自然数である。ｔｋは、ｋ回目にデータを取得した時刻である。システム同定装置１は、制御部１０によって、第１出力のデータセットを記憶部１２に格納してよい。
【００３７】
システム同定装置１は、ステップＳ１の処理とステップＳ２の処理とを並行して実行し、Ｎ個のデータを取得するまで処理を繰り返す。システム同定装置１は、Ｎ個のデータを取得後、ステップＳ３に進む。
【００３８】
システム同定装置１は、生成部１４によって第２出力（ｕ２（ｔ））を生成し、第２入力を対象システム２に印加する（ステップＳ３）。生成部１４は、ｕ１（ｔ）とｕ２（ｔ）との関係が式（６）を満たすように、ステップＳ１で生成した第１入力に基づいて、第２入力を生成する。生成部１４は、非線形要素２４に係る出力が消去されるように、第１入力と第２入力とを生成するともいえる。第１出力と第２出力との間に成り立つ所定の関係は、非線形要素２４に係る出力が消去されるように定められる関係であるともいえる。システム同定装置１は、制御部１０によって、第２入力を記憶部１２に格納してよい。
【００３９】
システム同定装置１は、取得部１６によって、第２入力に応じた対象システム２からの第２出力（ｙ２（ｔ））を取得する（ステップＳ４）。取得部１６は、第２出力を、所定のサンプリング時刻で、Ｎ回にわけて取得するものとする。Ｎは、第１出力を取得した回数と同じ数とされる。所定のサンプリング時刻は、第１出力を取得したサンプリング時刻に基づいて設定される。第１出力の各データを取得したサンプリング時刻がｔｋである場合、第２出力を取得するサンプリング時刻は、ｔｋ／αｄに設定される。第２出力のデータセット（Ｙ２）は、Ｙ２＝ｙ２（ｔｋ／αｄ）で表されるものとする。ｋは、Ｎ以下の自然数である。システム同定装置１は、制御部１０によって、第２出力のデータセットを記憶部１２に格納してよい。
【００４０】
システム同定装置１は、ステップＳ３の処理とステップＳ４の処理とを並行して実行し、Ｎ個のデータを取得するまで処理を繰り返す。システム同定装置１は、Ｎ個のデータを取得後、ステップＳ５に進む。
【００４１】
システム同定装置１は、制御部１０によって、第１出力のデータセット（Ｙ１）と第２出力のデータセット（Ｙ２）との差分信号を算出する（ステップＳ５）。制御部１０は、１〜Ｎの各ｋについて、ｙ１（ｔｋ）とｙ２（ｔｋ／αｄ）との差を計算し、ｅ（ｔｋ）で表される差分信号を算出する。つまり、制御部１０は、ｅ（ｔｋ）＝ｙ１（ｔｋ）−ｙ２（ｔｋ／αｄ）を計算して、差分信号を算出する。
【００４２】
システム同定装置１は、制御部１０によって、差分信号と第１出力とに基づいて伝達関数を算出し、対象システム２の線形要素２２を同定する（ステップＳ６）。伝達関数は、パラメータ（θ）によってパラメトライズされる場合、Ｇ＾（ｐ，θ）で表されるものとする。この場合、制御部１０は、以下の式（１２）及び式（１３）によって示される評価関数（Ｊ（θ））を最小化するパラメータ（θ）を同定することによって、伝達関数を同定しうる。
【数１１】
2018190288.tif 000013
【００４３】
システム同定装置１は、制御部１０によって、ステップＳ６の処理で同定した線形要素２２と、第１出力とに基づいて、非線形要素２４を同定する（ステップＳ７）。制御部１０は、第１出力のデータセットから、同定した線形要素２２に対応する応答成分を消去したデータセットを生成し、当該データセットと第１出力とに基づいて非線形要素２４を同定してよい。言い換えれば、制御部１０は、ｙ１（ｔ）−Ｇ＾（ｐ）ｕ１（ｔ）の計算結果とｕ１（ｔ）とに基づいて非線形要素２４を同定してよい。
【００４４】
制御部１０は、線形要素２２の同定結果と第２出力とに基づいて、非線形要素２４を同定してもよい。制御部１０は、第２出力のデータセットから、同定した線形要素２２に係る出力を消去したデータセットを生成してよい。制御部１０は、線形要素２２に係る出力を消去したデータセットと第２出力とに基づいて非線形要素２４を同定してよい。言い換えれば、制御部１０は、ｙ２（ｔ）−Ｇ＾（ｐ）ｕ２（ｔ）の計算結果とｕ２（ｔ）とに基づいて非線形要素２４を同定してよい。
【００４５】
システム同定装置１は、ステップＳ７の処理の後、図４のフローチャートの処理を終了してよいし、ステップＳ１に戻って処理を繰り返してもよい。
【００４６】
本実施形態に係るシステム同定装置１は、複数の入力信号を入力し、各入力信号に応じた出力信号に基づいて、システムを同定しうる。複数の入力信号が所定の関係を有することによって、システム同定の精度が向上されうる。対象システム２の非線形要素２４が入力に関する斉次多項式で表される場合、システム同定装置１は、対象システム２の出力から非線形要素２４に係る出力を消去しやすくなる。対象システム２の出力から非線形要素２４に係る出力が消去されるように所定の関係が定められる場合、システム同定装置１は、対象システム２の出力から非線形要素２４に係る出力を消去しやすくなる。
【００４７】
システム同定において、非線形要素２４を含むシステムの同定は、線形要素２２のみを含むシステムの同定よりも複雑となることがある。非線形要素２４を含むシステムの同定精度は、線形要素２２のみを含むシステムの同定精度よりも低くなることがある。システム同定装置１は、非線形要素２４に係る出力が消去された出力信号に基づいて線形要素２２を同定することによって、非線形要素２４を含む出力信号に基づいて同定する場合と比較して、線形要素２２の同定精度を向上しうる。システム同定装置１は、精度よく同定した線形要素２２に基づいて非線形要素２４を同定することによって、非線形要素２４の同定精度を向上しうる。
【００４８】
システムを同定する際に、システムの状態変数等の値が初期値として仮定されうる。状態変数の真値と仮定された初期値との間の差異が大きい場合、非線形要素２４に係る出力の誤差が大きくなりうる。非線形要素２４に係る出力の誤差が大きい場合、非線形要素２４を含むシステムの同定の精度は、さらに低くなりうる。システム同定装置１は、状態変数の真値と仮定された初期値との間の差異が大きい場合でも、非線形要素２４に係る出力を消去することによって、線形要素２２の同定精度を向上しうる。
【００４９】
［バッテリのパラメータ推定］
本実施形態に係るシステム同定装置１は、対象システム２がバッテリである場合、バッテリのパラメータ推定装置として機能しうる。以下、対象システム２がバッテリであるものとする。入力信号は、バッテリの充放電電流に対応する。この場合、生成部１４は、電流を出力可能な電源装置として構成されてよいし、電源装置の電流を制御するように構成されてよい。出力信号は、バッテリの端子電圧に対応する。この場合、取得部１６は、電圧センサとして構成されてよいし、電圧センサからバッテリの端子電圧の測定値を取得するように構成されてよい。
【００５０】
バッテリの内部状態は、バッテリのヒステリシス特性と、バッテリの開回路電圧と、バッテリに対する充放電電流に応じて発生する過電圧とを含むモデルによって表されうる。開回路電圧は、ＯＣＶ（Open Circuit Voltage）ともいう。ＯＣＶは、バッテリの電気化学的平衡状態における電極間の電位差である。ＯＣＶは、バッテリに充放電電流が流れない場合におけるバッテリの端子電圧に対応する。過電圧は、バッテリの内部インピーダンスで生じる電圧降下の大きさに対応する。
【００５１】
バッテリの内部状態を表すモデルは、図５及び図６で示されるようなバッテリ等価回路で近似されうる。バッテリ等価回路で近似されたモデルは、バッテリモデルともいう。図５において、バッテリの内部インピーダンスは、カウエル型回路で近似されている。図６において、バッテリの内部インピーダンスは、フォスタ型回路で近似されている。
【００５２】
図５及び図６において、バッテリ等価回路に入力される充放電電流は、ｕ（ｔ）として示される。ｕ（ｔ）が付された矢印は、バッテリを充電する電流の向きを表す。バッテリを充電する電流が流れる場合、ｕ（ｔ）は正の値となるものとする。バッテリから放電電流が流れる場合、ｕ（ｔ）は負の値となるものとする。ｕ（ｔ）に対応して内部インピーダンスで生じる過電圧は、η（ｔ）と表される。
【００５３】
図５及び図６において、バッテリ等価回路の端子電圧は、ｙ（ｔ）として示される。ｙ（ｔ）が付された矢印の先端側の端子は、バッテリの正極端子に対応するものとする。ｙ（ｔ）は、バッテリのＯＣＶと過電圧との和で表される。
【００５４】
バッテリのＯＣＶは、電圧源２０１で表される。電圧源２０１が出力する電圧は、ＯＣＶが時刻（ｔ）の関数であることを示すようにＯＣＶ（ｔ）で表される。ＯＣＶは、バッテリの充電率の関数として表されうる。バッテリの充電率は、ＳＯＣ（State Of Charge）ともいう。ＳＯＣとＯＣＶとの間の関係は、ＳＯＣ−ＯＣＶ特性といわれる。ＳＯＣ−ＯＣＶ特性は、例えば図７に示されるグラフで表されうる。図７の横軸及び縦軸はそれぞれ、ＳＯＣ及びＯＣＶを示す。ＳＯＣ−ＯＣＶ特性は、予め実験等によって求められうる。ＳＯＣがｘｃｈ（ｔ）と表される場合、図７に示されるような特性は、静的な非線形関数ｆＯＣＶ（・）を用いて、式（１４）のように表されうる。
【数１２】
2018190288.tif 000014
【００５５】
ｘｃｈ（ｔ）は、バッテリの満充電容量（Full Charge Capacity）を表すＣを用いて、式（１５）のように表される。
【数１３】
2018190288.tif 000015
【００５６】
ＳＯＣ−ＯＣＶ特性は、ヒステリシス特性を有することがある。ヒステリシス特性を有するＳＯＣ−ＯＣＶ特性は、充電時の特性と放電時の特性とが異なる。充電時のＯＣＶと放電時のＯＣＶとの差は、ヒステリシス電圧として表されうる。バッテリがヒステリシス特性を有する場合、電圧源２０１が出力する電圧は、ヒステリシス電圧を表すｈ（ｔ）とＯＣＶ（ｔ）との和となる。
【００５７】
バッテリがヒステリシス特性を有する場合、バッテリ等価回路の端子電圧を表すｙ（ｔ）は、式（１６）のように表される。
【数１４】
2018190288.tif 000016
【００５８】
ヒステリシス電圧は、例えばPlettによるヒステリシスモデルによれば、式（１７）のように表される。
【数１５】
2018190288.tif 000017
γは、ヒステリシス電圧の減衰速度を指定する正数である。ｍは、ヒステリシス電圧の最大値を表す。Plettによるヒステリシスモデルについては、例えば、以下の文献が参照されうる。
G. L. Plett: “Extended Kalman filtering for battery management systems of LiPB-based HEV battery packs Part 2. Modeling and identification”, Journal of Power Sources 134 (2004) 262-276
【００５９】
過電圧は、バッテリの内部インピーダンスとバッテリの充放電電流とで表される。内部インピーダンスが線形システムの伝達関数としてＧη（ｐ）で表される場合、過電圧は、以下の式（１８）のように表される。
【数１６】
2018190288.tif 000018
【００６０】
内部インピーダンスが図５に示されるカウエル型回路で近似される場合、Ｇη（ｐ）は、以下の式（１９）、式（２０）及び式（２１）で表される。
【数１７】
2018190288.tif 000019
ここで、ｎは自然数である。Ｒ０、Ｒｄ及びτｄは、バッテリのパラメータである。
【００６１】
図５及び図６に示されるバッテリ等価回路で近似されるバッテリのモデルは、システム同定装置１の同定対象となる対象システム２と対応づけられうる。バッテリのＯＣＶとヒステリシス電圧との和は、対象システム２の非線形要素２４に対応づけられうる。バッテリのモデルのうち、非線形要素２４に対応づけられる要素は、バッテリ非線形要素ともいう。バッテリの過電圧は、対象システム２の線形要素２２に対応づけられうる。バッテリのモデルのうち、線形要素２２に対応づけられる要素は、バッテリ線形要素ともいう。式（１９）〜（２１）に含まれるＲ０、Ｒｄ及びτｄは、線形要素２２に含まれる伝達関数をパラメトライズするパラメータ（θ）に含まれうる。
【００６２】
バッテリ非線形要素は、以下の式（２２）及び式（２３）のように表されうる。
【数１８】
2018190288.tif 000020
【００６３】
式（２２）によれば、バッテリ非線形要素は、α＞０においてｄ＝１の入力（ｕ（ｔ））に対して斉次なシステムであるといえる。本実施形態に係るシステム同定装置１は、例えば図４のフローチャートに示される処理によって、端子電圧のデータからバッテリ非線形要素の応答を消去しうる。バッテリ非線形要素の応答が消去されることによって、システム同定装置１は、バッテリ線形要素である過電圧の応答を用いて、Ｇη（ｐ）を同定しうる。
【００６４】
以下、バッテリのモデルのパラメータを同定するシミュレーションの例が説明される。
【００６５】
システム同定装置１は、生成部１４によってｕ１（ｔ）とｕ２（ｔ）とを生成し、バッテリのモデルに入力する。ｕ１（ｔ）は、図８（Ａ）のグラフに例示される信号であるものとする。ｕ２（ｔ）は、図８（Ａ）に例示されるｕ１（ｔ）に対して式（６）の関係が成り立つ信号であるものとする。
【００６６】
システム同定装置１は、取得部１６によってバッテリのモデルから出力を取得する。バッテリのモデルの出力は、図８（Ｂ）に示される過電圧（η）と、図８（Ｃ）に示されるＯＣＶと、ヒステリシス電圧とを含む。システム同定装置１は、実際には、出力から直接過電圧、ＯＣＶ、及びヒステリシス電圧を分けて取得できない。図８（Ｂ）及び図８（Ｃ）に示される過電圧及びＯＣＶの波形は、シミュレーションの結果として取得されうる。
【００６７】
図８（Ｂ）に示されるように、ｕ１（ｔ）に対応する出力であるｙ１（ｔ）に含まれる過電圧（η）の波形（実線）と、ｕ２（ｔ）に対応する出力であるｙ２（ｔ／αｄ）に含まれる過電圧（η）の波形（破線）とは、異なっている。一方で、図８（Ｃ）に示されるように、ｙ１（ｔ）に含まれるＯＣＶの波形と、ｙ２（ｔ／αｄ）に含まれるＯＣＶの波形とは一致している。図示されていないが、ｙ１（ｔ）に含まれるヒステリシス電圧の波形と、ｙ２（ｔ／αｄ）に含まれるヒステリシス電圧の波形とは、ＯＣＶの波形と同様に、一致している。つまり、非線形要素２４に係る出力は、ｙ１とｙ２とで一致する。
【００６８】
システム同定装置１は、図８（Ｂ）及び図８（Ｃ）に例示された波形を含む出力に基づいて、差分信号を算出しうる。非線形要素２４に係る出力がｙ１とｙ２とで一致することから、差分信号は、非線形要素２４に係る出力を含まない。つまり、本実施形態に係るシステム同定装置１及びシステム同定方法によれば、非線形要素２４に係る出力が消去されうることがシミュレーションによって確認された。
【００６９】
非線形要素２４に係る出力が消去された後、線形要素２２に係る出力に基づいて、過電圧に係るパラメータが同定されうる。シミュレーションでは、例えば、式（１９）に含まれるＲ０の推定値の平均がＲ０の真値と一致することが確かめられた。
【００７０】
本実施形態に係るシステム同定装置１は、バッテリのパラメータ推定装置として、電気自動車又はハイブリッド電気自動車等の車両に用いられてよい。
【００７１】
本開示に係る実施形態について諸図面及び実施例に基づき説明してきたが、当業者であれば本開示に基づき種々の変形又は修正を行うことが容易であることに注意されたい。従って、これらの変形又は修正は本開示の範囲に含まれることに留意されたい。例えば、各構成部、各ステップなどに含まれる機能などは論理的に矛盾しないように再配置可能であり、複数の構成部およびステップなどを１つに組み合わせたり、或いは分割したりすることが可能である。
【符号の説明】
【００７２】
１  システム同定装置
１０  制御部
１２  記憶部
１４  生成部
１６  取得部
２  対象システム
２２  線形要素
２４  非線形要素
【図１】
2018190288.tif 000021
【図２】
2018190288.tif 000022
【図３】
2018190288.tif 000023
【図４】
2018190288.tif 000024
【図５】
2018190288.tif 000025
【図６】
2018190288.tif 000026
【図７】
2018190288.tif 000027
【図８】
2018190288.tif 000028