发明名称 一种数据延拓式相关的相位差测量方法
摘要 本发明涉及信号处理领域,特别是一种数据延拓式相关的相位差测量方法。本发明的适用对象为任意同频正弦信号的相位差测量,本发明包括以下步骤:首先,对采样信号作周期预估计;然后,对数据长度是否整周期进行判断,对非整周期的数据长度作整周期数据延拓处理,使相关计算的长度恰好或者近似为整周期,从而达到提高相关法相位差测量精度的目的;最后,利用相关法相位差计算公式计算得到信号的相位差。本发明涉及的相位差测量方法实现简单、抗噪性能更强、精度更高,消除了相关法相位差测量精度受非整周期采样的影响,可结合其它方法共同提高相关法的相位差测量精度。
申请公布号 CN103926463A 申请公布日期 2014.07.16
申请号 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>=Ly<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;&theta;</mi><mo>=</mo><mi>arccos</mi><mrow><mo>(</mo><mfrac><mrow><msub><mi>R</mi><mi>xy</mi></msub><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></mrow><mrow><msqrt><msub><mi>R</mi><mi>xx</mi></msub><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msqrt><msqrt><msub><mi>R</mi><mi>yy</mi></msub><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msqrt></mrow></mfrac><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000101228840000011.GIF" wi="589" he="159" /></maths>
地址 401311 重庆市沙坪坝区大学城中国人民解放军后勤工程学院
您可能感兴趣的专利