(57)【要約】【課題】蓄電池の有効的な診断を行う。【解決手段】本発明の一態様としての蓄電池診断装置は、変換計算部と、残存価値計算部とを備える。前記変換計算部は、蓄電池の複数の状態量を、前記蓄電池の用途に応じた判別空間変換関数により、１次元以上の次元数の判別空間上の座標に写像する。前記残存価値計算部は、前記判別空間を蓄電池の使用可能領域と使用不可領域とに区分する電池性能判別関数の出力値を前記座標の値に基づき求め、前記電池性能判別関数の出力値に基づき、前記蓄電池の残存価値を計算する。【選択図】図１
2014059270.tif 000002
【特許請求の範囲】
【請求項１】
蓄電池の複数の状態量を、前記蓄電池の用途に応じた判別空間変換関数により、１次元以上の次元数の判別空間上の座標に写像する変換計算部と、
前記判別空間を蓄電池の使用可能領域と使用不可領域とに区分する電池性能判別関数の出力値を前記座標の値に基づき求め、前記電池性能判別関数の出力値に基づき、前記蓄電池の残存価値を計算する残存価値計算部と
を備えた蓄電池診断装置。
【請求項２】
複合残存価値計算部をさらに備え、
前記変換計算部は、前記蓄電池の複数の用途のそれぞれに応じた判別空間関数により、前記複数の状態量を、前記用途毎の判別空間上の座標に写像し、
前記残存価値計算部は、前記用途毎の判別空間を使用可能領域と使用不可領域とに区分する電池性能判別関数の出力値を前記座標の値に基づき計算し、前記電池性能判別関数の出力値に基づき、前記用途毎の前記蓄電池の残存価値を計算し、
前記複合残存価値計算部は、前記用途毎の前記残存価値を前記用途ごとの重み係数で重み付け総和することにより、前記蓄電池の複合残存価値を計算する
請求項1に記載の蓄電池診断装置。
【請求項３】
前記蓄電池を測定して前記蓄電池の複数の状態量を取得する電池状態測定部
をさらに備えた請求項1または2に記載の蓄電池診断装置。
【請求項４】
前記蓄電池の用途を1つ以上入力する蓄電池用途入力部
をさらに備えた請求項1ないし3のいずれか一項に記載の蓄電池診断装置。
【請求項５】
第1測定時点および前記第1測定時点より後の第2測定時点のそれぞれにおける蓄電池の複数の状態量を、前記蓄電池の用途に応じた判別空間変換関数により、１次元以上の次元数の判別空間上の第1座標および第2座標に写像する第1計算部と、
前記判別空間を蓄電池の使用可能領域と使用不可領域とに区分する電池性能判別関数について、前記第1座標および前記第2座標が前記使用可能領域に属する場合に、前記第1座標から前記電池性能判別関数上の最近傍点を求め、前記第1座標から前記第2座標へのベクトルについて、前記第1座標から前記最近傍点への方向の前記ベクトルの成分を計算する第2計算部と、
前記第1測定時点および前記第2測定時点間の時間長に対する前記成分の比率に基づいて、前記蓄電池の劣化進行を診断する劣化判定部と、
を備えた蓄電池診断装置。
【請求項６】
前記劣化判定部は、前記比率が閾値を超えた場合に、警告信号を出力する
請求項5に記載の蓄電池診断装置。
【請求項７】
前記劣化判定部は、前記第2座標と前記電池性能判別関数との距離が閾値以下になったら警告信号を出力する
請求項5または6に記載の蓄電池診断装置。
【請求項８】
第1測定時点および前記第1測定時点より後の第2測定時点のそれぞれにおける蓄電池の複数の状態量を、前記蓄電池の用途に応じた判別空間変換関数により、１次元以上の次元数の判別空間上の第1座標および第2座標にそれぞれ写像し、前記第1座標から前記第2座標へのベクトルを前記第2座標および第1座標の測定時点間の時間長で除算して第1ベクトルを取得する第1ベクトル取得部と、
前記蓄電池とは別の少なくとも1つの蓄電池について、複数の測定時点毎に当該蓄電池の複数の状態量を含むデータを格納する蓄電池状態量記憶部と、
前記格納部内の各前記データを読み出し、各前記データを前記判別空間変換関数により前記判別空間上の座標へ写像し、前記座標の中から前記第2座標との距離に応じて少なくとも1つの座標を第3座標として選択し、前記座標のうち前記第3座標の次の測定時点に対応する座標である第4座標に対する前記第3座標からのベクトルを前記第4座標および第3座標の測定時点間の時間長で除算して第2ベクトルを取得する第2ベクトル取得部と、
前記第1ベクトルと前記第2ベクトルを合成して予測ベクトルを得る予測ベクトル取得部と、
前記第2測定時点より後の第3測定時点と前記第2測定時点間の時間長を前記予測ベクトルに乗算し、乗算結果を前記第2座標に加算することにより、前記判別空間において前記第3測定時点の前記蓄電池の状態量の写像先を予測し、前記判別空間を蓄電池の使用可能領域と使用不可領域とに区分する判別関数と前記写像先とに基づいて、前記蓄電池の将来使用可否を診断する将来使用可否決定部と
を備えた蓄電池診断装置。
【請求項９】
前記将来使用可否決定部は、前記写像先が前記使用不可領域に属するときは、警告信号を出力する
ことを特徴とする請求項8に記載の蓄電池診断装置。
【請求項１０】
前記将来使用可否決定部は、前記写像先が前記使用可能領域に属し、前記写像先と前記判別関数との距離が閾値以下のときは、警告信号を出力する
ことを特徴とする請求項8または9に記載の蓄電池診断装置。
【請求項１１】
蓄電池の複数の状態量を、前記蓄電池の用途に応じた判別空間変換関数により、１次元以上の次元数の判別空間上の座標に写像する変換計算ステップと、
前記判別空間を蓄電池の使用可能領域と使用不可領域とに区分する電池性能判別関数の出力値を前記座標の値に基づき求め、前記電池性能判別関数の出力値に基づき、前記蓄電池の残存価値を計算する残存価値計算ステップと
を備えた蓄電池診断方法。
【請求項１２】
第1測定時点および前記第1測定時点より後の第2測定時点のそれぞれにおける蓄電池の複数の状態量を、前記蓄電池の用途に応じた判別空間変換関数により、１次元以上の次元数の判別空間上の第1座標および第2座標に写像する第1計算ステップと、
前記判別空間を蓄電池の使用可能領域と使用不可領域とに区分する電池性能判別関数について、前記第1座標および前記第2座標が前記使用可能領域に属する場合に、前記第1座標から前記電池性能判別関数上の最近傍点を求め、前記第1座標から前記第2座標へのベクトルについて、前記第1座標から前記最近傍点への方向の前記ベクトルの成分を計算する第2計算ステップと、
前記第1測定時点および前記第2測定時点間の時間長に対する前記成分の比率に基づいて、前記蓄電池の劣化進行を診断する劣化判定ステップと、
を備えた蓄電池診断方法。
【請求項１３】
第1測定時点および前記第1測定時点より後の第2測定時点のそれぞれにおける蓄電池の複数の状態量を、前記蓄電池の用途に応じた判別空間変換関数により、１次元以上の次元数の判別空間上の第1座標および第2座標にそれぞれ写像し、前記第1座標から前記第2座標へのベクトルを前記第2座標および第1座標の測定時点間の時間長で除算して第1ベクトルを取得する第1ベクトル取得ステップと、
前記蓄電池とは別の少なくとも1つの蓄電池について、複数の測定時点毎に当該蓄電池の複数の状態量を含むデータを格納する蓄電池状態量記憶部から前記データを読み出すステップと、
前記格納部内の各前記データを読み出し、各前記データを前記判別空間変換関数により前記判別空間上の座標へ写像し、前記座標の中から前記第2座標との距離に応じて少なくとも1つの座標を第3座標として選択し、前記座標のうち前記第3座標の次の測定時点に対応する座標である第4座標に対する前記第3座標からのベクトルを前記第4座標および第3座標の測定時点間の時間長で除算して第2ベクトルを取得する第2ベクトル取得ステップと、
前記第1ベクトルと前記第2ベクトルを合成して予測ベクトルを得る予測ベクトル取得ステップと、
前記第2測定時点より後の第3測定時点と前記第2測定時点間の時間長を前記予測ベクトルに乗算し、乗算結果を前記第2座標に加算することにより、前記判別空間において前記第3測定時点の前記蓄電池の状態量の写像先を予測し、前記判別空間を蓄電池の使用可能領域と使用不可領域とに区分する判別関数と前記写像先とに基づいて、前記蓄電池の将来使用可否を診断する将来使用可否決定ステップと
を備えた蓄電池診断方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施形態は、蓄電池診断装置およびその方法に関する。
【背景技術】
【０００２】
蓄電池はPHEV（Plug-in Hybrid Electric Vehicle）の動力源として、また再生可能エネルギーの出力安定化装置として需要が高まっており、高信頼化・長寿命化を実現するための研究開発が進んでいる。一方で、蓄電池性能向上の結果、蓄電池を組み込んだ製品側がの寿命を終えても、蓄電池が十分な性能保持している可能性が懸念されている。まだ十分な性能を持つ蓄電池を廃棄することは、環境負荷や経済性の観点から避けるべきである。解決策としては蓄電池の二次利用の仕組みを技術・制度の両面から進めることが検討されている。
【０００３】
蓄電池の二次利用に必要な技術の一つは、対象となった蓄電池がどれだけの経済的価値(残存価値)を持つかを計算することである。
【０００４】
従来、主に再利用(二次利用)の検査のプロセスに関する技術が知られている。ユーザが店舗等に持ち込んだ蓄電池を測定し、測定値に基づいて継続使用可能か否かを判定する。判定までの処理は外観検査や内部抵抗の測定等である。この技術では、継続使用可否は判定できるものの、蓄電池の使用者や回収業者は、現在の対象となる蓄電池の残存価値については知ることができない。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特許登録第4033130号
【発明の概要】
【発明が解決しようとする課題】
【０００６】
本発明の一側面は、上記問題点を解決するためになされたものであり、蓄電池を有効的に診断することが可能な蓄電池診断装置およびその方法を提供する。
【課題を解決するための手段】
【０００７】
本発明の一態様としての蓄電池診断装置は、変換計算部と、残存価値計算部とを備える。
【０００８】
前記変換計算部は、蓄電池の複数の状態量を、前記蓄電池の用途に応じた判別空間変換関数により、１次元以上の次元数の判別空間上の座標に写像する。
【０００９】
前記残存価値計算部は、前記判別空間を蓄電池の使用可能領域と使用不可領域とに区分する電池性能判別関数の出力値を前記座標の値に基づき求め、前記電池性能判別関数の出力値に基づき、前記蓄電池の残存価値を計算する。
【図面の簡単な説明】
【００１０】
【図１】本発明の実施形態に係る蓄電池診断装置のブロック図である。
【図２】図1の装置の動作フローチャートである。
【図３】蓄電池状態量履歴データを例示する図である。
【図４】電池状態量空間から電池性能判別空間への写像例を示す図である。
【図５】電池性能判別関数の例を示す図である。
【図６】電池変化量データ計算部の動作説明図である。
【図７】状態量予測値計算部の動作説明図である。
【図８】複合残存価値の計算例を示す図である。
【図９】判別空間変換関数および電池性能判別関数の学習例を示す図である。
【図１０】学習用蓄電池状態量および学習用蓄電池使用可／不可符号を例示する図である。
【発明を実施するための形態】
【００１１】
以下、図面を参照しながら、本発明の実施形態を説明する。
【００１２】
図1に本発明の実施形態に係る蓄電池診断装置のブロック図を示す。図2に、図1の装置による蓄電池の診断実行時の動作例のフローチャートを示す。本明細書中のステップ番号STはフローチャート中の番号に対応する。
【００１３】
（ST1）
ユーザが蓄電池を電池状態測定部1にセットする。ここでいう蓄電池は測定の構成単位を示す表現であり、セル(単一の電池)やモジュール(複数のセルを組み合わせて所定の電圧・電流を得られるようにしたもの)などの概念を含むものである。電池状態測定部1は、充放電器や交流測定器、X線装置など蓄電池の状態を計測できる装置であれば、特定の形態に限定されない。
【００１４】
（ST2）
次に対象蓄電池が再利用される可能性のある用途を蓄電池用途入力部3から入力する。再利用される可能性のある用途は、たとえば、「EV」「家庭用定置」「電動スクータ」等である。
【００１５】
この用途は、対象蓄電池のユーザが転売を望むものを選択すればよい。ユーザがEVとしての二次利用のみを望むのであれば、「EV」を入力すればよい。蓄電池用途入力部3からの入力値の範囲は、事前に設計者が設定し、ユーザに値を選択させるものとする。
【００１６】
蓄電池用途入力部3へ入力された用途情報は、ステップST5以降の処理で使用される。
【００１７】
（ST4）
電池ID入力部11から、対象蓄電池を特定するのに必要な情報を入力する。具体的には蓄電池側に具備されたIC等に記録された製造番号等を電子的に読み取ることが想定される。別の例として、同内容の情報をユーザがキーボードやバーコードリーダを用いて入力することも可能である。
【００１８】
（ST3）
蓄電池状態測定部1は、ステップST1でセットされた蓄電池のd次元状態量を計算する。
【００１９】
蓄電池状態測定部1の出力は、d次元の実数値の組である。これは、一種の多変量データである。1個の蓄電池の1回の測定からは、d個の実数値が得られる。
【００２０】
測定値として、充電時のデータの分析により得られる満充電までの所要時間や、内部抵抗値、電圧立ち上がりカーブの曲率の総和などの数値が想定される。また測定値は、X線等を用いた蓄電池内部の活物質状態の推定値であってもよい。
【００２１】
（ST13）
蓄電池状態量履歴データ記憶部27は、蓄電池状態計算部1が計算したd次元状態量を、電池ID入力部11から入力されたID情報と対応づけて内部に格納する。また、測定回数・測定日時（測定時間情報）も同時に記録する。測定時間情報は、診断履歴の変化傾向を計算するために必要な情報となる。図3に蓄電池状態量履歴データ記憶部27に格納されるデータ例を示す。
【００２２】
ID情報の数学的形式は限定されないが、特定蓄電池の一意性を確認できるものであることが必要である。この情報によって過去の測定データとの連続性を確保することが可能になる。
【００２３】
以下のST5、ST6、ST8、ST18、ST19、ST7、ST9の各ステップは、ステップST2で入力された用途数分、実行される。またステップST5から分岐される3つのパスは、並列実行しても逐次実行しても差し支えない。
【００２４】
（ST5）
電池性能判別空間変換計算部2は、電池性能判別空間変換関数記憶部14から当該用途に対応した判別空間変換関数を読み出し、当該判別空間変換関数により、電池状態測定部1で測定された電池状態量を、電池性能判別空間上の座標に写像する。電池性能判別空間変換関数記憶部14は、用途ごとに判別空間変換関数を記憶している。
【００２５】
電池状態量空間から電池性能判別空間への写像を模式的に図4に示す。電池状態量空間は、電池状態量の各状態量を軸とした空間である。電池性能判別空間は、蓄電池の使用可能/不可能に関する特徴を最もよく表した部分空間である。電池性能判別空間の次元数および電池性能判別空間の各軸を表す状態量（変数）は、使用する判別空間変換関数に応じて異なる。図中の○や×は、判別空間変換関数および後述する電池性能判別関数の学習の際に用いたデータを表している。○は当該用途で使用可能と判断されている蓄電池のデータ、×は使用不可と判断されている蓄電池のデータに対応する。
【００２６】
判別空間変換関数の数学的形式は限定されないが、多変量解析手法における線形判別分析の線形変換関数や、カーネル法と呼ばれる新しい多変量解析手法で用いられる非線形変換関数が想定される。
【００２７】
線型変換関数の場合、変換後の空間の次元が2次元とすると、線型変換関数はd×2の行列Aで与えることができる。したがって、電池状態量のベクトルデータをxとおくと、電池性能判別空間上の当該電池の2次元ベクトル値yは、y=Axの線型変換で得ることができる。行列Aを定める方法については後述する。
【００２８】
（ST6）
電池変化量データ計算部（第1計算部、第2計算部）9は、電池性能判別関数記憶部17から当該用途に対応する電池性能判別関数を読み出す。電池性能判別関数記憶部17は、用途ごとに電池性能判別関数を記憶している。
【００２９】
電池性能判別関数は、電池の使用可/不可を判別する関数である。電池性能判別関数を電池性能判別空間に描いた例を図5に示す。電池性能判別空間上では、電池性能判別関数は、使用可/不可を分割する超平面または曲面を表す。
【００３０】
電池変化量データ計算部9は、電池性能判別空間上の状態量(以降、変換後電池データと称する)を表す座標と、電池性能判別関数記憶部17から読み出した電池性能判別関数との距離を計算する。当該電池性能判別空間上の状態量は、今回測定された電池状態量を上記判別空間変換関数により変換して取得してもよいし、電池性能判別空間変換計算部2により計算された値をそのまま用いてもよい。
【００３１】
電池性能判別関数をf、変換後電池データをyとすると、関数値f(y)の値の絶対値を距離として採用してもよい。あるいは、fが線型平面であれば、点と平面の距離を計算する幾何学公式を用いてもよい。またfがカーネルサポートベクターマシンとよばれる非線型関数であれば、fを高次元の線形空間の関数とみなして、平面距離計算の公式を適用して距離計算することもできる。
【００３２】
サポートベクターマシンの詳細は渡辺他「学習システムの理論と実現」森北出版2005 ISBN:4-627-82941-8等に詳しい。
【００３３】
サポートベクターマシンの関数fは、以下の形式で与えられる。
【数１】
2014059270.tif 000003
yは変換後電池データであり、zは係数となる実数値である。Cは切片である。
【００３４】
またki(y)は、非線型関数であり、以下の数式で表される。
【数２】
2014059270.tif 000004
αはiごとに定まるベクトル定数であり、σは全てのki(y)で共通の定数である。
【００３５】
次に、電池変化量データ計算部9は、対象電池のIDデータを元に、蓄電池状態量履歴データ記憶部27から過去（たとえば前回）に測定した状態量データと、今回測定した状態量データを読み出し、当該用途に応じた電池性能判別空間変換関数を用いて、それぞれ変換後電池データに変換する。そして、これらの変換後電池データに基づき、過去（前回）の状態量を表す座標から電池性能判別関数面方向への単位時間あたりの移動量（移動成分）を計算する。
【００３６】
図6を用いて、移動成分の計算例を示す。図6における（2）の時点を今回の測定時点とし、（1）の前回の測定時点からの電池性能判別関数方向の移動成分の計算例を説明する。
【００３７】
前回測定時点（1）の点（座標）をA、今回測定時点（2）の点をBと表現する。TA、TBはそれぞれ、前回および今回の電池状態量の測定時点を表す。なお、fBは、点Bから電池性能判別関数までの距離を表す。
【００３８】
また点Aの電池性能判別関数上の再近傍点をCとする。
【００３９】
点Cは電池性能判別関数の関数形によっては解析的に定めることが可能であるが、解析的手段が不可能な場合は電池性能判別関数上の各点との距離を評価して探索的に定めることで対応する。
【００４０】
最初に点Aから点Cに向かうベクトルACを求める。
【００４１】
これは、AC=C-Aで計算できる。
【００４２】
次に点Aから点Bに向かうベクトルABを求める。
【００４３】
これは、AB=B-Aで計算できる。
【００４４】
さらにベクトルABとACの間の余弦cosθを計算する
cosθは以下の式で計算できる。
【数３】
2014059270.tif 000005
【００４５】
電池性能判別関数方向の移動成分MはベクトルABのAC上への射影の長さを測定期間で除したものとして計算される。測定期間はTB-TAで計算される。
【００４６】
これは単位時間当たりの使用不可領域への接近度を示すものであり、電池が急速に劣化することを検出する指標といえる。
【数４】
2014059270.tif 000006
【００４７】
（ST8）
急速劣化判定部10は、電池変化量データ計算部9で計算された電池性能判別関数面との距離が閾値以下になったか、または移動成分が閾値を超えたかを判断する。当該距離が閾値以下になった、または移動成分が閾値を超えた場合、あるいはこれらの両方が成立した場合、当該用途での使用可能限界が近づいたと判断する。このとき、急速劣化判定部10は、電池状態表示部4に警告信号を送信する。
【００４８】
図6の例では、距離fBが閾値以下になったとき、または移動成分Mが閾値を超えたとき、あるいはこれらの両方が成立した場合、当該用途での使用可能限界が近づいたと判断する。 
（ST18）
状態量予測値計算部（第1ベクトル取得部、第2ベクトル取得部、予測ベクトル取得部）25は、蓄電池状態量履歴データ記憶部27から対象蓄電池の状態量データの履歴を読み出す。また、電池性能判別空間変換関数記憶部14から用途に応じた判別空間変換関数を読み出す。履歴における対象蓄電池の今回と前回の状態を、当該判別空間変換関数によりそれぞれ変換後電池データに変換する。電池性能判別空間において、前回の変換後電池データを表す座標から、今回の変換後蓄電池データを表す座標への、変化方向および単位時間あたりの変化量を表すベクトル（第1ベクトルと呼ぶ）を計算する。
【００４９】
また、状態量予測値計算部25は、学習用蓄電池状態データ記憶部15の各データを該当用途に応じた判別空間変換関数により変換し、変換後蓄電池データの中から、今回測定時の対象蓄電池の変換後蓄電池データと類似したデータ（たとえば近傍に位置するデータ）を1個または複数個選択する。たとえば、各変換後蓄電池データを表す座標の中から、対象蓄電池の変換後蓄電池データを表す座標に対する距離に基づき、1個または複数個の座標を選択する。
【００５０】
ここで、学習用蓄電池状態データ記憶部15と学習用蓄電池使用可／不可符号記憶部16について説明する。学習用蓄電池状態データ記憶部15は、蓄電池状態量履歴データ記憶部27と同様の形式で、蓄電池の状態量、測定日時、測定回数、IDを記憶している。蓄電池状態量履歴データ記憶部27はユーザの蓄電池を測定したデータであったが、学習用蓄電池状態データ記憶部15のデータは、管理者が判別空間変換関数および電池性能判別関数の学習用データとして、ユーザの蓄電池とは別個の1つ以上の蓄電池から、複数の時点で測定して取得したものである。また、学習用蓄電池使用可／不可符号記憶部16は、学習用蓄電池状態データ記憶部15のデータと一対一に対応して、当該蓄電池が用途毎に使用可能な否かを示す真理値データを記憶する。学習用蓄電池状態データ記憶部15と学習用蓄電池使用可／不可符号記憶部16は1つの記憶部としてまとめても良い。図10に、学習用蓄電池状態データ記憶部15と学習用蓄電池使用可／不可符号記憶部16の1つの記憶部にまとめた場合のデータ例を示す。別々の記憶部とする場合は、少なくともIDと測定日時とを共通の項目として持てばよい。
【００５１】
状態量予測値計算部25は、選択した変換後蓄電池データと、当該選択した変換後蓄電池データの蓄電池の次の測定時点に対応する変換後蓄電池データ間で、電池性能判別空間上の変化方向と単位時間あたりの変化量を表すベクトル（第2ベクトルと呼ぶ）を計算する。つまり、第2ベクトルは、当該選択した変換後蓄電池データを表す座標から、当該次の測定時点に対応する変換後蓄電池データを表す座標への、変化方向と単位時間あたりの変化量を表す。
【００５２】
第1ベクトルと第2ベクトルを合成することで予測ベクトルを生成する。予測ベクトルは、対象蓄電池について、次回測定時に対する変化方向および単位時間あたりの変化量を表すベクトルである。
【００５３】
図7を用いて、具体例を説明する。
【００５４】
以下では（3）時点を今回の測定時点として、次回測定に対する予測ベクトルを得る例を示す。
【００５５】
最初に（2）の前回の測定時点から、(3)の今回の測定時点へのベクトルV2を計算する。
【００５６】
ベクトルV2は、（3）時点で測定した変換後電池データの座標から、前回測定した変換後電池データの座標の差分を取ることで得られる。すなわちV2=C-Bである。
【００５７】
そして、ベクトルV2を、TCおよびTB間の測定間隔時間で割って正規化することで、ベクトルV2’（第1ベクトル）を得ることができる。
【数５】
2014059270.tif 000007
【００５８】
次に状態量予測値計算部25は、学習用蓄電池状態量記憶部15から過去に測定された蓄電池の全状態量を読み出し、さらに電池性能判別空間変換関数記憶部14から当該用途に応じた判別空間変換関数を読み出す。各状態量を、判別空間変換関数により、それぞれ変換後電池データに変換する。
【００５９】
そして変換後電池データの集合から、（3）時点での変換後電池データC（今回測定した蓄電池の状態量の変換後電池データ）と最も類似するデータを選択する。類似の基準は電池性能判別空間上のユークリッド距離を想定するのが自然であるが、その他の距離関数を用いてもよい。
【００６０】
図7の例では、データXが、（3）時点のデータCの最近傍点に対応するため選択される。そして、選択したデータXと、当該データの次の測定時点のデータYとの間のベクトルV3を計算する。図7において、TxはデータXの測定時点、TyはデータYの測定時点を表す。
【００６１】
ベクトルV3は、選択したデータXと、次回測定時のデータYとの差分で計算される。すなわちV3=Y-Xである。
【００６２】
そして、ベクトルV3’（第2ベクトル）は、ベクトルV3を、TyおよびTx間の測定間隔の時間で割って正規化することで得られる。
【数６】
2014059270.tif 000008
【００６３】
V2’とV3’を合成することで、予測ベクトルV4を得ることができる。
【数７】
2014059270.tif 000009
【００６４】
（ST19）
将来使用不可決定部26は、状態量予測値計算部25で計算された予測ベクトルに、次回測定までの時間ΔTを乗じ、この結果を、今回測定された電池の変換後電池データに加算することで、次回測定時の変換後蓄電池データの予測値D（図7参照）を得る。次回の測定時点がTDと表すと、ΔT＝TD-TCである。この場合、次回測定時の予測値Dは以下のように計算される。
【数８】
2014059270.tif 000010
【００６５】
将来使用不可決定部26は、予測値Dが使用不可領域に入れば、警告信号を生成して、電池状態表示部4に出力する。あるいは、予測値Dが使用可能領域に属するが、前記予測値Dと電池性能判別関数との距離が閾値以下のときは、警告信号を生成して出力してもよい。
【００６６】
ステップS18で、類似する変換後蓄電池データを複数個読み出した場合は、複数の予測値が得られる。この場合、これらの予測値の平均、最悪値、または最良値が、使用不可領域に入れば、警告信号を生成する。
【００６７】
（ST7）
電池使用可能判別量計算部（残存価値計算部）7は、電池性能判別関数記憶部17から当該用途に対応する電池性能判別関数を読み出し、電池性能判別空間変換計算部2で計算された変換後蓄電池データに基づき、対象蓄電池について、当該用途の単純残存価値を計算する。図8の左側に、中古EV、定置蓄電池、スクータの3つの用途に応じてそれぞれ単純残存価値を計算する例を示す。この３用途は説明要の一例であり，本発明はこれら3用途に限定されるものではない．
電池性能判別関数の数学的な定義は限定されないが、電池性能判別空間上の状態量（変換後蓄電池データ）を入力とし、戻り値は実数値とする。戻り値の符号が正なら対象電池は使用可能、負なら使用不可能に対応する（図4参照）。電池性能判別関数fの出力を単純残存価値と定義してもよいが、0以下の負値は全て0とみなす等のフィルタ処理を加えてもよい。また電池性能判別関数fの出力が負（使用不可能）であれば、警告信号を生成して、電池状態表示部4に通知する。
【００６８】
（ST9）
電池状態表示部4が、電池使用可能判別量計算部7により計算された単純残存価値を表示する。
【００６９】
（ST11）
複合残存価値計算部8は、電池使用可能判別量計算部7で計算された各々の用途の各単純残存価値を、調整係数計算部6から入力された各調整係数を乗算して総和することで、複合残存価値を得る。図8の右側に、各単純残存価値を調整係数により重み付け総和して複合残存価値を得る例を示す。調整係数計算部6は、用途毎の調整係数の値を記憶している。調整係数計算部6は、複合残存価値計算部8からの要求に応じて用途毎の調整係数を出力する。
【００７０】
複合残存価値の計算式の一例を以下に示す。
【００７１】
最初に諸記号を定義する
【数９】
2014059270.tif 000011
【００７２】
各αは調整係数計算部6から入力される重み係数である。各αは、単純残存価値→金額への変換、または当該用途の相対的市場規模を反映するパラメタである。
【００７３】
（ST12）
電池状態表示部4は、複合残存価値計算部8で計算された複合残存価値を表示する。
【００７４】
以上の処理により、蓄電池のユーザは自分の蓄電池がどれだけの経済上の価値を持つのかを定量的に知ることができ、合わせて近い将来までの使用可能性を予測の形で知ることができるようになる。
【００７５】
次に、図9のフローチャートを用いて、判別空間変換関数及び電池性能判別関数の学習プロセスを示す。以下のST14、ST15、ST16のステップは設計者が設定した二次利用用途の数だけ実行される。
【００７６】
（ST13）
学習用データを学習用蓄電池状態データ記憶部15に格納し、学習用データに一対一に対応して使用可/不可を示す符号を学習用蓄電池使用可/不可符号記憶部16に格納し、学習を開始する。
【００７７】
学習用蓄電池状態データ記憶部15に記憶されるデータは蓄電池状態量履歴データ記憶部27と同形式であり、電池状態量・測定日時・測定回数・IDの4項目を含む(図10参照)。
【００７８】
また学習用蓄電池使用可/不可符号記憶部16は使用可/不可の符号を格納する(図10参照)。
【００７９】
使用可/不可符号は各二次利用用途における使用可/不可を表す真理値であり、前設定用途数ごとに定義される。したがって、一個の状態量の組に対して用途の数だけ真理値が対応付けされる。
【００８０】
（ST14）
電池性能判別空間変換関数計算部13は、学習用蓄電池状態データ記憶部15から学習用データ（蓄電池状態量）を学習に必要な所定の個数読み出し、さらに読み出した蓄電池状態量に対応する使用可/不可符号を学習用蓄電池使用可/不可符号記憶部16から読み出す。これらの読み出したデータを用いて、電池性能判別空間への変換関数を学習し、電池性能判別空間変換関数記憶部14に記憶する。
【００８１】
変換関数の例としては、以下に示すパターン認識の分野等次元数削減のために用いられる線型判別等が想定される。
【００８２】
同方法で用いられる変換関数の学習方法は、以下の文献に詳しい。  
石井他「わかりやすいパターン認識」1998オーム社 ISBN4-274-13149-1
【００８３】
ここでは上記文献の方法をもとに変換関数の学習方法の実施例の一つを説明する。
【００８４】
最初に各クラス内変動行列Sを以下で定義する。ここでのクラスは使用可/不可の各電池状態量の集合を指す。
【００８５】
Xは各クラス内の各電池状態量xの集合であり、d次元とする。mはクラス内の平均ベクトルである。
【数１０】
2014059270.tif 000012
【００８６】
ここから、各クラス内変動行列の総和として、クラス内変動行列SWを定めることができる。
【数１１】
2014059270.tif 000013
【００８７】
次にクラス間変動行列SBを、以下のように定める。niは各クラス内のデータ数、miは各クラスの平均ベクトル、mは全てのクラスのデータの平均である。
【００８８】
SBが大きければ大きいほど、2つの電池のクラスは空間上で離れていると考えられる。
【数１２】
2014059270.tif 000014
【００８９】
次に変換行列Aを導入する。ここでは(d、1)行列と仮定する。変換後電池データyは、以下で計算される。
【数１３】
2014059270.tif 000015
【００９０】
ここで全データを変換後電池データに変換した場合、変換後の空間で電池データの使用可/不可がもっとも明確に分離するように変換行列Aを定めればよい。
【００９１】
ここで変換後のクラス内変動行列
【数１４】
2014059270.tif 000016
を以下で定める。
【数１５】
2014059270.tif 000017
【００９２】
また変換後のクラス間変動行列
【数１６】
2014059270.tif 000018
を以下で定める。
【数１７】
2014059270.tif 000019
【００９３】
このとき、各クラスの分離度を示す指標として、次の関数
【数１８】
2014059270.tif 000020
を導入する。
【数１９】
2014059270.tif 000021
【００９４】
この関数値が大きければ、変換後電池データは使用可/不可ごとに小さくまとまり、かつ互いの集団間距離は乖離していると考えることができる。
【数２０】
2014059270.tif 000022
を最大にするAが、求める変換関数（変換行列）である。
【００９５】
Aを求める手順としては、最初に次の制約を導入する。
【数２１】
2014059270.tif 000023
【００９６】
そして上記制約のもとで
【数２２】
2014059270.tif 000024
を最大にするAを求める。
【００９７】
ここでラグランジュ乗数を導入すると、次のラグランジュ関数
【数２３】
2014059270.tif 000025
が得られる。
【数２４】
2014059270.tif 000026
【００９８】
これをAで微分して0とおくと以下が得られる。
【数２５】
2014059270.tif 000027
【００９９】
このとき
【数２６】
2014059270.tif 000028
の最大固有値に相当するλが、JS(A)を最大にするλである。
【数２７】
2014059270.tif 000029
の固有値問題を解き、対応する固有ベクトルを求めることでAが定まる。
【０１００】
さらに第二・第三位の固有値に対応する固有ベクトルにより、電池性能判別空間の複数の次元を定めることができる。
【０１０１】
これらの固有ベクトルを結合し、最終的な電池性能判別空間変換関数を得る。
【０１０２】
（ST15）
電池性能判別関数計算部8は、学習用蓄電池状態データ記憶部15から学習用データ（蓄電池状態量のデータ）を学習に必要な所定の個数読み出し、ステップS14で学習された判別空間変換関数を用いて、電池性能判別空間へ当該蓄電池状態量を変換する。
【０１０３】
（ST16）
電池性能判別関数計算部8は、ステップS15で読み出した蓄電池状態量に対応する符号を学習用蓄電池使用可/不可符号記憶部16から読み出す。使用可/不可符号の情報により、電池性能判別空間上の状態量（変換後蓄電池データ）は、使用可の蓄電池と不可の蓄電池を区別することが可能となる。
【０１０４】
電池性能判別関数は、電池性能判別空間上の状態量を入力とし、この使用可/不可符号を再現するように（たとえば使用可のものは正、使用不可のものは負の値が出力されるように）、所定の学習アルゴリズムを用いて計算される。
【０１０５】
たとえば電池性能判別関数の例として、前述のサポートベクターマシンを用いた場合、関数の計算に必要な係数行列は最適化問題をとき、必要な数値を計算すること手順で求めることができる。これらの詳細も前出の渡辺他「学習システムの理論と実現」森北出版2005 ISBN:4-627-82941-8等に詳しい。
【０１０６】
電池性能判別関数計算部8は、学習により得た電池性能判別関数を、電池性能判別関数記憶部17に格納する。
【０１０７】
前述のとおり、ステップST14、ST15、ST16を、設定された二次利用全てにおいて実施する。これにより、二次利用用途ごとに、判別空間変換関数と電池性能判別関数が得られる。
【０１０８】
（その他の実施形態）
上述したように、複数の二次利用用途を設定すれば、用途毎の単純残存価値を重み付け総和することで、複合残存価値を計算できる。一方で、二次利用用途を単一かつ一次利用用途と同じものに限定して、当該用途におけるリユースの可/不可判定、劣化診断ツールとして実現することも想定される。
【０１０９】
この際、単一の用途としてEVを設定した場合は、本実施形態をEVの蓄電池リユース判定ツールとして利用できる。
【０１１０】
なお、以上に説明した本実施形態における蓄電池診断装置は、例えば、汎用のコンピュータ装置を基本ハードウェアとして用いることで実現することが可能である。すなわち、蓄電池診断装置の各処理部は、上記のコンピュータ装置に搭載されたプロセッサにプログラムを実行させることにより実現することができる。このとき、蓄電池診断装置は、上記のプログラムをコンピュータ装置にあらかじめインストールすることで実現してもよいし、CD-ROMなどの記憶媒体に記憶して、あるいはネットワークを介して上記のプログラムを配布して、このプログラムをコンピュータ装置に適宜インストールすることで実現してもよい。また、残存価値診断装置内の記憶部は、装置内もしくは外付けのメモリ装置およびハードディスク、あるいは、CD-R, CD-RW, DVD-RAM, DVD-R 等の記録媒体によって構成されてもよい。
【図１】
2014059270.tif 000030
【図２】
2014059270.tif 000031
【図３】
2014059270.tif 000032
【図４】
2014059270.tif 000033
【図５】
2014059270.tif 000034
【図６】
2014059270.tif 000035
【図７】
2014059270.tif 000036
【図８】
2014059270.tif 000037
【図９】
2014059270.tif 000038
【図１０】
2014059270.tif 000039