法律状态公告日
法律状态信息
法律状态
2023-03-07
未缴年费专利权终止 IPC(主分类):G01S 3/00 专利号:ZL2014101080821 申请日:20140321 授权公告日:20160706
专利权的终止
2016-07-06
授权
授权
2014-07-30
实质审查的生效 IPC(主分类):G01S3/00 申请日:20140321
实质审查的生效
2014-07-02
公开
公开
技术领域
本发明涉及一种通过量子万有引力搜索机制和加权信号协方差极大似然方程的动态更新 共同来实现的考虑在冲击噪声环境下对动态目标的波达方向进行鲁棒跟踪的冲击噪声环境下 的量子万有引力搜索动态DOA估计方法。
背景技术
动态DOA估计技术作为无线通信系统、雷达和导弹制导的关键技术,是决定能否快速 跟踪信源信号,实现精确定位的关键。子空间更新算法是动态DOA估计领域中研究的重点。 其中子空间更新算法按其原理可分成三类:第一类是静态的DOA估计方法直接推广,如幂 迭代算法等基于特征分解批处理的子空间更新算法。第二大类是是将空间谱估计中的求解子 空间更新问题转化为最优化问题,然后求出最优解。如带约束的最优化问题梯度法和无约束 的最优化问题的投影逼近子空间跟踪。第三类是基于秩-1子空间更新算法,分别将协方差矩 阵的信号子空间及噪声子空间进行平均,从而得到一个秩等于2的矩阵,数据的更新就是对 这个矩阵的更新,典型算法如ROSA算法等。基于子空间的动态DOA估计虽具有计算量小, 处理速度快的特点,但步骤繁琐,在相干信源跟踪场合还需引入空间平滑等解相干算法,可 能会带来有效孔径的损失;而且在低信噪比情况下的跟踪性能并不理想。最大似然算法作为 一种经典的DOA估计方法,不仅具有优越的估计性能和解相干的能力。而且该算法直接能 处理数据协方差矩阵,在动态DOA估计中能避免子空间类跟踪方法需要不断重复的协方差 矩阵分解问题,可节省计算量。
经对现有技术文献的检索发现,刁鸣等在《系统工程与电子技术》(2009,Vol.29,No.12, pp.2046–2049)上发表的“一种新的基于粒子群算法的DOA跟踪方法”中使用了粒子群算法 和极大似然算法进行相干信号的动态DOA估计,在高斯噪声背景下效果较好,但在冲击噪 声背景下,方法性能恶化严重导致失效。赵大勇等在《山东大学学报》(2010,Vol.40,No.1, pp.133-138)上发表的“冲击噪声背景下的动态DOA跟踪”提出了粒子群算法和分数低阶矩的 极大似然方法求解弱冲击噪声环境下的动态DOA估计问题,使用分数低阶矩的动态更新能 在一定程度上实现动态DOA估计问题,但计算时间较长,依旧不能解决冲击噪声环境下的 动态DOA估计的收敛速度和精度问题,在特征指数小于1的强冲击噪声情况下基于粒子群 算法和分数低阶矩的极大似然方法就失效了。
已有文献表明,在恶劣噪声环境下,动态DOA估计问题依旧没有解决。针对冲击噪声 环境下的DOA估计是一个难题,尤其对相干信源的动态跟踪更难,故不能把高斯噪声环境 下的DOA估计算法直接进行移植。在冲击噪声下的鲁棒动态DOA估计首先应能建立高性能 的鲁棒动态跟踪的优化方程,冲击噪声环境下经典智能计算方法很难摆脱收敛速度和收敛性 能矛盾的制约,在现有计算条件下很难在有限的时间内搜索到最优解,需要设计新的智能算 法求解冲击噪声环境特别是强冲击噪声环境下的鲁棒动态DOA估计问题。
运用空间平滑等解相干技术对数据协方差矩阵进行预处理虽能解相干但损失阵列孔径。 而极大似然估计算法作为一种经典的DOA估计方法,在高斯噪声环境下不仅具有优越的估 计性能,还避免了子空间跟踪类方法需要不断重复协方差矩阵分解的问题。但极大似然算法 角度估计的另一缺点是搜索过程是一个复杂耗时的过程,运算量巨大。因此,本专利在冲击 噪声环境提出了一种基于加权信号协方差的极大似然算法的动态DOA估计方法。该方法先 将目标锁定在一个变化的搜索范围之内,再设计量子万有引力搜索机制在该搜索空间搜索加 权信号极大似然方程的所蕴含的最优角度。通过逐渐减小搜索的范围以及运用智能搜索机制, 有效解决了搜索方法的计算量问题。仿真结果表明这种冲击噪声环境下的动态目标DOA估 计方法能够保证所设计方法的实时性,而且具备解相干的能力和较好的跟踪精度,尤其在小 特征指数情况下,有非常明显的优势。
发明内容
本发明的目的在于提供一种有效性和鲁棒性更高的冲击噪声环境下的量子万有引力搜索 动态DOA估计方法。
本发明的目的这样实现的:
步骤一,对空间D个信源信号进行快拍采样和加权处理,由M个全向阵元构成的无模糊 线阵,D个中心频率相同的远场窄带信号从方向{θ1,θ2,…,θD}入射到阵列上,阵列接收第k 次快拍采样数据为x(k)=As(k)+n(k),式中x(k)=[x1(k),x2(k),…,xM(k)]T为阵列的接收数据 矢量;A(θ)=[a(θ1) a(θ2) … a(θD)]是信号导向矢量矩阵,θ=(θ1,θ2,…,θD)是角度矢量,θd(d=1,2,…,D)是第d个信源的来波方向;s(t)=[s1(t),s2(t),…,sD(t)]T为信号矢量;
n(k)=[n1(k),n2(k),…,nM(k)]T是独立同分布的满足SαS分布的冲击噪声向量;第d个导向矢 量为a(θd)=[exp(-jω0τ1d),exp(-jω0τ2d),…,exp(-jω0τMd)]Τ,d=1,2,…,D,式中, c为光速,f为入射信号的频率,λ为入射信号的波长,(l=1,2,…,M)表示第d个信号到达第l个阵元时相对于参考阵元的时延,El为第l个阵元的位 置,接收数据的加权信号为h(k)=x(k)/max{|x1(k)|0.7,|x2(k)|0.7,…,|xM(k)|0.7},max{} 为取最大值函数,定义阵元接收数据之间的加权信号协方差为C(k),第一次采样数据的加权 信号协方差为C(1)=h(1)hH(1),根据经典信息论准则估计信源个数D,t代表量子万有引力 搜索机制的迭代次数,初始时设t=0;
步骤二,量子万有引力搜索机制的所有物体初始状态确定,由N个物体组成的系统,每 个物体在D维搜索空间中运动,第i个物体的量子位置定义如下:
步骤三,第i个物体位置qi(t)=[qi1(t),qi2(t),…,qiD(t)]的加权信号协方差极大似然方程的 目标函数值为
步骤四,计算在t次迭代的每个物体质量:Mai=Mpi=Mii=Mi,i=1,2,…,N,
步骤五,量子位置yi(t+1)影射到定义区间qi(t+1),计算目标函数值,
fit[qi(t+1)]>fitq[pi(t)],则令局部最优量子位置否则,
pi(t+1)=pi(t),把最优的局部最优量;子位置设置成全局最优量子位置为
b(t+1)=[b1(t+1),b2(t+1),…,bD(t+1)],确定当前最好目标函数值best(t+1)和最差目标函数 值worst(t+1),
步骤六,判断是否达到最大迭代次数,若是,记录最优位置,执行步骤七;否则,t=t+1, 返回步骤四;
步骤七,快拍采样新数据x(k+1)=[x1(k+1),x2(k+1),…,xM(k+1)]T,对采样数据进行加权处 理:h(k+1)=x(k+1)/max{|x1(k+1)|0.7,|x2(k+1)|0.7,…,|xM(k+1)|0.7},计算当前采样 的加权信号协方差:
步骤八,如果达到最大跟踪次数,执行步骤九;否则,设k=k+1,t=0,返回步骤二 继续估计动态目标下一个时刻的方向;
步骤九,得到所有快拍采样下的全局最优位置就是检测的动态目标方向值,输出动态跟 踪结果。
本发明的有益效果在于:
(1)本发明解决了冲击噪声环境下的动态目标方向跟踪问题,使用所设计的量子万有引 力搜索机制作为搜索策略,所设计的方法具有快速和鲁棒性高的优点。
(2)相对于现有的弱冲击噪声环境下的动态目标跟踪方法,本发明可以同时利用量子万 有引力搜索机制和加权信号协方差极大似然优化方程的优势解决高斯噪声和冲击噪声环境下 的动态目标跟踪问题,且适用于特征指数小于1的恶劣冲击噪声环境,说明本方法的适用性 更广。
(3)实验结果表明,本发明所提出的量子万有引力DOA估计方法在对动态目标跟踪时能 够得到最优解,但时间开销远小于已有极大似然测向方法,估计精度却优于现有的极大似然 方法,说明了所提鲁棒动态跟踪方法的有效性。
附图说明
图1为量子万有引力动态DOA估计的跟踪方法示意图。
图2为量子万有引力搜索机制的结构示意图。
图3为特征指数α=1.4时两个独立信源的粒子群算法和分数低阶矩极大似然算法的动态 DOA估计跟踪情况。
图4为特征指数α=1.4时对于两个独立信源所提的量子万有引力动态DOA估计方法的 角度跟踪情况。
图5为特征指数α=1.4时两个相干信源粒子群算法和分数低阶矩极大似然算法的动态 DOA估计跟踪情况。
图6为特征指数α=1.4时对于两相干源所提的量子万有引力动态DOA估计方法的角度 跟踪情况。
图7为特征指数α=0.9时两个独立信源的基于粒子群算法和分数低阶矩极大似然算法的 动态跟踪情况。
图8为特征指数α=0.9时两个独立信源的所提的量子万有引力动态DOA估计方法的角 度跟踪情况。
图9为量子万有引力动态DOA估计方法的总体结构图。
具体实施方式
下面结合附图对本发明做进一步描述。
本发明针对现有冲击噪声环境下动态跟踪方法的不足,提出了一种可在复杂噪声环境和 恶劣测向背景达到鲁棒跟踪的动态目标跟踪方法。该方法先设计一种加权信号协方差的更新 方法,并提出解决连续优化问题的量子万有引力搜索算法去求解动态目标跟踪。在工程应用 中,当冲击噪声的特征指数为2时满足高斯噪声分布函数形式,故本发明所提出的方法也能 够解决高斯噪声环境的动态跟踪问题。
本发明是通过如下技术方案来实现的,主要包括以下步骤:
步骤一,对空间D个信源信号进行快拍采样和加权处理。由M个全向阵元构成的无模糊 线阵,D个中心频率相同的远场窄带信号从方向{θ1,θ2,…,θD}入射到阵列上,阵列接收第k 次快拍采样数据为x(k)=As(k)+n(k),式中x(k)=[x1(k),x2(k),…,xM(k)]T为阵列的接收数据 矢量;A(θ)=[a(θ1) a(θ2) … a(θD)]是信号导向矢量矩阵,θ=(θ1,θ2,…,θD)是角度矢量,θd(d=1,2,…,D)是第d个信源的来波方向;s(t)=[s1(t),s2(t),…,sD(t)]T为信号矢量; n(k)=[n1(k),n2(k),…,nM(k)]T是独立同分布的满足SαS分布的冲击噪声向量;第d个导向矢 量为a(θd)=[exp(-jω0τ1d),exp(-jω0τ2d),…,exp(-jω0τMd)]Τ,d=1,2,…,D,式中, c为光速,f为入射信号的频率,λ为入射信号的波长,(l=1,2,…,M)表示第d个信号到达第l个阵元时相对于参考阵元的时延,El为第l个阵元的位 置。接收数据的加权信号为h(k)=x(k)/max{|x1(k)|0.7,|x2(k)|0.7,…,|xM(k)|0.7},max{} 为取最大值函数。定义阵元接收数据之间的加权信号协方差为C(k),第一次采样数据的加权 信号协方差为C(1)=h(1)hH(1)。根据经典信息论准则估计信源个数D。t代表量子万有引力 搜索机制的迭代次数,初始时设t=0
步骤二,量子万有引力搜索机制的所有物体初始状态确定。考虑由N个物体组成的 系统,每个物体在D维搜索空间中运动,物体的位置代表优化问题的解。第i个物体的量子位 置定义如下:
qi(t)=[qi1(t),qi2(t),…,qiD(t)],对于所有的i=1,2,…,N,
qid(t)=yid(t)[ud(k)-ld(k)]+ld(k),d=1,2,…,D,第i个物体对应的速度为
vi(t)=[vi1(t),vi2(t),…,viD(t)],-0.2≤vid(t)≤0.2(d=1,2,…,D)表示第i个物体的第d维速度,t 代表迭代次数。为了使初始位置具有一定的分散性和均布性,量子位置在定义的量子域[0, 1]间随机初始化,在[-0.2,0.2]随机初始化速度的每一维。确定初始搜索空间
步骤三,第i个物体位置qi(t)=[qi1(t),qi2(t),…,qiD(t)]的加权信号协方差极大似然方程的 目标函数值为,式中,正交投影矩阵
步骤四,计算在t次迭代的每个物体质量:Mai=Mpi=Mii=Mi,i=1,2,…,N,
步骤五,量子位置yi(t+1)影射到定义区间qi(t+1),计算目标函数值,
步骤六,判断是否达到最大迭代次数,若是,记录最优位置,执行步骤七;否则,t=t+1, 返回步骤四。
步骤七,快拍采样新数据x(k+1)=[x1(k+1),x2(k+1),…,xM(k+1)]T,对采样数据进行加权处 理:h(k+1)=x(k+1)/max{|x1(k+1)|0.7,|x2(k+1)|0.7,…,|xM(k+1)|0.7},计算当前采样 的加权信号协方差:
步骤八,如果达到最大跟踪次数,执行步骤九;否则,设k=k+1,t=0,返回步骤二 继续估计动态目标下一个时刻的方向。
步骤九,得到所有快拍采样下的一系列全局最优位置就是检测的动态目标方向值,输出 动态跟踪结果。
本发明考虑到冲击噪声环境下完成动态跟踪的过程中能够同时考虑收敛精度,收敛速度 和去相干能力,使用量子引力搜索机制去搜索加权信号协方差极大似然方程中所包含的的波 达方向。所设计的动态跟踪系统还可以根据时间要求和性能要求自动确定迭代次数,从而使 所设计的量子万有引力动态跟踪方法满足更高性能要求。
实验中均使用阵元间隔为0.5倍波长的等距均匀线阵,阵元数M=8,搜索区间为[0,90°]。 所提的量子万有引力动态DOA估计方法,由于对接收数据使用了加权信号协方差处理方法, 所以记作QGSA-WSC,方法中参数G0=100,γ=20,ε=2.2204×10-16,加权信号协方差C更 新公式中的μ=0.95,方法中收敛因子取β=0.995,遗传因子δ=0.8,收敛半径取r=3。在 实验仿真过程中使用的用于比较的动态跟踪方法有:基于粒子群算法和分数低阶矩极大似然 算法(PSO-FLOM)的动态DOA估计(赵大勇等在《山东大学学报》(2010,Vol.40,No.1, pp.133-138)上发表的“冲击噪声背景下的动态DOA跟踪”)。为了考察两种动态DOA估计的收 敛速度和性能,量子万有引力动态DOA估计方法和粒子群算法的种群规模都设为30,最大迭 代次数设置方法相同,都取该次快拍第一维搜索边界区域差的整倍数,倍数值为zz=4。
假设冲击噪声环境下的广义信噪比GSNR=10dB,等功率两独立信源在不同冲击噪声环 境的两种方法跟踪结果如图3至图4,等功率两相干信源在不同冲击噪声环境的两种方法跟 踪结果如图5至图6。
从图3和图4可以看出,所提量子万有引力动态跟踪方法在冲击噪声背景下对独立信源 的跟踪精度要好于基于粒子群算法和分数低阶矩极大似然算法的动态DOA估计方法。从图5 和图6可以看出,所提动态跟踪算法在冲击噪声背景下对相干信源的跟踪速度和精度要好于 基于粒子群算法和分数低阶矩极大似然算法的动态DOA估计,有更好的收敛精度。
从图7和图8可以看出,在特征指数小于1的强冲击噪声环境下所提的方法对相干信源 的估计准确性要远小于基于粒子群算法和分数低阶矩极大似然算法,基于粒子群算法和分数 低阶矩极大似然算法的动态DOA跟踪完全失败,因为基于粒子群算法和分数低阶矩极大似 然算法的动态DOA跟踪仅适用于特征指数大于1的情况。
以图9的鲁棒动态跟踪装置为例,其它情况可以依此类推,在动态DOA估计系统存在M 个天线,开始对D个动态目标所发射的信号进行快拍采样,在冲击噪声环境下,使用已有信 源个数估计算法估计出D个信源数,计算初始加权信号协方差,使用量子万有引力搜索机制 估计波达方向的初始值。以波达方向估计值为中心确定下一个采样的搜索范围,对接受到的 当前快拍采样数据和前一快拍采样信号加权协方差用于更新加权信号协方差,使用量子万有 引力搜索机制去优化信号协方差极大似然方程,使用两种量子行为更新物体的速度和量子位 置,迭代求得物体的最优量子位置,根据映射态求得动态DOA估计的波达方向值,每个快 拍采样后的都循环迭代上述DOA估计过程,完成波达方向的动态跟踪过程。
冲击噪声环境下的量子万有引力搜索动态DOA估计方法,其综合特征是:冲击噪声环 境下基于加权信号协方差极大似然优化方程的建立;(2)解决连续优化问题的量子万有引力 搜索方法;(3)在量子万有引力搜索方法的速度和量子位置的更新过程中,所有物体分别使 用两种速度更新行为去迭代更新物体的速度。冲击噪声环境下加权信号协方差极大似然数据 处理方法的建立,根据新的加权信号,设计了加权信号协方差极大似然方程。充分考虑了冲 击噪声和相干信源的动态估计,可解决现有鲁棒动态跟踪方法在时间和估计精度之间的矛盾。 基于量子万有引力搜索算法的动态DOA估计方法,使用量子旋转门完成量子位置的更新, 可快速找到一个全局最优解,在演进的每一代物体进行经验交流,保留最优位置指导物体演 进。充分考虑了冲击噪声、相干信源的动态跟踪,可解决现有动态DOA估计方法在强冲击 噪声环境下也有好的鲁棒性和较高的收敛精度。两种速度更新方式。每个物体根据自己的已 有状态使用两种方式更新物体速度进而产生新的量子位置。基于量子万有引力演化机理的搜 索角度方法,在解决高斯噪声下的动态跟踪问题时同样具有收敛精度高、收敛速度快的优点。
机译: 相关天线噪声环境中相关多个信号的高分辨率DOA估计方法
机译: 无线位置系统中估计TDOA和FDOA的改进方法无线位置系统中估计TDOA和FDOA的改进方法
机译: 基于TDoA的TDoA步枪冲击点估计系统