基于CS理论的LFM信号DOA估计

第26卷第12期计算机应用研究VoL26No.122009年12月ApplicationR七searchofCompute玛Dec.2009

基于CS理论的LFM信号DOA估计水

杨璋,程旺宗

(西南交通大学信息科学与技术学院,成都610031)

摘要:针对常用IjF'M(1inear“qIIencymodulated)信号DOA(directionofarrival)估计算法采样数据量大,在低信噪比情况下估计效果不理想的问题,提出一种基于cs(comp瑚∞d靶nsing或c响pre鲳ives蛐pling)理论的LFM信号DOA估计新方法。实验仿真验证了CS理论的基本原理及运用其进行删信号DOA估计的可行性。

关键词:压缩感知;信号模型;阵列模型;原子库;波迭方向估计

中图分类号:TN91l文献标志码:A文章编号:1001-3695(2009)12—4642-03

doi:10.3969/j.issn.100l?3695.2009.12.067

DOAestimationofLFMsignalsbasedoncompressedsensing

YANGZII肌g,CHENCWang—zong

(S妇f矿删anMl面nS豳W母‰^nn锄,岛“如删砌0fD愕№泌j盼,c^m础6l003l,傀讹)Abstmct:ForthepmbleminDOAe8timationofLFMsi印alsthattheordin撕ly鲥thmeticoffenneedala玎91enuTIlberof

s锄-

ples帅dcan’tgainideale雎ctinconditj叩oflowSNR.ne

paperpmposedanewmethod

ofDOAesti眦tionofIinearh-

quencymodulation(LFM)signalba∞d∞CStheory.Experinlental8imul砒ionpmVidesaproof“t王le

b硇icide船ofCS山eory叽dthefe鹊ibiIityinLFMsignalDOAestimation.

Keywords:comp陀ssed鸵nsing;si印al删)del;阳可rnodel;m蛐library;DOA鹤timation

0引言

随着雷达技术和通信技术的发展。线性调频(LFM)信号得到了广泛的应用,利用阵列信号确定空间信号的波达方向(DOA)有着及其重要的意义。对于非平稳宽带信号的DOA估计,常用方法是基于空间时频分布(sTFD)的子空间法(如Wi印er-viuedistribution,wVD),wVD子空间算法就是利用信号的wVD良好的时频聚集性,先得到信号的空间时频域的数据协方差矩阵,然后对该矩阵进行特征值分解张成子空间(信号和噪声),这样就能利用子空间法得到DOA估计。但是L兀II信号的WVD时频聚集性在低信噪比的情况下下降明显,特别是对于多信源,wVD的交叉项严重影响算法效果,使得在低信噪比时基本不能进行准确的估计。

为了改进在低信号比和多信源情况下的估计效果,本文基于cs理论¨’21提出一种新的LFM信号的DOA估计方法,Do-noh0【3’提出:只要信号是稀疏的或可压缩的,就能使用很少的随机采样点高概率的还原信号。

1LFM信号模型和阵列模型

1.1信号模型

考虑典型的LFM信号:

巾hx—h(厶H等t2)】㈩.其中以为中心频率(起始频率);《为调频斜率,幅度为l。

LFM信号的瞬时频率表达式为

以t)=。岛+丘f(2)瞬时频率随时间呈线性关系。假设以采样频率,对信号进行采样,采样点为Ⅳ=1,2,3,…,Ⅳ,则信号表达式可以转换为

咖,=唧H争n+持)】㈤简化问题为

咖)=e印陋(厶n+茜)】(4)由奈奎斯特采样定律可知Z>2(厶+K),又由宽带信号的定义H],带宽和中心频率之比值介于10%~25%的称为宽带信号(wB),故厶∈(o,o.5],K∈[一O.1,0.1]。

均匀线阵模型如图l所示。

l23工

图l均匀线阵模型

1.2阵列模型

阵列模型为阵元数£的均匀线性阵列,阵元间距为d,假设有叮个LFM信号从不同方向入射到阵列上,则第m个阵元的输出为m=1,2,3,…,三,i=1,2,3,…,g。

‰(I)-耋Jf(I—r。i)+n。(I)(5)

_

r哪=÷(m—1)d8in(巩)(6)

收稿日期:2∞9.03.15:修回日期:2009.05.04.基金项目:国家自然科学基全资助项目(60772084);四川省重点科技计划项目(04GG021.O加一5,2006xl,一038);四J11省应用基础研究项目(04jY029?059-2t20(剐13?114)

作者简介:扬璋(1982一),男.四川达州人。硕士研究生,主要研究方向为信号与信息处理(声ⅡIid9999@163.咖);程旺宗(1984?),男,湖北黄冈人,项士研究生,主要研究方向为信号与信息处理.

万方数据

第12期杨璋,等:基于cs理论的LFM信号DOA估计?4643?

“牡exp【卢盯(胁等r2)】

以小exrh(厶n+精)】将式(7)代入式(5)得(7)

石],p越小矩阵性能越好。同时观测值肘需要满足的条件为

肘≥水l092(MK)(18)(8)

‰(I)=j}.e一伽恤“≯)7眦+M麻‘屯(I)+nm(£)(9)可以得到LFM信号的阵列模型:

z(f)=墨口。(口,f)sl(t)+Ⅳ(f)=A(伊,1)s(f)+Ⅳ(f)(Io)^(日.f)=【口l(口,t),。2(口。1),…,口。(日,‘)】(11)

dl(口,f)=【1,…,e。正1ub+‘一)7脚£+M,1ifi】T(12)考虑到延时fj:I.很小可以忽略掉该项,所以得到最终表示为q(峨,f)=【l,e一皿”蛐+卸)71.,….e一西妇+印)7毗】T(13)时间离散信号的表达式为

毛(巩,厅)=【l,…,e一正。妇+‘叫)由叫扣】T5i(n)(14)其中以和t为已知的参数,需要估计的是每个ⅡM信号毛的到达角吼‘引。

2压缩感知理论

2.1压缩感知理论的基本概念

一维离散时间信号s,可以看做是Ⅳ×l维的列向量虹乃],n=1,2,…,Ⅳo空间信号可以用Ⅳ×l维的基向量{嘭}兰。的线性组合表示。信号s可以表示为

Ⅳc是一个很小的常数[3’引。文献[9]指出随机高斯矩阵与大多数固定正交基构成的矩阵不相关。

最后是从观测值并中恢复原信号,即解式(16)中的0,由于观测值肘远小于信号长度,v,对求解欠定方程组式(16)是无解的。文献[3,6,lO]指出:针对可压缩信号s,求解欠定方程组的问题转换为最小岛范数问题:

min0怕lIos.Lz=蕊=巾怕(19)等价于:

minlIeos.L

z=籼=中怕(20)但是求解zo范数问题是一个NP问题¨o],Chen等人【111指出:求解更简单的z.优化问题会产生同等的解,上式变为

minlIe忆s.Lz=咖=甜旧’(21)这就使问题变成了一个凸优化问题,在有噪声的条件下:聋;=(①,,J)+‰;f=l,2,3,…,肘,吼E∥为噪声,最终求解式为

min¨z—o‘,,e0i+A8e4l(22)需要求的参数为9,即原信号s在基甲上的稀疏表达式可以化简为线性规划,解该类问题可以使用:局部最优的贪婪追踪算法,包括匹配追踪(MP)㈦、正交匹配追踪(OMP)【13J、树型匹配追踪(TMP)㈨;全局最优的基追踪(BP)[3?83等。

3基于CS理论的LFM信号DOA估计

?=。邑研嚷或,=怕(15’3.1建立过完备原子库

其中:9为系数组成的Ⅳxl维的列向量。假如把9的分量按从大到小顺序排列,其按一定量级呈现幂次(州er.1洲)速度衰减,明显大的系数个数为K个,在f0范数意义下可表示为K=0a0。则可以认为信号为K稀疏№】,并且可以用K项稀疏表示很好地逼近原信号。压缩感知理论的主要思想是:把原信号5投影到另一个与变换基甲不相关的朋×.7、r维观测矩阵①上,得到』If×1维的随机采样信号茹肘<<Ⅳ,实现了对原信号s的压缩采样。表达式为

茗=0l=o'f,e(16)并且可以凭借随机采样信号毒,通过将非凸问题转换为凸优化问题求解来精确重建原信号s。

2.2压缩感知理论的基本要求

首先原信号s必须是稀疏的或则可压缩的,并为其找到一个最佳的稀疏域甲,Mallat¨1提出了超完备库下的信号稀疏表示方法,说明了使用超完备冗余字典作为信号稀疏域的必要性。cS理论同时指出信号稀疏度越强(K值越小)重建原信号误差越小旧1。

其次是找到合理的观察矩阵o,根据cs理论,要实现信号s到∞上的投影,即原信号s的压缩采样过程(得到肘个观察值肘<<Ⅳ)同时保证能高概率地从肘次观测中重构出原信号s或信号s在基甲下的稀疏系数向量e。这就需要0甲满足有限等距性质(re8trictedi80metrypmperty,RIP)。文献[8]提出若m与甲满足不连续特性(不相干特性),则其很大概率上满足mP性质,当m与甲是∥空间的一组正交基时,定义。和甲之间的相关性因子:

p(o,甲)2√n×l?曝。I(妒‘,屯)I(17)其中:9。和甲;分别为矩阵中和甲的列向量。卢(0,甲)“l,

根据每个LFM信号的形式特点建立过完备原子库G,si(凡)为1×,v的离散LFM信号厶和峨为已知的参数,研为需要估计的DOA参数,取值是口,∈【一90。~900】(,=1,2,3,…,.,,.,为角度参数搜索个数),原子参数由参数啡单独决定,p,按照需要的搜索精度均匀地取值,构造出原子库G为.,×(Ⅳ+£)的矩阵:

gj2

e,【p(_-.皿霄(.兀田+后i托/^,)dsin吩/c)

唧(一.,21r(,0+t‘,I/^7)(£一1)d8in毋/c)

G=[g,,g:,…,gJ]T.LFM信号s可以表示为

5=∞(24)其中:e为信号s在原子库上投影的稀疏的系数表达式。3.2建立CS观察矩阵

一般情况下高斯随机矩阵可以作为通用的观察矩阵0,本文中使用近似伽分解¨纠的方法来改善高斯随机矩阵性能,同时对分解后的开矩阵对角线的值进行归一化处理,最后得到的矩阵作为观测矩阵①,可以得到更好的估计效果。

3.3由式(16)得到随机投影观测值x

本次投影完成了对信号,的压缩采样,得到肘个随机采样点(M<<Ⅳ),可以利用这肼个随机采样点来恢复原稀疏表达式曰,从而完成原信号DOA估计。

3.4LFM信号的DOA估计

使用式(22)作为解凸优化问题的目标函数,求解信号的稀疏表达式D,找到e中的最大值的坐标,在原子库G中查找此坐标位置的口值,即可得到DOA参数∥16’1引(对于多信源估计需查找最大几个值的坐标下对应的口值)。

万方数据

?4644?计算机应用研究第26卷

4计算机仿真

该仿真在MATLAB7.4a平台下进行,为了更好地描述方法性能.定义均方根误差RMSE。

R眦=以而(25)

其中:凡为估计次数;置表示第i次的估计值;名。。表示真值。

实验l单信源情况下基于cs理论的方法和WvD算法估计效果比较。

单个LFM信号从远场入射到14阵元均匀线阵,方位角200,阵元间距为信号最高频率对应波长的一半(A/2)。LFM信号的起始频率和调谐频率已知,分别为厶=0.35,||}=O.05,离散信号的长度为Ⅳ=256,观测矩阵中使用上面所述方法改进后的高斯随机矩阵,恢复算法为BP算法。压缩采样点的数量为肘=40,肘<<Ⅳ,在信号中加入均值为0方差为l的高斯白噪声,信噪比范围为一10—20dB,每隔5dB分别采用WvD算法和cs算法进行100次Monte.Cado实验,图2给出DOA估计的成功率曲线(wVD算法中认为估计出来的角度误差±0.15。为成功估计,cs方法使用原子库查找的方法,认为准确查找到相应信号模型为成功估计)。图3给出DOA估计的均方根误差(RMSE)曲线(由于观测矩阵采用的是随机矩阵识别概率上会存在小的波动)。

SN刚dBSN刚dB

图2单信源基于、vVD和CS方法图3单信源基于、nD和CS方法的识别概率与信噪比关系的RMSE与信噪比关系

由图2所示随着信噪比的提高基于cs理论的估计方法成功率迅速提高,成功估计概率明显高于WVD算法,由图3所示均方根误差(RMsE)的收敛速度也快于wVD算法。当信噪比大于O时估计成功率接近l。在低信噪比的情况下使用cs理论的估计方法与使用wVD算法相比有明显的优势。

实验2多信源情况下基于cs理论的方法和wVD算法估计效果比较。

两个LFM信号从远场入射到14阵元均匀线阵,信号l方位角20。,信号2方位角为lOo,阵元间距为信号最高频率对应波长的一半(A/2),信号1的起始频率和调谐频率分别为一曲=O.25,|i}.=o.066,信号的长度为川=256;信号2的起始频率和调谐频率分别为厶=0.25,如=O.066,信号的长度为M=256。选择改进后的高斯随机观测矩阵,恢复算法为BP算法,压缩采样点的数量为MI=坞=40,|jIf。<<Ⅳ1.肘2<<他,在两个信号中分别加入均值为O方差为1的高斯白噪声,信噪比范围同时从一10dB到20dB,每隔5dB分别采用wVD算法和cs算法进行100次Monte.cdo实验,图4和5分别给出信号l(20。)和信号2(100)DOA估计的成功率曲线(WVD算法中认为估计出来的角度误差±0.250为成功估计,cs方法使用原子库查找的方法,认为准确查找到相应信号模型为成功估计)。

从图4和5可以看出,随着信噪比的提高cs方法的估计成功率迅速提高,成功识别的概率明显高于wVD算法。由于交叉项的影响在低信噪比情况下wVD算法效果很差,而cS算法与WvD算法相比仍然有明显的优势(由于观测矩阵采用的是随机矩阵识别概率上会存在小的波动)。

鲁0.

|薹|;

SN刚dBSNR/dB图4信号l(20。)到达角的识别图5信号2(10。)到达角的识别概率与信噪比关系概率与信噪比关系以上实验说明了基于压缩感知(CS)理论的估计方法的有效性,在少量采样点和低信噪比的情况可以实现宽带LFM信号的DOA估计,并且获得较好的估计效果。

5结束语

本文阐述了CS理论的基本原理,并且在CS理论指导下,对原信号进行随机压缩采样,用采样得到的很少的点实现了LFM信号的DOA估计。实验仿真验证了cs理论实现LFM信号DOA估计的可行性,并在多种信噪比和多个信源的情况下与常用的WvD算法估计性能进行比较,体现了基于CS理论的DOA估计方法的可行性和优越性。

参考文献:

[1]cAND七sE.compressives删plillg[c]//P眦ofInt锄atioIlalc明?g陀8sofMalhem撕cs.Madrid:【s.n.],2006:1433-1452.

[2]DoNOH0D,TsAIGY.ExtenBion8ofcompr;essedsensing[J].Sig-怕IP啪essing.2006.86(3):533—548.

【3]DONOH0D.comp托s∞d∞聃ing[J].IEEET限nsonInfonTlationTheory,2006,52(4):1289—1306.

[4]王渡.LFM信号检测与参数估计技术研究[D].成都:电子科技大学。2005:l-11.

[5]吕雪.基于肿分解的宽带LFM信号参数估计[D】.西南交通大学,2003.

[6]cAND£sE,TAOT.Ne盯opIimalBi目rIalDeco代ryfroIn瑚domP睁jection8:universalencoding8t随te画es[J].IEEETransonInfor『T1a.tbnTh∞ry,2006,52(12):5406—5425.

[7]MAuATs.Matchingpu鹅ui协’订thtime—f呻encydicti佣a—es[J].1EEET陷n8onSignaIProce8s,1993,41(12):3397.3415.[8]o悄DEsE,RoMBERGJ,TAOT.R0b璐t佃。酣ainty一他ipI墨:∞c£sigrIalrec叽shlldionf而mhi圳y

i咖lpletef两uerIcyi—湘ti∞[J].IE匪T旧憾onl限'丌na【i∞Theory,2006,52(2):4朗.509.

[9】cANDEsE。wAKlNMB.AnintTDduc“ontocoInp佗鹪ive明巾pling[J].1EEESig怕IP代'∞ssingMagazine,2008,25(2):2l?30.[10]BARAN儿JKR.AIectureoncomp陀ssivesen8ing[J】.IEEESjgnaIPnocessjngMagazjne.2007,24(4):181—121.

[11]cHENss,DONOH0DL,sAuNDERsMA.At伽jcde咖poBi6鲫bybasi8Pmuit[J].SlAMR酬洲,200l,43(1):129-159.

[12]TROPPJA.GILBERTAc.si伊laJrecoveryfmmpanialinfbmlationbyⅢhogolIal眦tchiIIgPu舢it[J].IEEET限nsonInfo肿atjonTheory.2007,53(12):4655?4666.

[13]ELADM.oPIimizedpmjectio岫forcompreBsed蚰n8illg[J].IEEETransonSignaIPnocessing,2007,52(12):5695—5702.

[14]LAc,D0MN.si印且Jr即onstmctionusi“gspar眈旬ree他pre鲫tation[C]//Procofw删elte8ⅪmsP匝O—cs帅dPllotoIIi∞.s明DiegD:[8.n.],2005.

[15]傅迎华.可压缩传感重构算法与近似QR分解[J].计算机应用.2008。28(9):2300?2303.

[16】f1JcHsJJ.bnwp删inginspectml∞tjmti蛐:印pljcationt0帅yprocessing[c]//ProcoflcAssP.1996:316l-3164.

[17]cEVHERV,GuRBuzAc,MccLELLANH,eI耐.c伽p嘲ivewi托lesB删y8fhb谢nge出毗ti仰[c]//Pr∞0fICAsSP.2008:2497—2500.

[18]GuRBuzAC,MCcLELLANJH.Ao咖pre聃ivebe&mfomingmethod【c]//P啊ofJcAssP.2008:2617.2620.

万方数据

相关文档
最新文档