CN101692109B - 基于系统辨识的阻抗测量方法 - Google Patents
基于系统辨识的阻抗测量方法 Download PDFInfo
- Publication number
- CN101692109B CN101692109B CN2009100935893A CN200910093589A CN101692109B CN 101692109 B CN101692109 B CN 101692109B CN 2009100935893 A CN2009100935893 A CN 2009100935893A CN 200910093589 A CN200910093589 A CN 200910093589A CN 101692109 B CN101692109 B CN 101692109B
- Authority
- CN
- China
- Prior art keywords
- impedance
- measured object
- measured
- rpx
- model
- 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
Abstract
本发明公开了一种基于系统辨识的阻抗测量方法。该方法包括步骤:S1,利用阻抗测量仪器测得不同角频率下的阻抗值;S2,根据被测量对象及其与所述阻抗测量仪器探头的连接方式,建立含寄生参量的阻抗模型;S3,根据所测得的阻抗值和所建立的阻抗模型采用参数估计的方法进行系统辨识,得到所述被测试对象的实际阻抗值、电容、电感和各寄生参量的值。本发明的技术方案能够从被测试对象的阻抗测量值中分析出实际阻抗值,以消除寄生参量的影响;而且消除了由被测量对象与测量仪器探头接触的不确定性及各寄生参量的变化对最终测量结果的影响。
Description
技术领域
本发明涉及阻抗测量领域,尤其涉及一种基于系统辨识的阻抗测量方法。
背景技术
阻抗测量在各领域中都有着广泛的需求。在阻抗测量中,一个重要的问题是寄生参量的问题。虽然现有的测量仪器都采取了一定的方法(如屏蔽线法、N探针法等)来减小寄生参量的值,但寄生电阻、寄生电感和寄生电容还是不可避免,这会严重影响实际阻抗的测量。
如上所述,现有测量方法和测量仪器的最大缺陷就在于测得的阻抗值是包含了寄生参量的阻抗值,而不是被测试对象的实际阻抗值。而很多情况下,在测量阻抗、电容或者电感时,寄生参量的存在会造成阻抗的测量值与被测量对象的实际值之间存在很大的差别。而且,由于每次测量时,被测量对象与测量仪器探头的接触情况不确定,因此不同的接触会导致对同一测量对象的多次测量结果差别较大。
发明内容
本发明的目的在于针对现有技术的不足,提供一种能够获得被测量对象的实际阻抗、电容、电感及寄生参量值的方法。
为达到上述目的,本发明提供了一种基于系统辨识的阻抗测量方法,包括以下步骤:
S1,利用阻抗测量仪器测得不同角频率下的阻抗值;
S2,根据被测量对象和与所述阻抗测量仪器的探头的连接方式,建立含寄生参量的阻抗模型;
S3,建立目标函数,并根据所测得的阻抗值和所建立的阻抗模型采用参数估计的方法进行系统辨识,得到所述被测试对象的实际阻抗、电容、电感和各寄生参量的值。
其中,若被测量对象为阻抗,则所建立的阻抗模型的总阻抗Z可以为:
其中,Zx是被测量对象的阻抗,Zxr和Zxi分别是Zx的实部和虚部,Rp是并联寄生电阻,Cp是并联寄生电容,Rs是串联寄生电阻,Ls是串联寄生电感; w为角频率;
若被测量对象为电容,则所建立的阻抗模型的总阻抗Z可以为:
其中,Cx、Rpx和Rsx分别是被测量对象的电容值、并联电阻值和串联电阻值;
若被测对象为电感,则所建立的阻抗模型的总阻抗Z可以为:
其中,Lx、Rpx和Rsx分别是被测量对象的电感值、并联电阻值和串联电阻值。
若被测量对象为阻抗,则所建立的阻抗模型的总阻抗Z还可以为:
其中,Zxr和Zxi分别是被测量对象的阻抗Zx的实部和虚部,Cp是并联寄生电容,Rs是串联寄生电阻, w为角频率;
若被测量对象为电容,则所建立的阻抗模型的总阻抗Z还可以为:
其中,Cx、Rpx和Rsx分别是被测量对象的电容值、并联电阻值和串联电阻值;
若被测对象为电感,则所建立的阻抗模型的总阻抗Z还可以为:
其中,Lx、Rpx和Rsx分别是被测量对象的电感值、并联电阻值和串联电阻值。
在所述步骤S3中所采用的参数估计方法优选为最小二乘法。
若被测量对象为电容,还可以采用直流I-V法与参数估计方法结合进行系统辨识,得到所述被测试对象的实际电容、串并联电阻和各寄生参量的值。
本发明的技术方案通过在不同角频率下测量阻抗值,然后利用系统辨识的方法建立具有寄生参量的阻抗模型,并进行参数估计,从而得到被测量对象的实际阻抗值和其它各参量。本发明的方法能够从被测试对象的阻抗测量值中分析出实际阻抗值,以消除寄生参量的影响;而且消除了由被测量对象与测量仪器探头接触的不确定性及各寄生参量的变化对最终测量结果的影响。
附图说明
图1是本发明实施方式的方法流程图;
图2是利用本发明的方法测量阻抗时所建立的含寄生参量的阻抗模型;
图3是利用本发明的方法测量阻抗时所建立的含寄生参量的简化阻抗模型;
图4是利用本发明的方法测量电容时所建立的含寄生参量的阻抗模型;
图5是利用本发明的方法测量电容时所建立的含寄生参量的简化阻抗模型;
图6是利用本发明的方法测量电感时所建立的含寄生参量的阻抗模型;
图7是利用本发明的方法测量电感时所建立的含寄生参量的简化阻抗模型。
具体实施方式
下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明,但不用来限制本发明的范围。
本发明的核心思想是建立含寄生参量的阻抗模型,利用不同角频率点测得的阻抗值,基于阻抗模型进行系统辨识,从而获得被测量对象的实际阻抗和阻抗模型中各参量(包括电容、电感、串并联电阻、寄生电容和寄生电阻等)。图1是本发明实施方式的方法流程图,如图1所示,具体实施步骤包括:
S1:利用现有的阻抗测量仪器(可以采用例如射频的电流-电压法、网络分析法或自平衡电桥法)测量不同角频率点的阻抗值。
S2:建立含寄生参量的阻抗模型
(1)若被测量对象为阻抗
在各种阻抗测量仪器中,普遍存在着寄生电阻、寄生电容和寄生电感,它们与被测量对象的关系如图2所示。其中,Zx是被测量对象的实际阻抗,Rp是并联寄生电阻,Cp是并联寄生电容,Rs是串联寄生电阻,Ls是串联寄生电感。
由图1可以得到总输出阻抗Z为
其中, w为角频率。
在实际的阻抗测量仪器中,并联寄生电阻Rp很大,串联寄生电感Ls很小,对阻抗值测量的影响可以忽略,从而得到如图3所示的简化模型。将Zx表示成实部Zxr和虚部Zxi,Zx=Zxr+jZxi。
该简化模型的总阻抗Z为
(2)若被测量对象为电容
在实际情况中不存在纯电容,因为电容中常包含有电阻成分。当要测量电容对象的电容值时,所建立的含寄生参量的阻抗模型,如图4所示。图4中的寄生参量与图2相同,虚线框内是被测电容的本征模型,其中,Cx、Rpx和Rsx分别是被测量对象的实际电容值、并联电阻值和串联电阻值。
由图4可以得出总输出阻抗Z为
在实际的阻抗测量仪器中,寄生参量Rp很大,Ls很小,因此这两者对电容值测量的影响可以忽略,从而得到如图5所示的简化模型。
该简化模型的总阻抗Z为
(3)若被测量对象为电感
在实际情况中不存在纯电感,因为电感中亦包含有电阻成分。当要测量电感对象的电感值时,建立如图6所示的含寄生参量的模型。图6中的寄生参量与图2相同,虚线框内是电感的本征模型,其中,Lx、Rpx和Rsx分别是被测量对象的实际电感值、并联电阻值和串联电阻值。
总输出阻抗Z为
在实际的阻抗测量仪器中,寄生参量Rp很大,Ls很小,对电感值测量的影响可以忽略,从而得到如图7所示的简化模型。
该简化模型的总阻抗Z为
S3:参数估计
此步骤为系统辨识中的参数估计步骤。参数估计的方法很多,可以采用例如,最小二乘法、极大似然法等,本实施例采用最小二乘法来估计模型中的各个参量值。
由于阻抗值Z是一个复数,在进行最小二乘法辨识时,需将Z转换成如幅值-相角形式或实部-虚部形式,本实施例选择实部-虚部形式,从而得到系统辨识的目标函数为
其中,Zrdata,Zidata分别是所测得的阻抗值的实部和虚部,Zr,Zi分别是阻抗模型的输出阻抗Z的实部和虚部,N为在不同角频率下进行测量的测量点的总个数。
采用最小二乘法进行参数估计时,根据阻抗模型选取各参数的值,使得阻抗模型的输出阻抗Z满足(7)式,即残差平方和为最小。其中,阻抗模型的选择,在本实施例中,是根据测量对象分别为阻抗、电容、电感,分别选择图3、5、7所示模型,模型输出的阻抗值Z分别利用上述(2)、(4)、(6)式计算得到。参数估计的结果是被测试对象的实际阻抗值、电容、电感和各寄生参量的值。例如,当被测量对象为阻抗时,得到实际阻抗值Zx以及各寄生参量(并联寄生电容Cp,串联寄生电阻Rs)的值。
如果被测量对象是电容,仅仅采用上述参数估计方法进行系统辨识,所获取的并联电阻Rpx的值在很多情况下不准确。因此,可结合直流I-V法,较精确地估计出并联电阻Rpx的值。
直流I-V法的方法为:对被测量对象施加一直流电压U0(例如在图4中,是在并联寄生电容Cp的左端点与串联寄生电感Ls的右端点之间施加该直流电压U0),测得经过被测量对象的电流I0,得到直流电阻值
则可以得到Rpx的估计值如式(8)所示,其作为Rpx的系统辨识结果。
Rpx0=R0-Rsx0-Rs0 (8)
其中,Rsx0,Rs0分别是对Rsx,Rs利用上述的最小二乘法进行参数估计的系统辨识结果。
由以上实施方案可以看出,本发明通过在不同角频率下测量阻抗值,然后利用系统辨识的方法建立具有寄生参量的阻抗模型,并进行参数估计,从而得到被测量对象的实际阻抗值。本发明的方法能够从被测试对象的阻抗测量值中分析出实际阻抗值,以消除寄生参量的影响;而且消除了由被测量对象与测量仪器探头接触的不确定性及各寄生参量的变化对最终测量结果的影响。
以上实施方式仅用于说明本发明,而并非对本发明的限制,有关技术领域的普通技术人员,在不脱离本发明的精神和范围的情况下,还可以做出各种变化和变型,因此所有等同的技术方案也属于本发明的范畴,本发明的专利保护范围应由权利要求限定。
Claims (4)
1.一种基于系统辨识的阻抗测量方法,包括以下步骤:
S1,利用阻抗测量仪器测得不同角频率下的阻抗值;
S2,根据被测量对象及其与所述阻抗测量仪器探头的连接方式,建立含寄生参量的阻抗模型;
S3,根据所测得的阻抗值和所建立的阻抗模型采用参数估计方法进行系统辨识,得到所述被测试对象的实际阻抗、电容、电感和各寄生参量的值;
所述含寄生参量的阻抗模型包括:
若被测量对象为阻抗,则所建立的阻抗模型的总阻抗Z为:
若被测量对象为电容,则所建立的阻抗模型的总阻抗Z为:
其中,Cx、Rpx和Rsx分别是被测量对象的电容值、并联电阻值和串联电阻值;
若被测对象为电感,则所建立的阻抗模型的总阻抗Z为:
其中,Lx、Rpx和Rsx分别是被测量对象的电感值、并联电阻值和串联电阻值。
2.一种基于系统辨识的阻抗测量方法,包括以下步骤:
S1,利用阻抗测量仪器测得不同角频率下的阻抗值;
S2,根据被测量对象及其与所述阻抗测量仪器探头的连接方式,建立含寄生参量的阻抗模型;
S3,根据所测得的阻抗值和所建立的阻抗模型采用参数估计方法进行系统辨识,得到所述被测试对象的实际阻抗、电容、电感和各寄生参量的值;
所述含寄生参量的阻抗模型包括:
若被测量对象为阻抗,则所建立的阻抗模型的总阻抗Z为:
若被测量对象为电容,则所建立的阻抗模型的总阻抗Z为:
其中,Cx、Rpx和Rsx分别是被测量对象的电容值、并联电阻值和串联电阻值;
若被测对象为电感,则所建立的阻抗模型的总阻抗Z为:
其中,Lx、Rpx和Rsx分别是被测量对象的电感值、并联电阻值和串联电阻值。
3.如权利要求1或2所述的基于系统辨识的阻抗测量方法,其特征在于,在所述步骤S3中所采用的参数估计方法为最小二乘法。
4.如权利要求1或2所述的基于系统辨识的阻抗测量方法,其特征在于,所述步骤S3中的系统辨识,若被测量对象为电容,还采用直流I-V法与参数估计方法相结合进行系统辨识,得到所述被测试对象的实际电容、串并联电阻和各寄生参量的值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100935893A CN101692109B (zh) | 2009-10-13 | 2009-10-13 | 基于系统辨识的阻抗测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100935893A CN101692109B (zh) | 2009-10-13 | 2009-10-13 | 基于系统辨识的阻抗测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101692109A CN101692109A (zh) | 2010-04-07 |
CN101692109B true CN101692109B (zh) | 2011-08-17 |
Family
ID=42080817
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009100935893A Expired - Fee Related CN101692109B (zh) | 2009-10-13 | 2009-10-13 | 基于系统辨识的阻抗测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101692109B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5589988B2 (ja) * | 2011-07-28 | 2014-09-17 | 横河電機株式会社 | 電池インピーダンス測定装置 |
US9428071B2 (en) * | 2014-01-14 | 2016-08-30 | Ford Global Technologies, Llc | Impedance based battery parameter estimation |
CN104049149B (zh) * | 2014-06-20 | 2017-06-23 | 国家电网公司 | 基于二元线性回归模型的配电网双端线路段阻抗分析方法 |
CN105653848A (zh) * | 2015-12-25 | 2016-06-08 | 大连理工大学 | 一种利用单端口阻抗和灵敏度测量耦合线圈参数的方法 |
CN105717369B (zh) * | 2016-02-16 | 2019-07-12 | 常州同惠电子股份有限公司 | 一种用电容表测量电容的串联等效电感的方法 |
CN106532685B (zh) * | 2016-10-26 | 2019-02-01 | 浙江大学 | 用于并网逆变器稳定分析的广义阻抗判据计算方法及应用 |
CN106992765A (zh) * | 2017-04-18 | 2017-07-28 | 河北工业大学 | 谐波电流激励下降低感性电路阻抗值的方法 |
CN108470104A (zh) * | 2018-03-22 | 2018-08-31 | 常州同惠电子股份有限公司 | 用于阻抗分析仪的等效电路参数建模计算方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6118275A (en) * | 1998-06-25 | 2000-09-12 | Korea Kumho Petrochemical Co., Ltd. | Method and apparatus for measuring battery capacity using voltage response signal based on pulse current |
US6208147B1 (en) * | 1998-06-16 | 2001-03-27 | Korea Kumho Petrochenical Co., Ltd. | Method of and apparatus for measuring battery capacity by impedance spectrum analysis |
CN1543912A (zh) * | 2003-11-18 | 2004-11-10 | 华中科技大学 | 一种生物组织多频率阻抗测量方法及其装置 |
CN1761366A (zh) * | 2005-08-18 | 2006-04-19 | 南京大学 | 基于系统辨识的扬声器单元非线性参数测量方法 |
-
2009
- 2009-10-13 CN CN2009100935893A patent/CN101692109B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6208147B1 (en) * | 1998-06-16 | 2001-03-27 | Korea Kumho Petrochenical Co., Ltd. | Method of and apparatus for measuring battery capacity by impedance spectrum analysis |
US6118275A (en) * | 1998-06-25 | 2000-09-12 | Korea Kumho Petrochemical Co., Ltd. | Method and apparatus for measuring battery capacity using voltage response signal based on pulse current |
CN1543912A (zh) * | 2003-11-18 | 2004-11-10 | 华中科技大学 | 一种生物组织多频率阻抗测量方法及其装置 |
CN1761366A (zh) * | 2005-08-18 | 2006-04-19 | 南京大学 | 基于系统辨识的扬声器单元非线性参数测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101692109A (zh) | 2010-04-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101692109B (zh) | 基于系统辨识的阻抗测量方法 | |
CN101907437B (zh) | 一种基于小波差分算法的电缆故障测距方法 | |
Cataliotti et al. | Current transformers effects on the measurement of harmonic active power in LV and MV networks | |
CN101799494B (zh) | 基于盲信号处理技术的传导电磁干扰噪声软分离方法 | |
RU2009135814A (ru) | Измерение полного сопротивления линии электропередачи | |
CN103344828A (zh) | 一种新的谐波发射水平评估方法 | |
CN105353226A (zh) | 基于散射参数和智能算法的emi噪声源阻抗等效参数提取方法 | |
CN105223466B (zh) | 一种利用模极大值比的特高压直流输电线路单端测距方法 | |
CN105997072A (zh) | 一种基于测量精度的电阻抗检测用电极接触状态检测方法 | |
CN105974264A (zh) | 一种基于相电流暂态特征的故障选线方法 | |
CN106338237A (zh) | 基于频响阻抗法的变压器绕组变形检测方法 | |
US20090088995A1 (en) | Method for determining the linear electrical response of a transformer, generator or electrical motor | |
CN102590598B (zh) | 基于多阈值比较的超声波信号过零点预测方法 | |
CN109239533A (zh) | 一种基于人工神经网络的特高压直流输电线路的故障定位方法 | |
CN102129042B (zh) | 电池内阻检测的方法及装置 | |
CN102043084A (zh) | 一种检测避雷器阻性电流的检测方法 | |
CN103529350A (zh) | 一种大功率变压器匝间短路故障的检测方法及装置 | |
CN101153899A (zh) | 数字多用表自动测量系统及方法 | |
CN1996030A (zh) | 测定传导性电磁干扰噪声源内阻抗的装置 | |
CN104459335A (zh) | 一种用于配电网对地电容的检测装置及其检测方法 | |
CN101210943A (zh) | 测量电力直流系统对地绝缘电阻的方法 | |
CN102012462A (zh) | 用于检测电源电路火线-地线回路阻抗的装置及方法 | |
CN102608425B (zh) | 一种宽频阻抗测量系统及宽频阻抗测量方法 | |
CN103913654A (zh) | 高压直流输电系统接地极故障特征的识别方法以及设备 | |
CN104678339B (zh) | 一种用于探针式微波电压测量系统的校准装置、系统及方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110817 Termination date: 20171013 |