CN107239602B - 一种基于曲线拟合的探头天线模型快速计算方法 - Google Patents

一种基于曲线拟合的探头天线模型快速计算方法 Download PDF

Info

Publication number
CN107239602B
CN107239602B CN201710360973.XA CN201710360973A CN107239602B CN 107239602 B CN107239602 B CN 107239602B CN 201710360973 A CN201710360973 A CN 201710360973A CN 107239602 B CN107239602 B CN 107239602B
Authority
CN
China
Prior art keywords
probe
fitting
model
calculating
yaghjian
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
CN201710360973.XA
Other languages
English (en)
Other versions
CN107239602A (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.)
CETC 41 Institute
Original Assignee
CETC 41 Institute
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 CETC 41 Institute filed Critical CETC 41 Institute
Priority to CN201710360973.XA priority Critical patent/CN107239602B/zh
Publication of CN107239602A publication Critical patent/CN107239602A/zh
Application granted granted Critical
Publication of CN107239602B publication Critical patent/CN107239602B/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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling

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)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种基于曲线拟合的探头天线模型快速计算方法，该方法首先根据探头口径大小以及探头复反射系数计算出随频率变化的参数C0，然后根据C0与频率的分布，利用曲线拟合得到反映参数C0与频率关系的拟合函数式，将拟合系数作为探头参数存储，在近远场变换过程中，根据测量频率、拟合系数即可利用拟合函数式计算得到C0，然后根据新的Yaghjian模型计算得到探头H面方向图。本方法在进行近远场变换过程中探头方向图的计算复杂度显著低于求解二次复数方程，而且也避免了存储整个频段内的复反射系数数据。

Description

一种基于曲线拟合的探头天线模型快速计算方法
技术领域
本发明涉及一种基于曲线拟合的探头天线模型快速计算方法。
背景技术
自天线被发明以来，其在雷达、对抗、通信和导航等国防工业及民用领域的应用与日俱增，成为无线设备中不可或缺的部分。天线测量是伴着天线的设计出现的，是指导天线设计和验证检验天线性能的重要手段。对于天线来说，进行信号的辐射和接收是其核心价值所在，因此天线方向图特性也是其核心指标。为满足天线方向图特性测试的需求，天线自动化测试系统也应运而生，并伴随着天线设计制造技术的发展而逐步发展。
随着计算机技术的进步和现代测量手段的发展，天线测量技术也取得了重大进展，已形成多种测量方法，包括远场测量、近场测量和紧缩场测量等。
天线近场测量是用一个特性已知的探头，抽测天线近区某一表面上场的幅度及相位分布，通过严格的数学变换式确定天线的口面场和远场特性的过程。根据选取的测量面的不同，通常分为平面近场扫描、柱面近场扫描以及球面近场扫描。近场天线测试技术已经因其具有获得的信息量大、环境及电气随机干扰小、计算精度高、投资小、可全天候工作等优点得到广泛应用。近场测量作为一种间接测量方法，需要进行近场幅相数据至远场方向图的变换，而近场幅相数据的获取需要利用近场测量探头，探头的修正在近远场变换过程中必不可少。
目前近场测量中大多以标准开口矩形波导作为测量探头，探头补偿过程中需要探头的远场方向图，实际应用中不会对每个探头进行方向图测试然后修正，开口波导的形式简单，可由理论模型直接计算得到其方向图。
近场探头目前较为常用的理论模型是Stratton-Chu模型以及Yaghjian模型，两种模型计算探头的E面方向图相同，如图1所示：探头宽边长度为a，短边长度为b，以探头口面法线方向为Z轴建立坐标系，则其E面方向图可由式(1)计算：
Figure BDA0001300380300000011
其中，AE为方向图最大点(即θ＝0点)的幅值，计算归一化方向图时，可设置为1；
对于主模TE10模，归一化传播常数
Figure BDA0001300380300000012
k为波数(与探头工作频率f的关系为k＝2πc/f)；Γ为探头复反射系数，经验证，对于前半空间的E面方向图(近场测试时仅使用探头前半空间)可以取Γ＝0。
Stratton-Chu模型采用式(2)进行H面方向图的计算：
Figure BDA0001300380300000021
其中，AH＝-ik2abE0/8，该模型在|θ|＞30°时精度较差，根据文献“ApproximateFormulas for the Far Field and Gaionf Open-Ended Rectangular Waveguide”，在θ＝30°时，Stratton-Chu探头模型与实测方向图误差就达到了1dB左右。
Yaghjian在该文献中提出了电场积分方程计算H面方向图的方法，如式(3)所示：
Figure BDA0001300380300000022
该模型在|θ|＜60°的范围内达到与实测方向图差别在1dB之内的精度，超过60°后，其误差较大，特别是在接近90°时，该模型计算的方向图将迅速下降到0，因此Yaghjian又在Stratton-Chu模型基础上考虑了边缘电流的影响，提出了新的Yaghjian模型，即式(4)所示：
Figure BDA0001300380300000023
由于EH(θ)与式(1)中的EE(θ)在θ＝0时相同，所以可以得到：
Figure BDA0001300380300000024
其中，C0为正实数，与探头口径以及其复反射系数Γ有关，需要根据探头的端口输入功率与远场辐射功率关系式建立复数方程求解，探头输入功率如式(6)所示：
Figure BDA0001300380300000025
探头的辐射功率则可通过远场方向图积分得到，如式(7)所示：
Figure BDA0001300380300000031
然后根据式(8)建立关于C0的二次复数方程并进行求解：
P0＝Pr (8)
最后根据式(4)即可得到探头H面方向图。
由上述内容可以看出，现有探头模型的问题主要是H面方向图的计算，若根据式(2)或式(3)模型进行计算，其计算精度较差，特别是在毫米波频段，探头边缘尺寸相对探头口径已较大，其边缘电流影响更为明显，所以采用高精度的Yaghjian模型十分必要。
根据式(4)的Yaghjian模型进行方向图计算，需要先存储该探头的随频率变化的复反射系数测量数据，计算时根据探头工作频点进行插值，然后求解由式(6)和式(7)建立的二次复数方程求出C0，最后根据式(4)计算得到方向图。
然而，高精度的Yaghjian模型在计算过程中需要存储随频率变化的复反射系数测量数据，并进行二次复数方程求解，其过程复杂，模型计算效率低。
发明内容
本发明的目的在于提出一种基于曲线拟合的探头天线模型快速计算方法，以降低近远场变换过程中探头方向图的计算复杂度，同时避免存储整个频段内的复反射系数数据。
为了实现上述目的，本发明采用如下技术方案：
一种基于曲线拟合的探头天线模型快速计算方法，包括步骤：
a测量探头在其工作频段上的复反射系数Γ(fj)，其中，fj为工作频点；
b计算出每个工作频点fj对应的参数C0在探头工作频段上的分布即C0(fj)，其中，C0为根据探头口径大小以及探头复反射系数Γ(fj)计算出随频率变化的参数；
c利用C0(fj)及其工作频点fj进行曲线拟合，得到相应的拟合函数式，该拟合函数式中包含对应的拟合系数；
d将拟合系数以及探头口径作为探头固有参数进行存储；
e在近远场变换过程中，根据工作频率f以及步骤d中存储的拟合系数，按照步骤c中的拟合函数式计算得到参数C0；
f将步骤e中得到的参数C0代入新的Yaghjian模型计算得到探头H面方向图。
优选地，所述步骤c中拟合函数式如下：
Figure BDA0001300380300000041
其中，C0(f)表示参数C0与频率f的函数关系，ai表示拟合系数，fi表示频率f的i次幂，N表示曲线拟合的最高幂级数。
优选地，所述步骤c中的曲线拟合包括多项式拟合、指数拟合或三角函数拟合。
本发明具有如下优点：
本发明方法首先根据探头口径大小以及探头复反射系数计算出随频率变化的参数C0，然后根据C0与频率的分布利用曲线拟合得到反映C0与频率关系的拟合函数式，将拟合系数作为探头参数存储，在近远场变换过程中，根据测量频率、拟合系数即可利用拟合函数式计算得到C0，然后根据新的Yaghjian模型计算得到探头H面方向图。上述方法无需存储整个频段上的复反射系数Γ(fj)，只需存储拟合系数(适用于全频段)，减少探头固定存储数据。此外，探头模型计算过程中避免了二次复数方程的求解，只需计算多项式，计算简便易行。
附图说明
图1为探头及坐标系示意图；
图2为本发明中基于曲线拟合的探头天线模型快速计算方法的流程示意图。
具体实施方式
本发明的基本思想为：首先根据一个频段内探头天线的实测复反射系数，根据Yaghjian方法计算出有限频点内的C0参数，然后利用曲线拟合方法，拟合出该参数与工作频率的关系，在近远场变换的同时只需要进行曲线的计算即可，可大幅降低模型的计算复杂度。
下面结合附图以及具体实施方式对本发明作进一步详细说明：
如图2所示，一种基于曲线拟合的探头天线模型快速计算方法，包括步骤：
a测量探头在其工作频段上的复反射系数Γ(fj)，其中，fj为工作频点；
b计算出每个工作频点fj对应的参数C0在探头工作频段上的分布即C0(fj)，其中，C0为根据探头口径大小以及探头复反射系数Γ(fj)计算出随频率变化的参数；
具体的，参数C0的计算根据背景技术部分的公式(6)到公式(8)。
c利用C0(fj)及其工作频点fj进行曲线拟合，得到相应的拟合函数式，该拟合函数式中包含对应的拟合系数。曲线拟合包括多项式拟合、指数拟合或三角函数拟合等。
以N次多项式拟合为例，得到如式(9)所示的拟合函数式：
Figure BDA0001300380300000051
其中，C0(f)表示参数C0与频率f的函数关系，ai表示拟合系数，fi表示频率f的i次幂，N表示曲线拟合的最高幂级数。
d将拟合系数ai以及探头口径a，b作为探头固有参数进行存储；
e在近远场变换过程中，根据工作频率f以及步骤d中存储的拟合系数，按照步骤c中的拟合函数式(9)计算得到参数C0，C0值计算简单快速；
f将步骤e中得到的参数C0代入新的Yaghjian模型计算得到探头H面方向图；
其中，新的Yaghjian模型是Yaghjian在Stratton-Chu模型基础上考虑了边缘电流的影响提出的，新的Yaghjian模型公式具体可以参见本说明书背景技术部分公式(4)。
本发明利用曲线拟合公式求解探头方向图模型中未知参数，探头补偿时能够避免随频率变化的复反射系数的存储，并且避免复杂的二次复数方程的求解，仅需将几个简单的拟合系数作为与探头口径大小类似的固有参数存储，在计算时直接根据关系式得到方向图即可。
当然，以上说明仅仅为本发明的较佳实施例，本发明并不限于列举上述实施例，应当说明的是，任何熟悉本领域的技术人员在本说明书的教导下，所做出的所有等同替代、明显变形形式，均落在本说明书的实质范围之内，理应受到本发明的保护。

Claims (3)

1.一种基于曲线拟合的探头天线模型快速计算方法，其特征在于，包括步骤：
a测量探头在其工作频段上的复反射系数Γ(fj)，其中，fj为工作频点；
b计算出每个工作频点fj对应的参数C0在探头工作频段上的分布即C0(fj)，其中，C0为根据探头口径大小以及探头复反射系数Γ(fj)计算出随频率变化的参数；
c利用C0(fj)及其工作频点fj进行曲线拟合，得到相应的拟合函数式，该拟合函数式中包含对应的拟合系数；
d将拟合系数以及探头口径作为探头固有参数进行存储；
e在近远场变换过程中，根据工作频率f以及步骤d中存储的拟合系数，按照步骤c中的拟合函数式计算得到参数C0；
f将步骤e中得到的参数C0代入新的Yaghjian模型计算得到探头H面方向图，
新的Yaghjian模型是Yaghjian在Stratton-Chu模型基础上考虑了边缘电流的影响提出的，新的Yaghjian模型公式为：
Figure FDA0002426587220000011
2.根据权利要求1所述的一种基于曲线拟合的探头天线模型快速计算方法，其特征在于，所述步骤c中拟合函数式如下：
Figure FDA0002426587220000012
其中，C0(f)表示参数C0与频率f的函数关系，ai表示拟合系数，fi表示频率f的i次幂，N表示曲线拟合的最高幂级数。
3.根据权利要求1或2所述的一种基于曲线拟合的探头天线模型快速计算方法，其特征在于，所述步骤c中的曲线拟合包括多项式拟合、指数拟合或三角函数拟合。
CN201710360973.XA 2017-05-22 2017-05-22 一种基于曲线拟合的探头天线模型快速计算方法 Active CN107239602B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710360973.XA CN107239602B (zh) 2017-05-22 2017-05-22 一种基于曲线拟合的探头天线模型快速计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710360973.XA CN107239602B (zh) 2017-05-22 2017-05-22 一种基于曲线拟合的探头天线模型快速计算方法

Publications (2)

Publication Number Publication Date
CN107239602A CN107239602A (zh) 2017-10-10
CN107239602B true CN107239602B (zh) 2020-07-03

Family

ID=59985050

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710360973.XA Active CN107239602B (zh) 2017-05-22 2017-05-22 一种基于曲线拟合的探头天线模型快速计算方法

Country Status (1)

Country Link
CN (1) CN107239602B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112583502B (zh) 2019-09-27 2022-06-10 维沃移动通信有限公司 一种探头天线确定方法及装置
CN111025029A (zh) * 2019-11-18 2020-04-17 扬州船用电子仪器研究所(中国船舶重工集团公司第七二三研究所) 一种三倍频相控阵阵面近场方向图测试系统及方法
CN112417682B (zh) * 2020-11-20 2022-10-28 中国人民解放军63921部队 天线的远场辐射功率数据的参数拟合方法和装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1131643A1 (de) * 1998-11-16 2001-09-12 DeTeMobil Deutsche Telekom MobilNet GmbH Verfahren zur korrektur der frequenz- und längenabhängigen leitungsdämpfung bei fdr-messungen an hochfrequenzkabeln
CN101267062A (zh) * 2008-04-30 2008-09-17 西安电子科技大学 基于拟合变形反射面的天线电性能预测方法
CN101277278A (zh) * 2008-05-21 2008-10-01 清华大学 一种针对有轨交通的无线信道预测方法
CN102393513A (zh) * 2011-08-31 2012-03-28 中国电子科技集团公司第三十八研究所 一种基于自然分布场景和稀少定标器的极化定标技术

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1131643A1 (de) * 1998-11-16 2001-09-12 DeTeMobil Deutsche Telekom MobilNet GmbH Verfahren zur korrektur der frequenz- und längenabhängigen leitungsdämpfung bei fdr-messungen an hochfrequenzkabeln
CN101267062A (zh) * 2008-04-30 2008-09-17 西安电子科技大学 基于拟合变形反射面的天线电性能预测方法
CN101277278A (zh) * 2008-05-21 2008-10-01 清华大学 一种针对有轨交通的无线信道预测方法
CN102393513A (zh) * 2011-08-31 2012-03-28 中国电子科技集团公司第三十八研究所 一种基于自然分布场景和稀少定标器的极化定标技术

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
A NTFF transformation method in planar near-field antenna measurement derived from FDTD;Du Liuge 等;《ISAPE2012》;20121026;101-104 *
基于混合算法近远场变换中的探头补偿;陈玉林 等;《微波学报》;20140615;15-17 *
近远场变换中探头方向图的改进与比较;陈玉林 等;《雷达科学与技术》;20141015;第12卷(第5期);557-560,568 *

Also Published As

Publication number Publication date
CN107239602A (zh) 2017-10-10

Similar Documents

Publication Publication Date Title
CN113225147B (zh) 一种阵列天线总辐射功率的测量方法、装置和系统
CN111371513B (zh) 用于测试待测设备的天线阵列的方法及测试系统
CN107239602B (zh) 一种基于曲线拟合的探头天线模型快速计算方法
CN110632399B (zh) 用于球面近场测量数据的修正方法及天线方向图测量方法
CN104730503B (zh) 确定高分辨率sar参考目标rcs对定标影响的方法及补偿方法
CN107632208B (zh) 一种球面近场天线测量方法及系统
CN103616569A (zh) 一种毫米波平面近场测试相位修正方法
CN109374990B (zh) 一种天线相位中心校准方法
CN111366793B (zh) 一种减小截断误差的平面天线近场测量方法及系统
CN104391183A (zh) 近场测量天线远区场特性快速计算方法
CN110907904B (zh) 一种太阳法方位波束宽度修正方法
CN107247193B (zh) 天线近场测试方法及装置
CN103543425A (zh) 一种自动补偿网络分析仪测量面变化误差的方法
CN115792835A (zh) 基于探头补偿和相位中心补正的目标rcs近场测量方法
CN102608434B (zh) 一种毫米波黑体散射系数的测量方法
CN109241565B (zh) 天线球面近场测量中采样位置误差校正方法和装置
CN105116261B (zh) 一种测量波纹喇叭损耗的方法
CN112710903A (zh) 一种基于多探头球面场天线测试系统的相位中心测试方法
CN110445557B (zh) 深空测控干涉测量大口径天线指向标校方法及装置
CN105510365B (zh) 地面反射系数的现场测量方法
TWI511481B (zh) 移動終端總全向靈敏度量測之加速方法
CN113970670B (zh) 箔条空气混合介电常数测量方法、系统、设备、介质及终端
CN102508049A (zh) 一种基于Hankel函数外推的天线中场测试方法
CN115047257A (zh) 一种基于球面近场测量的天线自动化测量系统
CN113899956B (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