JP4153745B2 - Impedance parameter estimation device - Google Patents
Impedance parameter estimation device Download PDFInfo
- Publication number
- JP4153745B2 JP4153745B2 JP2002221768A JP2002221768A JP4153745B2 JP 4153745 B2 JP4153745 B2 JP 4153745B2 JP 2002221768 A JP2002221768 A JP 2002221768A JP 2002221768 A JP2002221768 A JP 2002221768A JP 4153745 B2 JP4153745 B2 JP 4153745B2
- Authority
- JP
- Japan
- Prior art keywords
- impedance
- equivalent circuit
- test
- parameter
- estimation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Landscapes
- Measurement Of Resistance Or Impedance (AREA)
Description
【0001】
【発明の属する技術分野】
本発明は、供試インピーダンスの等価回路およびパラメータ値を推定するインピーダンスパラメータの推定装置に関するものである。
【0002】
【従来の技術】
従来、例えば電池等の起電力を含む測定対象物や起電力を含まない測定対象物の内部インピーダンスを精度よく測定する場合、測定手段として交流法が一般に用いられる。
【0003】
図10は、交流法によるインピーダンス測定回路の概略構成を示している。図において、1は測定対象物の供試インピーダンスで、この供試インピーダンス1は、図11に示すように、例えば測定対象物が電池の場合、溶液抵抗Rs、反応抵抗r1、r2、r3、二重層容量C1、C2、C3、配線等のインダクタンス分Lsからなる等価回路で表される。なお、起電力は省略してある。
【0004】
供試インピーダンス1には、電池の放電電流を決定する直流電流源2が接続されるとともに、抵抗3、交流電流計4を介して交流信号源5が接続されている。交流信号源5は、供試インピーダンス1に印加される交流信号の周波数を可変可能にしている。
【0005】
そして、インピーダンス測定に際して、供試インピーダンス1に印加する交流信号の周波数を広い範囲に渡って変化させ、それぞれの周波数における供試インピーダンス1の両端電圧を交流電圧計6で測定するとともに、抵抗3を流れる電流を交流電流計4で測定し、これらの測定値の比である複素インピーダンスをインピーダンス計測部7で求める。
【0006】
次に、供試インピーダンス1の複素インピーダンスを図12に示すように横軸にインピーダンスの実数部、縦軸にインピーダンスの虚数部を表わす複素平面にプロットして複素平面インピーダンス図(コールコール・プロット等)を作成し、この複素平面インピーダンス図から供試インピーダンス1の等価回路およびパラメータを求める。
【0007】
次に、複素平面インピーダンス図から求めた供試インピーダンス1の等価回路およびパラメータを初期値として、これら初期値から複素インピーダンスを計算により求め、この計算結果を等価回路・パラメタ入力部9よりカーブフィット部8へ入力する。そして、カーブフィット部8において、等価回路およびパラメータの初期値から求められたインピーダンスと、実際の測定結果から求められたインピーダンス測定値との二乗誤差を求める。ここで、二乗誤差を求める方法としては、周知の最急降法(「アルゴリズム辞典」1994年、共立出版、pp.270参照)やシンプレックス法(「A Simplex Method for Function Minimization」、Computer Journal、Vol.7、p.308-313)などの非線型推定手法が知られており、二乗誤差関数が最小となるまで、カーブフィット部8での動作を繰り返し、実際に測定したインピーダンスと計算値との誤差が最小となるパラメータを求める。そして、このカーブフィット部8で求められた結果が、供試インピーダンス1の等価回路およびパラメータ値の推定結果として、表示/出力部10に出力される。
【0008】
従って、このような従来のインピーダンスパラメータ推定装置によれば、例えば、供試インピーダンス1の等価回路が、図13に示すように溶液抵抗Rsと、1組の反応抵抗rと二重層容量Cだけから構成される非常に単純なものの場合は、測定で得られるインピーダンスのコールコール・プロットから簡単に等価回路およびパラメータ値を推定することができる。
【0009】
ここで、図14に示す複素平面インピーダンス図(コールコール・プロット)から図13に示す等価回路のパラメータを求める場合を説明すると、このときの等価回路のインピーダンスは、各周波数によってコールコール・プロット上でほぼ半円の軌跡を画く単純な形となる。この場合、周波数が非常に高いときのインピーダンスは、虚数部がゼロとなり、このときの実数部が溶液抵抗Rsに相当し、周波数の低いときのインピーダンスも虚数部がゼロとなり、このときの実数部はRs+rに相当する。これにより、この実数部より先程求めた溶液抵抗Rs分を引き算すれば、反応抵抗rが求められる。また、インピーダンスの虚数部の最大値およびそのときの周波数ω’より1/(ω’C)=r/2の関係があるので、二重層容量Cも、C=2/(ω’r)より求められる。そして、このようにして求めたRs、r、Cのパラメータ値は、初期値として最急降法等を採用したカーブフィット部8に与えられ、実際の測定結果から求められたインピーダンスとの二乗誤差が求められ、これらの誤差が最小となるパラメータが供試インピーダンス1の等価回路およびパラメータとして推定される。
【0010】
このようにして、供試インピーダンス1の等価回路が非常に単純な場合には、実際の測定結果から求められたインピーダンスの結果から、簡単に供試インピーダンス1の等価回路およびパラメータを推定することができる。
【0011】
【発明が解決しようとする課題】
ところが、通常、供試インピーダンス1の等価回路は、図13で述べたような単純な構成をしているものは殆どなく、複数の反応抵抗や二重層容量およびインダクタンス分等からなる複雑な構成からなっている。そして、等価回路の構成が複雑になると、実際の測定結果から求められたインピーダンスのコールコール・プロットも図14に示すような単純な形とはならず、図12で述べたような複雑な形状となる。このことは、コールコール・プロットが単純な1つの半円とならないため、コールコール・プロットから等価回路やパラメータ値を求めることは非常に困難であり、勿論、自動化も事実上行うことができない。
【0012】
また、供試インピーダンス1は、通常、等価回路を構成する反応抵抗や二重層容量の組(rとCの並列接続の個数)などが何個含まれているかを事前に知ることはできず、これらのrとCの並列接続の個数は、実際の測定結果から求められたインピーダンスのコールコール・プロットから推定するしかない。しかし、上述したように複数の半円が明確に判別できるコールコール・プロットとなることは非常に稀で、半円同士が重なってしまうような状態では、半円の個数から等価回路を正しく判断することは極めて難しく、このため、自動化が不可能であるばかりか、かなりの熟練を要する作業となり、作業者に多大な負担を与えるという問題がある。
【0013】
さらに、パラメータを推定する手段として用いられる最急降法やシンプレックス法などの非線型推定手法は、等価回路が正しく判断され、パラメータ初期値が求めるべき真の値に近ければ、非常に精度よくパラメータ値を推定することができるが、等価回路の判断が不正確で、パラメータ初期値が真の値から大きく離れていると、真の値に近づくことはなく、全く異なったパラメータ値を推定してしまう、つまり発散してしまうという問題が生じる。これは、最急降法やシンプレックス法などの非線型推定手法が、測定インピーダンスとの二乗誤差関数を最小にするのではなく、二乗誤差関数の極小値を探索するという性質によるものである。
【0014】
本発明は上記事情に鑑みてなされたもので、供試インピーダンスの等価回路およびパラメータを高精度に、迅速かつ簡単に求めることができるインピーダンスパラメータの推定装置を提供することを目的とする。
【0015】
【課題を解決するための手段】
請求項1記載の発明は、供試インピーダンスと、前記供試インピーダンスに印加される、複数の異なる周波数の交流信号を発生する交流信号源を有し、前記交流信号を印加することにより前記供試インピーダンスのインピーダンス値を測定するインピーダンス測定手段と、前記インピーダンス測定手段により測定された各周波数毎のインピーダンス値より伝達関数形式の連立方程式を作成し、この連立方程式で求められた伝達関数を状態方程式に変換するとともに、該状態方程式を対角化変換して1次形と2次形の合成に分解し、これら分解された1次形と2次形の部分をそれぞれ相当する等価回路およびパラメータに変換して前記供試インピーダンスの等価回路およびパラメータを推定する等価回路推定手段と、前記等価回路推定手段により推定された等価回路が描くインピーダンス軌跡の面積に対し予め所定の閾値が設定されるとともに、該閾値を満たす等価回路のみを選択する等価回路選択手段とを具備したことを特徴としている。
【0016】
請求項2記載の発明は、供試インピーダンスと、前記供試インピーダンスに印加される、複数の異なる周波数の交流信号を発生する交流信号源を有し、前記交流信号を印加することにより前記供試インピーダンスのインピーダンス値を測定するインピーダンス測定手段と、前記インピーダンス測定手段により測定された各周波数毎のインピーダンス値より伝達関数形式の連立方程式を作成し、この連立方程式で求められた伝達関数を状態方程式に変換するとともに、該状態方程式を対角化変換して1次形と2次形の合成に分解し、これら分解された1次形と2次形の部分をそれぞれ相当する等価回路およびパラメータに変換して前記供試インピーダンスの等価回路およびパラメータを推定する等価回路推定手段と、前記等価回路推定手段により推定された等価回路が描くインピーダンス軌跡の面積に対し予め所定の個数が設定されるとともに、面積が大きい方から所定の個数の等価回路のみを選択する等価回路選択手段とを具備したことを特徴としている。
【0017】
請求項3記載の発明は、1又は2記載の発明において、前記等価回路推定手段または前記等価回路選択手段から取得された等価回路およびパラメータを初期値とし、これら初期値からインピーダンスを求めるとともに、実際に測定したインピーダンス計測結果との誤差を求め、該誤差が最小となるときの前記供試インピーダンスの等価回路を構成するパラメータを推定するパラメータ推定手段をさらに設けたことを特徴としている。
【0018】
請求項4記載の発明は、請求項1乃至3のいずれか一記載の発明において、前記インピーダンス測定手段の交流信号源は、交流信号として掃引正弦波、ランダムノイズ、多重正弦波あるいはインパルス信号を発生するものであることを特徴としている。
【0019】
請求項5記載の発明は、請求項1乃至4のいずれか一記載の発明において、前記インピーダンス測定手段により測定された前記供試インピーダンスのインピーダンス値は、通信回線または外部記憶手段を介して取得されるものであることを特徴としている。
【0020】
この結果、本発明によれば、カーブフィット方法のようにコールコール・プロットを目視しながら等価回路およびパラメータ初期値を予想して手作業で設定するような必要がなくなり、供試インピーダンスの等価回路およびパラメータを高精度に、迅速かつ簡単に求めることができ、これら処理の自動化も実現できる。
【0021】
また、本発明によれば、意味のない不用な等価回路を除去し、有意と予想される等価回路のみを選択することができる。
【0022】
また、本発明によれば、求められた等価回路およびパラメータに対し、さらにカーブフィット計算を行うことにより、誤差の少ない等価回路およびパラメータ値を推定することができる。
【0023】
また、本発明によれば、他の場所に用意された多様な測定対象物の供試インピーダンスについても、簡単に等価回路およびパラメータ値を推定することができる。
【0024】
【発明の実施の形態】
以下、本発明の実施の形態を図面に従い説明する。
【0025】
(第1の実施の形態)
図1は、本発明が適用されるインピーダンスパラメータの推定装置の概略構成を示している。
【0026】
図において、11は供試インピーダンスで、ここでの供試インピーダンス11は、起電力を発生する電池などからなっている。勿論、供試インピーダンス11は、起電力を含まないものでもよい。
【0027】
供試インピーダンス11には、インピーダンス測定手段12が接続されている。このインピーダンス測定手段12は、供試インピーダンス11に接続される負荷装置13を有している。負荷装置13は、供試インピーダンス11をどのくらいの電流で測定するかを決める通電電流を決定するものである。
【0028】
なお、負荷装置13は、供試インピーダンス11の通電電流を設定できるものであれば、増幅器などに置き代えることができる。
【0029】
負荷装置13には、交流信号源としての掃引正弦波発生器14が接続されている。この掃引正弦波発生器14は、負荷装置13の通電電流に比べて振幅の遥かに小さい掃引正弦波状の電流を発生するもので、負荷装置13の通電電流に周波数掃引された正弦波状の変化を与える。
【0030】
供試インピーダンス11には、電圧検出部15が並列に接続されるとともに、電流検出部16が直列に接続されている。電圧検出部15は、供試インピーダンス11の端子電圧を測定し、また、電流検出部16は、供試インピーダンス11に流れる電流を測定するものである。この電流検出部16は、供試インピーダンス11と電圧検出部15との間にあってもよい。
【0031】
電圧検出部15および電流検出部16には、インピーダンス計測部17が接続されている。インピーダンス計測部17は、掃引正弦波発生器14より負荷装置13の通電電流に周波数掃引された正弦波状の変化を与えている状態で、それぞれの周波数における電圧検出部15および電流検出部16の測定値から供試インピーダンス11のインピーダンス計測値を求める。
【0032】
インピーダンス計測部17には、等価回路推定手段としての等価回路推定部18が接続されている。等価回路推定部18は、インピーダンス計測部17で計測されたインピーダンス測定値から供試インピーダンスの等価回路およびパラメータ値を推定するもので、これら等価回路およびパラメータ値の推定のための処理として、図2に示すフローチャートが実行される。
【0033】
図2において、まず、ステップ18aでは、インピーダンス計測部17で求められた各周波数毎のインピーダンス測定値より伝達関数形式の連立方程式を作成する。この場合、周波数fでのインピーダンス計測値Z( ω )を、下記の伝達関数形式にフィットする。なお、式(1)中においてs=j2πfである。また、mおよびnは、次数であり、m≦nかつ予想される回路の次数より充分に大きな値とする。
【0034】
【数1】
【0035】
インピーダンス測定手段12の測定で得られた全周波数f1〜fk(s1〜sk)でのインピーダンス計測結果Z1〜Zkをまとめて行列形式で表記すると式(2)のように表せる。
【0036】
【数2】
【0037】
次に、式(2)の右辺の一部を左辺に移行することで、以下の連立方程式を得る。
【0038】
【数3】
【0039】
そして、上式を Ax=y とおき、x=A-1・yからベクトルx、すなわちa 0 〜a n 、b1〜bmを求める。
【0040】
なお、インピーダンス計測値を伝達関数にフィットする方法としては、上記に述べた以外にも各種(例えば、E.C.LEVY、「Complex-Curve Fitting」、IRE TRANSACTIONS ON AUTOMATIC CONTROL、1959、May、p.37-43)の方法があり、何れの方法を採用してもよい。
【0041】
次に、ステップ18bに進む。このステップ18bでは、ステップ18aで求められた伝達関数を状態方程式に変換する。具体的には、ステップ18aで求めた伝達関数を、周知の方法(例えば「入門現代制御理論」、啓学出版、p30〜)を用いて下記の状態方程式に変換する。
【0042】
【数4】
【0043】
次に、ステップ18cに進む。このステップ18cでは、状態方程式の解を単純な1次形と2次形の合成に分解する。この場合、式(4)の状態方程式の状態行列を対角要素とそれに隣接する要素以外は「0」となる形(下式(5))に対角化変換(例えば、白岩謙一「基礎課程 線形代数入門」、サイエンス社、1992、181ページ〜184ページ)する。このような対角化変換することにより、測定インピーダンスに当てはめた式(1)で示す伝達関数は、ブロック線図上では図3に示すような2次形回路(対角要素のうち、2×2行列の部分)と1次形回路の並列接続、電気回路では図4に示すような直列接続として表現できる。
【0044】
【数5】
【0045】
そして、ステップ18dに進む。このステップ18dでは、分解された1次形と2次形の部分をそれぞれ相当する等価回路およびパラメータに変換する。この場合、式(5)に示す対角化変換された状態方程式は、分解された状態行列An、Bn、Cn、Dnに応じて、1次形の部分は、式(6)で表わされる図5に示すようなcpとrpの並列接続の回路あるいは式(7)で表わされる図6に示すようなLpとrpの並列接続の回路に変換される。
【0046】
【数6】
【0047】
この場合、これら2つの回路は、状態行列Bn×Cnが0未満ならばLpとrpの並列接続の回路、0以上ならcpとrpの並列接続の回路と簡単に判別できる。
【0048】
なお、ここでは最も単純な2種類の1次形の回路を示したが、インピーダンス特性が1次形の回路であれば他の回路に変換してもよい。
【0049】
同様にして、2次形の部分の回路を求める。この場合、2次形の部分状態行列An、Bn、Cn、Dnを式(8)に示すような伝達関数形に変換する。
【0050】
【数7】
【0051】
ここで、2次形の部分は、図7に示すようなLp、cp、rpの並列接続の回路に相当するので、各々対角要素に相当する回路およびパラメータ値を式(9)のように求める。
【0052】
【数8】
【0053】
なお、ここでも最も単純な2次形の回路を示したが、インピーダンス特性が2次形の回路であれば他の回路に変換してもよい。
【0054】
このような等価回路推定部18での一連の処理により、インピーダンス計測部17より求められたインピーダンス計測値から供試インピーダンス11の等価回路およびパラメータ値を推定することができる。
【0055】
このときの処理は、従来のカーブフィット方法と異なり、外部から等価回路およびパラメータ値を全く与えることなく、インピーダンス計測値だけから等価回路およびパラメータ値を推定することができる。つまり、カーブフィット方法のようにコールコール・プロットを目視しながら等価回路およびパラメータ初期値を予想して手作業で設定するようなことが全く不要になるので、供試インピーダンス11の等価回路およびパラメータを高精度に、迅速かつ簡単に求めることができ、これら処理の自動化も実現できる。
【0056】
なお、このような等価回路推定部18での供試インピーダンスの等価回路およびパラメータ値を推定する一連の処理において、ステップ18bでの伝達関数を状態方程式に変換する処理は省略することができる。つまり、ステップ18aで、インピーダンス測定手段12により測定された各周波数毎のインピーダンスより伝達関数形式の連立方程式を作成した後、直ちにステップ18cに進んで、連立方程式の解を単純な1次形と2次形の合成に分解するようにしても、上述したと同様な効果を期待できる。
【0057】
図1に戻って、等価回路推定部18には、等価回路選択手段としての等価回路選択部19が接続されている。
【0058】
ところで、等価回路推定部18で求められる等価回路およびパラメータの個数は、図2に示すステップ18aで述べた式(1)の次数mとnに比例し、通常多数の等価回路が出力される。これら求められた全ての等価回路を出力結果として採用してもよいが、この中の多くの回路は、インピーダンス計測値に与える影響が殆どない、意味のない不用な等価回路である。
【0059】
そこで、等価回路選択部19を設けて、等価回路推定部18より求められる等価回路のうち有意と予想される等価回路のみを選択し、他の意味のない等価回路を削除することにより、最終的に得られる等価回路の規模を必要最小限に抑える。これにより、装置の取り扱いを簡単にできるとともに、計算時間などの節約を実現することができる。
【0060】
この場合、等価回路選択部19で用いられる等価回路を選択する方法は、コールコール・プロット上で等価回路が占める面積の大きさに着目している。つまり、等価回路推定部18より求められる等価回路は、例えば、図8(a)(b)(c)に示すようにコールコール・プロット上で半円あるいは円周のインピーダンス軌跡を描くが、この時の扇状の面積(図中の斜線部分)に着目し、全ての等価回路の面積のうち、予め設定した閾値(例えば1%)を満たすインピーダンス軌跡を描く等価回路を有意な等価回路として採用し、閾値を満たさないインピーダンス軌跡しか描かない等価回路は、インピーダンス測定値に与える影響が無視できると判断して削除する。このときの閾値は、ユーザが任意に設定できることが望ましい。
【0061】
なお、閾値は設けずに、大きな面積を描く等価回路から順番に所定の数だけ有意な等価回路として採用するようにしてもよい。また、他の等価回路を選択する手段として、扇状の面積以外にも例えば円弧の面積などを利用することも可能である。
【0062】
図1に戻って、等価回路選択部19には、パラメータ推定手段としてのパラメータ推定部20が接続されている。このパラメータ推定部20は、インピーダンス計測部17に接続されている。
【0063】
等価回路選択部19で選択した等価回路およびパラメータ値をそのままカーブフィット結果として出力してもよいが、さらに高精度なフィット結果を必要とすることがある。
【0064】
そこで、ここでは等価回路選択部19で選択した等価回路およびパラメータ値をパラメータ推定部20に入力して、さらにカーブフィット計算を行うようにしている。この場合、等価回路選択部19で選択した等価回路およびパラメータ値をパラメータ推定部20の初期値としてそのまま利用可能である。あるいは、等価回路選択部19で選択した等価回路およびパラメータ値をユーザが任意に編集できるようにしてもよい。
【0065】
パラメータ推定部20では、等価回路およびパラメータ値の初期値とインピーダンス計測部17で求められたインピーダンス計測値との二乗誤差を求める。この二乗誤差を求める方法としては、上述したと同様に、周知の最急降下法(「アルゴリズム辞典」1994年、共立出版、pp.270参照)やシンプレックス法(「A Simplex Method for Function Minimization」、Computer Journal、Vol.7、p.308-313)などの方法が用いられ、二乗誤差関数が最小となるまで、パラメータの初期値を変更していき、インピーダンス計測値と初期値の誤差が最小となるパラメータ値を求める。このカーブフィット計算においては、パラメータ値は変更するが等価回路は変更しない。
【0066】
そして、このようにしてパラメータ推定部20で求められた結果が、供試インピーダンス1の等価回路およびパラメータ値の推定結果として出力される。
【0067】
この場合、パラメータ推定部20より出力される等価回路やパラメータ値を任意に編集して、さらにパラメータ推定部20でカーブフィット計算を繰り返し行うようにすれば、さらに高精度なカーブフィットを行うことができる。
【0068】
なお、上述では、等価回路選択部19で選択した等価回路およびパラメータ値をパラメータ推定部20の初期値としてパラメータ推定部20に与えるようにしているが、等価回路推定部18で求めた等価回路およびパラメータ値を、パラメータ推定部20の初期値として利用することも可能である。この場合、等価回路推定部18で求めた等価回路およびパラメータ値をユーザが任意に編集できる構成を備えるようにしてもよい。
【0069】
パラメータ推定部20でカーブフィット計算を行なった結果は、表示/出力部21に出力され、供試インピーダンス11の等価回路とパラメータ推定値として表示される。
【0070】
従って、このようにすれば、通電電流およびそれに重畳した周波数を掃引した正弦波状の電流を供給したときの供試インピーダンス11の端子電圧および電流をそれぞれ測定し、これらの測定結果から供試インピーダンス11のインピーダンス計測値を求め、このインピーダンス測定値を等価回路推定部18に入力し、この等価回路推定部18での一連の処理により、供試インピーダンス11の等価回路およびパラメータ値を推定するようにしている。これにより、従来のカーブフィット方法と異なり、外部から等価回路およびパラメータ値を全く与えることなく、インピーダンス計測値だけから等価回路およびパラメータ値を推定することができる。つまり、従来のカーブフィット方法のように、コールコール・プロットを目視しながら等価回路およびパラメータ初期値を予想して手作業で設定するようなことが全く不要となるので、等価回路およびパラメータ値を高精度に、迅速かつ簡単に求めることができ、さらに、これら処理の自動化も実現できる。
【0071】
また、等価回路選択部19により、意味のない不用な等価回路を除去し、有意と予想される等価回路のみを選択することができるので、最終的に得られる等価回路の規模を必要最小限に抑えることができ、取り扱いを簡単にできるとともに、計算時間などの節約を図ることもできる。
【0072】
また、パラメータ推定部20を用いて、等価回路選択部19により選択された等価回路およびパラメータに対し、さらにカーブフィット計算を行うようにしているので、誤差の少ない等価回路およびパラメータ値を求めることができる。
【0073】
なお、上述した実施の形態では、インピーダンスを測定するため、供試インピーダンス11に掃引正弦波状の電流変化を生じさせたが、掃引正弦波以外にランダムノイズや多重正弦波、インパルス信号など、インピーダンスを計測するために使用される交流信号を用いても、同様の効果が得られる。
【0074】
(第2の実施の形態)
次に、本発明の第2の実施の形態を説明する。
【0075】
図9は、本発明の第2の実施の形態の概略構成を示すもので、図1と同一部分には同符号を付している。
【0076】
この場合、第1の実施の形態で述べたインピーダンス測定手段12は、別体になっており、このインピーダンス測定手段12で求められたインピーダンス計測値は、予め外部記録手段としての磁気ディスク装置22に保存されている。そして、この磁気ディスク装置22から読み出されたインピーダンス計測値を等価回路推定部18に入力する。この場合、等価回路推定部18以降の構成と動作は、第1の実施の形態で述べたと同様なので、ここでは省略する。
【0077】
このような構成としても、第1の実施の形態と同様な効果を期待でき、加えて、他の場所に用意された多様な測定対象物の供試インピーダンスについても、簡単に等価回路およびパラメータ値を推定することができる。
【0078】
なお、上述では、インピーダンス計測値はここでは磁気ディスク装置22に保存されたものを用いる例を述べたが、他の記録媒体や通信回線からデータを読み出すものでも全く同様の効果が得られることは勿論である。
【0079】
また、インピーダンスに着目して測定や解析や等価回路変換を行ったが、インピーダンスの逆数であるアドミタンスとして扱っても同様な効果があることは勿論である。アドミタンスとして扱った場合、等価回路は1次と2次の回路の並列接続として表現できる。
【0080】
その他、本発明は、上記の実施の形態に限定されるものでなく、実施段階では、その要旨を変更しない範囲で種々変形することが可能である。
【0081】
さらに、上記実施の形態には、種々の段階の発明が含まれており、開示されている複数の構成要件における適宜な組み合わせにより種々の発明が抽出できる。例えば、実施の形態に示されている全構成要件から幾つかの構成要件が削除されても、発明が解決しようとする課題の欄で述べた課題を解決でき、発明の効果の欄で述べられている効果が得られる場合には、この構成要件が削除された構成が発明として抽出できる。
【0082】
【発明の効果】
以上述べたように本発明によれば、供試インピーダンスの等価回路およびパラメータを高精度に、迅速かつ簡単に求めることができるインピーダンスパラメータの推定装置を提供できる。
【図面の簡単な説明】
【図1】本発明の第1の実施の形態のインピーダンスパラメータの推定装置の概略構成を示す図。
【図2】第1の実施の形態に用いられる等価回路推定部の動作を説明するためのフローチャート。
【図3】第1の実施の形態に用いられる等価回路推定部を説明するための図。
【図4】第1の実施の形態に用いられる等価回路推定部を説明するための図。
【図5】第1の実施の形態に用いられる等価回路推定部を説明するための図。
【図6】第1の実施の形態に用いられる等価回路推定部を説明するための図。
【図7】第1の実施の形態に用いられる等価回路推定部を説明するための図。
【図8】第1の実施の形態に用いられる等価回路選択部を説明するための図。
【図9】本発明の第2の実施の形態の概略構成を示す図。
【図10】従来のインピーダンスパラメータの推定装置の概略構成を示す図。
【図11】従来のインピーダンスパラメータの推定装置を説明するための図。
【図12】従来のインピーダンスパラメータの推定装置を説明するための図。
【図13】従来のインピーダンスパラメータの推定装置を説明するための図。
【図14】従来のインピーダンスパラメータの推定装置を説明するための図。
【符号の説明】
1…供試インピーダンス
2…直流電流源
3…抵抗
4…交流電流計
5…交流信号源
6…交流電圧計
7…インピーダンス計測部
8…カーブフィット部
9…等価回路・パラメータ入力部
10…表示/出力部
11…供試インピーダンス
12…インピーダンス測定手段
13…負荷装置
14…掃引正弦波発生器
15…電圧検出部
16…電流検出部
17…インピーダンス計測部
18…等価回路推定部
19…等価回路選択部
20…パラメータ推定部
21…表示/出力部
22…磁気ディスク装置[0001]
BACKGROUND OF THE INVENTION
The present invention relates to an equivalent circuit for a test impedance and an impedance parameter estimation device for estimating a parameter value.
[0002]
[Prior art]
Conventionally, for example, when the internal impedance of a measurement object including an electromotive force such as a battery or a measurement object not including an electromotive force is accurately measured, an alternating current method is generally used as a measurement unit.
[0003]
FIG. 10 shows a schematic configuration of an impedance measurement circuit using an alternating current method. In the figure, reference numeral 1 denotes a test impedance of a measurement object. As shown in FIG. 11, for example, when the measurement object is a battery, the test impedance 1 is a solution resistance Rs, a reaction resistance r1, r2, r3, two. It is represented by an equivalent circuit composed of an inductance Ls such as multilayer capacitances C1, C2, C3 and wiring. The electromotive force is omitted.
[0004]
A DC
[0005]
In the impedance measurement, the frequency of the AC signal applied to the test impedance 1 is changed over a wide range, the voltage across the test impedance 1 at each frequency is measured with the
[0006]
Next, the complex impedance of the test impedance 1 is plotted on a complex plane representing the real part of the impedance on the horizontal axis and the imaginary part of the impedance on the vertical axis as shown in FIG. ) And an equivalent circuit and parameters of the test impedance 1 are obtained from this complex plane impedance diagram.
[0007]
Next, using the equivalent circuit and parameters of the test impedance 1 obtained from the complex plane impedance diagram as initial values, the complex impedance is obtained by calculation from these initial values, and the calculation result is obtained from the equivalent circuit / parameter input unit 9 by the curve fitting unit. Input to 8. Then, the
[0008]
Therefore, according to such a conventional impedance parameter estimating apparatus, for example, an equivalent circuit of the test impedance 1 is composed only of the solution resistance Rs, the set of reaction resistance r, and the double layer capacitance C as shown in FIG. For a very simple configuration, the equivalent circuit and parameter values can be easily estimated from a Cole-Cole plot of the impedance obtained from the measurement.
[0009]
Here, the case where the parameters of the equivalent circuit shown in FIG. 13 are obtained from the complex plane impedance diagram (Cole-Cole plot) shown in FIG. 14 will be described. The impedance of the equivalent circuit at this time depends on each frequency on the Cole-Cole plot. It becomes a simple shape that draws a semicircular locus. In this case, the imaginary part of the impedance when the frequency is very high is zero, the real part is equivalent to the solution resistance Rs, and the imaginary part is also zero when the frequency is low. Corresponds to Rs + r. Thus, the reaction resistance r can be obtained by subtracting the solution resistance Rs obtained previously from the real part. Further, since there is a relationship of 1 / (ω′C) = r / 2 from the maximum value of the imaginary part of the impedance and the frequency ω ′ at that time, the double layer capacitance C is also from C = 2 / (ω′r). Desired. The parameter values of Rs, r, and C obtained in this way are given to the
[0010]
In this way, when the equivalent circuit of the test impedance 1 is very simple, the equivalent circuit and parameters of the test impedance 1 can be easily estimated from the result of the impedance obtained from the actual measurement result. it can.
[0011]
[Problems to be solved by the invention]
However, in general, the equivalent circuit of the test impedance 1 hardly has a simple configuration as described in FIG. 13, and has a complicated configuration including a plurality of reaction resistances, double layer capacitances, inductance components, and the like. It has become. When the configuration of the equivalent circuit becomes complicated, the Cole-Cole plot of the impedance obtained from the actual measurement result does not have a simple shape as shown in FIG. 14, but has a complicated shape as described in FIG. It becomes. This is because the Cole-Cole plot does not become a simple semicircle, so it is very difficult to obtain an equivalent circuit and parameter values from the Cole-Cole plot, and of course, automation cannot be practically performed.
[0012]
In addition, the test impedance 1 cannot usually know in advance how many reaction resistances or double layer capacitance pairs (number of parallel connections of r and C) that constitute an equivalent circuit are included, The number of r and C connected in parallel can only be estimated from a Cole-Cole plot of impedance obtained from actual measurement results. However, as described above, it is very rare that a plurality of semicircles can be clearly identified as a Cole-Cole plot. In a state where semicircles overlap each other, the equivalent circuit is correctly determined from the number of semicircles. It is extremely difficult to do this, and therefore, not only is automation impossible, but it requires a lot of skill, and it imposes a great burden on the operator.
[0013]
Furthermore, nonlinear estimation methods such as the steepest descent method and simplex method used as a means for estimating parameters are very accurate if the equivalent circuit is correctly judged and the initial parameter value is close to the true value to be obtained. The value can be estimated, but if the equivalent circuit judgment is inaccurate and the parameter initial value is far from the true value, the true value will not be approached, and a completely different parameter value will be estimated. That is, the problem of divergence occurs. This is because nonlinear estimation methods such as the steepest descent method and the simplex method do not minimize the square error function with the measured impedance, but search for the minimum value of the square error function.
[0014]
The present invention has been made in view of the above circumstances, and an object of the present invention is to provide an impedance parameter estimation device that can quickly and easily obtain an equivalent circuit and parameters of a test impedance with high accuracy.
[0015]
[Means for Solving the Problems]
The invention according to claim 1 has a test impedance and an AC signal source that generates an AC signal of a plurality of different frequencies applied to the test impedance, and the test signal is applied by applying the AC signal. An impedance measuring means for measuring the impedance value of the impedance, and a simultaneous equation in a transfer function form from the impedance value for each frequency measured by the impedance measuring means, and the transfer function obtained by the simultaneous equation is used as a state equation. In addition, the state equation is diagonalized and decomposed to be synthesized into a primary form and a secondary form, and the decomposed primary form and secondary form parts are converted into corresponding equivalent circuits and parameters, respectively. Equivalent circuit estimation means for estimating the equivalent circuit and parameters of the test impedance;A predetermined threshold is set in advance for the area of the impedance locus drawn by the equivalent circuit estimated by the equivalent circuit estimation means, and equivalent circuit selection means for selecting only an equivalent circuit satisfying the threshold;It is characterized by comprising.
[0016]
Claim2The described inventionAn impedance having a test impedance and an AC signal source for generating a plurality of AC signals of different frequencies applied to the test impedance, and measuring an impedance value of the test impedance by applying the AC signal A simultaneous equation in a transfer function format is created from the impedance value for each frequency measured by the measuring unit and the impedance measuring unit, and the transfer function obtained by the simultaneous equation is converted into a state equation, and the state equation is Diagonal transformation is performed to break down the composite of the primary form and the secondary form, and the parts of the broken down primary form and secondary form are converted into corresponding equivalent circuits and parameters, respectively. Equivalent circuit estimation means for estimating a circuit and parameters;A predetermined number is set in advance with respect to the area of the impedance locus drawn by the equivalent circuit estimated by the equivalent circuit estimation unit, and only a predetermined number of equivalent circuits are selected from the larger area.And equippedIt is characterized by that.
[0017]
Claim3The invention described is 1Or 2In the described invention, the equivalent circuit and the parameter acquired from the equivalent circuit estimation unit or the equivalent circuit selection unit are set as initial values, and an impedance is obtained from these initial values and an error from an actually measured impedance measurement result is obtained. Further, it is characterized by further comprising parameter estimation means for estimating a parameter constituting the equivalent circuit of the test impedance when the error is minimized.
[0018]
Claim4The invention described in claims 1 to3In the invention according to any one of the above, the AC signal source of the impedance measuring means generates a swept sine wave, random noise, multiple sine wave or impulse signal as an AC signal.
[0019]
Claim5The invention described in claims 1 to4In the invention according to any one of the above, the impedance value of the test impedance measured by the impedance measuring means is obtained via a communication line or an external storage means.
[0020]
As a result, according to the present invention, there is no need to manually set the equivalent circuit and parameter initial values while visually observing the Cole-Cole plot as in the curve fitting method. And parameters can be obtained quickly and easily with high accuracy, and automation of these processes can also be realized.
[0021]
Further, according to the present invention, it is possible to remove meaningless unnecessary equivalent circuits and select only equivalent circuits expected to be significant.
[0022]
Further, according to the present invention, it is possible to estimate an equivalent circuit and a parameter value with less error by further performing curve fitting calculation on the obtained equivalent circuit and parameter.
[0023]
Further, according to the present invention, it is possible to easily estimate the equivalent circuit and parameter values for the test impedances of various measurement objects prepared in other places.
[0024]
DETAILED DESCRIPTION OF THE INVENTION
Hereinafter, embodiments of the present invention will be described with reference to the drawings.
[0025]
(First embodiment)
FIG. 1 shows a schematic configuration of an impedance parameter estimation apparatus to which the present invention is applied.
[0026]
In the figure,
[0027]
Impedance measuring means 12 is connected to the
[0028]
Note that the
[0029]
A sweep
[0030]
A
[0031]
An
[0032]
An equivalent
[0033]
In FIG. 2, first, in
[0034]
[Expression 1]
[0035]
The total frequency f obtained by the measurement of the impedance measuring means 121~ Fk(S1~ Sk) Impedance measurement result Z1~ ZkCan be expressed as a formula (2).
[0036]
[Expression 2]
[0037]
Next, the following simultaneous equations are obtained by shifting a part of the right side of Expression (2) to the left side.
[0038]
[Equation 3]
[0039]
And the above equation is set as Ax = y, and x = A-1Y to vector x, ie a 0 ~ A n , B1~ BmAsk for.
[0040]
In addition to the methods described above, impedance measurement values can be fitted to transfer functions in various ways (for example, ECLEVY, “Complex-Curve Fitting”, IRE TRANSACTIONS ON AUTOMATIC CONTROL, 1959, May, p. 37- 43), and any method may be adopted.
[0041]
Next, the process proceeds to step 18b. In
[0042]
[Expression 4]
[0043]
Next, the process proceeds to step 18c. In
[0044]
[Equation 5]
[0045]
Then, the process proceeds to Step 18d. In
[0046]
[Formula 6]
[0047]
In this case, these two circuits can be easily identified as a circuit in which Lp and rp are connected in parallel if the state matrix Bn × Cn is less than 0, and a circuit in which cp and rp are connected in parallel if the state matrix Bn × Cn is 0 or more.
[0048]
Although the two simplest primary circuits are shown here, the circuits may be converted to other circuits as long as the impedance characteristic is a primary circuit.
[0049]
Similarly, a circuit of a quadratic part is obtained. In this case, the quadratic partial state matrices An, Bn, Cn, and Dn are converted into a transfer function form as shown in Expression (8).
[0050]
[Expression 7]
[0051]
Here, since the quadratic part corresponds to a parallel connection circuit of Lp, cp, and rp as shown in FIG. 7, the circuit and parameter values corresponding to the diagonal elements are expressed as in equation (9). Ask.
[0052]
[Equation 8]
[0053]
Although the simplest secondary circuit is shown here, it may be converted to another circuit as long as the impedance characteristic is a secondary circuit.
[0054]
Through such a series of processes in the equivalent
[0055]
In this process, unlike the conventional curve fitting method, the equivalent circuit and parameter values can be estimated from only the impedance measurement values without giving any equivalent circuit and parameter values from the outside. In other words, it is not necessary to manually set the equivalent circuit and parameter initial values while visually observing the Cole-Cole plot as in the curve fitting method, so that the equivalent circuit and parameters of the
[0056]
In the series of processes for estimating the equivalent circuit of the test impedance and the parameter values in the equivalent
[0057]
Returning to FIG. 1, the equivalent
[0058]
Incidentally, the number of equivalent circuits and parameters obtained by the equivalent
[0059]
Therefore, the equivalent
[0060]
In this case, the method of selecting an equivalent circuit used in the equivalent
[0061]
Instead of providing a threshold value, a predetermined number of significant circuits may be employed in order from an equivalent circuit that draws a large area. Further, as means for selecting another equivalent circuit, for example, the area of an arc can be used in addition to the fan-shaped area.
[0062]
Returning to FIG. 1, the equivalent
[0063]
The equivalent circuit and parameter values selected by the equivalent
[0064]
Therefore, here, the equivalent circuit selected by the equivalent
[0065]
The
[0066]
Then, the result obtained by the
[0067]
In this case, if the equivalent circuit or parameter value output from the
[0068]
In the above description, the equivalent circuit selected by the equivalent
[0069]
The result of the curve fitting calculation performed by the
[0070]
Accordingly, in this way, the terminal voltage and current of the
[0071]
In addition, since the equivalent
[0072]
Further, since the curve fitting calculation is further performed on the equivalent circuit and parameter selected by the equivalent
[0073]
In the above-described embodiment, in order to measure the impedance, a swept sine wave-like current change is generated in the
[0074]
(Second Embodiment)
Next, a second embodiment of the present invention will be described.
[0075]
FIG. 9 shows a schematic configuration of the second embodiment of the present invention, and the same parts as those in FIG.
[0076]
In this case, the impedance measuring means 12 described in the first embodiment is a separate body, and the impedance measurement value obtained by the impedance measuring means 12 is previously stored in the
[0077]
Even with such a configuration, the same effect as that of the first embodiment can be expected. In addition, for the test impedances of various measurement objects prepared in other places, an equivalent circuit and parameter values can be easily obtained. Can be estimated.
[0078]
In the above description, an example in which the impedance measurement value stored in the
[0079]
In addition, measurement, analysis, and equivalent circuit conversion are performed focusing on the impedance, but it is a matter of course that the same effect can be obtained by treating it as admittance that is the reciprocal of impedance. When treated as admittance, an equivalent circuit can be expressed as a parallel connection of a primary circuit and a secondary circuit.
[0080]
In addition, the present invention is not limited to the above-described embodiments, and various modifications can be made without departing from the scope of the invention in the implementation stage.
[0081]
Furthermore, the above embodiments include inventions at various stages, and various inventions can be extracted by appropriately combining a plurality of disclosed constituent elements. For example, even if some constituent requirements are deleted from all the constituent requirements shown in the embodiment, the problem described in the column of the problem to be solved by the invention can be solved, and is described in the column of the effect of the invention. If the above effect is obtained, a configuration from which this configuration requirement is deleted can be extracted as an invention.
[0082]
【The invention's effect】
As described above, according to the present invention, it is possible to provide an impedance parameter estimation device that can quickly and easily obtain an equivalent circuit and parameters of a test impedance with high accuracy.
[Brief description of the drawings]
FIG. 1 is a diagram showing a schematic configuration of an impedance parameter estimation apparatus according to a first embodiment of the present invention.
FIG. 2 is a flowchart for explaining the operation of an equivalent circuit estimation unit used in the first embodiment;
FIG. 3 is a diagram for explaining an equivalent circuit estimation unit used in the first embodiment;
FIG. 4 is a diagram for explaining an equivalent circuit estimation unit used in the first embodiment;
FIG. 5 is a diagram for explaining an equivalent circuit estimation unit used in the first embodiment;
FIG. 6 is a diagram for explaining an equivalent circuit estimation unit used in the first embodiment;
FIG. 7 is a diagram for explaining an equivalent circuit estimation unit used in the first embodiment;
FIG. 8 is a diagram for explaining an equivalent circuit selection unit used in the first embodiment;
FIG. 9 is a diagram showing a schematic configuration of a second embodiment of the present invention.
FIG. 10 is a diagram showing a schematic configuration of a conventional impedance parameter estimation apparatus.
FIG. 11 is a diagram for explaining a conventional impedance parameter estimation apparatus;
FIG. 12 is a diagram for explaining a conventional impedance parameter estimation apparatus;
FIG. 13 is a diagram for explaining a conventional impedance parameter estimation apparatus;
FIG. 14 is a diagram for explaining a conventional impedance parameter estimation apparatus;
[Explanation of symbols]
1 ... Test impedance
2 ... DC current source
3 ... resistance
4 ... AC ammeter
5 ... AC signal source
6 ... AC voltmeter
7: Impedance measurement unit
8 ... Curve fitting part
9 ... Equivalent circuit / parameter input section
10. Display / output unit
11 ... Test impedance
12: Impedance measuring means
13 ... Load device
14 ... Sweep sine wave generator
15 ... Voltage detector
16 ... Current detector
17 ... Impedance measurement section
18 ... Equivalent circuit estimation unit
19 ... Equivalent circuit selector
20: Parameter estimation unit
21 ... Display / output unit
22 ... Magnetic disk unit
Claims (5)
前記供試インピーダンスに印加される、複数の異なる周波数の交流信号を発生する交流信号源を有し、前記交流信号を印加することにより前記供試インピーダンスのインピーダンス値を測定するインピーダンス測定手段と、
前記インピーダンス測定手段により測定された各周波数毎のインピーダンス値より伝達関数形式の連立方程式を作成し、この連立方程式で求められた伝達関数を状態方程式に変換するとともに、該状態方程式を対角化変換して1次形と2次形の合成に分解し、これら分解された1次形と2次形の部分をそれぞれ相当する等価回路およびパラメータに変換して前記供試インピーダンスの等価回路およびパラメータを推定する等価回路推定手段と、
前記等価回路推定手段により推定された等価回路が描くインピーダンス軌跡の面積に対し予め所定の閾値が設定されるとともに、該閾値を満たす等価回路のみを選択する等価回路選択手段と
を具備したことを特徴とするインピーダンスパラメータの推定装置。The test impedance,
An impedance measurement means for measuring an impedance value of the test impedance by applying the AC signal, the AC signal source generating an AC signal of a plurality of different frequencies applied to the test impedance;
Create simultaneous equations in the form of transfer functions from the impedance values for each frequency measured by the impedance measuring means, convert the transfer functions obtained by these simultaneous equations into state equations, and diagonalize the state equations Then, it is decomposed into a composite of the primary form and the secondary form, and the parts of the decomposed primary form and secondary form are converted into equivalent circuits and parameters, respectively, and the equivalent circuit and parameters of the test impedance are changed. Equivalent circuit estimation means for estimation ;
A predetermined threshold is set in advance for the area of the impedance locus drawn by the equivalent circuit estimated by the equivalent circuit estimation means, and equivalent circuit selection means for selecting only an equivalent circuit satisfying the threshold is provided. Impedance parameter estimation device.
前記供試インピーダンスに印加される、複数の異なる周波数の交流信号を発生する交流信号源を有し、前記交流信号を印加することにより前記供試インピーダンスのインピーダンス値を測定するインピーダンス測定手段と、
前記インピーダンス測定手段により測定された各周波数毎のインピーダンス値より伝達関数形式の連立方程式を作成し、この連立方程式で求められた伝達関数を状態方程式に変換するとともに、該状態方程式を対角化変換して1次形と2次形の合成に分解し、これら分解された1次形と2次形の部分をそれぞれ相当する等価回路およびパラメータに変換して前記供試インピーダンスの等価回路およびパラメータを推定する等価回路推定手段と、
前記等価回路推定手段により推定された等価回路が描くインピーダンス軌跡の面積に対し予め所定の個数が設定されるとともに、面積が大きい方から所定の個数の等価回路のみを選択する等価回路選択手段と、
を具備したことを特徴とするインピーダンスパラメータの推定装置。 The test impedance,
An impedance measurement means for measuring an impedance value of the test impedance by applying the AC signal, the AC signal source generating an AC signal of a plurality of different frequencies applied to the test impedance;
Create simultaneous equations in the form of transfer functions from the impedance values for each frequency measured by the impedance measuring means, convert the transfer functions obtained by these simultaneous equations into state equations, and diagonalize the state equations Then, it is decomposed into a composite of the primary form and the secondary form, and the parts of the decomposed primary form and secondary form are converted into equivalent circuits and parameters, respectively, and the equivalent circuit and parameters of the test impedance are changed. Equivalent circuit estimation means for estimation;
A predetermined number is set in advance with respect to the area of the impedance locus drawn by the equivalent circuit estimated by the equivalent circuit estimation means, and equivalent circuit selection means for selecting only a predetermined number of equivalent circuits from the larger area ;
An impedance parameter estimation apparatus comprising:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002221768A JP4153745B2 (en) | 2002-07-30 | 2002-07-30 | Impedance parameter estimation device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002221768A JP4153745B2 (en) | 2002-07-30 | 2002-07-30 | Impedance parameter estimation device |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2004061372A JP2004061372A (en) | 2004-02-26 |
JP4153745B2 true JP4153745B2 (en) | 2008-09-24 |
Family
ID=31941992
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2002221768A Expired - Fee Related JP4153745B2 (en) | 2002-07-30 | 2002-07-30 | Impedance parameter estimation device |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4153745B2 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102901929A (en) * | 2011-07-28 | 2013-01-30 | 横河电机株式会社 | Device for calculating impedances of battery cell and battery impedance measuring system |
Families Citing this family (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4635249B2 (en) * | 2004-08-26 | 2011-02-23 | 学校法人 芝浦工業大学 | Winding test equipment |
JP4946131B2 (en) * | 2006-03-29 | 2012-06-06 | 横河電機株式会社 | Apparatus and method for measuring characteristics of fuel cell |
JP5004815B2 (en) * | 2008-02-01 | 2012-08-22 | 三井造船株式会社 | Evaluation method of silicon thin film |
JP5246030B2 (en) | 2008-09-26 | 2013-07-24 | 富士通株式会社 | Circuit automatic design program, method and apparatus |
KR101047113B1 (en) | 2008-11-11 | 2011-07-07 | 고려대학교 산학협력단 | Identity Parameter Extraction Method and Equivalent Circuit Determination Method in Electrical Circuit System |
JP5362421B2 (en) * | 2009-04-17 | 2013-12-11 | 株式会社エヌエフ回路設計ブロック | Equivalent circuit element constant estimation method, apparatus and equivalent circuit element constant estimation program, characteristic measurement method, apparatus and characteristic measurement program |
JP5267398B2 (en) | 2009-09-25 | 2013-08-21 | 富士通株式会社 | Pareto data generation program, method and apparatus for automatic circuit design, and automatic circuit design program, method and apparatus |
JP5626869B2 (en) * | 2010-01-26 | 2014-11-19 | 日置電機株式会社 | measuring device |
JP5353764B2 (en) | 2010-03-02 | 2013-11-27 | 富士通株式会社 | Automatic design support program, method and apparatus |
JP2013527613A (en) * | 2010-05-18 | 2013-06-27 | エスエムエー ソーラー テクノロジー アーゲー | Photovoltaic system and method for diagnosing contact of apparatus |
JP5934479B2 (en) * | 2011-07-25 | 2016-06-15 | 日置電機株式会社 | Equivalent circuit analysis apparatus and equivalent circuit analysis method |
JP5719723B2 (en) * | 2011-08-26 | 2015-05-20 | 日置電機株式会社 | Measuring apparatus, measuring method, and program |
JP5766551B2 (en) * | 2011-08-29 | 2015-08-19 | 日置電機株式会社 | Equivalent circuit parameter measuring apparatus and equivalent circuit parameter measuring method |
JP6046418B2 (en) * | 2012-08-27 | 2016-12-14 | 日置電機株式会社 | Equivalent circuit analysis apparatus and equivalent circuit analysis method |
JP2015117971A (en) * | 2013-12-17 | 2015-06-25 | 新電元工業株式会社 | Device and method for measuring impedance |
US9428071B2 (en) * | 2014-01-14 | 2016-08-30 | Ford Global Technologies, Llc | Impedance based battery parameter estimation |
CN107018474B (en) * | 2017-05-24 | 2024-04-12 | 上海传英信息技术有限公司 | Intelligent power amplifier calibration method and intelligent power amplifier calibration device |
KR102442474B1 (en) * | 2018-11-30 | 2022-09-14 | 주식회사 엘지에너지솔루션 | Apparatus and method for measuring internal resistance of battery cell |
-
2002
- 2002-07-30 JP JP2002221768A patent/JP4153745B2/en not_active Expired - Fee Related
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102901929A (en) * | 2011-07-28 | 2013-01-30 | 横河电机株式会社 | Device for calculating impedances of battery cell and battery impedance measuring system |
CN102901929B (en) * | 2011-07-28 | 2016-01-20 | 横河电机株式会社 | For calculating device and the battery impedance measuring system of cell impedance |
Also Published As
Publication number | Publication date |
---|---|
JP2004061372A (en) | 2004-02-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4153745B2 (en) | Impedance parameter estimation device | |
JP3190313B2 (en) | Method and apparatus for quantifying characteristic factors of power storage device using nonlinear equivalent circuit model | |
JP2782062B2 (en) | Pulse-based impedance meter and pulse-based method for measuring complex impedance | |
NL2002665C2 (en) | Time-domain reflectometer. | |
US7383140B2 (en) | Capacitance, inductance and impedance measurements using multi-tone stimulation and DSP algorithms | |
Woodington et al. | Behavioral model analysis of active harmonic load-pull measurements | |
JP6462948B2 (en) | Energy storage cell impedance measuring apparatus, method and related system | |
JP2009522572A (en) | Method and apparatus for determining non-linear behavior | |
CN108051661B (en) | Diode parasitic parameter extraction method based on ringing current time-frequency domain characteristics | |
Akçay et al. | Identification of power transformer models from frequency response data: A case study | |
JP3782026B2 (en) | Impedance parameter estimation method and apparatus | |
WO2000010017A1 (en) | Method and apparatus for differentiating contributions of harmonic pollution sources | |
CN114325057B (en) | PT primary voltage reconstruction method based on inverse black box and inverse electromagnetic dual model | |
US6519538B1 (en) | Method for predicting stability characteristics of power supplies | |
CN106199285A (en) | The arbitrarily capacitance characteristic under AC carrier measures equipment and measuring method thereof | |
Negri et al. | Inductively Coupled In-Circuit Measurement of Multiport Admittance Parameters | |
JP4845133B2 (en) | High-order harmonic resonance frequency characteristic estimation method | |
TWI281029B (en) | Automatic compensation method of AC leakage current, constituent circuits, and testing device using the circuit | |
JPS61267102A (en) | Plant modeling device | |
US6754616B1 (en) | Method of emulating an ideal transformer valid from DC to infinite frequency | |
JP4511272B2 (en) | RECORDING MEDIUM RECORDING EQUIVALENT CIRCUIT MODEL FOR STORAGE ELEMENT, RECORDING PROGRAM, RECORDING MEDIUM, derivation device, SIMULATION PROGRAM, RECORDING MEDIUM, SIMULATION DEVICE, DESIGN METHOD, GOOD / FAILURE JUDGING METHOD, AND GOOD / FAILURE JUDGING DEVICE | |
De Paula et al. | A simple and accurate cable modeling suitable for high-frequency phenomena analysis in PWM motor drives | |
JP7415089B1 (en) | Noise analysis device, noise analysis method, and program | |
JP7279489B2 (en) | Measuring device and its measuring method | |
JPH08285909A (en) | Nonlinear inductance model forming method for arrester |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20050520 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20051006 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20070420 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20070508 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20070706 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20071113 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20071226 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20080624 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20080704 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 4153745 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110711 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110711 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110711 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110711 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120711 Year of fee payment: 4 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120711 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120711 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120711 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130711 Year of fee payment: 5 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130711 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130711 Year of fee payment: 5 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
LAPS | Cancellation because of no payment of annual fees |