CN108959806B - 一种基于球面近场测量和球模式源的等效辐射建模方法 - Google Patents

一种基于球面近场测量和球模式源的等效辐射建模方法 Download PDF

Info

Publication number
CN108959806B
CN108959806B CN201810812658.0A CN201810812658A CN108959806B CN 108959806 B CN108959806 B CN 108959806B CN 201810812658 A CN201810812658 A CN 201810812658A CN 108959806 B CN108959806 B CN 108959806B
Authority
CN
China
Prior art keywords
spherical
measured
antenna
model
wave
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.)
Expired - Fee Related
Application number
CN201810812658.0A
Other languages
English (en)
Other versions
CN108959806A (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.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201810812658.0A priority Critical patent/CN108959806B/zh
Publication of CN108959806A publication Critical patent/CN108959806A/zh
Application granted granted Critical
Publication of CN108959806B publication Critical patent/CN108959806B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design 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)
  • Variable-Direction Aerials And Aerial Arrays (AREA)
  • Aerials With Secondary Devices (AREA)

Abstract

本发明属于天线测量技术领域,提供一种基于球面近场测量和球模式源的等效辐射建模方法,用以克服传统建模方法中理想偶极子的空间分布不变、缺乏灵活性、模型精度较低等问题。本发明以待测天线球面近场测量数据为基础,以球模式源为组成等效模型的基本单元,利用球面波展开和模式匹配法,建立等效模型,使得该模型能具备和待测天线相似的辐射特性,并且将球模式源所处位置的方位角和俯仰角进行变量分离,仅针对其方位角进行优化,通过这种一维优化的手段,可以在较短的时间内,提高模型精度,在工程中具有实用价值。

Description

一种基于球面近场测量和球模式源的等效辐射建模方法
技术领域
本发明属于天线测量技术领域,具体提供一种基于球面近场测量和球模式源的等效辐射建模方法。
背景技术
在天线测试领域,主要采用近场测试、远场测试两种方法,对于远场测试要求纯净的空间电磁环境及待测天线远场条件的测试距离,实施过程中存在诸多困难。近场测量在天线的近区范围内,不受远场测试中的距离效应和外界环境的影响,具有测试精度高、保密性强,全天候工作等优点;因此近场测量技术在天线近远场变换、天线诊断、电磁干扰预测和近场屏蔽效能计算都有重要的意义。
基于球面近场测量的等效辐射建模方法,是通过球面近场测量待测天线的近场数据,建立等效模型,使得等效模型能够产生和待测天线相似的辐射特性。现有的建模方法,其模型是由多个理想偶极子组成的,如文献“M.Serhir,P.Besnier,M.Drissim,"Anaccurate equi valent behavioral model of antenna radiation using a mode-matching technique based on sph erical near field measurements,"IEEETrans.Antennas Propag.,vol.56,no.1,pp.48-57,J an.2008”,其中每个理想偶极子的位置和激励被确定,等效模型就随之确定;但是,理想偶极子的位置是根据经验公式确定的,针对不同待测天线,其空间分布不变,缺乏灵活性,模型精度较低;并且由于该方法固有的缺陷,对理想偶极子的空间分布进行优化会耗费大量时间,无法在较短时间内提高精度。
发明内容
本发明的目的在于针对上述背景技术的缺陷,提供一种基于球面近场测量和球模式源的等效辐射建模方法;本发明以待测天线球面近场测量数据为基础,利用球面波展开和模式匹配法,建立等效模型,使得该模型能具备和待测天线相似的辐射特性,并且通过快速优化手段,在较短的时间内,提高模型精度。
为实现上述目的,本发明采用的技术方案为:
一种基于球面近场测量和球模式源的等效辐射建模方法,其特征在于,包括以下步骤:
步骤1.采用近场测量探头对待测天线的辐射场进行采样,测得的辐射场记为
Figure BDA0001739558970000021
并根据待测天线产生的辐射场在坐标原点的球面波展开式得到球面波函数的系数Qj,所述展开式为:
Figure BDA0001739558970000022
其中,
Figure BDA0001739558970000023
表示位置矢量,k为波数,η为波阻抗,Jt为球面波模式截断数,Qj为球面波函数的系数,
Figure BDA0001739558970000024
为归一化球面波函数,j=1,2,3,...,Jt
步骤2.在包围待测辐射体的最小虚拟球面上,采用LS个球模式源构建模型,
Figure BDA00017395589700000218
rmin为包围待测辐射体的最小虚拟球面的半径;第g个球模式源在坐标原点的球面波展开式为:
Figure BDA0001739558970000025
其中,
Figure BDA0001739558970000026
Figure BDA0001739558970000027
对应表示第g个球模式源的待定系数,
Figure BDA0001739558970000028
Figure BDA0001739558970000029
Figure BDA00017395589700000210
为第g、g=1,2,3,...,LS个球模式源的球面波函数转移系数;
Figure BDA00017395589700000211
则得到矩阵方程:
Figure BDA00017395589700000212
Figure BDA00017395589700000213
Figure BDA00017395589700000214
Figure BDA00017395589700000215
Figure BDA00017395589700000216
Figure BDA00017395589700000217
对等效模型的球模式源的空间分布进行优化,建立目标函数:
fitness=min{σmod}
其中:
Figure BDA0001739558970000031
Figure BDA0001739558970000032
为个体,采用遗传算法求解上述目标函数,得到最优解
Figure BDA0001739558970000033
即完成建模。
本发明的有益效果在于:
本发明提供一种基于球面近场测量和球模式源的等效辐射建模方法,提出了球模式源的概念,并对球模式源的空间分布进行优化,大幅提高模型精度;且其空间分布优化过程中,对球模式源的方位角
Figure BDA0001739558970000034
和俯仰角θ进行变量分离,仅需要对球模式源的方位角
Figure BDA0001739558970000035
进行优化,即一维优化,大大减少了计算时间。由此可知,相较于传统方法中理想偶极子的空间分布是根据经验公式给出的、无法进行快速优化,本发明模型精度高、建模时间短。
附图说明
图1为本发明基于球面近场测量和球模式源的等效辐射建模方法的原理示意图。
图2为实施例1中待测二元贴片天线阵结构示意图。
图3为实施例1中天线阵与其等效模型的辐射场对比图。
图4为实施例1中原天线、传统模型、本发明模型三者的辐射场对比图;其中,(a)(b)为待测天线辐射场,(c)(d)为传统等效模型辐射场,(e)(f)为本发明中的等效模型辐射场。
图5为实施例2中待测半波天线结构示意图。
图6为实施例2中天线与其等效模型的辐射场对比图。
图7为实施例2中原天线、传统模型、本发明模型三者的辐射场对比图;其中,(a)(b)为待测天线辐射场,(c)(d)为传统等效模型辐射场,(e)(f)为本发明中的等效模型辐射场。
具体实施方式
下面结合附图和实施例对本发明做进一步详细说明。
本发明提供一种基于球面近场测量和球模式源的等效辐射建模方法,其测量模型如图1所示,在本实施例中,将待测天线放置在球面近场测量转台上,在包围待测辐射体的半径为1米测量球面Ωm上,按照方位角步进
Figure BDA0001739558970000041
和俯仰角步进Δθ=5°,用连接矢量网络分析仪的近场测量探头对待测天线的辐射场进行采样,包括幅度和相位,测得的辐射场记为
Figure BDA0001739558970000042
根据球面波展开理论,待测天线产生的辐射场可以在坐标原点O进行球面波展开如下:
Figure BDA0001739558970000043
其中,
Figure BDA0001739558970000044
表示位置矢量,k为波数,η为波阻抗,Jt为球面波模式截断数,Qj为球面波函数的系数,
Figure BDA0001739558970000045
为归一化球面波函数;
本发明提出一种新的理想电磁辐射源—球模式源,球模式源的定义如下:
一种假想点源,且其辐射场
Figure BDA0001739558970000046
在所处位置进行球面波展开如下:
Figure BDA0001739558970000047
其中,α1、α2、α3、α4、α5,和α6为球模式源的待定系数;
本发明所建立的等效模型由球模式源组成,如图1所示,球模式源分布在包围待测辐射体的最小虚拟球面Ωf上,Ωf和Ωm是同心圆,等效模型的辐射场
Figure BDA0001739558970000048
的球面波展开如下:
Figure BDA0001739558970000049
其中,
Figure BDA00017395589700000416
为球模式源的个数,rmin为包围待测辐射体的最小虚拟球面Ωf的半径,
Figure BDA00017395589700000410
表示向上取整运算,
Figure BDA00017395589700000411
Figure BDA00017395589700000417
对应表示第g个球模式源的待定系数;
Figure BDA00017395589700000412
Figure BDA00017395589700000413
是球面波函数转移系数,其用于将球模式源在所处位置的球面波展开转换为在坐标原点O点的球面波展开;
至此,待测辐射体和等效模型在坐标原点O点的球面波展开都已给出,分别为式(1)和式(3),令式(1)等于式(3),得到:
Figure BDA00017395589700000414
根据球面波函数的正交性,式(4)左右两边的球面波函数系数应该分别相等,即:
Figure BDA00017395589700000415
将式(4)写成矩阵方程的形式:
Figure BDA0001739558970000051
其中:
Figure BDA0001739558970000052
Figure BDA0001739558970000053
Figure BDA0001739558970000054
Figure BDA0001739558970000055
Figure BDA0001739558970000056
其中,T为矩阵转置符号;
利用最小二乘法求解矩阵方程(6):
Figure BDA0001739558970000057
其中,
Figure BDA0001739558970000058
为等效模型中的未知系数;此时等效模型的球面波展开系数
Figure BDA0001739558970000059
为:
Figure BDA00017395589700000510
等效模型辐射场与待测辐射体辐射场的误差可以由球面波展开系数的差值决定,将该误差记作:
Figure BDA00017395589700000511
其中,||·||1为平均差算符;对组成等效模型的球模式源的空间分布进行优化,目标函数为:
fitness=min{σmod} (10)
Figure BDA00017395589700000512
中的元素
Figure BDA00017395589700000513
b=1,2,3,4,5,6都具有相似的形式:
Figure BDA0001739558970000061
其中,s,m,n与b有对应关系b=2{n(n+1)+m-1}+s,且s=1或2,n>0,-n≤m≤n;i为虚数单位
σ,l,v与j有对应关系j=2{v(v+1)+l-1}+σ,且σ=1或2,v>0,-v≤l≤v;
Figure BDA0001739558970000062
是第g个球模式源在球坐标系下的坐标;
Figure BDA0001739558970000063
Figure BDA0001739558970000064
是球面波旋转函数,
Figure BDA0001739558970000065
是球面波函数平移函数;
Figure BDA0001739558970000066
矩阵写成两个矩阵的哈达马积(Hadamard product):
Figure BDA0001739558970000067
其中,
Figure BDA0001739558970000068
Figure BDA0001739558970000069
矩阵元素形式分别是
Figure BDA00017395589700000610
Figure BDA00017395589700000611
矩阵
Figure BDA00017395589700000612
只与球模式源的方位角
Figure BDA00017395589700000613
有关,而矩阵
Figure BDA00017395589700000614
只与球模式源的俯仰角θ0有关;组成等效模型的LS个球模式源的俯仰角分别为
Figure BDA00017395589700000615
而方位角则以式(10)为目标函数,通过遗传算法进行优化求解;至此,输出最优解
Figure BDA00017395589700000616
即完成模型建立。
在上述优化过程中,仅对方位角进行优化,所以本发明称之为一维优化。
另外,需要说明的是:本发明提出一种新的理想电磁辐射源—球模式源,也可转换为现有理想偶极子进行应用,球模式源的系数与理想偶极子的偶极矩具有如下关系:
Figure BDA00017395589700000617
Figure BDA00017395589700000618
Figure BDA00017395589700000619
Figure BDA00017395589700000620
Figure BDA0001739558970000071
Figure BDA0001739558970000072
其中,de,x,de,y,de,z,dm,x,dm,y,和dm,z分别代表沿x轴放置的电偶极子偶极矩,沿y轴放置的电偶极子偶极矩,沿z轴放置的电偶极子偶极矩,沿x轴放置的磁偶极子偶极矩,沿y轴放置的磁偶极子偶极矩,和沿z轴放置的磁偶极子偶极矩;通过上述6个转换式,一个球模式源可以等价转换为六个理想偶极子,由球模式源组成的等效模型也可以等价转换成由理想偶极子组成的等效模型。
在具体实施方式中采用惯用的均方根误差对模型精度进行对比,均方根误差定义如下:
Figure BDA0001739558970000073
其中,
Figure BDA0001739558970000074
Figure BDA0001739558970000075
其中,Eθ
Figure BDA0001739558970000076
是待测辐射体的电场和磁场分量,Eθ_M
Figure BDA0001739558970000077
是等效模型的电场和磁场分量。
实施例1
本实施例以工作在3GHz的二元贴片天线阵为例进行建模;该天线介质基板的相对介电常数为2.2,介质基板厚度为2.87mm;天线尺寸为L1=80mm,L2=46.648mm,W1=100mm,W2=31.1807mm,d1=16.1mm,d2=40mm,d3=9.4096mm,d4=16.676mm;馈电相位差为九十度;通过本发明建立等效模型由12个球模式源组成,可以等效转换为72个理想偶极子,如图3所示为该天线阵与其等效模型的辐射场对比图,均方根误差为6.01%,计算机建模时间(不包括球面近场测试时间)为303秒;如图4所示为原天线、传统模型、本发明模型三者的辐射场对比图;传统模型的均方根误差为13.08%;由此可见,本发明模型精度大幅提高,且建模时间短。
实施例2
本实施例亦一个放置在金属块旁的半波天线为例进行建模,其工作频率为5GHz;本实施例旨在仿真受环境影响的天线辐射特性;通过本发明建立等效模型由12个球模式源组成,可以等效转换为72个理想偶极子,如图6所示为该辐射体与其等效模型辐射场的对比图,均方根误差为5.04%,计算机建模时间(不包括球面近场测试时间)为263秒;如图4所示为原天线、传统模型、本发明模型三者的辐射场对比图;传统模型的均方根误差为13.52%;由此可见,本发明模型精度大幅提高,且建模时间短。
以上所述,仅为本发明的具体实施方式,本说明书中所公开的任一特征,除非特别叙述,均可被其他等效或具有类似目的的替代特征加以替换;所公开的所有特征、或所有方法或过程中的步骤,除了互相排斥的特征和/或步骤以外,均可以任何方式组合。

Claims (1)

1.一种基于球面近场测量和球模式源的等效辐射建模方法,其特征在于,包括以下步骤:
步骤1.采用近场测量探头对待测天线的辐射场进行采样,测得的辐射场记为
Figure FDA0001739558960000011
并根据待测天线产生的辐射场在坐标原点的球面波展开式得到球面波函数的系数Qj,所述展开式为:
Figure FDA0001739558960000012
其中,
Figure FDA0001739558960000013
表示位置矢量,k为波数,η为波阻抗,Jt为球面波模式截断数,Qj为球面波函数的系数,
Figure FDA0001739558960000014
为归一化球面波函数,j=1,2,3,...,Jt
步骤2.在包围待测辐射体的最小虚拟球面上,采用LS个球模式源构建模型,
Figure FDA0001739558960000015
rmin为包围待测辐射体的最小虚拟球面的半径;第g个球模式源在坐标原点的球面波展开式为:
Figure FDA0001739558960000016
其中,
Figure FDA0001739558960000017
Figure FDA0001739558960000018
对应表示第g个球模式源的待定系数,
Figure FDA0001739558960000019
Figure FDA00017395589600000110
Figure FDA00017395589600000111
为第g、g=1,2,3,...,LS个球模式源的球面波函数转移系数;
Figure FDA00017395589600000112
则得到矩阵方程:
Figure FDA00017395589600000113
Figure FDA00017395589600000114
Figure FDA00017395589600000115
Figure FDA00017395589600000116
Figure FDA00017395589600000117
Figure FDA00017395589600000118
对等效模型的球模式源的空间分布进行优化,建立目标函数:
fitness=min{σmod}
其中:
Figure FDA0001739558960000021
Figure FDA0001739558960000022
为个体,采用遗传算法求解上述目标函数,得到最优解
Figure FDA0001739558960000023
即完成建模。
CN201810812658.0A 2018-07-23 2018-07-23 一种基于球面近场测量和球模式源的等效辐射建模方法 Expired - Fee Related CN108959806B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810812658.0A CN108959806B (zh) 2018-07-23 2018-07-23 一种基于球面近场测量和球模式源的等效辐射建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810812658.0A CN108959806B (zh) 2018-07-23 2018-07-23 一种基于球面近场测量和球模式源的等效辐射建模方法

Publications (2)

Publication Number Publication Date
CN108959806A CN108959806A (zh) 2018-12-07
CN108959806B true CN108959806B (zh) 2022-03-15

Family

ID=64464554

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810812658.0A Expired - Fee Related CN108959806B (zh) 2018-07-23 2018-07-23 一种基于球面近场测量和球模式源的等效辐射建模方法

Country Status (1)

Country Link
CN (1) CN108959806B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112100789A (zh) * 2019-05-28 2020-12-18 北京邮电大学 一种基于模式展开的球面近场采样点相位还原方法及系统
CN110632399B (zh) * 2019-09-19 2020-12-22 电子科技大学 用于球面近场测量数据的修正方法及天线方向图测量方法
CN111553046B (zh) * 2020-03-13 2021-03-30 电子科技大学 一种基于球面波展开和源重构的天线辐射计算方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2640821A1 (fr) * 1988-12-16 1990-06-22 Thomson Csf Antenne a couverture tridimensionnelle et balayage electronique, du type reseau volumique rarefie aleatoire
CN104063544A (zh) * 2014-06-25 2014-09-24 电子科技大学 近场效应误差分析方法
CN104408276A (zh) * 2014-09-15 2015-03-11 电子科技大学 一种采样远场方向图诊断阵列天线失效阵元的方法
CN105445566A (zh) * 2015-11-13 2016-03-30 西北工业大学 球面多探头天线测试数据处理方法
CN105808800A (zh) * 2014-12-30 2016-07-27 中国舰船研究设计中心 一种电子设备机柜泄露辐射仿真预测方法
CN106707037A (zh) * 2017-01-23 2017-05-24 电子科技大学 一种材料电磁特性参数无损反射测量方法及装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2640821A1 (fr) * 1988-12-16 1990-06-22 Thomson Csf Antenne a couverture tridimensionnelle et balayage electronique, du type reseau volumique rarefie aleatoire
CN104063544A (zh) * 2014-06-25 2014-09-24 电子科技大学 近场效应误差分析方法
CN104408276A (zh) * 2014-09-15 2015-03-11 电子科技大学 一种采样远场方向图诊断阵列天线失效阵元的方法
CN105808800A (zh) * 2014-12-30 2016-07-27 中国舰船研究设计中心 一种电子设备机柜泄露辐射仿真预测方法
CN105445566A (zh) * 2015-11-13 2016-03-30 西北工业大学 球面多探头天线测试数据处理方法
CN106707037A (zh) * 2017-01-23 2017-05-24 电子科技大学 一种材料电磁特性参数无损反射测量方法及装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
口径面天线辐射近区场强计算;张加坤;《电子科技大学学报》;20050225;19-21+52 *

Also Published As

Publication number Publication date
CN108959806A (zh) 2018-12-07

Similar Documents

Publication Publication Date Title
CN109061323B (zh) 一种采用球面幅度扫描的近场天线测量方法
CN108959806B (zh) 一种基于球面近场测量和球模式源的等效辐射建模方法
CN109581078B (zh) 一种适用于半空间环境中天线的方向图测量系统及方法
CN110988499B (zh) 一种基于无相位近场测量的天线辐射特性获取方法
CN112798876B (zh) 一种利用插值算法的天线近场测量方法及装置
CN107085171B (zh) 一种用于变电站局部放电测向的特高频相控阵的优化布置方法
Zhang et al. An efficient probe calibration based near-field-to-near-field transformation for EMI diagnosis
Hu et al. A near-field to far-field RCS measurement method for multiple-scattering target
CN113533867A (zh) 基于Fourier插值的远场方向图快速测量方法
Varela et al. Under-sampled spherical near-field antenna measurements with error estimation
CN109286069A (zh) 基于偶极子的阵列天线场源重构方法及装置
Neitz et al. 3-D monostatic RCS determination from multistatic near-field measurements by plane-wave field synthesis
Varela et al. Fast spherical near-field to far-field transformation for offset-mounted antenna measurements
CN115524541B (zh) 基于准远场模式滤波的任意极化线天线阵方向图测量方法
CN112068101B (zh) 基于模式滤波的目标散射分离方法
Zou et al. A field iterative method for efficient source reconstruction based on magnitude-only and single-plane near-field scanning
Yaccarino et al. A comparison of conventional and phaseless planar near-field antenna measurements: the effect of probe position errors
Neitz et al. Antenna gain determination in the near-field using propagating plane-wave expansions
Wang et al. Correction for probe-positioning errors in digital holographic reconstruction phase-less near-field measurement
Ying et al. Research on Undersampling Application on Narrowband Antenna in Planner Near-field Measurement
CN115356550A (zh) 基于多极子与幅度项补偿的线阵天线近场测量方法
CN117591778A (zh) 一种基于相位迭代的球面偶极子阵列的相位恢复方法
Nie et al. Study on Multi-radiation Source Equivalence and Prediction
Hirose et al. 3D Absolute Gain Pattern Measurements Within a Strip Region Using a New Multi-Cut Fresnel Field to Far-Field Cylindrical Transformation
Capozzoli et al. Prolate function expansion of circularly supportée aperture fields in near-field antenna characterization

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20220315