发明名称 基于平滑伪魏格纳分布的二相编码信号参数估计方法
摘要 本发明涉及一种基于平滑伪魏格纳分布的二相编码信号参数估计方法,主要解决现有方法不能兼顾低信噪比与小计算量的问题。其实现步骤是:1)产生13位巴克码序列编码的二相编码信号;2)对二相编码信号作FFT变换求出其功率谱;3)对功率谱进行平滑处理,并根据平滑后功率谱的3dB带宽功率谱重心来估计信号载频;4)对二相编码信号进行平滑伪魏格纳分布变换,得到载频处的幅度序列;5)搜索幅度序列最小负尖峰间距得到信号码元宽度,并计算相邻负尖峰所包含的子码个数;6)根据负尖峰位置以及子码个数,得到码元序列的估计。本发明在低信噪比情况下提高了对二相编码信号的参数的估计精度,减小了计算量,可用于对目标的准确侦察。
申请公布号 CN106501787A 申请公布日期 2017.03.15
申请号 CN201610952840.7 申请日期 2016.11.02
申请人 西安电子科技大学 发明人 张林让;蒙妍;许建华;彭炜杰;赵永志
分类号 G01S7/41(2006.01)I 主分类号 G01S7/41(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 基于平滑伪魏格纳分布的二相编码信号参数估计方法,包括:(1)产生由13位巴克码序列[1 1 1 1 1 ‑1 ‑1 1 1 ‑1 1 ‑1 1]编码的二相编码信号s(t);(2)对二相编码信号s(t)进行FFT变换,得到其功率谱Y(m);(3)对功率谱Y(m)进行平滑处理,得到平滑后的功率谱Y<sub>s</sub>(k);(4)根据平滑后功率谱Y<sub>s</sub>(k)的3dB带宽功率谱重心,估计出信号载频<img file="FDA0001140109430000011.GIF" wi="67" he="71" />(5)获得幅度序列<img file="FDA0001140109430000012.GIF" wi="315" he="78" />5a)对二相编码信号s(t)进行平滑伪魏格纳分布SPWVD变换,得到二维幅频矩阵SPWVD(t,f),5b)对二维幅频矩阵SPWVD(t,f)在频率维进行搜索,得到信号载频<img file="FDA0001140109430000013.GIF" wi="45" he="71" />处的幅度序列<img file="FDA0001140109430000014.GIF" wi="315" he="77" />(6)计算二相编码信号s(t)的码元宽度<img file="FDA0001140109430000015.GIF" wi="67" he="71" />6a)搜索幅度序列<img file="FDA0001140109430000016.GIF" wi="308" he="78" />得到各个负尖峰位置,计算出相邻负尖峰间隔Δ;6b)取最小相邻负尖峰间隔Δ记为Δ<sub>min</sub>,则信号码元宽度的估计值<img file="FDA0001140109430000017.GIF" wi="42" he="71" />为:<maths num="0001"><math><![CDATA[<mrow><msub><mover><mi>T</mi><mo>^</mo></mover><mi>b</mi></msub><mo>=</mo><msub><mi>&Delta;</mi><mi>min</mi></msub><mo>&CenterDot;</mo><msub><mi>T</mi><mi>s</mi></msub><mo>,</mo></mrow>]]></math><img file="FDA0001140109430000018.GIF" wi="262" he="70" /></maths>其中,T<sub>s</sub>为信号采样周期;(7)计算由负尖峰分隔开的每段数据所包含的子码个数l:l=Δ/Δ<sub>min</sub>;(8)计算二相编码信号的码元序列<img file="FDA0001140109430000019.GIF" wi="83" he="79" />8a)根据负尖峰位置以及由负尖峰分隔开的每段数据所包含的子码个数l,粗略估计出码元序列<img file="FDA00011401094300000110.GIF" wi="81" he="71" />8b)对粗略估计出的码元序列<img file="FDA00011401094300000111.GIF" wi="81" he="71" />按照如下公式器得到码元序列的精确估计<img file="FDA00011401094300000112.GIF" wi="81" he="69" /><maths num="0002"><math><![CDATA[<mrow><msub><mover><mi>P</mi><mo>^</mo></mover><mi>m</mi></msub><mo>=</mo><mi>arg</mi><mi>min</mi><mo>{</mo><mo>|</mo><mo>|</mo><mover><mi>s</mi><mo>^</mo></mover><mrow><mo>(</mo><msub><mover><mi>f</mi><mo>^</mo></mover><mn>0</mn></msub><mo>,</mo><msub><mover><mi>T</mi><mo>^</mo></mover><mi>b</mi></msub><mo>,</mo><msub><mover><mi>P</mi><mo>~</mo></mover><mi>m</mi></msub><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mi>s</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>|</mo><mo>|</mo><mo>}</mo><mo>,</mo></mrow>]]></math><img file="FDA00011401094300000113.GIF" wi="886" he="110" /></maths>其中,||·||是取范数运算,<img file="FDA00011401094300000114.GIF" wi="310" he="87" />是由估计出的信号载频<img file="FDA00011401094300000115.GIF" wi="75" he="77" />码元宽度<img file="FDA00011401094300000116.GIF" wi="74" he="71" />码元序列<img file="FDA00011401094300000117.GIF" wi="53" he="70" />以及时间t重新构造出的二相编码信号,arg min是取最优解。
地址 710071 陕西省西安市雁塔区太白南路2号
您可能感兴趣的专利