一种统计mimo雷达多目标检测仿真方法
【专利摘要】本发明公开了一种统计MIMO雷达多目标检测仿真方法，其特征在于，包括如下步骤：步骤1、构建MIMO雷达模型；步骤2、利用功率谱密度求解MIMO雷达的特征值和特征向量；步骤3、基于信号和噪声特征值大小差异，进行多目标迭代检测仿真。本发明针对不完全分集信道情况下，提出MIMO雷达多目标检测性能算法，解决了在真实的信道环境中，由于各种干扰存在以及复杂的电磁环境下，可能出现信道部分相关，提出了符合实际信道环境下的MIMO雷达多目标检测方法。本发明可以模拟真实任意信道环境，无论是完全分集还是不完全分集情况，用于MIMO雷达多目标检测。
【专利说明】
一种统计MI MO雷达多目标检测仿真方法
技术领域
[0001 ]本发明涉及一种统计MBTO雷达多目标检测仿真方法。
【背景技术】
[0002] ]\nMO(Multiple-Input Multiple-Output:多输入多输出技术）是指利用多发射和 多接收天线，将通信链路分解成为许多并行的子信道，通过多路空间信道的复用和分集技 术，从而大大提高系统容量和降低误码率。
[0003] MIMO技术已经成为无线通信领域的关键技术之一，通过近几年的持续发展，MMO 技术将越来越多地应用于各种无线通信系统。中外学者在MBTO无线通信技术的启发下，针 对现代雷达系统高性能的要求，提出了 MMO雷达新体系。MMO雷达定义为多天线协作发射 多重信号，并使用多接收天线联合进行目标检测与参数估计的雷达系统。目标检测要求对 由信号、噪声和其他干扰信号组成的接收机输出进行处理，使得以规定的高概率输出所希 望得到的信号，而以低概率输出噪声和其他干扰信号的随机虚警。
[0004] 在真实的信道环境中，由于各种干扰存在以及复杂的电磁环境下，可能出现信道 部分相关的情况，故常规检测算法中，MHTO雷达利用多个阵列同时发射可区分的信号，利用 单目标雷达截面积的闪烁特性引入参数上有差异的目标观测通道，基于NP准则获得检测和 估计增益，提高雷达对各种空间目标和隐身目标检测、识别和参数估计的性能。但常规检测 算法需要计算多维高斯协方差矩阵的行列式、逆矩阵、特征值和特征向量，计算复杂度过 高，而且不适用多目标的检测。
【发明内容】
[0005] 针对上述问题，本发明提供一种统计MIMO雷达多目标检测仿真方法，首先利用功 率谱密度来求解特征值和特征向量，然后基于特征门限理论进行多目标检测，针对真实信 道中可能存在的情况，提出了符合实际信道环境下的MMO雷达多目标检测方法。
[0006] 为实现上述技术目的，达到上述技术效果，本发明通过以下技术方案实现：
[0007] -种统计MBTO雷达多目标检测仿真方法，其特征在于，包括如下步骤：
[0008] 步骤1、构建MBTO雷达模型；
[0009] 步骤2、利用功率谱密度求解MMO雷达的特征值和特征向量；
[0010] 步骤3、基于信号和噪声特征值大小差异，进行多目标迭代检测仿真。
[0011] 优选，假设MMO雷达模型由M个发射天线、N个接收天线组成，则第j个接收天线接 收的信号yj表不为：
[0012]
[0013]式中，X1表示第i个发射天线所发射的信号，表示第j个接收天线上的总噪声，hu 表示通道增益，j = ^发射功率归一化因子；
[0014] 假设发射信号正交，接收端进行匹配滤波并转换成矩阵Y，则：
[0015]
[0016] 式中，Y=[yj]MNX1，H=[hij]MN，n=[nj]MNxi〇
[0017] 优选，利用功率谱密度求解特征值\和特征向量U1的方法如下：
[0018]
[0019] ⑵
[0020]式中，fi = i/L，i=0，l，'"，L_l，Pyy〇是功率谱密度函数，L是样本个数;exp是以自 然常数e为底的指数函数。
[0021] 优选当目标观测通道完全独立时，当Rr =(1/Μ>ΙΛ?ν ,，则共有MN个独立的目 标观测通道，有MN个相等的特征
[0022] 当目标观测通道不完全独立时又分为完全相关和部分相关两部分：
[0023] 若目标观测通道完全相关，当Kt = G / Λ?)ΙΛ?ν + d、，则共有MN个特征值，其中一个 特征值为為=W+ <:，其余特征值为4 ;
[0024]若目标观测通道部分相关，则假设接收信号表示为y=n+j γ，求得：
[0025]
[0026] 通过构建I AImn-Ry I =〇特征方程依次求得参数λ，对应为特征值A1，…，λΜΝ;
[0027] 式中，Ry是接收?目号协方差矩阵，Imn是丽维单位矩阵，是噪声方差，rI和y是接收 信号的实部和虚部随机矢量，EU是求解数学期望值。
[0028] 优选，基于信号和噪声特征值大小差异估计目标个数，设存在k个目标，假设噪声 特征值多样性为q，则信号的最小特征值AK = AM_q应超过所设定的检测门限，而Ν-k个噪声的 特征值Ak+1，…，λ Μ应低于检测门限，检测流程如下：
[0029] 步骤a、首先从q = l开始确定二元假设如下：
[0030] Ho = k<N-q
[0031] Hi = k = N-q
[0032] 步骤b、利用信号特征值大于噪声方差特征值，构建基于特征值门限的循环表达式 为：
[0033]
[0034]式中，为噪声特征值的上界，Ho表示目标不存在，出表示目标存在;将求解的特 征值A1排序代入循环表达式，其中，若M小于XuN-q，则满足Ho,则令q = q+l返回步骤a继续检 验直至满足Hi或者q=M;若λκ满足Hi，即AK>A uN-q，则停止检验，此时检测到的目标数目k = N_ q〇
[0035]本发明的有益效果是：本发明针对不完全分集信道情况下，提出MMO雷达多目标 检测性能算法，解决了在真实的信道环境中，由于各种干扰存在以及复杂的电磁环境下，可 能出现信道部分相关，提出了符合实际信道环境下的MMO雷达多目标检测方法。本发明可 以模拟真实任意信道环境，无论是完全分集还是不完全分集情况，用于MMO雷达多目标检 测 。
【附图说明】
[0036]图1是本发明一种统计MBTO雷达多目标检测仿真方法的流程图；
[0037]图2是本发明MBTO雷达模型的示意图；
[0038]图3本发明中MBTO雷达单目标检测曲线；
[0039]图4本发明中MBTO雷达双目标检测曲线。
【具体实施方式】
[0040]下面结合附图和具体的实施例对本发明技术方案作进一步的详细描述，以使本领 域的技术人员可以更好的理解本发明并能予以实施，但所举实施例不作为对本发明的限 定。
[0041] -种统计MBTO雷达多目标检测仿真方法，包括如下步骤：
[0042] 步骤1、构建MBTO雷达模型；
[0043] 步骤2、利用功率谱密度求解MMO雷达的特征值和特征向量；
[0044] 步骤3、基于信号和噪声特征值大小差异，进行多目标迭代检测仿真。
[0045] 在实际仿真过程中，在完全分集信道情况下和在部分相关信道情况下，检测器的 信噪比计算都要求计算多维高斯协方差矩阵的行列式、逆矩阵、特征值和特征向量，需要耗 费大量的资源。
[0046] 而本发明中，若向量Y1xl是源于零均值广义平稳高斯随机过程的随机矢量，当样本 个数L-〇〇时，接收信号协方差矩阵Ry的特征值A i和特征向量Ui可以利用功率谱密度来求 解，并且考虑到复杂的信道传输环境，首先利用功率谱函数推导出特征值和特征函数，然后 借助门限定理给出M頂0雷达多目标检测的统一仿真方法。下面结合图1进行详细的介绍。
[0047] 如图2所示，假设MMO雷达模型由M个发射天线、N个接收天线组成，则第j个接收天 线接收的信号w弄示为，
[0048]
[0049] 式中，X1表示第i个发射天线所发射的信号，如表示第j个接收天线上的总噪声，Iuj 表示通道增益，j = 1，…
％发射功率归一化因子；
[0050] 假设发射信号正交，接收端进行匹配滤波并转换成矩阵Y，则：
[0051
[0052」 式干，Y = Lyj」MNX1，H=[hij]MN，n=[nj]MNxi〇
[0053] 利用功率谱密度求解特征值M和特征向量Ul的方法如下：
[0054] Ai = Pyy(fi) (1)
[0055]
⑵
[0056] 式中，fi = i/L，i=0，l，'"，L_l，Pyy〇是功率谱密度函数，L是样本个数;exp是以自 然常数e为底的指数函数。
[0057]得到特征值后信号处理的全过程如下：
[0058] 完全分集下特征值分布情况:当目标观测通道完全独立时，当Ri. = (IZM)Iuv ， 即MIMO雷达完全满足分集条件，则共有MN个独立的目标观测通道，有MN个相等的特征值 為=丄+ σ..:。
[0059] 当目标观测通道不完全独立时又分为完全相关和部分相关两部分：
[0060] 完全相关下特征值分布情况：当Rr =(1/Μ)?,1?ν + 〇·χν，贝拱有MN个特征值，其中一 个特征值为4 =# + <，其余特征值为4 = σ,?。
[0061 ]信道部分相关情况下特征值分布情况：
[0062] 不完全分集信道情况下，根据通道矩阵矢量均为复高斯过程，假设接收信号表示 为y=n+i γ，龙得：
[0063]
[0064] 通过构建I AImn-Ry I =〇特征方程依次求得参数λ，对应为特征值A1，…，λΜΝ;
[0065]式中，Ry是接收?目号协方差矩阵，Imn是MN维单位矩阵，是噪声方差，η和γ是接收 信号的实部和虚部随机矢量，EU是求解数学期望值。
[0066]在部分相关信道情况下，此时各观测角度不再独立，Ry产生信号降维，其空间维数 介于1和MN之间，其特征值可能互不相等，也可能存在部分相同。
[0067] 本文算法检验的是噪声特征值的多样性，设存在k个目标，假设噪声特征值多样性 为q，则信号的最小特征值λκ = λ"ι应超过所设定的检测门限，而Ν-k个噪声的特征值 Ak+1，…，λΜ应低于检测门限，检测流程如下：
[0068] 步骤a、首先从q = l开始确定二元假设如下：
[0069] Ho = k<N-q
[0070] Hi = k = N-q
[0071] 步骤b、利用信号特征值大于噪声方差特征值，构建基于特征值门限的循环表达式 为：
[0072]
[0073]式中，AuN-q为噪声特征值的上界，Ho表示目标不存在，出表示目标存在;将求解的特 征值A1排序代入循环表达式，其中，若M小于XuN-q，则满足Ho,则令q = q+l返回步骤a继续检 验直至满足Hi或者q=M;若λκ满足Hi，即AK>A uN-q，则停止检验，此时检测到的目标数目k = N_ q〇
[0074]任意信道情况下，在分析特征值基础上，根据门限定理，得到统一目标检测表达 式：
[0075] 表达式含义为：从q= I，将求解的特征值11排序代入循环表达式，其中，若A1小于 A1Vq，则满足Ho,则令q = q+l返回步骤a继续检验直至满足H1或者q = M;若λκ满足H1，即λκ> AVq，则停止检验，此时检测到的目标数目k = N_q。
[0076] 图3是本发明中的0°单目标检测曲线，仿真所用的统计MBTO雷达系数为:M=2，N = 4。图4是本发明中0°、15°双目标检测曲线，仿真所用的统计M頂0雷达系数为:M=2，N=4。
[0077] 从图3和图4中可以看到，本文中提出的检测方法在单目标和双目标检测上都满足 检测性能曲线Pd vs SNR变化趋势，即满足单目标乃至多目标下目标检测应用，可以用于多 目标检测。在低信噪比情况下，各信道之间的相关性越强其检测性能越好，MIMO雷达的检测 性能随着观测通道相关性的增强而改善，这一点容易理解，在低信噪比条件下，相干处理增 益要高于空间分集对检测性能的影响，在高信噪比情况下，信道之间的相关性会引起一定 的检测性能损失，M頂0雷达的检测性能优于相控阵雷达，通道之间的相关性越小，MIMO雷达 的检测性能越好，随着相关性的增强，检测性能会随之下降。改变发射、接收阵元数目，上述 结论仍然成立，此时性能检测曲线随着发射、接收阵元数目稍微变化。
[0078]本专利在多目标散射模型的基础上，建立了回波信号的模型并分析了Mnro雷达满 足空间分集工作的条件，针对传统检测算法需要计算多维高斯协方差矩阵的行列式、逆矩 阵、特征值和特征向量，本发明提出了大数据情境下利用功率谱密度函数仿真计算的MIMO 雷达不同信道环境下的特征值和特征向量，在特征值的基础上，利用噪声特征值的多样性 估计目标个数，从而实现MMO雷达多目标检测的统一算法。
[0079] 解决了在真实的信道环境中，由于各种干扰存在以及复杂的电磁环境下，可能出 现信道部分相关，提出了符合实际信道环境下的MMO雷达多目标检测方法。本发明可以模 拟真实任意信道环境，无论是完全分集还是不完全分集情况，用于M頂0雷达多目标检测。
[0080] 以上仅为本发明的优选实施例，并非因此限制本发明的专利范围，凡是利用本发 明说明书及附图内容所作的等效结构或者等效流程变换，或者直接或间接运用在其他相关 的技术领域，均同理包括在本发明的专利保护范围内。
【主权项】
1. 一种统计ΜΙΜΟ雷达多目标检测仿真方法，其特征在于，包括如下步骤： 步骤1、构建ΜΙΜΟ雷达模型； 步骤2、利用功率谱密度求解ΜΙΜΟ雷达的特征值和特征向量； 步骤3、基于信号和噪声特征值大小差异，进行多目标迭代检测仿真。2. 根据权利要求1所述的一种统计ΜΙΜΟ雷达多目标检测仿真方法，其特征在于，假设 ΜΙΜΟ雷达模型由Μ个发射天线、Ν个接收天线组成，则第j个接收天线接收的信号yj表示为：式中，X康示第i个发射天线所发射的信号，η徒示第j个接收天线上的总噪声，hi康示 通道增益，j = 1，…，N，tpl为发射功率归一化因子； V Μ 假设发射信号正交，接收端进行匹配滤波并转换成矩阵Υ，则：式中，Υ二[yj ]ΜΝΧ1 ,Η二比ij]MN,n二[nj ]ΜΝΧ1。3. 根据权利要求2所述的一种统计ΜΙΜΟ雷达多目标检测仿真方法，其特征在于，利用功 率谱密度求解特征值λι和特征向量m的方法如下： Ai=Pyy(fi) (1)式中，fi = i/L，i = 0，l，…，L-l，Pyy()是功率谱密度函数，L是样本个数;e邱是W自然常 数e为底的指数函数。4. 根据权利要求3所述的一种统计ΜΙΜΟ雷达多目标检测仿真方法，其特征在于： 当目标观测通道完全独立时，当Ri- = (1/始)1,。, +诗，则共有ΜΝ个独立的目标观测通 道，有MN个相等的特征值4 =·^ +的; Μ 当目标观测通道不完全独立时又分为完全相关和部分相关两部分： 若目标观测通道完全相关，当R,- = (1 / +巧，则共有丽个特征值，其中一个特征 值为4 = W + σ,:，其余特征值为Λ = 若目标观测通道部分相关，则假设接收信号表示为y=n+j 丫，求得：通过构建I λ?ΜΝ-Κγ I = 0特征方程依次求得参数λ，对应为特征值λι，…，λΜΝ; 式中，Ry是接收信号协方差矩阵，Ιμν是ΜΝ维单位矩阵，的是噪声方差，η和丫是接收信号 的实部和虚部随机矢量，ΕΠ 是求解数学期望。5. 根据权利要求4所述的一种统计ΜΙΜΟ雷达多目标检测仿真方法，其特征在于，基于信 号和噪声特征值大小差异估计目标个数，设存在k个目标，假设噪声特征值多样性为q，则信 号的最小特征值λκ = λΜ-。应超过所设定的检测口限，而N-k个噪声的特征值Ak+l，…，λΜ应低于 检测口限，检测流程如下： 步骤a、首先从q= 1开始确定二元假设如下： Ho = k<N-q Hi = k = N-q 步骤b、利用信号特征值大于噪声方差特征值，构建基于特征值口限的循环表达式为：式中，AVq为噪声特征值的上界，册表示目标不存在，化表示目标存在;将求解的特征值 λι排序代入循环表达式，其中，若λι小于XUN-q，则满足化，则令q = q+l返回步骤a继续检验直 至满足化或者q=M;若λκ满足化，即λκ>λν。，则停止检验，此时检测到的目标数目k = N-q。
【文档编号】G01S7/41GK105842674SQ201610239715
【公开日】2016年8月10日
【申请日】2016年4月18日
【发明人】刘伟伟, 唐蕾, 沈卫康
【申请人】南京工程学院