一种基于单快拍MUSIC算法的距离向处理方法技术

本发明专利技术公开了一种基于单快拍MUSIC算法的距离向处理方法，它是基于运用单快拍MUSIC算法估计目标的位置，首先对原始回波信号进行去斜处理，然后根据单快拍MUSIC算法进行构造自相关矩阵、特征分解等过程；计算MUSIC谱并进行谱峰搜索，这样就能确定目标所在的距离单元格；最后使用最小二乘法计算目标所在距离单元格的值。本发明专利技术与传统的通过脉冲压缩方法进行距离向处理相比，本发明专利技术仅需要单个回波信号即可使用MUSIC算法进行距离向的处理，分辨率有较大的提高，能够实现距离向超分辨能力，本发明专利技术特别适用于目标较少，精度要求高的情况。

【技术实现步骤摘要】
一种基于单快拍MUSIC算法的距离向处理方法
本技术专利技术属于雷达
，它特别涉及了合成孔径雷达(SAR)成像

技术介绍
合成孔径雷达(SAR)是一种高分辨率微波成像雷达，具有全天时和全天候工作的优点，已被广泛应用各个领域，如地形测绘、制导、环境遥感和资源勘探等。SAR应用的重要前提和信号处理的主要目标是通过成像算法获取高分辨、高精度的微波图像。其获得的高分辨微波图像已经被广泛应用于诸多领域，如生成数字高程图、观测火山活动和洪灾情况、监测陆地和海洋交通等。距离向处理是SAR成像过程中对信号进行的第一步处理方式。通过方位向处理，获得目标的一维图像。距离向处理传统方法有匹配滤波、去斜处理等方法，详见“皮亦鸣.杨建宇等.合成孔径雷达成像原理[M].电子科技大学出版社,2007”。阵列信号处理(ASP)作为信号处理的一个重要分支是在空间不同位置上以一定的规律布觉传感器，形成传感器阵列。阵列的配置由两部分内容组成：一是每个阵元的天线方向图；二是阵列的几何结构(即阵元的物理位置)，通常我们将其分成线性阵列、平面阵列和立体阵列三类。阵列信号处理被广泛应用于雷达、通信、声纳、测向、地震探测、射电天文和电子医疗工程等领域。MUSIC算法是一种基于矩阵特征空间分解的方法。信号处理的观测空间分解为信号子空间和噪声子空间，显然这两个空间是正交的。信号子空间由阵列接收到的数据协方差矩阵中与信号对应的特征向量组成，噪声子空间则由协方差矩阵中所有最小特征值(噪声方差)对应的特征向量组成。通常使用MUSIC算法进行波达方向(DOA)估计。目前，用于SAR成像的MUSIC算法有部分研究，有的是需要先聚焦成像，再使用MUSIC算法提高分辨率，详见“张平,商建,杨汝良.一种有效的二维MUSIC超分辨SAR成像算法[J].系统仿真学报,2010,22(1):184-187”。有的是需要多个快拍，即多个距离向回波信号进行联合处理，且只计算出目标的位置，没有计算出信号的回波幅度，详见“张东浩.线阵三维SAR成像算法研究及仿真[D].电子科技大学,2010”。有的是计算出目标的位置及幅度，但是需要多个快拍才能完成计算，详见“成晨.三维合成孔径雷达超分辨成像方法研究[D].电子科技大学,2012”。
技术实现思路
本专利技术提出了一种基于单快拍MUSIC算法的距离向处理方法，它主要运用单快拍MUSIC算法估计目标的位置，找到目标位置以后，再通过最小二乘法估计目标的幅度。通过这两步处理，就完成了SAR成像过程中的距离向处理。为了方便描述本专利技术的内容，首先作以下术语定义：定义1、合成孔径雷达(SAR)合成孔径雷达是将雷达是将雷达固定于载荷运动平台上，结合运动平台的运动以合成线性阵列以达到运动向的分辨率，再利用雷达波束向回波延时实现距离一维成像，从而实现观测目标二维成像的一种合成孔径雷达技术。定义2、标准的MUSIC算法标准的MUSIC算法是一种基于矩阵特征空间分解的方法。信号处理的观测空间可以分解为信号子空间和噪声子空间，显然这两个空间是正交的。信号子空间由阵列接收到的数据协方差矩阵中与信号对应的特征向量组成，噪声子空间则由协方差矩阵中所有最小特征值(噪声方差)对应的特征向量组成。通过MUSIC算法得到的频谱叫做MUSIC谱。详见“何子述.现代数字信号处理及其应用[M].清华大学出版社,2009”。定义3、标准合成孔径雷达原始回波仿真方法标准合成孔径雷达原始回波仿真方法是指给定雷达系统参数、平台轨迹参数以及观测场景参数等所需的参数条件下，基于合成孔径雷达成像原理仿真得到具有SAR回波信号特性的原始回波信号的方法，详细内容可参考文献“LanG.CummingFrank,HWong.合成孔径雷达成像:算法与实现[M].电子工业出版社,2012”。将垂直于雷达平台运动的方向叫做距离向。对距离向进行离散化处理，离散化之后的距离向由多个小格组成，每个小格被称为距离单元。定义4、单快拍快拍数一般是指时域上的采样点数。单快拍是指全部阵元采样一次。协方差矩阵受到快拍数影响，如果快拍数少，则协方差矩阵估计不准确，那么噪声子空间不准，最终将导致MUSIC算法性能下降。定义5、标准的特征分解方法特征分解(Eigendecomposition)又被称为谱分解(Spectraldecomposition)。它是将矩阵分解为由其特征值和特征向量表示的矩阵之积的方法。详见“矩阵理论”，黄廷祝编，高等教育出版社。定义6、标准的去斜处理方法去斜处理就是将输入的线性调频信号与一个调频率大小相同但符号相反的参考信号相乘，去掉2次相位，得到频率与信号输入位置成正比的单频信号。由于单频信号的频率与输入信号位置相关，可通过FFT变换实现对线性调频信号的压缩。定义7、最小二乘法最小二乘法是一种数学优化技术。它通过最小化误差的平方和寻找数据的最佳函数匹配。利用最小二乘法可以简便地求得未知的数据，并使得这些求得的数据与实际数据之间误差的平方和为最小。最小二乘法还可用于曲线拟合。其他一些优化问题也可通过最小化能量或最大化熵用最小二乘法来表达。详见“SorensonHW.Least-squaresestimation:fromGausstoKalman[J].IEEEspectrum,1970,7(7):63-68”。定义8、标准的施密特正交化方法施密特正交化(Schmidtorthogonalization)是求欧氏空间正交基的一种方法。在线性代数中，如果内积空间上的一组向量能够组成一个子空间，那么这一组向量就称为这个子空间的一个基。施密特正交化提供了一种方法，能够通过这一子空间上的一个基得出子空间的一个正交基，并可进一步求出对应的标准正交基。详见“线性代数”，同济大学数学系编，高等教育出版社。定义9、共轭转置共轭转置即对复数矩阵A进行转置并取共轭，记为AH，可通过标准的共轭转置方法计算得出，。详见“线性代数”，同济大学数学系编，高等教育出版社。定义10、标准的矩阵求逆方法假设矩阵A和矩阵B，如果AB＝E，其中E为单位矩阵，则称矩阵B是矩阵A的右逆矩阵，通常将矩阵B写为A-1，根据矩阵A，通过标准的矩阵求逆方法计算出矩阵A-1。详见“矩阵理论”，黄廷祝，钟守铭等，高等教育出版社。定义11、共轭复数共轭复数，两个实部相等，虚部互为相反数的复数互为共轭复数。详见“线性代数”，同济大学数学系编，高等教育出版社。定义12、标准的谱峰搜索方法谱峰搜索即找出信号频谱中谱峰的位置，可以将谱峰视为极大值点，在谱峰搜索的过程中，通常将设置一个门限值，小于门限值的频谱会被忽略。MATLAB中已有谱峰搜索函数，即’findpeaks’，可以直接调用，实现对MUSIC谱的谱峰搜索，本文中的谱峰搜索方法即基于MATLAB中的’findpeaks’函数。详见MATLAB中的帮助文档。定义13、标准合成孔径雷达原始回波仿真方法标准合成孔径雷达原始回波仿真方法是指给定雷达系统参数、平台轨迹参数以及观测场景参数等所需的参数条件下，基于合成孔径雷达成像原理仿真得到具有SAR回波信号特性的原始回波信号的方法，详细内容可参考文献“LanG.CummingFrank,HWong.合成孔径雷达成像:算法与实现[M].电子工业出版本文档来自技高网...

【技术保护点】
一种基于单快拍MUSIC算法的距离向处理方法，其特征是它包括以下几个步骤：步骤1、初始化SAR系统参数：初始化SAR系统参数包括：雷达工作中心频率，记做fc；雷达载频波长，记做λ；雷达发射基带信号的信号带宽，记做Br；雷达发射信号脉冲宽度，记做Tr；雷达发射信号的调频斜率，记做Kr；雷达接收系统的采样频率，记做fs；电磁波在空气中的传播速度，记做C；距离向采样时间为TP，n为距离向时间序号，n＝1，2，…，Nr，n为自然数，Nr为距离向采样点总数；上述参数均为SAR系统标准参数，其中雷达中心频率fc，雷达载频波长λ，雷达发射基带信号的信号带宽Br，雷达发射信号脉冲宽度Tr，雷达发射信号调频斜率Kr，雷达接收系统的采样频率fs，波门延时为τ0，在SAR系统设计过程中已经确定；根据SAR成像系统方案和观测方案，距离向处理所需要的初始化成像系统参数均为已知；步骤2、初始化SAR回波信号：在第t个距离向采样得到的原始回波数据，记为SaE(t)，t＝1，2，…，TP，t为自然数；在实际线阵SAR系统中，原始回波数据SaE(t)可由雷达系统数据接收机提供；在仿真线阵SAR成像过程中，原始回波数据SaE(t)根据雷达系统参数，采用标准合成孔径雷达原始回波仿真方法产生得到；在雷达回波数据进行距离向处理之前，原始回波数据SaE(t)均已知；初始化场景大小为[‑X，X]，将原点设置场景中心，定义场景中有M个散射点，且目标都在距离向接收范围内，这M个散射点相对于场景中心位置为x1，x2，...，xM，这M个点的散射系数分别记为...

【技术特征摘要】
1.一种基于单快拍MUSIC算法的距离向处理方法，其特征是它包括以下几个步骤：步骤1、初始化SAR系统参数：初始化SAR系统参数包括：雷达工作中心频率，记做fc；雷达载频波长，记做λ；雷达发射基带信号的信号带宽，记做Br；雷达发射信号脉冲宽度，记做Tr；雷达发射信号的调频斜率，记做Kr；雷达接收系统的采样频率，记做fs；电磁波在空气中的传播速度，记做C；距离向采样时间为TP，n为距离向时间序号，n＝1，2，…，Nr，n为自然数，Nr为距离向采样点总数；上述参数均为SAR系统标准参数，其中雷达中心频率fc，雷达载频波长λ，雷达发射基带信号的信号带宽Br，雷达发射信号脉冲宽度Tr，雷达发射信号调频斜率Kr，雷达接收系统的采样频率fs，波门延时为τ0，在SAR系统设计过程中已经确定；根据SAR成像系统方案和观测方案，距离向处理所需要的初始化成像系统参数均为已知；步骤2、初始化SAR回波信号：在第t个距离向采样得到的原始回波数据，记为SaE(t)，t＝1，2，…，TP，t为自然数；在实际线阵SAR系统中，原始回波数据SaE(t)可由雷达系统数据接收机提供；在仿真线阵SAR成像过程中，原始回波数据SaE(t)根据雷达系统参数，采用标准合成孔径雷达原始回波仿真方法产生得到；在雷达回波数据进行距离向处理之前，原始回波数据SaE(t)均已知；初始化场景大小为[-X，X]，将原点设置场景中心，定义场景中有M个散射点，且目标都在距离向接收范围内，这M个散射点相对于场景中心位置为x1，x2，...，xM，这M个点的散射系数分别记为通过公式得到M个散射点所对应的回波时延，分别为其中，x为目标相对于场景中心的距离，C为光速，原始回波数据SaE(t)通过采用标准合成孔径雷达原始回波仿真方法产生得到，即其中，A为回波信号的增益，这里将A设为1，根据采样频率fs对连续回波信号SaE(t)进行采样，得到离散回波信号SAE(n)，即其中，Nr为距离向采样点数；采用标准的去斜处理方法，得到参考信号其中，n＝1，2，…，Nr，n为自然数；上式中所使用的参数均为步骤1中初始化的参数；步骤3、对原始回波数据进行去斜处理：通过标准的取共轭方法对参考信号href(n)进行处理，得到共轭参考信号通过将步骤2中的回波信号SAE(n)与共轭参考信号相乘，实现去除回波信号相位中二次项，得到处理后的信号SE(n)，这个过程用公式表示为其中，n＝1，2，…，Nr，n为自然数，为M个点的散射系数，为M个散射点所对应的回波时延分，Kr为雷达发射信号的调频斜率，fs为雷达接收系统的采样频率，fc为雷达工作中心频率，SAE(n)为步骤2中得到的离散回波信号；步骤4、生成噪声子空间矩阵：步骤4.1、构造自相关矩阵：步骤4.1.1初始化矩阵及信号向量：初始化一个Nr×Nr维零矩阵Ax，矩阵Ax的结构如下：矩阵Ax中的元素ax，y均为0，其中，x，y为自然数，x＝1，2，…，Nr，y＝1，2，…，Nr；回波信号SE(n)表示为向量形式SE＝[SE(1)，SE(2)，…，SE(Nr)]，其中，SE(n)为步骤3中得到的n时刻的回波信号，n＝1，2，…，Nr；步骤4.1.2、对矩阵Ax第m行处理的第1步：选择向量SE＝[SE(1)，SE(2)，…，SE(Nr)]中第1个到第Nr-m+1个元素，这Nr-m+1个元素定义为向量SE(m，1)＝[SE(1)，SE(2)，…，SE(Nr-m+1)]，然后将向量SE(m，1)填入矩阵Ax中的第m行，第m列到第Nr列的元素位置，得到矩阵Ax(m，1)，其中，Ax是步骤4.1.1中定义的矩阵，矩阵Ax(m，1)的结果如下其中，SE(n)为步骤4.1.1中向量的元素SE＝[SE(1)，SE(2)，…，SE(Nr)]，n＝1，2，…，Nr-m+1；ax，y为矩阵Ax中的元素，其中，x，y为自然数，x＝1，2，…，Nr，y＝1，2，…，Nr；步骤4.1.3、对矩阵Ax第m行处理的第2步：选择向量SE＝[SE(1)，SE(2)，…，SE(Nr)]中第2个到第m个元素，这m-1个元素定义为向量SE(m，2)＝[SE(2)，SE(3)，…，SE(m)]，然后将向量SE(m，2)中元素进行反向排序得到向量SE(m，-2)＝[SE(m)，SE(m-1)，…，SE(2)]，采用标准的取共轭的方法对向量SE(m，-2)进行处理得到共轭向量SE*(m，-2)＝[SE*(m)，SE*(m-1)，…，SE*(2)]，最后，将向量SE*(m，-2)填入矩阵Ax(m，1)的第m行，第1列到第m列的位置，得到矩阵Ax(m，2)，其中，Ax(m，1)是步骤4.1.2中的矩阵，矩阵Ax(m，2)的结果如下

【专利技术属性】
技术研发人员：张晓玲，李良，周灵杰，韦顺军，师君，
申请(专利权)人：电子科技大学，
类型：发明
国别省市：四川,51