CN113189403B - 一种自适应正交解调方法 - Google Patents
一种自适应正交解调方法 Download PDFInfo
- Publication number
- CN113189403B CN113189403B CN202110390745.3A CN202110390745A CN113189403B CN 113189403 B CN113189403 B CN 113189403B CN 202110390745 A CN202110390745 A CN 202110390745A CN 113189403 B CN113189403 B CN 113189403B
- Authority
- CN
- China
- Prior art keywords
- phase
- factor
- sequence
- orthogonal
- information
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R27/00—Arrangements for measuring resistance, reactance, impedance, or electric characteristics derived therefrom
- G01R27/02—Measuring real or complex resistance, reactance, impedance, or other two-pole characteristics derived therefrom, e.g. time constant
- G01R27/26—Measuring inductance or capacitance; Measuring quality factor, e.g. by using the resonance method; Measuring loss factor; Measuring dielectric constants ; Measuring impedance or related variables
- G01R27/2605—Measuring capacitance
Abstract
本发明公开了一种自适应正交解调方法,于包括:搭建C/V转换电路,采集响应信号形成离散序列;估算离散序列中被测信号的单周期序列长度,得到一组同相因子和正交因子;根据同相因子和正交因子获取一组被测信号的幅值信息和相位信息;对同相因子和正交因子进行补偿、得到补偿同相因子和补偿正交因子;根据补偿同相因子和正交因子可再次获取一组被测信号的幅值信息和相位信息;根据幅值信息和相位信息再次计算同相因子和正交因子的补偿余项,并得到补偿同相因子和补偿正交因子,进而再次得到一组幅值信息和相位信息,重复上述迭代过程,当相邻两次得到的幅值信息和相位信息相差小于设定阈值、则停止迭代并记录最后结果。
Description
技术领域
本发明涉及信号处理技术领域,尤其涉及一种自适应正交解调方法。
背景技术
基于交流式电容测量电路在测量仪器中有着广泛的应用,很多实验室高精度仪表就是基于交流式电容测量电路的,其等效电路图如图2所示,测量时,需要用合适频率的正弦信号作为激励信号VI,则相应信号为当电路反馈为阻性反馈时,该环节的响应快速性更好,阻性反馈意味着jωR12C36<<1,因此该环节的传递函数为进而可求得响应信号的幅值特性和相位特性为因此只要求得响应信号的幅值即可测出被测电容值,理想情况下,响应信号是与激励信号同频率的正弦信号。
可以通过使用解调算法从被测信号中计算得到幅值信息和相位信息。在实际应用上,解调运算可以分为模拟解调和数字解调(正交解调)。其中,模拟解调一般用一个模拟乘法器来实现,适用于被测信号具有恒定的相对相位,若相对相位发生变化,会导致解调结果精度下降,这是因为模拟解调输出的幅值信息中含有相位元素。因此在测量相对相位经常变化的信号时,正交解调具有更好的效果。在解调原理上,正交解调能够独立输出准确的幅值信息和相位信息,二者不会相互干扰。被测信号频率已知且恒定是正交解调能够正常使用的条件之一。但在测量实践中,被测信号的频率并非是恒定的,会在噪声信号的影响下产生一定的频率波动。
近期有研究提出使用滤波器或是快速傅里叶变换来解决被测信号频率波动的问题。而当频率小范围波动时,比如5%的波动幅度,这就意味着快速傅里叶变换以及滤波器等解决方法没有明显效果。但这样的频率波动足以对解调结果产生较为明显的影响,导致解调精度下降。解决解调精度下降问题的方法通常是采用高采样频率,这是一种根据经验的方法,高采样频率意味着使用高速AD芯片,虽然能够在一定程度上提高解调精度,但在大幅度增加成本的前提下精度提升依然有限,没有真正解决问题。
发明内容
根据现有技术存在的问题,本发明公开了一种自适应正交解调方法,具体包括如下步骤:
对于正弦信号的幅值和相位的测量,采集该信号获取离散序列;
估算离散序列中被测信号的单周期序列长度,并对比周期设定值求出频率波动系数,基于频率波动系数生成同相解调序列和正交解调序列,舍去离散序列多余部分并只保留一个周期长度。
将离散序列分别与同相解调序列和正交解调序列各项相乘并累加得到一组同相因子和正交因子;
根据同相因子和正交因子获取一组被测信号的幅值信息和相位信息;
根据幅值信息和相位信息计算同相因子和正交因子的补偿余项,对同相因子和正交因子进行补偿、得到补偿同相因子和补偿正交因子;
根据补偿同相因子和正交因子可再次获取一组被测信号的幅值信息和相位信息;通过补偿同相因子和补偿正交因子得到的幅值信息和相位信息,相比于通过同相因子和正交因子得到的幅值信息和相位信息误差较小;
根据幅值信息和相位信息再次计算同相因子和正交因子的补偿余项,并得到补偿同相因子和补偿正交因子,进而再次得到一组幅值信息和相位信息;相比于上一次得到的幅值信息和相位信息,误差进一步缩小;重复上述迭代过程,当相邻两次得到的幅值信息和相位信息相差小于设定阈值(此数值可根据实际情况需要调整设定,如0.01%),则停止迭代并记录最后结果。
进一步的,所述同相解调序列和正交解调序列采用如下方式获取:
获取离散序列的实际周期长度;
基于被测信号单周期序列长度、并对比设定的激励信号周期,求出频率波动系数δ,将被测信号序列写成:
所述补偿同相因子和补偿正交因子采用如下方式获取:
将(2.2)和(2.3)代入以下公式(2.4)和(2.5)中,计算出补偿同相因子R'和补偿正交因子I'
R'=R-c (2.4)
I'=I+s (2.5)
由于采用了上述技术方案,本发明提供的一种自适应正交解调方法,该方法从解调原理入手,从本质上降低频率波动对正交解调的影响,可以使解调结果达到更高的精度,本发明从算法上加以改进,不需要额外引入其他器件,不会提升系统复杂度,本发明能够在保证较高的精度的同时选取更低的采样频率,因此可以选择较为低速的AD芯片,可在一定程度上降低系统成本。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明中C/V转换电路等效电路图;
图2为本发明中解调算法的流程图;
图3为本发明中电容性反馈响应曲线示意图;
图4为本发明中电阻性反馈响应曲线示意图。
具体实施方式
为使本发明的技术方案和优点更加清楚,下面结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚完整的描述:
如图2所示的一种自适应正交解调方法,具体包括如下步骤:
S1:搭建C/V转换电路,用具有合适的幅值与频率的正弦激励信号激励该电路产生响应信号,使用合适的AD芯片采集响应信号形成离散序列。
S11:在测量电容领域中,交流法具有高信噪比、抗杂散能力强以及工作频带宽等优点,因此适合测量微小电容。搭建用于交流法测量微小电容的测量电路:C/V转换电路,由电路图可求得该环节的传递函数为:
由电路图可知,该电路的反馈电路是由反馈电容和反馈电阻并联而成的,同时存在电阻性和电容性两种反馈,反馈电阻Rf和反馈电容Cf的取值则可以决定哪一种反馈为主要部分,两种反馈的响应曲线有较大区别。
当jωRfCf>>1时,主要部分为电容性反馈,传递函数(3.1)化简为(3.2),此时节约响应曲线如图3所示。
当jωRfCf<<1时,主要部分为电阻性反馈,传递函数(3.1)化简为(3.3),此时节约响应曲线如图4所示。
从图中可以看出,在快速性方面电阻性反馈明显优于电阻性反馈。根据仿真所得结果,电阻性反馈电路的稳定时间在几百纳秒数量级,而电容性反馈电路的稳定时间在100微秒数量级。所以本设计选择阻性反馈为主导。并以此作为电阻和电容数值选择依据。
S12:对于激励频率的选择,根据算式(3.3)可知,激励信号频率越高,传递函数数值越大,对解调越有利;但频率越高,所需要的采集芯片频率就越高;在选择频率时应在阻性反馈的前提下选择。所以激励频率的选择应综合上述三点考虑。
S13:使用AD芯片采集S12中的响应信号生成响应信号的离散序列。上面曾提到,采集频率直接关系到离散序列长度。经本专利改进过解调算法后,采样频率可相对原有解调算法的采样频率大幅缩减,可选为激励频率的10到20倍。此外,由于频率波动的影响,在采集时无法得知响应信号具体周期,因此采集时间应适当加长,保证序列中包含一个完整的周期。
S2:基于S1中采集到的序列,估算被测信号单周期序列长度,并对比周期设定值求出频率波动系数。基于频率波动系数生成同相解调序列和正交解调序列,舍去离散序列多余部分并只保留一个周期长度。
S21:基于S1中采集到的被测信号离散序列,找出其中相邻且异号的离散点。对于正弦信号,在其零点附近可近似为线性。设Vx(m)和Vx(m+1)异号,则根据上述理论零点(穿越点)位置可近似估算为
求出离散序列中两个零点位置,间距即为半周期长度,进而可求得被测信号序列单周期实际长度Nx。
S22:基于S21得出的被测信号单周期序列长度Nx,并对比周期设定值N求出频率波动系数δ=N/Nx。对Nx取整,并仅保留序列的一个周期长度,多余部分可革除。被测信号序列也可以写成:
S3:将S1中得到的测信号离散序列和S2中的同相、正交解调序列各项相乘并累加,得到同相因子和正交因子。
S31:将S2中得到的被测信号离散序列(3.5)分别和S2中得到的同相解调序列和正交解调序列相乘累加,即:
称R、I分别为同相因子和正交因子。
化简同相、正交因子,得到:
S41:根据S3中得到的算式(3.8)和(3.9),可以得到:
S5:基于S4中得到的幅值信息和相位信息可计算S3中同相因子和正交因子的补偿余项,并对S3中的因子进行补偿,得到一组补偿同相因子和补偿正交因子。
S52:将S51的计算结果带入算式(3.14)和(3.15)中,计算出一组补偿同相因子R'和正交因子I'。
R'=R0-c (3.14)
I'=I0+s (3.15)
S7:比较两次解调结果,判断变化量是否低于阈值。若低于阈值则结束解调输出结果;若高于阈值则进行迭代计算直至变化量低于阈值。
S71:将S6得到的幅值信息K1和相位信息与S4得到的幅值信息K0和相位信息带入公式(3.16)、(3.17),得到ηK和并与幅值设定阈值ηK0和相位设定阈值进行比较。当迭代次数较少时(如小于20次),则应同时满足公式(3.18)和(3.19)时停止迭代;当迭代次数较多时(如大于等于20次),对于公式(3.18)和(3.19)满足其中之一时,停止迭代。若不符合上述两种情况,则解调继续,执行S72。
将S6的结果带回到S5中计算同相因子和正交因子的补偿余项,进而可再次获取一组补偿同相因子R'和正交因子I'。再次执行S6,可得到一组幅值信息K1和相位信息比较两次S6得到的幅值信息K和相位信息判断变化量是否低于阈值:若低于阈值,结束解调,输出最近一次S6结果;若高于阈值,重复上述过程,继续迭代,直至结果变化量低于阈值为止。
ηK≤ηK0(3.18)
通过实验来验证改进算法的可行性,利用计算机软件计算解调结果。其中,波动系数δ的范围为[0.95,1.05];相位角的取值范围为[15°,75°];幅值K设置为1000;对于原有的正交解调算法,解调序列长度取值范围为[100,1000];对于经本专利改进的算法,解调序列长度取值范围为[5,20]。经过大量数据验证,本专利的改进算法能够有效的提升精度的同时降低采样频率,列举电容测量中部分具有代表性的数据如表1表2中显示。
表1:在频率波动影响下原有的正交解调结果
表2:经本专利改进后的正交解调结果
基于交流法测量微小电容的实验,验证改进的算法的可行性。在实物环境下,本设计仅以交流法测量3pF电容为例,验证结果如表3所示。
表3实物解调结果测试
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,根据本发明的技术方案及其发明构思加以等同替换或改变,都应涵盖在本发明的保护范围之内。
Claims (4)
1.一种自适应正交解调方法,其特征在于包括:
对于正弦信号的幅值和相位的测量,采集该信号获取离散序列;
估算离散序列中被测信号的单周期序列长度,并对比周期设定值求出频率波动系数,基于频率波动系数生成同相解调序列和正交解调序列,舍去离散序列多余部分并只保留一个周期长度;
将离散序列分别与同相解调序列和正交解调序列各项相乘并累加得到一组同相因子和正交因子;
根据同相因子和正交因子获取一组被测信号的幅值信息和相位信息;
根据幅值信息和相位信息计算同相因子和正交因子的补偿余项,对同相因子和正交因子进行补偿、得到补偿同相因子和补偿正交因子;
根据补偿同相因子和正交因子可再次获取一组被测信号的幅值信息和相位信息;
根据幅值信息和相位信息再次计算同相因子和正交因子的补偿余项,并得到补偿同相因子和补偿正交因子,并采取同样的方式再次得到一组幅值信息和相位信息,重复上述迭代过程,当相邻两次得到的幅值信息和相位信息相差小于设定阈值、则停止迭代并记录最后结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110390745.3A CN113189403B (zh) | 2021-04-12 | 2021-04-12 | 一种自适应正交解调方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110390745.3A CN113189403B (zh) | 2021-04-12 | 2021-04-12 | 一种自适应正交解调方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113189403A CN113189403A (zh) | 2021-07-30 |
CN113189403B true CN113189403B (zh) | 2022-11-25 |
Family
ID=76975540
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110390745.3A Active CN113189403B (zh) | 2021-04-12 | 2021-04-12 | 一种自适应正交解调方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113189403B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5053983A (en) * | 1971-04-19 | 1991-10-01 | Hyatt Gilbert P | Filter system having an adaptive control for updating filter samples |
CN1155800A (zh) * | 1995-10-16 | 1997-07-30 | 罗拉尔航天公司 | 振幅相位键控数字通信系统用的自适应数字码元恢复 |
CN105302935A (zh) * | 2015-08-10 | 2016-02-03 | 工业和信息化部电信研究院 | 数字解调和测量分析方法 |
CN207928317U (zh) * | 2016-12-15 | 2018-10-02 | 意法半导体股份有限公司 | 阻抗测量设备 |
CN109459070A (zh) * | 2018-11-15 | 2019-03-12 | 浙江理工大学 | 一种pgc相位解调法中相位延迟提取与补偿方法 |
CN110307780A (zh) * | 2019-06-25 | 2019-10-08 | 浙江理工大学 | 基于迭代计算的pgc相位解调误差实时补偿方法 |
-
2021
- 2021-04-12 CN CN202110390745.3A patent/CN113189403B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5053983A (en) * | 1971-04-19 | 1991-10-01 | Hyatt Gilbert P | Filter system having an adaptive control for updating filter samples |
CN1155800A (zh) * | 1995-10-16 | 1997-07-30 | 罗拉尔航天公司 | 振幅相位键控数字通信系统用的自适应数字码元恢复 |
CN105302935A (zh) * | 2015-08-10 | 2016-02-03 | 工业和信息化部电信研究院 | 数字解调和测量分析方法 |
CN207928317U (zh) * | 2016-12-15 | 2018-10-02 | 意法半导体股份有限公司 | 阻抗测量设备 |
CN109459070A (zh) * | 2018-11-15 | 2019-03-12 | 浙江理工大学 | 一种pgc相位解调法中相位延迟提取与补偿方法 |
CN110307780A (zh) * | 2019-06-25 | 2019-10-08 | 浙江理工大学 | 基于迭代计算的pgc相位解调误差实时补偿方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113189403A (zh) | 2021-07-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR101294566B1 (ko) | 레졸버의 위치 오차를 적응적으로 보상하기 위한 장치 | |
JP3772121B2 (ja) | エンコーダの信号処理装置 | |
CN106506101B (zh) | 一种接收机幅频响应自动化校准方法及装置 | |
Janecki | Edge effect elimination in the recursive implementation of Gaussian filters | |
CN100353169C (zh) | 测试电子元件的方法及其仪器 | |
CN115616266B (zh) | 一种波形幅度值测量方法、装置、终端设备以及存储介质 | |
CN113484603B (zh) | 一种电能量计量系统及方法 | |
US10734927B2 (en) | Motor drive phase sector buffer, offset and gain error circuitry | |
CN106063132B (zh) | 嵌入式装置的芯片上模数转换器(adc)线性度测试 | |
CN109117816A (zh) | 基于六阶样条插值小波的信号奇异点检测方法 | |
CN102495280B (zh) | 一种抗噪音宽带频率测量方法及锁相频率计 | |
CN113189403B (zh) | 一种自适应正交解调方法 | |
US20070226602A1 (en) | Measuring device and method for measuring relative phase shifts of digital signals | |
CN112445357A (zh) | 校准方法、装置及存储介质 | |
US20110050615A1 (en) | Processing circuit for determining touch points of touch event on touch panel and related method | |
JP2011013018A (ja) | 高分解能信号処理装置および高分解能信号処理方法 | |
US8312327B2 (en) | Correcting apparatus, PDF measurement apparatus, jitter measurement apparatus, jitter separation apparatus, electric device, correcting method, program, and recording medium | |
CN114964401A (zh) | 一种基于fft的科氏流量计信号处理方法 | |
CN111624400B (zh) | 正弦信号频率测量方法 | |
CN104391175A (zh) | 具有宽频率范围揭示和保持相位信息的测频系统及其测频方法 | |
CN114325094B (zh) | 一种相位信息测量装置及方法 | |
Amira et al. | Quasi-balanced bridge method for the measurements of the impedances | |
CN112816752B (zh) | 电容互感器温度补偿方法及设备 | |
CN112556569B (zh) | 传感器温度补偿的方法、装置、设备及存储介质 | |
CN113822329B (zh) | 一种水电机组主轴摆度信号处理方法及装置 |
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 |