CN112269065B - 一种双激励的高精度绝对阻抗测量方法 - Google Patents

一种双激励的高精度绝对阻抗测量方法 Download PDF

Info

Publication number
CN112269065B
CN112269065B CN202011102344.5A CN202011102344A CN112269065B CN 112269065 B CN112269065 B CN 112269065B CN 202011102344 A CN202011102344 A CN 202011102344A CN 112269065 B CN112269065 B CN 112269065B
Authority
CN
China
Prior art keywords
amplitude
excitation
asa
excitation signals
measured
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
CN202011102344.5A
Other languages
English (en)
Other versions
CN112269065A (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.)
Hubei University of Technology
Original Assignee
Hubei University of Technology
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 Hubei University of Technology filed Critical Hubei University of Technology
Priority to CN202011102344.5A priority Critical patent/CN112269065B/zh
Publication of CN112269065A publication Critical patent/CN112269065A/zh
Application granted granted Critical
Publication of CN112269065B publication Critical patent/CN112269065B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R27/00Arrangements for measuring resistance, reactance, impedance, or electric characteristics derived therefrom
    • G01R27/02Measuring real or complex resistance, reactance, impedance, or other two-pole characteristics derived therefrom, e.g. time constant
    • G01R27/14Measuring resistance by measuring current or voltage obtained from a reference source

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of Resistance Or Impedance (AREA)

Abstract

本发明公开了一种双激励的高精度绝对阻抗测量方法,主要包括硬件电路设计与软件算法设计两部分。硬件部分包括含有至少2个通道数模转换(DAC)与3个通道模数转换(ADC)模块的单片机系统、2路DAC输出跟随电路、被测阻抗与标准电阻串联环电路及环电路中间结点传感信号的两级放大电路。2路DAC跟随电路的输出加载到串联环电路的两端,且1路DAC跟随电路与两级传感放大信号连接到3个ADC转换通道进行AD转换。软件部分包括设计DDS方法,通过DAC的输出2路正弦激励信号,设计ADC传感采集信号正弦拟合算法,并根据拟合结果设计调整激励信号的幅值,避免两级放大信号的饱和。软件部分还包括根据传感信号大小与激励信号的幅值,结合标准电阻进一步计算绝对阻抗的大小。

Description

一种双激励的高精度绝对阻抗测量方法
技术领域
本发明涉及测量电路技术与信号处理技术领域,尤其涉及一种双激励的高精度绝对阻抗测量方法。
背景技术
阻抗的测量被广泛应用于工业领域,如用于温度测量的热敏电阻、铂电阻,用于位移测量的各种绕线和非绕线的可变电阻器,以及其他如压敏、湿敏、光敏、气敏电阻等的测量。它们所涉及到的物理量在很多场合要求测量范围大、测量精度高,相应地对被测阻抗的测量范围与测量精度也提出了很高的要求。
惠更斯电桥测量方法可以通过加大差分放大倍数,满足测量精度高的要求,但这种情况下,一旦被测阻抗范围变化较大,电桥差分输出信号会很快饱和,导致测量结果失去意义。因此,这种测量方案只适合阻抗变化很小的相对测量,如材料应力测量中应变片电阻的相对测量等。
一般的电阻分压法,虽然可以实现电阻的大范围绝对测量,但由于这种方法本身输出信号较大,放大倍数也不可能很大,所以测量精度的大小取决于信号采样的分辨率,一般为模数转换(ADC)的分辨率。理论上ADC的转换位数越多,分辨率越高,但实际上由于电源及其他干扰噪声的影响,单纯提高ADC的转换分辨率并不能降低电路的信噪比,也就不能带来测量精度的真正提高。同时,ADC转换位数的提高必然带来测量速度的减慢与成本的提高。
这些矛盾决定了大范围与高精度的阻抗测量在传统测量方案中难以实现。
发明内容
为解决以上问题本发明提出一种双激励的高精度绝对阻抗测量电路与测量方法,在满足测量精度不降低的前提下、提高阻抗测量范围的方法。
本发明的技术方案如下:一种双激励的高精度绝对阻抗测量方法,包括由单片机内部的DAC模块输出幅值不同、相位相差180°的正弦激励信号Asinωt、-Bsinωt,使激励信号的幅值A>0、B>0,A+B=Const,Const为常数,再通过运放的同向跟随后,加载到由被测电阻Rx与标准电阻R串联电路的两端,被测电阻Rx与标准电阻R的中间结点输出电压通过同向一级放大k1倍,输出电压Amsinωt,再通过反向二级放大k2倍,输出电压AOsinωt,正激励信号Asinωt及输出电压Amsinωt和AOsinωt分别送入单片机内部进行AD转换及信号处理,得到被测阻抗Rx的大小,其中,在测量被测阻抗Rx的大小前,需根据以下步骤调整激励信号A和B的幅值:
步骤1.根据所述被测电阻Rx的原始标定值为R0计算激励信号A和B的理想幅值,
Figure GDA0003785636900000011
B=Const-A,Rx=R0
步骤2.根据步骤1得到的激励信号A和B的理想幅值,计算激励信号A和B幅值的上限:A和B,下限:A和B
Figure GDA0003785636900000012
B=Const-A,B=Const-A,所述单片机采用参考电压进行ADC转换,参考电压为2Asa,Asa为假设饱和幅值,k=k1*k2<0;
步骤3.计算当被测电阻Rx相对原始标定值为R0降低时,激励信号A和B的幅值应调整为AN=A,BN=B,对应的Rx的测量量程为
Figure GDA0003785636900000021
计算当被测电阻Rx相对原始标定值为R0升高时,激励信号A和B的幅值应调整为AN=A,BN=B,对应的Rx的测量量程为
Figure GDA0003785636900000022
步骤4.根据被测电阻Rx相对原始标定值的变化方向,按步骤3对应调整激励信号A和B的幅值生成正弦激励信号ANsinωt、-BNsinωt并加载到由被测电阻Rx与标准电阻R串联电路的两端,经同向一级放大和反向二级放大后输出电压的幅值为AO,单片机计算得到被测电阻Rx为
Figure GDA0003785636900000023
进一步地,若步骤4经同向一级放大和反向二级放大后输出电压的幅值AO>Asa,则单片机根据同向一级放大后输出点电压的幅值Am和正弦激励信号Asinωt调整前的幅值A,按下式计算激励信号A和B的幅值调整值ANx、BNx
Figure GDA0003785636900000024
Figure GDA0003785636900000025
BNx=Const-ANx
调整激励信号A和B的幅值为ANx、BNx,生成新的相应正弦激励信号后,加载到被测电阻Rx与标准电阻R串联电路的两端,得到新的电压AO,按照步骤4单片机计算得到被测电阻Rx为
Figure GDA0003785636900000026
优选地,步骤4之后还包括:步骤5.若步骤4经同向一级放大和反向二级放大后输出电压的幅值AO>Asa,当被测电阻Rx相对原始标定值为R0升高时,根据步骤3计算得到Rx的测量量程上限
Figure GDA0003785636900000027
重新计算激励信号A和B的幅值调整值AN=AN+1、BN=BN+1
Figure GDA0003785636900000028
BN+1=Const-AN+1,对应的Rx的测量量程为
Figure GDA0003785636900000029
重复步骤4,若单片机采集到的反向二级放大后输出电压幅值AO不大于Asa,则单片机计算得到被测电阻Rx为
Figure GDA00037856369000000210
步骤6.若经步骤5单片机采集到的反向二级放大后输出电压幅值AO大于Asa,则使BN=BN+1,AN=AN+1,重复步骤5计算得到
Figure GDA00037856369000000211
和激励信号A和B的幅值调整值
Figure GDA00037856369000000212
BN=BN+2=Const-AN+2,单片机继续比较AO,若AO不大于Asa,则单片机计算得到被测电阳
Figure GDA0003785636900000031
Figure GDA0003785636900000032
步骤7.若AO>Asa,则重复步骤6,直至AO不大于Asa,单片机计算得到被测电阻
Figure GDA0003785636900000033
优选地,步骤4之后还包括:步骤5.若步骤4经同向一级放大和反向二级放大后输出电压的幅值AO>Asa,当被测电阻Rx相对原始标定值为R0降低时,根据步骤3计算得到Rx的测量量程下限
Figure GDA0003785636900000034
重新计算激励信号A和B的幅值调整值AN=AN+1、BN=BN+1
Figure GDA0003785636900000035
BN+1=Const-AN+1,对应的Rx的测量量程为
Figure GDA0003785636900000036
重复步骤4,若单片机采集到的反向二级放大后输出电压幅值AO不大于Asa,则单片机计算得到被测电阻Rx为
Figure GDA0003785636900000037
步骤6.若经步骤5单片机采集到的反向二级放大后输出电压幅值Ao大于Asa,则使BN=BN+1,AN=AN+1,重复步骤5计算得到
Figure GDA0003785636900000038
和激励信号A和B的幅值调整值
Figure GDA0003785636900000039
BN=BN+2=Const-AN+2,单片机继续比较AO,若AO不大于Asa,则单片机计算得到被测电阻
Figure GDA00037856369000000310
Figure GDA00037856369000000311
步骤7.若AO>Asa,则重复步骤6,直至AO不大于Asa,单片机计算得到被测电阻
Figure GDA00037856369000000312
进一步地,同向一级放大时,放大倍数k1取1。
本发明的有益效果为:本发明在满足测量精度不降低的前提下、提高阻抗测量范围
附图说明
图1为本发明的一种双激励的高精度绝对阻抗测量电路与测量方法的测量实施方案框图;
图2为本发明的一种双激励的高精度绝对阻抗测量电路与测量方法的测量电路原理图;
图3为本发明的一种双激励的高精度绝对阻抗测量电路与测量方法的单片机主要引脚连接图;
图4为本发明的一种双激励的高精度绝对阻抗测量电路与测量方法的测量与放大电路图;
图5为本发明的一种双激励的高精度绝对阻抗测量电路与测量方法的测量电路仿真波形图;
图6为本发明的一种双激励的高精度绝对阻抗测量电路与测量方法的DMA传送生成正弦模拟量过程图;
图7为本发明的一种双激励的高精度绝对阻抗测量电路与测量方法的激励信号的最佳值与上下限分布图。
具体实施方式
下面结合具体实施例对本发明作进一步详述。虽然以下描述了本发明的优选实施方式,然而应该理解,可以以各种形式实现本发明而不应被这里阐述的实施方式所限制,可以对这些实施方式和实施例进行各种变化和更改,也就是说,这些变化和更改也属于本发明的构思范围。
实施例1
一种双激励的高精度绝对阻抗测量测量方法,如图2所示的测量原理图,由单片机内部的DAC模块输出幅值不同,相位相差180°的正弦激励信号Asinωt,-Bsinωt,使激励信号的幅值A>0、B>0,A+B=Const,Const为常数,加载到由被测电阻Rx与标准电阻R串联电路的两端。中间结点C电压VC(t)满足如下的公式:
Figure GDA0003785636900000041
中间结点C电压经过一级放大k1后,输出电压Amsinωt,输出电压Vm(t)=Amsinωt经过二级放大k2后,输出电压Vo(t)=AOsinωt:
Figure GDA0003785636900000042
Am、Ao分别为一级、二级放大后电压信号的幅值,为准确测量得到Am、Ao的值,要求两级放大后的信号均不出现饱和,即要求降低放大倍数或者选择适当的A、B值,使ARx-BR≈0,显然降低放大倍数一定会导致测量精度的降低,因此选择适当的A、B值,在保证放大倍数不降低的条件下,避免输出信号的饱和是本发明的关键核心。
若二级放大不饱和,可以通过公式的反推计算得到被测电阻为:
Figure GDA0003785636900000043
根据所述被测电阻Rx的原始标定值为R0计算激励信号A和B的理想幅值,
Figure GDA0003785636900000044
B=Const-A,Rx=R0
显然,若Ao=0,A=B,则Rx=R0=R,满足电路典型值的分析结果。
选择适当的A、B值的步骤按照如下进行:
第一次测量时,令A=B进行激励信号的加载,也就是电阻串联环两端加载幅值相等,相位相差180°的两路正弦激励信号,那么一级放大电压Vm(t)的幅值为:
Figure GDA0003785636900000045
因为一级放大倍数k1较小,此级放大不会出现Vm(t)饱和现象,根据Am的值,应用式(4)的公式计算Rx的值。由于此时放大倍数较小,计算得到的Rx值精度相对较低。
Figure GDA0003785636900000046
由式(2)可得
Figure GDA0003785636900000051
k1k2=k
(6)
Figure GDA0003785636900000052
其中,Asa为假设饱和幅值,k为总体放大倍数
假设Ao=Asa时,为正饱和;Ao=-Asa达到负饱和。考虑到单片机采用参考电压为2.5V进行ADC转换,那么信号电压转换范围为0-2.5V,可以认为这个范围的一半为信号的饱和幅值,即Asa=1.25V。
那么Rx的测量范围为:
Figure GDA0003785636900000053
若被测阻抗变化未超出Rx的测量范围,二级放大再次输出不饱和|Ao|<Asa,为提高Rx的测量精度,根据初始设置时设置的激励信号A和B的理想幅值,计算激励信号A和B幅值的上限:A和B,下限:A和B
Figure GDA0003785636900000054
B=Const-A,B=Const-A (9)
所述单片机采用参考电压进行ADC转换,参考电压为2Asa,Asa为假设饱和幅值,k=k1*k2<0,
计算当被测电阻Rx相对原始标定值为R0降低时,激励信号A和B的幅值应调整为AN=A,BN=B,对应的Rx的测量量程为
Figure GDA0003785636900000055
调整激励信号A和B的幅值生成正弦激励信号ANsinωt、-BNsinωt并加载到由被测电阻Rx与标准电阻R串联电路的两端,经同向一级放大和反向二级放大后输出电压的幅值为Ao,单片机计算得到被测电阻Rx为
Figure GDA0003785636900000056
计算当被测电阻Rx相对原始标定值为R0升高时,激励信号A和B的幅值应调整为AN=A,BN=B,对应的Rx的测量量程为
Figure GDA0003785636900000057
调整激励信号A和B的幅值生成正弦激励信号ANsinωt、-BNsinωt并加载到由被测电阻Rx与标准电阻R串联电路的两端,经同向一级放大和反向二级放大后输出电压的幅值为AO,单片机计算得到被测电阻Rx为
Figure GDA0003785636900000058
若被测阻抗变化超出Rx的测量范围,二级放大再次输出饱和|AO|>Asa,则单片机根据同向一级放大后输出点电压的幅值Am和正弦激励信号Asinωt调整前的幅值A,按下式计算激励信号A和B的幅值调整值ANx、BNx
Figure GDA0003785636900000059
Figure GDA00037856369000000510
BNx=Const-ANx
调整激励信号A和B的幅值为ANx、BNx,生成新的相应正弦激励信号后,加载到被测电阻Rx与标准电阻R串联电路的两端,得到新的电压AO,按照公式(3)重新计算精确的Rx值,
Figure GDA0003785636900000061
实施例2
一种双激励的高精度绝对阻抗测量电路的总体测量方案如下图1所示:由单片机内部的DAC模块输出一定频率的正弦激励信号,跟随输出后加载到测量电阻串联环上,中间结点输出经过两级放大后,送入单片机内部进行AD转换以及必要的信号处理,得到被测阻抗的大小。
按照图1的测量方案设计各模块的硬件电路图。
如图3、单片机主要连接硬件图
单片机选用STM32F4型的单片机,内部含有2通道的12位数模转换(DAC)模块以及3个独立模数转换(ADC)模块,可以3通道同步AD转换以及一共多达16通道的模数转换。主要电路设计原理图如下所示:这里端口PA4、PA5分别设置为DAC的模拟输出,端口PA0-PA2设置为3通道ADC输入。模拟输入输出参考电压2.5V,由稳压芯片ADR391的输出提供给VREF+引脚来实现。外部8M的晶振通过OSC_IN、OSC_OUT引脚与单片机相连,内部倍频后主频达到168M。内部还含有两个通道DMA,可以实现DAC的传送,也可以实现ADC的传输,有多达8个以上定时器,可以控制ADC与DAC的转换频率。
如图4、测量与放大电路
其次DAC跟随电路、电阻串联环以及二级放大电路设计如下:单片机通过DAC输出的激励信号通过运放的同向跟随后,加载到由标准电阻R与被测阻抗RX组成的串联电路两端,再通过同向与反向两级放大输出。
电阻串联环中的标准电阻R尽量选用与被测阻抗RX接近的电阻。若被测阻抗变化很大,R可以选择被测阻抗变化范围的中间值。
一级放大采用同向放大,由于输入阻抗大,不影响串联环中间结点电压的输出。此时放大倍数不宜过大,可以选用1-3倍,由电阻R1、R2决定。二级放大采用反向放大方式,可以采用较大的放大倍数,根据具体设计要求选用10-100倍。
电路的特点在于:(1)运放选择单电源运放,保证在单电源状态下运放能够工作,有利于测量系统的简化及对电源的要求;(2)为了适合单电源条件下电路的正常工作,输入信号应该是有直流偏置的交流激励信号,直流偏置为
Figure GDA0003785636900000062
VREF为ADC、DAC的参考电压信号。输入信号如图4左边上下量运放输出表达公式;(3)一级放大倍数为
Figure GDA0003785636900000063
因为一级放大倍数不能太大,因此R1一般远远大于R2;(4)二级放大倍数为
Figure GDA0003785636900000064
且这样的放大电路保证一、二级放大输出均含有
Figure GDA0003785636900000065
的偏置量,便于后期的ADC转换。
若取图4中A=0.29,B=0.71,R=10k,Rx=22.5k,R1=10k,R2=1k,R3=1k,R4=50k,激励信号频率为1kHz时,下图5仿真显示了四个运放的输出波形图。
软件程序与算法的设计。
1、DAC输出幅值可变的正弦信号的算法设计
利用DAC与直接数字合成(DDS)技术生成正弦激励信号的步骤包括:
(1)12位的DAC正弦数值量的计算
在确定周期点的数量n与幅值A、B的前提下,对于12位的DAC,按照如下的初始化与递推公式快速计算n个数字量计算每个点的值,并保存在数组sinArray中。
初始化公式如下:
Figure GDA0003785636900000071
如果是对幅值B对应的正弦信号初始化,这里的A换成B。δ为间隔相位。考虑到我们最终并不需要余弦量,这里的余弦值没有设置专门的数组存放。程序中后期计算阻抗结果时,注意要把A=2.5Amp/4096来替换原来的A。
递推公式如下:
Figure GDA0003785636900000072
如果是对幅值B对应的正弦信号,公式中的第三个式子的加号应改变为减号。另外,由于一个周期仅需要n个数据,因此递推后的数组中最后多余的一个值应当舍掉。
(2)DAC启动触发源的选择及初始化设置
若已知激励信号的频率,并且已知每个周期的点数为n的情况下,由TIMER定时器生成DAC触发频率。对于f=1kHz,即周期T=1000us的正弦信号,信号相邻两点的间隔时间为
Figure GDA0003785636900000073
对于168M频率的主频,DAC触发定时器分频值m为:
Figure GDA0003785636900000074
选定适当的n值,保证m为大于一定量的整数。例如n=200,则m=840。
(3)利用DMA传送方式实现2通道DAC的同步输出
为了生成连续的正弦激励信号输出,也为了不占用单片机的程序执行时间,利用DMA循环传送的方式传送存储器数组中的正弦量到DAC,输出对应的DAC模拟量。
DMA传送结构如图6所示:TIMER触发DAC转换,DAC转换完毕后向DMA申请传输数据,DMA收到请求后,从数组中取出DMA指针指向正弦数字量传送给DAC模块。在取出数字量的同时,DMA指针指向下一个数字存储空间。直到最后一个数字量传送完毕后,DMA指针再次指向第一个正弦值的存储空间,不断循环在波器上生成连续正弦信号。
利用该方法产生频率为1kHz,幅值为0.625V,直流偏置为1.25V的连续正弦信号,连接示波器。
2、ADC采样正弦信号幅值与相位拟合算法的设计
(1)三通道ADC同步采样初始化
图中端口PA0、PA1、PA2上的模拟输入量分别由单片机内部的三个ADC模块ADC1、ADC2、ADC3进行转换,并且采用同步转换的方式进行。
初始化还要求设置一定的ADC采样时间,并保证采样时间与转换时间之和小于总的采样周期。例如采样频率为100kHz,采样周期为10us,转换时间1us加上设置的采样时间3us,一共4us,小于采样周期10us,满足ADC转换要求。
(2)定时器设置采样频率
ADC的转换触发由定时器功能设置完成,保证采样时间间隔的准确性。定时器中累加器初值的设置由采样频率或采样周期决定,若采样周期为Ts(us),对于主频为168M的单片机,触发初值为168Ts。例如若采样频率为100k,采样周期10us,那么定时器累加器初值为1680.
(3)DMA设置采样数据的传输
同DAC输出相反,这里的DMA把ADC转换的数据传送到单片机内部的存储空间中去。DAC与ADC虽然都使用DMA传送,但由于使用的是不同通道的DMA,因此彼此之间不发生冲突。
(4)采样结果的正弦拟合,拟合方法及程序设计如下,
假设正弦信号假设正弦激励信号周期为T,ADC采样周期为TS,那么每个周期的采样点数
Figure GDA0003785636900000081
如T=1ms,TS=10ms,那么nT=100。
为了准确计算各采样信号幅值与相位,要求采样各个周期的信号,一般每次要采样多个周期,若采样r个周期,即一共采样点数为r·nT个。
对于采样信号Xi而言,以一级放大输出为例,理论上应满足如下的公式:
Figure GDA0003785636900000082
这里i从0开始,即i=0~r·nT-1。
展开可得:
Figure GDA0003785636900000083
对于整周期的正弦最小二乘拟合,可以推导得出如下的公式:
Figure GDA0003785636900000084
进一步,可以根据上式结果得到Am
Figure GDA0003785636900000085
令A点激励信号,一级放大输出信号,二级放大输出信号的幅值相位分别为Am1,Am2,Am3
Figure GDA0003785636900000086
对三路同步采集的ADC信号经过以上公式(11)处理后,计算出这些参数的具体大小。
在不考虑误差影响条件下,理论上应有:
Figure GDA0003785636900000087
(5)进一步计算出Am与Ao,根据a、b值计算结果如下:
实际电路计算中,由于电路分布电容及被测阻抗本身非纯电阻等的影响,且放大倍数相对图中A点激励信号不仅有大小,还有有符号的特征,所以,Am、Ao与Am2、Am3并不能直接赋值,
Figure GDA0003785636900000091
也不一定满足公式(12)所述的关系。
因此,需要运用公式(13)得到Am,AO
Figure GDA0003785636900000092
以上公式充分考虑了实际电路中输出信号与激励信号由同相向反相转化,相位连续变化的过程,相位由0向π变化过程。也就是,当相差为π/2时,可以认为电桥平衡,输出幅值为0。
实施例3、DAC两通道幅值的调整算法
(1)确定最佳A值与A的上下限值
根据一级放大信号拟合结果,确定DAC两通道模拟信号A与B的输出幅值满足A+B=1.2V,考虑到单片机采用参考电压为2.5V进行ADC转换,那么信号电压转换范围为0-2.5V,可以认为这个范围的一半为信号的饱和幅值,即=1.25V。
被测阻抗Rx对应最佳A值,表示此A值使得中间电位为0,AO=0,因此k(ARx-BR)=0,A+B=Const,对应
Figure GDA0003785636900000093
对应的上限A值,使得输出正饱和;
Figure GDA0003785636900000094
对应下限A值,使得输出负饱和。图7显示了若总体放大倍数k=10时,满足A+B=1.2V,Asa=1.25V,R在0.2R到5R变化范围下满足公式(18)的最佳A值与上下限值曲线。从图7中可以看出,若Rx=1.5R时,理想的A值处于图中点2对应纵坐标值,Rx值的测量变化范围为点1→3对应横坐标值;若能判断出Rx值将变小,可以设置A值为图中点4对应纵坐标值,即上限值AN=A+1.25/k,此时Rx的测量变化范围为对应点5→4的横坐标值,相对于1→2而言,变化范围扩大了;若判断Rx值将变大,可以设置A值为图中点6对应纵坐标值,即下限值AN=A-1.25/k,此时Rx的测量变化范围为对应点6→7的横坐标值,相对于2→3而言,变化范围扩大了。很明显,在已知阻值大小变化方向的前提下,这样的激励信号幅值设置具有更大的测量范围。
Figure GDA0003785636900000095
Figure GDA0003785636900000096
若目前模拟放大后的输出为0,并且预测Rx将变大,则可调整AN=A-1.25/k,若Rx变小,则可调整AN=A+1.25/k,这样保证正式测量前,A、B设置的初值使得被测量Rx具有最大的测量范围。
(2)根据A的上下限计算待测电阻Rx的测量范围
假设饱和幅值量为Asa,总体放大倍数为k,激励信号和为A与B,且A+B=Const。标准电阻为R,那么计算有:
Figure GDA0003785636900000101
那么Rx的范围为:
Figure GDA0003785636900000102
(3)调整算法示例
作为一个示例,电路中信号的幅值满足A+B=Const=1.2V,且满足0.2<A或B<1.0。电路中一级放大k1=1倍,二级放大k2=-50倍。
假设标准阻抗位R=10k,Rx从10k变化到9k,再变化到7k、6.9k时,A、B的加载过程为:
Rx的原始标准值为10k时,
Figure GDA0003785636900000103
时,理论上,Am=0,Ao=0,A=0.657,B=0.543;A=0.607,B=0.593,预估Rx降低时,Rx的量程为(7.595,8.987),预估Rx降低时,Rx的量程为(8.987,10.619),激励信号A和B应调整为A=A=0.657,B=B=0.543;当变化到Rx≈9k时,若未调整激励信号的幅值A和B,仍按A=B=0.6生成激励信号,按照式(2)验算得到Am≈-0.316,Ao≈1.58,显然此时|Ao|>1.25,实际输出已饱和,Rx的变化后的值不在Rx的量程内,即使按调整后的激励信号A和B加载,A=A=0.657,B=B=0.543,也得到A0=1.27>1.25,实际输出已饱和。
因此按照式(10),根据Am=-0.316,调整得到AN≈0.6316,BN≈0.5684;并用AN,BN替代A、B后,重新有输出Am≈0,Ao≈0。然后根据Ao值,按照式(3)精确计算Rx的值,精确计算后Rx的值为8.999k。
当被测电阻又变化到近似Rx≈7k时,Rx不在Rx的量程内,A的调整范围按式(18)计算理论值为(0.681,0.731),B的范围为(0.469,0.519),
若不调整初始激励信号的幅值A和B,仍按A=B=0.6生成激励信号,按照式(2)计算得到Am=-0.1059,按照式(10),调整得到AN=0.7059,BN=0.4941,落入A和B的理论范围内;并用AN,BN替代A、B后,重新有输出Am≈0,Ao≈0。然后根据Ao值,按照式(3)精确计算Rx的值,精确计算后Rx的值为6.999k。
若变化到Rx=6.9k时,A的调整范围按式(19)计算理论值为(0.685,0.735),B的范围为(0.464,0.515),若A=AN=0.7059,B=BN=0.4941,落入A和B的理论范围内,按照式(2)计算得到Am=-0.0042,Ao=0.21,显然此时|Ao|<1.25,输出未饱和,不需要调整A、B的值,按照式(3),可以计算出被测量Rx之值,Rx=0.6899k。

Claims (5)

1.一种双激励的高精度绝对阻抗测量方法,包括由单片机内部的DAC模块输出幅值不同、相位相差180°的正弦激励信号Asinωt、-Bsinωt,使激励信号的幅值A>0、B>0,A+B=Const,Const为常数,再通过运放的同向跟随后,加载到由被测电阻Rx与标准电阻R串联电路的两端,被测电阻Rx与标准电阻R的中间结点输出电压通过同向一级放大k1倍,输出电压Amsinωt,再通过反向二级放大k2倍,输出电压AOsinωt,正激励信号Asinωt及输出电压Amsinωt和AOsinωt分别送入单片机内部进行AD转换及信号处理,得到被测阻抗Rx的大小,其特征在于,在测量被测阻抗Rx的大小前,需根据以下步骤调整激励信号A和B的幅值:
步骤1.根据所述被测电阻Rx的原始标定值为R0计算激励信号A和B的理想幅值,
Figure FDA0003785636890000011
B=Const-A,Rx=R0
步骤2.根据步骤1得到的激励信号A和B的理想幅值,计算激励信号A和B幅值的上限:A和B,下限:A和B
Figure FDA0003785636890000012
B=Const-A,B=Const-A,所述单片机采用参考电压进行ADC转换,参考电压为2Asa,Asa为假设饱和幅值,k=k1*k2<0;
步骤3.计算当被测电阻Rx相对原始标定值为R0降低时,激励信号A和B的幅值应调整为AN=A,BN=B,对应的Rx的测量量程为
Figure FDA0003785636890000013
计算当被测电阻Rx相对原始标定值为R0升高时,激励信号A和B的幅值应调整为AN=A,BN=B,对应的Rx的测量量程为
Figure FDA0003785636890000014
步骤4.根据被测电阻Rx相对原始标定值的变化方向,按步骤3对应调整激励信号A和B的幅值生成正弦激励信号ANsinωt、-BNsinωt并加载到由被测电阻Rx与标准电阻R串联电路的两端,经同向一级放大和反向二级放大后输出电压的幅值为AO,单片机计算得到被测电阻Rx为
Figure FDA0003785636890000015
2.根据权利要求1所述的双激励的高精度绝对阻抗测量方法,其特征在于,若步骤4经同向一级放大和反向二级放大后输出电压的幅值AO>Asa,则单片机根据同向一级放大后输出点电压的幅值Am和正弦激励信号Asinωt调整前的幅值A,按下式计算激励信号A和B的幅值调整值ANx、BNx
Figure FDA0003785636890000016
Figure FDA0003785636890000017
BNx=Const-ANx
调整激励信号A和B的幅值为ANx、BNx,生成新的相应正弦激励信号后,加载到被测电阻Rx与标准电阻R串联电路的两端,得到新的电压AO,按照步骤4单片机计算得到被测电阻Rx为
Figure FDA0003785636890000021
3.根据权利要求1所述的双激励的高精度绝对阻抗测量方法,其特征在于,步骤4之后还包括:步骤5.若步骤4经同向一级放大和反向二级放大后输出电压的幅值AO>Asa,当被测电阻Rx相对原始标定值为R0升高时,根据步骤3计算得到Rx的测量量程上限
Figure FDA0003785636890000022
重新计算激励信号A和B的幅值调整值AN=AN+1、BN=BN+1
Figure FDA0003785636890000023
BN+1=Const-AN+1,对应的Rx的测量量程为
Figure FDA0003785636890000024
重复步骤4,若单片机采集到的反向二级放大后输出电压幅值AO不大于Asa,则单片机计算得到被测电阻Rx为
Figure FDA0003785636890000025
Figure FDA0003785636890000026
步骤6.若经步骤5单片机采集到的反向二级放大后输出电压幅值AO大于Asa,则使BN=BN+1,AN=AN+1,重复步骤5计算得到
Figure FDA0003785636890000027
和激励信号A和B的幅值调整值
Figure FDA0003785636890000028
BN=BN+2=Const-AN+2,单片机继续比较AO,若AO不大于Asa,则单片机计算得到被测电阻
Figure FDA00037856368900000217
Figure FDA0003785636890000029
步骤7.若AO>Asa,则重复步骤6,直至AO不大于Asa,单片机计算得到被测电阻
Figure FDA00037856368900000210
n>2。
4.根据权利要求1所述的双激励的高精度绝对阻抗测量方法,其特征在于,步骤4之后还包括:步骤5.若步骤4经同向一级放大和反向二级放大后输出电压的幅值AO>Asa,当被测电阻Rx相对原始标定值为R0降低时,根据步骤3计算得到Rx的测量量程下限
Figure FDA00037856368900000211
重新计算激励信号A和B的幅值调整值AN=AN+1、BN=BN+1
Figure FDA00037856368900000212
BN+1=Const-AN+1,对应的Rx的测量量程为
Figure FDA00037856368900000213
重复步骤4,若单片机采集到的反向二级放大后输出电压幅值AO不大于Asa,则单片机计算得到被测电阻Rx为
Figure FDA00037856368900000214
Figure FDA00037856368900000215
步骤6.若经步骤5单片机采集到的反向二级放大后输出电压幅值AO大于Asa,则使BN=BN+1,AN=AN+1,重复步骤5计算得到
Figure FDA00037856368900000216
和激励信号A和B的幅值调整值
Figure FDA0003785636890000031
BN=BN+2=Const-AN+2,单片机继续比较AO,若AO不大于Asa,则单片机计算得到被测电阻
Figure FDA0003785636890000032
Figure FDA0003785636890000033
步骤7.若AO>Asa,则重复步骤6,直至AO不大于Asa,单片机计算得到被测电阻
Figure FDA0003785636890000034
5.根据权利要求2所述的双激励的高精度绝对阻抗测量方法,其特征在于,同向一级放大时,放大倍数k1取1。
CN202011102344.5A 2020-10-15 2020-10-15 一种双激励的高精度绝对阻抗测量方法 Active CN112269065B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011102344.5A CN112269065B (zh) 2020-10-15 2020-10-15 一种双激励的高精度绝对阻抗测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011102344.5A CN112269065B (zh) 2020-10-15 2020-10-15 一种双激励的高精度绝对阻抗测量方法

Publications (2)

Publication Number Publication Date
CN112269065A CN112269065A (zh) 2021-01-26
CN112269065B true CN112269065B (zh) 2022-11-08

Family

ID=74338167

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011102344.5A Active CN112269065B (zh) 2020-10-15 2020-10-15 一种双激励的高精度绝对阻抗测量方法

Country Status (1)

Country Link
CN (1) CN112269065B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113074759A (zh) * 2021-03-31 2021-07-06 中国航发北京航科发动机控制系统科技有限公司 一种基于lvdt传感器的位移测量系统
CN114137305A (zh) * 2021-12-27 2022-03-04 苏州格里德医学传感技术有限公司 一种频响测试仪器和一种复阻抗测试仪器
CN115639408A (zh) * 2022-12-23 2023-01-24 中大智能科技股份有限公司 一种高精度土壤电导率自适应监测系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN2060894U (zh) * 1988-12-31 1990-08-22 浙江省送变电工程公司 光导微机型电介质测定仪
WO2015188779A1 (zh) * 2014-06-12 2015-12-17 中国人民解放军第四军医大学 一种用于电阻抗成像的高精度数据采集系统
CN106066425A (zh) * 2016-07-29 2016-11-02 中国电子科技集团公司第四十研究所 一种阻抗测量装置及其实现校准补偿的方法
CN110794321A (zh) * 2019-12-09 2020-02-14 江苏由甲申田新能源科技有限公司 一种电池储能装置的绝缘电阻监测方法及系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2544286A (en) * 2015-11-10 2017-05-17 Abb Ltd Method and apparatus for electrode impedance measurement

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN2060894U (zh) * 1988-12-31 1990-08-22 浙江省送变电工程公司 光导微机型电介质测定仪
WO2015188779A1 (zh) * 2014-06-12 2015-12-17 中国人民解放军第四军医大学 一种用于电阻抗成像的高精度数据采集系统
CN106066425A (zh) * 2016-07-29 2016-11-02 中国电子科技集团公司第四十研究所 一种阻抗测量装置及其实现校准补偿的方法
CN110794321A (zh) * 2019-12-09 2020-02-14 江苏由甲申田新能源科技有限公司 一种电池储能装置的绝缘电阻监测方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王选择 等.一种高信噪比的高阻抗测量处理电路.《仪表技术与传感器》.2017, *

Also Published As

Publication number Publication date
CN112269065A (zh) 2021-01-26

Similar Documents

Publication Publication Date Title
CN112269065B (zh) 一种双激励的高精度绝对阻抗测量方法
CN100472965C (zh) 具有最小化转换误差的a/d转换器
US9748969B1 (en) Method of operation for an oversampled data converter
CN107994903B (zh) 模数转换电路及流水线模数转换器
JPH0316320A (ja) アナログ・デジタル変換器
CN110082602B (zh) 一种全阻抗测量电路和测量装置
CN108680616B (zh) 一种数字化处理的湿度传感装置及方法
US11057047B2 (en) Ratiometric gain error calibration schemes for delta-sigma ADCs with capacitive gain input stages
CN102226825B (zh) 全数字差动电容检测装置
RU2400711C1 (ru) Тензометрический измерительный преобразователь
Linthish et al. Autonulling-based multichannel impedance measurement system for capacitive sensors
US6891430B1 (en) Techniques for signal measurement using a conditionally stable amplifier
JP3765915B2 (ja) 増幅器の温度零点補正装置
JPS59133728A (ja) A/d変換器
SU1755070A1 (ru) Устройство дл измерени и контрол температуры
JPS63224522A (ja) 増幅装置
Petrellis et al. Capacitive sensor estimation based on self-configurable reference capacitance
RU2731033C1 (ru) Тензопреобразователь давления мостового типа
RU2308688C1 (ru) Способ нормализации зарядочувствительной характеристики пьезоэлектрических преобразователей
SU754326A1 (ru) Цифровой фазометр 1
SU879602A1 (ru) Быстродействующий преобразователь отношени двух напр жений в код
Petrellis et al. A digitally configurable reference capacitance with mismatch compensation for biosensor readout circuits
CN117420359A (zh) 一种全动态范围高精度阻值测量结构及其测量方法
SU834550A1 (ru) Устройство дл измерени отношени дВуХ пЕРЕМЕННыХ НАпР жЕНий
SU809187A1 (ru) Цифрова измерительна система

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