发明名称 一种基于线性规划的近场目标定位方法及其系统
摘要 本发明提供一种基于线性规划的近场目标定位方法及其系统。利用压缩感知理论实现近场目标定位,在未知信号形式的情况下,选定参考阵元并使其工作在正常的采样频率下,其他阵元工作在远低于奈奎斯特采样频率下;将参考阵元输出信号作为参考目标信号,得到所有阵元的采样数据;生成延时表和稀疏基矩阵;产生随机映射矩阵,并得到系数矩阵;线性规划求解,获取稀疏向量;基于获得的稀疏向量估计结果,从预先设定的方位-距离集得到近场目标定位信息。根据本发明,传感器无需工作在奈奎斯特采样频率以上,可大幅降低采样率,节省传感器工作能耗,提高系统资源利用率;对目标带宽没有限制,可同时适用于窄、宽带目标定位,也适用于非高斯型目标的定位。
申请公布号 CN101644774B 申请公布日期 2012.08.15
申请号 CN200910088660.9 申请日期 2009.07.07
申请人 中国科学院声学研究所 发明人 田彪;李宇;黄海宁
分类号 G01S15/06(2006.01)I;G01S15/08(2006.01)I 主分类号 G01S15/06(2006.01)I
代理机构 北京法思腾知识产权代理有限公司 11318 代理人 杨小蓉
主权项 1.一种基于线性规划的近场目标定位方法,其特征在于,包括以下步骤:1)将含L个阵元线列阵中的某个阵元选定为参考阵元,并以参考阵元输出信号作为参考目标信号ζ<sub>rf</sub>(t),在未知信号形式的情况下,使参考阵元工作在依据奈奎斯特采样定律确定的正常的采样频率下,其他阵元工作在远低于奈奎斯特采样频率下,此时,参考阵元工作的采样频率为F<sub>s</sub>,该参考阵元输出信号作为参考目标信号ζ<sub>rf</sub>(t),采样点数设为K,其它阵元的采样频率为F<sub>cs</sub>,采样点数为M;在已知目标信号形式的情况下,包括参考阵元在内所有的阵元均工作在远低于奈奎斯特采样频率下,采样频率为F<sub>cs</sub>,此时,直接选定一个已知的目标信号作为参考目标信号;2)基于参考目标信号,得到所有阵元的采样数据β<sub>i</sub>,并且<img file="FDA0000149951970000011.GIF" wi="376" he="80" />LM×1维,上标T表示转置,这里,在未知目标信号形式的情况下,由于所述参考阵元的采样频率F<sub>s</sub>不同于压缩采样频率F<sub>cs</sub>,需要对该参考阵元的数据进行抽取,以使得其采样点数与其他阵元的采样点数相同;3)生成延时表和稀疏基矩阵Ψ,包括如下子步骤:31)对每一个阵元产生相应的稀疏基矩阵Ψ<sub>i</sub>;<img file="FDA0000149951970000012.GIF" wi="1500" he="301" />其中,i=1,2,…,L;n=1,2,…,N;p=1,2,…,P,式中,N表示探测的角度数,P表示探测的距离数,α<sub>n,p</sub>=[θ<sub>n</sub>,R<sub>p</sub>]表示延时参数对,θ<sub>n</sub>代表方位,R<sub>p</sub>代表距离,Δ<sub>i</sub>(α<sub>n,p</sub>)为其它阵元与参考阵元间的延时;32)将各阵元的稀疏基矩阵组合成阵列的稀疏基矩阵:<maths num="0001"><![CDATA[<math><mrow><mi>&Psi;</mi><mo>=</mo><msup><mrow><mo>[</mo><msubsup><mi>&Psi;</mi><mn>1</mn><mi>T</mi></msubsup><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msubsup><mi>&Psi;</mi><mi>L</mi><mi>T</mi></msubsup><mo>]</mo></mrow><mi>T</mi></msup><mo>,</mo></mrow></math>]]></maths>LNP×K维:4)产生随机映射矩阵Φ,并得到系数矩阵A=ΦΨ,包括如下子步骤:41)对每一个阵元产生相应的随机映射矩阵φ<sub>i</sub>:φ<sub>i</sub>为M×NP维随机产生矩阵;42)产生相应的随机映射矩阵Φ;Φ为LM×LNP维块对角矩阵,φ<sub>i</sub>为其块对角元素;43)依据随机映射矩阵Φ和稀疏基矩阵Ψ,计算系数矩阵A=ΦΨ;5)线性规划求解,获取稀疏向量b:<img file="FDA0000149951970000021.GIF" wi="316" he="71" />当||A<sup>T</sup>(β-Ab)||<sub>∞</sub><ε<sub>N</sub>σ时,得到稀疏向量b的估计,其中,参数ε<sub>N</sub>表示误差的容忍度,取值为<img file="FDA0000149951970000022.GIF" wi="194" he="56" />稀疏向量b=[b<sub>1,1</sub>,b<sub>1,2</sub>,…,b<sub>N,P</sub>]<sup>T</sup>;6)基于所获得的稀疏向量b的估计结果,得到近场目标定位信息:稀疏向量b中的非零正元素b<sub>n,p</sub>对应方位-距离对α<sub>n,p</sub>=[θ<sub>n</sub>,R<sub>p</sub>],从而从预先设定的方位-距离集中得到近场目标的位置信息:方位θ<sub>n</sub>和距离R<sub>p</sub>。
地址 100190 北京市海淀区北四环西路21号中国科学院声学研究所