[发明专利]基于机电耦合与最小二乘法的变形阵列天线电性能补偿方法

摘要

本发明公开了一种基于机电耦合与最小二乘法的变形阵列天线电性能补偿方法，包括1)确定阵列天线的结构参数和电磁工作参数；2)通过有限元分析得到天线阵面的结构变形量；3)计算理想情况和变形情况阵列天线的电性能；4)结合最小二乘法，计算得到补偿激励电流；5)判断补偿后的阵列天线电性能是否满足指标要求。本发明基于阵列天线的机电耦合模型，结合最小二乘法，对比理想和变形阵列天线的远区电场分布来求解激励电流的补偿量，然后调整电流实现对变形阵列天线电性能的补偿。该方法可在不增加天线结构重量和天线结构复杂度情况下，明显改善由环境载荷引起的天线电性能恶化，确保在各种服役环境下阵列天线都能正常工作。

主权项

基于机电耦合与最小二乘法的变形阵列天线电性能补偿方法，其特征在于，包括如下过程：(1)根据平面阵列天线的基本结构，确定平面阵列天线的几何模型参数、材料属性和电磁工作参数；(2)根据平面阵列天线的几何模型参数及材料属性，建立平面阵列天线结构有限元模型；根据平面阵列天线的安装形式确定天线有限元模型的约束位置及约束方式；对天线结构有限元模型施加随机振动加速度功率谱，计算平面阵列天线随机振动变形量，分别提取平面阵列天线的各个辐射单元中心节点在x,y,z方向上的位置偏移量(Δx,Δy,Δz)；(3)根据平面阵列天线的电磁工作参数及各个辐射单元中心节点的位置偏移量，利用阵列天线的机电耦合模型，分别计算理想情况和变形情况下的阵列天线的电性能；其中，机电耦合模型计算阵列天线电性能，按如下过程进行：(3a)设阵列天线中共有N个辐射单元，第n号辐射单元的位置矢量为观察点P相对于坐标系O‑xyz所在的方向(θ,φ)上的方位矢量为其中：cosαx=sinθcosφcosαy=sinθsinφcosαz=cosθ]]>式中，xn、yn、zn分别为位置矢量在x、y、z方向的投影；分别为x、y、z方向的单位矢量；αx、αy、αz分别为位置矢量与x、y、z方向的夹角；(3b)当天线阵面产生变形时，第n号辐射单元产生的位置误差为则变形后第n号辐射单元的位置矢量为：r→n′=(r→n+Δr→n)=(xn+Δxn)i^+(yn+Δyn)j^+(zn+Δzn)k^;]]>式中，Δxn、Δyn、Δzn分别为位置误差在x、y、z方向的投影；(3c)根据阵列天线电磁叠加原理，理想情况下阵列天线的方向图函数为：结合(3b)得到变形阵列天线方向图函数为：式中，为阵列天线的波常数，λ为阵列天线工作波长，fn(θ,φ)为辐射单元方向图函数，An是天线单元激励电流幅度，是天线单元激励电流相位；根据阵列天线方向图函数绘制的阵列天线方向图，计算得到阵列天线的增益、副瓣电平及波束指向等电性能参数指标；(4)由平面阵列天线机电耦合模型得到的理想阵列天线和变形阵列天线的方向图函数，结合最小二乘法，求解理想和变形天线方向图之间差值平方最小的方程，得到激励电流的补偿量；将补偿激励电流带入变形阵列天线的方向图函数表达式中，得到补偿后的天线方向图；其中，结合最小二乘法求解激励电流的补偿量按如下过程进行：(4a)利用机电耦合模型建立变形阵列天线方向图函数表达式，假定补偿激励电流为未知量，将变形阵列天线方向图中的激励电流替换为补偿激励电流得到补偿后的阵列天线方向图；(4b)理想天线方向图函数为E0(θ,φ)，补偿后的天线方向图为Ec(θ,φ)，则理想情况和补偿后天线方向图差值平方最小的方程为：Min∫∫|Ec(θ,φ)‑E0(θ,φ)|2·w(θ,φ)dθdφ式中，w(θ,φ)为与方位俯仰角(θ,φ)有关的权因子；(4c)分别将理想和补偿后的天线方向图在M个方位俯仰角上进行离散，理想天线方向图离散后的表达式为：补偿后天线方向图离散后的表达式为：式中，(θ1,φ1),(θ2,φ2)…(θM,φM)为M个方位俯仰角；为M个方位俯仰角分别对应的方位矢量；E0(θ1,φ1),E0(θ2,φ2)…E0(θM,φM)为M个方位俯仰角分别对应的理想天线方向图函数；Ec(θ1,φ1),Ec(θ2,φ2)…Ec(θM,φM)为M个方位俯仰角分别对应的补偿后的天线方向图函数；f(θ1,φ1),f(θ2,φ2)…f(θM,φM)为M个方位俯仰角分别对应的辐射单元的方向图函数；为阵列天线N个辐射单元分别对应的位置矢量；为阵列天线N个辐射单元分别对应产生的位置误差；A1,A2…AN为理想情况下阵列天线N个单元分别对应的激励电流幅度；为理想情况下阵列天线N个单元分别对应的激励电流相位；A′1,A′2…A′N为补偿后阵列天线N个单元分别对应的激励电流幅度，为补偿后阵列天线N个天线单元分别对应的激励电流相位；(4d)理想和补偿后天线方向图函数差值平方的方程写成离散形式：Min∫∫|Ec(θ,φ)-E0(θ,φ)|2·w(θ,φ)dθdφ=MinΣm=1M|Ec(θm,φm)-E0(θm,φm)|2w(θm,φm)=Min(W·|McI′-M0I0|2)]]>式中，Ec(θm,φm)为第m个方位俯仰角(θm,φm)对应的补偿后的天线方向图；E0(θm,φm)为第m个方位俯仰角(θm,φm)对应的理想阵列天线的方向图；w(θm,φm)为第m个方位俯仰角(θm,φm)对应的权因子；Mc、M0为M×N维的矩阵，I′、I0为N×1维的列向量，W为M×M维的对角阵，其各自具体的表达式如下：其中，w(θ1,φ1),w(θ2,φ2)…w(θM,φM)为M个方位俯仰角(θ1,φ1),(θ2,φ2)…(θM,φM)分别对应的权因子；M0为理想情况下天线方向图权系数矩阵；和Mc为存在位置误差情况下天线方向图权系数矩阵；(4e)求解理想和补偿后天线方向图差值平方的方程得到补偿后的激励电流为：I′=(McTWMc)-1(McTWM0)I0]]>式中，I0表示理想情况下阵列天线的激励电流，是N×1维的列向量；I′表示补偿之后阵列天线的激励电流，是N×1维的列向量；W是权系数矩阵，为M×M维的对角阵；M0和Mc分别为理想和存在位置误差情况下天线方向图权系数矩阵，均为M×N维的矩阵；将该补偿激励电流带入变形阵列天线的方向图函数中，绘制补偿后的变形阵列天线方向图，计算补偿后的变形阵列天线电性能，得到包括增益、副瓣电平以及波束指向等天线电性能参数；(5)判断补偿后的天线电性能是否满足指标要求，如果满足要求，则该方法能够实现对变形阵列天线电性能的补偿；否则，修改阵列天线的结构参数，并重复步骤(1)到步骤(5)，直至满足要求。