(57)【特許請求の範囲】
【請求項１】
供試インピーダンスに流れる通電電流がその基準値である直流電流に矩形波電流の加算により矩形波状の変化を有し、前記通電電流により前記供試インピーダンスから得られる電圧データを時系列的に求め、これら電圧データから前記供試インピーダンスの等価回路を構成するパラメータを推定することを特徴とするインピーダンスパラメータの推定方法。
【請求項２】
供試インピーダンスに流れる通電電流がその基準値である直流電流に矩形波電流の加算により矩形波状の変化を有し、前記通電電流を以て前記供試インピーダンスから得られる電圧データを時系列的に求め、得られた複数の電圧データの平均値を求め、この平均値から前記供試インピーダンスの等価回路を構成するパラメータを推定することを特徴とするインピーダンスパラメータの推定方法。
【請求項３】
供試インピーダンスに流れる通電電流がその基準値である直流電流に矩形波電流の加算により矩形波状の変化を有し、前記通電電流を以て前記供試インピーダンスから得られる電圧データを時系列的に求め、これら電圧データから前記供試インピーダンスの等価回路を構成する複数のパラメータを推定し、この推定により得られた複数のパラメータの平均値を求め、この平均値から前記供試インピーダンスの等価回路を構成するパラメータを推定することを特徴とするインピーダンスパラメータの推定方法。
【請求項４】
供試インピーダンスに流れる通電電流がその基準値である直流電流に矩形波電流の加算により矩形波状の変化を有し、前記通電電流により前記供試インピーダンスから得られる電圧データを時系列的に求め、
これら電圧データから推定される、前記供試インピーダンスの等価回路を構成するパラメータを求め、該パラメータを初期値とし、この初期値を前記等価回路を構成するパラメータに用いて前記供試インピーダンスの両端電圧を算出し、
この両端電圧と、前記供試インピーダンスから実際に時系列的に測定された両端電圧の電圧データとの誤差を算出し、
この誤差が最小となるときの前記等価回路を構成するパラメータを推定することを特徴とするインピーダンスパラメータの推定方法。
【請求項５】
供試インピーダンスに流れる通電電流がその基準値である直流電流に矩形波電流の加算により矩形波状の変化を有し、前記通電電流を以て前記供試インピーダンスから得られる電圧データを時系列的に求め、これら電圧データの平均値を求め、この平均値から前記供試インピーダンスの等価回路を構成するパラメータを推定し、この推定により得られたパラメータを初期値とし、この初期値を前記供試インピーダンスの等価回路を構成するパラメータに用いて前記供試インピーダンスの両端電圧を算出し、
この両端電圧と、前記供試インピーダンスから実際に時系列的に測定された両端電圧の電圧データとの誤差を算出し、
この誤差が最小となるときの前記供試インピーダンスの等価回路を構成するパラメータを推定することを特徴とするインピーダンスパラメータの推定方法。
【請求項６】
供試インピーダンスの等価回路を構成するパラメータを推定するインピーダンスパラメータの推定装置であって、
基準値である直流電流に矩形波電流の加算により予め定めた時間だけ変化させた通電電流を前記供試インピーダンスに供給する電流供給手段と、
前記通電電流に前記変化を与える直前の時刻、前記通電電流を変化させる前記時間内の少なくとも２点の時刻で、前記供試インピーダンスの両端電圧を測定する電圧測定手段と、
を備え、前記電圧測定手段で得られた電圧データから前記供試インピーダンスの等価回路を構成するパラメータを推定することを特徴とするインピーダンスパラメータの推定装置。
【請求項７】
供試インピーダンスの等価回路を構成するパラメータを推定するインピーダンスパラメータの推定装置であって、
基準値である直流電流に矩形波電流の加算により予め定めた時間だけ変化させた通電電流を前記供試インピーダンスに供給する電流供給手段と、
前記通電電流に前記変化を与える直前の時刻、前記通電電流に変化が与えられる前記時間内の少なくとも２点の時刻の各時刻で前記供試インピーダンスの両端電圧を測定する電圧測定手段と、
前記電圧測定手段により前記各時刻で得た電圧データを用いて算出された前記パラメータを初期値とし、この初期値を前記供試インピーダンスの等価回路を構成するパラメータに用いて前記供試インピーダンスの両端電圧を算出し、この両端電圧と前記電圧測定手段により得た前記電圧データとの誤差を最小とするカーブフィット手段と、
を備えることを特徴とするインピーダンスパラメータの推定装置。
【発明の詳細な説明】
【０００１】
【産業上の利用分野】
本発明は、等価回路で表わされた供試インピーダンスのパラメータ値を推定するインピーダンスパラメータの推定方法及び装置に関するものである。
【０００２】
【従来の技術】
従来、電池などの起電力を含む測定対象物や、あるいは起電力を含まない測定対象物の内部インピーダンスを測る場合、その内部インピーダンスを等価回路で表現し、その各パラメータを次のような方法で推定していた。
【０００３】
（Ａ）電流遮断法：
図７は電流遮断法によるインピーダンスパラメータ推定回路例である。供試インピーダンスは、図８に示すような起電力Ｖ１を含む等価回路で表わされ、例えば、溶液抵抗Ｒs、反応抵抗ｒ1及び二重層容量Ｃdから成り、等価回路中のパラメータＲsのみを測定できる。
【０００４】
図７において、直流電流源１０１からの電流は、スイッチ１０２を経由して供試インピーダンス１０３に供給され、この供試インピーダンス１０３の両端電圧が電圧測定部１０４で測定される。
【０００５】
測定時には、先ず、スイッチ１０２をオン状態として供試インピーダンス１０３に電流Iｍを流しておく。次に、キャパシタＣdと抵抗ｒ1による時定数の影響が無視できる程度の短時間だけスイッチ１０２をオフとし、その後、再びオンとする。このとき得られる供試インピーダンス１０３に流れる電流と両端の電圧の時間関係が図９に示されている。電流を遮断した前後における供試インピーダンス１０３の両端の電圧Ｖ0、Ｖ1を電圧測定部１０４で測定し、得られた値を用いて供試インピーダンスのパラメータＲsを推定する。
【０００６】
（Ｂ）交流法：
図１０は交流法による供試インピーダンス推定回路例である。供試インピーダンス２０５には、直流電流源２０１及び交流信号源２０４の交流電流が供給される。交流信号源２０４は、抵抗２０３と交流電流計２０２を介して供試インピーダンス２０５に接続されている。
【０００７】
測定に際しては、供試インピーダンス２０５に印加する交流の周波数を広い範囲に渡って変化させ、各周波数における供試インピーダンス２０５の両端電圧を交流電圧計２０６で測定し、その値と交流電流計２０２の測定値との比である複素インピーダンスを算出する。次に、この供試インピーダンスの複素インピーダンスを図１１に示すように複素平面にプロットして複素平面インピーダンス図（コールコール・プロット等）を作成し、このグラフから供試インピーダンス２０５のパラメータＲs、ｒ1、Ｃdを推定している。
【０００８】
【発明が解決しようとする課題】
しかしながら、上記に示す従来のインピーダンスパラメータの推定方法では次のような問題点があった。
【０００９】
上述（Ａ）の電流遮断法では、図８に示す等価回路中のパラメータＲsは推定できるものの、ｒ1やＣdは推定できない。供試インピーダンス中に起電力を含まない場合であれば、電流遮断の時間を増せば測定した電圧からｒ1も求まる。しかし、通常、供試インピーダンスは通電電流の大きさによってパラメータ値が変化することを考慮すると、電流遮断法では電流を大きく変化させるので、どの電流値のインピーダンスパラメータ値であるかが一意に定まらない。
【００１０】
また、電流遮断を行う際、スイッチ１０２の性能上の制約から図９の電圧波形の変化直後に過渡現象が生じるため、それが整定するまでの間はＶ1測定に支障が出る場合がある。
【００１１】
更に、（Ｂ）の交流法では、図８の等価回路のパラメータＲs、ｒ1、Ｃdをすべて推定できる。また、供試インピーダンスが非線形性を有し通電電流によってそのパラメータが変化する場合でも通電電流、すなわち直流電流源２０１の出力電流値は一意に確定できる。
【００１２】
しかし、一般的にコールコール・プロットを得るためには多数の周波数について測定を行わなければならないため、測定に時間がかかるという問題がある。
【００１３】
更に、得られたコールコール・プロットに対して等価回路を表す半円を手動で当てはめ、その半径とＸ軸上の中心座標の値からパラメータを推定しているが、このようにパラメータ推定において手動部分が多いため自動化が困難である。
【００１４】
そこで、本発明の目的は、通電電流に対する供試インピーダンスの等価回路パラメータを高精度且つ迅速に求めることができるインピーダンスパラメータの推定方法及び装置を提供することにある。
【００１５】
【課題を解決するための手段】
前述の課題を解決するため、本発明によるインピーダンスパラメータの推定方法及び装置は、次のような特徴的な構成を採用している。
【００１６】
（１）  供試インピーダンスに流れる通電電流がその基準値である直流電流に矩形波電流の加算により矩形波状の変化を有し、前記通電電流により前記供試インピーダンスから得られる電圧データを時系列的に求め、これら電圧データから前記供試インピーダンスの等価回路を構成するパラメータを推定することを特徴とするインピーダンスパラメータの推定方法。
（２）  供試インピーダンスに流れる通電電流がその基準値である直流電流に矩形波電流の加算により矩形波状の変化を有し、前記通電電流を以て前記供試インピーダンスから得られる電圧データを時系列的に求め、得られた複数の電圧データの平均値を求め、この平均値から前記供試インピーダンスの等価回路を構成するパラメータを推定することを特徴とするインピーダンスパラメータの推定方法。
（３）  供試インピーダンスに流れる通電電流がその基準値である直流電流に矩形波電流の加算により矩形波状の変化を有し、前記通電電流を以て前記供試インピーダンスから得られる電圧データを時系列的に求め、これら電圧データから前記供試インピーダンスの等価回路を構成する複数のパラメータを推定し、この推定により得られた複数のパラメータの平均値を求め、この平均値から前記供試インピーダンスの等価回路を構成するパラメータを推定することを特徴とするインピーダンスパラメータの推定方法。
（４）  供試インピーダンスに流れる通電電流がその基準値である直流電流に矩形波電流の加算により矩形波状の変化を有し、前記通電電流により前記供試インピーダンスから得られる電圧データを時系列的に求め、
これら電圧データから推定される、前記供試インピーダンスの等価回路を構成するパラメータを求め、該パラメータを初期値とし、この初期値を前記等価回路を構成するパラメータに用いて前記供試インピーダンスの両端電圧を算出し、
この両端電圧と、前記供試インピーダンスから実際に時系列的に測定された両端電圧の電圧データとの誤差を算出し、
この誤差が最小となるときの前記等価回路を構成するパラメータを推定することを特徴とするインピーダンスパラメータの推定方法。
（５）  供試インピーダンスに流れる通電電流がその基準値である直流電流に矩形波電流の加算により矩形波状の変化を有し、前記通電電流を以て前記供試インピーダンスから得られる電圧データを時系列的に求め、これら電圧データの平均値を求め、この平均値から前記供試インピーダンスの等価回路を構成するパラメータを推定し、この推定により得られたパラメータを初期値とし、この初期値を前記供試インピーダンスの等価回路を構成するパラメータに用いて前記供試インピーダンスの両端電圧を算出し、
この両端電圧と、前記供試インピーダンスから実際に時系列的に測定された両端電圧の電圧データとの誤差を算出し、
この誤差が最小となるときの前記供試インピーダンスの等価回路を構成するパラメータを推定することを特徴とするインピーダンスパラメータの推定方法。
（６）  供試インピーダンスの等価回路を構成するパラメータを推定するインピーダンスパラメータの推定装置であって、
基準値である直流電流に矩形波電流の加算により予め定めた時間だけ変化させた通電電流を前記供試インピーダンスに供給する電流供給手段と、
前記通電電流に前記変化を与える直前の時刻、前記通電電流を変化させる前記時間内の少なくとも２点の時刻で、前記供試インピーダンスの両端電圧を測定する電圧測定手段と、
を備え、前記電圧測定手段で得られた電圧データから前記供試インピーダンスの等価回路を構成するパラメータを推定することを特徴とするインピーダンスパラメータの推定装置。
（７）  供試インピーダンスの等価回路を構成するパラメータを推定するインピーダンスパラメータの推定装置であって、
基準値である直流電流に矩形波電流の加算により予め定めた時間だけ変化させた通電電流を前記供試インピーダンスに供給する電流供給手段と、
前記通電電流に前記変化を与える直前の時刻、前記通電電流に変化が与えられる前記時間内の少なくとも２点の時刻の各時刻で前記供試インピーダンスの両端電圧を測定する電圧測定手段と、
前記電圧測定手段により前記各時刻で得た電圧データを用いて算出された前記パラメータを初期値とし、この初期値を前記供試インピーダンスの等価回路を構成するパラメータに用いて前記供試インピーダンスの両端電圧を算出し、この両端電圧と前記電圧測定手段により得た前記電圧データとの誤差を最小とするカーブフィット手段と、
を備えることを特徴とするインピーダンスパラメータの推定装置。
【００２０】
【発明の実施の形態】
以下、本発明によるインピーダンスパラメータの推定方法及び装置の好適実施形態を添付図を参照しながら説明する。
【００２１】
＜第１の実施形態＞
図１は本発明の一実施形態による供試インピーダンスのパラメータ推定回路図である。
【００２２】
本実施形態による供試インピーダンスのパラメータ推定方法及び装置の概略を述べると次のようになる。すなわち、供試インピーダンスをどのくらいの電流で測定するかを決める通電電流を流すと同時に、これに比べて振幅のはるかに小さい矩形状電流を与え、このときの供試インピーダンス両端の電圧を時系列的に採取し、この採取データから供試インピーダンスの等価回路中のパラメータＲs、ｒ1、Ｃdの各初期値を算出する。測定データの採取はこの１回のみで良い。
【００２３】
次に、これら初期値を供試インピーダンスの等価回路パラメータに代入して供試インピーダンス両端の電圧を計算により求める。続いて、この計算値と、実際に測定した前記時系列電圧データとの二乗誤差を求める。
【００２４】
そして、周知の最急降法（「アルゴリズム辞典」１９９４年、共立出版、ｐｐ．２７０参照）により前記二乗誤差関数が最小となるまで上記操作を繰り返す。
【００２５】
図１において、供試インピーダンス３には、第１の直流電流源１からその測定ポイントを決める通電電流が供給され、更にこの通電電流に加算する形で第２の直流電流源５がスイッチ２を介して予め定めた時間（供試インピーダンスの値にもよるが、例えば０．１秒）だけ供給される。この第２の直流電流源５の電流値Imは、第１の電流源１による通電電流に比較して数％程度の十分小さい値である。これはImの値があまり大きいと測定時の通電電流の値が一意的に定まらないからである。
【００２６】
測定に際しては、直流電流源１によって測定ポイントの通電電流を設定した後、スイッチ２を短い時間だけオフし、その時点を挟んだ前後の供試インピーダンス両端の時系列電圧データを採取する。尚、スイッチ２をオフした直後の時系列電圧データのいくつかはスイッチ２の性能上の制約から生じる測定回路の過渡応答の影響を避けるため用いない。
【００２７】
供試インピーダンス３の両端電圧は電圧測定部４により測定され、時系列的電圧データとして収集してデジタルデータ化している。そして電圧測定部４で得られた時系列電圧データはカーブフィット部６に送出される。カーブフィット部６では、前記電圧測定部４による時系列電圧データと予め与えられた等価回路モデルとから、以下に述べるカーブフィット処理により供試インピーダンスの等価回路パラメータＲs、ｒ1、Ｃdを推定する。カーブフィット部６で得られた各パラメータ値は表示部７で表示される。
【００２８】
上記カーブフィット処理は以下の処理（１）〜（７）の手順に従って実行される。
【００２９】
（１）図２に示す、第２の直流電流源５とスイッチ２による矩形状電流を加算する直前の電圧Ｖ0、直後の電圧Ｖ1、測定用電流値Imから等価パラメータＲsの初期値Ｒsiを次式で算出する。ここで電圧Ｖ1採取においては、矩形状電流の加算直後の過渡現象を避けることは前述した通りである。
Ｒsi ＝（Ｖ1−Ｖ0）／Im
【００３０】
（２）次に、矩形状電流の加算後、加算終了前の電圧ＶLから等価パラメータｒ1の初期値ｒ1iを次式で算出する。ここで、ＶLの採取タイミングであるが、矩形状電流の終了間際まで待つ必要はない。図２に示すように、矩形状電流の印加時間の途中で十分である。これは初期値算出の過程は、あくまでも以後の計算の「タネ」となる数値を決定するだけなので厳密さは要求されないからである。
ｒ1i ＝｛( ＶL−Ｖ0 )／Im ｝−Ｒsi
【００３１】
（３）測定した時系列電圧データ内の中程の時刻T(M)と、そのときの電圧Ｖ(M)から等価回路パラメータＣdの初期値を次式で算出する。
すなわち、図８の等価回路において、時刻T(M)のときの両端電圧をＶ(M)、流れる電流をImとすると次式を得る。
Ｖ(M)＝Im・Ｒsi+Im・ｒ1i{１-exp(-T(M)/(Ｃdi・ｒ1i）)}
【００３２】
したがって、Ｃdi は次式で得られる。
Ｃdi ＝−T(M)／[ｒ1i・log{１−((Ｖ(M)/Im)−Ｒsi)／ｒ1i}]
尚、図２ではＶ(M)＜ＶLであるが、これは逆であっても良い。
【００３３】
（４）各パラメータの初期値Ｒsi、ｒ1i、Ｃdiを図８に示す等価回路に適用した場合、供試インピーダンス両端のモデル応答の電圧ｆ(Ｒs,ｒ1,Ｃd)は（３）のＶ(M)の式と同等であるので、次のようになる。
ｆ（Ｒs，ｒ1，Ｃd）＝［Ｒs +ｒ1{1-exp(-T(M)/(Ｃd・ｒ1))} ］・Im
ここで、モデル応答ｆ（Ｒs，ｒ1，Ｃd）と実測値Ｖ(M)との差を取り、これを次式に示す誤差関数erf（ Ｒs，ｒ1，Ｃd ）とする。
【数１】
0003782026.tif 000002
式を整理すると、誤差関数は次のようになる。
【数２】
0003782026.tif 000003
【００３４】
（５）上記で求めた等価パラメータの初期値Ｒsi、ｒ1i、Ｃdiから始めて、この二乗誤差関数erf ( Ｒs,ｒ1、Ｃd )が最小値を示すまで以下の（６）、（７）の処理を繰り返す。
【００３５】
（６）この二乗誤差関数式を各変数Ｒs、ｒ1、Ｃdについて偏微分し、これら変数から成る一次結合（Ｃ1・Ｒs＋Ｃ2・ｒ1＋Ｃ3・Ｃd：ここでＣ1〜Ｃ3は係数）を調べて最も傾きが急な方向を求める。
【００３６】
（７）上記処理（６）で決定された方向に変数Ｒs、ｒ1、Ｃdの値を変化させ、その方向で、二乗誤差が最少となる点を求め、再び処理（６）を行う。
【００３７】
上記の二乗誤差関数erf(Ｒs,ｒ1,Ｃd)が最小値を示したときのＲs、ｒ1、Ｃdが求める供試インピーダンスの各パラメータの推定値である。
【００３８】
図２に示す、時系列電圧データの採取期間であるが、これはスイッチ２をオフしている期間中のすべてのデータを採取する必要はなく、図２に示すように矩形状電流の途中までのデータ採取で良い。
【００３９】
本発明では、１回のデータ収集（スイッチ２による矩形状電流加算は１回）だけで供試インピーダンスのパラメータを推定できる。
【００４０】
尚、本実施形態では、第２の直流電流源５をスイッチ２を介して加算しているが、第１の直流電流源のみを用い、通電電流を与えるとともに矩形状の電流変化を与えても同様の効果が得られることは勿論である。
【００４１】
＜第２の実施形態＞
図３は第２の実施形態に使用する直流電流源５と２極双投スイッチ８であり、他の部分は図１の第１の実施形態と同じである。
【００４２】
第２の実施形態では直流電流源５が２極双投のスイッチ８に接続されており、２極双投スイッチ８をａ側、ｂ側に交互に切り換えることより矩形状電流の極性を正負交互に供試インピーダンス３に加算することができる。その様子を図４に示す。
【００４３】
測定の手順等は第１の実施形態と同様であり、矩形状電流を交互に印加して、このときの供試インピーダンス両端の時系列データを採取する。そしてこのデータを用いて各矩形状電流印加時のパラメータ推定を行い、これらを平均化して最終的に供試インピーダンスのパラメータを推定する。
【００４４】
このように正負交互に測定用電流を加算することにより通電電流の平均値を変化させること無く測定が行える。したがって、この技術を複数回のデータ収集に適用すれば、測定のＳ／Ｎ向上はもちろんのこと、通電電流の平均値を変化させることなく測定を行えるため、測定条件をより厳密に定めることができるというメリットがある。
【００４５】
ここで、第２の実施形態において、第２の直流電流源５及び２極双投スイッチ８を用いずに、第１の直流電流源の電流値を通電電流である基準値に対して増減しても同じ効果が得られることは勿論である。
【００４６】
図５は第３の実施形態を示す図であり、供試インピーダンス３に対して図２又図４に示す測定パルスを繰り返し印加して測定データの平均化を行うように構成されている。こうすることにより、より確度の高いインピーダンスパラメータの推定を可能とした例である。供試インピーダンス３の両端電圧を電圧測定部４で複数回採取し、これらのデータを平均化して測定時のＳＮ比を向上させる。
【００４７】
図５のブロック図を参照しながら本実施形態を説明する。供試インピーダンス３の両端電圧は、Ａ／Ｄ変換回路５０１でデジタル信号に変換される。このデジタル信号は、加算回路５０２及び波形メモリ５０３から成る累積回路に加えられる。加算回路には２つの入力端子があり、一方にはＡ／Ｄ変換回路５０１からのデータが入力され、他方には波形メモリの出力データが入力されている。スケーリング回路５０４は、累積回数に対応して、波形メモリ５０３の出力データをスケーリングする。制御回路５０５は、Ａ／Ｄ変換回路５０１のサンプリング・スタートや波形メモリのタイミングを決める。
【００４８】
具体的な動作を説明すると、先ず測定を開始するとき、すなわち初期状態では波形メモリはクリアされている。この状態で供試インピーダンス３に図２又は図４に示す電流変化を与えて１回目の測定を行い、Ａ／Ｄ変換回路５０１からデータdata1が出力される。すると、このデータdata1は、加算回路５０２の一方の入力に供給され、他方の入力には波形メモリ５０３の出力が供給されている。この波形メモリは、前述の通り、初期状態ではクリアされているので、この入力値はゼロである。従って、加算回路５０２の出力はデータdata1となり、これが波形メモリに蓄積される。
【００４９】
２回目の測定において、Ａ／Ｄ変換回路５０１からデータdata2が出力されたとする。すると、加算回路には、このデータdata２及び波形メモリ５０３からの出力（data1）とが入力されるので、加算回路５０２の出力はデータ（data1＋data2）となる。
【００５０】
以下、３回目の測定ではデータdata3が、４回目の測定ではデータdata4が、というように順次Ｎ回のデータ採取を行うと、波形メモリにはデータ（data1＋data2＋・・・＋data N）が蓄積される。この積算結果は、スケーリング回路５０４に入力され、その出力からは１／Ｎにスケールダウンされたデータ（data1＋data2＋・・・＋data N）／Ｎ、即ち平均値が出力される。そして、この測定データの平均値を用いて第１の実施例と同様に、供試インピーダンス３のパラメータ値を推定する。
【００５１】
以上述べたように平均化を行えば、測定値の分散、即ちバラツキは測定回数をＮとすると１／√Ｎに減少する。これは例えば、測定回数Ｎを１６とするとそのバラツキは１／４になることを意味する。この結果、測定の信頼度が大幅に向上し、供試インピーダンス３のパラメータ推定の確度も向上する。
【００５２】
本発明は第１の実施形態でも述べたように、平均化を行わない場合、わずか１回の測定（電流変化）を行うだけでインピーダンスパラメータ推定が行えるという本質的に高速な測定方法である。このため、平均化のために測定回数が増加したとしても、トータルな測定時間は従来の交流法など対しても十分な優位性を持っている。
【００５３】
また、本実施例では測定データの平均化を行い、その後にパラメータの推定を行う場合について述べたが、これとは違い、測定データを１回採取する毎にパラメータを推定し、これらの推定パラメータ値を平均化しても良い。
【００５４】
尚、本発明は、上述実施形態に限定されず種々の変形が可能である。例えば、上記実施形態では等価回路モデルとして図８の回路を用いたが、異なる回路構成、すなわち抵抗と容量の並列部分が複数存在する場合等でも本発明は適用可能である。また、図８の等価回路において、起電力Ｖ１を含まない場合でも本発明は適用可能である。
【００５５】
また、上記実施形態では、カーブフィットの手法として各測定点の推定等価回路の応答からの誤差の二乗和を最小とする手法を用いたが、他の評価関数によるカーブフィットを用いることもできる。また、等価回路モデルに合わせて、その他任意の数学的回帰分析手法を用いることもできる。
【００５６】
更に、上記実施形態ではカーブフィット手法を用いて等価回路パラメータＲs、ｒ1、Ｃdの値を推定している。しかし、Ｒsのみを求めたい場合は、図６に示すように印加する矩形状電流の幅をＣdとｒ1による時定数の影響が無視できる程度の短時間とし、実施形態１に示すカーブフィット処理（１）のみを行う。こうすることにより、測定時間を電流遮断法と同程度に短縮しながら、測定時の電流変化を極めて少なくできるため、通電電流を一意に定めることができる。
【００５７】
【発明の効果】
以上述べたように、本発明は、通電電流及びそれに重畳した矩形状電流を供給したときの供試インピーダンスの端子に発生する時系列電圧データを採取し、この採取データを用いてカーブフィット処理を行うことにより供試インピーダンスの等価回路パラメータを推定している。本発明によれば、供試インピーダンス中に起電力を有していても任意のパラメータを少なくとも１回の矩形状の電流変化で推定することができるという優れた特長を持っている。
【００５８】
また本発明によれば、測定に際しての矩形状の電流変化は微少であるため、供試インピーダンスが非線形の場合であっても直流通電電流を一意に定めることが可能である。したがって、パラメータ推定の測定条件の明確化が図れるという従来技術には無い優れた効果を持っている。
【００５９】
更に、通電電流に対して矩形状電流を正負交互に加算、すなわち通電電流に対して矩形状に電流を増減することにより、通電電流の平均値に影響を与えること無く測定が行えるという極めて優れた特長を持っている。
【図面の簡単な説明】
【図１】本発明の一実施形態による供試インピーダンスのパラメータ推定回路図である。
【図２】第１の実施形態における矩形状電流を加算する前後の電流と電圧の変化を示す図である。
【図３】本発明の第２の実施形態に使用する電流源と２極双投スイッチを示す図である。
【図４】図３に示す本発明の第２の実施形態による測定結果としての電流と電圧の変化を示す図である。
【図５】本発明の第３の実施形態による供試インピーダンスのパラメータ推定回路図である。
【図６】本発明において、Ｒsのみを求める場合の測定結果である。
【図７】従来の電流遮断法によるインピーダンスパラメータ推定回路例を示す図である。
【図８】供試インピーダンスの等価回路である。
【図９】電流遮断法における供試インピーダンスの両端電圧波形とそこに流れる電流の時間波形図である。
【図１０】交流法による供試インピーダンス推定回路例を示す図である。
【図１１】図９に示す供試インピーダンス推定回路で作成する複素平面インピーダンス図である。
【符号の説明】
１    第１の直流電流源
２    スイッチ
３    供試インピーダンス
４    電圧測定部
５    第２の直流電流源
６    カーブフィット部
７    表示部
８    ２極双投スイッチ
１０１  直流電流源
１０２  スイッチ
１０３  供試インピーダンス
１０４  電圧測定部
２０１  直流電流源
２０２  交流電流計
２０３  抵抗
２０４  交流信号源
２０５  供試インピーダンス
２０６  交流電圧計
５０１  Ａ／Ｄ変換回路
５０２  加算回路
５０３  波形メモリ
５０４  スケーリング回路
５０５  制御回路
【図１】
0003782026.tif 000004
【図２】
0003782026.tif 000005
【図３】
0003782026.tif 000006
【図４】
0003782026.tif 000007
【図５】
0003782026.tif 000008
【図６】
0003782026.tif 000009
【図７】
0003782026.tif 000010
【図８】
0003782026.tif 000011
【図９】
0003782026.tif 000012
【図１０】
0003782026.tif 000013
【図１１】
0003782026.tif 000014