• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
    • / 21
    • 下载费用:30 金币  

    重庆时时彩输钱了咋办: 基于凸优化算法的雷达阵列和差波束方向图优化方法.pdf

    关 键 词:
    基于 优化 算法 雷达 阵列 波束 方向 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201710206959.4

    申请日:

    2017.03.31

    公开号:

    CN106772260A

    公开日:

    2017.05.31

    当前法律状态:

    实审

    有效性:

    审中

    法律详情: 实质审查的生效IPC(主分类):G01S 7/02申请日:20170331|||公开
    IPC分类号: G01S7/02 主分类号: G01S7/02
    申请人: 西安电子科技大学
    发明人: 曹运合; 张奕; 凤宏晓; 苏洪涛; 周生华
    地址: 710071 陕西省西安市太白南路2号
    优先权:
    专利代理机构: 西安睿通知识产权代理事务所(特殊普通合伙) 61218 代理人: 惠文轩
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201710206959.4

    授权公告号:

    |||

    法律状态公告日:

    2017.06.23|||2017.05.31

    法律状态类型:

    实质审查的生效|||公开

    摘要

    本发明公开了一种基于凸优化算法的雷达阵列和差波束方向图优化方法,思路为:确定包含M个阵元的雷达圆形阵列均匀分布在半径为R的圆周上,然后确定和波束主瓣区域、和波束零陷区域、和波束旁瓣区域、和波束低旁瓣区域,以及差波束主瓣区域、差波束零陷区域、差波束旁瓣区域、差波束低旁瓣区域;分别得到约束和波束主瓣区域的代价函数和约束差波束主瓣区域的代价函数;进而分别得到雷达圆形阵列的和波束方向图综合优化模型和雷达圆形阵列的差波束方向图综合优化模型,计算雷达圆形阵列和波束最优权矢量;依次得到期望目标的角误差归一化曲线的斜率和待优化的雷达圆形阵列差和差波束方向图的角误差曲线斜率,并计算雷达圆形阵列差波束的最优权矢量。

    权利要求书

    1.一种基于凸优化算法的雷达阵列和差波束方向图优化方法,其特征在于,包括以下
    步骤:
    步骤1,确定包含M个阵元的雷达圆形阵列均匀分布在半径为R的圆周上,然后计算得到
    雷达圆形阵列的方向矢量;M、R分别为大于0的正整数;
    步骤2,根据雷达圆形阵列的方向矢量,分别计算得到雷达圆形阵列的和波束方向图和
    雷达圆形阵列的差波束方向图,进而分别得到和波束主瓣区域、和波束零陷区域、和波束旁
    瓣区域、和波束低旁瓣区域,以及差波束主瓣区域、差波束零陷区域、差波束旁瓣区域、差波
    束低旁瓣区域;
    步骤3,分别得到约束和波束主瓣区域的代价函数和约束差波束主瓣区域的代价函数;
    步骤4,根据约束和波束主瓣区域的代价函数和约束差波束主瓣区域的代价函数,分别
    得到雷达圆形阵列的和波束方向图综合优化模型和雷达圆形阵列的差波束方向图综合优
    化模型;
    步骤5,根据雷达圆形阵列的和波束方向图综合优化模型,计算得到雷达圆形阵列和波
    束最优权矢量;
    步骤6,将雷达圆形阵列包含的M个阵元各自方向分别指向第一设定方向和第二设定方
    向,所述第一设定方向和第二设定方向为不同的方向;进而分别形成两个波束方向图,记为
    第一波束方向图和第二波束方向图,然后分别计算得到雷达圆形阵列的原始和波束接收能
    量和雷达圆形阵列的原始差波束接收能量;
    步骤7,根据雷达圆形阵列的原始和波束接收能量和雷达圆形阵列的原始差波束接收
    能量,计算得到期望目标的角误差归一化曲线,并拟合该期望目标的角误差归一化曲线得
    到期望目标的角误差归一化曲线的斜率;
    步骤8,根据雷达圆形阵列和波束最优权矢量,计算得到待优化的雷达圆形阵列差和差
    波束方向图的角误差曲线,并拟合得到待优化的雷达圆形阵列差和差波束方向图的角误差
    曲线斜率;
    步骤9,根据雷达圆形阵列的差波束方向图综合优化模型和待优化的雷达圆形阵列差
    和差波束方向图的角误差曲线斜率,计算得到雷达圆形阵列差波束的最优权矢量。
    2.如权利要求1所述的一种基于凸优化算法的雷达阵列和差波束方向图优化方法,其
    特征在于,在步骤1中,所述确定包含M个阵元的雷达圆形阵列均匀分布在半径为R的圆周
    上,还包括:
    确定包含M个阵元的雷达圆形阵列按照阵元序号依次均匀分布在半径为R的圆周上,其
    中M个阵元分别为无方向性;并且距离雷达圆形阵列的设定位置处存在信号源,其中设定位
    置为到雷达圆形阵列的距离大于圆周半径R,所述信号源向雷达圆形阵列发射信号,分别将
    信号入射进雷达圆形阵列的方位角记为信号源的方位角θ,将信号入射进雷达圆形阵列的
    俯仰角记为信号源的俯仰角φ;同时确定雷达圆形阵列检测范围内存在期望目标,该期望
    目标的方位角为δ,期望目标的俯仰角为β;
    所述雷达圆形阵列的方向矢量a(θ,φ),其表达式为:
    <mrow> <mi>a</mi> <mrow> <mo>(</mo> <mi>&theta;</mi> <mo>,</mo> <mi>&phi;</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mi>exp</mi> <mrow> <mo>(</mo> <mi>j</mi> <mn>2</mn> <mi>&pi;</mi> <mi>R</mi> <mi> </mi> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&phi;</mi> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mo>(</mo> <mi>&theta;</mi> <mo>-</mo> <msub> <mi>&gamma;</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mo>/</mo> <mi>&lambda;</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>exp</mi> <mrow> <mo>(</mo> <mi>j</mi> <mn>2</mn> <mi>&pi;</mi> <mi>R</mi> <mi> </mi> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&phi;</mi> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mo>(</mo> <mi>&theta;</mi> <mo>-</mo> <msub> <mi>&gamma;</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>/</mo> <mi>&lambda;</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>exp</mi> <mrow> <mo>(</mo> <mi>j</mi> <mn>2</mn> <mi>&pi;</mi> <mi>R</mi> <mi> </mi> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&phi;</mi> <mi>cos</mi> <mo>(</mo> <mi>&theta;</mi> <mo>-</mo> <msub> <mi>&gamma;</mi> <mrow> <mi>M</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> <mo>/</mo> <mi>&lambda;</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
    其中,θ表示信号源的方位角,φ表示信号源的俯仰角;令γm表示第m个阵元到第1个阵
    元的角度,γm=2πm/M,m∈{0,1,…,M-1},λ表示信号源向雷达圆形阵列发射信号的波长,R
    表示圆周的半径,exp表示指数函数操作,j为虚数单位,sin表示求正弦操作,cos表示求余
    弦操作,M表示雷达圆形阵列包含的阵元个数,M、R分别为大于0的正整数。
    3.如权利要求1所述的一种基于凸优化算法的雷达阵列和差波束方向图优化方法,其
    特征在于,步骤2的子步骤为:
    2a)设定待优化的雷达圆形阵列和波束权矢量为wH,
    wH=[wH1,…,wHm,…,wHM]T,上标T表示转置,m∈{1,2,…,M},wHm表示待优化的雷达圆形
    阵列和波束权矢量内第m个阵元的和波束权值,M表示雷达圆形阵列包含的阵元个数,与待
    优化的雷达圆形阵列和波束权矢量内包含的阵元个数相同;并根据待优化的雷达圆形阵列
    和波束权矢量wH和雷达圆形阵列的方向矢量a(θ,φ),计算得到雷达圆形阵列的和波束方
    向图FH(θ),其表达式为:
    FH(θ)=wH×a(θ)
    其中,a(θ)表示雷达圆形阵列的方向矢量a(θ,φ)的方位维矢量;
    2b)设定待优化的雷达圆形阵列差波束权矢量为wD,
    wD=[wD1,…,wDm,…,wDM]T,wDm表示待优化的雷达圆形阵列差波束权矢量内第m个阵元的
    差波束权值,m∈{1,2,…,M},M表示雷达圆形阵列包含的阵元个数,与待优化的雷达圆形阵
    列差波束权矢量内包含的阵元个数相同;利用待优化的雷达圆形阵列差波束权矢量wD和雷
    达圆形阵列的方向矢量a(θ,φ),计算得到雷达圆形阵列的差波束方向图FD(θ),其表达式
    为:
    FD(θ)=wD×a(θ)
    2c)将雷达圆形阵列的和波束方向图FH(θ)和雷达圆形阵列的差波束方向图FD(θ)分别
    进行区域划分,将雷达圆形阵列的和波束方向图FH(θ)划分为和波束主瓣区域ΘHmain、和波
    束零陷区域ΘHnull、和波束旁瓣区域ΘHside、和波束低旁瓣区域ΘHlow;所述和波束主瓣区域
    ΘHmain为雷达圆形阵列的和波束方向图FH(θ)内接收功率在-13dB到0dB之间的区域,所述和
    波束零陷区域ΘHnull为雷达圆形阵列的和波束方向图FH(θ)内接收功率小于-80dB的区域,
    所述和波束旁瓣区域ΘHside为雷达圆形阵列的和波束方向图FH(θ)内接收功率小于-30dB的
    区域,所述和波束低旁瓣区域ΘHlow为雷达圆形阵列的和波束方向图FH(θ)内接收功率小于-
    40dB的区域;其中和波束主瓣区域ΘHmain内包含待优化的和波束的主瓣BH;
    将雷达圆形阵列的差波束方向图FD(θ)划分为:差波束主瓣区域ΘDmain、差波束零陷区
    域ΘDnull、差波束旁瓣区域ΘDside、差波束低旁瓣区域ΘDlow,所述差波束主瓣区域ΘDmain为
    雷达圆形阵列的差波束方向图FD(θ)内接收功率在-13dB到0dB之间的区域,所述差波束零
    陷区域ΘDnull为雷达圆形阵列的差波束方向图FD(θ)内接收功率小于-80dB的区域,所述差
    波束低旁瓣区域ΘDlow为雷达圆形阵列的差波束方向图FD(θ)内接收功率小于-30dB的区域,
    所述差波束低旁瓣区域ΘDlow为雷达圆形阵列的差波束方向图FD(θ)内接收功率小于-40dB
    的区域;其中差波束主瓣区域ΘDmain内包含待优化的差波束主瓣BD。
    4.如权利要求3所述的一种基于凸优化算法的雷达阵列和差波束方向图优化方法,其
    特征在于,步骤3的子步骤为:
    3a)分别设置雷达圆形阵列的和波束期望主瓣BH0,并在和主瓣区域ΘHmain内将待优化的
    和波束的主瓣BH和雷达圆形阵列的和波束期望主瓣BH0之差的向量2-范数最小化,进而得到
    约束和波束主瓣区域ΘHmain的代价函数:
    min||Fh(θHm”)-FH(θHm')||2θHm'∈ΘHmain m'=1,2,…,M'
    其中,m'∈{1,2,…,M'},M'表示和主瓣区域ΘHmain对应的方位角区间内包含的采样点
    总个数,θHm'表示和主瓣区域ΘHmain对应的方位角区间内第m'个采样点的有限近似采样角
    度,θHm'∈ΘHmain,Fh(θHm”)为雷达圆形阵列的和波束期望主瓣BH0内第m”个采样点有限近似采
    样角度对应的方向图,FH(θHm')为和主瓣区域ΘHmain对应的方位角区间内第m'个采样点的有
    限近似采样角度对应方向图,||·||2为向量2-范数,∈表示属于,min为求最小值操作,m”
    ∈{1,2,…,M”},M”为雷达圆形阵列的和波束期望主瓣BH0包含的采样点总个数;雷达圆形阵
    列的和波束期望主瓣BH0包含的采样点总个数与主瓣区域Θmain对应的方位角区间内包含的
    采样点总个数取值相等;
    3b)设置雷达圆形阵列的差波束期望主瓣BD0,并在差波束主瓣区域ΘDmain内将待优化的
    差波束主瓣BD与雷达圆形阵列的差波束期望主瓣BD0之差的向量2-范数最小化,得到约束差
    波束主瓣区域ΘDmain的代价函数:
    min||Fd(θDn')-FD(θDn)||2 θDn∈Θmain n=1,2,…,N
    其中,θDn∈ΘDmain,n∈{1,2,…,N},θDn表示差主瓣区域ΘDmain对应的方位角区间内第n
    个采样点的有限近似采样角度,N为差主瓣区域ΘDmain对应的方位角区间内包含的采样点总
    个数,FD(θDn)为差主瓣区域ΘDmain对应的方位角区间内第n个采样点的有限近似采样角度对
    应的方向图,Fd(θDn')为雷达圆形阵列的差波束期望主瓣BD0内第n'个采样点有限近似采样
    角度对应的方向图,
    n'∈{1,2,…,N'},N'为雷达圆形阵列的差波束期望主瓣BD0内包含的采样点总个数;雷
    达圆形阵列的差波束期望主瓣BD0内包含的采样点总个数与差主瓣区域ΘDmain对应的方位
    角区间内包含的采样点总个数取值相等,∈表示属于。
    5.如权利要求4所述的一种基于凸优化算法的雷达阵列和差波束方向图优化方法,其
    特征在于,步骤4的子步骤为:
    4a)根据约束和波束主瓣区域ΘHmain的代价函对雷达圆形阵列的和波束期望主瓣BH0进
    行∞-范数约束,进而得到雷达圆形阵列的和波束方向图综合优化模型:
    <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mi>min</mi> </mtd> <mtd> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>F</mi> <mi>h</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>&theta;</mi> <mrow> <msup> <mi>Hm</mi> <mrow> <mo>&prime;</mo> <mo>&prime;</mo> </mrow> </msup> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>F</mi> <mi>H</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>&theta;</mi> <mrow> <msup> <mi>Hm</mi> <mo>&prime;</mo> </msup> </mrow> </msub> <mo>)</mo> </mrow> <mo>|</mo> <msub> <mo>|</mo> <mn>2</mn> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>&theta;</mi> <mrow> <msup> <mi>Hm</mi> <mo>&prime;</mo> </msup> </mrow> </msub> <mo>&Element;</mo> <msub> <mi>&Theta;</mi> <mrow> <mi>H</mi> <mi>m</mi> <mi>a</mi> <mi>i</mi> <mi>n</mi> </mrow> </msub> </mrow> </mtd> <mtd> <mrow> <msup> <mi>m</mi> <mo>&prime;</mo> </msup> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mn>...</mn> <mo>,</mo> <msup> <mi>M</mi> <mo>&prime;</mo> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>s</mi> <mo>.</mo> <mi>t</mi> <mo>.</mo> </mrow> </mtd> <mtd> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>F</mi> <mi>h</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>&theta;</mi> <mrow> <mi>H</mi> <mi>s</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>|</mo> <msub> <mo>|</mo> <mi>&infin;</mi> </msub> <mo>&le;</mo> <msub> <mi>&xi;</mi> <mi>H</mi> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>&theta;</mi> <mrow> <mi>H</mi> <mi>s</mi> </mrow> </msub> <mo>&Element;</mo> <msub> <mi>&Theta;</mi> <mrow> <mi>H</mi> <mi>s</mi> <mi>i</mi> <mi>d</mi> <mi>e</mi> </mrow> </msub> </mrow> </mtd> <mtd> <mrow> <mi>s</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mn>...</mn> <mo>,</mo> <mi>S</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow></mrow> </mtd> <mtd> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>F</mi> <mi>h</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>&theta;</mi> <mrow> <mi>H</mi> <mi>e</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>|</mo> <msub> <mo>|</mo> <mi>&infin;</mi> </msub> <mo>&le;</mo> <msub> <mi>&epsiv;</mi> <mi>H</mi> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>&theta;</mi> <mrow> <mi>H</mi> <mi>e</mi> </mrow> </msub> <mo>&Element;</mo> <msub> <mi>&Theta;</mi> <mrow> <mi>H</mi> <mi>n</mi> <mi>u</mi> <mi>l</mi> <mi>l</mi> </mrow> </msub> </mrow> </mtd> <mtd> <mrow> <mi>e</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mn>...</mn> <mo>,</mo> <mi>E</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow></mrow> </mtd> <mtd> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>F</mi> <mi>h</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>&theta;</mi> <mrow> <mi>H</mi> <mi>l</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>|</mo> <msub> <mo>|</mo> <mi>&infin;</mi> </msub> <mo>&le;</mo> <msub> <mi>&eta;</mi> <mi>H</mi> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>&theta;</mi> <mrow> <mi>H</mi> <mi>l</mi> </mrow> </msub> <mo>&Element;</mo> <msub> <mi>&Theta;</mi> <mrow> <mi>H</mi> <mi>l</mi> <mi>o</mi> <mi>w</mi> </mrow> </msub> </mrow> </mtd> <mtd> <mrow> <mi>l</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mn>...</mn> <mo>,</mo> <mi>L</mi> </mrow> </mtd> </mtr> </mtable> </mfenced>
    其中,m'∈{1,2,…,M'},M'表示和主瓣区域ΘHmain对应的方位角区间内包含的采样点
    总个数,θHm'表示和主瓣区域ΘHmain对应的方位角区间内第m'个采样点的有限近似采样角
    度,Fh(θHm”)为雷达圆形阵列的和波束期望主瓣BH0内第m”个采样点有限近似采样角度对应
    的方向图,FH(θHm')为和主瓣区域ΘHmain对应的方位角区间内第m'个采样点的有限近似采样
    角度对应方向图,∈表示属于,m”∈{1,2,…,M”},M”为雷达圆形阵列的和波束期望主瓣BH0
    包含的采样点总个数,ξH为雷达圆形阵列的和波束方向图预先设定的最大旁瓣电平,εH为雷
    达圆形阵列的和波束方向图预先设定的零陷方向最大电平,ηH为雷达圆形阵列的和波束方
    向图预先设定的低旁瓣区域最大预设电平,θHs表示和波束旁瓣区域ΘHside对应的方位角区
    间内第s个采样点的有限近似采样角度,θHs∈ΘHside,s∈{1,2,…,S},S表示和波束旁瓣区
    域ΘHside对应的方位角区间内包含的采样点总个数;θHe表示和波束零陷区域ΘHnull对应的
    方位角区间内第e个采样点的有限近似采样角度,θHe∈ΘHnull,e∈{1,2,…,E},E表示和波
    束零陷区域ΘHnull对应的方位角区间内包含的采样点总个数;θHl表示和波束低旁瓣区域
    ΘHlow对应的方位角区间内第l个采样点的有限近似采样角度,θHl∈ΘHlow,l∈{1,2,…,L},
    L表示和波束低旁瓣区域ΘHlow对应的方位角区间内包含的采样点总个数;Fh(θHs)为雷达圆
    形阵列的和波束期望主瓣BH0内第s个采样点的有限近似采样角度对应的方向图,Fh(θHe)为
    雷达圆形阵列的和波束期望主瓣BH0内第e个采样点的有限近似采样角度对应的方向图,Fh
    (θHl)为雷达圆形阵列的和波束期望主瓣BH0内第l个采样点的有限近似采样角度对应的方
    向图,||·||∞为向量无穷∞-范数操作,s.t.表示约束条件,∈表示属于;
    4b)根据约束差波束主瓣区域ΘDmain的代价函数对雷达圆形阵列的差波束期望主瓣BD0
    进行∞-范数约束,得到雷达圆形阵列的差波束方向图综合优化模型:
    <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mi>min</mi> </mtd> <mtd> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>F</mi> <mi>d</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>&theta;</mi> <mrow> <msup> <mi>Dn</mi> <mo>&prime;</mo> </msup> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>F</mi> <mi>D</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>&theta;</mi> <mrow> <mi>D</mi> <mi>n</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>|</mo> <msub> <mo>|</mo> <mn>2</mn> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>&theta;</mi> <mrow> <mi>D</mi> <mi>n</mi> </mrow> </msub> <mo>&Element;</mo> <msub> <mi>&Theta;</mi> <mrow> <mi>D</mi> <mi>m</mi> <mi>a</mi> <mi>i</mi> <mi>n</mi> </mrow> </msub> </mrow> </mtd> <mtd> <mrow> <mi>n</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mn>...</mn> <mo>,</mo> <mi>N</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>s</mi> <mo>.</mo> <mi>t</mi> <mo>.</mo> </mrow> </mtd> <mtd> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>F</mi> <mi>d</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>&theta;</mi> <mrow> <msup> <mi>Ds</mi> <mo>&prime;</mo> </msup> </mrow> </msub> <mo>)</mo> </mrow> <mo>|</mo> <msub> <mo>|</mo> <mi>&infin;</mi> </msub> <mo>&le;</mo> <msub> <mi>&xi;</mi> <mi>D</mi> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>&theta;</mi> <mrow> <msup> <mi>Ds</mi> <mo>&prime;</mo> </msup> </mrow> </msub> <mo>&Element;</mo> <msub> <mi>&Theta;</mi> <mrow> <mi>D</mi> <mi>s</mi> <mi>i</mi> <mi>d</mi> <mi>e</mi> </mrow> </msub> </mrow> </mtd> <mtd> <mrow> <msup> <mi>s</mi> <mo>&prime;</mo> </msup> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mn>...</mn> <mo>,</mo> <msup> <mi>S</mi> <mo>&prime;</mo> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow></mrow> </mtd> <mtd> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>F</mi> <mi>d</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>&theta;</mi> <mrow> <msup> <mi>De</mi> <mo>&prime;</mo> </msup> </mrow> </msub> <mo>)</mo> </mrow> <mo>|</mo> <msub> <mo>|</mo> <mi>&infin;</mi> </msub> <mo>&le;</mo> <msub> <mi>&epsiv;</mi> <mi>D</mi> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>&theta;</mi> <mrow> <msup> <mi>De</mi> <mo>&prime;</mo> </msup> </mrow> </msub> <mo>&Element;</mo> <msub> <mi>&Theta;</mi> <mrow> <mi>D</mi> <mi>n</mi> <mi>u</mi> <mi>l</mi> <mi>l</mi> </mrow> </msub> </mrow> </mtd> <mtd> <mrow> <msup> <mi>e</mi> <mo>&prime;</mo> </msup> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mn>...</mn> <mo>,</mo> <msup> <mi>E</mi> <mo>&prime;</mo> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow></mrow> </mtd> <mtd> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>F</mi> <mi>d</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>&theta;</mi> <mrow> <msup> <mi>Dl</mi> <mo>&prime;</mo> </msup> </mrow> </msub> <mo>)</mo> </mrow> <mo>|</mo> <msub> <mo>|</mo> <mi>&infin;</mi> </msub> <mo>&le;</mo> <msub> <mi>&eta;</mi> <mi>D</mi> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>&theta;</mi> <mrow> <msup> <mi>Dl</mi> <mo>&prime;</mo> </msup> </mrow> </msub> <mo>&Element;</mo> <msub> <mi>&Theta;</mi> <mrow> <mi>D</mi> <mi>l</mi> <mi>o</mi> <mi>w</mi> </mrow> </msub> </mrow> </mtd> <mtd> <mrow> <msup> <mi>l</mi> <mo>&prime;</mo> </msup> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mn>...</mn> <mo>,</mo> <msup> <mi>L</mi> <mo>&prime;</mo> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow></mrow> </mtd> <mtd> <mrow> <msubsup> <mi>w</mi> <mi>D</mi> <mi>T</mi> </msubsup> <mo>&times;</mo> <mi>A</mi> <mrow> <mo>(</mo> <msub> <mi>&theta;</mi> <mi>B</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mo>=</mo> <mn>0</mn> </mrow> </mtd> <mtd> <mrow> <msub> <mi>&theta;</mi> <mi>B</mi> </msub> <mo>=</mo> <mi>&delta;</mi> </mrow> </mtd> <mtd> <mrow></mrow> </mtd> </mtr> </mtable> </mfenced>
    其中,n∈{1,2,…,N},θDn表示差主瓣区域ΘDmain对应的方位角区间内第n个采样点的有
    限近似采样角度,N为差主瓣区域ΘDmain对应的方位角区间内包含的采样点总个数,FD(θDn)
    为差主瓣区域ΘDmain对应的方位角区间内第n个采样点的有限近似采样角度对应的方向图,
    Fd(θDn')为雷达圆形阵列的差波束期望主瓣BD0内第n'个采样点有限近似采样角度对应的方
    向图,||·||2为向量2-范数,min为求最小值操作,n'∈{1,2,…,N'},N'为雷达圆形阵列的
    差波束期望主瓣BD0内包含的采样点总个数;ξD为雷达圆形阵列的差波束方向图预先设定的
    最大旁瓣电平,εD为雷达圆形阵列的差波束方向图预先设定的零陷方向最大电平,ηD为雷达
    圆形阵列的差波束方向图预先设定的低旁瓣区域最大预设电平,θDs'表示差波束旁瓣区域
    ΘDside对应的方位角区间内第s'个采样点的有限近似采样角度,θDs'∈ΘDside,s'∈{1,
    2,…,S'},S'表示差波束旁瓣区域ΘDside对应的方位角区间内内包含的采样点总个数;θDe'
    表示差波束零陷区域ΘDnull对应的方位角区间内第e'个采样点的有限近似采样角度,θDe'∈
    ΘDnull,e'∈{1,2,…,E'},E'表示差波束零陷区域ΘDnull对应的方位角区间内包含的采样
    点总个数;θDl'表示差波束低旁瓣区域ΘDlow对应的方位角区间内第l'个采样点的有限近似
    采样角度,θDl'∈ΘDlow,l'∈{1,2,…,L'},L'表示差波束低旁瓣区域ΘDlow对应的方位角区
    间内包含的采样点总个数;Fd(θDs')为雷达圆形阵列的差波束期望主瓣BD0内第s'个采样点
    的有限近似采样角度对应的方向图,Fd(θDe')为雷达圆形阵列的差波束期望主瓣BD0内第e'
    个采样点的有限近似采样角度对应的方向图,Fd(θDl')为雷达圆形阵列的差波束期望主瓣
    BD0内第l'个采样点的有限近似采样角度对应的方向图,上标T为矩阵转置操作,A(θB)为雷
    达圆形阵列的方位角导向矢量,wD为待优化的雷达圆形阵列差波束权矢量,θB为雷达圆形阵
    列的方位角,θB=δ表示雷达圆形阵列的方位角与期望目标的方位角取值相等。
    6.如权利要求5所述的一种基于凸优化算法的雷达阵列和差波束方向图优化方法,其
    特征在于,在步骤5中,所述雷达圆形阵列和波束最优权矢量,其得到过程为:
    根据雷达圆形阵列的和波束方向图综合优化模型,计算得到雷达圆形阵列的和波束方
    向图的凸优化代价函数和约束条件为:
    <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mi>min</mi> </mtd> <mtd> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>F</mi> <mi>h</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>&theta;</mi> <mrow> <msup> <mi>Hm</mi> <mrow> <mo>&prime;</mo> <mo>&prime;</mo> </mrow> </msup> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>F</mi> <mi>H</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>&theta;</mi> <mrow> <msup> <mi>Hm</mi> <mo>&prime;</mo> </msup> </mrow> </msub> <mo>)</mo> </mrow> <mo>|</mo> <msub> <mo>|</mo> <mn>2</mn> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>&theta;</mi> <mrow> <msup> <mi>Hm</mi> <mo>&prime;</mo> </msup> </mrow> </msub> <mo>&Element;</mo> <msub> <mi>&Theta;</mi> <mrow> <mi>H</mi> <mi>m</mi> <mi>a</mi> <mi>i</mi> <mi>n</mi> </mrow> </msub> </mrow> </mtd> <mtd> <mrow> <msup> <mi>m</mi> <mo>&prime;</mo> </msup> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mn>...</mn> <mo>,</mo> <msup> <mi>M</mi> <mo>&prime;</mo> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>s</mi> <mo>.</mo> <mi>t</mi> <mo>.</mo> </mrow> </mtd> <mtd> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>F</mi> <mi>h</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>&theta;</mi> <mrow> <mi>H</mi> <mi>s</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>|</mo> <msub> <mo>|</mo> <mi>&infin;</mi> </msub> <mo>&le;</mo> <msub> <mi>&xi;</mi> <mi>H</mi> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>&theta;</mi> <mrow> <mi>H</mi> <mi>s</mi> </mrow> </msub> <mo>&Element;</mo> <msub> <mi>&Theta;</mi> <mrow> <mi>H</mi> <mi>s</mi> <mi>i</mi> <mi>d</mi> <mi>e</mi> </mrow> </msub> </mrow> </mtd> <mtd> <mrow> <mi>s</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mn>...</mn> <mo>,</mo> <mi>S</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow></mrow> </mtd> <mtd> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>F</mi> <mi>h</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>&theta;</mi> <mrow> <mi>H</mi> <mi>e</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>|</mo> <msub> <mo>|</mo> <mi>&infin;</mi> </msub> <mo>&le;</mo> <msub> <mi>&epsiv;</mi> <mi>H</mi> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>&theta;</mi> <mrow> <mi>H</mi> <mi>e</mi> </mrow> </msub> <mo>&Element;</mo> <msub> <mi>&Theta;</mi> <mrow> <mi>H</mi> <mi>n</mi> <mi>u</mi> <mi>l</mi> <mi>l</mi> </mrow> </msub> </mrow> </mtd> <mtd> <mrow> <mi>e</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mn>...</mn> <mo>,</mo> <mi>E</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow></mrow> </mtd> <mtd> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>F</mi> <mi>h</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>&theta;</mi> <mrow> <mi>H</mi> <mi>l</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>|</mo> <msub> <mo>|</mo> <mi>&infin;</mi> </msub> <mo>&le;</mo> <msub> <mi>&eta;</mi> <mi>H</mi> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>&theta;</mi> <mrow> <mi>H</mi> <mi>l</mi> </mrow> </msub> <mo>&Element;</mo> <msub> <mi>&Theta;</mi> <mrow> <mi>H</mi> <mi>l</mi> <mi>o</mi> <mi>w</mi> </mrow> </msub> </mrow> </mtd> <mtd> <mrow> <mi>l</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mn>...</mn> <mo>,</mo> <mi>L</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow></mrow> </mtd> <mtd> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>w</mi> <mi>H</mi> </msub> <mo>|</mo> <msub> <mo>|</mo> <mn>2</mn> </msub> <mo>&le;</mo> <mi>&upsi;</mi> </mrow> </mtd> <mtd> <mrow></mrow> </mtd> <mtd> <mrow></mrow> </mtd> </mtr> </mtable> </mfenced>
    其中,υ为设定的加权系数约束值,wH为待优化的雷达圆形阵列和波束权矢量;求得雷达
    圆形阵列的和波束方向图的凸优化代价函数和约束条件的最佳解,该最佳解即为雷达圆形
    阵列和波束最优权矢量w′H。
    7.如权利要求1所述的一种基于凸优化算法的雷达阵列和差波束方向图优化方法,其
    特征在于,在步骤6中,所述雷达圆形阵列的原始和波束接收能量和雷达圆形阵列的原始差
    波束接收能量,其得到过程为:
    将雷达圆形阵列和包含的M个阵元各自方向分别指向(δ-θ')方向和(δ+θ')方向,进而
    分别形成两个波束方向图,记为第一波束方向图F(δ-θ')和第二波束方向图F(δ+θ'),δ为期
    望目标的方位角,θ'表示第一波束方向图与第二波束方向图的夹角,且
    λ表示信号源向雷达圆形阵列发射信号的波长,R表示圆周的半径,M表示雷达圆形阵列包含
    的阵元个数,δ为期望目标的方位角,cos表示取余弦操作;
    然后将第一波束方向图F(δ-θ')和第二波束方向图F(δ+θ')叠加,得到雷达圆形阵列的
    原始和波束方向图;将第一波束方向图F(δ-θ')和第二波束方向图F(δ+θ')相减,得到雷达
    圆形阵列的原始差波束方向图;
    从而分别计算得到雷达圆形阵列的原始和波束接收能量EΣ和雷达圆形阵列的原始差
    波束接收能量E△,其表达式分别为:
    EΣ=FΣ2(δ)=[F(δ-θ)+F(δ+θ)]2
    E△=FΣ(δ)F△(δ)=[F(δ-θ)-F(δ+θ)][F(δ-θ)+F(δ+θ)]
    其中,FΣ(δ)为δ方向上雷达圆形阵列的原始和波束方向图,
    FΣ(δ)=F(δ-θ)+F(δ+θ);F△(δ)为δ方向上雷达圆形阵列的原始差波束方向图,
    F△(δ)=[F(δ-θ)-F(δ+θ)]。
    8.如权利要求1所述的一种基于凸优化算法的雷达阵列和差波束方向图优化方法,其
    特征在于,在步骤7中,所述期望目标的角误差归一化曲线的斜率,其得到过程为:
    根据公式:
    <mrow> <mi>&epsiv;</mi> <mo>=</mo> <mfrac> <mrow> <mo>|</mo> <msub> <mi>E</mi> <mi>&Delta;</mi> </msub> <mo>|</mo> </mrow> <mrow> <mo>|</mo> <msub> <mi>E</mi> <mi>&Sigma;</mi> </msub> <mo>|</mo> </mrow> </mfrac> <mo>=</mo> <mfrac> <mrow> <msub> <mi>F</mi> <mi>&Sigma;</mi> </msub> <mrow> <mo>(</mo> <mi>&delta;</mi> <mo>)</mo> </mrow> <msub> <mi>F</mi> <mi>&Delta;</mi> </msub> <mrow> <mo>(</mo> <mi>&delta;</mi> <mo>)</mo> </mrow> </mrow> <mrow> <msup> <msub> <mi>F</mi> <mi>&Sigma;</mi> </msub> <mn>2</mn> </msup> <mrow> <mo>(</mo> <mi>&delta;</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>=</mo> <mfrac> <mrow> <msub> <mi>F</mi> <mi>&Delta;</mi> </msub> <mrow> <mo>(</mo> <mi>&delta;</mi> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>F</mi> <mi>&Sigma;</mi> </msub> <mrow> <mo>(</mo> <mi>&delta;</mi> <mo>)</mo> </mrow> </mrow> </mfrac> </mrow>
    计算得到期望目标的角误差归一化曲线ε;并对期望目标的角误差归一化曲线ε进行拟
    合,计算得到期望目标的角误差归一化曲线的斜率
    9.如权利要求1所述的一种基于凸优化算法的雷达阵列和差波束方向图优化方法,其
    特征在于,在步骤8中,所述待优化的雷达圆形阵列差和差波束方向图的角误差曲线斜率,
    其得到过程为:
    计算得到待优化的雷达圆形阵列差和差波束方向图的角误差曲线ε′:
    <mrow> <msup> <mi>&epsiv;</mi> <mo>&prime;</mo> </msup> <mo>=</mo> <mfrac> <mrow> <mi>i</mi> <mi>m</mi> <mi>a</mi> <mi>g</mi> <mo>&lsqb;</mo> <msup> <mrow> <mo>(</mo> <msubsup> <mi>w</mi> <mi>H</mi> <mrow> <mo>&prime;</mo> <mi>T</mi> </mrow> </msubsup> <mi>A</mi> <mo>(</mo> <mi>&theta;</mi> <mo>)</mo> <mo>)</mo> </mrow> <mo>*</mo> </msup> <mrow> <mo>(</mo> <msubsup> <mi>w</mi> <mi>D</mi> <mi>T</mi> </msubsup> <mi>A</mi> <mo>(</mo> <mi>&theta;</mi> <mo>)</mo> </mrow> <mo>)</mo> <mo>&rsqb;</mo> </mrow> <mrow> <mo>|</mo> <msubsup> <mi>w</mi> <mi>H</mi> <mrow> <mo>&prime;</mo> <mi>T</mi> </mrow> </msubsup> <mi>A</mi> <mrow> <mo>(</mo> <mi>&theta;</mi> <mo>)</mo> </mrow> <msup> <mo>|</mo> <mn>2</mn> </msup> </mrow> </mfrac> </mrow>
    其中,w′H为雷达圆形阵列和波束最优权矢量,wD为待优化的雷达圆形阵列差波束权矢
    量,上标T为矩阵转置操作,上标*为矩阵取共轭操作,Imag表示取虚部操作,A(θ)为雷达圆
    形阵列的方位角导向矢量;
    然后对待优化的雷达圆形阵列差和差波束方向图的角误差曲线ε′进行拟合,进而得到
    待优化的雷达圆形阵列差和差波束方向图的角误差曲线斜率
    10.如权利要求5所述的一种基于凸优化算法的雷达阵列和差波束方向图优化方法,其
    特征在于,在步骤9中,所述雷达圆形阵列差波束的最优权矢量w′D,其得到过程为:
    首先计算雷达圆形阵列的差波束方向图优化模型:

    其中,ξD为雷达圆形阵列的差波束方向图预先设定的最大旁瓣电平,εD为雷达圆形阵列
    的差波束方向图预先设定的零陷方向最大电平,ηD为雷达圆形阵列的差波束方向图预先设
    定的低旁瓣区域最大预设电平,θDs'表示差波束旁瓣区域ΘDside对应的方位角区间内第s'个
    采样点的有限近似采样角度,θDs'∈ΘDside,s'∈{1,2,…,S'},S'表示差波束旁瓣区域
    ΘDside对应的方位角区间内内包含的采样点总个数;θDe'表示差波束零陷区域ΘDnull对应的
    方位角区间内第e'个采样点的有限近似采样角度,θDe'∈ΘDnull,e'∈{1,2,…,E'},E'表示
    差波束零陷区域ΘDnull对应的方位角区间内包含的采样点总个数;θDl'表示差波束低旁瓣区
    域ΘDlow对应的方位角区间内第l'个采样点的有限近似采样角度,θDl'∈ΘDlow,l'∈{1,
    2,…,L'},L'表示差波束低旁瓣区域ΘDlow对应的方位角区间内包含的采样点总个数;Fd
    (θDs')为雷达圆形阵列的差波束期望主瓣BD0内第s'个采样点的有限近似采样角度对应的方
    向图,Fd(θDe')为雷达圆形阵列的差波束期望主瓣BD0内第e'个采样点的有限近似采样角度
    对应的方向图,Fd(θDl')为雷达圆形阵列的差波束期望主瓣BD0内第l'个采样点的有限近似
    采样角度对应的方向图,||·||∞为向量无穷∞-范数操作,s.t.表示约束条件,A(θB)为雷
    达圆形阵列的方位角导向矢量,wD为待优化的雷达圆形阵列差波束权矢量,θB为雷达圆形阵
    列的方位角,θB=δ表示雷达圆形阵列的方位角与期望目标的方位角取值相等,上标T为矩
    阵转置操作,s.t.表示约束条件;表示待优化的雷达圆形阵列差和差波束方向图
    的角误差曲线斜率与期望目标的角误差归一化曲线的斜率取值相等;
    然后求得雷达圆形阵列的方向图优化模型的最佳解,所述雷达圆形阵列的方向图的优
    化模型的最佳解为雷达圆形阵列差波束的最优权矢量w′D。

    说明书

    基于凸优化算法的雷达阵列和差波束方向图优化方法

    技术领域

    本发明属于雷达技术领域,特别涉及一种基于凸优化算法的雷达阵列和差波束方
    向图优化方法,适用于雷达圆形阵列和波束和差波束的方向图优化。

    背景技术

    随着雷达技术的发展,相控阵雷达越来越多地应用到了实际的工程项目中。然而
    目前的相控阵雷达天线基本上还是线阵和传统平面阵列天线。线阵列的结构相对来说比较
    简单,易于处理,但实际中它只能有大约120度左右的方位角覆盖,其提供的增益和方向图
    等特性随扫描角的不同而改变,尤其是扫描角偏离阵列法线方向较大是,阵列性能急剧下
    降,所以这些大大限制了其使用范围。传统的平面阵列同样也存在着一些类似缺点,如波束
    扫描范围窄(局限在120度以内);波束宽度随着扫描角的增加而增加;阵列单元之间的互耦
    效应是扫描角的函数,难于保持平衡等,圆形阵列能避免上述阵列的一些缺点,它不仅可以
    提供360度的方位角,而且通过循环移动阵列激励,简单而灵活的操纵波束的方位,而在俯
    仰方向上也有一个理想的方向特性。同时,圆结构的对称性使其基本维持等性能的波束形
    状和天线增益,并大体上保持互耦平衡。

    得益于计算机技术的飞速发展,高效的阵列方向图综合方法已成为研究热点。目
    前对于线阵及平面阵的波束方向图的综合主要方法是对和差波束进行加窗处理,来达到控
    制副瓣电平的目的。而对于圆形阵列的波束方向图优化主要应用的综合算法有遗传算法、
    模拟退火算法、粒子群算法,以及蚁群算法,这些算法从本质上来说都是基于随机性的自然
    算法,存在综合过程中需要进行大规模的搜索,导致此类算法的收敛速度慢,并容易陷入局
    部最优解等问题。

    发明内容

    针对上述现有技术存在的不足,本发明的目的在于提出一种基于凸优化算法的雷
    达阵列和差波束方向图优化方法,该种基于凸优化算法的雷达阵列和差波束方向图优化方
    法最突出的优势在于,利用凸优化算法能够稳定、有效地求取方向图的最优权值,易于工程
    实现。

    实现本发明目的的技术思路是:将雷达圆形阵列和差波束方向图优化为一个凸优
    化问题,并建立合适的目标函数和约束条件,然后利用MATLAB的凸优化工具包,求得雷达圆
    形阵列和差波束方向图的最优权值。

    为达到上述技术目的,本发明采用如下技术方案予以实现。

    一种基于凸优化算法的雷达阵列和差波束方向图优化方法,包括以下步骤:

    步骤1,确定包含M个阵元的雷达圆形阵列均匀分布在半径为R的圆周上,然后计算
    得到雷达圆形阵列的方向矢量;M、R分别为大于0的正整数;

    步骤2,根据雷达圆形阵列的方向矢量,分别计算得到雷达圆形阵列的和波束方向
    图和雷达圆形阵列的差波束方向图,进而分别得到和波束主瓣区域、和波束零陷区域、和波
    束旁瓣区域、和波束低旁瓣区域,以及差波束主瓣区域、差波束零陷区域、差波束旁瓣区域、
    差波束低旁瓣区域;

    步骤3,分别得到约束和波束主瓣区域的代价函数和约束差波束主瓣区域的代价
    函数;

    步骤4,根据约束和波束主瓣区域的代价函数和约束差波束主瓣区域的代价函数,
    分别得到雷达圆形阵列的和波束方向图综合优化模型和雷达圆形阵列的差波束方向图综
    合优化模型;

    步骤5,根据雷达圆形阵列的和波束方向图综合优化模型,计算得到雷达圆形阵列
    和波束最优权矢量;

    步骤6,将雷达圆形阵列包含的M个阵元各自方向分别指向第一设定方向和第二设
    定方向,所述第一设定方向和第二设定方向为不同的方向;进而分别形成两个波束方向图,
    记为第一波束方向图和第二波束方向图,然后分别计算得到雷达圆形阵列的原始和波束接
    收能量和雷达圆形阵列的原始差波束接收能量;

    步骤7,根据雷达圆形阵列的原始和波束接收能量和雷达圆形阵列的原始差波束
    接收能量,计算得到期望目标的角误差归一化曲线,并拟合该期望目标的角误差归一化曲
    线得到期望目标的角误差归一化曲线的斜率;

    步骤8,根据雷达圆形阵列和波束最优权矢量,计算得到待优化的雷达圆形阵列差
    和差波束方向图的角误差曲线,并拟合得到待优化的雷达圆形阵列差和差波束方向图的角
    误差曲线斜率;

    步骤9,根据雷达圆形阵列的差波束方向图综合优化模型和待优化的雷达圆形阵
    列差和差波束方向图的角误差曲线斜率,计算得到雷达圆形阵列差波束的最优权矢量。

    本发明与现有技术相比具有以下优点:

    现有的雷达圆形阵列的波束方向图优化主要应用的综合算法包含遗传算法、模拟
    退火算法、粒子群算法,以及蚁群算法,这些算法均在综合过程中分别需要进行大规模的搜
    索,导致这些算法的收敛速度都很慢,并都容易陷入局部最优解的问题;而本发明方法利用
    凸优化算法对雷达圆形阵列的波束方向图进行优化,能够稳定、有效、快速地求解出方向图
    优化的最优权值,且本发明方法步骤简单,易于工程实现;另外,本发明方法能有效降低雷
    达圆形阵列和差波束的方向图的副瓣,且能够抑制干扰,并能够增强雷达在有干扰环境下
    的杂波抑制性能。

    附图说明

    下面结合附图和具体实施方式对本发明作进一步详细说明。

    图1为本发明的一种基于凸优化算法的雷达阵列和差波束方向图优化方法实现流
    程图;

    图2为雷达圆形阵列所在场景的几何示意图;

    图3为使用本发明方法得到的雷达阵列和差波束方向图:

    图4为雷达圆形阵列的原始归一化误差曲线对比曲线图:

    图5为使用本发明方法后得到的雷达圆形阵列归一化误差曲图。

    具体实施方式

    参照图1,为本发明的一种基于凸优化算法的雷达阵列和差波束方向图优化方法
    实现流程图;其中所述基于凸优化算法的雷达阵列和差波束方向图优化方法,包括以下步
    骤:

    步骤1,确定包含M个阵元的雷达圆形阵列均匀分布在半径为R的圆周上,然后计算
    得到雷达圆形阵列的方向矢量;M、R分别为大于0的正整数。

    具体地,确定包含M个阵元的雷达圆形阵列按照阵元序号依次均匀分布在半径为R
    的圆周上,其中M个阵元分别为无方向性;并且距离雷达圆形阵列的设定位置处存在信号
    源,其中设定位置为到雷达圆形阵列的距离大于圆周半径R,所述信号源向雷达圆形阵列发
    射信号,分别将信号入射进雷达圆形阵列的方位角记为信号源的方位角θ,将信号入射进雷
    达圆形阵列的俯仰角记为信号源的俯仰角φ;同时确定雷达圆形阵列检测范围内存在期望
    目标,该期望目标的方位角为δ,期望目标的俯仰角为β。

    如图2所示,建立雷达圆形阵列所在场景的几何示意图,即以雷达圆形阵列的中心
    为原点o建立三维坐标系xoyz,在三维坐标系xoyz中,首先确定信号源发射信号的入射方向
    和期望目标方向,所述期望目标方向由期望目标的方位角和期望目标的俯仰角确定,其中
    期望目标的方位角是期望目标方向与z轴的夹角,所述期望目标的俯仰角是期望目标方向
    在xoy平面的投影与x轴的夹角;进而计算得到雷达圆形阵列的方向矢量a(θ,φ),其表达式
    为:


    其中,θ表示信号源的方位角,φ表示信号源的俯仰角;令γm表示第m个阵元到第1
    个阵元的角度,γm=2πm/M,m∈{0,1,…,M-1},λ表示信号源向雷达圆形阵列发射信号的波
    长,R表示圆周的半径,exp表示指数函数操作,j为虚数单位,sin表示求正弦操作,cos表示
    求余弦操作,M表示雷达圆形阵列包含的阵元个数,M、R分别为大于0的正整数。

    步骤2,根据雷达圆形阵列的方向矢量a(θ,φ),分别计算得到雷达圆形阵列的和
    波束方向图FH(θ)和雷达圆形阵列的差波束方向图FD(θ),进而分别得到和波束主瓣区域
    ΘHmain、和波束零陷区域ΘHnull、和波束旁瓣区域ΘHside、和波束低旁瓣区域ΘHlow,以及差波
    束主瓣区域ΘDmain、差波束零陷区域ΘDnull、差波束旁瓣区域ΘDside、差波束低旁瓣区域
    ΘDlow。

    步骤2的子步骤为:

    2a)设定待优化的雷达圆形阵列和波束权矢量为wH,wH=[wH1,…,wHm,…,wHM]T,上
    标T表示转置,m∈{1,2,…,M},wHm表示待优化的雷达圆形阵列和波束权矢量内第m个阵元的
    和波束权值,M表示雷达圆形阵列包含的阵元个数,与待优化的雷达圆形阵列和波束权矢量
    内包含的阵元个数相同;并根据待优化的雷达圆形阵列和波束权矢量wH和雷达圆形阵列的
    方向矢量a(θ,φ),计算得到雷达圆形阵列的和波束方向图FH(θ),其表达式为:

    FH(θ)=wH×a(θ)

    其中,a(θ)表示雷达圆形阵列的方向矢量a(θ,φ)的方位维矢量。

    2b)设定待优化的雷达圆形阵列差波束权矢量为wD,wD=[wD1,…,wDm,…,wDM]T,上
    标T表示转置,wDm表示待优化的雷达圆形阵列差波束权矢量内第m个阵元的差波束权值,m∈
    {1,2,…,M},M表示雷达圆形阵列包含的阵元个数,与待优化的雷达圆形阵列差波束权矢量
    内包含的阵元个数相同;利用待优化的雷达圆形阵列差波束权矢量wD和雷达圆形阵列的方
    向矢量a(θ,φ),计算得到雷达圆形阵列的差波束方向图FD(θ),其表达式为:

    FD(θ)=wD×a(θ)

    其中,a(θ)表示雷达圆形阵列的方向矢量a(θ,φ)的方位维矢量。

    2c)为了得到更好的波束优化结果,将雷达圆形阵列的和波束方向图FH(θ)和雷达
    圆形阵列的差波束方向图FD(θ)分别进行区域划分,将雷达圆形阵列的和波束方向图FH(θ)
    划分为和波束主瓣区域ΘHmain、和波束零陷区域ΘHnull、和波束旁瓣区域ΘHside、和波束低旁
    瓣区域ΘHlow;所述和波束主瓣区域ΘHmain为雷达圆形阵列的和波束方向图FH(θ)内接收功
    率在-13dB到0dB之间的区域,所述和波束零陷区域ΘHnull为雷达圆形阵列的和波束方向图
    FH(θ)内接收功率小于-80dB的区域,所述和波束旁瓣区域ΘHside为雷达圆形阵列的和波束
    方向图FH(θ)内接收功率小于-30dB的区域,所述和波束低旁瓣区域ΘHlow为雷达圆形阵列的
    和波束方向图FH(θ)内接收功率小于-40dB的区域;其中和波束主瓣区域ΘHmain内包含待优
    化的和波束的主瓣BH。

    将雷达圆形阵列的差波束方向图FD(θ)划分为:差波束主瓣区域ΘDmain、差波束零
    陷区域ΘDnull、差波束旁瓣区域ΘDside、差波束低旁瓣区域ΘDlow,所述差波束主瓣区域
    ΘDmain为雷达圆形阵列的差波束方向图FD(θ)内接收功率在-13dB到0dB之间的区域,所述差
    波束零陷区域ΘDnull为雷达圆形阵列的差波束方向图FD(θ)内接收功率小于-80dB的区域,
    所述差波束低旁瓣区域ΘDlow为雷达圆形阵列的差波束方向图FD(θ)内接收功率小于-30dB
    的区域,所述差波束低旁瓣区域ΘDlow为雷达圆形阵列的差波束方向图FD(θ)内接收功率小
    于-40dB的区域;其中差波束主瓣区域ΘDmain内包含待优化的差波束主瓣BD。

    步骤3,分别得到约束和波束主瓣区域ΘHmain的代价函数和约束差波束主瓣区域
    ΘDmain的代价函数。

    3a)范数逼近是一个可解的凸问题,利用范数逼近的方法将雷达圆形阵列的和波
    束和雷达圆形阵列的差波束分别优化转化为凸优化问题;设置雷达圆形阵列的和波束期望
    主瓣BH0,并在和主瓣区域ΘHmain内将待优化的和波束的主瓣BH和雷达圆形阵列的和波束期
    望主瓣BH0之差的向量2-范数最小化,进而得到约束和波束主瓣区域ΘHmain的代价函数:

    min||Fh(θHm”)-FH(θHm')||2θHm'∈ΘHmain m'=1,2,…,M'

    其中,m'∈{1,2,…,M'},M'表示和主瓣区域ΘHmain对应的方位角区间内包含的采
    样点总个数,θHm'表示和主瓣区域ΘHmain对应的方位角区间内第m'个采样点的有限近似采样
    角度,θHm'∈ΘHmain,Fh(θHm”)为雷达圆形阵列的和波束期望主瓣BH0内第m”个采样点有限近似
    采样角度对应的方向图,FH(θHm')为和主瓣区域ΘHmain对应的方位角区间内第m'个采样点的
    有限近似采样角度对应方向图,∈表示属于,||·||2为向量2-范数,min为求最小值操作,
    m”∈{1,2,…,M”},M”为雷达圆形阵列的和波束期望主瓣BH0包含的采样点总个数;雷达圆形
    阵列的和波束期望主瓣BH0包含的采样点总个数与主瓣区域Θmain对应的方位角区间内包含
    的采样点总个数取值相等,∈表示属于。

    3b)确定雷达圆形阵列的差波束方向图FD(θ)的零点,在该雷达圆形阵列的差波束
    方向图FD(θ)的零点方向上雷达圆形阵列的差波束接收功率接近0。

    首先确定雷达圆形阵列的差波束方向图上的半主瓣宽度R为圆
    阵半径,λ为信号源向雷达圆形阵列发射信号的波长;然后将雷达圆形阵列的差波束方向图
    中从至θ区域,以及从θ至区域分别设置为雷达圆形阵列差波束方向图的
    主瓣区域,θ为信号源的方位角,δ为期望目标的方位角。

    然后设置雷达圆形阵列的差波束期望主瓣BD0,并在差波束主瓣区域ΘDmain内将待
    优化的差波束主瓣BD与雷达圆形阵列的差波束期望主瓣BD0之差的向量2-范数最小化,得到
    约束差波束主瓣区域ΘDmain的代价函数:

    min||Fd(θDn')-FD(θDn)||2θDn∈Θmain n=1,2,…,N

    其中,θDn∈ΘDmain,n∈{1,2,…,N},θDn表示差主瓣区域ΘDmain对应的方位角区间内
    第n个采样点的有限近似采样角度,N为差主瓣区域ΘDmain对应的方位角区间内包含的采样
    点总个数,FD(θDn)为差主瓣区域ΘDmain对应的方位角区间内第n个采样点的有限近似采样角
    度对应的方向图,Fd(θDn')为雷达圆形阵列的差波束期望主瓣BD0内第n'个采样点有限近似
    采样角度对应的方向图,||·||2为向量2-范数,∈表示属于,min为求最小值操作,n'∈{1,
    2,…,N'},N'为雷达圆形阵列的差波束期望主瓣BD0内包含的采样点总个数;雷达圆形阵列
    的差波束期望主瓣BD0内包含的采样点总个数与差主瓣区域ΘDmain对应的方位角区间内包
    含的采样点总个数取值相等。

    步骤4,根据约束和波束主瓣区域ΘHmain的代价函数和约束差波束主瓣区域ΘDmain
    的代价函数,分别得到雷达圆形阵列的和波束方向图综合优化模型和雷达圆形阵列的差波
    束方向图综合优化模型。

    4a)根据约束和波束主瓣区域ΘHmain的代价函数,控制和波束零陷区域ΘHnull、和
    波束旁瓣区域ΘHside、和波束低旁瓣区域ΘHlow分别在设定电平之下,对雷达圆形阵列的和
    波束期望主瓣BH0进行∞-范数约束,进而得到雷达圆形阵列的和波束方向图综合优化模型:


    其中,m'∈{1,2,…,M'},M'表示和主瓣区域ΘHmain对应的方位角区间内包含的采
    样点总个数,θHm'表示和主瓣区域ΘHmain对应的方位角区间内第m'个采样点的有限近似采样
    角度,Fh(θHm”)为雷达圆形阵列的和波束期望主瓣BH0内第m”个采样点有限近似采样角度对
    应的方向图,FH(θHm')为和主瓣区域ΘHmain对应的方位角区间内第m'个采样点的有限近似采
    样角度对应方向图,∈表示属于,m”∈{1,2,…,M”},M”为雷达圆形阵列的和波束期望主瓣
    BH0包含的采样点总个数,ξH为雷达圆形阵列的和波束方向图预先设定的最大旁瓣电平,通
    常取值为不大于-30dB;εH为雷达圆形阵列的和波束方向图预先设定的零陷方向最大电平,
    通常取值为不大于-80dB;ηH为雷达圆形阵列的和波束方向图预先设定的低旁瓣区域最大
    预设电平,通常取值为不大于-40dB;θHs表示和波束旁瓣区域ΘHside对应的方位角区间内第
    s个采样点的有限近似采样角度,θHs∈ΘHside,s∈{1,2,…,S},S表示和波束旁瓣区域ΘHside
    对应的方位角区间内包含的采样点总个数;θHe表示和波束零陷区域ΘHnull对应的方位角区
    间内第e个采样点的有限近似采样角度,θHe∈ΘHnull,e∈{1,2,…,E},E表示和波束零陷区
    域ΘHnull对应的方位角区间内包含的采样点总个数;θHl表示和波束低旁瓣区域ΘHlow对应的
    方位角区间内第l个采样点的有限近似采样角度,θHl∈ΘHlow,l∈{1,2,…,L},L表示和波束
    低旁瓣区域ΘHlow对应的方位角区间内包含的采样点总个数;Fh(θHs)为雷达圆形阵列的和
    波束期望主瓣BH0内第s个采样点的有限近似采样角度对应的方向图,Fh(θHe)为雷达圆形阵
    列的和波束期望主瓣BH0内第e个采样点的有限近似采样角度对应的方向图,Fh(θHl)为雷达
    圆形阵列的和波束期望主瓣BH0内第l个采样点的有限近似采样角度对应的方向图,||·||∞
    为向量无穷∞-范数操作,s.t.表示约束条件,∈表示属于。

    4b)根据约束差波束主瓣区域ΘDmain的代价函数,控制差波束零陷区域ΘDnull、差
    波束旁瓣区域ΘDside、差波束低旁瓣区域ΘDlow分别在某个设计的电平之下,对雷达圆形阵
    列的差波束期望主瓣BD0进行∞-范数约束,得到雷达圆形阵列的差波束方向图综合优化模
    型:


    其中,n∈{1,2,…,N},θDn表示差主瓣区域ΘDmain对应的方位角区间内第n个采样
    点的有限近似采样角度,N为差主瓣区域ΘDmain对应的方位角区间内包含的采样点总个数,
    FD(θDn)为差主瓣区域ΘDmain对应的方位角区间内第n个采样点的有限近似采样角度对应的
    方向图,Fd(θDn')为雷达圆形阵列的差波束期望主瓣BD0内第n'个采样点有限近似采样角度
    对应的方向图,||·||2为向量2-范数,min为求最小值操作,n'∈{1,2,…,N'},N'为雷达圆
    形阵列的差波束期望主瓣BD0内包含的采样点总个数;ξD为雷达圆形阵列的差波束方向图预
    先设定的最大旁瓣电平,通常取值为不大于-30dB;εD为雷达圆形阵列的差波束方向图预先
    设定的零陷方向最大电平,通常取值为不大于-80dB;ηD为雷达圆形阵列的差波束方向图预
    先设定的低旁瓣区域最大预设电平,通常取值为不大于-40dB;θDs'表示差波束旁瓣区域
    ΘDside对应的方位角区间内第s'个采样点的有限近似采样角度,θDs'∈ΘDside,s'∈{1,
    2,…,S'},S'表示差波束旁瓣区域ΘDside对应的方位角区间内内包含的采样点总个数;θDe'
    表示差波束零陷区域ΘDnull对应的方位角区间内第e'个采样点的有限近似采样角度,θDe'∈
    ΘDnull,e'∈{1,2,…,E'},E'表示差波束零陷区域ΘDnull对应的方位角区间内包含的采样
    点总个数;θDl'表示差波束低旁瓣区域ΘDlow对应的方位角区间内第l'个采样点的有限近似
    采样角度,θDl'∈ΘDlow,l'∈{1,2,…,L'},L'表示差波束低旁瓣区域ΘDlow对应的方位角区
    间内包含的采样点总个数;Fd(θDs')为雷达圆形阵列的差波束期望主瓣BD0内第s'个采样点
    的有限近似采样角度对应的方向图,Fd(θDe')为雷达圆形阵列的差波束期望主瓣BD0内第e'
    个采样点的有限近似采样角度对应的方向图,Fd(θDl')为雷达圆形阵列的差波束期望主瓣
    BD0内第l'个采样点的有限近似采样角度对应的方向图,上标T为矩阵转置操作,||·||∞为
    向量无穷∞-范数操作,s.t.表示约束条件,∈表示属于,A(θB)为雷达圆形阵列的方位角导
    向矢量,wD为待优化的雷达圆形阵列差波束权矢量,θB为雷达圆形阵列的方位角,θB=δ表示
    雷达圆形阵列的方位角与期望目标的方位角取值相等。

    步骤5,根据雷达圆形阵列的和波束方向图综合优化模型,计算得到雷达圆形阵列
    和波束最优权矢量w′H。

    具体地,为了提高雷达圆形阵列的辐射效率,在方向设计算法的最优化过程中增
    加雷达圆形阵列的功率约束,


    wHm表示待优化的雷达圆形阵列和波束权矢量内第m个阵元的和波束权值,M表示雷
    达圆形阵列包含的阵元个数,υ为设定的加权系数约束值,设定的加权系数约束值取值范围
    为0到1之间,通常取值为0.5;用以控制权矢量的稳健性和增益损失;根据雷达圆形阵列的
    和波束方向图综合优化模型,计算得到雷达圆形阵列的和波束方向图的凸优化代价函数和
    约束条件为:


    其中,ξH为雷达圆形阵列的和波束方向图预先设定的最大旁瓣电平,通常取值为
    不大于-30dB;εH为雷达圆形阵列的和波束方向图预先设定的零陷方向最大电平,通常取值
    为不大于-80dB;ηH为雷达圆形阵列的和波束方向图预先设定的低旁瓣区域最大预设电平,
    通常取值为不大于-40dB;θHs表示和波束旁瓣区域ΘHside对应的方位角区间内第s个采样点
    的有限近似采样角度,θHs∈ΘHside,s∈{1,2,…,S},S表示和波束旁瓣区域ΘHside对应的方
    位角区间内包含的采样点总个数;θHe表示和波束零陷区域ΘHnull对应的方位角区间内第e
    个采样点的有限近似采样角度,θHe∈ΘHnull,e∈{1,2,…,E},E表示和波束零陷区域ΘHnull
    对应的方位角区间内包含的采样点总个数;θHl表示和波束低旁瓣区域ΘHlow对应的方位角
    区间内第l个采样点的有限近似采样角度,θHl∈ΘHlow,l∈{1,2,…,L},L表示和波束低旁瓣
    区域ΘHlow对应的方位角区间内包含的采样点总个数;Fh(θHs)为雷达圆形阵列的和波束期
    望主瓣BH0内第s个采样点的有限近似采样角度对应的方向图,Fh(θHe)为雷达圆形阵列的和
    波束期望主瓣BH0内第e个采样点的有限近似采样角度对应的方向图,Fh(θHl)为雷达圆形阵
    列的和波束期望主瓣BH0内第l个采样点的有限近似采样角度对应的方向图,||·||∞为向量
    无穷∞-范数操作,s.t.表示约束条件,υ为设定的加权系数约束值,wH为待优化的雷达圆形
    阵列和波束权矢量。

    然后利用MATLAB的CVX(凸优化)工具包,求得雷达圆形阵列的和波束方向图的凸
    优化代价函数和约束条件的最佳解,该最佳解即为雷达圆形阵列和波束最优权矢量w′H。

    步骤6,将雷达圆形阵列和包含的M个阵元各自方向分别指向第一设定方向和第二
    设定方向,所述第一设定方向和第二设定方向为不同的方向;进而分别形成两个波束方向
    图,记为第一波束方向图和第二波束方向图,然后分别计算得到雷达圆形阵列的原始和波
    束接收能量EΣ和雷达圆形阵列的原始差波束接收能量E△。

    具体地,利用双波束指向法,将雷达圆形阵列和包含的M个阵元各自方向分别指向
    (δ-θ')方向和(δ+θ')方向,进而分别形成两个波束方向图,记为第一波束方向图F(δ-θ')和
    第二波束方向图F(δ+θ'),δ为期望目标的方位角,θ'表示第一波束方向图与第二波束方向
    图的夹角,且λ表示信号源向雷达圆形阵列发射信号的波长,R表示圆周
    的半径,M表示雷达圆形阵列包含的阵元个数,δ为期望目标的方位角,cos表示取余弦操作。

    然后将第一波束方向图F(δ-θ')和第二波束方向图F(δ+θ')叠加,得到雷达圆形阵
    列的原始和波束方向图;将第一波束方向图F(δ-θ')和第二波束方向图F(δ+θ')相减,得到
    雷达圆形阵列的原始差波束方向图。

    从而分别计算得到雷达圆形阵列的原始和波束接收能量EΣ和雷达圆形阵列的原
    始差波束接收能量E△,其表达式分别为:

    EΣ=FΣ2(δ)=[F(δ-θ)+F(δ+θ)]2

    E△=FΣ(δ)F△(δ)=[F(δ-θ)-F(δ+θ)][F(δ-θ)+F(δ+θ)]

    其中,FΣ(δ)为δ方向上雷达圆形阵列的原始和波束方向图,

    FΣ(δ)=F(δ-θ)+F(δ+θ);F△(δ)为δ方向上雷达圆形阵列的原始差波束方向图,

    F△(δ)=[F(δ-θ)-F(δ+θ)]。

    步骤7,根据雷达圆形阵列的原始和波束接收能量EΣ和雷达圆形阵列的原始差波
    束接收能量E△,计算得到期望目标的角误差归一化曲线,并拟合该期望目标的角误差归一
    化曲线得到期望目标的角误差归一化曲线的斜率

    具体地,根据公式:


    计算得到期望目标的角误差归一化曲线ε;由于期望目标的角误差归一化曲线ε近
    似线性关系,在工程上近似为一条直线表示,并利用MATLAB中的polyfit()函数,对期望目
    标的角误差归一化曲线ε进行拟合,计算得到期望目标的角误差归一化曲线的斜率

    步骤8,根据雷达圆形阵列和波束最优权矢量w′H,计算得到待优化的雷达圆形阵
    列差和差波束方向图的角误差曲线,并拟合得到待优化的雷达圆形阵列差和差波束方向图
    的角误差曲线斜率。

    具体地,运用公式,计算得到待优化的雷达圆形阵列差和差波束方向图的角误差
    曲线ε′:


    其中,w′H为雷达圆形阵列和波束最优权矢量,wD为待优化的雷达圆形阵列差波束
    权矢量,上标T为矩阵转置操作,上标*为矩阵取共轭操作,Imag表示取虚部操作,A(θ)为雷
    达圆形阵列的方位角导向矢量。

    由于待优化的雷达圆形阵列差和差波束方向图的角误差曲线ε′近似线性关系,在
    工程上近似用一条直线表示,并利用MATLAB中的polyfit()函数,对待优化的雷达圆形阵
    列差和差波束方向图的角误差曲线ε′进行拟合,进而得到待优化的雷达圆形阵列差和差波
    束方向图的角误差曲线斜率

    步骤9,根据雷达圆形阵列的差波束方向图综合优化模型和待优化的雷达圆形阵
    列差和差波束方向图的角误差曲线斜率计算得到雷达圆形阵列差波束的最优权矢量
    w′D。

    具体地,对优化的和差波束的角误差归一化曲线进行约束,使得它和原始的和差
    波束产生的角误差归一化曲线相一致,使得待优化的雷达圆形阵列差和差波束方向图的角
    误差曲线斜率与期望目标的角误差归一化曲线的斜率取值相等,即进而得
    到雷达圆形阵列的差波束方向图优化模型:


    其中,ξD为雷达圆形阵列的差波束方向图预先设定的最大旁瓣电平,通常取值为
    不大于-30dB;εD为雷达圆形阵列的差波束方向图预先设定的零陷方向最大电平,通常取值
    为不大于-80dB;ηD为雷达圆形阵列的差波束方向图预先设定的低旁瓣区域最大预设电平,
    通常取值为不大于-40dB;θDs'表示差波束旁瓣区域ΘDside对应的方位角区间内第s'个采样
    点的有限近似采样角度,θDs'∈ΘDside,s'∈{1,2,…,S'},S'表示差波束旁瓣区域ΘDside对
    应的方位角区间内内包含的采样点总个数;θDe'表示差波束零陷区域ΘDnull对应的方位角区
    间内第e'个采样点的有限近似采样角度,θDe'∈ΘDnull,e'∈{1,2,…,E'},E'表示差波束零
    陷区域ΘDnull对应的方位角区间内包含的采样点总个数;θDl'表示差波束低旁瓣区域ΘDlow
    对应的方位角区间内第l'个采样点的有限近似采样角度,θDl'∈ΘDlow,l'∈{1,2,…,L'},
    L'表示差波束低旁瓣区域ΘDlow对应的方位角区间内包含的采样点总个数;Fd(θDs')为雷达
    圆形阵列的差波束期望主瓣BD0内第s'个采样点的有限近似采样角度对应的方向图,Fd
    (θDe')为雷达圆形阵列的差波束期望主瓣BD0内第e'个采样点的有限近似采样角度对应的方
    向图,Fd(θDl')为雷达圆形阵列的差波束期望主瓣BD0内第l'个采样点的有限近似采样角度
    对应的方向图,||·||∞为向量无穷∞-范数操作,s.t.表示约束条件,A(θB)为雷达圆形阵
    列的方位角导向矢量,wD为待优化的雷达圆形阵列差波束权矢量,θB为雷达圆形阵列的方位
    角,θB=δ表示雷达圆形阵列的方位角与期望目标的方位角取值相等,上标T为矩阵转置操
    作,s.t.表示约束条件。

    然后利用MATLAB的CVX(凸优化)工具包,求得雷达圆形阵列的方向图优化模型的
    最佳解,所述雷达圆形阵列的方向图的优化模型的最佳解为雷达圆形阵列差波束的最优权
    矢量w′D。

    至此求得雷达圆形阵列和波束最优权矢量w′H和雷达圆形阵列差波束的最优权矢
    量w′D,完成对雷达圆形阵列的和差波束方向图优化,所述雷达圆形阵列的和差波束方向图
    为雷达圆形阵列的和波束方向图FH(θ)和雷达圆形阵列的差波束方向图FD(θ)。

    通过以下仿真试验对本发明效果作进一步验证说明。

    (一)仿真条件:

    雷达圆形阵列为均匀圆阵,该雷达圆形阵列包含34个阵元,阵元间隔角度为5°,波
    长半径比2.6;设置和波束低旁瓣区域ΘHlow和差波束低旁瓣区域ΘDlow都为(35°,40°),设置
    和波束零陷区域ΘHnull和差波束零陷区域ΘDnull都为40±2°,设置雷达圆形阵列的和波束
    方向图预先设定的零陷方向最大电平和雷达圆形阵列的差波束方向图预先设定的零陷方
    向最大电平都为-80dB,设置雷达圆形阵列的和波束方向图预先设定的最大旁瓣电平和雷
    达圆形阵列的差波束方向图预先设定的最大旁瓣电平都为-30dB,设置雷达圆形阵列的和
    波束方向图预先设定的低旁瓣区域最大预设电平和雷达圆形阵列的差波束方向图预先设
    定的低旁瓣区域最大预设电平都为-40dB,设定的加权系数约束值为0.5。

    (二)仿真内容:采用本发明优化的和差波束方向图,仿真结果如图3所示;由原始
    和差波束方向图得到的角误差曲线,仿真结果如图4所示;由优化后和差波束方向图得到的
    角误差曲线,仿真结果如图5所示。

    (三)仿真分析

    从图3可以看出,采用本发明优化得到圆形阵列和差方向图符合预期设置的要求。

    从图4和图5可以看出,优化后的和差方向图角误差曲线和原始和差方向图角误差
    曲线相一致,对该阵列测角性能没有造成损失。

    综上所述,仿真实验验证了本发明的正确性,有效性和可靠性。

    显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精
    神和范围;这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围
    之内,则本发明也意图包含这些改动和变型在内。

    关于本文
    本文标题:基于凸优化算法的雷达阵列和差波束方向图优化方法.pdf
    链接地址://www.4mum.com.cn/p-6001483.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

    [email protected] 2017-2018 www.4mum.com.cn网站版权所有
    经营许可证编号:粤ICP备17046363号-1 
     


    收起
    展开
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
  • 北京pk10苹果计划软件 全天北京pk赛车计划人工 开奖公告快乐时时10分钟赚10万 幸运飞艇必中计划软件 快3大小单双玩法 双色球胆拖预测 龙虎游戏500元倍投方案稳赚 北京pk拾全天连中计划 定位胆规律 江西时时提前知道 山东时时开奖结果 pk10两期计划群官方 极速pk10官网下载 彩票开奖查询七星彩 网上购彩何时恢复 快三大小单双怎么玩稳赚