一种相位差估计的相频匹配方法组成比例

技术编号:12389322 阅读:91 留言:0更新日期:2015-11-25 22:32
本发明专利技术涉及信号处理领域,特别是一种相位差估计的相频匹配方法。本发明专利技术的适用对象为两路正弦信号的相位差估计,包括以下步骤:首先,根据正弦信号的和差角公式,利用两路正弦信号和信号频率实现两路正弦信号的90°相移,获得两路正弦信号的正交分量;其次,对两路正弦信号及其正交分量进行滑动互相关,获得互相关信号;然后,对两路正弦信号及其正交分量进行滑动自相关,获得自相关信号,实现自相关信号和互相关信号的频率匹配;最后,对互相关信号和自相关信号进行互相关,得到互相关函数,利用互相关函数相位的加权平均获得相位差估计值。本发明专利技术涉及的两路正弦信号相位差估计方法能够充分运用两路正弦信号的信息,提高正弦信号相位差估计精度,改善相位差估计的抗噪性能。

【技术实现步骤摘要】

本专利技术涉及信号处理领域,特别是相位差估计方法。
技术介绍
正弦信号的相位差估计广泛应用于战场检测、设备故障诊断、雷达、工业测量等领域。例如,科里奥利质量流量计就是通过测量两路正弦信号之间的相位差获得流体的质量流量;发电机出厂试验需要检测发电机三相输出电压的相位差,对相位差的精确估计有很高的要求,所以正弦信号相位差估计方法具有重要的研究意义和应用价值。近年来涌现出了多种相位差估计方法,主要有:过零检测法、相关法、希尔伯特变换法、正交时延估计方法等。(1)过零检测法(参考文献[1]:任海东,尹文庆,胡飞.基于LabVIEW的三种相位差测量法的对比分析[J].科学技术与工程,2010,1:263-268.):该方法通过计算两路正弦信号过零时刻的时间差,将时间差转换为相位差。该方法计算量小,测量速度快,但抗噪能力较差。(2)相关法(参考文献[2]:沈廷鳌,涂亚庆,李明,张海涛.基于相关原理的相位差测量改进算法及应用[J].振动与冲击,2014,33(21):177-182.):该方法通过两路正弦信号的互相关和自相关获得信号的相位差。该方法无需预知信号的频率即可获得两路信号的相位差,但是其抗噪能力较差,且信号为非整周期采样时误差较大。(3)希尔伯特变换法(参考文献[3]:刘维来,赵璐,王克逸,冯志华,龙潜.基于希尔伯特变换的科式流量计信号处理[J].计量学报,2013,34(5):446-451.):该方法通过对两路正弦信r>号进行Hilbert变换,得到解析信号,通过解析信号的互相关获得信号的相位差。该方法原理简单,但是信号非整周期采样时会出现端点效应,导致相位差误差较大,并且仅采用零序号互相关值,抗噪性能还有待提高。(4)正交时延估计方法(参考文献[4]:DouglasL.Maskell,GrahamS.Woods.Thediscrete-timequadraturesubsampleestimationofdelay[J].IEEEtransactionsonInstrumentationandMeasurement,2002,51(1):133-137.):通过信号时移获得正弦信号的正交分量,利用两路正弦信号的同相分量及其正交分量的互相关获得相位差估计值。该方法计算原理简单,计算量小,实时性强,但要求设置适当的采样频率才能实现正弦信号的90°相移,且由于互相关信号误差项的存在,导致正交时延估计方法为有偏估计。综上所述,正弦信号的相位差估计方法具有重要研究意义和应用价值,但现有方法存在估计精度不高、抗噪能力较差等问题,需要提出一种估计精度高、抗噪能力强的相位差估计方法。
技术实现思路
本专利技术的目的是提出一种相位差估计的相频匹配方法,适用于两路正弦信号的相位差估计,提高相位差估计精度,改善相位差估计的抗噪性能。本专利技术提出一种相位差估计的相频匹配方法。本专利技术的基本思想:首先,根据正弦信号的和差角公式,利用两路正弦信号和信号频率实现两路正弦信号的90°相移,获得两路正弦信号的正交分量;其次,对两路正弦信号及其正交分量进行滑动互相关,获得互相关信号;然后,对两路正弦信号及其正交分量进行滑动自相关,获得自相关信号;最后,对互相关信号和自相关信号进行互相关,得到互相关函数,利用互相关函数相位的加权平均获得相位差估计值。设有两路正弦信号为:x(t)=sx(t)+zx(t)y(t)=sy(t)+zy(t)---(1)]]>式中sx(t)=Acos(ωt+θ1)和sy(t)=Bcos(ωt+θ2)为不含噪声信号,ω为两路信号的圆周频率,A和B分别为两路信号的幅值,θ1和θ2分别为两路信号的初相位,zx(t)和zy(t)分别为两路信号的噪声。对于两路正弦信号x和y,设在N时刻采样得到N0点信号为:xN={x(1),x(2)...x(N0)本文档来自技高网...
一种相位差估计的相频匹配方法

【技术保护点】
一种相位差估计的相频匹配方法,其特征在于:适用对象为两路正弦信号的相位差估计。该方法包括以下步骤:第一步:根据正弦信号的和差角公式,利用两路正弦信号xN(n)、yN(n)和信号频率ω计算得到两路正弦信号的正交分量和为:x^N(n)=xN(n-b)-xN(n)cos(bω)sin(bω)y^N(n)=yN(n-b)-yN(n)cos(bω)sin(bω)---(1)]]>第二步:对两路正弦信号xN(n)、yN(n)及其正交分量进行互相关,获得互相关信号r1,N(k)和r2,N(k)为:r1,N(k)=1N0-b-kΣn=1+bN0-k[xN(n)yN(n+k)+x^N(n)y^N(n+k)]r2,N(k)=1N0-b-kΣn=1+bN0-k[xN(n)y^N(n+k)-x^N(n)yN(n+k)]k=0,1...N0-b-1---(2)]]>式中N0为信号长度。第三步:对两路正弦信号及其正交分量进行自相关,得到自相关信号为:h1,N(k)=1N0-b-kΣn=1+bN0-k[xN(n)xN(n+k)+x^N(n)x^N(n+k)+yN(n)yN(n+k)+y^N(n)y^N(n+k)]h2,N(k)=1N0-b-kΣn=1+bN0-k[xN(n)x^N(n+k)-x^N(n)xN(n+k)+xN(n)x^N(n+k)-x^N(n)xN(n+k)]---(3)]]>其中h1,N(0)=1,h2,N(0)=0,k=0,1…N0‑b‑1。第四步:对自相关信号h1,N(k)、h2,N(k)和互相关信号r1,N(k)、r2,N(k)进行互相关,得到互相关函数,利用互相关函数相位的加权平均获得相位差估计值Δθ为:Δθ=θ2-θ1=Σk=0N0-b-1w(k)r2,N(k)h1,N(k)-r1,N(k)h2,N(k)r1,N(k)h1,N(k)+r2,N(k)h2,N(k)---(4)]]>式中w(k)表示权重值,w(k)=N0-b-kΣk=0N0-b-1(N0-b-k),k=0,1,2...N0-b-1.]]>第五步:利用式(5)和(6)计算得到两路正弦信号xN+1和yN+1的正交分量和为x^N+1(N0)=xN+1(N0-b)-xN+1(N0)cos(bω)sin(bω)y^N+1(N0)=yN+1(N0-b)-yN+1(N0)cos(bω)sin(bω)---(5)]]>x^N+1(1:N0-1)=x^N(2:N0)y^N+1(1:N0-1)=y^N(2:N0)---(6)]]>利用式(7)计算得到第N+1时刻两路正弦信号及其正交分量的第k点互相关信号r1,N+1(k)和r2,N+1(k)为:r1,N+1(k)=r1,N(k)-1N0-b-k[a11-a12]r2,N+1(k)=r2,N(k)-1N0-b-k[a21-a22]---(7)]]>式中a11=xN(1+b)yN(1+b+k)+x^N(1+b)y^N(1+b+k),]]>a12=xN+1(N0-k)yN+1(N0)+x^N+1(N0-k)y^N+1(N0),]]>a21=xN(1+b)y^N(1+b+k)-x^N(1+b)yN(1+b+k),]]>a22=xN+1(N0-k)y^N+1(N0)-x^N+1(N0-k)yN+1(N0).]]>计算得到第N+1时刻第k点自相关信号h1,N+1(k)和h2,N+1(k)为:h1,N+1(k)=h1,N (k)-1N0-b-k[b11-b12]h2,N+1(k)=h2,N(k)-1N0-b-k[b21-b22]---(8)]]>式中b11=xN(1+b)xN(1+b+k)+x^N(1+b)x^N(1+b+k)+yN(1+b)yN(1+b+k)+y^N(1+b)y^N(1+b+k),]]>b12=xN+1(N0-k)xN+1(N0)+x^N+1(N0-k)x^N+1(N0)+yN+1(N0-k)yN+1(N0)+y^N+1(N0-k)y^N+1(N0),]]>b21=xN(1+b)x^N(1+b+k)-x^N(1+b)xN(1+b+k)+yN(1+b)y^N(1+b+k)-y^N(1+b)yN(1+b+k),]]>b22=xN+1(N0-k)x^N+1(N0)-x^N+1(N0-k)xN+1(N0)+yN+1(N0-k)y^N(N0)-y^N+1(N0-k)yN+1(N0).]]>取N=N+1,并重复步骤第四步到第五步,直到计算到最后一点采样信号。...

【技术特征摘要】
1.一种相位差估计的相频匹配方法,其特征在于:适用对象为两路正弦信号的相位差估
计。
该方法包括以下步骤:
第一步:根据正弦信号的和差角公式,利用两路正弦信号xN(n)、yN(n)和信号频率ω计
算得到两路正弦信号的正交分量和为:
x^N(n)=xN(n-b)-xN(n)cos(bω)sin(bω)y^N(n)=yN(n-b)-yN(n)cos(bω)sin(bω)---(1)]]>第二步:对两路正弦信号xN(n)、yN(n)及其正交分量进行互相关,获得互
相关信号r1,N(k)和r2,N(k)为:
r1,N(k)=1N0-b-kΣn=1+bN0-k[xN(n)yN(n+k)+x^N(n)y^N(n+k)]r2,N(k)=1N0-b-kΣn=1+bN0-k[xN(n)y^N(n+k)-x^N(n)yN(n+k)]k=0,1...N0-b-1---(2)]]>式中N0为信号长度。
第三步:对两路正弦信号及其正交分量进行自相关,得到自相关信号为:
h1,N(k)=1N0-b-kΣn=1+bN0-k[xN(n)xN(n+k)+x^N(n)x^N(n+k)+yN(n)yN(n+k)+y^N(n)y^N(n+k)]h2,N(k)=1N0-b-kΣn=1+bN0-k[xN(n)x^N(n+k)-x^N(n)xN(n+k)+xN(n)x^N(n+k)-x^N(n)xN(n+k)]---(3)]]>其中h1,N(0)=1,h2,N(0)=0,k=0,1…N0-b-1。
第四步:对自相关信号h1,N(k)、h2,N(k)和互相关信号r1,N(k)、r2,N(k)进行互相关,得到
互相关函数,利用互相关函数相位的加权平均获得相位差估计值Δθ为:
Δθ=θ2-θ1=Σk=0N0-b-1w(k)r2,N(k)h1,N(k)-r1,N(k)h2,N(k)r1,N(k)h1,N(k)+r2,N(k)h2,N(k)---(4)]]>式中w(k)表示权重值,w(k)=N0-b-kΣk=0N0-b-1(N0-b-k),k=0,1,2...N0-b-1.]]>第五步:利用式(5)和(6)计算得到两路...

【专利技术属性】
技术研发人员:涂亚庆沈艳林张海涛李明
申请(专利权)人:中国人民解放军后勤工程学院
类型:发明
国别省市:重庆;85

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1