CN105092969B - 一种相位差估计的相频匹配方法 - Google Patents
一种相位差估计的相频匹配方法 Download PDFInfo
- Publication number
- CN105092969B CN105092969B CN201510448265.2A CN201510448265A CN105092969B CN 105092969 B CN105092969 B CN 105092969B CN 201510448265 A CN201510448265 A CN 201510448265A CN 105092969 B CN105092969 B CN 105092969B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mover
- mtr
- mtd
- 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.)
- Expired - Fee Related
Links
Abstract
本发明涉及信号处理领域,特别是一种相位差估计的相频匹配方法。本发明的适用对象为两路正弦信号的相位差估计,包括以下步骤:首先,根据正弦信号的和差角公式,利用两路正弦信号和信号频率实现两路正弦信号的90°相移,获得两路正弦信号的正交分量;其次,对两路正弦信号及其正交分量进行滑动互相关,获得互相关信号;然后,对两路正弦信号及其正交分量进行滑动自相关,获得自相关信号,实现自相关信号和互相关信号的频率匹配;最后,对互相关信号和自相关信号进行互相关,得到互相关函数,利用互相关函数相位的加权平均获得相位差估计值。本发明涉及的两路正弦信号相位差估计方法能够充分运用两路正弦信号的信息,提高正弦信号相位差估计精度,改善相位差估计的抗噪性能。
Description
技术领域
本发明涉及信号处理领域,特别是相位差估计方法。
背景技术
正弦信号的相位差估计广泛应用于战场检测、设备故障诊断、雷达、工业测量等领域。例如,科里奥利质量流量计就是通过测量两路正弦信号之间的相位差获得流体的质量流量;发电机出厂试验需要检测发电机三相输出电压的相位差,对相位差的精确估计有很高的要求,所以正弦信号相位差估计方法具有重要的研究意义和应用价值。
近年来涌现出了多种相位差估计方法,主要有:过零检测法、相关法、希尔伯特变换法、正交时延估计方法等。
(1)过零检测法(参考文献[1]:任海东,尹文庆,胡飞.基于LabVIEW的三种相位差测量法的对比分析[J].科学技术与工程,2010,1:263-268.):该方法通过计算两路正弦信号过零时刻的时间差,将时间差转换为相位差。该方法计算量小,测量速度快,但抗噪能力较差。
(2)相关法(参考文献[2]:沈廷鳌,涂亚庆,李明,张海涛.基于相关原理的相位差测量改进算法及应用[J].振动与冲击,2014,33(21):177-182.):该方法通过两路正弦信号的互相关和自相关获得信号的相位差。该方法无需预知信号的频率即可获得两路信号的相位差,但是其抗噪能力较差,且信号为非整周期采样时误差较大。
(3)希尔伯特变换法(参考文献[3]:刘维来,赵璐,王克逸,冯志华,龙潜.基于希尔伯特变换的科式流量计信号处理[J].计量学报,2013,34(5):446-451.):该方法通过对两路正弦信号进行Hilbert变换,得到解析信号,通过解析信号的互相关获得信号的相位差。该方法原理简单,但是信号非整周期采样时会出现端点效应,导致相位差误差较大,并且仅采用零序号互相关值,抗噪性能还有待提高。
(4)正交时延估计方法(参考文献[4]:Douglas L.Maskell,Graham S.Woods.Thediscrete-time quadrature subsample estimation of delay[J].IEEE transactionson Instrumentation and Measurement,2002,51(1):133-137.):通过信号时移获得正弦信号的正交分量,利用两路正弦信号的同相分量及其正交分量的互相关获得相位差估计值。该方法计算原理简单,计算量小,实时性强,但要求设置适当的采样频率才能实现正弦信号的90°相移,且由于互相关信号误差项的存在,导致正交时延估计方法为有偏估计。
综上所述,正弦信号的相位差估计方法具有重要研究意义和应用价值,但现有方法存在估计精度不高、抗噪能力较差等问题,需要提出一种估计精度高、抗噪能力强的相位差估计方法。
发明内容
本发明的目的是提出一种相位差估计的相频匹配方法,适用于两路正弦信号的相位差估计,提高相位差估计精度,改善相位差估计的抗噪性能。
本发明提出一种相位差估计的相频匹配方法。
本发明的基本思想:
首先,根据正弦信号的和差角公式,利用两路正弦信号和信号频率实现两路正弦信号的90°相移,获得两路正弦信号的正交分量;其次,对两路正弦信号及其正交分量进行滑动互相关,获得互相关信号;然后,对两路正弦信号及其正交分量进行滑动自相关,获得自相关信号;最后,对互相关信号和自相关信号进行互相关,得到互相关函数,利用互相关函数相位的加权平均获得相位差估计值。
设有两路正弦信号为:
式中sx(t)=Acos(ωt+θ1)和sy(t)=Bcos(ωt+θ2)为不含噪声信号,ω为两路信号的圆周频率,A和B分别为两路信号的幅值,θ1和θ2分别为两路信号的初相位,zx(t)和zy(t)分别为两路信号的噪声。
对于两路正弦信号x和y,设在N时刻采样得到N0点信号为:
式中xN(n)=sx,N(n)+zx,N(n),yN(n)=sy,N(n)+zy,N(n),N0为信号的长度。
第一步:根据正弦信号的和差角公式,利用两路正弦信号和信号频率计算两路正弦信号的正交分量。
对于正弦信号s(n-b)=Acos(ω(n-b)+θ),根据正弦信号的和差角公式可得:
根据式(3),正弦信号s(n)的正交分量为:
为获得两路正弦信号xN(n)和yN(n)的正交分量,根据正弦信号的和差角公式,利用两路正弦信号和信号频率计算得到两路正弦信号的正交分量为:
和的期望值为:
第二步:对两路正弦信号及其正交分量进行互相关,得到互相关信号。
为获得两路正弦信号的相位差,对两路正弦信号及其正交分量进行互相关,得到互相关信号为:
互相关信号r1,N(k)和r2,N(k)的期望值为:
由式(8)可知,互相关信号r1,N(k)和r2,N(k)的相位均含有相位差信息,所以互相关信号的有效利用有助于提高两路信号相位差的估计精度。
第三步:对两路正弦信号及其正交分量进行自相关,得到自相关信号。
由于互相关信号的相位除包括相位差和频率两部分,利用互相关信号的相位求解相位差需要克服频率对相位差估计的影响。为克服互相关信号相位中的频率对相位差估计的影响,对两路正弦信号及其正交分量进行自相关,得到自相关信号为:
式中h1,N(0)=1,h2,N(0)=0,k=0,1…N0-b-1。
自相关信号h1,N(k)和h2,N(k)的期望值为:
由式(10)可知,自相关信号h1,N(k)和h2,N(k)的初相位为零,与互相关信号频率相同,自相关信号与互相关信号的相位差即为两路正弦信号的相位差。
第四步:对互相关信号和自相关信号的互相关,得到互相关函数,利用互相关函数相位的加权平均获得相位差估计值。
式中w(k)表示权重值,
第五步:在N+1时刻采样得到两路信号的第N0+1个采样点x(N0+1)和y(N0+1),此时取N0点采样信号为:
计算得到两路正弦信号xN+1和yN+1的正交分量为:
计算得到第N+1时刻第k点互相关信号r1,N+1(k)和r2,N+1(k)为:
式中
计算得到第N+1时刻第k点自相关信号h1,N+1(k)和h2,N+1(k)为:
式中
取N=N+1,并重复步骤第四步到第五步,直到计算到最后一点采样信号。
附图说明
下面根据附图和具体实施方式对本发明进一步阐述。以两路正弦信号为例进行说明。
图1为方法的实现流程图。
图中:1表示两路正弦信号;2表示两路正弦信号的正交分量;3表示互相关信号;4表示自相关信号;5表示互相关函数;6表示相位差;7表示90°相移;8表示滑动互相关;9表示滑动自相关;10表示互相关;11表示相位加权平均。
图2为两路正弦信号。
图中:12表示第一路正弦信号;13表示第二路正弦信号。
图3为两路正弦信号的正交分量。
图中:14表示第一路正弦信号的正交分量;15表示第二路正弦信号的正交分量。
由上述分析可知,本文方法基本流程如下:
①利用式(5)计算得到两路正弦信号xN(n)和yN(n)的正交分量和
②利用式(7)计算得到两路正弦信号及其正交分量的互相关信号r1,N(k)和r2,N(k);
③利用式(9)计算得到两路正弦信号及其正交分量的自相关信号h1,N(k)和h2,N(k);
④利用式(11)计算得到相位差估计值Δθ;
⑤利用式(13)和(14)计算得到两路正弦信号xN+1和yN+1的正交分量和利用式(15)计算得到两路正弦信号及其正交分量的互相关信号r1,N+1(k)和r2,N+1(k),利用式(16)计算得到两路正弦信号及其正交分量的自相关信号h1,N+1(k)和h2,N+1(k),取N=N+1,并重复步骤④-⑤,直到计算到最后一点采样信号。
Claims (1)
1.一种相位差估计的相频匹配方法,其特征在于:适用对象为两路正弦信号的相位差估计;
该方法包括以下步骤:
第一步:根据正弦信号的和差角公式,利用两路正弦信号xN(n)、yN(n)和信号频率ω计算得到两路正弦信号的正交分量和为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>-</mo>
<mi>b</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>b</mi>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>b</mi>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>y</mi>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>-</mo>
<mi>b</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>y</mi>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>b</mi>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>b</mi>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
第二步:对两路正弦信号xN(n)、yN(n)及其正交分量进行互相关,获得互相关信号r1,N(k)和r2,N(k)为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>r</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>N</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>b</mi>
<mo>-</mo>
<mi>k</mi>
</mrow>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
<mo>+</mo>
<mi>b</mi>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>k</mi>
</mrow>
</munderover>
<mo>&lsqb;</mo>
<msub>
<mi>x</mi>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>y</mi>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>r</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>N</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>b</mi>
<mo>-</mo>
<mi>k</mi>
</mrow>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
<mo>+</mo>
<mi>b</mi>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>k</mi>
</mrow>
</munderover>
<mo>&lsqb;</mo>
<msub>
<mi>x</mi>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>y</mi>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>k</mi>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
<mn>1...</mn>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>b</mi>
<mo>-</mo>
<mn>1</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
式中N0为信号长度;
第三步:对两路正弦信号及其正交分量进行自相关,得到自相关信号为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>h</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>N</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>b</mi>
<mo>-</mo>
<mi>k</mi>
</mrow>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
<mo>+</mo>
<mi>b</mi>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>k</mi>
</mrow>
</munderover>
<mo>&lsqb;</mo>
<msub>
<mi>x</mi>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>x</mi>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>y</mi>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>y</mi>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>h</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>N</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>b</mi>
<mo>-</mo>
<mi>k</mi>
</mrow>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
<mo>+</mo>
<mi>b</mi>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>k</mi>
</mrow>
</munderover>
<mo>&lsqb;</mo>
<msub>
<mi>x</mi>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>x</mi>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>y</mi>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>y</mi>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
其中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)进行互相关,得到互相关函数,利用互相关函数相位的加权平均获得相位差估计值Δθ为:
<mrow>
<mi>&Delta;</mi>
<mi>&theta;</mi>
<mo>=</mo>
<msub>
<mi>&theta;</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&theta;</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>b</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<mi>w</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<msub>
<mi>r</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>N</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>h</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>N</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>r</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>N</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>h</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>N</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>r</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>N</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>h</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>N</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>r</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>N</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>h</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>N</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
式中w(k)表示权重值,θ1表示正弦信号xN(n)的初相位,θ2表示正弦信号yN(n)的初相位;
第五步:利用式(5)和(6)计算得到两路正弦信号xN+1和yN+1的正交分量和为
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>b</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>b</mi>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>b</mi>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>y</mi>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>b</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>y</mi>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>b</mi>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>b</mi>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>:</mo>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>:</mo>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>:</mo>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>:</mo>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
利用式(7)计算得到第N+1时刻两路正弦信号及其正交分量的第k点互相关信号r1,N+1(k)和r2,N+1(k)为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>r</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>r</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>N</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>b</mi>
<mo>-</mo>
<mi>k</mi>
</mrow>
</mfrac>
<mo>&lsqb;</mo>
<msub>
<mi>a</mi>
<mn>11</mn>
</msub>
<mo>-</mo>
<msub>
<mi>a</mi>
<mn>12</mn>
</msub>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>r</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>r</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>N</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>b</mi>
<mo>-</mo>
<mi>k</mi>
</mrow>
</mfrac>
<mo>&lsqb;</mo>
<msub>
<mi>a</mi>
<mn>21</mn>
</msub>
<mo>-</mo>
<msub>
<mi>a</mi>
<mn>22</mn>
</msub>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
式中
<mrow>
<msub>
<mi>a</mi>
<mn>12</mn>
</msub>
<mo>=</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>y</mi>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>a</mi>
<mn>21</mn>
</msub>
<mo>=</mo>
<msub>
<mi>x</mi>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>b</mi>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>b</mi>
<mo>+</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>b</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>y</mi>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>b</mi>
<mo>+</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>a</mi>
<mn>22</mn>
</msub>
<mo>=</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>y</mi>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
计算得到第N+1时刻第k点自相关信号h1,N+1(k)和h2,N+1(k)为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>h</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>h</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>N</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>b</mi>
<mo>-</mo>
<mi>k</mi>
</mrow>
</mfrac>
<mo>&lsqb;</mo>
<msub>
<mi>b</mi>
<mn>11</mn>
</msub>
<mo>-</mo>
<msub>
<mi>b</mi>
<mn>12</mn>
</msub>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>h</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>h</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>N</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>b</mi>
<mo>-</mo>
<mi>k</mi>
</mrow>
</mfrac>
<mo>&lsqb;</mo>
<msub>
<mi>b</mi>
<mn>21</mn>
</msub>
<mo>-</mo>
<msub>
<mi>b</mi>
<mn>22</mn>
</msub>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
式中
<mrow>
<msub>
<mi>b</mi>
<mn>11</mn>
</msub>
<mo>=</mo>
<msub>
<mi>x</mi>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>b</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>x</mi>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>b</mi>
<mo>+</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>b</mi>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>b</mi>
<mo>+</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>y</mi>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>b</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>y</mi>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>b</mi>
<mo>+</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>b</mi>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>b</mi>
<mo>+</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>b</mi>
<mn>12</mn>
</msub>
<mo>=</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>x</mi>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>y</mi>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>y</mi>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>b</mi>
<mn>21</mn>
</msub>
<mo>=</mo>
<msub>
<mi>x</mi>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>b</mi>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>b</mi>
<mo>+</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>b</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>x</mi>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>b</mi>
<mo>+</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>y</mi>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>b</mi>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>b</mi>
<mo>+</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>b</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>y</mi>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>b</mi>
<mo>+</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>b</mi>
<mn>22</mn>
</msub>
<mo>=</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>x</mi>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>y</mi>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mi>N</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>y</mi>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
取N=N+1,并重复步骤第四步到第五步,直到计算到最后一点采样信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510448265.2A CN105092969B (zh) | 2015-07-23 | 2015-07-23 | 一种相位差估计的相频匹配方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510448265.2A CN105092969B (zh) | 2015-07-23 | 2015-07-23 | 一种相位差估计的相频匹配方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105092969A CN105092969A (zh) | 2015-11-25 |
CN105092969B true CN105092969B (zh) | 2018-04-06 |
Family
ID=54573859
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510448265.2A Expired - Fee Related CN105092969B (zh) | 2015-07-23 | 2015-07-23 | 一种相位差估计的相频匹配方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105092969B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106771600B (zh) * | 2017-01-07 | 2019-07-30 | 中国人民解放军陆军勤务学院 | 频率失配条件下的自适应无偏相位差估计方法 |
CN109186722B (zh) * | 2018-11-14 | 2020-07-31 | 合肥工业大学 | 一种内嵌旋涡发生体的永磁式钠流量计的在线校准方法和系统 |
RU205068U1 (ru) * | 2021-04-05 | 2021-06-25 | Евгений Борисович Колесников | Устройство сдвига фазы на 90 градусов |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101315397A (zh) * | 2008-06-25 | 2008-12-03 | 中国海洋石油总公司 | 一种幅相测量方法 |
CN104142425A (zh) * | 2014-07-13 | 2014-11-12 | 中国人民解放军后勤工程学院 | 一种正弦信号频率估计的相位匹配方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100111217A1 (en) * | 2002-12-30 | 2010-05-06 | Fuente Vicente D | Method and system for reducing phase difference and doppler effects in detection and communication systems |
WO2007117525A2 (en) * | 2006-04-03 | 2007-10-18 | Wionics Research | Frequency offset correction for an ultrawideband communication system |
JP6017399B2 (ja) * | 2013-10-23 | 2016-11-02 | オリンパス株式会社 | 撮像装置及び位相差検出方法 |
-
2015
- 2015-07-23 CN CN201510448265.2A patent/CN105092969B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101315397A (zh) * | 2008-06-25 | 2008-12-03 | 中国海洋石油总公司 | 一种幅相测量方法 |
CN104142425A (zh) * | 2014-07-13 | 2014-11-12 | 中国人民解放军后勤工程学院 | 一种正弦信号频率估计的相位匹配方法 |
Non-Patent Citations (2)
Title |
---|
Phase-shift correlation method for accurate phase difference estimation in range finder;YU-ZHOU LIU等;《Applied Optics》;20150410;第54卷(第11期);3470-3477 * |
基于相关原理的相位差测量改进算法及应用;沈廷鳌等;《振动与冲击》;20141231;第33卷(第21期);177-182 * |
Also Published As
Publication number | Publication date |
---|---|
CN105092969A (zh) | 2015-11-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105556332B (zh) | 管道或管中的fmcw雷达的频散校正 | |
CN100480706C (zh) | 利用固定相移来测量同频信号相位差的方法及电路 | |
CN104697623B (zh) | 变速激励下叶片异步振动参数识别方法 | |
CN1996986A (zh) | 全相位时移相位差频谱校正法 | |
CN109029598B (zh) | 一种互相关法气体超声波流量计中基于主次包络结合处的参考波形选取方法 | |
CN105092969B (zh) | 一种相位差估计的相频匹配方法 | |
CN107843740B (zh) | 一种融合振动与声音信号频谱特征的转速测量方法 | |
CN103760418A (zh) | 一种基于希尔伯特变换的互相关相位差测量方法 | |
CN104076200A (zh) | 一种双子段相位差频率估计方法及其装置 | |
CN105022036B (zh) | 风廓线雷达风速测定方法 | |
CN111272193B (zh) | 一种基于噪声功率谱估计的mems陀螺仪在线频差识别方法 | |
CN107305223B (zh) | 一种改进的相位差频率估计方法 | |
CN106461437A (zh) | 扩展范围adc流量计 | |
CN103926463B (zh) | 一种数据延拓式相关的相位差测量方法 | |
CN104502703A (zh) | 基于频域动态模型的电力信号同步相量测量方法 | |
CN102866261A (zh) | 检测超声波在测流速中飞行时间的方法 | |
CN104061970B (zh) | 一种电磁流量信号检测方法 | |
CN102565541B (zh) | 一种用于电学层析成像系统的递推解调方法 | |
CN109030944B (zh) | 一种基于频移和dft的实正弦信号相位差估计方法 | |
CN104914439B (zh) | 一种超声波测距的双相位测量方法 | |
CN106772214A (zh) | 塔康信标模拟器方位角精度标校方法 | |
Gong et al. | On-line calibration method of SINS/DVL integrated navigation system | |
CN104483549A (zh) | 基于数据延拓和Hilbert变换的相位差测量方法 | |
CN102944360B (zh) | 一种微速差双转子系统的不平衡信号提取方法 | |
CN107064546B (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180406 Termination date: 20180723 |
|
CF01 | Termination of patent right due to non-payment of annual fee |