发明名称 一种数据延拓式相关的相位差测量方法
摘要 本发明涉及信号处理领域,特别是一种数据延拓式相关的相位差测量方法。本发明的适用对象为任意同频正弦信号的相位差测量,本发明包括以下步骤:首先,对采样信号作周期预估计;然后,对数据长度是否整周期进行判断,对非整周期的数据长度作整周期数据延拓处理,使相关计算的长度恰好或者近似为整周期,从而达到提高相关法相位差测量精度的目的;最后,利用相关法相位差计算公式计算得到信号的相位差。本发明涉及的相位差测量方法实现简单、抗噪性能更强、精度更高,消除了相关法相位差测量精度受非整周期采样的影响,可结合其它方法共同提高相关法的相位差测量精度。
申请公布号 CN103926463B 申请公布日期 2016.08.17
申请号 CN201410055705.3 申请日期 2014.02.16
申请人 中国人民解放军后勤工程学院 发明人 涂亚庆;沈廷鳌;李明
分类号 G01R25/00(2006.01)I 主分类号 G01R25/00(2006.01)I
代理机构 代理人
主权项 一种数据延拓式相关的相位差测量方法,其特征在于:适用对象为任意同频正弦信号的相位差测量;该方法包括以下步骤:第一步:对信号的频率或周期进行预估计,假定获得两路短时信号的采样数据分别为x=[x<sub>1</sub>,x<sub>2</sub>,…,x<sub>N‑1</sub>,x<sub>N</sub>],y=[y<sub>1</sub>,y<sub>2</sub>,…,y<sub>N‑1</sub>,y<sub>N</sub>],P∈Z<sup>+</sup>为一个周期的采样点数,周期估计结果设为P;第二步:对数据是否整周期进行判断,假定N除以P的商是k,余数为m,若m=0,则相关长度与信号的整周期相匹配,成整数倍关系,此时数据无需处理,若m>0,说明相关长度超过k个周期但未达到k+1个周期;第三步:进行数据处理,可分两种情况进行数据处理:1)进行数据的剔除,即根据信号的周期性,从采样数据x,y中剔除整周期数据外的m个数据[x<sub>kP+1</sub>,…,x<sub>kP+m</sub>],[y<sub>kP+1</sub>,…,y<sub>kP+m</sub>],使相关长度为整周期,剔除数据后的信号序列为:x<sub>e</sub>=[x<sub>1</sub>,x<sub>2</sub>,…,x<sub>kP</sub>],y<sub>e</sub>=[y<sub>1</sub>,y<sub>2</sub>,…,y<sub>kP</sub>];2)进行数据的延拓,即根据信号的周期性,从采样数据x,y中查找P‑m个数据[x<sub>(k‑1)P+m+1</sub>,…,x<sub>kP</sub>],[y<sub>(k‑1)P+m+1</sub>,…,y<sub>kP</sub>],将其合并到原数据序列x,y之后,使相关长度为整周期,数据延拓后的信号序列为:x<sub>e</sub>=[x<sub>1</sub>,x<sub>2</sub>,…,x<sub>N</sub>,x<sub>(k‑1)P+m+1</sub>,…,x<sub>kP</sub>],y<sub>e</sub>=[y<sub>1</sub>,y<sub>2</sub>,…,y<sub>N</sub>,y<sub>(k‑1)P+m+1</sub>,…,y<sub>kP</sub>];第四步:利用相关法计算相位差,把经数据延拓处理后的信号x<sub>e</sub>,y<sub>e</sub>分别作自相关和互相关运算,再利用互相关函数计算出相位差,假定Δθ为两路信号的相位差,R<sub>xy</sub>(0)、R<sub>xx</sub>(0)、R<sub>yy</sub>(0)分别为x<sub>e</sub>,y<sub>e</sub>的互相关和自相关函数,则<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>&Delta;</mi><mi>&theta;</mi><mo>=</mo><mi>a</mi><mi>r</mi><mi>c</mi><mi>c</mi><mi>o</mi><mi>s</mi><mrow><mo>(</mo><mfrac><mrow><msub><mi>R</mi><mrow><mi>x</mi><mi>y</mi></mrow></msub><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></mrow><mrow><msqrt><mrow><msub><mi>R</mi><mrow><mi>x</mi><mi>x</mi></mrow></msub><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></mrow></msqrt><msqrt><mrow><msub><mi>R</mi><mrow><mi>y</mi><mi>y</mi></mrow></msub><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></mrow></msqrt></mrow></mfrac><mo>)</mo></mrow><mo>.</mo></mrow>]]></math><img file="FSB0000152313530000011.GIF" wi="595" he="160" /></maths>
地址 401311 重庆市沙坪坝区大学城中国人民解放军后勤工程学院