(57)【要約】  本発明は、電気化学デバイスから受け取った少なくとも１つの変数（Ｓ）によって決まる少なくとも１つの特性値（Ｖａｌｉ）を取得するステップ（１０３）と、取得した特性値から前記デバイスの少なくとも１つの欠陥（Ｄｉ）を決定するステップ（１０５）とを含む、電気化学デバイスの欠陥を検出する方法に関する。この方法はウェーブレット変換を含む数学的操作を含み、その操作は受け取った変数（Ｓ）から特性値（Ｖａｌｉ）を取得するように行われる。本発明は、さらに、１つのそのような方法を行うデバイス、ならびに対応するコンピュータプログラムに関する。
2012530925.tif 000002
【特許請求の範囲】
【請求項１】
電気化学デバイスの欠陥を検出する欠陥検出方法であって、前記電気化学デバイスから受け取った少なくとも１つの変数（Ｓ）から少なくとも１つの特性値（Ｖａｌｉ）を取得するステップ（１０３）と、前記取得した値から前記電気化学デバイスの少なくとも１つの欠陥（Ｄｉ）を決定するステップ（１０５）とを含み、ウェーブレット変換を含む数学的操作が、受け取った前記変数（Ｓ）から前記特性値（Ｖａｌｉ）を取得するように行われる、欠陥検出方法。
【請求項２】
前記ウェーブレット変換は、前記取得した特性値が、値のスケール変数ａおよび値の平行移動変数ｂによって決まる少なくとも１つのウェーブレット係数Ｓａ，ｂを含む離散ウェーブレット変換である請求項１に記載の欠陥検出方法。
【請求項３】
複数の特性値が、スケールレベル変数ｊがａ未満である１組のウェーブレット係数ｗｊ，ｐを複数の組のウェーブレット係数ｗｊ＋１，ｐに分解することによって取得される、請求項２に記載の欠陥検出方法。
【請求項４】
前記値のスケール変数ａがａ＝ａ０ｊを満足し、ここで、ｊはスケールレベルであり、ａ０はスケールパラメータであり、１組の特性値が、０から所与の分解レベルまでの値ｊごとに、各組のウェーブレット係数ｗｊ，ｐをａ０組のウェーブレット係数ｗｊ＋１，ｐに連続的に分解することによって取得される、請求項３に記載の欠陥検出方法。
【請求項５】
前記所与の分解レベルが最大分解レベルＭに対応する、請求項４に記載の欠陥検出方法。
【請求項６】
前記決定するステップが、前記特性値を、少なくとも１つの第１の欠陥クラス（Ｃ１）を第２の欠陥クラス（Ｃ２）から分離する少なくとも１つの決定要素（Ｈ）と比較するステップを含む、請求項１から５のいずれか一項に記載の欠陥検出方法。
【請求項７】
前記決定要素（Ｈ）が、複数の特性値を複数の欠陥クラス（Ｃ１、Ｃ２）に事前に分類することによって定義される、請求項６に記載の欠陥検出方法。
【請求項８】
前記複数の特性値を取得するステップ（１０３）と前記欠陥を決定するステップ（１０５）との間に、前記取得した複数の特性値の中から少なくとも１つの関連値（Ｖａｌｉ’）を選択するステップ（１０４）があり、前記欠陥の決定が前記関連値（Ｖａｌｉ’）からなされる、請求項１から７のいずれか一項に記載の欠陥検出方法。
【請求項９】
前記電気化学デバイスから受け取った前記変数を予備処理するステップ（１０１）を含む、請求項１から８のいずれか一項に記載の欠陥検出方法。
【請求項１０】
前記予備処理するステップが、前記電気化学デバイスから受け取った前記変数の少なくとも１つの周波数成分を除去するステップを含む、請求項９に記載の欠陥検出方法。
【請求項１１】
請求項１から１０のいずれか一項に記載の欠陥検出方法のステップを実行するための命令を含む、電気通信ネットワークを介してダウンロード可能であり、かつ／またはコンピュータのメモリに記憶され、かつ／または前記コンピュータのリーダと協同するように意図された記憶媒体に記憶されるコンピュータプログラム。
【請求項１２】
電気化学デバイスのための欠陥検出デバイスであって、
少なくとも１つの変数（Ｓ）を前記電気化学デバイスから受け取り、ウェーブレット変換を含む数学的操作を行うことによって前記変数（Ｓ）から少なくとも１つの特性値（Ｖａｌｉ）を生成するように構成された処理モジュール（２０３）と、
前記処理モジュールから受け取った少なくとも１つの値から前記電気化学デバイスの少なくとも１つの欠陥を決定するように構成された決定モジュール（２０５）と
を含む、欠陥検出デバイス。
【請求項１３】
前記決定モジュールが、前記処理モジュールから受け取った前記特性値に応じて欠陥を示す変数（Ｄｉ）を供給するように構成された欠陥分類手段を含む、請求項１２に記載の欠陥検出デバイス。
【請求項１４】
前記欠陥分類手段がニューラルネットワークを利用する、請求項１３に記載の欠陥検出デバイス。
【請求項１５】
前記ウェーブレット変換が、値のスケール変数ａおよび値の平行移動変数ｂによって決まるウェーブレット係数Ｓａ，ｂへの離散変換である、前記変数（Ｓ）の複数のウェーブレット係数への分解を含む、請求項１２から１４のいずれか一項に記載の欠陥検出デバイス。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は電気化学デバイスの欠陥を検出する分野に関する。
【背景技術】
【０００２】
現在、化学エネルギーの電気エネルギーへの変換、または逆の変換によることを意味する「電気化学」デバイスとして知られる様々なデバイスがある。
【０００３】
このタイプのデバイスの第１のカテゴリは、電気エネルギーを電気デバイスに供給するかまたはそれを後の使用のために貯蔵する目的で化学エネルギーを電気エネルギーに変換するデバイスに関係する。そのようなデバイスの例はバッテリ、燃料電池、またはスーパーキャパシタである。
【０００４】
このタイプのデバイスの第２のカテゴリは、化学反応を行うか、または生成物もしくは反応物を分離するのに電気による様々な方法を使用するデバイスに関する。そのようなデバイスは、一般に、電着、放電機械加工、または電解浮選などの「電気化学的」方法を使用する。
【０００５】
これらの電気化学デバイスの信頼性および耐用寿命は様々な現象によって制限される。化学エネルギーを電気エネルギーに変換するデバイス、例えば燃料電池などでは、２つの主要な現象が性能の低下またはデバイスの完全な故障さえもたらす。
【０００６】
これらのデバイスの耐用寿命は、充電／放電のサイクルによって、あるいは度重なる遮断および起動または電力需要の変動による間欠的動作によって低下する。
【０００７】
さらに、例えば、使用される電気化学的方法のいくつかのパラメータの制御における故障などのいくつかの事象（反応物供給の中断、反応生成物および副次的生成物の不十分な管理）、媒体の毒作用、または構成要素もしくはモジュールの故障が、動作中に生じることがある。
【０００８】
これらの有害な現象は、それらを検出し、それらを修正できると思われる診断方法の使用を必要とする。
【０００９】
従来の電気化学デバイスでは、従来の診断方法は通常、いくつかのパラメータを知ることに基づき、それらのパラメータは、これらのシステムの外部または内部のものである場合があり、電気化学デバイス自体に挿入された内部センサなどの特定の計測器が必要となる。
【００１０】
そのような計測器は、高価であったり、センサ、特にデバイスの内部にあるものを取り付けるにはあまり適していない電気化学システムの形状のため実装するのが困難であったりすることが多いので、必ずしも望ましいとは限らない。
【００１１】
さらに、それら内部センサの組み込みのためには、デバイスが改造され、それにより、欠陥の可能性が増加し、診断が不正確になることがある。さらに、移動式用途で電気化学システムを使用する場合、デバイスのサイズはできるだけ小さくなければならず、したがって、診断計測器もそうでなければならず、そのため、従来の方法を使用するのは不可能である。
【００１２】
プロトン交換膜燃料電池（ＰＥＭＦＣ）の特定であるが排他的でない場合には、劣化機構を理解し、信頼性および耐用寿命を改善するために研究が行われてきている。いくつかの物理的モデルが開発されている。しかし、それらは、それらを使用するために測定することが困難であるか、または不可能でさえあるいくつかのパラメータを必要とする。さらに、それらの複雑さは、一般に、かなりの計算時間を必要とし、そのため、リアルタイム診断でそれらを適用するのは困難である。
【発明の概要】
【発明が解決しようとする課題】
【００１３】
本発明はこれらの欠点を克服することを目的とする。
【００１４】
したがって、本発明の目的は非侵入的方法で電気化学デバイスの欠陥を検出する方法を提供することである。
【００１５】
本発明の別の目的は最少の計測器を使用して欠陥を検出する方法を提供することである。
【００１６】
本発明のさらなる別の目的は、タイプ、形状、サイズ、または用途と無関係に様々なシステムで使用することができる、電気化学デバイスの欠陥を検出する一般的な方法を提供することである。
【００１７】
最後に、本発明の別の目的はリアルタイムで使用することができる、欠陥を検出する方法を提供することである。
【課題を解決するための手段】
【００１８】
このために、本発明は、電気化学デバイスの欠陥を検出する方法を提案し、その方法は、前記電気化学デバイスから受け取った少なくとも１つの変数から少なくとも１つの特性値を取得するステップと、この取得した値に基づいて前記電気化学デバイスの少なくとも１つの欠陥を決定するステップとを含み、ウェーブレット変換を含む数学的操作が、受け取った前記変数から特性値を取得するように行われる。
【００１９】
有利には、ウェーブレット変換は、取得した特性値がスケール変数ａおよび平行移動変数ｂによって決まる少なくとも１つのウェーブレット係数Ｓａ，ｂを含む離散ウェーブレット変換である。この離散化はウェーブレットへの分解に必要とされる計算時間を改善する。
【００２０】
特に、複数の特性値は、スケールレベル変数ｊがａ未満である１組のウェーブレット係数ｗｊ，ｐを複数の組のウェーブレット係数ｗｊ＋１，ｐに分解することによって取得される。
【００２１】
好ましくは、値のスケール変数ａはａ＝ａ０ｊを満足し、ここで、ｊはスケールレベルであり、ａ０はスケールパラメータであり、１組の特性値は、０から所与の分解レベルまでの値ｊごとに、各組のウェーブレット係数ｗｊ，ｐをａ０組のウェーブレット係数ｗｊ＋１，ｐに連続的に分解することによって取得される。これはウェーブレット分解の間に高レベルの細部をもたらす。
【００２２】
好ましくは、所与の分解レベルは最大分解レベルに対応し、その結果、最高レベルの細部がウェーブレット分解の間に取得される。
【００２３】
好ましい実施形態では、決定するステップは、特性値を、少なくとも１つの第１の欠陥クラスを第２の欠陥クラスから分離する少なくとも１つの決定要素と比較するステップを含む。有利には、決定要素は、複数の特性値を複数の欠陥クラスに事前に分類することによって定義される。
【００２４】
好ましくは、複数の特性値を取得するステップと欠陥を決定するステップとの間に、取得した複数の特性値の中から少なくとも１つの関連値（Ｖａｌｉ’）を選択するステップがあり、欠陥の決定が前記関連値からなされる。これにより計算時間が速くなる。
【００２５】
それは、この方法が電気化学デバイスから受け取った変数を予備処理するステップを含む場合に特に有利である。特に、この予備処理するステップは、必要とされる計算時間を最適化するために、電気化学デバイスから受け取った変数の少なくとも１つの周波数成分を除去するステップを含む。
【００２６】
本発明は、さらに、上述の方法のステップを実行するための命令を含むコンピュータプログラムを提案する。
【００２７】
本発明は、さらに、電気化学デバイスの欠陥を検出するためのデバイスを提案し、そのデバイスは、少なくとも１つの変数をこの電気化学デバイスから受け取り、ウェーブレット変換を含む数学的操作を行うことによってこの変数から少なくとも１つの特性値を生成するように構成された処理モジュール、ならびに処理モジュールから受け取った少なくとも１つの値から電気化学デバイスの少なくとも１つの欠陥を決定するように構成された決定モジュールを含む。
【００２８】
本発明の目的である方法、コンピュータプログラム、および欠陥検出デバイスは、以下の説明を読み、添付図面を検討することにより一層よく理解されるであろう。
【図面の簡単な説明】
【００２９】
【図１】本発明による、電気化学デバイスの欠陥を検出する方法のステップを示す図である。
【図２】離散ウェーブレット変換の使用に由来する第１のタイプのツリー構造、すなわち完全ツリーを示す図である。
【図３】離散ウェーブレット変換の使用に由来する第２のタイプのツリー構造、すなわち部分ツリーを示す図である。
【図４】事前分類法で定義されるようなマージン、サポートベクトル、および分離超平面の概念を示す図である。
【図５】燃料電池の特定の例における本発明による事前欠陥分類の一例を示す図である。
【図６】本発明による電気化学デバイス中の欠陥検出デバイスを概略的に示す図である。
【発明を実施するための形態】
【００３０】
まず、本発明の電気化学デバイスの欠陥を検出する方法のステップを示す図１を参照する。
【００３１】
ここで使用される「電気化学デバイス」という用語は、例えば、化学反応を達成するか、または生成物もしくは反応物を分離するために、化学エネルギーの変換によって電気エネルギーを生成し、それを供給する（直接または一時的にそれを貯蔵することによって）ことができる任意のデバイス、ならびに電気エネルギーの化学エネルギーへの変換を使用することができる任意のデバイスを包含することを思い起こされたい。
【００３２】
そのようなデバイスは、バッテリ、燃料電池、またはスーパーキャパシタからなることができる。または、そのようなデバイスは、電気めっき、放電機械加工、電気合成、電気精製、電気濃縮、または電解浮選のための電池などの電解槽からなることができる。そのようなデバイスは、さらに、電気透析器からなることができる。
【００３３】
電気化学デバイスのあり得る欠陥を検出するために、本発明の方法は、電気化学デバイスから受け取った変数Ｓに行われるいくつかの連続操作を含むことになり、それにより、発生源の内部へセンサを挿入する必要なく非侵入診断を行うことができるようになる。
【００３４】
発生源から受け取った変数Ｓは、電気化学デバイスの特性評価を可能にする任意のタイプの信号からなることができる。
【００３５】
電気を生成するデバイスの場合には、この変数Ｓは、単に、デバイスから出力として送出される電圧、電流、または電力などの任意の信号とすることができる。
【００３６】
化学変換を行うために電気を使用するデバイスの場合には、この変数Ｓは、そのようなデバイスに入力される指定のパラメータに対するデバイスの応答とすることができる。例えば、指定の電流が入力される場合、変数Ｓはこのデバイスの電圧応答とすることができる。逆に、指定の電圧が入力される場合、変数はこのデバイスの電流応答とすることができる。最後に、指定の電力が入力される場合、変数はデバイスの電流または電圧応答とすることができる。
【００３７】
以下で詳述する非限定の例では、この変数Ｓは指定の電流で動作するバッテリの端子で測定された出力電圧であるが、電圧または電力が指定される場合バッテリからの電流を、電圧または電流が指定される場合バッテリからの電力を、あるいは、任意の動作モードの場合、生成物もしくは反応物の圧力もしくは濃度、反応物もしくは生成物の流量、温度、またはこれらの変数の任意の時間的もしくは空間的変化を使用することを容易に考慮に入れることができる。
【００３８】
第１の主ステップ１０３の間、第１の処理が、分析されるべき電気化学デバイスから受け取った変数Ｓに適用され、電気化学デバイスの１つまたは複数の欠陥を特徴づける１つまたは複数の値Ｖａｌｉ（ここで１≦ｉ≦ｎ）が取得されることになる。
【００３９】
特に、変数Ｓが現実にはアナログである場合、取得した特性値Ｖａｌｉは後続のデジタル処理で使用することができるデジタル変数になることになる。
【００４０】
次に、これらの取得した特性値Ｖａｌｉから、電気化学デバイスの１つまたは複数の欠陥を第２の主ステップ１０５の間に決定することができる。
【００４１】
ステップ１０５は、定義された基準に従って特性Ｖａｌｉをグループ化することによって特性Ｖａｌｉがある程度組織化された構造に分割される場合に管理しないか、または１組の既に分類されたデータに基づいて管理することができる。
【００４２】
例えば、管理された場合にこれを行うために、取得した特性値Ｖａｌｉは、各々が電気化学デバイスの特定の状態、例えば、特定のタイプの欠陥が存在する状態に関連する一連の事前に分類された値と比較される。この比較から、電気化学デバイスの少なくとも１つのあり得る欠陥Ｄｉを導き出すことができる。
【００４３】
したがって、本発明の方法は、受け取った変数Ｓから値Ｖａｌｉを取得するために、第１の主ステップ１０３の間にウェーブレット変換を含む数学的操作を使用することによって最初に特徴付けられる。
【００４４】
ウェーブレットは、中心時間のまわりに局所化された限定期間の数学関数Ψである。その名称（ウェーブレット）はそのコンパクトでかつ振動する性質を反映している。任意の数学関数は、それが振動し、有限のエネルギーのものであり、０に等しい平均を有する性質を有する場合、ウェーブレットと考えることができる。
【００４５】
変数を分析する他の方法に対するウェーブレット解析の第１の利点は、「マザーウェーブレット」として使用可能な多くの関数があることである。
【００４６】
マザーウェーブレットとしてしばしば使用される関数はメキシカンハット関数である。その数式は以下のとおりである。
【００４７】
【数１】
2012530925.tif 000003
【００４８】
このマザーウェーブレットから、一群のウェーブレット（Ψａ，ｂ（ｔ））ａ，ｂは以下の式による時間平行移動および膨張（またはウェーブレット圧縮）によって定義される。
【００４９】
【数２】
2012530925.tif 000004
【００５０】
分解のための基底として働くこの一群のウェーブレットにより、所与の変数Ｓを分析することが可能になる。「平行移動変数」と呼ばれる変数ｂは時間局所化パラメータであり、一方、「スケール変数」と呼ばれる変数ａはスケール因子に対応する。大きいスケールは信号の全体図に対応し、小さいスケールは詳細の記述に対応する。ウェーブレット変換の使用により、変数はすべてのスケールで分析することができ、多重解像度分析をもたらす。
【００５１】
したがって、この信号に含まれる異なるスケールで生じる異なる現象（したがって、欠陥）に関する情報を取得することができる。分解の各レベルで、信号は異なるスケールで取得され、それにより、あるスケールの分解から次の（より細密な）スケールに進む場合に現象を局所化させることが可能になる。
【００５２】
ウェーブレット変換を行うために、スケール変数ａごとに、ウェーブレットは分析されるべき変数だけ時間軸の原点からシフトされ（平行移動変数ｂを変更することによって）、２つの間の一連の相関が計算される。
【００５３】
これらの相関の結果は１組の「ウェーブレット係数」Ｓａｂに対応し、ウェーブレット係数は、ウェーブレットの形状が分析されるべき変数の形状に近づくとき最大になり、以下の式を満足する。
【００５４】
【数３】
2012530925.tif 000005
【００５５】
ここで、＊は共役を示し、Ｒは実数の組である。次に、変数Ｓはこれらの係数Ｓａ，ｂによって記述され、次に、これらの係数はデバイスの欠陥を決定するための特性値Ｖａｌｉとして働くことができる。
【００５６】
これは、変数ａおよびｂが連続的に変化する場合、変数Ｓの連続ウェーブレット変換（ＣＷＴ）と呼ばれる。そのような連続ウェーブレット変換は信号のＳ（ｇ）の完全な記述を与えるが、その代償が、式（２）を実施することの困難さ、および結果として生じる高い冗長性である。
【００５７】
この困難を克服するために、変数ａおよびｂにいくつかの離散値を使用するように限定する離散ウェーブレット変換と呼ばれるウェーブレット変換のタイプを使用することが有利である。多くの用途では、ａおよびｂの値は以下のように定義され、
（４）            ａ＝ａ０ｊ、
ここで、ａ０＞１であり、
（５）            ｂ＝ｋ．ｂ０．ａ０ｊ、
ここで、ｂ０＞０であり、ｊおよびｋは整数である。
【００５８】
変数ｊおよびｋはそれぞれスケールレベルおよび平行移動レベルである。取得される結果は一連の離散値であり、これはウェーブレット系列分解と呼ばれる。全く例示的な目的のために、ここで選ばれた値はｂ０＝１およびａ０＝２であり、それらはスケール変数が連続の値１、２、４、８などを仮定するダイアディック分解に対応する。
【００５９】
この特定の場合には、信号を分解するために使用されるウェーブレットは次式の
【００６０】
【数４】
2012530925.tif 000006
【００６１】
を満足し、ウェーブレット係数は、
【００６２】
【数５】
2012530925.tif 000007
【００６３】
によって定義され、オリジナルの変数Ｓは、
【００６４】
【数６】
2012530925.tif 000008
【００６５】
に従った対応するウェーブレットに基づいて定義される。
【００６６】
信号処理の単純な観点から、２つのタイプの要素のみ、すなわち、解析されるべき変数と、それを解析またはフィルタ処理する関数とが考慮されるべきである。この観点から、信号の離散ウェーブレット変換は、この信号がフィルタバンクを通過させられることと見なすことができる。
【００６７】
所与のスケールにおいて、離散ウェーブレット変換は、前のスケールからの係数がａ０フィルタからなるバンクを通過させられることからなる。因子ａ０が２に等しい例では、ローパスフィルタが信号の粗い画像を与え、ハイパスフィルタが細部を示す。これらの２つのフィルタは相補的であり、一方によって除去された周波数は他方によって回復される。フィルタバンクの反復使用により、図２および３に示されるツリー構造がもたらされる。
【００６８】
スケールレベルｊごとに、係数の組ｗｊ，ｐが取得され、ここで、パラメータｐはツリーにおける位置を示し、０と２ｊ−１との間で変化し、１組の係数ｗｊ，ｐに対応する各ノードに対して左側にあるノードの数に等しい。それは周波数インデックスと見なすことができる。組ｗｊ，ｐは一連の係数Ｓｊ，ｋを含み、ここで、ｋは０から２Ｍ−ｊ−１まで変化し、パラメータＭは分解されるべき信号の分解の最高レベルに対応し、それは、例えば、この信号の長さが２の整数のべき乗である場合に、この信号の長さの自然対数に対応することができる。
【００６９】
したがって、固定のスケールレベルｊでは、係数の組ｗｊ，ｐは、
【００７０】
【数７】
2012530925.tif 000009
【００７１】
を満足する。
【００７２】
図２は、３つの連続するレベルのフィルタ処理が適用される変数Ｓを示す。スケールレベルｊごとに、０から最大分解レベルＭまで、連続する係数の組ｗｊ，ｐが取得され、それは、前のスケールレベルｊ−１の係数ｗｊ−１，ｐの組の各々へのローパスフィルタ（「Ｌｏ」で記号化された）およびハイパスフィルタ（「Ｈｉ」で記号化された）の適用に対応する。
【００７３】
ウェーブレットパケット変換として知られているそのような変換は、全部の分解レベルの各々において変数Ｓを完全に特徴付けることができるという意味で完全である。レベルｊでは、２ｊ組の係数（またはノード）が取得される。信号は各分解レベルで完全に示されるので、いくつかのレベルを有する完全な「ツリー」による変数Ｓのこの表現は冗長である。そのようなツリー構造を用いて、欠陥を識別するために所与の欠陥の「重要な」パケットのみを選択し、これらのパケットのみを使用することが可能である。
【００７４】
図３は、３つの連続するレベルのフィルタ処理を用いる部分ツリーを示す別の例を示す。この特別の場合には、ハイパスフィルタおよびローパスフィルタが、ゼロの「ツリー位置」変数ｐを有する（変数Ｓの「低周波数」成分を示す）係数ｗｊ，ｐの組にのみ適用され、この場合、係数ｗｊ，ｐはそれぞれスケールａ＝０、ａ＝２、およびａ＝４に対応する組ｗ０，０、ｗ１，０、およびｗ２，０である。
【００７５】
この場合、分解は、いかなるｊについても係数ｗｊ，０の組に限定される。変数Ｓの「高周波数」成分はもはや分解されず、したがって、低周波数成分ほど詳細には解析されない。
【００７６】
そのような分解は、取得されるべき係数が選択される場合、図２の分解ほど完全でないが、欠陥を決定するために分解されるべき変数Ｓの範囲が前もって分かっている場合に有用となり得る。この場合、この分解は計算時間がより速く、より効率的であり、特別なタイプの欠陥に直接集中することができる。
【００７７】
より一般的には、ウェーブレットまたはウェーブレットパケットへの分解の後に取得された係数によって、これらの信号の周波数コンテンツは利用することが可能になる。所与の欠陥と関連する分解信号のあらゆる変化が、離散ウェーブレット変換の１つまたは複数の分解レベルまたはウェーブレットパケット分解の１つまたは複数のパケットにおいて分かることになる。
【００７８】
そのような分解により、取得した係数ｗｊ，ｐの異なる組を使用して、基準を満たす事象のエネルギー、エントロピー、平均、最大、最小、標準偏差、数などのような１つまたは複数の特性を特性評価することが可能になる。次に、これらの特性（係数ｗｊ，ｐの組と同様に）は特性値Ｖａｌｉに対応することができ、それを使用して、あり得る欠陥を第２の主ステップ１０５の間に決定することになる。
【００７９】
これを行うために、取得した値Ｖａｌｉは事前に分類された一連の値と比較され、事前に分類された値の各々は電気化学デバイスの特定の状態、例えば、正常状態Ｄ０、またはあるタイプの欠陥に対応する状態Ｄｉに関連する。電気化学デバイスのあり得る欠陥はこの比較から導き出すことができる。
【００８０】
事前分類で使用された値は、各々が特別なタイプの欠陥に対応する１つまたは複数の欠陥クラスＣ１、Ｃ２に分類されるステップ１０３で取得された特性値と本質的に同等な値である。欠陥と値とのこの関連性は、分析されるべきデバイスの製造業者からのデータを使用して、または訓練およびフィードバックによって行うことができる。
【００８１】
値の欠陥クラスへの事前分類により、クラスの区別のために１つまたは複数の決定要素を定義することができることになる。取得した特性値Ｖａｌｉはステップ１０５で決定要素と比較されて、値Ｖａｌｉが欠陥クラスに属するかどうかが決定される。
【００８２】
そのような決定要素の形態は、考慮している次元の数に依存する。事前分類が単一の決定軸に関連してなされる場合、これらの決定要素は、値Ｖａｌｉが比較される閾値になることになる。
【００８３】
２つの値Ｖａｌｉ間の相関が観測される２次元分類では、決定要素は例えば直線になることになる。一般に、決定要素は、次元Ｎの空間の分離面に、例えば、線形の場合には分離超平面になることになる。
【００８４】
値の事前分類は様々な方法を使用して行うことができる。特に有利な１つの方法はサポートベクターマシンを使用することからなる。
【００８５】
サポートベクターマシン（すなわちＳＶＭ）は管理された学習に基づいた識別技法である。
【００８６】
これらのサポートベクターマシンは、高次元データで働くことができ、堅固な理論的基盤を有し、実際に良好な結果を与えるという利点を有する。さらに、適用モデルに関係なく、サポートベクターマシンの性能は他の分類法と同様またはそれよりも良好である。
【００８７】
サポートベクターマシンは以下の２つの重要な概念に基づく。
【００８８】
１）訓練集合の境界までの最小距離を最大にすることができる、クラスを分離する「最適の」境界の構築。これは、既知のアルゴリズムが適用される二次最適化問題としてその問題を定式化することによって行われる。
【００８９】
２）データが線形に分離可能でない場合を、適切な空間における線形分離のより簡単な場合にカーネル関数を使用して変形することができるように、サポートベクターマシンが入力データを示す空間を、より高い次元の空間、おそらく無限次元の空間に変換する。
【００９０】
この方法により、最初に変数を２つのクラスに分類することができる。しかし、より多くのクラスに分類するための拡張が存在する。
【００９１】
分離されるべき２つのクラス「−１」および「＋１」、ならびに以下の訓練集合の場合を考える。
【００９２】
【数８】
2012530925.tif 000010
【００９３】
２つの異なるクラスに属するデータを分離する最適超平面を構築するための２つの場合、すなわち、データが線形に分離可能であるか、またはデータが線形的に分離可能でない場合がある。
【００９４】
データが線形に分離可能である第１の場合、最適超平面Ｈは、
【００９５】
【数９】
2012530925.tif 000011
【００９６】
を満足し、それは、
【００９７】
【数１０】
2012530925.tif 000012
【００９８】
のように書くこともできる。
【００９９】
次に、超平面におけるポイントｘからの距離は、次式の
【０１００】
【数１１】
2012530925.tif 000013
【０１０１】
に従って、超平面上へのこのポイントの正射影によって与えられる。
【０１０２】
次に、２つのクラスの観測点（ｏｂｓｅｒｖａｔｉｏｎ）と超平面との間の最小距離に対応するマージンＭａを定義することができる。
【０１０３】
【数１２】
2012530925.tif 000014
【０１０４】
決定境界を意味する最適分離超平面ＨはこのマージンＭａを最大にするものであり、超平面に対する２つのクラスの距離の和を最大にし、したがって、式（７）の制約条件下で‖ｗ‖を最小にすることと同じものである。しかし、‖ｗ‖よりも‖ｗ‖２を最小にすることが容易である場合がある。
【０１０５】
したがって、最小化の問題は、以下の線形制約により二次関数を最小にする問題として定式化することができる。
【０１０６】
【数１３】
2012530925.tif 000015
【０１０７】
図４は、特定の２次元の場合のマージン、サポートベクトル、および分離超平面Ｈのこれらの概念を示す。
【０１０８】
この図において、値の２つのグループは、それぞれ欠陥Ｄ１および欠陥Ｄ２を示す２つのクラスＣ１およびＣ２に分類される。これらの２つのクラスＣ１とＣ２との間の境界を示す分離超平面Ｈは、「サポートベクトル」と呼ばれる各クラスＣｉのそれぞれの極限値Ｖｓ１およびＶｓ２に対して定義されたマージンＭａを最小にするものである。
【０１０９】
この分離超平面Ｈが事前訓練によって定義された後、決定ステップ１０５は、取得した特性値をこの分離超平面Ｈに対して位置決めすることからなり、それにより、特性値をクラスＣ１およびＣ２の一方に分類し、関連する欠陥Ｄｉを導き出すことができる。
【０１１０】
（１５）で明示した最小化問題の解決は、例えば、制約条件ごとにラグランジュ乗数を使用して行われる。この場合、以下の式が得られる。
【０１１１】
【数１４】
2012530925.tif 000016
【０１１２】
ラグランジュ関数は、ｗおよびｂに対して最小にされ、αに対して最大にされなければならない。
【０１１３】
キューン・タッカー条件下でラグランジュ関数の偏導関数を相殺することによって、以下の系が得られる。
【０１１４】
【数１５】
2012530925.tif 000017
【０１１５】
したがって、最適ラグランジュ乗数はαｉ＊である。αｉ＊がこの系の解である場合、
【０１１６】
【数１６】
2012530925.tif 000018
【０１１７】
であることが示される。
【０１１８】
最接近ポイントに対応するαｉ＊のみが非ゼロである（図４のサポートベクトルＶｓ１およびＶｓ２）。
【０１１９】
したがって、線形の場合の分離超平面はｗその最適値ｗ＊と置き換えることによって得られる。次に、以下の式が得られる。
【０１２０】
【数１７】
2012530925.tif 000019
【０１２１】
データが線形に分離可能でない他の場合、最適超平面は以下条件を満足するものである。
【０１２２】
１）適切に分類されたベクトルと最適超平面との間の最大距離は最大でなければならない、
【０１２３】
２）不適切に分類されたベクトルと最適超平面との間の距離は最小でなければならない。
【０１２４】
より形式的には、非線形変換φが入力ベクトルｘに適用される。その結果得られる空間は特徴空間と呼ばれる。この空間は、訓練集合のすべてのポイントについてｙｉｈ（ｘｉ）＞０であるとして
【０１２５】
【数１８】
2012530925.tif 000020
【０１２６】
を満足する超平面、すなわち、特徴空間における分離超平面を求めて探索される。
【０１２７】
次に、線形分離の上述の場合と同じ手順を使用して、最適化問題は以下のように書くことができる。
【０１２８】
【数１９】
2012530925.tif 000021
【０１２９】
したがって、Ｋ（ｘｉ，ｘｊ）＝φ（ｘｉ）．φ（ｘｊ），で定義されるカーネル関数を導入することによって、非線形の場合の分離超平面は次式を有する。
【０１３０】
【数２０】
2012530925.tif 000022
【０１３１】
カーネル関数を使用すると、以前に説明した線形の場合に移行する。線形、多項式、ガウシアン、ラプラシアンカーネルなどの多数のカーネル関数がある。
【０１３２】
ニューラルネットワークまたはｋ個の最近傍を使用する方法などの他の先行する分類法を使用することができる。
【０１３３】
有利な一実施形態では、第１のオプションの予備処理ステップ１０１が第１の主ステップ１０３の前に行われ、変数Ｓを前処理して特性検出方法が最適化される。
【０１３４】
前処理の一例は、一次ステップ１０３がウェーブレット分解を使用する場合には、重要でない周波数を有する成分を変数Ｓから除去することからなる。これにより、次に重要成分のみが分解されるのでこの分解が最適化される。
【０１３５】
これを行うために、閾値パラメータとして働くカットオフ周波数を有するフィルタを使用することができる。閾値の値は、例えば、フィードバック、システムの知見、または重要周波数帯域から経験的に決定される。
【０１３６】
特に有利な実施形態では、第１の主ステップ１０３の間に取得されたものの中からいくつかの最適値Ｖａｌｉを選択するために第１の主ステップ１０３と欠陥決定の第２の主ステップ１０５との間に最適化ステップ１０４がある。
【０１３７】
実際には、欠陥決定ステップで値Ｖａｌｉを直接使用するのは、ステップ１０３の間に生成された特性値Ｖａｌｉの組が非常に大きく、ノイズまたは望ましくない成分で汚染されている場合、問題を引き起こすことがある。
【０１３８】
さらに、特性値Ｖａｌｉはしばしば冗長な情報を含み、したがって、処理をすべて必要とするわけではない。
【０１３９】
さらに、ステップ１０５の決定プロセスを最適化するには、選択ステップ１０４の間にできるだけ処理されるべき値Ｖａｌｉの数を減少させて、決定ステップ１０５にとって最良と考えられる最も関連する値Ｖａｌ’ｉ（ここで、１≦ｉ≦ｍ、ｍ＜ｎ）のみを保持することが望ましい。これは電気化学デバイスの診断の頑強性の改善および計算時間の低減に寄与する。
【０１４０】
これを行うために、１つの特定の実施形態は、最良ウェーブレット基底を選択する方法を使用することができる。この方法は、「最良基底」と呼ばれる基底を選択する基準を使用することに基づく。
【０１４１】
この方法は以下の２つのステップ、すなわち、
−  ステップ１０３においてウェーブレットへの分解の間に取得された特性値Ｖａｌｉの組に選んだ基準を適用するステップと、
−  「重要でない」と考えられる特性値Ｖａｌｉを除去するために、選んだ基準に応じて、重要性の増加または減少の順序によって前のステップの間に見いだされた特性を並び替えるステップと
を含む。
【０１４２】
次に、残りの特性値Ｖａｌｉは決定ステップ１０５の間に使用できることになる。
【０１４３】
検出のための最適基底の一例は、様々な周波数と時間情報との間の分離可能性を最大にする基底である。決定のための最適基底は、様々な欠陥のクラス間の分離可能性、言い換えれば区別を最大にする基底である。
【０１４４】
この状況において、いくつかの基準を定義することができる。
【０１４５】
１）第１の例では、以下の式、
【０１４６】
【数２１】
2012530925.tif 000023
【０１４７】
に従って、２つのシーケンスｘとｙとの時間−周波数エネルギー分布の間の距離を測定することからなる「クロスエントロピー」が使用される。
【０１４８】
この値は、２つの異なるクラスを示す分布ｘｉとｙｉとの間のカルバック・ライブラー発散に対応する。
【０１４９】
次に、
【０１５０】
【数２２】
2012530925.tif 000024
【０１５１】
のように最適化されるべき基準Ｄ（ｘ，ｙ）を定義することができる。
【０１５２】
この例では、訓練集合の各クラスは、まず、各ノードがクラスの要素に対して係数の２乗の平均シーケンスを含むツリーによって示される。
【０１５３】
上述で定義された基準は２値であるので、すべてのクラスにペアで適用され、最終基準は結果として生じる２値の基準の合計である。
【０１５４】
２）第２の例では、ステップ１０５の間に様々なクラスを分離するのに向けてステップ１０３の間に取得された係数またはパケットのキャパシティ（ｃａｐａｃｉｔｙ）を最大にするために、ここでの基準は、互いに最も遠く離れているクラス間の分散を意味する「クラス間慣性」を最大にするとともに、互いにできるだけ近いクラスの分散を意味する「クラス内慣性」を最小にすることである。したがって、基準は、全慣性に対するクラス内慣性の比からなることができる。クラス間慣性に対するクラス内慣性の比の使用を考えることもできる。
【０１５５】
欠陥のｋ個のクラスＣ１、Ｃ２、…、Ｃｋと、それぞれｎ１、ｎ２、…、ｎｋ要素を含むそれぞれの重心ｇ１、ｇ２、…、ｇｋとがある一例では、全ポイントクラウドの重心はｇで表すことができる。
【０１５６】
次に、クラス間慣性は以下の式で定義される。
【０１５７】
【数２３】
2012530925.tif 000025
【０１５８】
ここで、ｄは定義された距離、例えばユークリッド距離である。
【０１５９】
クラス内慣性は以下の式で定義される。
【０１６０】
【数２４】
2012530925.tif 000026
【０１６１】
全慣性に関しては、以下の式で定義される。
【０１６２】
【数２５】
2012530925.tif 000027
【０１６３】
したがって、最終基準Ｒは、例えば、以下の関係で与えられる。
【０１６４】
【数２６】
2012530925.tif 000028
【０１６５】
分離されるべきクラスは事前に定義される。したがって、すべての欠陥クラス間を同時に区別するか、または２つずつクラスを分離するか、または所与のクラスを他のクラスすべてから分離することができる。
【０１６６】
別の特定の実施形態では、ステップ１９４における次元数の低減は特異値分解を使用する。
【０１６７】
ｍ行ｎ列の行列Ｍを特異値に分解するのは、それを
【０１６８】
【数２７】
2012530925.tif 000029
【０１６９】
の形態で書くことと同じであり、ここで、Σは、例えば、減少順に行列Ｍの特異値λｉ，ｉ＝１，…ｍを含む対角行列であり、
【０１７０】
【数２８】
2012530925.tif 000030
【０１７１】
であることを思い出されたい。
【０１７２】
２つの他の行列ＵおよびＶは特異値λｉ，ｉ＝１，…ｍに対応する特異ベクトル（右および左特異ベクトル）を含む。
【０１７３】
特異値は慣性または代表性の程度を反映すると解釈され、特異ベクトルは、初期データ（行列Ｍ）の変化が最も高い軸である。特異値が減少順に順序づけられる場合、最終値はデータの最小変化を含むものである。したがって、ｍからｐの特異値（ｐ＜ｎ）の変化に基づく次元数低減では、ｍ−ｐ＋１の最終特異値λｉに含まれる情報は無視できることを仮定する。
【０１７４】
例えば、主成分分析などの次元数低減のための他の方法を考えることができる。
【０１７５】
次に、図５は、燃料電池の非限定の場合の本発明による欠陥分類の一例を示す。
【０１７６】
この図５において、１組の特性値が２つの別個のウェーブレットパケットの関数としてグラフに示されている。
【０１７７】
これらの特性値の位置は、訓練によって、分析されるべき燃料電池の特定の動作状態と関連づけられている。
【０１７８】
グラフの中心に位置する特性値の第１のグループは、燃料電池の正常動作状態に対応する特性値のクラスＣ０を定義する。
【０１７９】
グラフに左側に位置する特性値の第２のグループは、燃料電池がドライアウト欠陥を有する異常動作状態に対応する特性値のクラスＣ１を定義する。
【０１８０】
最後に、グラフに右側に位置する特性値の第３のグループは、燃料電池がフラディング欠陥を有する異常動作状態に対応する特性値のクラスＣ２を定義する。
【０１８１】
これらのクラスＣ０、Ｃ１、Ｃ２は、訓練によって、および当該の様々な状態を有する電池の特性値の測定によって定義される。例えば、上述のように計算された決定要素によって定義されたこれらのクラス値の境界は、モジュール２０５に格納され、それぞれ正常状態、ドライアウト欠陥による異常状態、およびフラディング欠陥による異常状態を示す状態変数Ｄ０、Ｄ１、Ｄ２に関連づけられる。
【０１８２】
新しい診断が、訓練中に管理された分類で使用されたものと同じタイプの燃料電池で行われる場合、同じウェーブレットパケットが観測され、特性値は区域Ｚ０からＺ２のうちの１つに位置づけられる。特定の状態を示す状態変数Ｄ０からＤ２は、測定された特性値が位置する区域の関数として決定モジュール２０５によって生成されることになる。
【０１８３】
最後に、図６は、本発明の電気化学デバイスの欠陥を検出するためのデバイスの略図を示す。
【０１８４】
この図において、上述で示したタイプのいずれの電気化学デバイス２００も変数Ｓを検出デバイス２０１に供給する。
【０１８５】
この検出デバイス２０１は、決定モジュール２０５に接続された処理モジュール２０３を含み、決定モジュール２０５はそれ自体インターフェースモジュール２０７に接続される。
【０１８６】
処理モジュール２０３は、第１の主変換ステップ１０３ならびに上述のように関連特性値Ｖａｌｉ’の予備処理および選択のあり得るオプションのステップ１０１および１０４を行うように構成される。
【０１８７】
そのような処理モジュール２０３は、この目的のために、デジタル値を使用して計算を行うか、またはコンピュータプログラムを実行することができるプロセッサ、マイクロプロセッサ、または他の構成要素、例えば集積回路からなることができる。
【０１８８】
分析されるべき変数Ｓが現実にはアナログである場合、処理モジュール２０３は、処理することができるデジタル値に変数Ｓを変換するためのアナログ−デジタル変換器を含むことができる。
【０１８９】
処理モジュール２０３は、１つまたは複数の値Ｖａｌｉを取得する主ステップ１０３を実行した後、前記値を欠陥決定モジュール２０５に供給する。オプションの選択ステップ１０４が、さらに処理モジュール２０３によって行われる特定の場合には、関連特性値Ｖａｌｉ’が欠陥決定モジュール２０５に供給される。
【０１９０】
図６の例では、値Ｖａｌｉがいくつかの並列接続によって供給されるように示されているが、これらの値が直列的に送出される場合には単一接続を使用することができる。第１の並列の実施形態は値をより迅速に移送するが、第２の直列の実施形態は簡単であり、モジュール２０３と２０５との間の接続のコストを低減する。
【０１９１】
決定モジュール２０５は、処理モジュール２０５から受け取った値Ｖａｌｉから電気化学デバイスの１つまたは複数の特性Ｄｉを決定するように構成される。これを達成するために、決定モジュール２０５は、特性がこれらの値の関数として分類される分類手段を含むことができる。
【０１９２】
そのような分類手段は、例えば、ファジー論理または入力として受け取った値に基づいて異なる欠陥を分類するように訓練されたニューラルネットワークに基づく方法を使用することができる。
【０１９３】
この分類手段は、サポートベクターマシン、主成分分析、またはｋ個の最近傍を決定することなどの統計的な方法を使用することもできる。
【０１９４】
したがって、いくつかの値Ｖａｌｉに応じて、決定モジュール２０５は、電気化学デバイスの特性を示す１つまたは複数の変数Ｄｉ、例えば正常状態（Ｄ０）あるいは１つまたは複数の欠陥の特性を出力するのにその分類手段を利用することになる。
【０１９５】
次に、これらの値Ｄｉはインターフェースモジュール２０７によって受け取られ、インターフェースモジュール２０７は、電気化学デバイスの作動状態を、受け取った変数Ｄｉに応じて、検出デバイス２０１のユーザに、または検出デバイスから下流に位置する制御システムに示すことになる。
【０１９６】
これは、表示装置（変数Ｄｉに応じて状態または欠陥のタイプを明示することができる）、警報音、または分析されるべき電気化学デバイスの正常または異常作動状態をユーザまたは下流に配置された調整および制御システムに通知する任意の他の信号として行うことができる。
【０１９７】
もちろん、本発明は、上述で説明および示した例の特定の詳細に限定されず、本発明の範囲から逸脱することなく本発明から他の実施形態を考えることができる。
【図１】
2012530925.tif 000031
【図２】
2012530925.tif 000032
【図３】
2012530925.tif 000033
【図４】
2012530925.tif 000034
【図５】
2012530925.tif 000035
【図６】
2012530925.tif 000036
【国際調査報告】
2012530925.tif 000037
2012530925.tif 000038
2012530925.tif 000039
2012530925.tif 000040
2012530925.tif 000041
2012530925.tif 000042