专利名称：：一种适用于稀疏微波成像的随机噪声雷达的信号处理方法
技术领域：
：本发明涉及微波成像
技术领域：
，特别是一种适用于稀疏微波成像的随机噪声雷达的信号处理方法。
背景技术：
：传统的高分辨率检测雷达、合成孔径雷达和逆合成孔径雷达等是通过发射天线对被观测区域发射电磁脉冲，并利用接收天线接收雷达的回波。回波通常被转化成数字信号并存储起来，计算机通过对这些数据的处理实现目标检测和场景成像。在将回波的模拟信号转化成数字信号的时候必须采用模数转换器(ADC)。这个采样过程可以发生在匹配滤波之前或者匹配滤波之后(MerrillI.Skolnik,"RadarHandbook(SecondEdition)”，McGraw-Hill,1990；andN.LevanonandΕ.Mozeson,"RadarSignals，，，2004)。*艮据雷达分辨率的公式，分辨率和发射信号的带宽成反比。为了获得高分辨的雷达检测和成像结果，必须发射带宽很宽的雷达信号。而根据经典的奈奎斯特采样定理，模拟信号到数字信号的采样率(ADC速率)至少必须为带宽的两倍(对于复信号，采样率至少应该等于信号带宽)。这对ADC提出了很高的要求，而且对数据的存储、处理和传输也都带来很大的负担。传统的雷达检测和成像方法是基于匹配滤波实现的，即通过求发射信号和接收信号自相关实现的。自相关结果在有回波的时刻会出现理想的冲击脉冲。当ADC速率低于奈奎斯特采样率时，采用传统方法的成像结果会出现模糊。为了提高雷达的截获概率性能和电子反对抗能力，二十世纪60年代美国和欧洲就对以微波噪声源作为其发射信号或者调制形式的随机信号雷达给予了广泛的关注，60年代末期。美国的Purdue大学最早研制了一部试验型随机信号雷达。法国的Carpenter教授讨论了米波波段的相关法随机信号雷达。70年代中期美国Minnesota大学的Kaveh教授做了交错脉冲随机信号雷达的研究。同时，英国London大学的Rarest等人做了固态随机信号雷达的研究。80年代以来，电子技术的发展给随机信号雷达的运用提供了可能，已有随机信号雷达运用于地质探测和微波成像。南京理工大学的刘国岁教授从70年代开始研究随机信号雷达，在相关法、反相关法和频率法实现探测方面取得了一定的研究成果(刘国岁、顾红和苏为民，“随机信号雷达”，国防工业出版社，200。高分辨率的随机噪声探测雷达和成像雷达同样碰到了ADC数据率高，数据存储量、处理量和传输量大的问题，有必要研究随机噪声雷达数据率减低的方法。压缩感知(Compressedknsing，简称CS)是二i^一世纪发展起来的一种稀疏信号处理技术，它从观测信号中获取尽可能少的观测数据，同时尽可能地保留信号中所含的信息，再通过优化算法对原始信号进行有效的逼近和恢复。2006年，Donoho系统阐述了11(Ip表示P阶范数，数学计算的表达式为II·IIp)优化和信号稀疏性之间的关联，同时他还首次使用了压缩感知的概念，针对待估计信号在某个空间可稀疏表征的特性、采用特定的降维压缩采样、利用优化方法实现信号重建，将信号的采样、恢复及信息提取直接建立在信号稀疏特性表征的基础上。同年，Canc^s和Tao指出了I1与Itl优化问题等价的充分条件(RestrictedIsometryftOperty)，阐明了I1优化以及稀疏信号重建之间的关系，为设计高效的重建算法提供了理论指导。压缩感知正是基于下面的一系列文章的[1]D.LDonoho,“Compressedsensing，，，IEEETrans.Inform.Theory,vol.52，no.4，pp.1289-1306，Apr.2006.[2]E.CandesandT.Tao,"Decodingbylinearprogramming,，，IEEETrans.Inform.Theory,vol.51，no.12，pp.4203-4215，Dec.2005.[3]E.Candes,J.Romberg，andT.Tao,“Robustuncertaintyprinciples:Exactsignalreconstructionfromhighlyincompletefrequencyinformation,，，IEEETrans.Inform.Theory,vol.52，no.2，pp.489-509，Feb.2006.[4]E.Candes，J.Romberg，andT.Tao,"Stablesignalrecoveryfromincompleteandinaccuratemeasurements，，，Comm.PureApp1.Math.，vol.59,no.8,pp.1207-1223，2006.[5]E.CandesandT.Tao,“Near-optimalsignalrecoveryfromrandomprojectionsanduniversalencodingstrategies？“IEEETrans.Inform.Theory,vol.52，no.12，pp.5406-5245，Dec.2006.[6]E.CandesandJ.Romberg，"Sparsityandincoherenceincompressivesampling，，，InverseProb.，vol.23，no.3，pp.969-986，June2007.压缩感知是基于这样的假设对于信号、图像和其他的数据，往往存在一个能够稀疏表示它们的基。也就是说对于长度为N的信号χ(η)，存在一个基Ψ，只要Κ(Κ<Ν)个非零的基向量就能表示它X=m=^eiIPi=jz9^,(1)i=lI=I1其中，Φ=[約φ2·■·W7J为χ(η)的稀疏基；θ是一个NXl的列向量，它只有k个元素非零。这个稀疏基可以是小波基、离散余弦基和傅里叶基等。对于具有这样性质的信号x(n)，根据压缩感知理论，可以只需获取M=0(Klog(Ν/Κ))观测量就能很好的恢复出信号。也就是说，我们只需x(n)的M个线性映射_νΟ)=Φχ=〈χ，《〉，πι=1，…，M就能很好的恢复信号。Φ是一个MXN的观测矩阵。由于Μ<Ν，信号χ(η)的恢复是一个欠定的病态问题，根据压缩感知理论，只要观测矩阵Φ和稀疏基矩阵Ψ的乘积满足RIP(R.G.Baraniuk,Μ.Davenport,R.A.Devore,andΜ.B.Wakin.“Asimplifiedproofoftherestrictedisometrypropertyforrandommatrices，”ConstructiveApproximation,2007)，I1优化就能很好的给出值最大的K个稀疏基系数θ—argminθ^st.y=ΦΦ0因此，恢复的结果可以表示为x=W(3)I1优化是数值求解问题，其解法主要分两大类。如凸优化和线性规划等，该类方法及其软件实现都已较成熟并有着广泛应用。另一类是基于贪心算法(GreedyAlgorithm)的方法，它起始于“匹配追踪(MatchingRirsuit)”方法。(2)-类是基于迭代寻优技术的方法，方法在效率方面普遍高于基于迭代寻优技术的方法。目前，优化问题的数值求解方法及其软件实现已较成熟，而且计算效率也在逐步提高。
发明内容本发明的目的是提出一种适用于稀疏微波成像的随机噪声雷达的信号处理方法，对于宽带和高分辨率的随机噪声雷达系统，能够有效的降低ADC的速率，缓解ADC压力；能有效减少信号的存储、处理和传输的数据量。为达到上述目的，本发明的技术解决方案是一种适用于稀疏微波成像的随机噪声雷达的信号处理方法，其包括步骤Si.对具有稀疏特性的目标场景，雷达系统的发射信号是带限的高斯随机白噪声；S2.接收时用低速均勻采样的方法获得比奈奎斯特采样定理要求所需的观测量更少的观测数据；S3.建立信号的回波模型；S4.结合发射信号形式和数据获取方式建立观测矩阵；S5.通过稀疏信号处理理论的压缩感知的信号处理方法，优化求解获得场景目标的后向散射系数；S6.经后续处理，得到高分辨率的目标检测和雷达图像。所述的信号处理方法，其所述Sl步中的带限高斯随机白噪声，系统的热噪声就是其信号源，具有信号实现简单、低截获概率性能和强电子反对抗能力；信号的随机性提高了观测矩阵的RIP，增加了算法准确恢复的概率。所述的信号处理方法，其所述S2步中的观测数据，是由被观测场景的稀疏性决定的；观测数据量决定了随机噪声雷达系统中ADC的信号采样率小于奈奎斯特采样率，大大减小了ADC的压力，以及后续的信号存贮、处理和传输等的压力。所述的信号处理方法，其所述S3步中，回波形式为权利要求1.一种适用于稀疏微波成像的随机噪声雷达的信号处理方法，其特征在于，包括步骤51.对具有稀疏特性的目标场景，雷达系统的发射信号是带限的高斯随机白噪声；52.接收时用低速均勻采样的方法获得比奈奎斯特采样定理要求所需的观测量更少的观测数据；53.建立信号的回波模型；54.结合发射信号形式和数据获取方式建立观测矩阵；55.通过稀疏信号处理理论的压缩感知的信号处理方法，优化求解获得场景目标的后向散射系数；56.经后续处理，得到高分辨率的目标检测和雷达图像。2.根据权利要求1所述的信号处理方法，其特征在于，所述Sl步中的带限高斯随机白噪声，系统的热噪声就是其信号源，具有信号实现简单、低截获概率性能和强电子反对抗能力；信号的随机性提高了观测矩阵的RIP，增加了算法准确恢复的概率。3.根据权利要求1所述的信号处理方法，其特征在于，所述S2步中的观测数据，是由被观测场景的稀疏性决定的；观测数据量决定了随机噪声雷达系统中ADC的信号采样率小于奈奎斯特采样率，大大减小了ADC的压力，以及后续的信号存贮、处理和传输等的压力。4.根据权利要求1所述的信号处理方法，其特征在于，所述S3步中，回波形式为Γ/\272.Rs(t)=Φ一~)ρexpj(W0+6w){t一—)+φ■cIlI式中，t为时间，C为光速，R为目标到天线的距离，a(t)是发射信号的幅度，W0是发射信号的中心频率，Sw是均勻分布的随机频率分量，ρ和φ分别为目标的幅度和相位。5.根据权利要求1所述的信号处理方法，其特征在于，所述S4步中的观测矩阵，是根据发射信号的形式和系统的观测几何构成，该矩阵是带状矩阵，只有矩阵的对角线附近的元素不为零，而且矩阵规模远小于根据满足奈奎斯特采样定理构建的观测矩阵。6.根据权利要求1所述的信号处理方法，其特征在于，所述S5步中的压缩感知的信号处理方法，是基于迭代寻优技术中的凸优化或线性规划方法；或是基于贪心算法的方法。7.根据权利要求1所述的信号处理方法，其特征在于，将观测数据和基于压缩感知的随机噪声雷达的观测矩阵分块分别处理，分块方法减小了信号处理器的处理压力，并提高了处理速度。8.根据权利要求7所述的信号处理方法，其特征在于，所述分块分别处理，包括步骤51.将观测数据根据计算效率的需求分为Q块，记第q块数据为Yes(q)，它对应的观测场景记为X(q)；52.根据和第q块观测数据向量列位置选取观测矩阵中位置相同的部分，对应的分块观测矩阵记为Ocs(q)；53.根据优化算法，恢复出第q块场景；S4将Q块场景的恢复结果组合起来，得到完整区域的观测结果。9.根据权利要求8所述的信号处理方法，其特征在于，所述步骤Sl中分为Q块，其Q块可大小不等。10.根据权利要求1、5或7所述的信号处理方法，其特征在于，所述建立观测矩阵，在观测矩阵实际的构建中，在发射随机信号的同时，通过线性延时线提取出采样率和回波数据获取采样率一致的发射信号,作为观测矩阵的列向量。全文摘要一种适用于稀疏微波成像的随机噪声雷达的信号处理方法，涉及微波成像技术，对具有稀疏特性的目标场景，系统的发射信号是带限的高斯随机白噪声，接收时用低速均匀采样的方法获得比奈奎斯特采样定理要求所需的观测量更少的观测数据。在结合发射信号形式和数据获取方式建立观测矩阵之后，通过稀疏信号处理理论的压缩感知优化求解获得场景目标的后向散射系数，得高分辨率的目标检测和成像。为了提高计算效率采用适用于稀疏微波成像的随机噪声雷达的分块信号处理方法，分块处理时，结合数据采用的分块形式建立相应的分块观测矩阵。与传统的雷达系统相比，本发明在获得相同的分辨率时只需很少的观测数据；在采用相同的观测数据量时可获得更高的分辨率。文档编号G01S13/89GK102207547SQ201010139169公开日2011年10月5日申请日期2010年3月31日优先权日2010年3月31日发明者吴一戎,张冰尘,林月冠,洪文申请人:中国科学院电子学研究所