发明名称 一种大动态高精度同步连续频率测量方法
摘要 本发明属于频率测量技术领域,具体涉及一种大动态高精度同步连续频率测量方法。本发明的方法包括以下步骤:步骤1、采用被测信号对同步后的采样时间闸门信号进行计数;步骤2、采用高频脉冲计数;步骤3、采用延时单元内插计数;步骤4、误差修正,确定被测信号频率。本发明的方法细分多周期同步法中高频时钟无法识别的量化误差,提高误差分辨率,实现大动态多路信号短时间同步连续高精度测量。
申请公布号 CN103698602B 申请公布日期 2016.05.18
申请号 CN201310689682.7 申请日期 2013.12.16
申请人 北京自动化控制设备研究所 发明人 赵振涌;何杰;李婷婷;王宁;张伟斌;袁寰;邢娜
分类号 G01R23/10(2006.01)I 主分类号 G01R23/10(2006.01)I
代理机构 核工业专利中心 11007 代理人 包海燕
主权项 一种大动态高精度同步连续频率测量方法,其特征在于:包括以下步骤:步骤1、采用被测信号对同步后的采样时间闸门信号进行计数;步骤2、采用高频脉冲计数;步骤3、采用延时单元内插计数;步骤4、误差修正,确定被测信号频率;步骤1包括以下步骤:设被测信号频率为f<sub>x</sub>、周期为T<sub>x</sub>;采样时间闸门信号的频率为f<sub>s</sub>、周期为T<sub>s</sub>;将被测信号和采样时间闸门信号进行同步处理,用被测信号对同步后的采样时间闸门信号进行计数,得到采样时间闸门信号一个周期内被测信号的完整周期数:<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><msub><mi>T</mi><mi>s</mi></msub><mo>=</mo><msub><mi>N</mi><mi>x</mi></msub><msub><mi>T</mi><mi>x</mi></msub><mo>+</mo><msub><mi>&Delta;</mi><mn>0</mn></msub></mtd></mtr><mtr><mtd><msub><mi>&Delta;</mi><mn>0</mn></msub><mo>=</mo><msub><mi>T</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>T</mi><mi>i</mi></msub></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000929895900000011.GIF" wi="326" he="159" /></maths>式中,N<sub>x</sub>为同步后的采样时间闸门信号一个周期内被测信号的完整周期数;△<sub>0</sub>为第i次采样时间闸门计数后的总误差;T<sub>i‑1</sub>为第i次采样时间闸门上升沿到被测信号相邻上升沿的时间;T<sub>i</sub>为第i+1次采样时间闸门上升沿到被测信号相邻上升沿的时间;步骤2包括以下步骤:设高频脉冲频率为f<sub>0</sub>、周期为T<sub>0</sub>,采用高频脉冲对T<sub>i‑1</sub>、T<sub>i</sub>、T<sub>xi‑1</sub>、T<sub>xi</sub>进行计数:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>T</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>=</mo><msub><mi>n</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><msub><mi>f</mi><mn>0</mn></msub><mo>+</mo><msub><mi>&Delta;</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>T</mi><mi>i</mi></msub><mo>=</mo><msub><mi>n</mi><mi>i</mi></msub><msub><mi>f</mi><mn>0</mn></msub><mo>+</mo><msub><mi>&Delta;</mi><mi>i</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>T</mi><mrow><mi>x</mi><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>=</mo><msub><mi>n</mi><mrow><mi>x</mi><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><msub><mi>f</mi><mn>0</mn></msub><mo>+</mo><msub><mi>&Delta;</mi><mrow><mi>x</mi><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>T</mi><mrow><mi>x</mi><mi>i</mi></mrow></msub><mo>=</mo><msub><mi>n</mi><mrow><mi>x</mi><mi>i</mi></mrow></msub><msub><mi>f</mi><mn>0</mn></msub><mo>+</mo><msub><mi>&Delta;</mi><mrow><mi>x</mi><mi>i</mi></mrow></msub></mrow></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000929895900000021.GIF" wi="450" he="309" /></maths>其中,<maths num="0003" id="cmaths0003"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>&Delta;</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>=</mo><msub><mi>T</mi><mrow><mi>d</mi><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>T</mi><mrow><mi>s</mi><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>&Delta;</mi><mi>i</mi></msub><mo>=</mo><msub><mi>T</mi><mrow><mi>d</mi><mi>i</mi></mrow></msub><mo>-</mo><msub><mi>T</mi><mrow><mi>s</mi><mi>i</mi></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>&Delta;</mi><mrow><mi>x</mi><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>=</mo><msub><mi>T</mi><mrow><mi>x</mi><mi>d</mi><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>T</mi><mrow><mi>s</mi><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>&Delta;</mi><mrow><mi>x</mi><mi>i</mi></mrow></msub><mo>=</mo><msub><mi>T</mi><mrow><mi>x</mi><mi>d</mi><mi>i</mi></mrow></msub><mo>-</mo><msub><mi>T</mi><mrow><mi>s</mi><mi>i</mi></mrow></msub></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000929895900000022.GIF" wi="389" he="310" /></maths>式中,T<sub>xi‑1</sub>为跨第i‑1次采样时间闸门上升沿的被测信号整周期;T<sub>xi</sub>为跨第i次采样时间闸门上升沿的被测信号整周期;n<sub>i‑1</sub>、n<sub>i</sub>、n<sub>xi‑1</sub>、n<sub>xi</sub>依次为T<sub>i‑1</sub>、T<sub>i</sub>、T<sub>xi‑1</sub>、T<sub>xi</sub>内高频脉冲的个数;△<sub>i‑1</sub>为T<sub>i‑1</sub>内高频脉冲计数后形成的相对高频脉冲非整周期误差;T<sub>di‑1</sub>为跨第i‑1次采样时间闸门上升沿的高频脉冲非整周期误差;T<sub>si‑1</sub>为跨第i‑1次采样时间闸门上升沿的被测信号所在周期的下一个周期上升沿所在高频脉冲周期的非整周期误差;△<sub>i</sub>为T<sub>i</sub>内高频脉冲计数后形成的相对高频脉冲非整周期误差;T<sub>di</sub>为跨第i次采样时间闸门上升沿的高频脉冲非整周期误差;T<sub>si</sub>为跨第i次采样时间闸门上升沿的被测信号上升沿所在高频脉冲非整周期误差;△<sub>xi‑1</sub>为T<sub>xi‑1</sub>内高频脉冲计数后形成的相对高频脉冲非整周期误差;T<sub>xdi‑1</sub>为跨第i‑1次采样时间闸门上升沿的被测信号所在周期上升沿所在高频脉冲周期的非整周期误差;△<sub>xi</sub>为T<sub>xi</sub>高频脉冲计数后形成的相对高频脉冲非整周期误差;T<sub>xdi</sub>为跨第i次采样时间闸门上升沿的被测信号所在周期上升沿所在高频脉冲周期的非整周期误差;步骤3包括以下步骤:设延时单元频率为f<sub>d</sub>、周期为T<sub>d</sub>,采用延时单元对T<sub>di‑1</sub>、T<sub>si‑1</sub>、T<sub>di</sub>、T<sub>si</sub>、T<sub>xdi‑1</sub>、T<sub>xdi</sub>进行内插计数:<maths num="0004" id="cmaths0004"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><msub><mi>T</mi><mrow><mi>d</mi><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>=</mo><msub><mi>m</mi><mrow><mi>d</mi><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><msub><mi>T</mi><mi>d</mi></msub></mtd></mtr><mtr><mtd><mrow><msub><mi>T</mi><mrow><mi>s</mi><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>=</mo><msub><mi>m</mi><mrow><mi>s</mi><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><msub><mi>T</mi><mi>d</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>T</mi><mrow><mi>d</mi><mi>i</mi></mrow></msub><mo>=</mo><msub><mi>m</mi><mrow><mi>d</mi><mi>i</mi></mrow></msub><msub><mi>T</mi><mi>d</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>T</mi><mrow><mi>s</mi><mi>i</mi></mrow></msub><mo>=</mo><msub><mi>m</mi><mrow><mi>s</mi><mi>i</mi></mrow></msub><msub><mi>T</mi><mi>d</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>T</mi><mrow><mi>x</mi><mi>d</mi><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>=</mo><msub><mi>m</mi><mrow><mi>x</mi><mi>d</mi><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><msub><mi>T</mi><mi>d</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>T</mi><mrow><mi>x</mi><mi>d</mi><mi>i</mi></mrow></msub><mo>=</mo><msub><mi>m</mi><mrow><mi>x</mi><mi>d</mi><mi>i</mi></mrow></msub><msub><mi>T</mi><mi>d</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>T</mi><mn>0</mn></msub><mo>=</mo><msub><mi>m</mi><mn>0</mn></msub><msub><mi>T</mi><mi>d</mi></msub></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000929895900000031.GIF" wi="332" he="539" /></maths>式中,m<sub>di‑1</sub>、m<sub>si‑1</sub>、m<sub>di</sub>、m<sub>si</sub>、m<sub>xdi‑1</sub>、m<sub>xdi</sub>、m<sub>0</sub>依次为T<sub>di‑1</sub>、T<sub>si‑1</sub>、T<sub>di</sub>、T<sub>si</sub>、T<sub>xdi‑1</sub>、T<sub>xdi</sub>、T<sub>0</sub>内延时单元个数;步骤4包括以下步骤:采用下式计算被测信号频率:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>f</mi><mrow><mi>x</mi><mi>i</mi></mrow></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>N</mi><mrow><mi>x</mi><mi>i</mi></mrow></msub><mo>+</mo><mfrac><mrow><msub><mi>m</mi><mn>0</mn></msub><msub><mi>n</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><mrow><mo>(</mo><msub><mi>m</mi><mrow><mi>d</mi><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>m</mi><mrow><mi>s</mi><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow></mrow><mrow><msub><mi>m</mi><mn>0</mn></msub><msub><mi>n</mi><mrow><mi>x</mi><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><mrow><mo>(</mo><msub><mi>m</mi><mrow><mi>x</mi><mi>d</mi><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>m</mi><mrow><mi>s</mi><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mfrac><mrow><msub><mi>m</mi><mn>0</mn></msub><msub><mi>n</mi><mi>i</mi></msub><mo>+</mo><mrow><mo>(</mo><msub><mi>m</mi><mrow><mi>d</mi><mi>i</mi></mrow></msub><mo>-</mo><msub><mi>m</mi><mrow><mi>s</mi><mi>i</mi></mrow></msub><mo>)</mo></mrow></mrow><mrow><msub><mi>m</mi><mn>0</mn></msub><msub><mi>n</mi><mrow><mi>x</mi><mi>i</mi></mrow></msub><mo>+</mo><mrow><mo>(</mo><msub><mi>m</mi><mrow><mi>x</mi><mi>d</mi><mi>i</mi></mrow></msub><mo>-</mo><msub><mi>m</mi><mrow><mi>s</mi><mi>i</mi></mrow></msub><mo>)</mo></mrow></mrow></mfrac><mo>)</mo></mrow><mfrac><msub><mi>f</mi><mn>0</mn></msub><msub><mi>f</mi><mi>s</mi></msub></mfrac></mrow>]]></math><img file="FDA0000929895900000032.GIF" wi="1213" he="156" /></maths>N<sub>xi</sub>为第i个采样时间闸门信号内被测信号的完整周期数;f<sub>xi</sub>为第i个采样时间闸门信号内被测信号的频率。
地址 100074 北京市丰台区云岗北区西里1号院