基于qiwo_ift组合算法的阵列天线唯相位变换波束赋形方法
【技术领域】
[0001] 本发明涉及一种基于阵列天线唯相位变换的波束赋形方法。涉及在电流幅度不变 仅相位变化的情况下实现对远场方向图变化的技术。
【背景技术】
[0002] 阵列天线唯相位变换波束赋形技术是一种仅通过相位加权来改变方向图赋形的 方法。不同的辐射方向图是为了实现不同的作用。在很多实际应用中，如空中交通管制雷 达、航天器、卫星以及无线通信等，需要产生不同的辐射方向图。唯相位变换波束赋形方法 使得天线系统仅使用一个天线阵即可完成不同的任务，这就需要通过一些算法寻求一组合 理的激励权值来控制天线波束的指向、宽度、形状、副瓣电平、零陷点的位置。采用只改变馈 电相位分布的仅相位加权方法可使其不改变原有功率分配馈电网络和不增加新设备的情 况下，利用计算机控制移相器值的改变实现波束赋形，是非常经济的可行方法，由于相位控 制比功率控制更加容易实现且拥有更小的功率损耗，因此可以通过唯相位变换技术减少天 线增益损失。因此，研究阵列天线的唯相位变换技术是非常必要的。
[0003] 常规的方法是以单元的激励相位作为优化变量，通过优化算法来完成，如Davis 将最小二乘法应用于唯相加权的阵列波束赋形；Haupt用遗传算法实现唯相加权的自适应 零点生成；Khzmalyan采用单坐标最小化的方法实现唯相加权波束赋形和自适应零点生 成，常规的方法当单元的数量较多时，优化变量也就很多，这样导致优化过程非常耗时，优 化的结果无法逼近理想方向图。
[0004] 本发明利用量化入侵杂草算法（QIWO)快速获得理想方向图的幅相值作为经典傅 里叶（IFFT)迭代算法的初值，结合了量化入侵杂草算法良好的平衡全局收敛和局部收敛 能力以及迭代傅里叶算法成熟和运算速度快的优点，回避了初值选择不当后迭代傅里叶算 法收敛速度慢或无法收敛的缺点，最终实现在功率保持一定的情况下，仅通过相位加权完 成方向图赋形的变化。
【发明内容】
[0005] 实现本发明的目的在于提供一种线性阵列唯相位变化方向图综合方法。
[0006] 实现本发明目的的技术解决方案为：将线性阵列的每一单元电流幅度与相位幅度 做等间隔分段处理，取分段后每段的中间值作为量化值，算法运行中计算可能的初值，中间 迭代值以及满足唯相位变换方向图的终值均量化到选择分段后每段的中间值；选定智能算 法中的初始种群数，最大种群数标准差，标准差最终值，非线性调节因子，缩放影子，适应度 函数，循环次数等固定参数和计算公式，将随机选取的量化值作为每个阵列单元幅度相位 的初始迭代值，取得满足条件的单一方向图综合的幅相值。每个单元的幅相值作为迭代傅 里叶变换的初值，叶逆变换（IFFT)得出阵列方向图，然后与期望方向图比较，得出新的方 向图，再应用二维快速傅里叶变换（FFT)反算出新的单元的激励。
[0007] 本发明关于一种基于QIW0_IFT组合算法的阵列天线唯相位变换波束赋形方法效 果体现在下述几个方面：
[0008] (1)量化的相位值，便于应用在实际工程。
[0009] (2)量化入侵杂草算法将每个单元的电流和相位进行量化处理，将无限可能的值 简化为有限个可能的值，通过缩小搜索范围来获得较强的局部搜索能力以增强算法的整体 收敛性能。
[0010] (3)将量化入侵杂草算法得出的单元幅相值作为迭代傅里叶算法的初值，提高经 典迭代傅里叶算法收敛速度的同时保证了算法最终的收敛。
【附图说明】
[0011] 图1 :共电流幅度仅相位加权下平顶波束赋形效果图。
[0012] 图2 :共电流幅度仅相位加权下余割平方赋形效果图。
[0013] 图3 :计算得出的共电流幅度值。
[0014] 图4 :计算得出的各方向图对应的相位值。
【具体实施方式】
[0015] 本发明的技术方案具体实施过程如下：
[0016] 步骤一：将每个单元的幅相参数量化，确定幂数N，对线性阵列的每一个单元的电 流幅度和相位幅值分别按照上下限进行等间隔分段，分段数量为2 N，取每个分段的中间值 作为量化值，并对每一个量化后电流幅度和相位幅度分别度按1~2N进行编号；假设QITO 算法的单元数为M，则初始种群数初值设定为2M，最大种群数设定10M，可产生的种子区间 为[1，Μ/2]，解搜索空间为[1，2 N] ;QIW0算法中的适应度函数选择为：
[0017] 0. 8 X IMSLL-SLVL|+0. 2|LL_MAX-NLVL
[0018] 其中：MSLL为方向图旁瓣最高值；SLVL为理想方向图旁瓣最高值；NLVL为所有指 定角度的模值，LL_MAX为理想方向图的波纹的上下限；
[0019] 步骤二：将随机产生的每个单元的幅相量化值带入QIWO算法，迭代后进而得出每 一个理想方向图对应的单元幅相值，流程包括；
[0020] (a)设QIWO算法的量化初值为随机产生的整数值，对于单元数为M，则初始解为X =[X1, X2,…X2J ;
[0021] (b)算法中杂草生长的过程，杂草每次成长的最新值的前M个值即为每个单元的 幅度，后M个值即为每个单元的相位值，每个杂草个体可产生的种子数为：
[0022]
[0023] 其中$_和F_为此次进化中的最大最小适应度值，SniaJP S_为可产生的最大种 子数和最小种子数，f (Xi)为第i个杂草个体的适应度值，Fflcra为函数向下取整；
[0024] (C)让每个单元对应的幅度值，即种子个体在D维搜索空间，以正态分布随机分散 在父代杂草附近，即第i个单元的第S个幅度或相位种子的位置为；
[0025] Xlis=XJNCoiO2)
[0026] ⑷当在迭代过程中，当种群数超过最大种群数的限制时，则按照方向图的适应度 计算结果排序，对排序后的个体按适应度值由大到小依次选出IOM个，作为该次迭代最终 保留下来的种群；
[0027] 经过上述流程，满足最大迭代次数和达到理想方向图要求两个条件中的任意一个 后，完成QIWO算法的计算过程，最终每个方向图的单元幅相值作为迭代傅里叶变换运算的 初值；
[0028] 步骤三：每个单元的幅相值作为迭代傅里叶变换的初值，傅里叶逆变换得出阵列 方向图，然后与期望方向图比较，得出新的方向图；
[0029] 步骤四：将新的方向图应用FFT反算出新的单元激励；
[0030] (a)在IFFT得出阵列方向图的过程中，每一个单元的模值为上次求出多个方向图 对应模值的均值，相位幅值则为上次求出方向图各自对应的相位值，并用〇补齐余值；
[0031] (b)在FFT反算各单元激励的过程中，对于方向图与理想方向图的模值不一致性， 采用保留原方向图的相位值，用理想方向图的模值直接覆盖的方法，迭代后取前M值作为M 个单元的幅相值；
[0032] (C)在迭代傅里叶算法的过程中当迭代次数超过预定值后自动停止迭代。
[0033] 实例1 :线阵天线由29个单元组成，单元间距0. 5 λ，约束条件为平顶波束主瓣为 sin θ I彡0.41，波纹为0. 1，副瓣为-25dB，余割平方赋形左副瓣范围-1彡sin Θ彡-〇. 15， 右副瓣范围为0. 7彡sin Θ彡1，主瓣范围-0.034彡sin Θ彡0.6,波纹为0. 1，副瓣 为-25dB。表1中为各方向图同功率不同相位值。
[0034] 表 1
[0035]
【主权项】
1. 一种基于QIWO_IFT组合算法的阵列天线唯相位变换波束赋形方法，其特征在于包 括以下步骤： 步骤一：将每个单元的幅相参数量化，确定幂数N，对线性阵列的每一个单元的电流幅 度和相位幅值分别按照上下限进行等间隔分段，分段数量为2N，取每个分段的中间值作为 量化值，并对每一个量化后电流幅度和相位幅度分别度按1~2N进行编号；假设QIWO算法 的单元数为M，则初始种群数初值设定为2M，最大种群数设定10M，可产生的种子区间为[1， M/2]，解搜索空间为[1，2N]，QIWO算法中的适应度函数选择为： 0?8XIMSLL-SLVLI+0? 2ILL_MAX-NLVLI， 其中：MSLL为方向图旁瓣最高值；SLVL为理想方向图旁瓣最高值；NLVL为所有指定角 度的模值，LL_MAX为理想方向图的波纹的上下限； 步骤二：将随机产生的每个单元的幅相量化值带入QIWO算法，迭代后进而得出每一个 理想方向图对应的单元幅相值，流程包括； (a) 设QIWO算法的量化初值为随机产生的整数值，对于单元数为M，则初始解为X= [Xi，乂2，? ? ?XgM]; (b) 算法中杂草生长的过程，杂草每次成长的最新值的前M个值即为每个单元的幅度， 后M个值即为每个单元的相位值，每个杂草个体可产生的种子数为：其中$_和F_为此次进化中的最大最小适应度值，SniaJPS_为可产生的最大种子数 和最小种子数，f(Xi)为第i个杂草个体的适应度值，Fflcra为函数向下取整； (c) 让每个单元对应的幅度值，即种子个体在D维搜索空间，以正态分布随机分散在父 代杂草附近，即第i个单元的第S个幅度或相位种子的位置为； Xlis=X^(0, 〇2) (d) 当在迭代过程中，当种群数超过最大种群数的限制时，则按照方向图的适应度计算 结果排序，对排序后的个体按适应度值由大到小依次选出IOM个，作为该次迭代最终保留 下来的种群； 经过上述流程，满足最大迭代次数和达到理想方向图要求两个条件中的任意一个后， 完成QIWO算法的计算过程，最终每个方向图的单元幅相值作为迭代傅里叶变换运算的初 值； 步骤三：每个单元的幅相值作为迭代傅里叶变换的初值，傅里叶逆变换得出阵列方向 图，然后与期望方向图比较，得出新的方向图； 步骤四：将新的方向图应用FFT反算出新的单元激励： (a) 在IFFT得出阵列方向图的过程中，每一个单元的模值为上次求出多个方向图对应 模值的均值，相位幅值则为上次求出方向图各自对应的相位值，并用〇补齐余值； (b) 在FFT反算各单元激励的过程中，对于方向图与理想方向图的模值不一致性，采用 保留原方向图的相位值，用理想方向图的模值直接覆盖的方法，迭代后取前M值作为M个单 元的幅相值； (c) 在迭代傅里叶算法的过程中当迭代次数超过预定值后自动停止迭代。
【专利摘要】本发明涉及一种基于QIWO_IFT组合算法的阵列天线唯相位变换波束赋形方法，属于无线通信、信号处理技术领域。本发明涉及一种阵列天线唯相位变换波束赋形方法，即在功率保持不变的情况下仅通过相位加权来实现不同辐射要求的波束，该发明包括以下部分：使用量化入侵杂草算法得出每一个理想方向图对应的单元幅相值作为迭代傅里叶变换的初值，应用离散傅里叶逆变换(IFFT)得出数组方向图，然后与期望方向图比较，得出新的方向图，再应用快速傅里叶变换(FFT)反算出新的单元激励。本发明结合量化入侵杂草算法方向图赋形的快速性、有效性以及迭代傅里叶算法能最终收敛和收敛速度快的特点，提高仅相位变化各方向图的空域内波束覆盖的能力，因而具有较好的工程实用性。
【IPC分类】H01Q25/00
【公开号】CN105048112
【申请号】CN201510362210
【发明人】祁峥东, 高国明, 袁洪, 简玲
【申请人】中国船舶重工集团公司第七二四研究所
【公开日】2015年11月11日
【申请日】2015年6月26日