CN109408956B - 基于单元节点位移的反射面天线平均功率方向图建模方法 - Google Patents

基于单元节点位移的反射面天线平均功率方向图建模方法 Download PDF

Info

Publication number
CN109408956B
CN109408956B CN201811239147.0A CN201811239147A CN109408956B CN 109408956 B CN109408956 B CN 109408956B CN 201811239147 A CN201811239147 A CN 201811239147A CN 109408956 B CN109408956 B CN 109408956B
Authority
CN
China
Prior art keywords
unit
displacement
node
antenna
electric field
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.)
Active
Application number
CN201811239147.0A
Other languages
English (en)
Other versions
CN109408956A (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 CN201811239147.0A priority Critical patent/CN109408956B/zh
Publication of CN109408956A publication Critical patent/CN109408956A/zh
Application granted granted Critical
Publication of CN109408956B publication Critical patent/CN109408956B/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
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation

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)
  • Aerials With Secondary Devices (AREA)

Abstract

本发明公开了一种基于单元节点位移的反射面天线平均功率方向图建模方法，包括：输入反射面天线几何参数与电参数；计算理想天线远区辐射电场；计算理想天线远区辐射功率；进行反射面三角形网格划分；计算单元节点位移与单元中心点位移的节点转换矩阵；组集单元节点位移与单元中心点位移的总体转换矩阵；计算电场相对单元节点位移的单元一阶导数、二阶Hessian阵；组集电场相对单元节点位移的总体一阶导数、二阶Hessian阵；计算电场相对单元中心点位移的总体一阶导数、二阶Hessian阵；输入表面随机误差均方根值；计算天线远区辐射功率平均值；判断电性能是否满足要求；输出辐射功率方向图；更新表面随机误差均方根值。本发明可指导反射面天线面板加工与制造。

Description

基于单元节点位移的反射面天线平均功率方向图建模方法
技术领域
本发明属于雷达天线技术领域，具体涉及雷达天线领域中的一种基于单元节点位移的反射面天线平均功率方向图建模方法。
背景技术
反射面天线广泛应用于射电天文、雷达、通信、探测等领域。在外部工作环境下，反射面天线很容易受到自重、温度、风荷等外部载荷影响；在加工与制造过程，反射面天线面板也受到加工制造与安装误差影响；针对外部载荷引入的系统误差与加工制造引入的随机误差一直是反射面天线设计领域研究的热点。
Y.Rahmat-Samii在文献“An efficient computational method forcharacterizing the effects of random surface errors on the average powerpattern of reflectors”(IEEE Trans.Antennas and Propagation,1983年第31卷第1期，92-98)公开了一种基于概率方法分析天线表面随机误差对电性能影响的分析方法。王猛、段宝岩、王伟等人在文献“反射面天线表面误差对平均功率方向图的影响”(西安电子科技大学学报，2014年第41卷第6期188-194)中提出了一种表面随机误差与系统误差同时存在下的平均功率方向图计算方法。由于现有方法均是将表面误差以相位误差的形式引入到电性能计算中，导致在公式推导中的繁琐与耗时，难以形成快速分析的目的。因此，本发明为了适应提高分析效率的需要，创新性地将天线辐射电场描述为三角形单元节点位移的二阶近似公式，提出一种基于单元节点位移的反射面天线平均功率方向图建模方法。
发明内容
本发明的目的是克服上述现有技术的不足，提供一种基于单元节点位移的反射面天线平均功率方向图建模方法。该方法基于单元节点位移二阶近似公式，利用转换矩阵，提出了分析随机误差对反射面天线平均功率方向图影响的建模方法，可指导反射面天线面板加工与制造。
本发明的技术方案是：基于单元节点位移的反射面天线平均功率方向图建模方法，其特征在于，包括如下步骤：
(1)输入反射面天线几何参数与电参数
输入用户提供的反射面天线几何参数与电参数；其中几何参数包括半径、焦距；电参数包括工作波长、自由空间波常数、口径场幅度分布函数、锥销电平、口径场形状指数以及包括天线增益、波瓣宽度、副瓣电平、指向精度在内的电性能要求；
(2)计算理想天线远区辐射电场
根据反射面天线几何参数与电参数，利用几何光学法计算理想天线远区辐射电场；
(3)计算理想天线远区辐射功率
根据理想天线远区辐射电场信息，通过下式计算理想天线远区辐射功率
Figure BDA0001838848450000021
其中，G0表示理想天线远区辐射功率，E0表示理想天线远区辐射电场，下标0代表理想天线，上标*表示取共轭运算；
(4)进行反射面三角形网格划分
根据用户提供的天线几何参数与工作波长，将反射面口径面划分为一系列等边三角形单元，等边三角形单元边长满足以下关系式
Figure BDA0001838848450000031
其中，λ为工作波长，l为口径面等边三角形边长；
(5)计算单元节点位移与单元中心点位移的节点转换矩阵
根据反射面三角形网格划分后三角形单元信息，通过下式计算单元节点位移与单元中心点位移的节点转换矩阵
Figure BDA0001838848450000032
其中，Ki表示第i个单元节点位移与其环绕的单元中心点位移之间的节点转换矩阵，mi为反射面三角形网格划分后环绕第i个单元节点的三角形单元数目，Ki的维数为mi行1列，上标T表示转置运算；
(6)组集单元节点位移与单元中心点位移的总体转换矩阵
根据反射面三角形网格划分后三角形单元信息，将单元节点位移与单元中心点位移的节点转换矩阵，按照三角形单元信息进行组集，通过下式组集单元节点位移与单元中心点位移的总体转换矩阵
Figure BDA0001838848450000033
其中，K表示单元节点位移与单元中心点位移的总体转换矩阵，Ki为第i个单元节点位移与单元中心点位移的节点转换矩阵，N为单元节点总数；
(7)计算电场相对单元节点位移的单元一阶导数、二阶Hessian阵；
(8)组集电场相对单元节点位移的总体一阶导数、二阶Hessian阵；
(9)计算电场相对单元中心点位移的总体一阶导数、二阶Hessian阵；
(10)输入表面随机误差均方根值
根据反射面天线面板加工制造误差，输入天线表面随机误差均方根值；
(11)计算天线远区辐射功率平均值；
(12)判断电性能是否满足要求
判断天线远区辐射功率平均值是否满足天线增益、波瓣宽度、副瓣电平、指向精度在内的电性能要求，如果满足要求则转至步骤(13)，否则转至步骤(14)；
(13)输出辐射功率方向图
当天线远区辐射功率平均值满足天线增益、波瓣宽度、副瓣电平、指向精度在内的电性能要求时，输出辐射功率方向图；
(14)更新表面随机误差均方根值
当天线远区辐射功率平均值不满足天线增益、波瓣宽度、副瓣电平、指向精度在内的电性能要求时，更新表面随机误差均方根值，转至步骤(10)。
所述步骤(2)中的几何光学法是一种基于口径场分布的高频近似方法，计算公式如下：
Figure BDA0001838848450000041
其中，E0表示理想天线远区辐射电场，∫表示积分运算，ρ表示反射面上点在口径面内的极坐标分量，Q(ρ)表示用户输入的反射面上点对应的口径场幅度分布函数，exp表示自然对数的指数运算，j为虚数单位，k为自由空间波常数，
Figure BDA0001838848450000051
为反射面上点在口径面内的位置矢量，
Figure BDA0001838848450000052
为远区观察点单位位置矢量，s表示天线口径面，ds表示在口径面内进行积分运算。
所述步骤(7)计算电场相对单元节点位移的单元一阶导数、二阶Hessian阵是根据用户提供的天线几何参数与电参数，结合三角形单元，通过下式电场相对单元节点位移的单元一阶导数
Gi＝[Gi,1,Gi,2,Gi,3]T
Figure BDA0001838848450000053
其中，Gi为电场相对单元节点位移的单元一阶导数，下标i表示第i个单元，Gi,l表示电场相对第i个单元上第l个节点的单元一阶导数分量，l＝1,2,3,∫表示积分运算，Q(ρ)表示用户输入的口径场幅度分布函数，ρ表示单元内部点在口径面内的极坐标分量，exp表示自然对数的指数运算，j为虚数单位，k为自由空间波常数，
Figure BDA0001838848450000054
为单元内部点在口径面内的位置矢量，
Figure BDA0001838848450000055
为远区观察点单位位置矢量，γl表示单元内部相对于第l个节点的单元形函数，ξ表示单元内部点在在馈源坐标系下的角度分量，si表示第i个单元在口径面内的投影面积，dsi表示在第i个单元投影面积内进行积分运算；
根据用户提供的天线几何参数与电参数，结合三角形单元，通过下式电场相对单元节点位移的单元二阶Hessian阵
Figure BDA0001838848450000056
Figure BDA0001838848450000057
其中，Hi为电场相对单元节点位移的单元二阶Hessian阵，下标i表示第i个单元，Hi,uv表示电场相对第i个单元上第u、v个节点的单元二阶Hessian阵元素分量，u,v＝1,2,3,∫表示积分运算，Q(ρ)表示用户输入的口径场幅度分布函数，ρ表示单元内部点在口径面内的极坐标分量，exp表示自然对数的指数运算，j为虚数单位，k为自由空间波常数，
Figure BDA0001838848450000061
为单元内部点在口径面内的位置矢量，
Figure BDA0001838848450000062
为远区观察点单位位置矢量，γu、γv分别表示单元内部相对于第u、v个节点的单元形函数，ξ表示单元内部点在在馈源坐标系下的角度分量，si表示第i个单元在口径面内的投影面积，dsi表示在第i个单元投影面积内进行积分运算。
所述步骤(8)组集电场相对单元节点位移的总体一阶导数、二阶Hessian阵是根据反射面三角形网格划分后三角形单元信息，通过下式组集电场相对单元节点位移的总体一阶导数、二阶Hessian阵
Figure BDA0001838848450000063
Figure BDA0001838848450000064
其中，GE表示电场相对单元节点位移的总体一阶导数，Gi为电场相对于单元节点位移的单元一阶导数，M表示三角形单元总数，A代表组集运算，HE表示电场相对单元节点位移的总体二阶Hessian阵，Hi为电场相对于单元节点位移的单元二阶Hessian阵。
所述步骤(9)计算电场相对单元中心点位移的总体一阶导数、二阶Hessian阵是根据单元节点位移与单元中心点位移的总体转换矩阵，结合电场相对单元节点位移的总体一阶导数、二阶Hessian阵，通过下式计算电场相对于单元中心点位移的总体一阶导数、二阶Hessian阵
G＝KGE
H＝KHEKT
其中，G为电场相对于单元中心点位移的总体一阶导数，H为电场相对于单元中心点位移的二阶Hessian阵，GE表示电场相对单元节点位移的总体一阶导数，HE表示电场相对单元节点位移的总体二阶Hessian阵，K表示单元节点位移与单元中心点位移的总体转换矩阵，上标T表示转置运算。
所述步骤(11)计算天线远区辐射功率平均值是根据理想天线远区辐射电场、远区辐射功率、电场相对单元中心点位移的单元一阶、二阶系数与表面随机误差均方根值，通过下式计算天线远区辐射功率平均值
Figure BDA0001838848450000071
其中，
Figure BDA0001838848450000072
表示天线远区辐射功率平均值，G0表示理想天线远区辐射功率，E0表示理想天线远区辐射电场，下标0代表理想天线，上标*表示取共轭运算，M为反射面上三角形单元总数，Hii、Hij、Hjj分别表示电场相对于单元中心点位移的总体二阶Hessian阵的第i行第i列分量、第i行第j列分量，第j行第j列分量，Gi电场相对于单元中心点位移的总体一阶导数第i个分量，σ表示用户输入的天线表面随机误差均方根值。
本发明的有益效果：本发明首先输入天线几何参数与电参数信息，计算理想天线远区辐射电场与远区辐射功率；其次，进行反射面三角形网格划分，分别计算单元节点位移与单元中心点位移的节点转换矩阵、总体转换矩阵；然后，计算电场相对单元节点位移的单元一阶导数、二阶Hessian阵、组集电场相对单元节点位移的总体一阶导数、二阶Hessian阵；再次，计算电场相对单元中心点位移的总体一阶导数、二阶Hessian阵；然后，根据用户输入的表面随机误差均方根值，计算天线远区辐射功率平均值；最后，判断辐射功率平均值是否满足电性能要求，并输出辐射功率方向图，以此指导反射面面板加工与制造。
与现有技术相比，本发明具有以下优点：
1.本发明基于单元节点位移二阶近似公式获得随机误差影响下反射面天线远区辐射功率平均值，在保证计算精度的前提下，避免了繁琐的公式推导，提高分析效率；
2.本发明从概率的角度出发，获得反射面天线远区辐射功率平均值，采用二阶近似公式保证了分析的准确性。
以下将结合附图对本发明做进一步详细说明。
附图说明
图1为本发明的流程图；
图2为传统方法与本发明方法在增益平均值的比较曲线；
图3为传统方法与本发明方法在平均功率方向图的比较曲线。
具体实施方式
下面结合附图对本发明具体实施方式作进一步的详细描述：
实施例1
基于单元节点位移的反射面天线平均功率方向图建模方法，包括如下步骤：
步骤1，输入反射面天线几何参数与电参数
输入用户提供的反射面天线几何参数与电参数；其中几何参数包括半径、焦距；电参数包括工作波长、自由空间波常数、口径场幅度分布函数、锥销电平、口径场形状指数以及包括天线增益、波瓣宽度、副瓣电平、指向精度在内的电性能要求；
步骤2，计算理想天线远区辐射电场
根据反射面天线几何参数与电参数，利用几何光学法计算理想天线远区辐射电场；
步骤3，计算理想天线远区辐射功率
根据理想天线远区辐射电场信息，通过下式计算理想天线远区辐射功率
Figure BDA0001838848450000091
其中，G0表示理想天线远区辐射功率，E0表示理想天线远区辐射电场，下标0代表理想天线，上标*表示取共轭运算；
步骤4，进行反射面三角形网格划分
根据用户提供的天线几何参数与工作波长，将反射面口径面划分为一系列等边三角形单元，等边三角形单元边长满足以下关系式
Figure BDA0001838848450000092
其中，λ为工作波长，l为口径面等边三角形边长；
步骤5，计算单元节点位移与单元中心点位移的节点转换矩阵
根据反射面三角形网格划分后三角形单元信息，通过下式计算单元节点位移与单元中心点位移的节点转换矩阵
Figure BDA0001838848450000093
其中，Ki表示第i个单元节点位移与其环绕的单元中心点位移之间的节点转换矩阵，mi为反射面三角形网格划分后环绕第i个单元节点的三角形单元数目，Ki的维数为mi行1列，上标T表示转置运算；
步骤6，组集单元节点位移与单元中心点位移的总体转换矩阵
根据反射面三角形网格划分后三角形单元信息，将单元节点位移与单元中心点位移的节点转换矩阵，按照三角形单元信息进行组集，通过下式组集单元节点位移与单元中心点位移的总体转换矩阵
Figure BDA0001838848450000101
其中，K表示单元节点位移与单元中心点位移的总体转换矩阵，Ki为第i个单元节点位移与单元中心点位移的节点转换矩阵，N为单元节点总数；
步骤7，计算电场相对单元节点位移的单元一阶导数、二阶Hessian阵；
步骤8，组集电场相对单元节点位移的总体一阶导数、二阶Hessian阵；
步骤9，计算电场相对单元中心点位移的总体一阶导数、二阶Hessian阵；
步骤10，输入表面随机误差均方根值
根据反射面天线面板加工制造误差，输入天线表面随机误差均方根值；
步骤11，计算天线远区辐射功率平均值；
步骤12，判断电性能是否满足要求
判断天线远区辐射功率平均值是否满足天线增益、波瓣宽度、副瓣电平、指向精度在内的电性能要求，如果满足要求则转至步骤13，否则转至步骤14；
步骤13，输出辐射功率方向图
当天线远区辐射功率平均值满足天线增益、波瓣宽度、副瓣电平、指向精度在内的电性能要求时，输出辐射功率方向图；
步骤14，更新表面随机误差均方根值
当天线远区辐射功率平均值不满足天线增益、波瓣宽度、副瓣电平、指向精度在内的电性能要求时，更新表面随机误差均方根值，转至步骤10。
实施例2
如图1所示，本发明提供了一种基于单元节点位移的反射面天线平均功率方向图建模方法，包括如下步骤：
步骤1，输入反射面天线几何参数与电参数
输入用户提供的反射面天线几何参数与电参数；其中几何参数包括半径、焦距；电参数包括工作波长、自由空间波常数、口径场幅度分布函数、锥销电平、口径场形状指数以及包括天线增益、波瓣宽度、副瓣电平、指向精度在内的电性能要求；
步骤2，计算理想天线远区辐射电场
根据反射面天线几何参数与电参数，通过下式计算理想天线远区辐射电场；
Figure BDA0001838848450000111
其中，E0表示理想天线远区辐射电场，∫表示积分运算，ρ表示反射面上点在口径面内的极坐标分量，Q(ρ)表示用户输入的反射面上点对应的口径场幅度分布函数，exp表示自然对数的指数运算，j为虚数单位，k为自由空间波常数，
Figure BDA0001838848450000112
为反射面上点在口径面内的位置矢量，
Figure BDA0001838848450000113
为远区观察点单位位置矢量，s表示天线口径面，ds表示在口径面内进行积分运算；
步骤3，计算理想天线远区辐射功率
根据理想天线远区辐射电场信息，通过下式计算理想天线远区辐射功率
Figure BDA0001838848450000121
其中，G0表示理想天线远区辐射功率，E0表示理想天线远区辐射电场，下标0代表理想天线，上标*表示取共轭运算；
步骤4，进行反射面三角形网格划分
根据用户提供的天线几何参数与工作波长，将反射面口径面划分为一系列等边三角形单元，等边三角形单元边长满足以下关系式
Figure BDA0001838848450000122
其中，λ为工作波长，l为口径面等边三角形边长。
步骤5，计算单元节点位移与单元中心点位移的节点转换矩阵
根据反射面三角形网格划分后三角形单元信息，通过下式计算单元节点位移与单元中心点位移的节点转换矩阵
Figure BDA0001838848450000123
其中，Ki表示第i个单元节点位移与其环绕的单元中心点位移之间的节点转换矩阵，mi为反射面三角形网格划分后环绕第i个单元节点的三角形单元数目，Ki的维数为mi行1列，上标T表示转置运算；
步骤6，组集单元节点位移与单元中心点位移的总体转换矩阵
根据反射面三角形网格划分后三角形单元信息，将单元节点位移与单元中心点位移的节点转换矩阵，按照三角形单元信息进行组集，通过下式组集单元节点位移与单元中心点位移的总体转换矩阵
Figure BDA0001838848450000131
其中，K表示单元节点位移与单元中心点位移的总体转换矩阵，Ki为第i个单元节点位移与单元中心点位移的节点转换矩阵，N为单元节点总数；
步骤7，计算电场相对单元节点位移的单元一阶导数、二阶Hessian阵
根据用户提供的天线几何参数与电参数，结合三角形单元，通过下式电场相对单元节点位移的单元一阶导数
Gi＝[Gi,1,Gi,2,Gi,3]T
Figure BDA0001838848450000132
其中，Gi为电场相对单元节点位移的单元一阶导数，下标i表示第i个单元，Gi,l表示电场相对第i个单元上第l个节点的单元一阶导数分量，l＝1,2,3,∫表示积分运算，Q(ρ)表示用户输入的口径场幅度分布函数，ρ表示单元内部点在口径面内的极坐标分量，exp表示自然对数的指数运算，j为虚数单位，k为自由空间波常数，
Figure BDA0001838848450000133
为单元内部点在口径面内的位置矢量，
Figure BDA0001838848450000134
为远区观察点单位位置矢量，γl表示单元内部相对于第l个节点的单元形函数，ξ表示单元内部点在在馈源坐标系下的角度分量，si表示第i个单元在口径面内的投影面积，dsi表示在第i个单元投影面积内进行积分运算；
根据用户提供的天线几何参数与电参数，结合三角形单元，通过下式电场相对单元节点位移的单元二阶Hessian阵
Figure BDA0001838848450000141
Figure BDA0001838848450000142
其中，Hi为电场相对单元节点位移的单元二阶Hessian阵，下标i表示第i个单元，Hi,uv表示电场相对第i个单元上第u、v个节点的单元二阶Hessian阵元素分量，u,v＝1,2,3,∫表示积分运算，Q(ρ)表示用户输入的口径场幅度分布函数，ρ表示单元内部点在口径面内的极坐标分量，exp表示自然对数的指数运算，j为虚数单位，k为自由空间波常数，
Figure BDA0001838848450000143
为单元内部点在口径面内的位置矢量，
Figure BDA0001838848450000144
为远区观察点单位位置矢量，γu、γv分别表示单元内部相对于第u、v个节点的单元形函数，ξ表示单元内部点在在馈源坐标系下的角度分量，si表示第i个单元在口径面内的投影面积，dsi表示在第i个单元投影面积内进行积分运算；
步骤8，组集电场相对单元节点位移的总体一阶导数、二阶Hessian阵
根据反射面三角形网格划分后三角形单元信息，通过下式组集电场相对单元节点位移的总体一阶导数、二阶Hessian阵
Figure BDA0001838848450000145
Figure BDA0001838848450000146
其中，GE表示电场相对单元节点位移的总体一阶导数，Gi为电场相对于单元节点位移的单元一阶导数，M表示三角形单元总数，A代表组集运算，HE表示电场相对单元节点位移的总体二阶Hessian阵，Hi为电场相对于单元节点位移的单元二阶Hessian阵；
步骤9，计算电场相对单元中心点位移的总体一阶导数、二阶Hessian阵
根据单元节点位移与单元中心点位移的总体转换矩阵，结合电场相对单元节点位移的总体一阶导数、二阶Hessian阵，通过下式计算电场相对于单元中心点位移的总体一阶导数、二阶Hessian阵
G＝KGE
H＝KHEKT
其中，G为电场相对于单元中心点位移的总体一阶导数，H为电场相对于单元中心点位移的二阶Hessian阵，GE表示电场相对单元节点位移的总体一阶导数，HE表示电场相对单元节点位移的总体二阶Hessian阵，K表示单元节点位移与单元中心点位移的总体转换矩阵，上标T表示转置运算；
步骤10，输入表面随机误差均方根值
根据反射面天线面板加工制造误差，输入天线表面随机误差均方根值；
步骤11，计算天线远区辐射功率平均值
根据理想天线远区辐射电场、远区辐射功率、电场相对单元中心点位移的单元一阶、二阶系数与表面随机误差均方根值，通过下式计算天线远区辐射功率平均值
Figure BDA0001838848450000151
其中，
Figure BDA0001838848450000152
表示天线远区辐射功率平均值，G0表示理想天线远区辐射功率，E0表示理想天线远区辐射电场，下标0代表理想天线，上标*表示取共轭运算，M为反射面上三角形单元总数，Hii、Hij、Hjj分别表示电场相对于单元中心点位移的总体二阶Hessian阵的第i行第i列分量、第i行第j列分量，第j行第j列分量，Gi电场相对于单元中心点位移的总体一阶导数第i个分量，σ表示用户输入的天线表面随机误差均方根值；
步骤12，判断电性能是否满足要求
判断天线远区辐射功率平均值是否满足天线增益、波瓣宽度、副瓣电平、指向精度在内的电性能要求，如果满足要求则转至步骤13，否则转至步骤14；
步骤13，输出辐射功率方向图
当天线远区辐射功率平均值满足天线增益、波瓣宽度、副瓣电平、指向精度在内的电性能要求时，输出辐射功率方向图；
步骤14，更新表面随机误差均方根值
当天线远区辐射功率平均值不满足天线增益、波瓣宽度、副瓣电平、指向精度在内的电性能要求时，更新表面随机误差均方根值，转至步骤10。
本发明的优点可通过以下仿真实验进一步说明：
1.仿真条件：
反射面天线口径100λ，焦距100λ，λ为工作波长，口径场幅度分布函数为
Figure BDA0001838848450000161
其中，Q为口径场幅度分布函数，为口径面内节点极坐标分量，a为口径面半径，P为口径场幅度分布函数控制参数，B+C＝1，ET＝20lgB，ET为口径场幅度分布函数在反射面边缘的锥销，ET＝-20dB，P＝1。分别分析表面随机误差均方根值ε为λ/20～λ/90的天线功率方向图计算结果。
2.仿真结果：
采用本发明的方法进行表面随机误差存在下的辐射功率方向图计算，并与传统方法进行比较。图2为采用传统方法与本发明方法得到的天线增益平均值随表面随机误差均方根值的变化曲线。图3为采用传统方法与本发明方法在表面随机误差均方值ε为λ/30时的天线平均功率方向图曲线。可以看出在表面随机误差均方值小于λ/30时，本发明方法与传统方法在天线增益上具有较好的吻合性，同时主瓣与近副瓣区域吻合性也非常好。
综上所述，本发明首先输入天线几何参数与电参数信息，计算理想天线远区辐射电场与远区辐射功率；其次，进行反射面三角形网格划分，分别计算单元节点位移与单元中心点位移的节点转换矩阵、总体转换矩阵；然后，计算电场相对单元节点位移的单元一阶导数、二阶Hessian阵、组集电场相对单元节点位移的总体一阶导数、二阶Hessian阵；再次，计算电场相对单元中心点位移的总体一阶导数、二阶Hessian阵；然后，根据用户输入的表面随机误差均方根值，计算天线远区辐射功率平均值；最后，判断辐射功率平均值是否满足电性能要求，并输出辐射功率方向图，以此指导反射面面板加工与制造。
与现有技术相比，本发明具有以下优点：
1.本发明基于单元节点位移二阶近似公式获得随机误差影响下反射面天线远区辐射功率平均值，在保证计算精度的前提下，避免了繁琐的公式推导，提高分析效率；
2.本发明从概率的角度出发，获得反射面天线远区辐射功率平均值，采用二阶近似公式保证了分析的准确性。
本实施方式中没有详细叙述的部分属本行业的公知的常用手段，这里不一一叙述。以上例举仅仅是对本发明的举例说明，并不构成对本发明的保护范围的限制，凡是与本发明相同或相似的设计均属于本发明的保护范围之内。

Claims (6)

1.基于单元节点位移的反射面天线平均功率方向图建模方法，其特征在于，包括如下步骤：
(1)输入反射面天线几何参数与电参数
输入用户提供的反射面天线几何参数与电参数；其中几何参数包括半径、焦距；电参数包括工作波长、自由空间波常数、口径场幅度分布函数、锥销电平、口径场形状指数以及包括天线增益、波瓣宽度、副瓣电平、指向精度在内的电性能要求；
(2)计算理想天线远区辐射电场
根据反射面天线几何参数与电参数，利用几何光学法计算理想天线远区辐射电场；
(3)计算理想天线远区辐射功率
根据理想天线远区辐射电场信息，通过下式计算理想天线远区辐射功率
Figure FDA0001838848440000011
其中，G0表示理想天线远区辐射功率，E0表示理想天线远区辐射电场，下标0代表理想天线，上标*表示取共轭运算；
(4)进行反射面三角形网格划分
根据用户提供的天线几何参数与工作波长，将反射面口径面划分为一系列等边三角形单元，等边三角形单元边长满足以下关系式
Figure FDA0001838848440000012
其中，λ为工作波长，l为口径面等边三角形边长；
(5)计算单元节点位移与单元中心点位移的节点转换矩阵
根据反射面三角形网格划分后三角形单元信息，通过下式计算单元节点位移与单元中心点位移的节点转换矩阵
Figure FDA0001838848440000021
其中，Ki表示第i个单元节点位移与其环绕的单元中心点位移之间的节点转换矩阵，mi为反射面三角形网格划分后环绕第i个单元节点的三角形单元数目，Ki的维数为mi行1列，上标T表示转置运算；
(6)组集单元节点位移与单元中心点位移的总体转换矩阵
根据反射面三角形网格划分后三角形单元信息，将单元节点位移与单元中心点位移的节点转换矩阵，按照三角形单元信息进行组集，通过下式组集单元节点位移与单元中心点位移的总体转换矩阵
Figure FDA0001838848440000022
其中，K表示单元节点位移与单元中心点位移的总体转换矩阵，Ki为第i个单元节点位移与单元中心点位移的节点转换矩阵，N为单元节点总数；
(7)计算电场相对单元节点位移的单元一阶导数、二阶Hessian阵；
(8)组集电场相对单元节点位移的总体一阶导数、二阶Hessian阵；
(9)计算电场相对单元中心点位移的总体一阶导数、二阶Hessian阵；
(10)输入表面随机误差均方根值
根据反射面天线面板加工制造误差，输入天线表面随机误差均方根值；
(11)计算天线远区辐射功率平均值；
(12)判断电性能是否满足要求
判断天线远区辐射功率平均值是否满足天线增益、波瓣宽度、副瓣电平、指向精度在内的电性能要求，如果满足要求则转至步骤(13)，否则转至步骤(14)；
(13)输出辐射功率方向图
当天线远区辐射功率平均值满足天线增益、波瓣宽度、副瓣电平、指向精度在内的电性能要求时，输出辐射功率方向图；
(14)更新表面随机误差均方根值
当天线远区辐射功率平均值不满足天线增益、波瓣宽度、副瓣电平、指向精度在内的电性能要求时，更新表面随机误差均方根值，转至步骤(10)。
2.根据权利要求1所述的基于单元节点位移的反射面天线平均功率方向图建模方法，其特征在于：所述步骤(2)中的几何光学法是一种基于口径场分布的高频近似方法，计算公式如下：
Figure FDA0001838848440000031
其中，E0表示理想天线远区辐射电场，∫表示积分运算，ρ表示反射面上点在口径面内的极坐标分量，Q(ρ)表示用户输入的反射面上点对应的口径场幅度分布函数，exp表示自然对数的指数运算，j为虚数单位，k为自由空间波常数，
Figure FDA0001838848440000032
为反射面上点在口径面内的位置矢量，
Figure FDA0001838848440000033
为远区观察点单位位置矢量，s表示天线口径面，ds表示在口径面内进行积分运算。
3.根据权利要求1所述的基于单元节点位移的反射面天线平均功率方向图建模方法，其特征在于：所述步骤(7)计算电场相对单元节点位移的单元一阶导数、二阶Hessian阵是根据用户提供的天线几何参数与电参数，结合三角形单元，通过下式电场相对单元节点位移的单元一阶导数
Gi＝[Gi,1,Gi,2,Gi,3]T
Figure FDA0001838848440000041
其中，Gi为电场相对单元节点位移的单元一阶导数，下标i表示第i个单元，Gi,l表示电场相对第i个单元上第l个节点的单元一阶导数分量，l＝1,2,3,∫表示积分运算，Q(ρ)表示用户输入的口径场幅度分布函数，ρ表示单元内部点在口径面内的极坐标分量，exp表示自然对数的指数运算，j为虚数单位，k为自由空间波常数，
Figure FDA0001838848440000042
为单元内部点在口径面内的位置矢量，
Figure FDA0001838848440000043
为远区观察点单位位置矢量，γl表示单元内部相对于第l个节点的单元形函数，ξ表示单元内部点在馈源坐标系下的角度分量，si表示第i个单元在口径面内的投影面积，dsi表示在第i个单元投影面积内进行积分运算；
根据用户提供的天线几何参数与电参数，结合三角形单元，通过下式电场相对单元节点位移的单元二阶Hessian阵
Figure FDA0001838848440000044
Figure FDA0001838848440000045
其中，Hi为电场相对单元节点位移的单元二阶Hessian阵，下标i表示第i个单元，Hi,uv表示电场相对第i个单元上第u、v个节点的单元二阶Hessian阵元素分量，u,v＝1,2,3,∫表示积分运算，Q(ρ)表示用户输入的口径场幅度分布函数，ρ表示单元内部点在口径面内的极坐标分量，exp表示自然对数的指数运算，j为虚数单位，k为自由空间波常数，
Figure FDA0001838848440000046
为单元内部点在口径面内的位置矢量，
Figure FDA0001838848440000047
为远区观察点单位位置矢量，γu、γv分别表示单元内部相对于第u、v个节点的单元形函数，ξ表示单元内部点在馈源坐标系下的角度分量，si表示第i个单元在口径面内的投影面积，dsi表示在第i个单元投影面积内进行积分运算。
4.根据权利要求1所述的基于单元节点位移的反射面天线平均功率方向图建模方法，其特征在于：所述步骤(8)组集电场相对单元节点位移的总体一阶导数、二阶Hessian阵是根据反射面三角形网格划分后三角形单元信息，通过下式组集电场相对单元节点位移的总体一阶导数、二阶Hessian阵
Figure FDA0001838848440000051
Figure FDA0001838848440000052
其中，GE表示电场相对单元节点位移的总体一阶导数，Gi为电场相对于单元节点位移的单元一阶导数，M表示三角形单元总数，A代表组集运算，HE表示电场相对单元节点位移的总体二阶Hessian阵，Hi为电场相对于单元节点位移的单元二阶Hessian阵。
5.根据权利要求1所述的基于单元节点位移的反射面天线平均功率方向图建模方法，其特征在于：所述步骤(9)计算电场相对单元中心点位移的总体一阶导数、二阶Hessian阵是根据单元节点位移与单元中心点位移的总体转换矩阵，结合电场相对单元节点位移的总体一阶导数、二阶Hessian阵，通过下式计算电场相对于单元中心点位移的总体一阶导数、二阶Hessian阵
G＝KGE
H＝KHEKT
其中，G为电场相对于单元中心点位移的总体一阶导数，H为电场相对于单元中心点位移的二阶Hessian阵，GE表示电场相对单元节点位移的总体一阶导数，HE表示电场相对单元节点位移的总体二阶Hessian阵，K表示单元节点位移与单元中心点位移的总体转换矩阵，上标T表示转置运算。
6.根据权利要求1所述的基于单元节点位移的反射面天线平均功率方向图建模方法，其特征在于：所述步骤(11)计算天线远区辐射功率平均值是根据理想天线远区辐射电场、远区辐射功率、电场相对单元中心点位移的单元一阶、二阶系数与表面随机误差均方根值，通过下式计算天线远区辐射功率平均值
Figure FDA0001838848440000061
其中，
Figure FDA0001838848440000062
表示天线远区辐射功率平均值，G0表示理想天线远区辐射功率，E0表示理想天线远区辐射电场，下标0代表理想天线，上标*表示取共轭运算，M为反射面上三角形单元总数，Hii、Hij、Hjj分别表示电场相对于单元中心点位移的总体二阶Hessian阵的第i行第i列分量、第i行第j列分量，第j行第j列分量，Gi电场相对于单元中心点位移的总体一阶导数第i个分量，σ表示用户输入的天线表面随机误差均方根值。
CN201811239147.0A 2018-10-23 2018-10-23 基于单元节点位移的反射面天线平均功率方向图建模方法 Active CN109408956B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811239147.0A CN109408956B (zh) 2018-10-23 2018-10-23 基于单元节点位移的反射面天线平均功率方向图建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811239147.0A CN109408956B (zh) 2018-10-23 2018-10-23 基于单元节点位移的反射面天线平均功率方向图建模方法

Publications (2)

Publication Number Publication Date
CN109408956A CN109408956A (zh) 2019-03-01
CN109408956B true CN109408956B (zh) 2022-12-30

Family

ID=65469406

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811239147.0A Active CN109408956B (zh) 2018-10-23 2018-10-23 基于单元节点位移的反射面天线平均功率方向图建模方法

Country Status (1)

Country Link
CN (1) CN109408956B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110110450B (zh) * 2019-05-11 2023-02-21 西安电子科技大学 基于解析机电耦合模型的刚性可展开天线装配误差分析方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7026989B1 (en) * 2004-01-23 2006-04-11 Itt Manufacturing Enterprises, Inc. Methods and apparatus for shaping antenna beam patterns of phased array antennas
CN107679336B (zh) * 2017-10-20 2021-02-05 西安电子科技大学 基于二阶近似公式的反射面天线表面随机误差分析方法

Also Published As

Publication number Publication date
CN109408956A (zh) 2019-03-01

Similar Documents

Publication Publication Date Title
CN109472066B (zh) 基于单元中心点位移的反射面天线随机误差分析方法
CN106886656B (zh) 一种三维阵列天线方向图旁瓣抑制方法
CN107679336B (zh) 基于二阶近似公式的反射面天线表面随机误差分析方法
CN105302962B (zh) 基于结构-电磁混合单元的反射面天线机电集成优化设计方法
CN112364467B (zh) 一种放松反射面天线远场分析电磁网格尺寸的方法
CN109063266A (zh) 一种基于并行贝叶斯优化的蒙皮天线机电协同设计方法
CN110059422B (zh) 频率选择表面曲面天线罩的电磁散射特性仿真方法
CN107357962A (zh) 一种基于自适应代理模型的天线罩肋截面尺寸优化方法
CN103513225B (zh) 一种基于空间增益的稀疏平面阵形优化方法
CN105160115A (zh) 基于近似与灵敏度分析的反射面天线机电集成优化设计方法
CN109408956B (zh) 基于单元节点位移的反射面天线平均功率方向图建模方法
CN109408958B (zh) 考虑面片拼合误差的伞状天线平均功率方向图建模方法
CN105024166B (zh) 基于子阵的平面阵列天线方向图综合方法
CN111276822A (zh) 一种天线方向图主瓣可控的天线阵列设计方法
CN108153954B (zh) 一种基于结构热类比的反射面天线温度快速重构方法
Guney et al. Input resistance calculation for circular microstrip antennas using adaptive neuro-fuzzy inference system
CN114169201B (zh) 面向电性能的反射面天线结构加权优化方法
CN109408957B (zh) 一种考虑确定性误差的伞状天线随机误差分析方法
CN109301500A (zh) 切比雪夫微带阵列天线的设计方法
CN110110450B (zh) 基于解析机电耦合模型的刚性可展开天线装配误差分析方法
Oğuzer et al. E‐polarized beam scattering by an open cylindrical PEC strip having an arbitrary “conical‐section” profile
CN106785470B (zh) 一种构架式天线索网机电集成制作方法
CN108268747B (zh) 基于力密度法的低栅瓣索网天线机电综合设计方法
Mutonkole Modelling of antenna responses
Chen et al. Machine-learning Assisted Synthesis for Series-Fed Microstrip Array

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