(57)【要約】【課題】二次電池のインピーダンスを測定する。【解決手段】本発明による評価方法は、（１）二次電池に定電流Ｉｃを印加して、前記二次電池を充電または放電するステップを有し、さらに、前記（１）のステップ中に、（２）任意のタイミングで、前記二次電池に電流入力信号Ｉ（ｔ）を印加し、電圧応答信号Ｖ（ｔ）を取得するステップであって、前記電流入力信号Ｉ（ｔ）は、前記定電流Ｉｃに矩形波Ｉｐ（ｔ）を重畳することにより得られるステップと、（３）電流入力信号Ｉ（ｔ）および前記電圧応答信号Ｖ（ｔ）をフーリエ変換して、それぞれ、第１の周波数領域データＰＩ（ω）および第２の周波数領域データＰＶ（ω）を取得するステップと、（４）前記第１および第２の周波数領域データから、前記二次電池の所望の充電深度または放電深度でのインピーダンススペクトルを計算するステップと、を有する。【選択図】図１
2014238948.tif 000002
【特許請求の範囲】
【請求項１】
二次電池のインピーダンスの評価方法であって、
（１）二次電池に定電流Ｉｃを印加して、前記二次電池を充電または放電するステップを有し、
さらに、前記（１）のステップ中に、
（２）任意のタイミングで、前記二次電池に電流入力信号Ｉ（ｔ）を印加し、電圧応答信号Ｖ（ｔ）を取得するステップであって、前記電流入力信号Ｉ（ｔ）は、前記定電流Ｉｃに矩形波Ｉｐ（ｔ）を重畳することにより得られるステップと、
（３）電流入力信号Ｉ（ｔ）をフーリエ変換して、第１の周波数領域データＰＩ（ω）を取得し、前記電圧応答信号Ｖ（ｔ）をフーリエ変換して、第２の周波数領域データＰＶ（ω）を取得するステップと、
（４）前記第１および第２の周波数領域データから、以下の式により、前記二次電池の所望の充電深度または放電深度でのインピーダンススペクトルを計算するステップと、
Ｚ（ω）＝ＰＶ（ω）／ＰＩ（ω）      （１）式
を有することを特徴とする評価方法。
【請求項２】
当該評価方法は、前記二次電池に過充電または過放電が生じないように実施されることを特徴とする請求項１に記載の評価方法。
【請求項３】
前記定電流Ｉｃが正の場合、前記矩形波Ｉｐ（ｔ）は、前記重畳の間、前記電流入力信号Ｉ（ｔ）が前記定電流Ｉｃを上回らないように選定され、
前記定電流Ｉｃが負の場合、前記矩形波Ｉｐ（ｔ）は、前記重畳の間、前記電流入力信号Ｉ（ｔ）が前記定電流Ｉｃを下回らないように選定されることを特徴とする請求項１または２に記載の評価方法。
【請求項４】
前記矩形波Ｉｐ（ｔ）は、前記定電流Ｉｃの印加を一時的に停止することにより形成されることを特徴とする請求項１乃至３のいずれか一つに記載の評価方法。
【請求項５】
前記矩形波Ｉｐ（ｔ）は、ダブルパルスであることを特徴とする請求項１乃至３のいずれか一つに記載の評価方法。
【請求項６】
前記矩形波Ｉｐ（ｔ）は、単位波形の複数の繰り返しを有することを特徴とする請求項１乃至５のいずれか一つに記載の評価方法。
【請求項７】
さらに、
（５）前記インピーダンススペクトルＺ（ω）から得られる前記二次電池の一つの特性パラメータを、所定の値と比較するステップ
を有することを特徴とする請求項１乃至６のいずれか一つに記載の評価方法。
【請求項８】
前記所定の値は、前記二次電池の使用前に得られたインピーダンススペクトルＺ０（ω）から得られる前記二次電池の一つの特性パラメータ、または前記二次電池の一つの特性パラメータの閾値であることを特徴とする請求項７に記載の評価方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、二次電池のインピーダンスの評価方法に関する。
【背景技術】
【０００２】
充電による繰り返し使用が可能な蓄電池、すなわち二次電池は、各種用途において広く使用されている。特に、リチウムイオン二次電池は、小型、大容量で起電力が高い等の特徴があり、例えば、通信機器（例えば携帯電話、スマートフォン）、ＡＶ機器（例えば携帯型ゲーム機）、情報機器（例えばパーソナルコンピュータ）、および自動車（例えばハイブリッド車、電気自動車）など、幅広い分野で広く使用されている。
【発明の概要】
【発明が解決しようとする課題】
【０００３】
二次電池は、充電／放電を繰り返すうちに特性が劣化する。また、劣化した二次電池を使用し続けると、短絡や発熱などにより、大きな事故につながる場合がある。このため、二次電池に対しては、開発段階または定期的に、充放電曲線などの電池特性を評価することの重要性が指摘されている。また、最近では、電池特性として、単なる充放電曲線に留まらず、二次電池の内部抵抗（例えばインピーダンス）を評価することも検討されている。
【０００４】
ここで、各種電池特性の評価は、二次電池の充電過程中または放電過程中に、一度に実施することが効率的である。
【０００５】
しかしながら、これまで、充電特性（または放電特性）と内部抵抗など、複数の電池特性を一度に合わせて評価する方法はあまり知られていない。このため、電池の内部抵抗を評価する際には、いったん二次電池の充電または放電を中断してから、測定を実施しなければならないという問題がある。
【０００６】
また、この場合、内部抵抗の評価結果のばらつきを少なくし、測定精度を高めるためには、二次電池の充電または放電を停止してから二次電池が安定するまで、十分に時間を空けてから、内部抵抗の測定を実施する必要がある（例えば３０分〜１時間など）。しかしながら、このような対応では、二次電池の充電または放電の完了までに、相当の時間がかかってしまう。これは、例えば、電気自動車など、充電に長い時間をかけることが難しい場合に、大きな問題となり得る。
【０００７】
また、前述のような方法で測定された内部抵抗の測定結果は、二次電池のいわば静的な状態、すなわち平衡状態での挙動に相当する。従って、得られた結果は、二次電池が実際に使用される状態（充電過程または放電過程）における内部抵抗を正確に反映していない可能性がある。
【０００８】
このような背景から、迅速かつ正確に、二次電池の内部抵抗を評価することが可能な方法が要望されている。
【０００９】
本発明は、このような背景に鑑みなされたものであり、本発明では、従来に比べてより迅速に、実際の使用環境に近い状態で、二次電池のインピーダンスを測定することが可能な評価方法を提供することを目的とする。
【課題を解決するための手段】
【００１０】
本発明では、二次電池のインピーダンスの評価方法であって、
（１）二次電池に定電流Ｉｃを印加して、前記二次電池を充電または放電するステップを有し、
さらに、前記（１）のステップ中に、
（２）任意のタイミングで、前記二次電池に電流入力信号Ｉ（ｔ）を印加し、電圧応答信号Ｖ（ｔ）を取得するステップであって、前記電流入力信号Ｉ（ｔ）は、前記定電流Ｉｃに矩形波Ｉｐ（ｔ）を重畳することにより得られるステップと、
（３）電流入力信号Ｉ（ｔ）をフーリエ変換して、第１の周波数領域データＰＩ（ω）を取得し、前記電圧応答信号Ｖ（ｔ）をフーリエ変換して、第２の周波数領域データＰＶ（ω）を取得するステップと、
（４）前記第１および第２の周波数領域データから、以下の式により、前記二次電池の所望の充電深度または放電深度でのインピーダンススペクトルを計算するステップと、
Ｚ（ω）＝ＰＶ（ω）／ＰＩ（ω）      （１）式
を有することを特徴とする評価方法が提供される。
【００１１】
ここで、本発明による評価方法は、前記二次電池に過充電または過放電が生じないように実施されても良い。
【００１２】
また、本発明による評価方法において、前記定電流Ｉｃが正の場合、前記矩形波Ｉｐ（ｔ）は、前記重畳の間、前記電流入力信号Ｉ（ｔ）が前記定電流Ｉｃを上回らないように選定され、
前記定電流Ｉｃが負の場合、前記矩形波Ｉｐ（ｔ）は、前記重畳の間、前記電流入力信号Ｉ（ｔ）が前記定電流Ｉｃを下回らないように選定されても良い。
【００１３】
また、本発明による評価方法において、前記矩形波Ｉｐ（ｔ）は、前記定電流Ｉｃの印加を一時的に停止することにより形成されても良い。
【００１４】
また、本発明による評価方法において、前記矩形波Ｉｐ（ｔ）は、ダブルパルスであっても良い。
【００１５】
また、本発明による評価方法において、前記矩形波Ｉｐ（ｔ）は、単位波形の複数の繰り返しを有しても良い。
【００１６】
また、本発明による評価方法は、さらに、
（５）前記インピーダンススペクトルＺ（ω）から得られる前記二次電池の一つの特性パラメータを、所定の値と比較するステップ
を有しても良い。
【００１７】
この場合、前記所定の値は、前記二次電池の使用前に得られたインピーダンススペクトルＺ０（ω）から得られる前記二次電池の一つの特性パラメータ、または前記二次電池の一つの特性パラメータの閾値であっても良い。
【発明の効果】
【００１８】
本発明では、従来に比べてより迅速に、実際の使用環境に近い状態で、二次電池のインピーダンスを測定することが可能な評価方法を提供することができる。
【図面の簡単な説明】
【００１９】
【図１】本発明による第１の評価方法のフローの一例を模式的に示した図である。
【図２】二次電池に定電流Ｉｃを印加した際の電流Ｉおよび電圧Ｖの経時変化を模式的に示した図である。
【図３】本発明による第１の評価方法の一工程において使用されるデータおよび得られるデータの一例を模式的に示した図である。
【図４】本発明による第１の評価方法において使用され得る矩形波の一例を模式的に示した図である。
【図５】本発明による第１の評価方法において使用され得る矩形波の一例を模式的に示した図である。
【図６】本発明による第１の評価方法において使用され得る矩形波の一例を模式的に示した図である。
【図７】本発明による第１の評価方法において使用され得る矩形波の一例を模式的に示した図である。
【図８】本発明による第１の評価方法において使用され得る矩形波の一例を模式的に示した図である。
【図９】本発明による第１の評価方法の一工程において得られるデータの一例を模式的に示した図である。
【図１０】本発明による第１の評価方法の一工程において得られるデータの一例を模式的に示した図である。
【図１１】本発明による第２の評価方法のフローの一例を模式的に示した図である。
【図１２】二次電池に対して定電流Ｉｃを通電させた際の電流Ｉおよび電圧Ｖの経時変化を模式的に示した図である。
【図１３】本発明による第２の評価方法の一工程において使用されるデータおよび得られるデータの一例を模式的に示した図である。
【図１４】本発明による第２の評価方法の一工程において得られるデータの一例を模式的に示した図である。
【図１５】本発明による第２の評価方法の一工程において得られるデータの一例を模式的に示した図である。
【図１６】実施例において想定した二次電池の等価回路を示した図である。
【図１７】二次電池に印加される電流入力信号Ｉ（ｔ）（＞０）の波形を示したグラフである。
【図１８】二次電池に図１７に示した電流入力信号Ｉ（ｔ）を入力した際に得られる電圧応答信号Ｖ（ｔ）の波形を示したグラフである。
【図１９】電流入力信号Ｉ（ｔ）のフーリエ変換によって得られた第１の周波数データＰＩ（ω）（図１９（ａ））、および電圧応答信号Ｖ（ｔ）のフーリエ変換によって得られた第２の周波数データＰＶ（ω）（図１９（ｂ））である。
【図２０】第１の周波数データＰＩ（ω）および第２の周波数データＰＶ（ω）から計算されたインピーダンススペクトル（ナイキストプロット）である。
【発明を実施するための形態】
【００２０】
以下、図面を参照して、本発明について詳しく説明する。
【００２１】
本発明では、二次電池のインピーダンスの評価方法であって、
（１）二次電池に定電流Ｉｃを印加して、前記二次電池を充電または放電するステップを有し、
さらに、前記（１）のステップ中に、
（２）任意のタイミングで、前記二次電池に電流入力信号Ｉ（ｔ）を印加し、電圧応答信号Ｖ（ｔ）を取得するステップであって、前記電流入力信号Ｉ（ｔ）は、前記定電流Ｉｃに矩形波Ｉｐ（ｔ）を重畳することにより得られるステップと、
（３）電流入力信号Ｉ（ｔ）をフーリエ変換して、第１の周波数領域データＰＩ（ω）を取得し、前記電圧応答信号Ｖ（ｔ）をフーリエ変換して、第２の周波数領域データＰＶ（ω）を取得するステップと、
（４）前記第１および第２の周波数領域データから、以下の式により、前記二次電池の所望の充電深度または放電深度でのインピーダンススペクトルを計算するステップと、
Ｚ（ω）＝ＰＶ（ω）／ＰＩ（ω）      （１）式
を有することを特徴とする評価方法が提供される。
【００２２】
なお、各周波数ωは、周波数ｆとω＝２πｆの関係にある。従って、ＰＶ（ω）およびＰＩ（ω）は、それぞれ、ＰＶ（ｆ）およびＰＩ（ｆ）とも表記でき、その場合、（１）式で計算されるＺ（ω）は、Ｚ（ｆ）と表記できることに留意する必要がある。
【００２３】
前述のように、これまで、二次電池の充電または放電過程中に、二次電池の内部抵抗を同時に測定することは難しく、このため、二次電池の内部抵抗測定は、充電または放電をいったん中断させてから、実施する必要があった。
【００２４】
しかしながら、このような内部抵抗の測定方法では、内部抵抗の測定前に、二次電池の電圧状態が安定するまで十分に時間を空ける必要がある（例えば３０分〜１時間など）。
このため、二次電池の充電または放電をいったん中断させて内部抵抗の測定を実施すると、二次電池の充電または放電が完了するまでに、相当の時間がかかってしまうという問題がある。
【００２５】
また、このような方法で測定された内部抵抗の測定結果は、二次電池の平衡状態での挙動に相当する。従って、得られた結果は、二次電池が実際に使用される状態（充電過程または放電過程）における内部抵抗を正確に反映していない可能性がある。
【００２６】
これに対して、本発明では、二次電池に定電流Ｉｃを印加して、二次電池を充電または放電する過程中に、いわばｉｎ−ｓｉｔｕで、二次電池のインピーダンスが測定される。このため、本発明では、例えば、二次電池の充電曲線または放電曲線の測定と、内部抵抗の測定を、一度に行うことができる。
【００２７】
従って、本発明では、内部抵抗の測定を実施することにより、二次電池の充電が完了するまでに必要な時間が延びてしまうという問題を有意に抑制することができる。また、本発明では、内部抵抗の測定を実施するため、二次電池の放電、すなわち使用を中断する必要性も有意に回避することができる。
【００２８】
また、本発明では、インピーダンス測定は、充電過程中または放電過程中の二次電池に対して実施される。このため、従来のような平衡状態で得られる内部抵抗値とは異なり、測定された結果は、被測定対象である二次電池の実際の使用状態（充電過程または放電過程）における挙動を、より正確に反映したものとなる。従って、本発明では、二次電池の内部抵抗をより正確に評価することができる。
【００２９】
なお、本発明では、二次電池の内部抵抗を表す指標として、インピーダンススペクトルを採用する。
【００３０】
インピーダンススペクトルを内部抵抗の測定指標とすることにより、二次電池の内部抵抗を非破壊で測定することができる。また、二次電池の内部抵抗を、電解液抵抗、電荷移動抵抗、および界面抵抗など、各パラメータ特性毎に分離評価することが可能になる。また、内部抵抗を、アノード側成分およびカソード側成分に分離して評価することも可能となる。
【００３１】
なお、所望の充電深度または放電深度でのインピーダンスを求めることは重要である。ここで、充電深度とは、完全に放電した二次電池を完全に充電するために必要な電気量を１００％としたとき、二次電池の充電過程において、実際に二次電池に流した電気量の割合を意味する。一方、放電深度とは、完全に充電した二次電池を完全に放電するために必要な電気量を１００％としたとき、二次電池の放電過程において、実際に二次電池に流した電気量の割合を意味する。二次電池のインピーダンスは、充電深度および放電深度により変化するため、二次電池の性能評価またはモニタリングのためには、所定の充電深度または放電深度でのインピーダンスを求め、インピーダンススペクトルから得られるパラメータを相対的に評価することが重要である。
【００３２】
（本発明の一実施例による二次電池のインピーダンスの評価方法について）
次に、図１〜図１０を参照して、本発明の一実施例による二次電池のインピーダンスの評価方法（以下、単に「（本発明による）第１の評価方法」と称する）について説明する。
【００３３】
図１には、本発明による第１の評価方法の模式的なフロー図を示す。また、図２、図３、図９、ならびに図１０には、本発明による第１の評価方法の一工程において使用され、および／または得られるデータの一例を模式的に示す。さらに、図４〜図８には、本発明による第１の評価方法において使用され得る矩形波の一例を模式的に示す。
【００３４】
図１に示すように、本発明による第１の評価方法は、
（１）二次電池に定電流Ｉｃを印加して、前記二次電池を充電するステップ（Ｓ１１０）と、
（２）前記（１）中の任意のタイミングで、前記二次電池に電流入力信号Ｉ（ｔ）を印加し、電圧応答信号Ｖ（ｔ）を取得するステップであって、前記電流入力信号Ｉ（ｔ）は、前記定電流Ｉｃに矩形波Ｉｐ（ｔ）を重畳することにより得られるステップ（Ｓ１２０）と、
（３）電流入力信号Ｉ（ｔ）をフーリエ変換して、第１の周波数領域データＰＩ（ω）を取得し、前記電圧応答信号Ｖ（ｔ）をフーリエ変換して、第２の周波数領域データＰＶ（ω）を取得するステップ（Ｓ１３０）と、
（４）前記第１および第２の周波数領域データから、以下の式により、前記二次電池の所望の充電深度でのインピーダンススペクトルを計算するステップ（Ｓ１４０）と、
Ｚ（ω）＝ＰＶ（ω）／ＰＩ（ω）      （１）式
を有する。
【００３５】
以下、各ステップについて詳しく説明する。
【００３６】
（ステップＳ１１０）
まず、評価対象となる二次電池が準備される。二次電池の種類は、特に限られない。二次電池は、例えば、鉛蓄電池、ニッケルカドミウム二次電池、ニッケル亜鉛二次電池、ニッケル水素二次電池、リチウム空気電池、リチウムイオンキャパシタ、およびリチウムイオン二次電池等であっても良い。なお、リチウムイオン二次電池には、リチウムポリマー二次電池、リチウム硫黄二次電池等が含まれる。
【００３７】
なお、二次電池は、被給電装置から取り外した状態であっても良く、あるいは被給電装置に装着した状態のままであっても良い。
【００３８】
次に、二次電池に対して、定電流Ｉｃが印加され、二次電池が充電される。
【００３９】
図２には、二次電池に対して、定電流Ｉｃを印加した際の電流Ｉおよび電圧Ｖの経時変化を模式的に示す。図２（ａ）は、二次電池の電流変化挙動を示し、図２（ｂ）は、二次電池の電圧変化挙動を示している。
【００４０】
なお、この図２の例では、定電流Ｉｃは正の値であると仮定している（Ｉｃ＞０）。
【００４１】
この場合、定電流Ｉｃの印加により二次電池が充電されるため、図２（ｂ）に示すように、二次電池の電圧Ｖは、時間とともに徐々に上昇する。
【００４２】
定電流Ｉｃの値は、特に限られず、使用される二次電池の種類および容量等に応じて、適宜決定される。例えば、携帯電話等に使用される小型のリチウムイオン二次電池の場合、定電流Ｉｃの値は、数ｍＡ〜数Ａの範囲であっても良い。
【００４３】
また、このステップＳ１１０の際に、図２（ｂ）に示すような二次電池の充電曲線（電圧−時間曲線）が取得されても良い。なお、二次電池の充電は、過充電とならない時間で完了される。
【００４４】
（ステップＳ１２０）
前述のステップＳ１１０の期間中の任意のタイミングにおいて、定電流Ｉｃに矩形波Ｉｐ（ｔ）が重畳される。すなわち、二次電池に対して、定電流Ｉｃに矩形波Ｉｐ（ｔ）が重畳された電流入力信号Ｉ（ｔ）が印加される。また、電流入力信号Ｉ（ｔ）印加により生じる電圧応答信号Ｖ（ｔ）が測定される。
【００４５】
図３には、二次電池に対して、電流入力信号Ｉ（ｔ）を印加した際の電圧Ｖの変化（電圧応答信号Ｖ（ｔ））を模式的に示す。図３（ａ）は、二次電池への電流入力信号Ｉ（ｔ）を示しており、図３（ｂ）は、二次電池の電圧応答信号Ｖ（ｔ）を示している。
【００４６】
図３（ａ）に示す例では、充電過程中の二次電池に対して、時間ｔｐにおいて、矩形波Ｉｐ（ｔ）が重畳されている。矩形波Ｉｐ（ｔ）は、図３（ａ）の丸枠内に拡大して示すように、マイナス側のステップ状波形を有する。従って、電流入力信号Ｉ（ｔ）は、ほとんどの時間において、電流値がＩｃで一定であるものの、時間ｔｐにおいて一時的に電流値が低下する。
【００４７】
このような電流入力信号Ｉ（ｔ）の印加によって、二次電池には、図３（ｂ）に示すような電圧応答信号Ｖ（ｔ）が生じる。すなわち、電圧応答信号Ｖ（ｔ）は、時間ｔｐにおいて、電圧値が一時的に不連続的に変化（低下）する。ステップＳ１２０では、このような電圧応答信号Ｖ（ｔ）が取得される。
【００４８】
ここで、本願において、「矩形波」とは、必ずしも波形が時間に対して正確な長方形で表される電流波形のみを意味するものではない。
【００４９】
図４〜図７には、本願において使用され得る「矩形波」の一例を模式的に示す。
【００５０】
図４において、矩形波１０は、時間ｔｐ１における時間軸（横軸）に対して電流値が実質的に垂直に低下する第１の部分（「初期部分」ともいう）１２と、時間軸に対して電流値が実質的に一定となる第２の部分（「中間部分」とも言う）１４と、時間ｔｐ２における時間軸に対して電流値が実質的に垂直に上昇する第３の部分（「最終部分」とも言う）１６とを有する。
【００５１】
しかしながら、高精度の波形発生装置（パルスジェネレータ等）を使用しても、実際の測定系において、このような「理想的な」矩形波１０を発生させることは難しい場合が多い。
【００５２】
次に、図５に示す矩形波２０は、時間ｔｐ１から時間ｔｐ３までの範囲の電流値が時間に対して実質的に低下する初期部分２２と、時間ｔｐ３から時間ｔｐ２までの範囲の時間軸に対して電流値が実質的に一定となる中間部分２４と、時間ｔｐ２における時間軸に対して電流値が実質的に垂直に上昇する最終部分２６とを有する。
【００５３】
また、図６に示す矩形波３０は、時間ｔｐ１における時間軸に対して電流値が実質的に垂直に低下する初期部分３２と、時間ｔｐ１から時間ｔｐ３までの範囲の時間軸に対して電流値が実質的に一定となる中間部分３４と、時間ｔｐ３から時間ｔｐ２までの範囲の電流値が時間に対して実質的に上昇する最終部分３６とを有する。
【００５４】
さらに、図７に示す矩形波４０は、時間ｔｐ１から時間ｔｐ３までの範囲の電流値が時間に対して実質的に低下する初期部分４２と、時間ｔｐ３から時間ｔｐ４までの範囲の時間軸に対して電流値が実質的に一定となる中間部分４４と、時間ｔｐ４から時間ｔｐ２までの範囲の電流値が時間に対して実質的に上昇する最終部分４６とを有する。
【００５５】
この他にも、矩形波Ｉｐ（ｔ）として、各種形態を選定し得ることは、当業者には明らかである。例えば、図４〜図７において、矩形波１０〜４０は、各部分が直線状に示されている。しかしながら、本願において、矩形波Ｉｐ（ｔ）は、初期部分および／または最終部分が時間に対して曲線状に変化しても良い。また、矩形波Ｉｐ（ｔ）は、中間部分の一部に、時間に対して変動する部分を含んでも良い。
【００５６】
従って、本願において「矩形波」とは、少なくとも、時間軸に対して実質的に水平な中間部分（例えば、中間部分１４、２４、３４、４４）を有する波形の総称を意味する。
【００５７】
なお、矩形波Ｉｐ（ｔ）において、パルス幅（すなわち、時間ｔｐ１〜ｔｐ２の期間の長さ）は、特に限られない。矩形波Ｉｐ（ｔ）のパルス幅は、例えば、１μｓ〜１ｋｓの範囲であっても良く、１０μｓ〜１００ｓの範囲であることが好ましい。
【００５８】
また、矩形波Ｉｐ（ｔ）の中間部分１４、２４、３４、４４における電流値は、０（ゼロ）であっても良い。このような中間部分１４、２４、３４、４４の電流値が０となるような矩形波Ｉｐ（ｔ）は、例えば、充電中の二次電池に通電される定電流Ｉｃを一時的に遮断することにより簡単に形成することができ、すなわち、波形発生装置を使用せずに発生させることができる。このため、このような矩形波Ｉｐ（ｔ）を使用する場合、第１の評価方法を実施するための装置を簡略化することが可能となる。
【００５９】
さらに、第１の評価方法において、定電流Ｉｃに矩形波Ｉｐ（ｔ）が重畳されるタイミング、すなわち時間ｔｐは、二次電池の充電期間中であれば、特に限られない。すなわち、この第１の評価方法では、二次電池のインピーダンス測定を、充電期間中の任意のタイミングで実施することができる。
【００６０】
さらに、第１の評価方法において、図３〜図７に示す例では、各矩形波Ｉｐ（ｔ）は、マイナス側（各図の下側）に変化する波形を有する。
【００６１】
矩形波Ｉｐ（ｔ）がこのような形状を有する場合、二次電池の充電中に、定電流Ｉｃに矩形波Ｉｐ（ｔ）を重畳し、この入力電流を二次電池に印加しても、これにより二次電池に所定の値を超える電気量が供給されることはない。従って、この第１の評価方法では、インピーダンス評価のために、二次電池に過充電が生じることを有意に抑制することができる。
【００６２】
なお、同様の観点から言えば、矩形波Ｉｐ（ｔ）は、いわゆる「ダブルパルス」波形を有しても良い。
【００６３】
図８には、そのようなダブルパルス波形を有する矩形波Ｉｐ（ｔ）の一例を模式的に示す。
【００６４】
図８に示すように、この矩形波６０は、マイナス側（下向き）およびプラス側（上向き）の、２つのステップ状部分を有する。より具体的には、矩形波６０は、時間ｔｐ１における時間軸に対して電流値が実質的に垂直に低下する初期部分６２と、時間ｔｐ１から時間ｔｐ３までの範囲の時間軸に対して電流値が実質的に一定となる第１の中間部分６４と、時間ｔｐ３における時間軸に対して電流値が実質的に垂直に上昇する上昇部分６６と、時間ｔｐ３から時間ｔｐ２までの範囲の時間軸に対して電流値が実質的に一定となる第２の中間部分６７と、時間ｔｐ２における時間軸に対して電流値が実質的に垂直に低下する最終部分６９とを有する。
【００６５】
なお、前述のように、中間部分６４、６７を除く各部分６２、６６、６９は、時間軸に対して非垂直であっても良く、各部分６２、６６、６９は、曲線状であっても良い。また、第１の中間部分６４と第２の中間部分６７の時間長さは、実質的に同一であっても、異なっていても良い。さらに、初期部分６２の変化量（減少量）は、最終部分６９の変化量（減少量）と同一であっても異なっていても良い。
【００６６】
矩形波Ｉｐ（ｔ）がこのようなダブルパルス波形を有する場合も、第１の中間部分６４の期間および電流値（高さ）、ならびに第２の中間部分６７の期間および電流値（高さ）等を適正に選定することにより、二次電池に過充電が生じることを有意に抑制することができる。例えば、図８において、時間ｔｐ３〜時間ｔｐ２の期間に二次電池に供給される電気量Ｑ２が、時間ｔｐ１〜時間ｔｐ３の期間における「不足する」電気量Ｑ１を超えないように矩形波６０を形状を定めることにより、二次電池に過充電が生じることを有意に抑制することができる。
【００６７】
また、矩形波Ｉｐ（ｔ）が図４〜図７に示すような「マイナス側」のパルス波形１０、２０、３０、４０を有する場合、厳密には、インピーダンスの測定時には、二次電池に対する定電流Ｉｃのみの印加による充電時に比べて、供給される電気量は、幾分小さくなる。例えば、図４に示した矩形波１０の例では、矩形波１０の面積に相当する電気量が、通常の充電過程に比べて「減少」することになる。
【００６８】
このため、二次電池の充電過程においてインピーダンスの測定を実施した場合、厳密には、インピーダンス測定を実施せずに、単に充電のみを実施した場合に比べて、ごく僅かではあるが、充電完了までにより長い時間が必要となる。
【００６９】
しかしながら、矩形波Ｉｐ（ｔ）としてダブルパルス波形６０を適用した場合、図８に示した各波形部分６２、６４、６６、６７、６９の値を適正に選定することにより、そのような「減少」する電気量を少なくすることができ、さらには「減少」する電気量を実質的にゼロにすることが可能となる。
【００７０】
このため、ダブルパルス波形６０では、インピーダンス測定によって、二次電池の充電時間が長くなることをより確実に回避することが可能となる。
【００７１】
以上の記載では、定電流Ｉｃが印加された二次電池に対して、単一のサイクルを有する（１サイクルで構成された）矩形波Ｉｐ（ｔ）を重畳、印加する場合を例に、本発明による第１の評価方法の特徴について説明した。
【００７２】
しかしながら、本発明による第１の評価方法において、インピーダンス測定の際に二次電池に印加される矩形波Ｉｐ（ｔ）は、１サイクルのみで構成されたものに限られない。
【００７３】
例えば、矩形波Ｉｐ（ｔ）は、複数の繰り返し波形で構成されても良い。すなわち、例えば、図４〜図８において示した矩形波１０、２０、３０、４０、または６０を波形の一サイクル（一単位）としたとき、矩形波Ｉｐ（ｔ）は、そのような一サイクルの複数回の繰り返しを有しても良い。
【００７４】
そのような複数回のサイクルの繰り返しを有する矩形波Ｉｐ（ｔ）を使用した場合、積分効果により、インピーダンス測定結果の精度が向上する。
【００７５】
なお、サイクルの繰り返し回数は、測定精度の観点からは多い方が好ましいが、あまりに繰り返し数が多くなると、インピーダンスの測定および評価に必要な時間が長くなる。このため、サイクルの繰り返し数は、２回〜１０回の範囲であることが好ましい。
【００７６】
（ステップＳ１３０）
次に、電流入力信号Ｉ（ｔ）が、以下の（２）式によってフーリエ変換され、第１の周波数領域データＰＩ（ω）が取得される。
【００７７】
【数１】
2014238948.tif 000003
【００７８】
ここで、jは、虚数単位であり、ωは、角周波数であり、ｔは、時間である。なお、ＰＩ（ω）は複素数であり、図９（ａ）には強度の絶対値のωに対するプロットが模式的に示されている。
【００７９】
また、電圧応答信号Ｖ（ｔ）が、以下の（３）式によってフーリエ変換され、第２の周波数領域データＰＶ（ω）が取得される。なお、ＰＶ（ω）は複素数であり、図９（ｂ）には強度の絶対値のωに対するプロットが模式的に示されている。
【００８０】
【数２】
2014238948.tif 000004
【００８１】
各信号のフーリエ変換には、市販のコンピュータ内に保管されたフーリエ変換ソフトウェアまたは高速フーリエ変換（ＦＦＴ）ソフトウェア等が使用されても良い。あるいは、各信号のフーリエ変換には、測定系に別途設けられた高速フーリエ変換（ＦＦＴ）装置などを使用しても良い。
【００８２】
図９には、フーリエ変換によって得られる第１の周波数領域データＰＩ（ω）（図９（ａ））、および第２の周波数領域データＰＶ（ω）（図９（ｂ））の一例を模式的に示す。
【００８３】
（ステップＳ１４０）
次に、ステップＳ１３０で取得された第１の周波数領域データＰＩ（ω）、および第２の周波数領域データＰＶ（ω）から、以下の（１）式を用いて、二次電池の所望の充電深度でのインピーダンススペクトルが計算される：
Ｚ（ω）＝ＰＶ（ω）／ＰＩ（ω）      （１）式
図１０には、（１）式によって得られる二次電池のインピーダンススペクトルの一例を模式的に示す。図１０において、横軸Ｚ'は、インピーダンス実数部を示し、縦軸Ｚ"は、インピーダンス虚数部を示している。従って、インピーダンスは、
Ｚ＝Ｚ'−ｊＺ"    （４）式
で示される複素数である。
【００８４】
以上の工程により、二次電池のインピーダンスを測定することができる。また、インピーダンスの評価後には、ステップＳ１１０における二次電池の充電過程が継続され、二次電池が充電される。
【００８５】
このような本発明による第１の評価方法では、二次電池に定電流Ｉｃを印加して二次電池を充電する過程中に、いわばｉｎ−ｓｉｔｕで、二次電池のインピーダンスを測定することができる。すなわち、本発明による第１の評価方法では、内部抵抗の測定を実施する際に、二次電池の充電を停止する必要がない。このため、第１の評価方法では、二次電池の充電または放電が完了するまでに必要な時間が延びてしまうという問題を有意に抑制することができる。
【００８６】
また、第１の評価方法では、インピーダンス測定は、充電過程中の二次電池に対して実施される。このため、従来のような平衡状態で得られる内部抵抗値とは異なり、測定されたインピーダンス結果は、二次電池の実際の使用状態における挙動を、より正確に反映したものとなる。従って、第１の評価方法では、二次電池の内部抵抗をより正確に評価することができる。
【００８７】
（追加のステップ）
以上のステップＳ１１０〜ステップＳ１４０により、二次電池のインピーダンスを評価することができる。ただし、本発明による第１の評価方法では、ステップＳ１４０の後に、さらに、
（１）式により得られたインピーダンススペクトルＺ（ω）から得られる二次電池の少なくとも一つの特性パラメータを、所定の値と比較するステップ（ステップＳ１５０）
を実施しても良い。
【００８８】
このようなステップＳ１５０を実施することにより、二次電池の劣化の状態を把握したり、二次電池の劣化の状態をモニタリングしたりすることができる。
【００８９】
例えば、ステップＳ１４０で測定されたインピーダンススペクトルＺ（ω）は、未使用の二次電池のインピーダンススペクトルと比較されても良い。
【００９０】
これにより、二次電池の劣化状態、さらには二次電池の交換のタイミングを把握することができる。例えば、ステップＳ１４０で測定されたインピーダンススペクトルＺ（ω）から得られる各種パラメータ（例えば、電解液抵抗、電荷移動抵抗、界面抵抗、アノード側抵抗、および／またはカソード側抵抗など）が、未使用の二次電池の値と大きくかけ離れた場合、この二次電池を備える装置を使用する使用者に対して警告を発することにより、二次電池の使用を中断させたり、二次電池の交換を促したりすることが可能となる。
【００９１】
あるいは、測定されたインピーダンススペクトルＺ（ω）から得られる特性データは、予め定められた閾値と比較されても良い。この場合も、例えば、ステップＳ１４０で測定されたインピーダンススペクトルＺ（ω）から得られる特性データが閾値を超えた（または下回った）際に、使用者に対して警告を発することで、二次電池の使用を中断させたり、二次電池の交換を促したりすることが可能となる。
【００９２】
（本発明の一実施例による二次電池のインピーダンスの別の評価方法について）
次に、図１１を参照して、本発明の一実施例による二次電池のインピーダンスの別の評価方法（以下、「（本発明による）第２の評価方法」と称する）について説明する。
【００９３】
図１１には、本発明による第２の評価方法の模式的なフロー図を示す。
【００９４】
図１１に示すように、本発明による第２の評価方法は、
（１）二次電池に定電流Ｉｃを通電して、前記二次電池を放電するステップ（Ｓ２１０）と、
（２）前記（１）中の任意のタイミングで、前記二次電池に電流入力信号Ｉ（ｔ）を印加し、電圧応答信号Ｖ（ｔ）を取得するステップであって、前記電流入力信号Ｉ（ｔ）は、前記定電流Ｉｃに矩形波Ｉｐ（ｔ）を重畳することにより得られるステップ（Ｓ２２０）と、
（３）電流入力信号Ｉ（ｔ）をフーリエ変換して、第１の周波数領域データＰＩ（ω）を取得し、前記電圧応答信号Ｖ（ｔ）をフーリエ変換して、第２の周波数領域データＰＶ（ω）を取得するステップ（Ｓ２３０）と、
（４）前記第１および第２の周波数領域データから、以下の式により、前記二次電池の所望の放電深度でのインピーダンススペクトルを計算するステップ（Ｓ２４０）と、
Ｚ（ω）＝ＰＶ（ω）／ＰＩ（ω）      （１）式
を有する。
【００９５】
以下、各ステップについて詳しく説明する。
【００９６】
（ステップＳ２１０）
まず、二次電池に対して定電流Ｉｃが通電され、二次電池が放電される。
【００９７】
図１２には、二次電池に対して、定電流Ｉｃが通電された際の電流Ｉおよび電圧Ｖの経時変化を模式的に示す。図１２（ａ）は、二次電池の電流変化挙動を示し、図１２（ｂ）は、二次電池の電圧変化挙動を示している。
【００９８】
なお、この図１２の例では、定電流Ｉｃは負の値であると仮定している（Ｉｃ＜０）。
【００９９】
この場合、定電流Ｉｃの通電により二次電池が放電されるため、図２（ｂ）に示すように、二次電池の電圧Ｖは、時間とともに徐々に低下する。
【０１００】
定電流Ｉｃの値は、特に限られず、使用される二次電池の種類、容量、および二次電池に接続された負荷等に応じて、適宜決定される。例えば、携帯電話等に使用される小型のリチウムイオン二次電池の場合、定電流Ｉｃの絶対値は、数ｍＡ〜数Ａの範囲であっても良い。
【０１０１】
また、このステップＳ２１０の際に、図１２（ｂ）に示すような二次電池の放電曲線（電圧−時間曲線）が取得されても良い。なお、二次電池の放電は、過放電とならない時間で完了される。
【０１０２】
（ステップＳ２２０）
前述のステップＳ２１０の期間中の任意のタイミングにおいて、定電流Ｉｃに矩形波Ｉｐ（ｔ）が重畳される。すなわち、二次電池に対して、定電流Ｉｃに矩形波Ｉｐ（ｔ）が重畳された電流入力信号Ｉ（ｔ）が印加される。また、電流入力信号Ｉ（ｔ）印加により生じる電圧応答信号Ｖ（ｔ）が測定される。
【０１０３】
図１３には、二次電池に対して、電流入力信号Ｉ（ｔ）を印加した際の電圧Ｖの変化（電圧応答信号Ｖ（ｔ））を模式的に示す。図１３（ａ）は、二次電池への電流入力信号Ｉ（ｔ）を示しており、図１３（ｂ）は、二次電池の電圧応答信号Ｖ（ｔ）を示している。
【０１０４】
図１３（ａ）に示す例では、放電過程中の二次電池に対して、時間ｔｐにおいて、矩形波Ｉｐ（ｔ）が重畳されている。矩形波Ｉｐ（ｔ）は、図１３（ａ）の丸枠内に拡大して示すように、プラス側のステップ状波形を有する。従って、電流入力信号Ｉ（ｔ）は、ほとんどの時間において、電流値がＩｃで一定であるものの、時間ｔｐにおいて一時的に電流値が増加する（電流の絶対値が減少する）。
【０１０５】
このような電流入力信号Ｉ（ｔ）の印加によって、二次電池には、図１３（ｂ）に示すような電圧応答信号Ｖ（ｔ）が生じる。すなわち、電圧応答信号Ｖ（ｔ）は、時間ｔｐにおいて、電圧値が一時的に不連続的に変化（低下）する。ステップＳ２２０では、このような電圧応答信号Ｖ（ｔ）が取得される。
【０１０６】
第２の評価方法においても、定電流Ｉｃに矩形波Ｉｐ（ｔ）が重畳されるタイミング、すなわち時間ｔｐは、二次電池の放電期間中であれば、特に限られず、このため二次電池のインピーダンス測定は、二次電池の放電期間中の任意のタイミングで実施することができる。
【０１０７】
また、第２の評価方法において、矩形波Ｉｐ（ｔ）は、図１３（ａ）に示すように、プラス側（図１３の上側）に変化する波形を有しても良い。
【０１０８】
この場合、二次電池の放電中に、定電流Ｉｃに矩形波Ｉｐ（ｔ）を重畳しても、これにより二次電池から所定の値を超える電気量が放出されることはない。従って、この第２の評価方法では、インピーダンス評価のために、二次電池に過放電が生じることを有意に抑制することができる。
【０１０９】
このような第２の評価方法においても、例えば前述の図４〜図８等に示した各種波形を有する矩形波を適用することができる。ただし、この第２の評価方法の場合、矩形波Ｉｐ（ｔ）の波形は、図４〜図８に示した各波形が上下に反転された形態であることが好ましい。
【０１１０】
なお、前述のように、インピーダンス測定の際に二次電池に印加される矩形波Ｉｐ（ｔ）は、１サイクルのみで構成されたものに限られず、矩形波Ｉｐ（ｔ）は、複数の繰り返し波形で構成されても良い。
【０１１１】
（ステップＳ２３０）
次に、電流入力信号Ｉ（ｔ）が、以下の（２）式によってフーリエ変換され、第１の周波数領域データＰＩ（ω）が取得される。
【０１１２】
【数３】
2014238948.tif 000005
【０１１３】
ここで、ωは、角周波数であり、ｔは、時間である。
【０１１４】
また、電圧応答信号Ｖ（ｔ）が、以下の（３）式によってフーリエ変換され、第２の周波数領域データＰＶ（ω）が取得される。
【０１１５】
【数４】
2014238948.tif 000006
【０１１６】
各信号のフーリエ変換には、市販のコンピュータ内に保管されたフーリエ変換ソフトウェアまたは高速フーリエ変換（ＦＦＴ）ソフトウェア等が使用されても良い。あるいは、各信号のフーリエ変換には、測定系に別途設けられた高速フーリエ変換（ＦＦＴ）装置などを使用しても良い。
【０１１７】
図１４には、フーリエ変換によって得られた第１の周波数領域データＰＩ（ω）（図１４（ａ））、および第２の周波数領域データＰＶ（ω）（図１４（ｂ））の一例を模式的に示す。なお、ＰＩ（ω）とＰＶ（ω）は共に複素数であり、図１４にはそれぞれの強度の絶対値のωに対するプロットが模式的に示されている。
【０１１８】
（ステップＳ２４０）
次に、ステップＳ２３０で取得された第１の周波数領域データＰＩ（ω）、および第２の周波数領域データＰＶ（ω）から、以下の（１）式を用いて、所望の放電深度における二次電池のインピーダンススペクトルが計算される：
Ｚ（ω）＝ＰＶ（ω）／ＰＩ（ω）      （１）式
図１５には、（１）式によって得られる二次電池のインピーダンススペクトルの一例を模式的に示す。図１５において、横軸Ｚ'は、インピーダンス実数部を示し、縦軸Ｚ"は、インピーダンス虚数部を示している。従って、インピーダンスは、
Ｚ＝Ｚ'−ｊＺ"    （４）式
で示される複素数である。
【０１１９】
以上の工程により、放電過程中の二次電池のインピーダンスを測定することができる。また、インピーダンスの評価後には、ステップＳ２１０における二次電池の放電過程が継続され、二次電池が最終段階まで放電される。
【０１２０】
第２の評価方法では、二次電池に定電流Ｉｃを通電して二次電池を放電する過程中に、いわばｉｎ−ｓｉｔｕで、二次電池のインピーダンスを測定することができる。すなわち、第２の評価方法では、内部抵抗の測定を実施する際に、二次電池の放電を停止する必要がない。このため、第２の評価方法では、二次電池の放電曲線を取得するために必要な時間が長くなってしまうという問題を有意に抑制することができる。
【０１２１】
また、第２の評価方法では、インピーダンス測定は、放電過程中の二次電池に対して実施される。このため、従来のような平衡状態で得られる内部抵抗値とは異なり、測定されたインピーダンス結果は、二次電池の実際の使用状態における挙動を、より正確に反映したものとなる。従って、第２の評価方法においても、二次電池の内部抵抗をより正確に評価することが可能となる。
【０１２２】
（追加のステップ）
なお、図１１には示されていないが、第２の評価方法においても、ステップＳ２４０の後に、さらに、
（１）式により得られたインピーダンススペクトルＺ（ω）から得られる二次電池の少なくとも一つの特性パラメータを、所定の値と比較するステップ
を実施しても良い。
【実施例】
【０１２３】
以下、本発明の実施例について説明する。
【０１２４】
二次電池を模擬した等価回路を用いて、前述の第１の評価方法と同様の手順でインピーダンス計算を実施した。
【０１２５】
インピーダンスの評価は、数値計算ソフトウェア（ＭＡＴＬＡＢ：ＭａｔｈＷｏｒｋｓ社製）を用いて実施した。
【０１２６】
図１６には、計算に使用した二次電池の等価回路を示す。図１６に示すように、この等価回路８０は、電解液抵抗８２と、電荷移動抵抗８４と、電気二重層容量８６とで構成されている。電解液抵抗８２は０．１Ωとし、電荷移動抵抗８４は１Ωとし、電気二重層容量８６は１×１０−３Ｆと仮定した。
【０１２７】
フーリエ変換の条件は、以下の通りである：
・窓関数    なし
・サンプリング時間    １μｓ
図１７には、二次電池に印加される電流入力信号Ｉ（ｔ）（＞０）の形状を示す。この図に示すように、電流入力信号Ｉ（ｔ）は、２Ａの定電流Ｉｃに、マイナス側（−０．５Ａ）の矩形波（パルス幅＝５秒）を重畳して構成されたものである。
【０１２８】
図１８には、二次電池にこのような電流入力信号Ｉ（ｔ）を入力した際に得られる出力計算結果、すなわち電圧応答信号Ｖ（ｔ）の波形の拡大図を示す。
【０１２９】
次に、電流入力信号Ｉ（ｔ）をフーリエ変換し、第１の周波数データＰＩ（ｆ）を取得した。同様に、電圧応答信号Ｖ（ｔ）をフーリエ変換し、第２の周波数データＰＶ（ｆ）を取得した。第１の周波数データＰＩ（ｆ）および第２の周波数データＰＶ（ｆ）は複素数であり、それぞれの強度の絶対値を、図１９（ａ）および図１９（ｂ）に示した。
【０１３０】
さらに、第１の周波数データＰＩ（ｆ）および第２の周波数データＰＶ（ｆ）を用いて、（１）'式から、二次電池のインピーダンススペクトルを計算した。
Ｚ（ｆ）＝ＰＶ（ｆ）／ＰＩ（ｆ）      （１）'式
得られた結果（ナイキストプロット）を、図２０に示す。
【０１３１】
この図に示すように、インピーダンススペクトルは、略半円状の軌跡を示し、円の直径から、二次電池の電荷移動抵抗は、１Ωであることが計算された。
【産業上の利用可能性】
【０１３２】
本発明は、例えば、通信機器、ＡＶ機器、情報機器、および自動車などに使用されている二次電池等に適用することができる。
【符号の説明】
【０１３３】
１０      矩形波
１２      第１の部分（初期部分）
１４      第２の部分（中間部分）
１６      第３の部分（最終部分）
２０      矩形波
２２      初期部分
２４      中間部分
２６      最終部分
３０      矩形波
３２      初期部分
３４      中間部分
３６      最終部分
４０      矩形波
４２      初期部分
４４      中間部分
４６      最終部分
６０      矩形波
６２      初期部分
６４      第１の中間部分
６６      上昇部分
６７      第２の中間部分
６９      最終部分
８０      等価回路
８２      電解液抵抗
８４      電荷移動抵抗
８６      電気二重層容量
【図１】
2014238948.tif 000007
【図２】
2014238948.tif 000008
【図３】
2014238948.tif 000009
【図４】
2014238948.tif 000010
【図５】
2014238948.tif 000011
【図６】
2014238948.tif 000012
【図７】
2014238948.tif 000013
【図８】
2014238948.tif 000014
【図９】
2014238948.tif 000015
【図１０】
2014238948.tif 000016
【図１１】
2014238948.tif 000017
【図１２】
2014238948.tif 000018
【図１３】
2014238948.tif 000019
【図１４】
2014238948.tif 000020
【図１５】
2014238948.tif 000021
【図１６】
2014238948.tif 000022
【図１７】
2014238948.tif 000023
【図１８】
2014238948.tif 000024
【図１９】
2014238948.tif 000025
【図２０】
2014238948.tif 000026