CN105044683B - 多径条件下共形阵mimo雷达系统的发射波形设计方法 - Google Patents
多径条件下共形阵mimo雷达系统的发射波形设计方法 Download PDFInfo
- Publication number
- CN105044683B CN105044683B CN201510477718.4A CN201510477718A CN105044683B CN 105044683 B CN105044683 B CN 105044683B CN 201510477718 A CN201510477718 A CN 201510477718A CN 105044683 B CN105044683 B CN 105044683B
- Authority
- CN
- China
- Prior art keywords
- conformal
- radar system
- mimo radar
- conformal array
- matrix
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000013461 design Methods 0.000 title claims abstract description 24
- 238000000034 method Methods 0.000 title claims abstract description 18
- 239000011159 matrix material Substances 0.000 claims abstract description 71
- 238000001514 detection method Methods 0.000 claims abstract description 30
- 238000010586 diagram Methods 0.000 claims description 27
- 230000005540 biological transmission Effects 0.000 claims description 19
- 239000013598 vector Substances 0.000 claims description 18
- 230000005684 electric field Effects 0.000 claims description 13
- 230000005855 radiation Effects 0.000 claims description 9
- 238000005259 measurement Methods 0.000 claims description 6
- 230000021615 conjugation Effects 0.000 claims description 2
- 239000004576 sand Substances 0.000 claims description 2
- 238000005457 optimization Methods 0.000 description 6
- 238000004088 simulation Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 4
- 239000000969 carrier Substances 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 239000007787 solid Substances 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000000137 annealing Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000035899 viability Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
- G01S7/282—Transmitters
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明公开了一种多径条件下共形阵MIMO雷达系统的发射波形设计方法,包括以下步骤:(1)给定共形阵MIMO雷达系统的方位角θ的探测范围为[θmin,θmax],以及俯仰角的探测范围为设定共形阵天线的阵元个数N;(2)设定共形阵天线的第i个阵元包含的多径点个数为Mi;(3)设定共形阵MIMO雷达系统的期望的探测方向并设定共形阵MIMO雷达系统的期望的发射方向图在方位维的主瓣宽度wθ和俯仰维的主瓣宽度计算共形阵MIMO雷达系统的期望的发射方向图Bp;(4)根据共形阵MIMO雷达系统的期望的发射方向图Bp,构建共形阵MIMO雷达系统的波形矩阵S的设计准则;(5)根据所述设计准则,求解得到最终的共形阵MIMO雷达系统的波形矩阵
Description
技术领域
本发明属于雷达技术领域,特别涉及一种多径条件下共形阵MIMO雷达系统的发射波形设计方法,用于在多径条件下设计共形阵MIMO雷达的发射方向图。
背景技术
共形阵天线是指阵列表面与载体平台外形相吻合的阵列天线。与传统的平面阵天线相比,共形阵天线的特殊阵列流形使其各个阵元的朝向与阵元载体平台的曲面法向保持一致,从而使共形阵天线的波束扫描范围从平面阵天线通常能扫描到的±60°立体角的范围扩展到±90°立体角的范围,甚至扩展到全空间覆盖,可以显著提高雷达的视场范围;另外,共形阵天线能够与高速飞行载体,诸如战斗机、导弹、卫星等充分共形,而不会影响高速飞行载体的外形设计;因此,共形阵天线不仅能够大幅度降低飞行载体的雷达散射截面积(RCS),提高雷达在现代复杂电磁环境中的抗电磁干扰能力,而且还能满足飞机、导弹等飞行载体在高速飞行时的空气动力外形结构要求,能够降低飞行载体的空气阻力,极大提高了飞行载体的生存力与战斗力。鉴于共形阵天线的上述诸多优点,各国科研人员已经对共形阵雷达系统进行了大量研究。以色列研制的费尔康预警机曾被公认为是运用共形阵天线最成功的实例之一;以色列在此基础上研制了费尔康预警机的升级版-海雕预警机,该预警机同样采用了共形阵天线技术;美国的F-22和F-35战斗机也都使用了共形阵天线技术;王小谟院士在谈到我国的预警机发展趋势时表示:以后我国的预警机将采用共形阵天线技术。
但是,目前针对共形阵雷达系统的研究主要是针对共形相控阵雷达的研究,研究内容主要包括共形阵相控雷达的阵元布阵、自适应波束形成和共形阵相控雷达的空域信号处理等内容。由于共形阵天线主要安装在飞行载体上,而飞行载体一般都具有较大的飞行翼等突出结构,因此,共形阵天线在发射信号时会存在严重的多径问题。由于多径问题的存在,使共形阵天线的发射方向图会出现旁瓣升高和主瓣变形等问题,严重降低共形阵雷达系统的探测性能。虽然,共形相控阵雷达可以通过设计、发射波束形成导向矢量来缓解多径问题造成的影响,但是,效果仍然有限。
发明内容
针对上述现有技术的不足,本发明的目的在提出一种多径条件下共形阵MIMO雷达系统的发射波形设计方法,能够在多径条件下设计共形阵MIMO雷达的发射方向图。
实现本发明的技术思路是:首先,给定共形阵MIMO雷达系统的探测参数,并给定共形阵天线的阵元的多径点参数;然后,根据实际探测的需要,确定共形阵MIMO雷达系统的期望的发射方向图,并根据该发射方向图构建共形阵MIMO雷达系统的波形矩阵S的设计准则;最后根据该设计准则,使用优化算法求解得到最终的MIMO雷达系统的波形矩阵。
为实现上述技术目的,本发明采取以下技术方案予以实现。
一种多径条件下共形阵MIMO雷达系统的发射波形设计方法,其特征在于,包括以下具体步骤:
步骤1,给定共形阵MIMO雷达系统的方位角θ的探测范围为[θmin,θmax],以及俯仰角的探测范围为设定共形阵天线的阵元个数N,其中,第i个阵元相对于共形阵天线相位中心的相位补偿为第i个阵元的电场强度为i=1,2,…,N;所述共形阵MIMO雷达系统中包含所述共形阵天线;
步骤2,对于共形阵天线的第i个阵元,设定其包含的多径点个数为Mi,其中,第m个多径点的辐射强度为fi,m,第m个多径点的散射幅度为βi,m,第m个多径点相对于共形阵天线相位中心的相位补偿为i=1,2,…,N,m=1,2,…,Mi;
步骤3,设定共形阵MIMO雷达系统的期望的探测方向并设定共形阵MIMO雷达系统的期望的发射方向图在方位维的主瓣宽度wθ和俯仰维的主瓣宽度再根据步骤1设定的共形阵MIMO雷达系统的探测范围和阵元参数、步骤2设定的共形阵天线的阵元的多径点参数,计算共形阵MIMO雷达系统的期望的发射方向图Bp;
所述共形阵MIMO雷达系统的探测范围,包括共形阵MIMO雷达系统的方位角的探测范围以及俯仰角的探测范围;所述共形阵MIMO雷达系统的阵元参数,包括共形阵天线的阵元个数、每个阵元相对于共形阵天线相位中心的相位补偿以及每个阵元的电场强度;所述共形阵天线的阵元的多径点参数包括共形阵天线的每个阵元包含的多径点个数、每个多径点的辐射强度、每个多径点的散射幅度以及每个多径点相对于共形阵天线相位中心的相位补偿;
步骤4,根据共形阵MIMO雷达系统的期望的发射方向图Bp,构建共形阵MIMO雷达系统的波形矩阵S的设计准则;
步骤5,根据共形阵MIMO雷达系统的波形矩阵S的设计准则,使用优化算法求解得到最终的共形阵MIMO雷达系统的波形矩阵
本发明的有益效果为:本发明的多径条件下的共形阵MIMO雷达的波形设计方法,该方法在多径条件下通过优化共形阵MIMO雷达系统的波形来设计发射方向图,在优化过程中,相比传统的共形相控阵雷达系统,本发明方法中共形阵天线的每个阵元都可以设计独立的波形,因而具有更多的自由度,得到了旁瓣更低且主瓣指向性更好的发射方向图。
附图说明
下面结合附图说明和具体实施方式对本发明作进一步详细说明。
图1是本发明的流程图;
图2是序列二次规划算法求解最终的共形阵MIMO雷达系统的波形矩阵的流程图;
图3是仿真实验的共形阵MIMO雷达系统的共形阵天线的示意图;
图4是共形阵MIMO雷达系统的发射波形的发射方向图、优化发射波束形成导向向量得到的发射方向图、存在多径时共形天线的发射方向图以及没有多径时共形阵天线的发射方向图的的对比图,横坐标为方位角,单位为度(°),纵坐标为发射方向图。
具体实施方式
参照图1,本发明的多径条件下共形阵MIMO雷达系统的发射波形设计方法,包括以下具体步骤:
步骤1,给定共形阵MIMO雷达系统的方位角θ的探测范围为[θmin,θmax],以及俯仰角的探测范围为设定共形阵天线的阵元个数N,其中,第i个阵元相对于共形阵天线相位中心的相位补偿为第i个阵元的电场强度为i=1,2,…,N;所述共形阵MIMO雷达系统中包含所述共形阵天线。
步骤2,对于共形阵天线的第i个阵元,设定其包含的多径点个数为Mi,其中,第m个多径点的辐射强度为fi,m,第m个多径点的散射幅度为βi,m,第m个多径点相对于共形阵天线相位中心的相位补偿为i=1,2,…,N,m=1,2,…,Mi。
步骤3,设定共形阵MIMO雷达系统的期望的探测方向并设定共形阵MIMO雷达系统的期望的发射方向图在方位维的主瓣宽度wθ和俯仰维的主瓣宽度再根据步骤1设定的共形阵MIMO雷达系统的探测范围和阵元参数、步骤2设定的共形阵天线的阵元的多径点参数,计算共形阵MIMO雷达系统的期望的发射方向图Bp;
所述共形阵MIMO雷达系统的探测范围,包括共形阵MIMO雷达系统的方位角的探测范围以及俯仰角的探测范围;所述共形阵MIMO雷达系统的阵元参数,包括共形阵天线的阵元个数、每个阵元相对于共形阵天线相位中心的相位补偿以及每个阵元的电场强度;所述共形阵天线的阵元的多径点参数包括共形阵天线的每个阵元包含的多径点个数、每个多径点的辐射强度、每个多径点的散射幅度以及每个多径点相对于共形阵天线相位中心的相位补偿;
步骤3的具体子步骤为:
3.1将共形阵MIMO雷达系统的方位角θ的探测范围[θmin,θmax]等间隔均匀离散为Nθ维的方位角向量即将共形阵MIMO雷达系统的俯仰角的探测范围等间隔均匀离散为维的俯仰角向量即为保证获得共形阵MIMO雷达系统的期望的发射方向图,所述方位角向量中相邻元素的间隔(方位角离散角度间隔)、所述俯仰角向量中相邻元素的间隔(俯仰角离散角度间隔)均应小于等于1°,在本发明实例中,方位角离散角度间隔和俯仰角离散角度间隔均取0.5°;
3.2设定共形阵MIMO雷达系统的期望的探测方向θd和分别为共形阵MIMO雷达系统的期望的方向角和俯仰角,简称期望方向角和期望俯仰角;并设定共形阵MIMO雷达系统的期望的发射方向图Bp在方位维的主瓣宽度wθ和俯仰维的主瓣宽度
3.3根据共形阵天线的第i个阵元相对于共形阵天线相位中心的相位补偿以及共形阵天线的第i个阵元的电场强度计算没有多径时共形阵MIMO雷达系统的发射方向图B1,B1中位于第n1行第n2列的元素B1(n1,n2)为:
其中,n1=1,2,…,Nθ, 和分别表示共形阵MIMO雷达系统的方位角取且俯仰角取时,共形阵天线的第i个阵元的电场强度和第i个阵元相对于共形阵天线相位中心的相位补偿;表示期望方位角为θd且期望俯仰角为时,共形阵天线的第i个阵元相对于共形阵天线相位中心的相位补偿;i=1,2,…,N,N为共形阵天线的阵元个数;|·|表示取模,(·)*表示取共轭;
3.4如果以下两个不等式成立,则共形阵MIMO雷达系统的期望的发射方向图Bp中位于第n1行第n2列的元素Bp(n1,n2)=B1(n1,n2),否则,Bp(n1,n2)=0,所述两个不等式为:
步骤4,根据共形阵MIMO雷达系统的期望的发射方向图Bp,构建共形阵MIMO雷达系统的波形矩阵S的设计准则;
步骤4的具体子步骤为:
4.1设定共形阵MIMO雷达系统的发射波形的码元长度Ns,则共形阵MIMO雷达系统的波形矩阵S为Ns×N维的复矩阵,S=[s1,…,si,…,sN],其中,si为共形阵天线的第i个阵元的发射波形,i=1,2,…,N,N为共形阵天线的阵元个数,(·)T表示转置;
4.2计算存在多径时共形阵MIMO雷达系统的发射方向图Bs,Bs中位于第n1行第n2列的元素Bs(n1,n2)为:
其中,n1=1,2,…,Nθ,Nθ和分别为方向角向量和俯仰角向量的维数;和分别表示共形阵MIMO雷达系统的方位角取且俯仰角取时,共形阵天线的第i个阵元的电场强度和第i个阵元相对于共形阵天线相位中心的相位补偿;fi,m和βi,m分别表示共形阵天线的第i个阵元的第m个多径点的辐射强度和散射幅度;表示共形阵MIMO雷达系统的方位角取且俯仰角取时,共形阵天线的第i个阵元的第m个多径点相对于共形阵天线相位中心的相位补偿;sik为共形阵MIMO雷达系统的波形矩阵S中位于第k行第i列的元素,k=1,2,…,Ns,Ns为共形阵MIMO雷达系统的发射波形的码元长度;i=1,2,…,N,N为共形阵天线的阵元个数,m=1,2,…,Mi,Mi为共形阵天线的第i个阵元包含的多径点个数;|·|表示取模;
4.3根据共形阵MIMO雷达系统的期望的发射方向图Bp、存在多径时共形阵MIMO雷达系统的发射方向图Bs,构建共形阵MIMO雷达系统的波形矩阵S的设计准则为:
s.t.α>O
其中,Bp(n1,n2)为共形阵MIMO雷达系统的期望的发射方向图Bp中位于第n1行第n2列的元素,min表示最小化,max表示取最大值,s.t.表示约束条件;Φ为共形阵MIMO雷达系统的波形矩阵S的相位矩阵,S=exp(jΦ),j表示虚数单位,α为系数变量,n1=1,2,…,Nθ,|·|表示取模。
步骤5,根据共形阵MIMO雷达系统的波形矩阵S的设计准则,使用优化算法求解得到最终的共形阵MIMO雷达系统的波形矩阵
共形阵MIMO雷达系统的波形矩阵S的设计准则是一个极小极大优化问题,根据该设计准则,能够使用现有的优化算法,例如:模拟退火算法、遗传算法和序列二次规划算法等,求解得到最终的共形阵MIMO雷达系统的波形矩阵。本发明实例中,使用序列二次规划算法求解最终的共形阵MIMO雷达系统的波形矩阵,参照图2,具体步骤为:
5.1设置迭代的总次数Nc,设置最小误差值emin的初始值为无穷大,设置迭代次数nc=1,设置最优相位矩阵的初始值为全零矩阵;
5.2初始化共形阵MIMO雷达系统的第nc次迭代的波形矩阵的相位矩阵Φ(nc),并初始化第nc次迭代的系数变量α(nc),即将Φ(nc)中的每个元素设置为在0~2π间取值的一个随机值,将α(nc)设置为大于零的随机值;
5.3根据共形阵MIMO雷达系统的波形矩阵S的设计准则,运用序列二次规划算法对共形阵MIMO雷达系统的第nc次迭代的波形矩阵的相位矩阵Φ(nc)进行优化,求解得到第nc次迭代的优化后的相位矩阵Φ′(nc);
5.4比较第nc次迭代的优化后的相位矩阵Φ′(nc)的目标函数值与最小误差值emin的大小,如果小于emin,则令并将Φ′(nc)赋值给最优相位矩阵否则,忽略第nc次循环的优化后的相位矩阵Φ′(nc);
其中,为存在多径时共形阵MIMO雷达系统的第nc次迭代的发射方向图中位于第n1行第n2列的元素,的表达式为:
其中,n1=1,2,…,Nθ,Nθ和分别为方向角向量和俯仰角向量的维数;和分别表示共形阵MIMO雷达系统的方位角取且俯仰角取时,共形阵天线的第i个阵元的电场强度和第i个阵元相对于共形阵天线相位中心的相位补偿;fi,m和βi,m分别表示共形阵天线的第i个阵元的第m个多径点的辐射强度和散射幅度;表示共形阵MIMO雷达系统的方位角取且俯仰角取时,共形阵天线的第i个阵元的第m个多径点相对于共形阵天线相位中心的相位补偿;为共形阵MIMO雷达系统的第nc次迭代的波形矩阵中位于第k行第i列的元素,k=1,2,…,Ns,Ns为共形阵MIMO雷达系统的发射波形的码元长度;i=1,2,…,N,N为共形阵天线的阵元个数,m=1,2,…,Mi,Mi为共形阵天线的第i个阵元包含的多径点个数;|·|表示取模;
5.5若迭代次数nc小于Nc,令nc增加1,返回步骤5.2;否则,迭代结束,将最优相位矩阵作为最终的共形阵MIMO雷达系统的波形矩阵的相位矩阵,进而得到最终的共形阵MIMO雷达系统的波形矩阵
本发明的效果可通过以下仿真实验作进一步说明:
1)仿真条件:
共形阵MIMO雷达系统参数:参照图3,共形阵天线为圆弧阵,圆弧半径为57.3cm,包含8个阵元,每个阵元均采用微带天线,阵元弧长间距为5cm,阵元间的角度间隔为5°;共形阵MIMO雷达系统的载频为3GHz,波长为10cm,方位维的探测范围为[-72°,72°];共形阵天线的每个阵元包含两个多径点,每个多径点的位置为随机产生,每个多径点的散射幅度为在0~1间取值的随机值。
其他仿真参数:共形阵MIMO雷达系统的发射波形的码元长度为128,将方位维的探测范围[-72°,72°]以0.5°为间隔均匀离散,循环的总次数为10。
2)仿真内容及分析:分别绘制并比较以下四种发射方向图:共形阵MIMO雷达系统的发射波形的发射方向图,如图4中的实线所示;优化发射波束形成导向向量得到的发射方向图,如图4中的点线所示;存在多径时共形天线的发射方向图,如图4中的点划线所示;没有多径时共形阵天线的发射方向图,如图4中的虚线所示。
从图4中可以看出,与没有多径时共形阵天线的发射方向图相比,存在多径时共形阵天线的发射方向图的旁瓣明显升高;在存在多径的条件下,通过优化发射波束形成导向向量,可以降低发射方向图的旁瓣,但是该旁瓣仍然高于没有多径时共形阵天线的发射方向图的旁瓣,并主瓣位置发生了偏移;在存在多径的条件下,通过优化共形阵MIMO雷达系统的发射波形,可以显著降低发射方向图的旁瓣,并且该旁瓣远低于没有多径时共形阵天线的发射方向图的旁瓣,同时,保持主瓣位置不变,并使主瓣宽度变窄,主瓣指向性更好。
综上所述,在存在多径的条件下,运用本发明方法得到的共形阵MIMO雷达系统的发射波形,其发射方向图具有更低的旁瓣和更好的主瓣指向性,表明了本发明方法的有效性。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (2)
1.一种多径条件下共形阵MIMO雷达系统的发射波形设计方法,其特征在于,包括以下具体步骤:
步骤1,给定共形阵MIMO雷达系统的方位角θ的探测范围为[θmin,θmax],以及俯仰角的探测范围为设定共形阵天线的阵元个数N,其中,第i个阵元相对于共形阵天线相位中心的相位补偿为第i个阵元的电场强度为i=1,2,…,N;所述共形阵MIMO雷达系统中包含所述共形阵天线;其中,θmin表示共形阵MIMO雷达系统可探测的最小方位角;θmax表示共形阵MIMO雷达系统可探测的最大方位角;表示共形阵MIMO雷达系统可探测的最小俯仰角;表示共形阵MIMO雷达系统可探测的最大俯仰角;
步骤2,对于共形阵天线的第i个阵元,设定其包含的多径点个数为Mi,其中,第m个多径点的辐射强度为fi,m,第m个多径点的散射幅度为βi,m,第m个多径点相对于共形阵天线相位中心的相位补偿为i=1,2,…,N,m=1,2,…,Mi;
步骤3,设定共形阵MIMO雷达系统的期望的探测方向并设定共形阵MIMO雷达系统的期望的发射方向图在方位维的主瓣宽度wθ和俯仰维的主瓣宽度计算共形阵MIMO雷达系统的期望的发射方向图Bp;其中,θd和分别为共形阵MIMO雷达系统的期望的方位角和俯仰角,简称期望方位角和期望俯仰角;
其中,所述步骤3的具体子步骤为:
3.1将共形阵MIMO雷达系统的方位角θ的探测范围[θmin,θmax]等间隔均匀离散为Nθ维的方位角向量即将共形阵MIMO雷达系统的俯仰角的探测范围等间隔均匀离散为维的俯仰角向量即
3.2设定共形阵MIMO雷达系统的期望的探测方向并设定共形阵MIMO雷达系统的期望的发射方向图Bp在方位维的主瓣宽度wθ和俯仰维的主瓣宽度
3.3根据共形阵天线的第i个阵元相对于共形阵天线相位中心的相位补偿以及共形阵天线的第i个阵元的电场强度计算没有多径时共形阵MIMO雷达系统的发射方向图B1,B1中位于第n1行第n2列的元素B1(n1,n2)为:
其中,n1=1,2,…,Nθ, 和分别表示共形阵MIMO雷达系统的方位角取且俯仰角取时,共形阵天线的第i个阵元的电场强度和第i个阵元相对于共形阵天线相位中心的相位补偿;表示期望方位角为θd且期望俯仰角为时,共形阵天线的第i个阵元相对于共形阵天线相位中心的相位补偿;i=1,2,…,N,N为共形阵天线的阵元个数;|·|表示取模,(·)*表示取共轭;
3.4如果以下两个不等式成立,则共形阵MIMO雷达系统的期望的发射方向图Bp中位于第n1行第n2列的元素Bp(n1,n2)=B1(n1,n2),否则,Bp(n1,n2)=0,所述两个不等式为:
步骤4,根据共形阵MIMO雷达系统的期望的发射方向图Bp,构建共形阵MIMO雷达系统的波形矩阵S的设计准则;
其中,所述步骤4的具体子步骤为:
4.1设定共形阵MIMO雷达系统的发射波形的码元长度Ns,则共形阵MIMO雷达系统的波形矩阵S为Ns×N维的复矩阵,S=[s1,...,si,...,sN],其中,si为共形阵天线的第i个阵元的发射波形,N为共形阵天线的阵元个数,(·)T表示转置;sik为共形阵MIMO雷达系统的波形矩阵S中位于第k行第i列的元素,k=1,2,…,Ns;
4.2计算存在多径时共形阵MIMO雷达系统的发射方向图Bs,Bs中位于第n1行第n2列的元素Bs(n1,n2)为:
其中,n1=1,2,…,Nθ,Nθ和分别为方位角向量和俯仰角向量的维数;和分别表示共形阵MIMO雷达系统的方位角取且俯仰角取时,共形阵天线的第i个阵元的电场强度和第i个阵元相对于共形阵天线相位中心的相位补偿;fi,m和βi,m分别表示共形阵天线的第i个阵元的第m个多径点的辐射强度和散射幅度;表示共形阵MIMO雷达系统的方位角取且俯仰角取时,共形阵天线的第i个阵元的第m个多径点相对于共形阵天线相位中心的相位补偿;i=1,2,…,N,N为共形阵天线的阵元个数,m=1,2,…,Mi,Mi为共形阵天线的第i个阵元包含的多径点个数;|·|表示取模;
4.3根据共形阵MIMO雷达系统的期望的发射方向图Bp、存在多径时共形阵MIMO雷达系统的发射方向图Bs,构建共形阵MIMO雷达系统的波形矩阵S的设计准则为:
s.t.α>0
其中,Bp(n1,n2)为共形阵MIMO雷达系统的期望的发射方向图Bp中位于第n1行第n2列的元素,min表示最小化,max表示取最大值,s.t.表示约束条件;Φ为共形阵MIMO雷达系统的波形矩阵S的相位矩阵,S=exp(jΦ),j表示虚数单位,α为系数变量;
步骤5,根据共形阵MIMO雷达系统的波形矩阵S的设计准则,求解得到最终的共形阵MIMO雷达系统的波形矩阵。
2.如权利要求1所述的多径条件下共形阵MIMO雷达系统的发射波形设计方法,其特征在于,所述步骤5中,使用序列二次规划算法求解最终的共形阵MIMO雷达系统的波形矩阵,具体步骤为:
5.1设置迭代的总次数Nc,设置最小误差值emin的初始值为无穷大,设置迭代次数nc=1,设置最优相位矩阵的初始值为全零矩阵;
5.2初始化共形阵MIMO雷达系统的第nc次迭代的波形矩阵的相位矩阵Φ(nc),并初始化第nc次迭代的系数变量α(nc),即将Φ(nc)中的每个元素设置为在0~2π间取值的一个随机值,将α(nc)设置为大于零的随机值;
5.3根据共形阵MIMO雷达系统的波形矩阵S的设计准则,运用序列二次规划算法对共形阵MIMO雷达系统的第nc次迭代的波形矩阵的相位矩阵Φ(nc)进行优化,求解得到第nc次迭代的优化后的相位矩阵Φ′(nc);
5.4比较第nc次迭代的优化后的相位矩阵Φ′(nc)的目标函数值与最小误差值emin的大小,如果小于emin,则令并将Φ′(nc)赋值给最优相位矩阵否则,忽略第nc次迭代的优化后的相位矩阵Φ′(nc);
其中,Bp(n1,n2)为共形阵MIMO雷达系统的期望的发射方向图Bp中位于第n1行第n2列的元素;
其中,为存在多径时共形阵MIMO雷达系统的第nc次迭代的发射方向图中位于第n1行第n2列的元素,的表达式为:
其中,n1=1,2,…,Nθ,Nθ和分别为方位角向量和俯仰角向量的维数;和分别表示共形阵MIMO雷达系统的方位角取且俯仰角取时,共形阵天线的第i个阵元的电场强度和第i个阵元相对于共形阵天线相位中心的相位补偿;fi,m和βi,m分别表示共形阵天线的第i个阵元的第m个多径点的辐射强度和散射幅度;表示共形阵MIMO雷达系统的方位角取且俯仰角取时,共形阵天线的第i个阵元的第m个多径点相对于共形阵天线相位中心的相位补偿;为共形阵MIMO雷达系统的第nc次迭代的波形矩阵中位于第k行第i列的元素,k=1,2,…,Ns,Ns为共形阵MIMO雷达系统的发射波形的码元长度;i=1,2,…,N,N为共形阵天线的阵元个数,m=1,2,…,Mi,Mi为共形阵天线的第i个阵元包含的多径点个数;|·|表示取模;
5.5若迭代次数nc小于Nc,令nc增加1,返回步骤5.2;否则,迭代结束,将最优相位矩阵作为最终的共形阵MIMO雷达系统的波形矩阵的相位矩阵,进而得到最终的共形阵MIMO雷达系统的波形矩阵
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510477718.4A CN105044683B (zh) | 2015-08-06 | 2015-08-06 | 多径条件下共形阵mimo雷达系统的发射波形设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510477718.4A CN105044683B (zh) | 2015-08-06 | 2015-08-06 | 多径条件下共形阵mimo雷达系统的发射波形设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105044683A CN105044683A (zh) | 2015-11-11 |
CN105044683B true CN105044683B (zh) | 2017-05-24 |
Family
ID=54451360
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510477718.4A Active CN105044683B (zh) | 2015-08-06 | 2015-08-06 | 多径条件下共形阵mimo雷达系统的发射波形设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105044683B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105515627B (zh) * | 2015-12-07 | 2018-07-27 | 东南大学 | 一种大规模mimo检测方法及检测装置 |
CN106291474B (zh) * | 2016-08-22 | 2019-01-11 | 西安电子科技大学 | 基于圆柱阵的集中式mimo雷达波形优化方法 |
CN106646420B (zh) * | 2016-09-28 | 2019-08-06 | 西安电子科技大学 | 基于lfm信号的mimo雷达发射方向图设计方法 |
CN106842134B (zh) * | 2016-12-19 | 2019-07-09 | 西安电子科技大学 | 一种基于fft的平面阵mimo雷达波形设计方法 |
CN108459301B (zh) * | 2018-01-18 | 2021-10-08 | 西安电子科技大学 | 一种基于异构阵的mimo雷达波形设计方法 |
CN109507664B (zh) * | 2019-01-22 | 2020-05-22 | 中国人民解放军空军工程大学 | 压缩感知mimo雷达认知波形获取方法和装置 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101950014A (zh) * | 2010-09-02 | 2011-01-19 | 西安电子科技大学 | Mimo雷达波形在线设计方法 |
CN103076596A (zh) * | 2013-01-17 | 2013-05-01 | 西安电子科技大学 | 基于先验信息的mimo雷达发射方向图设计方法 |
CN103267963A (zh) * | 2013-05-07 | 2013-08-28 | 西安电子科技大学 | 基于相位调制的多准则雷达发射信号设计方法 |
CN103592642A (zh) * | 2013-10-22 | 2014-02-19 | 西安电子科技大学 | Mimo雷达波形的设计方法 |
CN103605113A (zh) * | 2013-12-03 | 2014-02-26 | 西安电子科技大学 | 多发多收干涉合成孔径雷达空时二维信号波形设计方法 |
CN103852751A (zh) * | 2014-03-26 | 2014-06-11 | 西安电子科技大学 | 基于接收波束形成的集中式mimo雷达波形设计方法 |
CN104076334A (zh) * | 2014-07-08 | 2014-10-01 | 西安电子科技大学 | 一种设计mimo雷达波形和发射天线阵列的方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006053161A (ja) * | 2005-10-17 | 2006-02-23 | Toshiba Corp | パルスデータ生成方法、シェイプデータ生成方法、パルスデータ生成装置、シェイプデータ生成装置及び送信パルス信号生成装置 |
JP2014115160A (ja) * | 2012-12-07 | 2014-06-26 | Furuno Electric Co Ltd | 探知装置、探知方法及び探知プログラム |
-
2015
- 2015-08-06 CN CN201510477718.4A patent/CN105044683B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101950014A (zh) * | 2010-09-02 | 2011-01-19 | 西安电子科技大学 | Mimo雷达波形在线设计方法 |
CN103076596A (zh) * | 2013-01-17 | 2013-05-01 | 西安电子科技大学 | 基于先验信息的mimo雷达发射方向图设计方法 |
CN103267963A (zh) * | 2013-05-07 | 2013-08-28 | 西安电子科技大学 | 基于相位调制的多准则雷达发射信号设计方法 |
CN103592642A (zh) * | 2013-10-22 | 2014-02-19 | 西安电子科技大学 | Mimo雷达波形的设计方法 |
CN103605113A (zh) * | 2013-12-03 | 2014-02-26 | 西安电子科技大学 | 多发多收干涉合成孔径雷达空时二维信号波形设计方法 |
CN103852751A (zh) * | 2014-03-26 | 2014-06-11 | 西安电子科技大学 | 基于接收波束形成的集中式mimo雷达波形设计方法 |
CN104076334A (zh) * | 2014-07-08 | 2014-10-01 | 西安电子科技大学 | 一种设计mimo雷达波形和发射天线阵列的方法 |
Non-Patent Citations (3)
Title |
---|
一种MIMO雷达多模式波形优化设计方法;王旭等;《电子与信息学报》;20150630;第37卷(第6期);1416-1423 * |
一种基于先验信息的MIMO雷达发射方向图设计方法;王旭等;《电子与信息学报》;20131231;第35卷(第12期);2802-2808 * |
低旁瓣MIMO雷达发射方向图设计;罗涛等;《电子与信息学报》;20131231;第35卷(第12期);2815-2822 * |
Also Published As
Publication number | Publication date |
---|---|
CN105044683A (zh) | 2015-11-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105044683B (zh) | 多径条件下共形阵mimo雷达系统的发射波形设计方法 | |
CN106291474B (zh) | 基于圆柱阵的集中式mimo雷达波形优化方法 | |
CN101369014B (zh) | 应用于多输入多输出雷达的双边约束自适应波束形成方法 | |
CN106772260B (zh) | 基于凸优化算法的雷达阵列和差波束方向图优化方法 | |
CN105785328B (zh) | 基于子阵划分的fda距离-角度解耦合波束形成方法 | |
CN103020363B (zh) | 一种通过改善阵列波束方向图旁瓣特性设计天线的方法 | |
CN105044684B (zh) | 基于射频隐身的mimo跟踪雷达发射波束的形成方法 | |
CN102565790A (zh) | 平面相控阵自适应的和差测角方法 | |
CN111400919B (zh) | 阵列天线中的低旁瓣波束设计方法 | |
CN103076596A (zh) | 基于先验信息的mimo雷达发射方向图设计方法 | |
CN102608580A (zh) | 数字阵列超低副瓣自适应数字波束形成方法 | |
CN106646420B (zh) | 基于lfm信号的mimo雷达发射方向图设计方法 | |
CN106646388A (zh) | 基于嵌套阵列的mimo雷达抗干扰方法 | |
CN104345299A (zh) | 基于简化ec的机载mimo雷达空时自适应处理方法 | |
CN105372633A (zh) | 一种相控阵雷达降维四通道抗主副瓣干扰的方法 | |
CN114488027A (zh) | 一种共形阵天线的波束调零展宽方法 | |
CN103513225B (zh) | 一种基于空间增益的稀疏平面阵形优化方法 | |
CN106646421B (zh) | 基于三维非均匀阵列的mimo雷达波形联合设计方法 | |
CN104678362A (zh) | Mimo天波超视距雷达波形优化方法 | |
CN110850383B (zh) | 一种基于共形阵的mimo雷达信号处理方法 | |
CN109633600B (zh) | 一种最小冗余线阵mimo-othr的doa估计方法 | |
CN113985358B (zh) | 一种基于相控阵mimo的接收adbf处理方法 | |
CN115015846B (zh) | 一种基于共形阵天线的广义自适应单脉冲测角方法 | |
CN109459729A (zh) | 一种基于恒模约束的mimo雷达非对称方向图设计方法 | |
Finden et al. | A dual polarized crossed bowtie array antenna for DVB-T based passive radar |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |