申请/专利权人:南京航空航天大学
申请日:2021-05-20
公开(公告)日:2022-07-26
公开(公告)号:CN113376605B
主分类号:G01S7/41
分类号:G01S7/41
优先权:
专利状态码:有效-授权
法律状态:2022.07.26#授权;2021.09.28#实质审查的生效;2021.09.10#公开
摘要:本发明公开了一种基于原子范数最小化的mmwave‑MIMO雷达DOA估计方法,首先根据mmwave‑MIMO雷达发射天线个数M,接收天线个数N,信源个数J,获取p×q×M×N的三维数据包D;其次,对三维数据包D先在对距离门做q点的傅里叶变换,给出目标所在距离门单元;依次从D中取出目标所在距离门单元下的二维数组Bi,将Bi按子脉冲进行脉冲累积得到单测量矢量bi;将单测量矢量bi等价于阵元数为M×N的均匀线阵下的DOA估计;给出基于时分mmwave‑MIMO雷达的超分辨DOA估计函数表达式,并转换成凸函数形式,最后求解原子范数最小化问题;得到目标角度参数及其导向矢量与噪声子空间的正交性。
主权项:1.一种基于原子范数最小化的mmwave-MIMO雷达DOA估计方法,其特征在于,包括以下步骤:S1、假设已知mmwave-MIMO雷达发射天线个数为M,接收天线个数为N,信源个数为J,其中JM×N;对发射天线端采用时分方式发射正交波形,依次从第1根发射天线开始到第M根发射天线结束,每间隔Tc秒发射线性调频信号,Tc为一个脉冲周期,在一个脉冲周期内,发射p个子脉冲进行脉冲积累;并且对子脉冲做距离门上的q次离散采样,最终得到p×q×M×N的三维数据包D;S2:对三维数据包D先取出距离门数据d=Dp,:,M×N,对d做q点的傅里叶变换,找出相对应的I个峰值点,即为目标所在距离门单元[q1,q2,…,qI],I≤J;不妨令循环变量i=1,…,J,qi为第i个目标所在距离门单元;若每个目标均对应不同的距离门单元,则I=J;反之,IJ;随后开始循环,初始化i=1;S3、判断循环变量i是否小于等于I,若是,则执行步骤S4;否则,则退出循环,执行步骤S9;S4、从三维数据包D中取出目标所在距离门单元下的二维数组Bi,其中Bi=D1:p,qi,1:M×N,将二维数组Bi按子脉冲进行脉冲累积得到单测量矢量bi;S5、将单测量矢量bi等价于阵元数为M×N的均匀线阵下的DOA估计;S6、根据原子范数理论,首先给出基于时分mmwave-MIMO雷达的超分辨DOA估计函数表达式,其次根据半正定规划理论将函数表达式转换成凸函数形式,最后通过CVX凸优化工具包来求解原子范数最小化问题;S7、对根据CVX凸优化工具包所求解得出的凸优化转换生成的中间矩阵Tμ的估计进行特征值分解EVD,最终通过谱峰搜索得到目标角度参数及其导向矢量与噪声子空间的正交性;S8、令i=i+1,返回执行步骤S3;S9、结束循环,此时已经遍历目标所在距离门,整合得出Θ=[Θ1,...,ΘI]及其每个距离门下的Ppower。
全文数据:
权利要求:
百度查询: 南京航空航天大学 一种基于原子范数最小化的mmwave-MIMO雷达超分辨DOA估计方法
免责声明
1、本报告根据公开、合法渠道获得相关数据和信息,力求客观、公正,但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解,仅供参考使用,不能作为本公司承担任何法律责任的依据或者凭证。