CN104990616B - 基于级联自适应陷波器的多路不同步脉冲计数补偿方法 - Google Patents
基于级联自适应陷波器的多路不同步脉冲计数补偿方法 Download PDFInfo
- Publication number
- CN104990616B CN104990616B CN201510366719.1A CN201510366719A CN104990616B CN 104990616 B CN104990616 B CN 104990616B CN 201510366719 A CN201510366719 A CN 201510366719A CN 104990616 B CN104990616 B CN 104990616B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msup
- theta
- lanf
- 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
Landscapes
- Measuring Frequencies, Analyzing Spectra (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明公开了一种基于级联自适应陷波器的多路不同步脉冲计数补偿方法,所述方法包括:在前端固定一自适应因子为0的Liu型自适应陷波器LANF,然后将LANF输出作为Mohsen型自适应陷波器MANF的输入,再利用MANF估计频率实时调整LANF并采用级联形式固定构成级联自适应陷波器SANF,最后将SANF应用于多路不同步脉冲计数补偿,可求得周期脉冲信号任一稳态时刻相位和相应补偿计数值。本发明解决了LANF仅能处理正弦信号,MANF处理周期脉冲输入信号状态变量时偏离周期轨道、暂态响应受输入信号幅值影响、频率估计稳态误差大等问题,不仅能处理非正弦周期信号并进行多路不同步脉冲计数补偿,同时也有效提高计数精度。
Description
技术领域
本发明涉及流量计脉冲计数补偿领域,尤其涉及一种基于级联自适应陷波器的多路不同步脉冲计数补偿方法。
背景技术
随着现代工业发展,流量计需求不断增加,尤其是脉冲输出式流量计(在流场稳定情况下,脉冲输出式流量计瞬时流量正比于脉冲频率)得到空前应用,对其快速、准确检定意义重大。通用的检定装置,工控机根据流量计脉冲信号的锁相信号驱动换向器,可对流量计的脉冲信号实现整周期截取。然而,在台位装置中,并不能保证多台位中各台位流量计脉冲计数的测量精度要求。
专利CN 103176045A提出了一种基于重合脉冲计数的异频双相位重合检测方法。该专利构建了异频信号相位重合预检测及相位重合脉冲群产生电路、双相位重合检测相位偏差修正电路和门时产生电路,声称解决了传统相位重合检测易受相位噪声和触发误差的影响,存在相位重合检测相位偏差的缺点。事实上,该专利给出的相位修正公式(其中,N2、N1对应两路脉冲的计数值,T2、T1对应两路脉冲的周期),不适合多台位流量计检定装置脉冲变频率、变脉宽的应用需求。
其它相位重合检测法,例如CN 102680808A、CN 102323739 B、CN103472299 A等,用于多台位流量计检定装置都存在检定时间不确定、需额外配套复杂的硬件装置、不利于基于工控机的检测系统集成等问题。
发明内容
为解决上述存在的问题与缺陷,本发明提供一种基于级联自适应陷波器的多路不同步脉冲计数补偿方法,该方法解决了LANF仅能处理正弦信号,MANF处理周期脉冲输入信号状态变量时偏离周期轨道、暂态响应受输入信号幅值影响、频率估计稳态误差大等问题,不仅能处理非正弦周期信号并进行多路不同步脉冲计数补偿,同时也有效提高计数精度。
本发明的目的通过以下的技术方案来实现:
一种基于级联自适应陷波器的多路不同步脉冲计数补偿方法,其特征在于,所述方法包括:
A在前端固定一自适应因子为0的Liu型自适应陷波器LANF;
B将LANF输出作为Mohsen型自适应陷波器MANF的输入;
C利用MANF估计频率实时调整LANF并采用级联形式固定构成级联自适应陷波器SANF;
D将SANF应用于多路不同步脉冲计数补偿,可求得周期脉冲信号任一稳态时刻相位和相应补偿计数值。
本发明有益效果是:
本发明解决了LANF仅能处理正弦信号,MANF处理周期脉冲输入信号状态变量时偏离周期轨道、暂态响应受输入信号幅值影响、频率估计稳态误差大等问题,不仅能处理非正弦周期信号并进行多路不同步脉冲计数补偿,同时也有效提高计数精度。
附图说明
图1是本发明所述的一种基于级联自适应陷波器的多路不同步脉冲计数补偿方法流程图;
图2是本发明所述的级联自适应陷波器SANF原理框图;
图3是脉冲计数补偿原理图;
图4是MANF状态变量稳态曲线图;
图5是MANF角频率估计曲线图;
图6是SANF状态变量x1,x210Hz脉冲响应曲线;
图7是SANF角频率估计曲线图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合实施例及附图对本发明作进一步详细的描述。
下面结合实施例及附图1对本发明作进一步详细的描述:
一种基于级联自适应陷波器的多路不同步脉冲计数补偿方法具体步骤包括:
步骤10、在前端固定一自适应因子为0的Liu型自适应陷波器LANF;其中,LANF的构建步骤为:
设二阶系统状态变量为x3,x4,阻尼比为ξ,自适应增益为γ。将作为输入,则可设计LANF为:
其中,参数设定阻尼比ξ=0.7,自适应增益γ=1。
步骤20、将LANF输出作为Mohsen型自适应陷波器MANF的输入;其中,MANF的构建步骤为:
设周期脉冲幅值为A,初相位为角频率为ω,则根据傅立叶级数展开式可得信号在任意t时刻测量值f为:
设二阶系统状态变量为x1,x2,阻尼比为ξ,自适应增益为γ,角频率ω瞬态估计值为θ,那么将上式作为输入可设计MANF为:
令基波信号幅值k=2A/π,理论上,对于周期输入信号,动态系统具有唯一周期轨道:
若不考虑式自适应律,则传递函数为:
以上参数设定中,幅值A=10V,频率f=10Hz,初相阻尼比ξ=0.7,自适应增益γ=1,θ0=100。
如图4和图5分别是MANF状态变量稳态曲线图和MANF角频率估计曲线图。
步骤30、利用MANF估计频率实时调整LANF并采用级联形式固定构成级联自适应陷波器SANF;其中,SANF的构建步骤为:
根据LANF原理,状态变量x3,x4稳态周期轨道为:
由稳态周期轨道Γ2可看出,x4与输入基波信号同频、同相。按缓慢积分流形原理设置γ=0,LANF可写成:
上式实质为一共振器,如果采用级联形式固定θ于MANF频率估计值,则对于θ同频输入周期脉冲信号,动态系统可滤除θ的高次谐波。将x3作为MANF输入,可循序减弱周期方波信号所带来的强谐波干扰,提升MANF信号处理能力。图2为根据以上原理设计的SANF原理框图。
如图6和图7分别为SANF状态变量x1,x210Hz脉冲响应曲线图和SANF角频率估计曲线图。
步骤40、将SANF应用于多路不同步脉冲计数补偿,可求得周期脉冲信号任一稳态时刻相位和相应补偿计数值;具体实施方法为:
将SANF应用于多路不同步脉冲计数补偿,则可根据稳态周期轨道Γ2求得周期脉冲信号任一稳态时刻相位:
令第i路脉冲信号在计数起始时刻t0相位为计数停止时刻t1相位为上升沿检测法脉冲计数值为Ni,则补偿计数值为:
图3为计数补偿原理示意图。
虽然本发明所揭露的实施方式如上,但所述的内容只是为了便于理解本发明而采用的实施方式,并非用以限定本发明。任何本发明所属技术领域内的技术人员,在不脱离本发明所揭露的精神和范围的前提下,可以在实施的形式上及细节上作任何的修改与变化,但本发明的专利保护范围,仍须以所附的权利要求书所界定的范围为准。
Claims (2)
1.基于级联自适应陷波器的多路不同步脉冲计数补偿方法,其特征在于,所述方法包括:
A在前端固定一自适应因子为0的Liu型自适应陷波器LANF;
B将LANF输出作为Mohsen型自适应陷波器MANF的输入;
C利用MANF估计频率实时调整LANF并采用级联形式固定构成级联自适应陷波器SANF;
D将SANF应用于多路不同步脉冲计数补偿,求得周期脉冲信号任一稳态时刻相位和相应补偿计数值;
所述LANF的构建步骤为:
设二阶系统状态变量为x3,x4,阻尼比为ξ,自适应增益为γ,初相位为角频率为ω;将作为输入,则可设计LANF为:
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>3</mn>
</msub>
<mo>=</mo>
<msub>
<mi>x</mi>
<mn>4</mn>
</msub>
<mo>,</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>4</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>&xi;&theta;x</mi>
<mn>4</mn>
</msub>
<mo>-</mo>
<msup>
<mi>&theta;</mi>
<mn>2</mn>
</msup>
<msub>
<mi>x</mi>
<mn>3</mn>
</msub>
<mo>+</mo>
<msup>
<mi>&theta;</mi>
<mn>2</mn>
</msup>
<mi>y</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mover>
<mi>&theta;</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>&gamma;x</mi>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>&theta;</mi>
<mn>2</mn>
</msup>
<mi>y</mi>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>&zeta;&theta;x</mi>
<mn>4</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,参数设定阻尼比ξ=0.7,自适应增益γ=1;
所述MANF的构建步骤为:
设周期脉冲幅值为A,初相位为角频率为ω,则根据傅立叶级数展开式可得信号在任意t时刻测量值f为:
设二阶系统状态变量为x1,x2,阻尼比为ξ,自适应增益为γ,角频率ω瞬态估计值为θ,那么将上式(3)作为输入可设计MANF为:
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</msub>
<mo>=</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>2</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>&xi;&theta;x</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msup>
<mi>&theta;</mi>
<mn>2</mn>
</msup>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<mn>2</mn>
<mi>&xi;</mi>
<mi>&theta;</mi>
<mi>f</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mover>
<mi>&theta;</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>&gamma;x</mi>
<mn>1</mn>
</msub>
<mi>&theta;</mi>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
令基波信号幅值k=2A/π,理论上,对于周期输入信号,动态系统(4)、(5)具有唯一周期轨道:
若不考虑式(5)自适应律,公式(4)传递函数为:
<mrow>
<msub>
<mi>H</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mi>&xi;</mi>
<mi>&theta;</mi>
</mrow>
<mrow>
<msup>
<mi>s</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<mn>2</mn>
<mi>&xi;</mi>
<mi>&theta;</mi>
<mi>s</mi>
<mo>+</mo>
<msup>
<mi>&theta;</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
以上参数设定中,幅值A=10V,频率f=10Hz,初相阻尼比ξ=0.7,自适应增益γ=1,θ0=100;
所述SANF的构建步骤为:
根据LANF原理,状态变量x3,x4稳态周期轨道为:
由稳态周期轨道Γ2可看出,x4与输入基波信号同频、同相;按缓慢积分流形原理设置γ=0,LANF可写成:
<mrow>
<msub>
<mi>H</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<msup>
<mi>&theta;</mi>
<mn>2</mn>
</msup>
<mrow>
<msup>
<mi>s</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<mn>2</mn>
<mi>&xi;</mi>
<mi>&theta;</mi>
<mi>s</mi>
<mo>+</mo>
<msup>
<mi>&theta;</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,ω为角频率、θ为瞬态估计值、ξ为阻尼比。
2.根据权利要求1所述的基于级联自适应陷波器的多路不同步脉冲计数补偿方法,其特征在于,所述步骤D中,将SANF应用于多路不同步脉冲计数补偿方法为:
将SANF应用于多路不同步脉冲计数补偿,则根据稳态周期轨道Γ2求得周期脉冲信号任一稳态时刻相位:
令第i路脉冲信号在计数起始时刻t0相位为计数停止时刻t1相位为上升沿检测法脉冲计数值为Ni,则补偿计数值为:
其中,x3,x4为二阶系统状态变量,Nci为周期脉冲信号任一稳态时刻相应补偿计数值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510366719.1A CN104990616B (zh) | 2015-06-26 | 2015-06-26 | 基于级联自适应陷波器的多路不同步脉冲计数补偿方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510366719.1A CN104990616B (zh) | 2015-06-26 | 2015-06-26 | 基于级联自适应陷波器的多路不同步脉冲计数补偿方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104990616A CN104990616A (zh) | 2015-10-21 |
CN104990616B true CN104990616B (zh) | 2018-01-19 |
Family
ID=54302458
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510366719.1A Active CN104990616B (zh) | 2015-06-26 | 2015-06-26 | 基于级联自适应陷波器的多路不同步脉冲计数补偿方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104990616B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
PL3491629T3 (pl) * | 2016-07-26 | 2020-12-28 | Alert Systems Aps | Sposób, aparat i system do wykrywania metalowych obiektów w strefie detekcji |
CN107782984A (zh) * | 2017-12-08 | 2018-03-09 | 国家电网公司 | 一种测量电网中各序分量的方法 |
CN112307907B (zh) * | 2020-10-09 | 2023-09-19 | 中国人民解放军陆军勤务学院 | 自适应陷波器偏差补偿的无偏频率估计方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5555190A (en) * | 1995-07-12 | 1996-09-10 | Micro Motion, Inc. | Method and apparatus for adaptive line enhancement in Coriolis mass flow meter measurement |
CN1151637A (zh) * | 1995-08-08 | 1997-06-11 | 三菱电机株式会社 | 频率合成器 |
CN1382995A (zh) * | 2002-04-24 | 2002-12-04 | 清华大学 | 一种频率和相位的数字测量方法 |
CN103248356A (zh) * | 2013-05-20 | 2013-08-14 | 上海理工大学 | 一种基于采用锁相环脉冲插值技术的计数器及实现方法 |
CN104154916A (zh) * | 2013-08-30 | 2014-11-19 | 北京航天发射技术研究所 | 一种基于激光陀螺捷联惯组的车载定位设备 |
CN104266732A (zh) * | 2014-10-24 | 2015-01-07 | 广州能源检测研究院 | 双台位流量计检定装置脉冲计数精度补偿方法 |
-
2015
- 2015-06-26 CN CN201510366719.1A patent/CN104990616B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5555190A (en) * | 1995-07-12 | 1996-09-10 | Micro Motion, Inc. | Method and apparatus for adaptive line enhancement in Coriolis mass flow meter measurement |
CN1151637A (zh) * | 1995-08-08 | 1997-06-11 | 三菱电机株式会社 | 频率合成器 |
CN1382995A (zh) * | 2002-04-24 | 2002-12-04 | 清华大学 | 一种频率和相位的数字测量方法 |
CN103248356A (zh) * | 2013-05-20 | 2013-08-14 | 上海理工大学 | 一种基于采用锁相环脉冲插值技术的计数器及实现方法 |
CN104154916A (zh) * | 2013-08-30 | 2014-11-19 | 北京航天发射技术研究所 | 一种基于激光陀螺捷联惯组的车载定位设备 |
CN104266732A (zh) * | 2014-10-24 | 2015-01-07 | 广州能源检测研究院 | 双台位流量计检定装置脉冲计数精度补偿方法 |
Non-Patent Citations (2)
Title |
---|
基于级联自适应陷波器的正选频率估计;张磊 等;《信号处理》;20060630;第22卷(第3期);全文 * |
自适应陷波器的科氏流量计信号频率跟踪方法;涂亚庆 等;《重庆大学学报》;20111031;第34卷(第10期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN104990616A (zh) | 2015-10-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101806832B (zh) | 一种低频率信号的频率测量方法 | |
CN104990616B (zh) | 基于级联自适应陷波器的多路不同步脉冲计数补偿方法 | |
CN101833036B (zh) | 一种交流电的瞬时相位测量方法 | |
US10429418B2 (en) | Frequency measurement device, frequency measurement method and under frequency load shedding device | |
CN103616814A (zh) | 一种基于fpga的同步采样时钟闭环校正方法和系统 | |
CN105224811B (zh) | 一种基于反馈迭代式频率跟踪的pmu动态数据处理方法 | |
CN103018555B (zh) | 一种高精度的电力参数软件同步采样方法 | |
CN105487034A (zh) | 一种0.05级电子式互感器校验方法及系统 | |
CN104897250B (zh) | 一种抗强谐波干扰的多台位流量计脉冲计数补偿方法 | |
CN107561476A (zh) | 一种计算cvt等效电路的电路参数的方法与装置 | |
CN113675850B (zh) | 一种基于非线性鲁棒估计的电网信息快速精准感知方法 | |
JP5000441B2 (ja) | 電力測定方法および電力測定装置 | |
CN102540252A (zh) | 基于互相关的高精度中值叠加方法 | |
CN104155521A (zh) | 相位差的确定方法和装置 | |
CN102495285B (zh) | 对称窗函数功率重心估计电力谐波参数的方法 | |
CN106597097A (zh) | 一种高精度频率测量方法 | |
CN103528634A (zh) | 一种科里奥利质量流量计云传输数字信号处理装置及方法 | |
CN101556173A (zh) | 一种基于小波降噪的科里奥利质量流量计数字解算装置及方法 | |
CN103543331B (zh) | 一种计算电信号谐波和间谐波的方法 | |
CN104266732B (zh) | 双台位流量计检定装置脉冲计数精度补偿方法 | |
CN103605904B (zh) | 基于误差估算的自补偿电力系统幅值算法 | |
CN103634003B (zh) | Osc频率自动校准电路及自动校准方法 | |
CN104020350B (zh) | 一种克服频率摄动的电压基波分量检测方法 | |
CN106291102A (zh) | 一种频标比对装置及方法 | |
CN103884910B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |