CN106027117B - 一种时延与多普勒频移联合估计方法 - Google Patents
一种时延与多普勒频移联合估计方法 Download PDFInfo
- Publication number
- CN106027117B CN106027117B CN201610287197.0A CN201610287197A CN106027117B CN 106027117 B CN106027117 B CN 106027117B CN 201610287197 A CN201610287197 A CN 201610287197A CN 106027117 B CN106027117 B CN 106027117B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msup
- tau
- msubsup
- sigmoid
- 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
- 238000000034 method Methods 0.000 title claims abstract description 14
- 230000004087 circulation Effects 0.000 claims abstract description 49
- 238000006243 chemical reaction Methods 0.000 claims abstract description 31
- 239000000654 additive Substances 0.000 claims description 10
- 230000000996 additive effect Effects 0.000 claims description 10
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 2
- 230000005540 biological transmission Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 abstract description 3
- 230000000694 effects Effects 0.000 description 2
- 125000004122 cyclic group Chemical group 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000000151 deposition Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000116 mitigating effect Effects 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
- H04B1/69—Spread spectrum techniques
- H04B1/707—Spread spectrum techniques using direct sequence modulation
- H04B1/7097—Interference-related aspects
- H04B1/71—Interference-related aspects the interference being narrowband interference
- H04B1/7101—Interference-related aspects the interference being narrowband interference with estimation filters
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
- H04B1/69—Spread spectrum techniques
- H04B1/7163—Spread spectrum techniques using impulse radio
- H04B1/719—Interference-related aspects
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Radar Systems Or Details Thereof (AREA)
- Organic Low-Molecular-Weight Compounds And Preparation Thereof (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
Abstract
本发明属于循环平稳信号处理技术领域,提供一种基于广义Sigmoid变换循环模糊函数的时延和多普勒频移联合估计方法。其特征是分别求取参考信号的可调Sigmoid变换循环自相关,待测信号的可调Sigmoid变换循环自相关,参考信号与待测信号之间的可调Sigmoid变换循环互相关和待测信号与参考信号之间的可调Sigmoid变换循环互相关,然后求取广义Sigmoid变换循环模糊函数和其绝对值,搜索其绝对值最大值所对应的横纵坐标,获得时延和多普勒频移的估计值。实验证明,本发明能够在脉冲噪声和同频干扰同时存在的情况下对时延和多普勒频移进行准确估计,且具有较强的抗脉冲噪声能力。
Description
技术领域
本发明属于循环平稳信号处理技术领域,涉及到时延和多普勒频移联合估计的方法,特别涉及一种使用广义Sigmoid变换循环模糊函数的时延与多普勒频移联合估计的方法。
背景技术
大多数通信、雷达信号都具有循环平稳特性,利用这一特性可以有效去除同频带干扰信号的影响。循环模糊函数就是应用了循环平稳理论的一个经典方法,用于同频干扰下的时延和多普勒频移联合估计具有很好的效果。但是该方法是基于噪声服从高斯分布的假设下的,当噪声中包含脉冲特性时算法性能将严重退化或失效。
基于分数低阶统计量理论的分数低阶循环模糊函数,可以同时在脉冲噪声和同频干扰下工作,但也有其固有的缺点:分数低阶的阶数需要噪声的先验知识,否则阶数的不当选择会影响方法的估计效果,另外,分数低阶统计量在脉冲性较强时,对脉冲噪声抑制能力不足,从而导致估计误差增大。故本发明提出一种可调Sigmoid变换循环相关,在可调Sigmoid循环相关的基础上,应用广义Sigmoid变换循环模糊函数在脉冲噪声和同频干扰共存条件下进行时延和多普勒频移的联合估计,并且可以通过调节可调Sigmoid变换循环相关的倾斜参数,调整对不同信号的衰减程度,以提高方法的适用范围。
发明内容
针对现有技术的不足,本发明提供一种利用广义Sigmoid变换循环模糊函数进行时延和多普勒频移的联合估计方法。该方法采用可调Sigmoid变换可以通过倾斜参数的调整有针对性地对噪声中的脉冲性进行抑制,利用同载频不同波特率的信号具有不同循环频率的特性,可以通过循环频率有效区分待测信号和同频干扰信号。故本方法可以在脉冲噪声和同频干扰同时存在时对时延和多普勒频移进行联合估计。
本发明的技术方案为:
一种基于广义Sigmoid变换循环模糊函数的时延与多普勒频移联合估计的方法,主要包括以下步骤:
第一步,计算广义Sigmoid变换循环模糊函数
1.1)利用可调Sigmoid循环自相关公式计算参考信号的可调Sigmoid循环自相关;
1.2)利用可调Sigmoid循环自相关公式计算待测信号的可调Sigmoid循环自相关;
1.3)利用可调Sigmoid循环互相关公式计算参考信号与待测信号之间的可调Sigmoid循环互相关。
1.4)利用可调Sigmoid循环互相关公式计算待测信号与参考信号之间的可调Sigmoid循环互相关。
1.5)根据广义Sigmoid变换循环模糊函数公式计算获得广义Sigmoid变换循环模糊函数。
第二步,参数估计
2.1计算1.5)结果中广义Sigmoid变换循环模糊函数的绝对值;
2.2搜索2.1)结果所对应的时延值和多普勒频移值,即为所求的估计值。
本方法能够在脉冲噪声和同频干扰同时存在的条件下对时延和多普勒频移正确估计,且抗脉冲噪声能力较强,满足实际需求。
附图说明
图1是本发明的算法流程图;
图2是本发明中广义Sigmoid变换循环模糊函数以BPSK信号为例的三维图;
图3(a)是本发明中广义Sigmoid变换循环模糊函数在频率等于真实频移时的时间轴截面;
图3(b)是本发明中广义Sigmoid变换循环模糊函数在时延等于真实时延时的频率轴截面。
具体实施方式
为使本发明实施例的目的、技术方案及其优点更加清楚,下面结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚完整的描述,整体算法流程图如图1所示:
第一步,计算广义Sigmoid变换循环模糊函数;
1.1利用可调Sigmoid变换循环自相关公式计算参考信号x(t)在时间区间 [-T/2,T/2]上的可调Sigmoid变换循环自相关如公式(1)所示:
采集数据为有限长度时的可调Sigmoid变换循环自相关可由公式(2)估计得到;
其中,ε表示循环频率,由s(t)的循环频率确定;τ表示时延;x(t)为接收到的参考信号,x(t)=s(t)+w1(t)+sI(t);w1(t)为加性噪声,设加性噪声为服从位置参数a=0、对称参数β=0的Alpha稳定分布噪声,若加性噪声为高斯分布或无干扰信号,不影响后续的估计步骤及估计结果;s(t)为具有循环平稳特性的发射信号;sI(t)为干扰信号,sI(t)和s(t)的载频相同,波特率不同;加性噪声项w1(t)和干扰项sI(t)为脉冲噪声与同频干扰并存下的信号模型; Sigmoid[x(t)]表示如公式(3)所示的对x(t)进行可调Sigmoid变换;
1.2)与1.1)类似,利用可调Sigmoid变换循环自相关公式计算待测信号y(t) 在时间区间[-T/2,T/2]上的可调Sigmoid变换循环自相关如公式(4) 所示:
采集数据为有限长度时的可调Sigmoid变换循环自相关可由公式(5)估计得到;
其中,为待测信号;D为待估计的时间延迟;fd为待估计的多普勒频移;w2(t)为加性噪声;
1.3利用可调Sigmoid循环互相关公式计算参考信号与待测信号之间的可调Sigmoid循环互相关
采集数据为有限长度时的估计公式为:
1.4利用可调Sigmoid循环互相关公式计算待测信号与参考信号之间的可调Sigmoid循环互相关
采集数据为有限长度时的估计公式为
1.5根据广义Sigmoid循环模糊函数公式,计算得到广义Sigmoid循环模糊函数u和f分别代表时延和多普勒频移;
第二步,参数估计
2.1)计算公式(10)结果中广义Sigmoid循环模糊函数的绝对值;
2.2)搜索广义Sigmoid循环模糊函数的绝对值最大值的一组时延值u和多普勒频移f,为待测信号的估计值如公式(11)所示。
Claims (1)
1.一种基于广义Sigmoid变换循环模糊函数的时延与多普勒频移联合估计的方法,其特征在于,包括以下步骤:
第一步,计算广义Sigmoid变换循环模糊函数
1.1)利用公式(1)可调Sigmoid变换循环自相关公式计算参考信号x(t)在时间区间[-T/2,T/2]上的可调Sigmoid变换循环自相关
<mrow>
<msubsup>
<mi>R</mi>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>T</mi>
<mi>S</mi>
</mrow>
<mi>&epsiv;</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mover>
<mo>=</mo>
<mi>&Delta;</mi>
</mover>
<munder>
<mrow>
<mi>l</mi>
<mi>i</mi>
<mi>m</mi>
</mrow>
<mrow>
<mi>T</mi>
<mo>&RightArrow;</mo>
<mi>&infin;</mi>
</mrow>
</munder>
<mfrac>
<mn>1</mn>
<mi>T</mi>
</mfrac>
<munderover>
<mo>&Integral;</mo>
<mrow>
<mo>-</mo>
<mi>T</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
<mrow>
<mi>T</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
</munderover>
<mi>S</mi>
<mi>i</mi>
<mi>g</mi>
<mi>m</mi>
<mi>o</mi>
<mi>i</mi>
<mi>d</mi>
<mo>&lsqb;</mo>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>+</mo>
<mi>&tau;</mi>
<mo>/</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<msup>
<mi>Sigmoid</mi>
<mo>*</mo>
</msup>
<mo>&lsqb;</mo>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>&tau;</mi>
<mo>/</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mi>&epsiv;</mi>
<mi>t</mi>
</mrow>
</msup>
<mi>d</mi>
<mi>t</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
由公式(2)估计得到采集数据为有限长度时的可调Sigmoid变换循环自相关
<mrow>
<msubsup>
<mover>
<mi>R</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>T</mi>
<mi>S</mi>
</mrow>
<mi>&epsiv;</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>T</mi>
</mfrac>
<munderover>
<mo>&Integral;</mo>
<mrow>
<mo>-</mo>
<mi>T</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
<mrow>
<mi>T</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
</munderover>
<mi>S</mi>
<mi>i</mi>
<mi>m</mi>
<mi>o</mi>
<mi>i</mi>
<mi>d</mi>
<mo>&lsqb;</mo>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>+</mo>
<mi>&tau;</mi>
<mo>/</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<msup>
<mi>Simoid</mi>
<mo>*</mo>
</msup>
<mo>&lsqb;</mo>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>&tau;</mi>
<mo>/</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mi>&epsiv;</mi>
<mi>t</mi>
</mrow>
</msup>
<mi>d</mi>
<mi>t</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,ε表示循环频率,由s(t)的循环频率确定;τ表示时延;x(t)为接收到的参考信号,x(t)=s(t)+w1(t)+sI(t);w1(t)为加性噪声,设加性噪声为服从位置参数a=0、对称参数β=0的Alpha稳定分布噪声,若加性噪声为高斯分布或无干扰信号,不影响后续的估计步骤及估计结果;s(t)为具有循环平稳特性的发射信号;sI(t)为干扰信号,sI(t)和s(t)的载频相同,波特率不同;加性噪声项w1(t)和干扰项sI(t)为脉冲噪声与同频干扰并存下的信号模型;公式(3)表示对x(t)进行可调Sigmoid变换;
<mrow>
<mi>S</mi>
<mi>i</mi>
<mi>g</mi>
<mi>m</mi>
<mi>o</mi>
<mi>i</mi>
<mi>d</mi>
<mo>&lsqb;</mo>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>=</mo>
<mfrac>
<mn>2</mn>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mi>&lambda;</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mfrac>
<mo>-</mo>
<mn>1</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
1.2)利用公式(4)可调Sigmoid变换循环自相关公式计算待测信号y(t)在时间区间[-T/2,T/2]上的可调Sigmoid变换循环自相关
<mrow>
<msubsup>
<mi>R</mi>
<mrow>
<mi>y</mi>
<mo>,</mo>
<mi>T</mi>
<mi>S</mi>
</mrow>
<mi>&epsiv;</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mover>
<mo>=</mo>
<mi>&Delta;</mi>
</mover>
<munder>
<mi>lim</mi>
<mrow>
<mi>T</mi>
<mo>&RightArrow;</mo>
<mi>&infin;</mi>
</mrow>
</munder>
<mfrac>
<mn>1</mn>
<mi>T</mi>
</mfrac>
<munderover>
<mo>&Integral;</mo>
<mrow>
<mo>-</mo>
<mi>T</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
<mrow>
<mi>T</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
</munderover>
<mi>S</mi>
<mi>i</mi>
<mi>g</mi>
<mi>m</mi>
<mi>o</mi>
<mi>i</mi>
<mi>d</mi>
<mo>&lsqb;</mo>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>+</mo>
<mi>&tau;</mi>
<mo>/</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<msup>
<mi>Sigmoid</mi>
<mo>*</mo>
</msup>
<mo>&lsqb;</mo>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>&tau;</mi>
<mo>/</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mi>&epsiv;</mi>
<mi>t</mi>
</mrow>
</msup>
<mi>d</mi>
<mi>t</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
采集数据为有限长度时的可调Sigmoid变换循环自相关可由公式(5)估计得到;
<mrow>
<msubsup>
<mover>
<mi>R</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>y</mi>
<mo>,</mo>
<mi>T</mi>
<mi>S</mi>
</mrow>
<mi>&epsiv;</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>T</mi>
</mfrac>
<munderover>
<mo>&Integral;</mo>
<mrow>
<mo>-</mo>
<mi>T</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
<mrow>
<mi>T</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
</munderover>
<mi>S</mi>
<mi>i</mi>
<mi>g</mi>
<mi>m</mi>
<mi>o</mi>
<mi>i</mi>
<mi>d</mi>
<mo>&lsqb;</mo>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>+</mo>
<mi>&tau;</mi>
<mo>/</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<msup>
<mi>Sigmoid</mi>
<mo>*</mo>
</msup>
<mo>&lsqb;</mo>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>&tau;</mi>
<mo>/</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mi>&epsiv;</mi>
<mi>t</mi>
</mrow>
</msup>
<mi>d</mi>
<mi>t</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,为待测信号;D为待估计的时间延迟;fd为待估计的多普勒频移;w2(t)为加性噪声;
1.3)利用可调Sigmoid循环互相关公式计算参考信号与待测信号之间的可调Sigmoid循环互相关
<mrow>
<msubsup>
<mi>R</mi>
<mrow>
<mi>x</mi>
<mi>y</mi>
<mo>,</mo>
<mi>T</mi>
<mi>S</mi>
</mrow>
<mi>&epsiv;</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mover>
<mo>=</mo>
<mi>&Delta;</mi>
</mover>
<munder>
<mrow>
<mi>l</mi>
<mi>i</mi>
<mi>m</mi>
</mrow>
<mrow>
<mi>T</mi>
<mo>&RightArrow;</mo>
<mi>&infin;</mi>
</mrow>
</munder>
<mfrac>
<mn>1</mn>
<mi>T</mi>
</mfrac>
<munderover>
<mo>&Integral;</mo>
<mrow>
<mo>-</mo>
<mi>T</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
<mrow>
<mi>T</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
</munderover>
<mi>S</mi>
<mi>i</mi>
<mi>g</mi>
<mi>m</mi>
<mi>o</mi>
<mi>i</mi>
<mi>d</mi>
<mo>&lsqb;</mo>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>+</mo>
<mi>&tau;</mi>
<mo>/</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<msup>
<mi>Sigmoid</mi>
<mo>*</mo>
</msup>
<mo>&lsqb;</mo>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>&tau;</mi>
<mo>/</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mi>&epsiv;</mi>
<mi>t</mi>
</mrow>
</msup>
<mi>d</mi>
<mi>t</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
采集数据为有限长度时的估计公式为:
<mrow>
<msubsup>
<mover>
<mi>R</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>x</mi>
<mi>y</mi>
<mo>,</mo>
<mi>T</mi>
<mi>S</mi>
</mrow>
<mi>&epsiv;</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>T</mi>
</mfrac>
<munderover>
<mo>&Integral;</mo>
<mrow>
<mo>-</mo>
<mi>T</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
<mrow>
<mi>T</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
</munderover>
<mi>S</mi>
<mi>i</mi>
<mi>g</mi>
<mi>m</mi>
<mi>o</mi>
<mi>i</mi>
<mi>d</mi>
<mo>&lsqb;</mo>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>+</mo>
<mi>&tau;</mi>
<mo>/</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<msup>
<mi>Sigmoid</mi>
<mo>*</mo>
</msup>
<mo>&lsqb;</mo>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>&tau;</mi>
<mo>/</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mi>&epsiv;</mi>
<mi>t</mi>
</mrow>
</msup>
<mi>d</mi>
<mi>t</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
1.4)利用可调Sigmoid循环互相关公式计算待测信号与参考信号之间的可调Sigmoid循环互相关
<mrow>
<msubsup>
<mi>R</mi>
<mrow>
<mi>y</mi>
<mi>x</mi>
<mo>,</mo>
<mi>T</mi>
<mi>S</mi>
</mrow>
<mi>&epsiv;</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mi>&Delta;</mi>
<mo>=</mo>
<munder>
<mrow>
<mi>l</mi>
<mi>i</mi>
<mi>m</mi>
</mrow>
<mrow>
<mi>T</mi>
<mo>&RightArrow;</mo>
<mi>&infin;</mi>
</mrow>
</munder>
<mfrac>
<mn>1</mn>
<mi>T</mi>
</mfrac>
<munderover>
<mo>&Integral;</mo>
<mrow>
<mo>-</mo>
<mi>T</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
<mrow>
<mi>T</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
</munderover>
<mi>S</mi>
<mi>i</mi>
<mi>g</mi>
<mi>m</mi>
<mi>o</mi>
<mi>i</mi>
<mi>d</mi>
<mo>&lsqb;</mo>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>+</mo>
<mi>&tau;</mi>
<mo>/</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<msup>
<mi>Sigmoid</mi>
<mo>*</mo>
</msup>
<mo>&lsqb;</mo>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>&tau;</mi>
<mo>/</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mi>&epsiv;</mi>
<mi>t</mi>
</mrow>
</msup>
<mi>d</mi>
<mi>t</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
采集数据为有限长度时的估计公式为
<mrow>
<msubsup>
<mover>
<mi>R</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>y</mi>
<mi>x</mi>
<mo>,</mo>
<mi>T</mi>
<mi>S</mi>
</mrow>
<mi>&epsiv;</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>T</mi>
</mfrac>
<munderover>
<mo>&Integral;</mo>
<mrow>
<mo>-</mo>
<mi>T</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
<mrow>
<mi>T</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
</munderover>
<mi>S</mi>
<mi>i</mi>
<mi>g</mi>
<mi>m</mi>
<mi>o</mi>
<mi>i</mi>
<mi>d</mi>
<mo>&lsqb;</mo>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>+</mo>
<mi>&tau;</mi>
<mo>/</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<msup>
<mi>Sigmoid</mi>
<mo>*</mo>
</msup>
<mo>&lsqb;</mo>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>&tau;</mi>
<mo>/</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mi>&epsiv;</mi>
<mi>t</mi>
</mrow>
</msup>
<mi>d</mi>
<mi>t</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
1.5)根据广义Sigmoid循环模糊函数公式,计算得到广义Sigmoid循环模糊函数u和f分别代表时延和多普勒频移;
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>C</mi>
<mrow>
<mi>y</mi>
<mi>x</mi>
<mo>,</mo>
<mi>G</mi>
<mi>S</mi>
</mrow>
<mi>&epsiv;</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>,</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mover>
<mo>=</mo>
<mi>&Delta;</mi>
</mover>
<mo>&Integral;</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>S</mi>
<mi>i</mi>
<mi>g</mi>
<mi>mo</mi>
<mi>i</mi>
<mi>d</mi>
</mrow>
<mi>&epsiv;</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>+</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>y</mi>
<mi>x</mi>
<mo>,</mo>
<mi>S</mi>
<mi>i</mi>
<mi>g</mi>
<mi>mo</mi>
<mi>i</mi>
<mi>d</mi>
</mrow>
<mi>&epsiv;</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mo>*</mo>
</msup>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mi>&pi;</mi>
<mi>f</mi>
<mi>&tau;</mi>
</mrow>
</msup>
<mi>d</mi>
<mi>&tau;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mo>&Integral;</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>y</mi>
<mo>,</mo>
<mi>S</mi>
<mi>i</mi>
<mi>g</mi>
<mi>mo</mi>
<mi>i</mi>
<mi>d</mi>
</mrow>
<mi>&epsiv;</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>+</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>y</mi>
<mi>x</mi>
<mo>,</mo>
<mi>S</mi>
<mi>i</mi>
<mi>g</mi>
<mi>mo</mi>
<mi>i</mi>
<mi>d</mi>
</mrow>
<mrow>
<mi>&epsiv;</mi>
<mo>-</mo>
<mi>f</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mo>*</mo>
</msup>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mi>&pi;</mi>
<mi>f</mi>
<mi>&tau;</mi>
</mrow>
</msup>
<mi>d</mi>
<mi>&tau;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mo>&Integral;</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>x</mi>
<mi>y</mi>
<mo>,</mo>
<mi>S</mi>
<mi>i</mi>
<mi>g</mi>
<mi>mo</mi>
<mi>i</mi>
<mi>d</mi>
</mrow>
<mi>&epsiv;</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>+</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>S</mi>
<mi>i</mi>
<mi>g</mi>
<mi>mo</mi>
<mi>i</mi>
<mi>d</mi>
</mrow>
<mrow>
<mi>&epsiv;</mi>
<mo>-</mo>
<mi>f</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mo>*</mo>
</msup>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mi>&pi;</mi>
<mi>f</mi>
<mi>&tau;</mi>
</mrow>
</msup>
<mi>d</mi>
<mi>&tau;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mo>&Integral;</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>x</mi>
<mi>y</mi>
<mo>,</mo>
<mi>S</mi>
<mi>i</mi>
<mi>g</mi>
<mi>mo</mi>
<mi>i</mi>
<mi>d</mi>
</mrow>
<mi>&epsiv;</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>+</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msubsup>
<mi>R</mi>
<mrow>
<mi>y</mi>
<mo>,</mo>
<mi>S</mi>
<mi>i</mi>
<mi>g</mi>
<mi>mo</mi>
<mi>i</mi>
<mi>d</mi>
</mrow>
<mrow>
<mi>&epsiv;</mi>
<mo>-</mo>
<mi>f</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mo>*</mo>
</msup>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mi>&pi;</mi>
<mi>f</mi>
<mi>&tau;</mi>
</mrow>
</msup>
<mi>d</mi>
<mi>&tau;</mi>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
第二步,参数估计
计算公式(10)广义Sigmoid循环模糊函数的绝对值;绝对值最大值的一组时延值u和多普勒频移f,由公式(11)计算待测信号的估计值
<mrow>
<mo>(</mo>
<mover>
<mi>D</mi>
<mo>^</mo>
</mover>
<mo>,</mo>
<msub>
<mover>
<mi>f</mi>
<mo>^</mo>
</mover>
<mi>d</mi>
</msub>
<mo>)</mo>
<mo>=</mo>
<mi>argmax</mi>
<mo>|</mo>
<msubsup>
<mi>C</mi>
<mrow>
<mi>y</mi>
<mi>x</mi>
<mo>,</mo>
<mi>G</mi>
<mi>S</mi>
</mrow>
<mi>&epsiv;</mi>
</msubsup>
<mo>(</mo>
<mi>u</mi>
<mo>,</mo>
<mi>f</mi>
<mo>)</mo>
<mo>|</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
<mo>.</mo>
</mrow>
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610287197.0A CN106027117B (zh) | 2016-05-04 | 2016-05-04 | 一种时延与多普勒频移联合估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610287197.0A CN106027117B (zh) | 2016-05-04 | 2016-05-04 | 一种时延与多普勒频移联合估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106027117A CN106027117A (zh) | 2016-10-12 |
CN106027117B true CN106027117B (zh) | 2018-04-10 |
Family
ID=57082234
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610287197.0A Expired - Fee Related CN106027117B (zh) | 2016-05-04 | 2016-05-04 | 一种时延与多普勒频移联合估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106027117B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107517173B (zh) * | 2017-09-21 | 2019-08-13 | 大连理工大学 | 基于广义循环相关熵的无线电信号载波频率估计方法 |
CN108037494B (zh) * | 2017-12-05 | 2021-05-14 | 大连大学 | 一种脉冲噪声环境下的雷达目标参数估计方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5247303A (en) * | 1992-07-20 | 1993-09-21 | University Corporation For Atmospheric Research | Data quality and ambiguity resolution in a doppler radar system |
CN1339924A (zh) * | 2000-08-21 | 2002-03-13 | 华为技术有限公司 | 自适应格型加权信道估计方法 |
US6366236B1 (en) * | 1999-08-12 | 2002-04-02 | Automotive Systems Laboratory, Inc. | Neural network radar processor |
-
2016
- 2016-05-04 CN CN201610287197.0A patent/CN106027117B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5247303A (en) * | 1992-07-20 | 1993-09-21 | University Corporation For Atmospheric Research | Data quality and ambiguity resolution in a doppler radar system |
US6366236B1 (en) * | 1999-08-12 | 2002-04-02 | Automotive Systems Laboratory, Inc. | Neural network radar processor |
CN1339924A (zh) * | 2000-08-21 | 2002-03-13 | 华为技术有限公司 | 自适应格型加权信道估计方法 |
Non-Patent Citations (3)
Title |
---|
《基于非线性变换的自适应时间延迟估计》;刘文红等;《上海电机学院学报》;20141231;第17卷(第5期);全文 * |
《稳定分布噪声下循环平稳信号时延与多普勒频移估计方法研究》;刘洋;《博士论文》;20131231;全文 * |
《脉冲噪声环境下的SCOT 加权时间延迟估计新方法》;孙永梅等;《通信学报》;20051231;第26卷(第12期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN106027117A (zh) | 2016-10-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103944174B (zh) | 基于互相关函数滤噪算法的低频振荡在线辨识方法 | |
CN106027117B (zh) | 一种时延与多普勒频移联合估计方法 | |
CN104038249B (zh) | 周期长码直扩信号伪随机码估计方法 | |
JP2020519162A (ja) | アップリンク同期タイミング偏差を決定するための方法および装置 | |
CN101977169A (zh) | Ofdm信号的时域参数盲估计方法 | |
CN105954769B (zh) | 一种时延与多普勒频移联合估计方法 | |
JP2012182627A5 (zh) | ||
CN106059973A (zh) | 频偏估计方法和系统 | |
CN103023831A (zh) | 一种适用于突发波形的载波频偏估计方法 | |
CN106878205A (zh) | 一种定时偏差估计方法、装置及终端 | |
CN105516051A (zh) | 数字预失真分数时延估计与信号对齐算法及系统 | |
CN109391573B (zh) | 一种基于线性调频信号的时频二维同步方法 | |
CN102932123B (zh) | 一种探测参考信号的检测方法与系统 | |
CN103152294A (zh) | 基于信号消除进行噪声估计的方法及系统 | |
EP2413552B1 (en) | Method and device for determining smooth window length in channel estimation | |
CN103354538A (zh) | 一种对水声通信中的接收信号进行多普勒补偿的方法 | |
CN101741773A (zh) | 时分同步码分多址接入系统中终端信道估计的方法及装置 | |
CN104184688B (zh) | 一种基于模糊函数的ofdm信号参数估计方法 | |
CN101982945B (zh) | 一种基于混沌理论的频谱感知方法 | |
CN106034094A (zh) | 一种信道估计方法及装置 | |
CN102238117B (zh) | 一种基于信道估计的速度测量方法和设备 | |
CN103825848A (zh) | 一种电力线通信系统中信号干扰噪声比的盲估计方法 | |
CN103634245B (zh) | 一种信号带宽估计方法及装置 | |
CN113660183B (zh) | 一种通信方法及装置 | |
CN110336764B (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180410 |