CN107436970A - 基于序列二次规划的三坐标雷达赋形方向图设计方法 - Google Patents
基于序列二次规划的三坐标雷达赋形方向图设计方法 Download PDFInfo
- Publication number
- CN107436970A CN107436970A CN201710544762.1A CN201710544762A CN107436970A CN 107436970 A CN107436970 A CN 107436970A CN 201710544762 A CN201710544762 A CN 201710544762A CN 107436970 A CN107436970 A CN 107436970A
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- mtd
- main lobe
- represent
- 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
Links
Classifications
-
- 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
-
- 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)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及一种基于序列二次规划的三坐标雷达赋形方向图设计方法，主要解决现有赋形方向图设计方法没有同时兼顾主瓣波动与旁瓣抑制、设计时间较长等问题。其实现过程是：1.设定期望方向图，设定方向图波动上限与下限；2.设定主瓣权重与旁瓣权重；3.产生多组随机相位权，并利用交替投影算法优化得到初始相位权；4.利用序列二次规划得到多组优化相位权；5.从多组优化相位权中选择最优权。本发明可以在一定程度上控制主瓣波动电平，同时可以抑制旁瓣电平，优化效率较高，具有工程实用性。
Description
技术领域
本发明属于雷达技术领域，涉及赋形方向图设计，可用于三坐标雷达不同工作模式下的发射赋形设计，并满足工程要求。
背景技术
阵列天线可以通过调整幅度、相位实现高增益窄波束或者特定的赋形波束，而且波束扫描方式灵活，因此在雷达、通信等领域具有广泛应用。在雷达领域中，发射方向图直接关系到目标信号的强度、杂波与干扰的抑制，也决定了雷达的威力，从而影响雷达***的探测性能和资源利用率。另外，在不同场景中，通常需要不同形状的波束，如空情三坐标雷达的搜索模式，通常需要余割平方的赋形波束，跟踪模式通常采用窄波束。因此，需要根据工作模式、场景的不同，设计不同的赋形方向图。
在雷达***中，要求发射方向图具有特定的主瓣形状和低旁瓣电平等。对于赋形方向图综合，目前主要包括：解析方法、智能优化算法和迭代优化算法等。解析方法，具有计算量小的特点，但仅局限于某些特定形式的方向图或性能有限。某些方向图综合问题可以转化为凸优化问题，从而可以基于优化理论进行求解，获得全局最优解。以上方法一般需要对权系数的幅度和相位同时进行优化，然而为了最大化发射机效率，需要保证权系数幅度恒定，即只对其相位进行优化。智能优化方法，如遗传算法等通过构造优化代价函数，对方向图的权向量进行优化以逼近期望的赋形方向图，通常计算量很大，而且需要合适的参数设定，对于大规模阵列难以应用。此外，交替投影、迭代FFT方法采用迭代的方式进行方向图综合，可以对方向图逐渐调整，达到较理想的结果。针对低旁瓣方向图综合，Dolph-Chebyshev方法和Taylor综合法可以获得较理想的旁瓣电平，但方向图主瓣形状较单一。
三坐标雷达方向图的主瓣波动，不仅会影响雷达的探测威力，也会影响目标的测量精度及跟踪性能，而过高的旁瓣电平会产生不必要的杂波及干扰，造成虚警或掩盖小目标。现有针对相位优化的方向图设计算法对主瓣形状控制及旁瓣抑制无法兼顾，因此设计的发射方向图主瓣通常存在较大波动，或者旁瓣电平不理想。
发明内容
要解决的技术问题
为了避免现有技术的不足之处，本发明提出一种基于序列二次规划的三坐标雷达赋形方向图设计方法对阵列相位权进行优化，以保证所设计方向图逼近期望方向图，而且具有较小的主瓣波动和旁瓣电平。
技术方案
一种基于序列二次规划的三坐标雷达赋形方向图设计方法，其特征在于步骤如下：
步骤1：根据要求设定期望方向图Pd(θn)，设定方向图波动的上限P1(θn)＝(1+δ1)Pd(θn)，n＝1,…,Ns与下限P2(θn)＝(1-δ2)Pd(θn)，n＝1,…,Ns，其中δ1表示上波纹波动参数，δ2表示下波纹波动参数，为归一化空间频率域均匀采样的离散角度，Ns表示空间采样点数；
步骤2：设定主瓣区域Ωm＝[θMin,θMax]，旁瓣区域Ωs＝[-90,θMin)∪(θMax,90]，其中θMax与θMin分别表示为主瓣的最大角度与最小角度；
根据主瓣区域方向图设定主瓣权重μm(θn)为
根据旁瓣区域方向图设定旁瓣权重μs(θn)为
步骤3：产生K组相位在[0,2π]上均匀分布的随机相位权k＝1,…,K，利用交替投影算法得到K组初始相位权k＝1,…,K；
步骤4：以主瓣波动和旁瓣电平抑制为准则，建立关于相位权的非线性优化模型：
s.t.αL≤α≤αH
其中相位权在θn方向上的方向图表示为a(θn)表示θn方向上的导向矢量，表示以相位权中元素为相位的恒模向量，符号j表示虚数单位，符号exp(·)表示以自然对数e为底的指数函数，(·)H表示共轭转置，η表示旁瓣抑制权重，α表示尺度因子，αL与αH分别表示尺度因子的上界与下界，是方向图主瓣逼近的代价函数，其反映了相位权所形成的方向图在αP1(θn)与αP2(θn)之外波纹大小，具体表示为：
分别以每组初始相位权k＝1,…,K为初始解，利用序列二次规划得到K组优化相位权k＝1,…,K；
步骤5：从K组优化相位权k＝1,…,K中选择使步骤4中优化模型的代价函数最小的相位权作为最优相位权。
步骤1所述的根据要求设定期望方向图，是按照如下方式确定的：
将期望方向图主瓣划分为等高度探测区域、远距离探测区域和低空探测区域，主瓣外的区域为旁瓣区域，设期望的探测高度为Hs，最远距离为R，天线上仰角度为γ，主瓣的最大角度与最小角度分别表示为θMax与θMin，远距离探测区域的起始角度为θE，主瓣区域最小增益为Pmin，旁瓣区域最小增益与最大增益分别为PSL1和PSL2，因此可设定期望方向图为Pd(θn)，具体表示如下
其中等高度探测区域为[θC-γ,θMax-γ]，远距离探测区域为[θE-γ,θC-γ]，低空探测区域为[θMin-γ,θE-γ]，csc(·)表示余割函数，衰减参数Q∈[2,4]表示余割函数的幂次，用于控制衰减，为等高度探测区域的起始角度，arcsin(·)表示反正弦函数，参数ξ1、ξ2表示两个旁瓣区域方向图斜率，参数ξ表示低空区域方向图斜率，具体定义为
ξ1＝[10log(PSL2)-10log(PSL1)]/[θMin-γ+90]
ξ2＝[10log(PSL2)-10log(PSL1)]/[θMax-γ-90]
ξ＝[10log(PMin)-10log(P0(θE))]/[θMin-θE]
其中log(·)表示10为底的对数，数值P0(θn)、P0(θC)和P0(θE)分别表示窄波束P0(θ)在θn、θC和θE处的取值，窄波束P0(θ)是角度θ的函数，定义为P0(θ)＝|bT(θ)b*(θC/2+θE/2-γ)|2，向量表示窄波束P0(θ)对应的导向矢量，|·|表示取模，表示形成窄波束P0(θn)所需的阵元个数，表示向下取整，(·)T表示转置，(·)*表示共轭，表示归一化空间频率域均匀采样的离散角度，Ns表示空间采样点数。
有益效果
本发明提出的一种基于序列二次规划的三坐标雷达赋形方向图设计方法，可以在一定程度上控制主瓣波动电平，同时可以抑制旁瓣电平；以交替投影算法的优化结果为初始相位权，并采用多组相位权对序列二次规划进行初始化，可以防止相位权收敛于性能较差的局部解，而且优化效率较高；只以相位作为优化变量，而各阵元激励信号的幅度相等，可以保证各阵元放大器工作于饱和放大状态，符合现有雷达***的工作方式，具有工程应用性。
附图说明
图1是本发明的实现流程图
图2是本发明与交替投影算法所设计的高仰角探测方向图比较
图3是本发明与交替投影算法所设计的下视探测方向图比较
具体实施方式
现结合实施例、附图对本发明作进一步描述：
参照图1，本实施例的具体实现步骤如下：
步骤1.设定期望方向图，设定方向图波动上限与下限。
假设雷达***是由M个阵元构成的均匀线阵，定义相位权为
其中φm表示对第m个阵元发射信号所需调节的相位，(·)T表示转置，则相位权在远场θ方向处的方向图表示为：
其中a(θ)＝[1,ej2πdsinθ/λ,…,ej2πd(M-1)sinθ/λ]T表示导向矢量，θ表示俯仰角，d表示阵元间距，λ表示信号波长，(·)H表示共轭转置，(·)*表示共轭，sin(·)表示正弦函数，表示以相位权中元素为相位的恒模向量，符号j表示虚数单位，符号exp(·)表示以自然对数e为底的指数函数；
确定归一化空间频率域均匀采样的离散角度具体定义如下：
其中Ns表示空间采样点数，arcsin(·)表示反正弦函数；
期望方向图主瓣包括等高度探测区域、远距离探测区域和低空探测区域，主瓣外的区域为旁瓣区域，设期望的探测高度为Hs，最远距离为R，天线上仰角度为γ，主瓣的最大角度与最小角度分别表示为θMax与θMin，远距离探测区域的起始角度为θE，主瓣区域最小增益为PMin，旁瓣区域最小增益与最大增益分别为PSL1和PSL2，因此可设定期望方向图为Pd(θn)，具体表示如下：
其中等高度探测区域为[θC-γ,θMax-γ]，远距离探测区域为[θE-γ,θC-γ]，低空探测区域为[θMin-γ,θE-γ]，csc(·)表示余割函数，衰减参数Q∈[2,4]表示余割函数的指数，用于控制衰减，为等高度探测区域的起始角度，参数ξ1、ξ2表示两个旁瓣区域方向图斜率，参数ξ表示低空区域方向图斜率，具体定义为：
ξ1＝[10log(PSL2)-10log(PSL1)]/[θMin-γ+90]
ξ2＝[10log(PSL2)-10log(PSL1)]/[θMax-γ-90]
ξ＝[10log(PMin)-10log(P0(θE))]/[θMin-θE]
其中log(·)表示10为底的对数，数值P0(θn)、P0(θC)和P0(θE)分别表示窄波束P0(θ)在θn、θC和θE处的取值，窄波束P0(θ)是角度θ的函数，定义为P0(θ)＝|bT(θ)b*(θC/2+θE/2-γ)|2，向量表示窄波束P0(θ)对应的导向矢量，|·|表示取模，表示形成窄波束P0(θn)所需的阵元个数，表示向下取整；
设定方向图波动的上限P1(θn)与下限P2(θn)，具体表示如下
P1(θn)＝(1+δ1)Pd(θn)，n＝1,…,Ns
P2(θn)＝(1-δ2)Pd(θn)，n＝1,…,Ns
其中δ1表示上波纹波动参数，δ2表示下波纹波动参数。
步骤2.设定主瓣权重与旁瓣权重。
设定主瓣区域Ωm＝[θMin,θMax]，旁瓣区域Ωs＝[-90,θMin)∪(θMax,90]；
根据主瓣区域方向图可设定主瓣权重为
根据旁瓣区域方向图可设定旁瓣权重为
步骤3.产生K组随机相位权，并利用交替投影算法优化得到初始相位权。
随机产生K组相位权k＝1,…,K，其中为M×1的向量，各元素服从[0,2π]上的均匀分布；
根据随机相位权利用交替投影算法计算得到K组初始相位权，其中第k组初始相位权按如下步骤计算，k＝1,…,K：
(3.1)设定交替投影算法的初始解为设定终止阈值ε，计算初始空域信号z(0)＝Fexp(jψ(0))，令t＝0；
其中表示Ns点IFFT矩阵，其第n列定义为
则Ns点FFT矩阵为FH
(3.2)计算期望空域信号y(t)，其维数为Ns×1，第n个元素为
其中arg(·)表示取相位；
(3.3)计算新解ψ(t+1)，其第m个元素为
(3.4)计算新解ψ(t+1)对应的空域信号z(t+1)＝Fexp(jψ(t+1))，判断终止条件是否成立，是则停止并令初始相位权否则令t＝t+1并重复步骤(3.2)-(3.4)。
步骤4.利用序列二次规划得到优化相位权。
(4.1)以主瓣波动与旁瓣电平抑制为准则，建立关于相位权的非线性优化模型：
根据步骤2设定的主瓣区域与旁瓣区域建立如下优化模型：
s.t.αL≤α≤αH <5>
其中相位权在θn方向上的方向图定义如式<1>所示，η表示旁瓣抑制权重，α表示尺度因子，αL与αH分别表示尺度因子的上界与下界，是方向图主瓣逼近的代价函数，其反映了相位权所形成的方向图在αP1(θn)与αP2(θn)之外的波纹大小，具体表示为：
n＝1,…,Ns，
其中P1(θn)表示主瓣波动的上界，P2(θn)表示主瓣波动的下界；
(4.2)以每组初始相位权为初始解，利用序列二次规划求解式<5>中非线性优化模型，得到K组优化相位权k＝1,…,K；
其中序列二次规划见[Antoniou A.,Lu W.,Practical Optimization Algorithmand Engineering Applications,Springer,2007.]。
步骤5.选择最优权。
从K组优化相位权k＝1,…,K中选择使式<5>中代价函数最小的相位权，作为最优相位权。
本发明的效果通过以下仿真对比试验进一步说明：
1.仿真条件：假设雷达***俯仰维包括M＝48个阵元，阵元间距为半波长，天线上仰角度为γ＝10°，空间采样点数Ns＝256，上波纹波动参数δ1＝0.02，下波纹波动参数δ2＝0.02，主瓣区域最小增益为PMin＝-15dB，旁瓣区域最小增益与最大增益为PSL1＝-40dB和PSL2＝-22dB，交替投影算法终止阈值ε＝10-3，相位权个数为K＝30，式<5>中旁瓣抑制权重μ＝1，尺度因子的上界与下界分别为 设置序列二次规划迭代次数为200。
2.仿真内容：
实验1：高仰角探测方向图
主瓣最大角度θMax＝50°，主瓣最小角度θMin＝0°，最大探测距离为R＝100km，最大探测高度为Hs＝20km，远距离探测起始角度为θE＝5°，衰减参数Q＝3，期望方向图如图2中点线所示，交替投影算法与本发明所设计方向图如图2中虚线和实线所示。
实验2：下视探测方向图
主瓣最大角度θMax＝30°，主瓣最小角度θMin＝-10°，最大探测距离为R＝100km，最大探测高度为Hs＝15km，远距离探测起始角度为θE＝2°，衰减参数Q＝3，期望方向图如图3中点线所示，交替投影算法与本发明所设计方向图如图3中虚线和实线所示。
3.仿真结果分析：
从图2中虚线可以看出，所设定期望方向图在高角度区域增益逐渐衰减，具有高仰角探测的能力。对比图2中实线与点线可以看出，相比交替投影方法，本发明方法的主瓣波动明显减小，且更加逼近期望方向图，虽然本发明方法设计的方向图旁瓣电平有一定抬高，考但虑实际中单个阵元方向图在偏离法向方向存在一定衰减，因此仍然可以对旁瓣有较好的抑制效果。
从图3中虚线可以看出，所设定期望方向图不仅在高角度区域存在一定增益，在低角度区域增益也是逐渐衰减的，具有下视探测的能力。对比图3中实线与点线可以看出，与交替投影方法相比，本发明方法的旁瓣电平有一定抬高，但主瓣波动明显减小，且更逼近期望方向图。
Claims (2)
1.一种基于序列二次规划的三坐标雷达赋形方向图设计方法，其特征在于步骤如下：
步骤1：根据要求设定期望方向图Pd(θn)，设定方向图波动的上限P1(θn)＝(1+δ1)Pd(θn)，n＝1,…,Ns与下限P2(θn)＝(1-δ2)Pd(θn)，n＝1,…,Ns，其中δ1表示上波纹波动参数，δ2表示下波纹波动参数，为归一化空间频率域均匀采样的离散角度，Ns表示空间采样点数；
步骤2：设定主瓣区域Ωm＝[θMin,θMax]，旁瓣区域Ωs＝[-90,θMin)∪(θMax,90]，其中θMax与θMin分别表示为主瓣的最大角度与最小角度；
根据主瓣区域方向图设定主瓣权重μm(θn)为
<mrow>
<msub>
<mi>&mu;</mi>
<mi>m</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>P</mi>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>&theta;</mi>
<mi>n</mi>
</msub>
<mo>&Element;</mo>
<msub>
<mi>&Omega;</mi>
<mi>m</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>&theta;</mi>
<mi>n</mi>
</msub>
<mo>&NotElement;</mo>
<msub>
<mi>&Omega;</mi>
<mi>m</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>N</mi>
<mi>s</mi>
</msub>
<mo>;</mo>
</mrow>
根据旁瓣区域方向图设定旁瓣权重μs(θn)为
<mrow>
<msub>
<mi>&mu;</mi>
<mi>s</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>P</mi>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mn>4</mn>
</mrow>
</msup>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>&theta;</mi>
<mi>n</mi>
</msub>
<mo>&Element;</mo>
<msub>
<mi>&Omega;</mi>
<mi>s</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>&theta;</mi>
<mi>n</mi>
</msub>
<mo>&NotElement;</mo>
<msub>
<mi>&Omega;</mi>
<mi>s</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>N</mi>
<mi>s</mi>
</msub>
<mo>;</mo>
</mrow>
步骤3：产生K组相位在[0,2π]上均匀分布的随机相位权利用交替投影算法得到K组初始相位权
步骤4：以主瓣波动和旁瓣电平抑制为准则，建立关于相位权的非线性优化模型：
s.t.αL≤α≤αH
其中相位权在θn方向上的方向图表示为a(θn)表示θn方向上的导向矢量，表示以相位权中元素为相位的恒模向量，符号j表示虚数单位，符号exp(·)表示以自然对数e为底的指数函数，(·)H表示共轭转置，η表示旁瓣抑制权重，α表示尺度因子，αL与αH分别表示尺度因子的上界与下界，是方向图主瓣逼近的代价函数，其反映了相位权所形成的方向图在αP1(θn)与αP2(θn)之外波纹大小，具体表示为：
分别以每组初始相位权为初始解，利用序列二次规划得到K组优化相位权
步骤5：从K组优化相位权中选择使步骤4中优化模型的代价函数最小的相位权作为最优相位权。
2.根据权利要求1所述的一种基于序列二次规划的三坐标雷达赋形方向图设计方法，其特征在于步骤1所述的根据要求设定期望方向图，是按照如下方式确定的：
将期望方向图主瓣划分为等高度探测区域、远距离探测区域和低空探测区域，主瓣外的区域为旁瓣区域，设期望的探测高度为Hs，最远距离为R，天线上仰角度为γ，主瓣的最大角度与最小角度分别表示为θMax与θMin，远距离探测区域的起始角度为θE，主瓣区域最小增益为Pmin，旁瓣区域最小增益与最大增益分别为PSL1和PSL2，因此可设定期望方向图为Pd(θn)，具体表示如下
其中等高度探测区域为[θC-γ,θMax-γ]，远距离探测区域为[θE-γ,θC-γ]，低空探测区域为[θMin-γ,θE-γ]，csc(·)表示余割函数，衰减参数Q∈[2,4]表示余割函数的幂次，用于控制衰减，为等高度探测区域的起始角度，arcsin(·)表示反正弦函数，参数ξ1、ξ2表示两个旁瓣区域方向图斜率，参数ξ表示低空区域方向图斜率，具体定义为
ξ1＝[10log(PSL2)-10log(PSL1)]/[θMin-γ+90]
ξ2＝[10log(PSL2)-10log(PSL1)]/[θMax-γ-90]
ξ＝[10log(PMin)-10log(P0(θE))]/[θMin-θE]
其中log(·)表示10为底的对数，数值P0(θn)、P0(θC)和P0(θE)分别表示窄波束P0(θ)在θn、θC和θE处的取值，窄波束P0(θ)是角度θ的函数，定义为P0(θ)＝|bT(θ)b*(θC/2+θE/2-γ)|2，向量表示窄波束P0(θ)对应的导向矢量，|·|表示取模，表示形成窄波束P0(θn)所需的阵元个数，表示向下取整，(·)T表示转置，(·)*表示共轭，表示归一化空间频率域均匀采样的离散角度，Ns表示空间采样点数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710544762.1A CN107436970B (zh) | 2017-07-06 | 2017-07-06 | 基于序列二次规划的三坐标雷达赋形方向图设计方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710544762.1A CN107436970B (zh) | 2017-07-06 | 2017-07-06 | 基于序列二次规划的三坐标雷达赋形方向图设计方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107436970A true CN107436970A (zh) | 2017-12-05 |
| CN107436970B CN107436970B (zh) | 2020-04-28 |
Family
ID=60459750
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710544762.1A Active CN107436970B (zh) | 2017-07-06 | 2017-07-06 | 基于序列二次规划的三坐标雷达赋形方向图设计方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107436970B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110501675A (zh) * | 2019-07-16 | 2019-11-26 | 北京工业大学 | 一种基于mimo雷达低旁瓣发射方向图设计方法 |
| CN115034075A (zh) * | 2022-06-22 | 2022-09-09 | 电子科技大学 | 非规则天线阵列矢量增益方向图的快速精准预测方法 |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4199168B2 (ja) * | 2004-07-14 | 2008-12-17 | 株式会社東芝 | アンテナ装置 |
| JP4564407B2 (ja) * | 2005-05-30 | 2010-10-20 | 株式会社東芝 | アンテナ装置 |
| CN101950014A (zh) * | 2010-09-02 | 2011-01-19 | 西安电子科技大学 | Mimo雷达波形在线设计方法 |
| CN102841335A (zh) * | 2012-09-10 | 2012-12-26 | 西安电子科技大学 | 基于迭代fft的mimo雷达快速波形合成方法 |
| CN103018732A (zh) * | 2013-01-17 | 2013-04-03 | 西安电子科技大学 | 基于空时联合优化的mimo雷达波形合成方法 |
| CN103810325A (zh) * | 2014-01-08 | 2014-05-21 | 西安电子科技大学 | 基于序列二次规划的低旁瓣线性稀布阵列天线优化方法 |
| CN104992000A (zh) * | 2015-06-18 | 2015-10-21 | 哈尔滨工业大学 | 一种基于l型阵列天线的波束形成及波束图优化方法 |
| CN106712864A (zh) * | 2017-01-20 | 2017-05-24 | 京信通信技术(广州)有限公司 | 一种智能天线性能测试及优化的方法及装置 |
| CN106772260A (zh) * | 2017-03-31 | 2017-05-31 | 西安电子科技大学 | 基于凸优化算法的雷达阵列和差波束方向图优化方法 |
| CN106842237A (zh) * | 2017-01-18 | 2017-06-13 | 南京理工大学 | 快速任意形状方向图主瓣保形自适应波束形成方法 |
-
2017
- 2017-07-06 CN CN201710544762.1A patent/CN107436970B/zh active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4199168B2 (ja) * | 2004-07-14 | 2008-12-17 | 株式会社東芝 | アンテナ装置 |
| JP4564407B2 (ja) * | 2005-05-30 | 2010-10-20 | 株式会社東芝 | アンテナ装置 |
| CN101950014A (zh) * | 2010-09-02 | 2011-01-19 | 西安电子科技大学 | Mimo雷达波形在线设计方法 |
| CN102841335A (zh) * | 2012-09-10 | 2012-12-26 | 西安电子科技大学 | 基于迭代fft的mimo雷达快速波形合成方法 |
| CN103018732A (zh) * | 2013-01-17 | 2013-04-03 | 西安电子科技大学 | 基于空时联合优化的mimo雷达波形合成方法 |
| CN103810325A (zh) * | 2014-01-08 | 2014-05-21 | 西安电子科技大学 | 基于序列二次规划的低旁瓣线性稀布阵列天线优化方法 |
| CN104992000A (zh) * | 2015-06-18 | 2015-10-21 | 哈尔滨工业大学 | 一种基于l型阵列天线的波束形成及波束图优化方法 |
| CN106842237A (zh) * | 2017-01-18 | 2017-06-13 | 南京理工大学 | 快速任意形状方向图主瓣保形自适应波束形成方法 |
| CN106712864A (zh) * | 2017-01-20 | 2017-05-24 | 京信通信技术(广州)有限公司 | 一种智能天线性能测试及优化的方法及装置 |
| CN106772260A (zh) * | 2017-03-31 | 2017-05-31 | 西安电子科技大学 | 基于凸优化算法的雷达阵列和差波束方向图优化方法 |
Non-Patent Citations (6)
| Title |
|---|
| A.HAMMAMI ET AL: "Planar array antenna pattern nulling based on sequential quadratic programming (SQP) algorithm", 《 EIGHTH INTERNATIONAL MULTI-CONFERENCE ON SYSTEMS, SIGNALS & DEVICES》 * |
| FRANCESCO BELFIORI ET AL: "Side-lobe suppression techniques for a uniform circular array", 《 THE 7TH EUROPEAN RADAR CONFERENCE》 * |
| 王旭: "MIMO雷达发射方向图与波形设计方法研究", 《中国博士学位论文全文数据库信息科技辑》 * |
| 王旭等: "一种MIMO雷达多模式波形优化设计方法", 《电子与信息学报》 * |
| 王旭等: "针对主瓣波动的阵列赋形方向图设计", 《火控雷达技术》 * |
| 莫海生: "MIMO雷达发射方向图综合与DOA估计方法研究", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110501675A (zh) * | 2019-07-16 | 2019-11-26 | 北京工业大学 | 一种基于mimo雷达低旁瓣发射方向图设计方法 |
| CN115034075A (zh) * | 2022-06-22 | 2022-09-09 | 电子科技大学 | 非规则天线阵列矢量增益方向图的快速精准预测方法 |
| CN115034075B (zh) * | 2022-06-22 | 2023-11-24 | 电子科技大学 | 非规则天线阵列矢量增益方向图的快速精准预测方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107436970B (zh) | 2020-04-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108920767B (zh) | 余割平方和合成相位双约束波瓣阵列天线优化设计方法 | |
| Salam et al. | Smart underground antenna arrays: A soil moisture adaptive beamforming approach | |
| CN107038299B (zh) | 一种考虑互耦效应的变形阵列天线远场方向图补偿方法 | |
| CN103076596B (zh) | 基于先验信息的mimo雷达发射方向图设计方法 | |
| CN105426578A (zh) | 一种基于遗传算法的mimo-sar面阵天线阵元位置优化方法 | |
| CN104615854B (zh) | 一种基于稀疏约束的波束展宽和旁瓣抑制方法 | |
| CN104537185B (zh) | 基于量子蝙蝠搜索的天线阵稀疏构建与方向图综合方法 | |
| CN106772260A (zh) | 基于凸优化算法的雷达阵列和差波束方向图优化方法 | |
| CN102608580B (zh) | 数字阵列超低副瓣自适应数字波束形成方法 | |
| CN103308904B (zh) | 一种兼顾近场和远场性能的三维摄像声纳***换能器阵列的稀疏优化方法 | |
| CN105158756A (zh) | 集中式mimo雷达射频隐身时多目标跟踪波束指向方法 | |
| CN104166136A (zh) | 一种基于干扰子空间跟踪的高效自适应单脉冲测角方法 | |
| CN105445709A (zh) | 一种稀布阵列近场无源定位幅相误差校正方法 | |
| Li et al. | Improved GA and PSO culled hybrid algorithm for antenna array pattern synthesis | |
| CN104391183A (zh) | 近场测量天线远区场特性快速计算方法 | |
| CN110501675A (zh) | 一种基于mimo雷达低旁瓣发射方向图设计方法 | |
| CN101916904A (zh) | 移动卫星通信多子阵平板天线阵及其优化方法 | |
| CN103513225B (zh) | 一种基于空间增益的稀疏平面阵形优化方法 | |
| CN110427590A (zh) | 基于自适应概率学习的大型稀疏阵列天线高效综合方法 | |
| CN110535518A (zh) | 一种宽波束发射波束形成优化设计方法 | |
| CN107290732A (zh) | 一种量子大***的单基地mimo雷达测向方法 | |
| CN107135023A (zh) | 用于毫米波通信***的三维训练码书设计方法及波束对准方法 | |
| CN109639329A (zh) | 唯相位加权波束快速赋形方法 | |
| CN107436970A (zh) | 基于序列二次规划的三坐标雷达赋形方向图设计方法 | |
| CN109271735B (zh) | 基于量子启发式引力搜索算法的阵列方向图综合方法 |
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 |