CN111984911A - 基于电磁对消的阵列天线辐射与散射综合优化方法 - Google Patents

基于电磁对消的阵列天线辐射与散射综合优化方法 Download PDF

Info

Publication number
CN111984911A
CN111984911A CN202010825165.8A CN202010825165A CN111984911A CN 111984911 A CN111984911 A CN 111984911A CN 202010825165 A CN202010825165 A CN 202010825165A CN 111984911 A CN111984911 A CN 111984911A
Authority
CN
China
Prior art keywords
array antenna
radiation
antenna
field
array
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.)
Granted
Application number
CN202010825165.8A
Other languages
English (en)
Other versions
CN111984911B (zh
Inventor
张帅
黄楠
白婵
张倩铢
刘雨果
柏文泉
周肖
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xidian University
Original Assignee
Xidian University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xidian University filed Critical Xidian University
Priority to CN202010825165.8A priority Critical patent/CN111984911B/zh
Publication of CN111984911A publication Critical patent/CN111984911A/zh
Application granted granted Critical
Publication of CN111984911B publication Critical patent/CN111984911B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/12—Computing arrangements based on biological models using genetic models
    • G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
    • H—ELECTRICITY
    • H01—ELECTRIC ELEMENTS
    • H01Q—ANTENNAS, i.e. RADIO AERIALS
    • H01Q21/00—Antenna arrays or systems
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Data Mining & Analysis (AREA)
  • Biophysics (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Artificial Intelligence (AREA)
  • Geometry (AREA)
  • Mathematical Analysis (AREA)
  • Computer Hardware Design (AREA)
  • Algebra (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Computational Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Medical Informatics (AREA)
  • Physiology (AREA)
  • Genetics & Genomics (AREA)
  • Biomedical Technology (AREA)
  • Computational Linguistics (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • Variable-Direction Aerials And Aerial Arrays (AREA)

Abstract

本发明提出了一种基于电磁对消的阵列天线辐射与散射综合优化方法，旨在实现良好辐射特性的同时，有效降低阵列天线的RCS，实现步骤为：构建电磁对消场景模型；获取阵列天线的辐射场；获取阵列天线的散射场；计算阵列天线总场；定义优化算法适应度函数fit(x)；利用优化算法获取阵列天线辐射和散射特性的综合优化结果。本发明解决了阵列天线辐射和散射特性综合优化的问题，实现了良好辐射特性的同时，有效降低阵列天线的RCS。

Description

基于电磁对消的阵列天线辐射与散射综合优化方法
技术领域
本发明属于天线技术领域，涉及一种阵列天线辐射与散射综合优化方法，具体涉及一种基于电磁对消的阵列天线辐射与散射综合优化方法，可用于阵列天线综合优化问题。
背景技术
随着探测技术的发展，目标雷达截面的减缩无疑具有重大的意义。天线作为发射和接收电磁波的散射体，随着其在军事平台的广泛应用，降低其RCS即有效降低设备被探测到的几率，因而研究阵列天线的RCS是一个关键课题。到目前为止，国内外仍然没有找到一种完全理想的途径：既可显著降低阵列天线的RCS又完全不影响其辐射性能。因此，针对阵列天线的辐射与散射综合优化有着重要的意义和广阔的应用前景。
近年来提出的RCS减缩技术主要有：外形技术、吸波材料技术、频率选择表面和阻抗加载技术等。外形技术通过合理设计目标外形和改善目标内部结构等方式减缩目标RCS；吸波材料技术通过对入射波电磁能量的消耗从而有效减少目标RCS；当双方雷达波处于不同工作频段时，可利用频率选择表面来实现对带外威胁雷达波的RCS减缩。这些传统的技术属于一种被动的技术，虽然发展迅速，应用范围很广，但也存在很多弊端。如吸波材料研制复杂、不便维护且对天线的辐射性能有负面影响；采用外形设计加工制造难度大，影响飞行器气动性能；频率选择表面对同频率同极化的威胁雷达波没有RCS减缩效果；阻抗加载技术只能在窄频带内降低天线的RCS。
由于阵列天线的辐射特性和低散射特性往往是一对难以调和的矛盾，现阶段关于阵列天线的辐射与散射综合优化尚处于初始研究阶段，相关理论与技术也相对匮乏，实现阵列天线良好辐射特性的同时有效降低阵列天线的RCS的研究报道极为有限。更进一步地，考虑到阵列天线的天然优势以及其在相控阵雷达系统的广泛应用，开展阵列天线的辐射与散射性能综合优化有着重要意义。例如，申请公布号为CN109950704A，名称为“一种用于强耦合宽带相控阵天线的带内RCS控制方法”的专利申请，公开了一种用于强耦合宽带相控阵天线的带内RCS控制方法，该方法主要基于阵列综合思想中的馈电相位优化的方法，在强耦合宽带相控阵天线各个单元端口与馈电网络之间接入一组不等长传输延时线的移相相位为优化变量，采用优化算法同时对阵列天线的辐射性能和散射性能进行优化，最终在辐射特性不被明显恶化的情况下，实现强耦合宽带相控阵天线带内RCS减缩。但存在的缺陷是接入了一段优化传输延时线，利用传输延时线相位贡献的阵列天线模式项RCS预估表达式对阵列天线模式项RCS进行预估近似。但是，阵列天线的RCS由天线模式项RCS和结构模式项RCS两部分组成，因阵列天线的各阵元馈电端口均匹配良好，天线模式项RCS的量级往往远低于结构模式项RCS。因此，该方法中的RCS控制方法仅能对阵列天线的天线模式项RCS进行优化控制，并不能降低结构模式项RCS，故而不能有效降低阵列天线的RCS。
发明内容
本发明的目的在于克服上述现有技术存在的缺陷，提出了一种基于电磁对消的阵列天线辐射与散射综合优化方法，旨在实现良好辐射特性的同时，有效降低阵列天线的RCS。
为实现上述目的，本发明采取的技术方案包括如下步骤：
(1)构建电磁对消场景模型：
在三维坐标系XYZ下构建放置在XOY平面内且包括N个周期性排布的天线单元的阵列天线A＝[A1,A2,…,An,…,AN]；在远区放置M个雷达T＝[T1,T2,…,Tm,…,TM]；设三维坐标系XYZ下一远区观察点P到坐标原点的距离为r，俯仰角为θ，方位角为
Figure BDA0002635895240000021
其中，N≥2，第n个天线单元An的馈电幅度、馈电相位和激励电流分别为αn、βn和In，
Figure BDA0002635895240000022
M≥1，第m个雷达Tm的方向为
Figure BDA0002635895240000023
(2)获取阵列天线A的辐射场
Figure BDA0002635895240000024
(2a)令N个天线单元的馈电相位均为0，天线单元An的馈电幅度αn＝1，其余N-1个天线单元的馈电幅度为0，并对An在远区观察点
Figure BDA0002635895240000025
处的阵中辐射场
Figure BDA0002635895240000026
进行测量：
Figure BDA0002635895240000027
其中，a表示辐射，B表示与天线单元An形式有关的比例系数，
Figure BDA0002635895240000028
表示天线单元An的归一化场强方向函数，
Figure BDA0002635895240000029
表示天线单元An相对于天线单元A1在空间产生场的相对相位；
(2b)对N个天线单元的阵中辐射场进行叠加，得到阵列天线A的辐射场
Figure BDA0002635895240000031
Figure BDA0002635895240000032
其中，∑表示求和；
(3)获取阵列天线A的散射场
Figure BDA0002635895240000033
(3a)令N个天线单元的馈电幅度与馈电相位均为0，当M＝1时，雷达T1在
Figure BDA0002635895240000034
方向以单位幅度均匀平面波对阵列天线A进行照射，则雷达T1探测到的阵列天线A的单站散射场为
Figure BDA0002635895240000035
Figure BDA0002635895240000036
其中，s表示散射，o表示单站，
Figure BDA0002635895240000037
表示阵列环境中An的远区散射场；
(3b)令N个天线单元的馈电幅度与馈电相位均为0，当M≥2时，雷达T1在
Figure BDA0002635895240000038
方向以单位幅度均匀平面波对阵列天线A进行照射，则其余M-1个雷达[T2，…,Tm,…,TM]在
Figure BDA0002635895240000039
方向探测到的阵列天线A的双站散射场为
Figure BDA00026358952400000310
Figure BDA00026358952400000311
其中，d表示双站，
Figure BDA00026358952400000312
表示阵列环境中An的远区散射场；
(4)计算阵列天线A的总场
Figure BDA00026358952400000313
Figure BDA00026358952400000314
其中，当M＝1时,
Figure BDA00026358952400000315
表示
Figure BDA00026358952400000316
当M≥2时，
Figure BDA00026358952400000317
表示
Figure BDA00026358952400000318
(5)定义优化算法的适应度函数fit(x)：
Figure BDA00026358952400000319
其中，x表示优化变量，x＝[α1,α2,…,αn,…,αN；β1,β2,…,βn,…,βN]T，[]T表示转置，max[Gain(x)]表示A的最大辐射方向增益，
Figure BDA0002635895240000041
表示A扫描过程中主波束指向，Psll表示峰值副瓣电平，Psll0表示对Psll进行抑制的期望值，
Figure BDA0002635895240000042
表示A的散射方向图在探测雷达方向上的雷达截面，σ0表示对A在探测雷达方向的雷达截面水平进行抑制的期望值，ω1，ω2，ω3和ω4表示权值系数；
(6)利用优化算法获取阵列天线A的辐射特性与散射特性的综合优化结果：
(6a)利用优化算法对优化变量x进行优化，得到优化后的优化变量为x*：
x*＝x′+x″
其中，x′＝[α′1,α′2,…,α′n,…,α′N；β′1,β′2,…,β′n,…,β′N]T，x″＝[α″1,α″2,…,α″n,…,α″N；β″1,β″2,…,β″n,…,β″N]T，α′n和α″n的取值范围为[0,1]，β′n和β″n的取值范围为[0,2π]；
(6b)将x′代入An的辐射场
Figure BDA0002635895240000043
得到优化后An的辐射场
Figure BDA0002635895240000044
并对所有天线单元优化后的辐射场进行叠加，得到阵列天线A优化后的辐射场
Figure BDA0002635895240000045
Figure BDA0002635895240000046
(6c)设定适应度函数fit(x*)中权值系数ω1，ω2，ω3和ω4的值，并利用优化算法，通过辐射场
Figure BDA0002635895240000047
对阵列天线A的散射场
Figure BDA0002635895240000048
进行电磁对消，实现对阵列天线A的单站RCS减缩、双站RCS减缩、全角域内的峰值RCS减缩以及整个工作频带内的RCS减缩；
(6d)将x″代入An的辐射场
Figure BDA0002635895240000049
得到An优化后的辐射场
Figure BDA00026358952400000410
对An优化后的辐射场
Figure BDA00026358952400000411
进行叠加，得到阵列天线A优化后的辐射场
Figure BDA00026358952400000412
Figure BDA00026358952400000413
(6e)设定适应度函数fit(x*)中权值系数ω1，ω2，ω3和ω4的值，并利用优化算法对辐射场
Figure BDA0002635895240000051
进行优化，实现但不限于最大辐射方向增益损失最小、峰值副瓣电平抑制在较低水平和扫描角最大程度的拓宽。
本发明与现有技术相比，具有如下优点：
本发明采用优化算法对每个天线单元的辐射场进行优化加权，产生一组阵列天线的辐射总场，该辐射总场分为两部分作用，其中一部分辐射场用来对阵列天线的散射场进行电磁对消，显著降低了阵列天线的RCS，另一部分用于优化阵列天线的辐射场，完成阵列天线的良好辐射特性，与现有技术相比，在实现良好辐射特性的同时，有效降低了阵列天线的RCS。
附图说明
图1为本发明的实现流程图；
图2为本发明的电磁对消场景模型图；
图3为本发明实施例的平面阵列天线的结构示意图；
图4为本发明实施例电磁对消前后沿
Figure BDA0002635895240000052
方向的照射雷达波入射时，
Figure BDA0002635895240000053
面内的双站RCS曲线图；
图5为本发明实施例电磁对消前后沿
Figure BDA0002635895240000054
方向的照射雷达波入射时，阵列天线全角域内的双站RCS峰值减缩曲线图；
图6为本发明实施例电磁对消前后沿
Figure BDA0002635895240000055
方向的照射雷达波入射时，阵列天线整个工作频带内的RCS减缩曲线图；
图7为本发明实施例的优化前后阵列天线主波束指向
Figure BDA0002635895240000056
处的辐射方向图；
图8为本发明实施例的优化前后阵列天线主波束指向
Figure BDA0002635895240000057
处的辐射方向图。
具体实施方式
以下结合附图和具体实施例，对本发明作进一步详细描述。
参照图1、本发明包括如下步骤：
步骤1)构建电磁对消场景模型：
参照图2，在三维坐标系XYZ下构建放置在XOY平面内且包括N个周期性排布的天线单元的阵列天线A＝[A1,A2,…,An,…,AN]；在远区放置M个雷达T＝[T1,T2,…,Tm,…,TM]；设三维坐标系XYZ下一远区观察点P到坐标原点的距离为r，俯仰角为θ，方位角为
Figure BDA0002635895240000061
其中，N≥2，本实施例中N＝80，第n个天线单元An的馈电幅度、馈电相位和激励电流分别为αn、βn和In，
Figure BDA0002635895240000062
M≥1，第m个雷达Tm的方向为
Figure BDA0002635895240000063
本发明适用于直线阵列、平面阵列和共形阵列等多种阵列形式，本实施例采用平面阵列天线，如图3所示，包含80个天线单元，从坐标原点处开始放置第一个天线单元，以x轴为列，y轴为行，每列等间距均匀放置5个单元，每行等间距均匀放置16个单元。
步骤2)获取阵列天线A的辐射场
Figure BDA0002635895240000064
(2a)令N个天线单元的馈电相位均为0，天线单元An的馈电幅度αn＝1，其余N-1个天线单元的馈电幅度为0，并对An在远区观察点
Figure BDA0002635895240000065
处的阵中辐射场
Figure BDA0002635895240000066
进行测量：
Figure BDA0002635895240000067
其中，a表示辐射，B表示与天线单元An形式有关的比例系数，
Figure BDA0002635895240000068
表示天线单元An的归一化场强方向函数，
Figure BDA0002635895240000069
表示天线单元An相对于天线单元A1在空间产生场的相对相位；
(2b)对N个天线单元的阵中辐射场进行叠加，得到阵列天线A的辐射场
Figure BDA00026358952400000610
Figure BDA00026358952400000611
其中，∑表示求和；
步骤3)获取阵列天线A的散射场
Figure BDA00026358952400000612
(3a)令N个天线单元的馈电幅度与馈电相位均为0，当M＝1时，雷达T1在
Figure BDA0002635895240000071
方向以单位幅度均匀平面波对阵列天线A进行照射，则雷达T1探测到的阵列天线A的单站散射场为
Figure BDA0002635895240000072
Figure BDA0002635895240000073
其中，s表示散射，o表示单站，
Figure BDA0002635895240000074
表示阵列环境中An的远区散射场；
(3b)令N个天线单元的馈电幅度与馈电相位均为0，当M≥2时，雷达T1在
Figure BDA0002635895240000075
方向以单位幅度均匀平面波对阵列天线A进行照射，则其余M-1个雷达[T2，…,Tm,…,TM]在
Figure BDA0002635895240000076
方向探测到的阵列天线A的双站散射场为
Figure BDA0002635895240000077
Figure BDA0002635895240000078
其中，d表示双站，
Figure BDA0002635895240000079
表示阵列环境中An的远区散射场；
步骤4)计算阵列天线A的总场
Figure BDA00026358952400000710
Figure BDA00026358952400000711
其中，当M＝1时,
Figure BDA00026358952400000712
表示
Figure BDA00026358952400000713
当M≥2时，
Figure BDA00026358952400000714
表示
Figure BDA00026358952400000715
步骤5)定义优化算法的适应度函数fit(x)：
Figure BDA00026358952400000716
其中，x表示优化变量，x＝[α1,α2,…,αn,…,αN；β1,β2,…,βn,…,βN]T，[]T表示转置，max[Gain(x)]表示A的最大辐射方向增益，
Figure BDA00026358952400000717
表示A扫描过程中主波束指向，Psll表示峰值副瓣电平，Psll0表示对Psll进行抑制的期望值，
Figure BDA00026358952400000718
表示A的散射方向图在探测雷达方向上的雷达截面，σ0表示对A在探测雷达方向的雷达截面水平进行抑制的期望值，ω1，ω2，ω3和ω4表示权值系数；优化算法可以为粒子群算法、蚁群算法、差分进化算法、空间映射算法和遗传算法等，本实施例采用遗传算法；
步骤6)利用遗传算法获取阵列天线A的辐射特性与散射特性的综合优化结果：
(6a)利用遗传算法对优化变量x进行优化，得到优化后的优化变量为x*：
x*＝x′+x″
其中，x′＝[α′1,α′2,…,α′n,…,α′N；β′1,β′2,…,β′n,…,β′N]T，x″＝[α″1,α″2,…,α″n,…,α″N；β″1,β″2,…,β″n,…,β″N]T，α′n和α″n的取值范围为[0,1]，β′n和β″n的取值范围为[0,2π]；本实施例利用遗传算法随机产生取值范围为[0,1]的[α′1,α′2,,α′n,,α′N]T、[α″1,α″2,…,α″n,…,α″N]T和取值范围为[0,2π]的[β′1,β′2,…,β′n,…,β′N]T、[β″1,β″2,…,β″n,…,β″NN]T，并通过选择、交叉、变异优化算子进一步优化，即可得到所有天线单元优化后的馈电幅度和馈电相位；
(6b)将x′代入An的辐射场
Figure BDA0002635895240000081
得到优化后An的辐射场
Figure BDA0002635895240000082
并对所有天线单元优化后的辐射场进行叠加，得到阵列天线A优化后的辐射场
Figure BDA0002635895240000083
Figure BDA0002635895240000084
本实施例产生的辐射场
Figure BDA0002635895240000085
与散射场
Figure BDA0002635895240000086
的幅度相同，相位相反，从而形成电磁对消；
(6c)设定适应度函数fit(x*)中权值系数ω1，ω2，ω3和ω4的值，并利用遗传算法，通过辐射场
Figure BDA0002635895240000087
对阵列天线A的散射场
Figure BDA0002635895240000088
进行电磁对消，实现对阵列天线A的单站RCS减缩、双站RCS减缩、全角域内的峰值RCS减缩以及整个工作频带内的RCS减缩；本实施例中权值系数的设定值为ω1＝1，ω2＝3，ω3＝2ω4＝1；
(6d)将x″代入An的辐射场
Figure BDA0002635895240000089
得到An优化后的辐射场
Figure BDA00026358952400000810
对An优化后的辐射场
Figure BDA00026358952400000811
进行叠加，得到阵列天线A优化后的辐射场
Figure BDA00026358952400000812
Figure BDA00026358952400000813
(6e)设定适应度函数fit(x*)中权值系数ω1，ω2，ω3和ω4的值，并利用遗传算法对辐射场
Figure BDA0002635895240000091
进行优化，实现了阵列天线A最大辐射方向增益损失在0.6dB以内、峰值副瓣电平抑制至-20dB以下、零波瓣宽度展宽小于2°和扫描角正负60°；本实施例权值系数的设定值为ω1＝1，ω2＝3，ω3＝2，ω4＝1。
以下结合实验结果，对本发明的技术效果作进一步说明：
1、实验条件和内容：
实验条件：Ansys 2019，Fortran 2013，Intel(R)Core(TM)i7-8700 CPU 3.20GHz，Windows 10(Pro)；
实验内容：沿
Figure BDA0002635895240000092
方向入射的单位幅度平面波对图3所示阵列天线进行照射，对阵列天线的辐射场与散射场进行测量，并利用遗传算法对每个天线单元辐射场进行优化加权，实现对阵列天线辐射与散射特性的综合优化。
2、实验结果分析：
参照图4，给出了电磁对消前后阵列天线的双站RCS(Bistatic RCS，BCS)结果对比图，可见通过电磁对消，在θ＝0°方向的单站RCS得到了完全抑制，实现了超过40dB的减缩。
参照图5，给出了电磁对消前后阵列天线全角域内RCS峰值减缩结果图，可见阵列的总双站RCS峰值出现在θ＝0°方向，幅值为-7.4dBsm，电磁对消后，阵列天线的双站RCS曲线基本保持在同一水平，幅值为-21.3dBsm，相比电磁对消前，实现了阵列天线双站RCS峰值减缩13.9dB，阵列天线在全角域范围内的RCS峰值得到了显著减缩。
参照图6，给出了电磁对消前后阵列天线在工作频带内的RCS减缩结果对比图，表明阵列天线相比电磁对消前可以在绝大部分工作频带内实现良好的减缩效果，且在20GHz频率入射波处获得了40.8dB的RCS减缩。
参照图7，给出了优化前后阵列天线主波束指向
Figure BDA0002635895240000093
的辐射特性结果对比图，可见通过优化，阵列天线的最大辐射方向保持不变、峰值副瓣电平抑制在-25dB以下，最大辐射方向增益仅损失0.5dB，具备良好的辐射特性。
参照图8，给出了优化前后阵列天线主波束指向
Figure BDA0002635895240000094
的辐射特性结果对比图，可见通过优化，阵列天线的最大辐射方向保持不变、峰值副瓣电平抑制在-25dB以下，最大辐射方向增益仅损失0.59dB，实现了阵列天线具备良好辐射特性的同时能够大角度扫描。
参考图4到图8的实验结果，可以看出，与现有技术相比本发明在实现良好辐射特性的同时，有效降低了阵列天线的RCS。

Claims (2)

1.一种基于电磁对消的阵列天线辐射与散射综合优化方法，其特征在于，包括如下步骤：
(1)构建电磁对消场景模型：
在三维坐标系XYZ下构建放置在XOY平面内且包括N个周期性排布的天线单元的阵列天线A＝[A1,A2,…,An,…,AN]；在远区放置M个雷达T＝[T1,T2,…,Tm,…,TM]；设三维坐标系XYZ下一远区观察点P到坐标原点的距离为r，俯仰角为θ，方位角为
Figure FDA0002635895230000011
其中，N≥2，第n个天线单元An的馈电幅度、馈电相位和激励电流分别为αn、βn和In，
Figure FDA0002635895230000012
第m个雷达Tm的方向为
Figure FDA0002635895230000013
(2)获取阵列天线A的辐射场
Figure FDA0002635895230000014
(2a)令N个天线单元的馈电相位均为0，天线单元An的馈电幅度αn＝1，其余N-1个天线单元的馈电幅度为0，并对An在远区观察点
Figure FDA0002635895230000015
处的阵中辐射场
Figure FDA0002635895230000016
进行测量：
Figure FDA0002635895230000017
其中，a表示辐射，B表示与天线单元An形式有关的比例系数，
Figure FDA0002635895230000018
表示天线单元An的归一化场强方向函数，
Figure FDA0002635895230000019
表示天线单元An相对于天线单元A1在空间产生场的相对相位；
(2b)对N个天线单元的阵中辐射场进行叠加，得到阵列天线A的辐射场
Figure FDA00026358952300000110
Figure FDA00026358952300000111
其中，∑表示求和；
(3)获取阵列天线A的散射场
Figure FDA0002635895230000021
(3a)令N个天线单元的馈电幅度与馈电相位均为0，当M＝1时，雷达T1在
Figure FDA0002635895230000022
方向以单位幅度均匀平面波对阵列天线A进行照射，则雷达T1探测到的阵列天线A的单站散射场为
Figure FDA0002635895230000023
Figure FDA0002635895230000024
其中，s表示散射，o表示单站，
Figure FDA0002635895230000025
表示阵列环境中An的远区散射场；
(3b)令N个天线单元的馈电幅度与馈电相位均为0，当M≥2时，雷达T1在
Figure FDA0002635895230000026
方向以单位幅度均匀平面波对阵列天线A进行照射，则其余M-1个雷达[T2，…,Tm,…,TM]在
Figure FDA0002635895230000027
方向探测到的阵列天线A的双站散射场为
Figure FDA0002635895230000028
Figure FDA0002635895230000029
其中，d表示双站，
Figure FDA00026358952300000210
表示阵列环境中An的远区散射场；
(4)计算阵列天线A的总场
Figure FDA00026358952300000211
Figure FDA00026358952300000212
其中，当M＝1时,
Figure FDA00026358952300000213
表示
Figure FDA00026358952300000214
当M≥2时，
Figure FDA00026358952300000215
表示
Figure FDA00026358952300000216
(5)定义优化算法的适应度函数fit(x)：
Figure FDA00026358952300000217
其中，x表示优化变量，x＝[α1,α2,…,αn,…,αN；β1,β2,…,βn,…,βN]T，[]T表示转置，max[Gain(x)]表示A的最大辐射方向增益，
Figure FDA0002635895230000031
表示A扫描过程中主波束指向，Psll表示峰值副瓣电平，Psll0表示对Psll进行抑制的期望值，
Figure FDA0002635895230000032
表示A的散射方向图在探测雷达方向上的雷达截面，σ0表示对A在探测雷达方向的雷达截面水平进行抑制的期望值，ω1，ω2，ω3和ω4表示权值系数；
(6)利用优化算法获取阵列天线A的辐射特性与散射特性的综合优化结果：
(6a)利用优化算法对优化变量x进行优化，得到优化后的优化变量为x*：
x*＝x′+x″
其中，x'＝[α′1,α′2,…,α′n,…,α′N；β′1,β′2,…,β′n,…,β′N]T，x″＝[α″1,α″2,…,α″n,…,α″N；β″1,β″2,…,β″n,…,β″N]T，α′n和α″n的取值范围为[0,1]，β′n和β″n的取值范围为[0,2π]；
(6b)将x'代入An的辐射场
Figure FDA0002635895230000033
得到优化后An的辐射场
Figure FDA0002635895230000034
并对所有天线单元优化后的辐射场进行叠加，得到阵列天线A优化后的辐射场
Figure FDA0002635895230000035
Figure FDA0002635895230000036
(6c)设定适应度函数fit(x*)中权值系数ω1，ω2，ω3和ω4的值，并利用优化算法，通过辐射场
Figure FDA0002635895230000037
对阵列天线A的散射场
Figure FDA0002635895230000038
进行电磁对消，实现对阵列天线A的单站RCS减缩、双站RCS减缩、全角域内的峰值RCS减缩以及整个工作频带内的RCS减缩；
(6d)将x”代入An的辐射场
Figure FDA0002635895230000039
得到An优化后的辐射场
Figure FDA00026358952300000310
对An优化后的辐射场
Figure FDA00026358952300000311
进行叠加，得到阵列天线A优化后的辐射场
Figure FDA00026358952300000312
Figure FDA0002635895230000041
(6e)设定适应度函数fit(x*)中权值系数ω1，ω2，ω3和ω4的值，并利用优化算法对辐射场
Figure FDA0002635895230000042
进行优化，实现但不限于最大辐射方向增益损失最小、峰值副瓣电平抑制在较低水平和扫描角最大程度的拓宽。
2.根据权利要求1所述的基于电磁对消的阵列天线辐射与散射综合优化方法，其特征在于，步骤(5)中所述的A的散射方向图在探测雷达方向上的雷达截面
Figure FDA0002635895230000043
计算公式为：
Figure FDA0002635895230000044
其中，
Figure FDA0002635895230000045
表示照射到阵列天线A上的入射波强度，
Figure FDA0002635895230000046
表示阵列天线A在该入射波照射下产生的散射强度，r→∞表示在远场区。
CN202010825165.8A 2020-08-17 2020-08-17 基于电磁对消的阵列天线辐射与散射综合优化方法 Active CN111984911B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010825165.8A CN111984911B (zh) 2020-08-17 2020-08-17 基于电磁对消的阵列天线辐射与散射综合优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010825165.8A CN111984911B (zh) 2020-08-17 2020-08-17 基于电磁对消的阵列天线辐射与散射综合优化方法

Publications (2)

Publication Number Publication Date
CN111984911A true CN111984911A (zh) 2020-11-24
CN111984911B CN111984911B (zh) 2023-11-14

Family

ID=73434498

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010825165.8A Active CN111984911B (zh) 2020-08-17 2020-08-17 基于电磁对消的阵列天线辐射与散射综合优化方法

Country Status (1)

Country Link
CN (1) CN111984911B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113258991A (zh) * 2021-04-13 2021-08-13 西北大学 一种散射通信系统、微带天线、阵列及通信距离提升方法
CN113964549A (zh) * 2021-12-22 2022-01-21 中国人民解放军海军工程大学 一种基于干扰对消的空间取样天线的设计方法及装置
CN114880616A (zh) * 2022-05-17 2022-08-09 中国人民解放军63660部队 一种获取瞬态电磁脉冲阵列天线辐射场的方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5353035A (en) * 1990-04-20 1994-10-04 Consejo Superior De Investigaciones Cientificas Microstrip radiator for circular polarization free of welds and floating potentials
CN104182636A (zh) * 2014-08-22 2014-12-03 西安电子科技大学 一种阵列天线辐射场和散射场综合低副瓣快速实现方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5353035A (en) * 1990-04-20 1994-10-04 Consejo Superior De Investigaciones Cientificas Microstrip radiator for circular polarization free of welds and floating potentials
CN104182636A (zh) * 2014-08-22 2014-12-03 西安电子科技大学 一种阵列天线辐射场和散射场综合低副瓣快速实现方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张帅;龚书喜;关莹;龚琦;: "考虑互耦的半波振子线阵辐射和散射方向图综合", 计算物理, no. 03 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113258991A (zh) * 2021-04-13 2021-08-13 西北大学 一种散射通信系统、微带天线、阵列及通信距离提升方法
CN113258991B (zh) * 2021-04-13 2023-02-14 西北大学 一种散射通信系统、微带天线、阵列及通信距离提升方法
CN113964549A (zh) * 2021-12-22 2022-01-21 中国人民解放军海军工程大学 一种基于干扰对消的空间取样天线的设计方法及装置
CN114880616A (zh) * 2022-05-17 2022-08-09 中国人民解放军63660部队 一种获取瞬态电磁脉冲阵列天线辐射场的方法
CN114880616B (zh) * 2022-05-17 2024-04-05 中国人民解放军63660部队 一种获取瞬态电磁脉冲阵列天线辐射场的方法

Also Published As

Publication number Publication date
CN111984911B (zh) 2023-11-14

Similar Documents

Publication Publication Date Title
CN111984911B (zh) 基于电磁对消的阵列天线辐射与散射综合优化方法
CN109950704B (zh) 一种用于强耦合宽带相控阵天线的带内rcs控制方法
CN111352080B (zh) Papr和相似性约束下低截获频控阵mimo雷达系统的设计方法
CN112151969B (zh) 基于广义散射矩阵的强耦合宽带相控阵带内rcs控制方法
Hao et al. Small-size broadband coding metasurface for RCS reduction based on particle swarm optimization algorithm
Cui et al. Beam-space reduced-dimension 3D-STAP for nonside-looking airborne radar
CN109560389B (zh) 基于新型散射对消方法的宽带rcs减缩超表面
Panduro et al. Energy Patterns of UWB antenna arrays with low side lobe level during beam-scanning
CN113126041B (zh) 基于惩罚因子可变的雷达分布式干扰信号产生方法
CN111710992B (zh) 一种基于等效线阵的弧形阵列天线的波束扫描方法及装置
CN113964548A (zh) 一种基于四维天线阵的带内散射缩减结构及其控制方法
CN111257871B (zh) 可用于微波凝视关联成像的单天线辐射源设计方法
Zhou et al. Efficient optimization and realization of a shaped-beam planar array for very large array application
CN106934102B (zh) 基于迭代傅里叶算法设计宽带随机表面的方法
CN112526510A (zh) 基于方向图分集的单通道角度超分辨方法及系统
Rocca et al. Synthesis of arbitrary sidelobes sum and difference patterns with common excitation weights
Liu et al. Recent development in nonuniformly spaced array synthesis methods
Atamanyuk et al. Nonuniform absorbing coating as an effective way to reduce radar visibility of an object with the surface formed by flat conductive plates
Zuo et al. Calculation of the extinction coefficient of dipoles cloud
CN113721210B (zh) 基于吸波-对消的深度rcs减缩超表面设计方法及超表面
CN116130980B (zh) 一种用于混合模式涡旋电磁波的相控电磁表面设计方法
CN113126039B (zh) 基于tch分解的stap雷达分布式干扰信号产生方法
Chopde et al. A Beam Summation Scheme for Ultra-wideband RCS Calculations in the High–frequency Regime
Li et al. Research on RCS Reduction of the Phased Array Antennas
Pu et al. A Non-adjacent Feed Selection Method for Beamforming Performance Enhancement for Multi-beam Satellite Antenna

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant