[go: up one dir, main page]

CN111539171B - 一种不依赖于初始点的时间常数估算方法 - Google Patents

一种不依赖于初始点的时间常数估算方法 Download PDF

Info

Publication number
CN111539171B
CN111539171B CN202010316181.4A CN202010316181A CN111539171B CN 111539171 B CN111539171 B CN 111539171B CN 202010316181 A CN202010316181 A CN 202010316181A CN 111539171 B CN111539171 B CN 111539171B
Authority
CN
China
Prior art keywords
formula
initial point
time constant
estimation method
following formula
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.)
Active
Application number
CN202010316181.4A
Other languages
English (en)
Other versions
CN111539171A (zh
Inventor
徐伟专
甘伟
刘国福
胡能钢
田武刚
文诚
范亚
熊艳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Yinhe Electric Co ltd
Original Assignee
Yinhe Electric Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Yinhe Electric Co ltd filed Critical Yinhe Electric Co ltd
Priority to CN202010316181.4A priority Critical patent/CN111539171B/zh
Publication of CN111539171A publication Critical patent/CN111539171A/zh
Application granted granted Critical
Publication of CN111539171B publication Critical patent/CN111539171B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/30Circuit design
    • G06F30/36Circuit design at the analogue level
    • G06F30/367Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of Resistance Or Impedance (AREA)

Abstract

本发明涉及测量系统中时间常数的估算方法,具体说是一种不依赖于初始点的时间常数估算方法,包括由等效电阻Rx和等效电容Cx组成的被测对象等效电路和该等效电路的测量电路,所述测量电路由理想恒压源Vin(t)、继电器KV、保护电阻R0和测量电阻Rm组成,并依据KCL和KVL定律和式(1)‑(8)得到τ、Rx、Cx和t0值。本发明的方法测试时间大幅缩短,测试系统效率明显提高,不依赖时间初始点,测试方案可行性强,并通过滤波消除了噪声和纹波干扰,提高了测量精度。且本发明方法不仅适用于测量一阶系统,也适用于仅有一个闭环主导极点的高阶系统。

Description

一种不依赖于初始点的时间常数估算方法
技术领域
本发明涉及测量系统中时间常数的估算方法,具体说是一种不依赖于初始点的时间常数估算方法。
背景技术
在测量一阶系统中,不可避免地会引入噪声和误差,若系统的时间常数较大,测量系统需要花费几分钟甚至十几分钟才能达到稳定,这显然无法满足工程实测需求。例如,在电机工作间隙测量绝缘电阻时,由于电机的绕组之间、绕组与地之间的等效电容值较大导致时间常数较大,若待测量系统达到稳态后求取时间常数,则会导致测量时间较长,降低测量系统和电机的工作效率。针对该现象,有学者提出了根据一阶电路的暂态阶段数学模型估算时间常数的方法,该方法在明确时间初始点的情况下能够准确估算得到时间常数,由于在复杂工况下一阶电路的时间初始点较难精准定位,因此该方法的可行性不强,无法满足工程需求。
发明内容
针对上述技术问题,本发明提供一种不依赖于初始点的时间常数估算方法,该方法通过建立测量系统和被测对象的数学模型,然后取模型中的任意时间间隔的三点,计算得到时间常数和时间初始点。该方法能够大幅缩短测量时间,显著提高测量系统的效率,而且不依赖于时间初始点,具备很强的可实施性。
本发明所采用的技术方案为:一种不依赖于初始点的时间常数估算方法,包括由等效电阻Rx和等效电容Cx组成的被测对象等效电路和该等效电路的测量电路,所述测量电路由理想恒压源Vin(t)、继电器KV、保护电阻R0和测量电阻Rm组成,并依据KCL和KVL定律得到以下式(1):
Figure BDA0002459672790000021
上述式(1)中,VRm(t)为测量电压,τ为时间常数,Rall=R0+Rm,t0为继电器KV完全闭合后的时间初始点,由式(1)得到以下式(2):
Figure BDA0002459672790000022
上述式(2)中,n(t)为测量电路产生的噪声,对K(t)进行滤波并引入常量A后,K(t)的表达式为以下式(3):
Figure BDA0002459672790000023
根据式(3)并引入变量Kf(t)得到以下式(4):
Figure BDA0002459672790000024
上述式(4)中t3-t2=t2-t1=Δt,由此解出τ。
进一步地,将上述解出的τ值代入上述式(3)得到以下式(5):
Figure BDA0002459672790000031
根据上述式(5)中的任意两个等式,得到以下式(6):
Figure BDA0002459672790000032
由上述式(6)得到Rx和A。
进一步地,由上述式(6)和式(1)得到以下式(7):
Figure BDA0002459672790000033
由上4.述式(7)得到Cx
进一步的,由上述式(6)和式(3)得到以下式(8):
Figure BDA0002459672790000034
由上述式(8)得到t0
由以上技术方案可知,本发明的方法测试时间大幅缩短,测试系统效率明显提高,不依赖时间初始点,测试方案可行性强,并通过滤波消除了噪声和纹波干扰,提高了测量精度。且本发明方法不仅适用于测量一阶系统,也适用于仅有一个闭环主导极点的高阶系统。
附图说明
图1是本发明的测量电路和被测对象等效电路的组成电路图。
具体实施方式
下面结合图1对本发明进一步详细说明:
本发明提供一种不依赖于初始点的时间常数估算方法,如图1中虚线框内为被测对象等效电路,它由等效电阻Rx和等效电容Cx组成。虚线框外为测量电路,它由理想恒压源Vin(t)、继电器KV、保护电阻R0和测量电阻Rm组成,R0通过对测量电压VRm(t)的限幅实现对后级测量电路的保护。对图1所示电路采用KCL和KVL定律可得以下式(1):
Figure BDA0002459672790000041
式(1)中,τ为时间常数Rall=R0+Rm,t0为继电器KV完全闭合后的开始测量点,也即时间初始点。由式(1)可得以下式(2):
Figure BDA0002459672790000042
式(2)中,n(t)为测量电路产生的噪声,为避免该噪声影响测量精度,需对K(t)进行滤波,滤波并引入常量A后,K(t)的表达式为以下式(3):
Figure BDA0002459672790000043
根据式(3)并引入变量Kf(t)可得以下式(4):
Figure BDA0002459672790000051
式(4)中t3-t2=t2-t1=Δt,由此解出τ。
若需要进一步求取被测对象的等效电阻、等效电容以及时间初始点,可按以下步骤进行:
将τ代入式(3)可得式(5):
Figure BDA0002459672790000052
根据上述式(5)中任意两个等式(建议选择t1和t3相关等式),按以下式(6)求解得到Rx和A:
Figure BDA0002459672790000053
由式(6)和(1)得到式(7),进一步求得Cx,;由式(6)和(3)得到式(8),进一步求得t0
Figure BDA0002459672790000061
Figure BDA0002459672790000062
为验证该上述估算方法的有效性,将测试系统和被测对象等效为理想元器件的组合,各个器件的取值为:R0=2MΩ,Rm=2kΩ,t0=5.5s,Rx=100MΩ,Cx=0.5uF。然后通过式(3)求得在各个时间点处的VRm。通过采用该估算方法求得的结果如下表。需指出,表中第n组待求量是通过第(n-2)、(n-1)、n组已知量求得。例如,第3组待求量是通过第1、2、3组已知量求得。表中显示,无论通过哪三组已知量所求得的τ、Rx、Cx和t0值都与预设值一致,由此验证了该估算方法的有效性。
Figure BDA0002459672790000063
上述实施方式仅供说明发明之用,而并非是对本发明的限制,所以有关技术领域的普通技术人员,在不脱离本发明精神和范围的情况下,还可以作出各种变化和变型,因此所有等同的技术方案也应属于本发明的范畴。

Claims (4)

1.一种不依赖于初始点的时间常数估算方法,包括由等效电阻Rx和等效电容Cx组成的被测对象等效电路和该等效电路的测量电路,其特征在于:所述测量电路由理想恒压源Vin(t)、继电器KV、保护电阻R0和测量电阻Rm组成,并依据KCL和KVL定律得到以下式(1):
Figure FDA0002459672780000011
上述式(1)中,VRm(t)为测量电压,τ为时间常数,Rall=R0+Rm,t0为继电器KV完全闭合后的时间初始点,由式(1)得到以下式(2):
Figure FDA0002459672780000012
上述式(2)中,n(t)为测量电路产生的噪声,对K(t)进行滤波并引入常量A后,K(t)的表达式为以下式(3):
Figure FDA0002459672780000013
根据式(3)并引入变量Kf(t)得到以下式(4):
Figure FDA0002459672780000014
上述式(4)中t3-t2=t2-t1=Δt,由此解出τ。
2.如权利要求1所述不依赖于初始点的时间常数估算方法,其特征在于:将上述解出的τ值代入上述式(3)得到以下式(5):
Figure FDA0002459672780000021
根据上述式(5)中的任意两个等式,得到以下式(6):
Figure FDA0002459672780000022
由上述式(6)得到Rx和A。
3.如权利要求2所述不依赖于初始点的时间常数估算方法,其特征在于:由上述式(6)和式(1)得到以下式(7):
Figure FDA0002459672780000023
由上述式(7)得到Cx
4.如权利要求2或3所述不依赖于初始点的时间常数估算方法,其特征在于:由上述式(6)和式(3)得到以下式(8):
Figure FDA0002459672780000031
由上述式(8)得到t0
CN202010316181.4A 2020-04-21 2020-04-21 一种不依赖于初始点的时间常数估算方法 Active CN111539171B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010316181.4A CN111539171B (zh) 2020-04-21 2020-04-21 一种不依赖于初始点的时间常数估算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010316181.4A CN111539171B (zh) 2020-04-21 2020-04-21 一种不依赖于初始点的时间常数估算方法

Publications (2)

Publication Number Publication Date
CN111539171A CN111539171A (zh) 2020-08-14
CN111539171B true CN111539171B (zh) 2023-05-12

Family

ID=71953519

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010316181.4A Active CN111539171B (zh) 2020-04-21 2020-04-21 一种不依赖于初始点的时间常数估算方法

Country Status (1)

Country Link
CN (1) CN111539171B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4952808A (en) * 1988-06-07 1990-08-28 U.S. Philips Corp. Thermal radiation detection apparatus
CN1543912A (zh) * 2003-11-18 2004-11-10 华中科技大学 一种生物组织多频率阻抗测量方法及其装置
CN1551501A (zh) * 2003-05-12 2004-12-01 因芬尼昂技术股份公司 校准电阻/电容滤波器电路之装置及方法
CN104198813A (zh) * 2014-05-09 2014-12-10 杭州电子科技大学 一种正交相关法测量超声换能器阻抗角的装置及方法
CN107656140A (zh) * 2017-08-11 2018-02-02 中国科学院地质与地球物理研究所 一种利用嵌入式处理器数字端口测量电阻的方法及电路
CN110096780A (zh) * 2019-04-23 2019-08-06 西安交通大学 一种超级电容一阶rc网络等效电路及参数确定方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4952808A (en) * 1988-06-07 1990-08-28 U.S. Philips Corp. Thermal radiation detection apparatus
CN1551501A (zh) * 2003-05-12 2004-12-01 因芬尼昂技术股份公司 校准电阻/电容滤波器电路之装置及方法
CN1543912A (zh) * 2003-11-18 2004-11-10 华中科技大学 一种生物组织多频率阻抗测量方法及其装置
CN104198813A (zh) * 2014-05-09 2014-12-10 杭州电子科技大学 一种正交相关法测量超声换能器阻抗角的装置及方法
CN107656140A (zh) * 2017-08-11 2018-02-02 中国科学院地质与地球物理研究所 一种利用嵌入式处理器数字端口测量电阻的方法及电路
CN110096780A (zh) * 2019-04-23 2019-08-06 西安交通大学 一种超级电容一阶rc网络等效电路及参数确定方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Anjali Chauhan ; Bahgat Sammakia ; Furat Afram ; Kanad Ghose ; Gamal Refai-Ahmed ; Dereje Agonafer.Analysing real time power of the microprocessor to estimate thermal time constant of the hotspot.13th InterSociety Conference on Thermal and Thermomechanical Phenomena in Electronic Systems.2012,全文. *
王德诚, 李军伟, 高松,孙海波,王鹏.基于转子时间常数在线辨识的车用异步电机转速估计.广西大学学报.第43卷(第2期),498-507. *

Also Published As

Publication number Publication date
CN111539171A (zh) 2020-08-14

Similar Documents

Publication Publication Date Title
WO2015073498A1 (en) Systems and methods for improved accuracy
CN106777545B (zh) 一种电阻闪烁噪声模型的建立方法及系统
CN110888070A (zh) 一种电池温度的估计方法、装置、设备及介质
CN113325352B (zh) 一种针对暂态电流测量用传感器的校准方法及系统
CN107782938A (zh) 电能表自热的误差补偿方法
CN111539171B (zh) 一种不依赖于初始点的时间常数估算方法
CN101308176A (zh) 一种在线测量电解电容的电路
CN112649652B (zh) 一种用于冲击电流量值的溯源方法及系统
Ruba et al. Detailed design of second order model of lithium-ion battery simulator based on experimental measurements
CN113900032B (zh) 利用瞬态响应校正放电深度的方法及系统
Petit Comparison of PD amplitudes of stator bars taken with different instruments
JPH08146057A (ja) 絶縁抵抗の測定方法
CN113361874B (zh) 一种汽车零部件eol测试评价方法
CN107342599A (zh) 一种机电暂态仿真中控制系统稳定性自动诊断和参数调整方法
CN108023576B (zh) 一种用于快沿脉冲发生器上升时间校准的方法
CN105811830A (zh) 基于数据统计和数值优化的永磁同步电机模型建模方法
CN114417583B (zh) 风电机组尾流速度损失的评估方法、装置及存储介质
US10184967B2 (en) Method of determining capacitance value of capacitor while taking applied alternating voltage into consideration, and program
CN110618301B (zh) 一种冲击电压波形调节方法及装置
CN116468000A (zh) 一种pcb设备的关键部位蠕变行为影响的预测方法、装置
Coenen et al. Statistical evaluation of manufacturing tolerances in electrical machines by simulation and measurement
DE102004045151A1 (de) Verfahren und Vorrichtung zur Bestimmung eines Brennraumdruckes in einem Verbrennungsmotor sowie zugehöriges Motorsteuerverfahren und zugehörige Motorsteuerung
CN110657833A (zh) 一种用于高精度源表一体化测量设备的新型校准方法
CN111130542A (zh) 一种模拟信号温度补偿方法
Mihai Metrological triangulation rules in ratio measurements of high standard resistance bridges

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant