ホーム > 特許ランキング > 本田技研工業株式会社
知財求人 - 知財ポータルサイト「IP Force」
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2023-11-14
(45)【発行日】2023-11-22
(54)【発明の名称】測定装置、測定方法、及びプログラム
(51)【国際特許分類】
G01R 27/02 20060101AFI20231115BHJP
G01R 31/389 20190101ALI20231115BHJP
G01R 31/367 20190101ALI20231115BHJP
H01M 10/48 20060101ALI20231115BHJP
【ＦＩ】
G01R27/02 A
G01R31/389
G01R31/367
H01M10/48 P
【請求項の数】
3
(21)【出願番号】P 2022046656
(22)【出願日】2022-03-23
(65)【公開番号】P2023140688
(43)【公開日】2023-10-05
【審査請求日】2022-11-29
(73)【特許権者】
【識別番号】000005326
【氏名又は名称】本田技研工業株式会社
(74)【代理人】
【識別番号】100165179
【弁理士】
【氏名又は名称】田▲崎▼ 聡
(74)【代理人】
【識別番号】100126664
【弁理士】
【氏名又は名称】鈴木 慎吾
(74)【代理人】
【識別番号】100154852
【弁理士】
【氏名又は名称】酒井 太一
(74)【代理人】
【識別番号】100194087
【弁理士】
【氏名又は名称】渡辺 伸一
(72)【発明者】
【氏名】冨永 由騎
(72)【発明者】
【氏名】西本 有里佳
【審査官】永井 皓喜
(56)【参考文献】
【文献】国際公開第２０２０／２６１７９９（ＷＯ，Ａ１）
【文献】登録実用新案第３２０７９９８（ＪＰ，Ｕ）
【文献】特開２０１５－１９７３６３（ＪＰ，Ａ）
【文献】特開２０１２－７９５８２（ＪＰ，Ａ）
【文献】特表２００８－５２２３５０（ＪＰ，Ａ）
【文献】中国特許出願公開第１０４９１４３１２（ＣＮ，Ａ）
(58)【調査した分野】(Int.Cl.，ＤＢ名)
Ｇ０１Ｒ ２７／０２
Ｇ０１Ｒ ３１／３８９
Ｇ０１Ｒ ３１／３６７
Ｈ０１Ｍ １０／４８
Ｈ０２Ｊ ７／００
(57)【特許請求の範囲】
【請求項1】
二次電池の交流インピーダンスの測定装置であって、インピーダンス演算器と、前記インピーダンス演算器と前記インピーダンス演算器に接続された機器と信号を入出力するための信号入出力インターフェイスと、を備え、前記インピーダンス演算器は次の（Ｓ１）～（Ｓ５）の処理を実行する測定装置。（Ｓ１）前記インピーダンス演算器は、横軸に周波数ｆ（Ｈｚ）、縦軸に位相θ（度）をプロットしてボード線図を作成する。ボード線図の作成時に、前記インピーダンス演算器は、インダクタンス成分と拡散成分のフィッティングを行い、測定結果のボード線図から除去する。
（Ｓ２）前記インピーダンス演算器は、作成したボード線図上で頂点を抽出して、Ｒ／Ｃ並列回路数として＋１をカウントし、前記頂点を中心にＲ／Ｃ並列回路で表すことができる山の形状をフィッティングする。ただし、測定ノイズと頂点とを見分けるため、棄却するθ値（θ０）が設定されており、前記インピーダンス演算器は、位相（θ）の値が予め定めたθ０以上である点から頂点を抽出する。
（Ｓ３）前記インピーダンス演算器は、フィッティングしたＲ／Ｃ並列回路で表現できるボード線図成分を、測定結果のボード線図から減ずる。
（Ｓ４）前記インピーダンス演算器が、ボード線図上で頂点を抽出しなくなればＲ／Ｃ並列回路数の同定を終了するが、頂点を抽出すれば、前記（Ｓ２）及び前記（Ｓ３）を繰り返す。ただし、頂点の位相（θ）の大きさは、予め設定した位相（θ）の大きさが予め定めたθ１以上である。
（Ｓ５）前記インピーダンス演算器が、Ｒ／Ｃ並列回路数の総カウント数を報告する。
【請求項2】
二次電池の交流インピーダンスの測定方法であって、プロセッサが、横軸に周波数ｆ（Ｈｚ）（対数軸）、縦軸に位相θ（度）をプロットしてボード線図を作成するステップＳ１と、
プロセッサが、作成したボード線図上で頂点を抽出して、Ｒ／Ｃ並列回路数として＋１をカウントし、前記頂点を中心にＲ／Ｃ並列回路で表すことができる山の形状をフィッティングするステップ２と、
プロセッサが、フィッティングしたＲ／Ｃ並列回路で表現できるボード線図成分を、作成したボード線図から減ずるステップＳ３と、
プロセッサが、ボード線図上で頂点を抽出しなければステップＳ５に進み、さらに頂点を抽出すれば、ステップＳ２及びステップＳ３を繰り返すように処理を選択するステップＳ４と、
プロセッサが、Ｒ／Ｃ並列回路数の総カウント数を報告して、Ｒ／Ｃ並列回路数の同定を終了するステップＳ５と
を含む、測定方法。
ただし、ステップＳ１において、ボード線図を作成する際に、インダクタンス成分と拡散成分のフィッティングを行い、測定結果のボード線図から除去する。
また、Ｓ２において、測定ノイズと頂点とを見分けるため、頂点の位相（θ）の大きさの下限値（θ０）を予め設定しておき、プロセッサが、θ値がθ０以上である頂点候補から頂点を抽出する。
さらに、ステップＳ４において、頂点の位相（θ）の大きさの下限値（θ１）を予め設定しておき、プロセッサが、θ値がθ１以上である頂点候補を頂点として抽出する。
【請求項3】
二次電池の交流インピーダンスの測定のためのプログラムであって、インピーダンス演算器と前記インピーダンス演算器と前記インピーダンス演算器に接続された機器との間で信号を入出力するための信号入出力インターフェイスと、を備えるインピーダンス測定装置において、前記インピーダンス演算器に搭載されたプロセッサにお以下の（Ｓ１）～（Ｓ５）の処理を命じるプログラム。（Ｓ１）横軸に周波数ｆ（Ｈｚ）（対数軸）、縦軸に位相θ（度）をプロットしてボード線図を作成する。
（Ｓ２）作成したボード線図上で頂点を抽出し、Ｒ／Ｃ並列回路数として＋１をカウントし、前記頂点を中心にＲ／Ｃ並列回路で表すことができる山の形状をフィッティングする。
（Ｓ３）フィッティングしたＲ／Ｃ並列回路で表現できるボード線図成分を、測定結果のボード線図から減じる。
（Ｓ４）ボード線図上で頂点を抽出しなくなればＲ／Ｃ並列回路の総カウント数を報告してＲ／Ｃ並列回路数の同定を終了し、頂点をさらに抽出すれば、（Ｓ２）及び（Ｓ３）の処理を繰り返す。
（Ｓ５）Ｒ／Ｃ並列回路数の総カウント数を報告する。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、測定装置、測定方法、及びプログラムに関する。
【背景技術】
【0002】
近年、リチウムイオン二次電池等の二次電池が搭載された電動車両の普及が進んでいる。二次電池は、充放電の繰り返しとともに、又は時間の経過とともに劣化する。二次電池の劣化の進行度合いを診断する手法として、交流インピーダンス測定法が公知である。
【0003】
交流インピーダンス測定法では、電極系に印加した交流電圧（電流）に対する応答電流（電圧）を測定し、その比がインピーダンスとなる。その表記法としては、ボード線図とナイキスト図がよく用いられる。ボード線図は、印加した交流電圧の周波数ｆ（Ｈｚ）の対数を横軸に、インピーダンスの絶対値｜Ｚ｜（Ω）の対数と位相差θ（ｄｅｇ）を、それぞれ縦軸にプロットした図である。｜Ｚ｜は印加した交流電圧と応答電流の振幅の比、θは電圧と電流の位相の差である。ナイキスト図では、横軸にインピーダンスの実数部を、縦軸に虚数部をとりインピーダンスを複素平面上に表記する。実数部Ｒｅ［Ｚ］（抵抗）は、式：Ｒｅ［Ｚ］＝｜Ｚ｜ｃｏｓθにより、虚数部Ｉｍ［Ｚ］（リアクタンス）は、式：Ｉｍ［Ｚ］＝｜Ｚ｜ｓｉｎθにより、それぞれ計算できる。
【0004】
特許文献１には、複数の二次電池を含む電池モジュールの特性に関する情報を処理する電池情報処理システムであって、前記電池モジュールの交流インピーダンスの測定結果を解析する解析装置と、前記交流インピーダンスを表現する等価回路モデルに含まれる複数の回路定数と、前記特性との間の相関関係が記憶された記憶装置とを備え、前記解析装置は、前記交流インピーダンスの実数成分に関するボード線図である第１の周波数特性図上と、前記交流インピーダンスの虚数成分に関するボード線図である第２の周波数特性図上とに、前記交流インピーダンスの測定結果をプロットし、前記第１の周波数特性図上のプロット結果に対するフィッティング処理により第１の多項式曲線を取得するとともに、前記第２の周波数特性図上のプロット結果に対するフィッティング処理により第２の多項式曲線を取得し、前記第１及び第２の多項式曲線を、前記交流インピーダンスの実数成分と虚数成分とに関するナイキスト図上のインピーダンス曲線に変換し、前記インピーダンス曲線から前記複数の回路定数を抽出し、前記相関関係を参照することによって、抽出された前記複数の回路定数から前記特性を評価する、電池情報処理システムが記載されている。
【0005】
特許文献２には、複数個の電池セルが直列に接続され、実負荷を高周波域を含む負荷変動を生じる状態で駆動する電池モジュールをリアルタイムで測定監視する電池監視装置に用いられるインピーダンス測定装置であって、前記各セルの電圧波形データ及び電流波形データを離散フーリエ変換し、電圧波形データの離散フーリエ変換結果を電流波形データの離散フーリエ変換結果で除算することによりインピーダンスを演算するＤＦＴ演算部と、このＤＦＴ演算部で演算されたインピーダンスデータの特徴に基づいて最適な等価回路モデルを選択する回路モデル選択部と、この回路モデル選択部で選択された等価回路モデルにおいて定数フィッティングを行う回路定数推定演算部、とで構成されている電池インピーダンス測定装置が記載され、さらに前記電池インピーダンス測定装置において、前記回路モデル選択部は、等価回路モデルの具体的な構成について、１）ワールブルグ素子の有無、２）ＬＲ並列回路の有無、３）ＲＣ並列回路の段数の手順で順次決定することが記載されている。
【先行技術文献】
【特許文献】
【0006】
【文献】特開２０１９－１９１０３０号公報
【文献】特開２０１３－０２９４１２号公報
【発明の概要】
【発明が解決しようとする課題】
【0007】
交流インピーダンス測定では、交流の入力信号（電圧又は電流）に対する出力信号（電圧又は電流）の比｜Ｚ｜と位相遅れθを測定する。入力信号の周波数ｆを変えることによって異なる律速過程の｜Ｚ｜及びθが得られる。すなわち、（１）微弱な交流信号を入力する、（２）応答を時間場から周波数場へ高速フーリエ変換する、（３）出力／入力インピーダンス比｜Ｚ｜及び位相遅れθを算出する（ボード線図）、（４）インピーダンス実部、虚部を複素平面上に表示する（ナイキスト図）、という手順で測定が行われる。インピーダンススペクトルは、横軸に実部（Ｚ’）、縦軸に虚部（Ｚ”）をプロットしたナイキスト図で表現される。二次電池の等価回路として、抵抗ＲとキャパシタＣを並列接続したＲ／Ｃ並列回路が用いられる。二次電池の交流インピーダンス測定においては、測定対象の二次電池の等価回路を決める必要があるが、理想的には、１つのＲ／Ｃ並列回路が、ナイキスト図上で独立した半円弧として現れる。しかし、実際の二次電池では複数の円弧が重なり合って現れるため、正しく分離することが難しい。特に、他と比較して小さい円弧成分は目視確認が難しく、詳細な抵抗成分を見落としていた。すなわち、Ｒ／Ｃ並列回路で表現される円弧成分が何個あるかを見極めるのが困難であった。
【0008】
本発明は、このような事情を考慮してなされたものであり、二次電池の交流インピーダンス測定において、等価回路におけるＲ／Ｃ並列回路数をより精確にカウントできる、測定装置、測定方法、及びプログラムを提供することを目的の一つとする。
【課題を解決するための手段】
【0009】
この発明に係る測定装置、測定方法、及びプログラムは、以下の構成を採用した。
（１）：この発明の一態様に係る測定装置は、インピーダンス演算器と、前記インピーダンス演算器と信号を入出力するための信号入出力インターフェイスと、を備え、前記インピーダンス演算器は次の（Ｓ１）～（Ｓ５）の処理を実行する測定装置である。
（Ｓ１）前記インピーダンス演算器は、横軸に周波数ｆ（Ｈｚ）、縦軸に位相θ（度）をプロットしてボード線図を作成する。
（Ｓ２）前記インピーダンス演算器は、作成したボード線図上で頂点を抽出して、Ｒ／Ｃ並列回路数として＋１をカウントし、前記頂点を中心にＲ／Ｃ並列回路で表すことができる山の形状をフィッティングする。
（Ｓ３）前記インピーダンス演算器は、フィッティングしたＲ／Ｃ並列回路で表現できるボード線図成分を、測定結果のボード線図から減ずる。
（Ｓ４）前記インピーダンス演算器が、ボード線図上で頂点を抽出しなくなればＲ／Ｃ並列回路数の同定を終了するが、頂点を抽出すれば、前記（Ｓ２）及び前記（Ｓ３）を繰り返す。ただし、頂点の位相（θ）の大きさは、予め設定した位相（θ）の大きさが予め定めたθ１以上である。
（Ｓ５）前記インピーダンス演算器が、Ｒ／Ｃ並列回路数の総カウント数を報告する。
ただし、（Ｓ１）において、前記インピーダンス演算器は、ボード線図を作成する際に、インダクタンス成分と拡散成分のフィッティングを行い、測定結果のボード線図から除去する。
また、（Ｓ２）において、測定ノイズと頂点とを見分けるため、頂点の位相（θ）の大きさの下限値（θ０）を予め設定しておき、前記インピーダンス演算器は、θ値がθ０以上である頂点候補から頂点を抽出する。
さらに、ステップＳ４において、頂点の位相（θ）の大きさの下限値（θ１）を予め設定しておき、前記インピーダンス演算器は、θ値がθ１以上である頂点候補を頂点として抽出する。
（１）：この発明の一態様に係る測定装置は、インピーダンス演算器と、前記インピーダンス演算器と信号を入出力するための信号入出力インターフェイスと、を備え、前記インピーダンス演算器は次の（Ｓ１）～（Ｓ５）の処理を実行する測定装置である。
（Ｓ１）前記インピーダンス演算器は、横軸に周波数ｆ（Ｈｚ）、縦軸に位相θ（度）をプロットしてボード線図を作成する。
（Ｓ２）前記インピーダンス演算器は、作成したボード線図上で頂点を抽出して、Ｒ／Ｃ並列回路数として＋１をカウントし、前記頂点を中心にＲ／Ｃ並列回路で表すことができる山の形状をフィッティングする。
（Ｓ３）前記インピーダンス演算器は、フィッティングしたＲ／Ｃ並列回路で表現できるボード線図成分を、測定結果のボード線図から減ずる。
（Ｓ４）前記インピーダンス演算器が、ボード線図上で頂点を抽出しなくなればＲ／Ｃ並列回路数の同定を終了するが、頂点を抽出すれば、前記（Ｓ２）及び前記（Ｓ３）を繰り返す。ただし、頂点の位相（θ）の大きさは、予め設定した位相（θ）の大きさが予め定めたθ１以上である。
（Ｓ５）前記インピーダンス演算器が、Ｒ／Ｃ並列回路数の総カウント数を報告する。
ただし、（Ｓ１）において、前記インピーダンス演算器は、ボード線図を作成する際に、インダクタンス成分と拡散成分のフィッティングを行い、測定結果のボード線図から除去する。
また、（Ｓ２）において、測定ノイズと頂点とを見分けるため、頂点の位相（θ）の大きさの下限値（θ０）を予め設定しておき、前記インピーダンス演算器は、θ値がθ０以上である頂点候補から頂点を抽出する。
さらに、ステップＳ４において、頂点の位相（θ）の大きさの下限値（θ１）を予め設定しておき、前記インピーダンス演算器は、θ値がθ１以上である頂点候補を頂点として抽出する。
【0010】
（２）：この発明の一態様に係る測定方法は、
プロセッサが、横軸に周波数ｆ（Ｈｚ）（対数軸）、縦軸に位相θ（度）をプロットしてボード線図を作成するステップＳ１と、
プロセッサが、作成したボード線図上で頂点を抽出して、Ｒ／Ｃ並列回路数として＋１をカウントし、前記頂点を中心にＲ／Ｃ並列回路で表すことができる山の形状をフィッティングするステップ２と、
プロセッサが、フィッティングしたＲ／Ｃ並列回路で表現できるボード線図成分を、作成したボード線図から減ずるステップＳ３と、
プロセッサが、ボード線図上で頂点を抽出しなければステップＳ５に進み、さらに頂点を抽出すれば、ステップＳ２及びステップＳ３を繰り返すように処理を選択するステップＳ４と、
プロセッサが、Ｒ／Ｃ並列回路数の総カウント数を報告して、Ｒ／Ｃ並列回路数の同定を終了するステップＳ５と
を含む。
ただし、ステップＳ１において、ボード線図を作成する際に、インダクタンス成分と拡散成分のフィッティングを行い、測定結果のボード線図から除去する。
また、Ｓ２において、測定ノイズと頂点とを見分けるため、頂点の位相（θ）の大きさの下限値（θ０）を予め設定しておき、プロセッサが、θ値がθ０以上である頂点候補から頂点を抽出する。
さらに、ステップＳ４において、頂点の位相（θ）の大きさの下限値（θ１）を予め設定しておき、プロセッサが、θ値がθ１以上である頂点候補を頂点として抽出する。
プロセッサが、横軸に周波数ｆ（Ｈｚ）（対数軸）、縦軸に位相θ（度）をプロットしてボード線図を作成するステップＳ１と、
プロセッサが、作成したボード線図上で頂点を抽出して、Ｒ／Ｃ並列回路数として＋１をカウントし、前記頂点を中心にＲ／Ｃ並列回路で表すことができる山の形状をフィッティングするステップ２と、
プロセッサが、フィッティングしたＲ／Ｃ並列回路で表現できるボード線図成分を、作成したボード線図から減ずるステップＳ３と、
プロセッサが、ボード線図上で頂点を抽出しなければステップＳ５に進み、さらに頂点を抽出すれば、ステップＳ２及びステップＳ３を繰り返すように処理を選択するステップＳ４と、
プロセッサが、Ｒ／Ｃ並列回路数の総カウント数を報告して、Ｒ／Ｃ並列回路数の同定を終了するステップＳ５と
を含む。
ただし、ステップＳ１において、ボード線図を作成する際に、インダクタンス成分と拡散成分のフィッティングを行い、測定結果のボード線図から除去する。
また、Ｓ２において、測定ノイズと頂点とを見分けるため、頂点の位相（θ）の大きさの下限値（θ０）を予め設定しておき、プロセッサが、θ値がθ０以上である頂点候補から頂点を抽出する。
さらに、ステップＳ４において、頂点の位相（θ）の大きさの下限値（θ１）を予め設定しておき、プロセッサが、θ値がθ１以上である頂点候補を頂点として抽出する。
【0011】
（３）：この発明の一態様に係るプログラムは、
インピーダンス測定装置が備えるインピーダンス演算器に搭載されたプロセッサが、
横軸に周波数ｆ（Ｈｚ）（対数軸）、縦軸に位相θ（度）をプロットしてボード線図を作成して、作成したボード線図上で頂点を抽出し、Ｒ／Ｃ並列回路数として＋１をカウントし、前記頂点を中心にＲ／Ｃ並列回路で表すことができる山の形状をフィッティングし、
フィッティングしたＲ／Ｃ並列回路で表現できるボード線図成分を、測定結果のボード線図から減じ、
ボード線図上で頂点を抽出しなくなればＲ／Ｃ並列回路の総カウント数を報告してＲ／Ｃ並列回路数の同定を終了し、頂点をさらに抽出すれば、Ｒ／Ｃ並列回路数として＋１をカウントし、前記頂点を中心にＲ／Ｃ並列回路で表すことができる山の形状をフィッティングし、フィッティングしたＲ／Ｃ並列回路で表現できるボード線図成分を、測定結果のボード線図から減じることを繰り返す、
プログラムである。
インピーダンス測定装置が備えるインピーダンス演算器に搭載されたプロセッサが、
横軸に周波数ｆ（Ｈｚ）（対数軸）、縦軸に位相θ（度）をプロットしてボード線図を作成して、作成したボード線図上で頂点を抽出し、Ｒ／Ｃ並列回路数として＋１をカウントし、前記頂点を中心にＲ／Ｃ並列回路で表すことができる山の形状をフィッティングし、
フィッティングしたＲ／Ｃ並列回路で表現できるボード線図成分を、測定結果のボード線図から減じ、
ボード線図上で頂点を抽出しなくなればＲ／Ｃ並列回路の総カウント数を報告してＲ／Ｃ並列回路数の同定を終了し、頂点をさらに抽出すれば、Ｒ／Ｃ並列回路数として＋１をカウントし、前記頂点を中心にＲ／Ｃ並列回路で表すことができる山の形状をフィッティングし、フィッティングしたＲ／Ｃ並列回路で表現できるボード線図成分を、測定結果のボード線図から減じることを繰り返す、
プログラムである。
【発明の効果】
【0012】
（１）～（３）の態様によれば、合成インピーダンスを表すボード線図では区別できなかった円弧成分を区別することが可能となり、等価回路におけるＲ／Ｃ並列回路数をより精確にカウントできるようになる。
【図面の簡単な説明】
【0013】
【図1】測定装置の使用形態を説明するブロック図である。
【図2】本実施形態のインピーダンス測定装置によるＲ／Ｃ等価回路数のカウントを説明するためのボード線図である。
【図3】本実施形態のインピーダンス測定装置によるＲ／Ｃ等価回路数のカウントを説明するためのボード線図である。
【図4】本実施形態のインピーダンス測定装置によるＲ／Ｃ等価回路数のカウントを説明するためのボード線図である。
【図5】本実施形態のインピーダンス測定装置によるＲ／Ｃ等価回路数のカウントを説明するためのボード線図である。
【図6】本実施形態のインピーダンス測定装置によるＲ／Ｃ等価回路数のカウントを説明するためのボード線図である。
【図7】本実施形態のインピーダンス測定装置によるＲ／Ｃ等価回路数のカウントを説明するためのボード線図である。
【図8】本実施形態のインピーダンス測定装置によるＲ／Ｃ等価回路数のカウントを説明するためのボード線図である。
【図9】従来のインピーダンス測定装置によるＲ／Ｃ等価回路数のカウントを説明するためのボード線図である。
【図10】本実施形態のインピーダンス測定方法を説明するフローチャートである。
【発明を実施するための形態】
【0014】
以下、図面を参照し、本発明の測定装置、測定方法、及びプログラムの実施形態について説明する。
本発明において、二次電池は特に限定されないが、リチウムイオン二次電池が好ましい態様である。
本発明において、二次電池は特に限定されないが、リチウムイオン二次電池が好ましい態様である。
【0015】
［測定装置］
本実施形態の測定装置１００は、後述するインピーダンス演算器８０及びインピーダンス演算器８０と信号を入出力するための信号入出力インターフェイス（図示せず）を備える。インピーダンス演算器８０は、また、交流電圧モニタ６０及び交流電流モニタ７０からの信号を入力するためのインターフェイスを備える。測定装置１００はインピーダンス演算器８０からなるものでもよいが、交流電圧モニタ６０、交流電流モニタ、電流センサ１２、掃引信号発生器５０を統合していてもよい。
本実施形態の測定装置１００は、後述するインピーダンス演算器８０及びインピーダンス演算器８０と信号を入出力するための信号入出力インターフェイス（図示せず）を備える。インピーダンス演算器８０は、また、交流電圧モニタ６０及び交流電流モニタ７０からの信号を入力するためのインターフェイスを備える。測定装置１００はインピーダンス演算器８０からなるものでもよいが、交流電圧モニタ６０、交流電流モニタ、電流センサ１２、掃引信号発生器５０を統合していてもよい。
【0016】
測定装置１００の使用形態を示す図１では、二次電池１０と電流センサ１２の直列回路の両端には、掃引信号発生器５０が接続されている。掃引信号発生器５０は、出力周波数が掃引変化する交流信号を、二次電池１０と電流センサ１２の直列回路に出力する。
【0017】
交流電圧モニタ６０は、二次電池１０の両端の交流電圧を測定してインピーダンス演算器８０に入力する。交流電流モニタ７０は、電流センサ１２に流れる交流電流を測定してインピーダンス演算器８０に入力する。
【0018】
インピーダンス演算器８０は、掃引信号発生器５０の出力信号の各周波数における交流電圧モニタ６０の測定電圧と交流電流モニタ７０の測定電流との比である二次電池１０の複素インピーダンスを算出する。
【0019】
インピーダンス演算器８０におけるＲ／Ｃ並列回路数の同定は以下のように行われる。
（１）インピーダンス演算器８０は、横軸に掃引信号発生器５０の出力信号の周波数ｆ（Ｈｚ）（対数軸）、縦軸に位相θ（度）をプロットしてボード線図を作成する（図２）。この際、インダクタンス（Ｌ）成分と拡散（例えば、ＣＰＥ）成分のフィッティングを行い、測定結果のボード線図から除去する。
（１）インピーダンス演算器８０は、横軸に掃引信号発生器５０の出力信号の周波数ｆ（Ｈｚ）（対数軸）、縦軸に位相θ（度）をプロットしてボード線図を作成する（図２）。この際、インダクタンス（Ｌ）成分と拡散（例えば、ＣＰＥ）成分のフィッティングを行い、測定結果のボード線図から除去する。
【0020】
（２）次に、インピーダンス演算器８０は、作成したボード線図上で極大かつ最大の点（以下「頂点」ともいう。）Ｔ１を抽出し、Ｒ／Ｃ並列回路数として１をカウントする。インピーダンス演算器８０は、点Ｔ１を中心にＲ／Ｃ並列回路で表すことができる山の形状をフィッティングする（図３）。なお、測定ノイズと極大かつ最大の点とを見分けるため、棄却するθ値（θ０）が設定されており、インピーダンス演算器８０は、位相（θ）の値が予め定めたθ０以上である点から極大かつ最大の点を抽出する。
【0021】
（３）インピーダンス演算器８０は、次に、フィッティングしたＲ／Ｃ並列回路で表現できるボード線図成分を、測定結果のボード線図から減ずる（図４）。すなわち、インピーダンス演算器８０は、同じ周波数（Ｈｚ）の位相（θ）の値を引き算する。
【0022】
（４）インピーダンス演算器８０が、ボード線図上で極大かつ最大の点（頂点）を抽出しなくなればインピーダンス演算器８０によるＲ／Ｃ並列回路数の同定を終了するが、まだ、極大かつ最大の点（頂点）を抽出すれば、（２）及び（３）を繰り返す。
ボード線図上の頂点の高さが予め設定した位相（θ）の大きさが予め定めたθ１（例えば、０．５）以下になったところでＲ／Ｃ並列回路数のカウント増加を止める。
ボード線図上の頂点の高さが予め設定した位相（θ）の大きさが予め定めたθ１（例えば、０．５）以下になったところでＲ／Ｃ並列回路数のカウント増加を止める。
【0023】
（２－２）インピーダンス演算器８０は、ボード線図上でさらに頂点Ｔ２を抽出し、Ｒ／Ｃ並列回路数として２をカウントする。頂点Ｔ２を中心にＲ／Ｃ並列回路で表すことができる山の形状をフィッティングする（図５）。
【0024】
（３－２）インピーダンス演算器８０は、フィッティングしたＲ／Ｃ並列回路で表現できるボード線図成分を、測定結果のボード線図から減ずる（図６）。
【0025】
（４－２）インピーダンス演算器８０が、ボード線図上で極大かつ最大の点（頂点）を抽出しなくなればインピーダンス演算器８０によるＲ／Ｃ並列回路数の同定を終了するが、まだ、極大かつ最大の点を抽出すれば、さらに（２）及び（３）を繰り返す。
【0026】
（２－３）インピーダンス演算器８０は、ボード線図上で頂点Ｔ３を抽出し、Ｒ／Ｃ並列回路数として３をカウントする。頂点Ｔ３を中心にＲ／Ｃ並列回路で表すことができる山の形状をフィッティングする（図７）。
【0027】
（３－３）インピーダンス演算器８０は、フィッティングしたＲ／Ｃ並列回路で表現できるボード線図成分を、測定結果のボード線図から減ずる（図８）。
【0028】
（４－３）インピーダンス演算器８０が、ボード線図上で極大かつ最大の点（頂点）を抽出しなくなればインピーダンス演算器８０によるＲ／Ｃ並列回路数の同定を終了するが、まだ、極大かつ最大の点を抽出すれば、さらに（２）及び（３）を繰り返す。
図８では、インピーダンス演算器８０はさらに頂点を抽出しなくなるので、ここでＲ／Ｃ並列回路数の同定を終了する。
図８では、インピーダンス演算器８０はさらに頂点を抽出しなくなるので、ここでＲ／Ｃ並列回路数の同定を終了する。
【0029】
（５）インピーダンス演算器８０は、Ｒ／Ｃ並列回路数の総カウント数を報告する。本実施形態のインピーダンス演算器８０は、図２に示すボード線図から３つの頂点を抽出できる。すなわち、本発明のインピーダンス測定装置１００は、図２に示すボード線図から、Ｒ／Ｃ並列回路数を３と同定することができる。
【0030】
従来のインピーダンス測定装置によるＲ／Ｃ並列回路数の同定結果を図９に示す。従来のインピーダンス測定装置は、図２と同一のボード線図から、Ｒ／Ｃ並列回路数を２と同定することができる。
したがって、本発明のインピーダンス測定装置は、従来、合成されたインピーダンスデータでは抽出できず見落としていた円弧成分（Ｒ／Ｃ並列回路）を抽出することができ、より精確なインピーダンス測定が可能となる。
したがって、本発明のインピーダンス測定装置は、従来、合成されたインピーダンスデータでは抽出できず見落としていた円弧成分（Ｒ／Ｃ並列回路）を抽出することができ、より精確なインピーダンス測定が可能となる。
【0031】
本実施形態のインピーダンス測定装置は、インピーダンス演算器及び信号入出力インターフェイスを備える。この構成要素は、例えば、ＣＰＵ（Central Processing Unit）などのハードウェアプロセッサがプログラム（ソフトウェア）を実行することにより実現される。これらの構成要素のうち一部又は全部は、ＬＳＩ（Large Scale Integration）やＡＳＩＣ（Application Specific Integrated Circuit）、ＦＰＧＡ（Field-Programmable Gate Array）、ＧＰＵ（Graphics Processing Unit）などのハードウェア（回路部；circuitryを含む）によって実現されてもよいし、ソフトウェアとハードウェアの協働によって実現されてもよい。プログラムは、予めＨＤＤ（Hard Disk Drive）やフラッシュメモリなどの記憶装置（非一過性の記憶媒体を備える記憶装置）に格納されていてもよいし、ＤＶＤやＣＤ－ＲＯＭなどの着脱可能な記憶媒体（非一過性の記憶媒体）に格納されており、記憶媒体がドライブ装置に装着されることでインストールされてもよい。
【0032】
上記説明した実施形態は、以下のように表現することができる。
コンピュータによって読み込み可能な命令（computer-readable instructions）を格納する記憶媒体（storage medium）と、
前記記憶媒体に接続されたプロセッサと、を備え、
前記プロセッサは、前記コンピュータによって読み込み可能な命令を実行することにより（the processor executing the computer-readable instructions to:）
横軸に周波数ｆ（Ｈｚ）（対数軸）、縦軸に位相θ（度）をプロットしてボード線図を作成し、
作成したボード線図上で極大かつ最大の点（以下「頂点」ともいう。）Ｔ１を抽出して、Ｒ／Ｃ並列回路数として１をカウントし、点Ｔ１を中心にＲ／Ｃ並列回路で表すことができる山の形状をフィッティングし、
フィッティングしたＲ／Ｃ並列回路で表現できるボード線図成分を、測定結果のボード線図から減じ、
ボード線図上で極大かつ最大の点（頂点）を抽出しなくなればＲ／Ｃ並列回路数の同定を終了してＲ／Ｃ並列回路数の総カウント数を報告し、さらに極大かつ最大の点（頂点）を抽出すれば、（ステップＳ２）及び（ステップＳ３）を繰り返す、
装置。
コンピュータによって読み込み可能な命令（computer-readable instructions）を格納する記憶媒体（storage medium）と、
前記記憶媒体に接続されたプロセッサと、を備え、
前記プロセッサは、前記コンピュータによって読み込み可能な命令を実行することにより（the processor executing the computer-readable instructions to:）
横軸に周波数ｆ（Ｈｚ）（対数軸）、縦軸に位相θ（度）をプロットしてボード線図を作成し、
作成したボード線図上で極大かつ最大の点（以下「頂点」ともいう。）Ｔ１を抽出して、Ｒ／Ｃ並列回路数として１をカウントし、点Ｔ１を中心にＲ／Ｃ並列回路で表すことができる山の形状をフィッティングし、
フィッティングしたＲ／Ｃ並列回路で表現できるボード線図成分を、測定結果のボード線図から減じ、
ボード線図上で極大かつ最大の点（頂点）を抽出しなくなればＲ／Ｃ並列回路数の同定を終了してＲ／Ｃ並列回路数の総カウント数を報告し、さらに極大かつ最大の点（頂点）を抽出すれば、（ステップＳ２）及び（ステップＳ３）を繰り返す、
装置。
【0033】
［測定方法］
本実施形態の測定方法は、以下のステップＳ１～Ｓ５を備える。
本実施形態の測定方法の説明では、図１０に示すフローチャート、図１の構成図、及び図２～４のボード線図を参照する。また、交流電圧モニタ６０及び交流電流モニタ７０は、測定データをインピーダンス演算器８０に入力しているものとする。
本実施形態の測定方法は、以下のステップＳ１～Ｓ５を備える。
本実施形態の測定方法の説明では、図１０に示すフローチャート、図１の構成図、及び図２～４のボード線図を参照する。また、交流電圧モニタ６０及び交流電流モニタ７０は、測定データをインピーダンス演算器８０に入力しているものとする。
【0034】
（ステップＳ１）インピーダンス演算器８０に搭載されたプロセッサ（以下、単に「プロセッサ」という。）が、横軸に周波数ｆ（Ｈｚ）（対数軸）、縦軸に位相θ（度）をプロットしてボード線図を作成する（図２）。この際、インダクタンス（Ｌ）成分と拡散（例えば、ＣＰＥ）成分のフィッティングを行い、測定結果のボード線図から除去する。
【0035】
（ステップＳ２）プロセッサが、作成したボード線図上で極大かつ最大の点（以下「頂点」ともいう。）Ｔ１を抽出し、Ｒ／Ｃ並列回路数として＋１をカウントする。プロセッサが、点Ｔ１を中心にＲ／Ｃ並列回路で表すことができる山の形状をフィッティングする（図３）。なお、測定ノイズと極大かつ最大の点とを見分けるため、棄却するθ値（θ０）が設定されており、プロセッサは、θ値がθ０以上である点から極大かつ最大の点（頂点）を抽出する。
【0036】
（ステップＳ３）プロセッサが、次に、フィッティングしたＲ／Ｃ並列回路で表現できるボード線図成分を、測定結果のボード線図から減ずる（図４）。すなわち、プロセッサが、同じ周波数（Ｈｚ）の位相（θ）の値を引き算する。
【0037】
（ステップＳ４）プロセッサが、ボード線図上で極大かつ最大の点（頂点）を抽出しなくなればステップＳ５に進む。インピーダンス演算器８０が、さらに極大かつ最大の点（頂点）を抽出すれば、（ステップＳ２）及び（ステップＳ３）を繰り返す。
【0038】
（ステップＳ５）プロセッサが、Ｒ／Ｃ並列回路数の総カウント数を報告して、Ｒ／Ｃ並列回路数の同定を終了する。
【0039】
［プログラム］
本実施形態のプログラムは、上述したインピーダンス演算器に搭載されたプロセッサに以下の（Ｓ１）～（Ｓ５）の処理を命じるプログラムである。
（Ｓ１）横軸に周波数ｆ（Ｈｚ）（対数軸）、縦軸に位相θ（度）をプロットしてボード線図を作成する。
（Ｓ２）作成したボード線図上で頂点を抽出し、Ｒ／Ｃ並列回路数として＋１をカウントし、前記頂点を中心にＲ／Ｃ並列回路で表すことができる山の形状をフィッティングする。
（Ｓ３）フィッティングしたＲ／Ｃ並列回路で表現できるボード線図成分を、測定結果のボード線図から減じる。
（Ｓ４）ボード線図上で頂点を抽出しなくなればＲ／Ｃ並列回路の総カウント数を報告してＲ／Ｃ並列回路数の同定を終了し、頂点をさらに抽出すれば、（Ｓ２）及び（Ｓ３）の処理を繰り返す。
（Ｓ５）Ｒ／Ｃ並列回路数の総カウント数を報告する。
本実施形態のプログラムは、上述したインピーダンス演算器に搭載されたプロセッサに以下の（Ｓ１）～（Ｓ５）の処理を命じるプログラムである。
（Ｓ１）横軸に周波数ｆ（Ｈｚ）（対数軸）、縦軸に位相θ（度）をプロットしてボード線図を作成する。
（Ｓ２）作成したボード線図上で頂点を抽出し、Ｒ／Ｃ並列回路数として＋１をカウントし、前記頂点を中心にＲ／Ｃ並列回路で表すことができる山の形状をフィッティングする。
（Ｓ３）フィッティングしたＲ／Ｃ並列回路で表現できるボード線図成分を、測定結果のボード線図から減じる。
（Ｓ４）ボード線図上で頂点を抽出しなくなればＲ／Ｃ並列回路の総カウント数を報告してＲ／Ｃ並列回路数の同定を終了し、頂点をさらに抽出すれば、（Ｓ２）及び（Ｓ３）の処理を繰り返す。
（Ｓ５）Ｒ／Ｃ並列回路数の総カウント数を報告する。
【0040】
以上、本発明を実施するための形態について実施形態を用いて説明したが、本発明はこうした実施形態に何等限定されるものではなく、本発明の要旨を逸脱しない範囲内において種々の変形及び置換を加えることができる。
【符号の説明】
【0041】
１０…二次電池
１２…電流センサ
５０…掃引信号発生器
６０…交流電圧モニタ
７０…交流電流モニタ
８０…インピーダンス演算器
１００…インピーダンス測定装置
１２…電流センサ
５０…掃引信号発生器
６０…交流電圧モニタ
７０…交流電流モニタ
８０…インピーダンス演算器
１００…インピーダンス測定装置
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】