发明名称 一种基于干扰抵消的短波多径信号同步方法
摘要 本发明提出了一种基于干扰抵消的短波多径信号同步方法,所述方法针对多径区分不明显的DS-SS系统,对基带信号码片匹配滤波并过抽样,与本地扩频码相关后通过串并变换、幅值平方、平滑,得到多径识别判决变量和最大值位置,进而得到路径同步位置与幅度,再对当前路径信号进行重构并抵消对后续路径的干扰,估计后续路径的位置,进而确定多径时延。本发明方法估计了短波多径信号的路径位置和路径幅度,重构多径信号,抵消了对后续路径的影响,重新识别信号和多径时延,从而实现多径信号的同步;在虚警概率基本不变的情况下增加对后续路径信号识别的概率,提高多径时延估计的准确性,提高了多径信号的检测概率。
申请公布号 CN103501187B 申请公布日期 2015.06.03
申请号 CN201310470951.0 申请日期 2013.10.10
申请人 中国人民解放军理工大学 发明人 郑学强;程云鹏;沈良;盛雁鸣;张玉明
分类号 H04B1/7115(2011.01)I;H04B1/7113(2011.01)I 主分类号 H04B1/7115(2011.01)I
代理机构 江苏永衡昭辉律师事务所 32250 代理人 王斌
主权项 一种基于干扰抵消的短波多径信号同步方法,其特征在于,其步骤如下:步骤1,将基带信号r(t)通过码片匹配滤波器g<sub>R</sub>(t),输出得到信号<img file="FDA0000647484620000011.GIF" wi="124" he="85" />信号<img file="FDA0000647484620000012.GIF" wi="100" he="86" />的表达式为:<img file="FDA0000647484620000013.GIF" wi="397" he="86" />其中,<img file="FDA0000647484620000014.GIF" wi="54" he="45" />表示卷积,对<img file="FDA0000647484620000015.GIF" wi="102" he="88" />以抽样速率Q/T<sub>c</sub>进行抽样得到信号<img file="FDA0000647484620000016.GIF" wi="135" he="87" />且令<img file="FDA0000647484620000017.GIF" wi="303" he="82" />Q是过抽样因子,T<sub>c</sub>是码片周期;将多径信号的路径号k初始化,并置为1;步骤2,用本地扩频码c=[c(0),...c(l),...c(L‑1)]对y<sub>k</sub>(n)进行相关运算,得到相关值e(m):<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>e</mi><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>l</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>L</mi><mo>-</mo><mn>1</mn></mrow></munderover><msup><mi>c</mi><mo>*</mo></msup><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow><msub><mi>y</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>lQ</mi><mo>+</mo><mi>m</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000647484620000018.GIF" wi="565" he="146" /></maths>其中,c<sup>*</sup>(l)为c(l)的复数共轭;l表示扩频码码片序号,且0≤l≤L‑1,L为扩频码长度;m=0,…,LQW‑1,W是平滑因子;步骤3,将相关值e(m)进行平方、平滑累加得到z(h),表示为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>z</mi><mrow><mo>(</mo><mi>h</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>w</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>W</mi><mo>-</mo><mn>1</mn></mrow></munderover><msup><mrow><mo>|</mo><mi>e</mi><mrow><mo>(</mo><mi>wLQ</mi><mo>+</mo><mi>h</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup></mrow>]]></math><img file="FDA0000647484620000019.GIF" wi="494" he="148" /></maths>其中h=0,…,LQ‑1;步骤4,计算z(h)的最大值z<sub>max</sub>(h)及对应位置<img file="FDA00006474846200000110.GIF" wi="131" he="83" /><img file="FDA00006474846200000111.GIF" wi="447" he="89" />z<sub>max</sub>(h)表达式如下:z<sub>max</sub>(h)=max{z(0),z(1),…,z(LQ‑1)}将z<sub>max</sub>(h)及其两侧各M个点剔除,M为自然数,计算剔除后的平均值<img file="FDA00006474846200000112.GIF" wi="145" he="82" />表示为:<img file="FDA0000647484620000021.GIF" wi="1511" he="518" />将最大值z<sub>max</sub>(h)与平均值<img file="FDA0000647484620000022.GIF" wi="116" he="87" />的比值作为判决变量<img file="FDA0000647484620000023.GIF" wi="465" he="89" />然后将判决变量λ(h)与预定的路径存在判决门限λ<sub>0</sub>进行比较:如果λ(h)&lt;λ<sub>0</sub>,说明当前路径不存在;若k=1,则返回步骤1,继续执行同步过程;否则,则判定已检测了全部k‑1条路径,退出同步过程;如果λ(h)≥λ<sub>0</sub>,说明当前路径存在,执行步骤5;步骤5,将最大值位置<img file="FDA0000647484620000024.GIF" wi="106" he="85" />作为第k条路径的捕捉位置,即<img file="FDA0000647484620000025.GIF" wi="317" he="92" />计算第k条路径幅度估计值α<sub>k</sub>:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>&alpha;</mi><mi>k</mi></msub><mo>=</mo><msqrt><msub><mi>z</mi><mi>max</mi></msub><mrow><mo>(</mo><mi>h</mi><mo>)</mo></mrow><mo>/</mo><mrow><mo>(</mo><msup><mi>WL</mi><mn>2</mn></msup><mo>)</mo></mrow></msqrt></mrow>]]></math><img file="FDA0000647484620000026.GIF" wi="455" he="110" /></maths>步骤6,利用第k条路径位置δ<sub>k</sub>(h)及幅度估计值α<sub>k</sub>,重构第k条路径信号x<sub>k</sub>(n),表示为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>x</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>&alpha;</mi><mi>k</mi></msub><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>WL</mi><mo>-</mo><mn>1</mn></mrow></munderover><mo>[</mo><mover><mi>c</mi><mo>~</mo></mover><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mi>g</mi><mrow><mo>(</mo><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mi>iQ</mi><mo>-</mo><msub><mi>&delta;</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>h</mi><mo>)</mo></mrow><mo>)</mo></mrow><mfrac><msub><mi>T</mi><mi>c</mi></msub><mi>Q</mi></mfrac><mo>)</mo></mrow><mo>]</mo></mrow>]]></math><img file="FDA0000647484620000027.GIF" wi="904" he="164" /></maths>其中,<img file="FDA0000647484620000028.GIF" wi="49" he="66" />是W个连续本地扩频码的级联,g(t)是码片发送滤波器g<sub>T</sub>(t)与匹配滤波器g<sub>R</sub>(t)的卷积,接收端已知g<sub>T</sub>(t),且g<sub>T</sub>(t)=g<sub>R</sub>(t);步骤7,将第k条路径对后续路径的干扰消除,即从当前信号y<sub>k</sub>(n)中减去重构的第k径信号x<sub>k</sub>(n),得到干扰消除后的信号y<sub>k+1</sub>(n),用于检测识别第k+1条路径,表示为:y<sub>k+1</sub>(n)=y<sub>k</sub>(n)‑x<sub>k</sub>(n)k值增加1后,返回步骤2继续执行。
地址 210007 江苏省南京市御道街标营2号