CN112417682A - 天线的远场辐射功率数据的参数拟合方法和装置 - Google Patents

天线的远场辐射功率数据的参数拟合方法和装置 Download PDF

Info

Publication number
CN112417682A
CN112417682A CN202011308932.4A CN202011308932A CN112417682A CN 112417682 A CN112417682 A CN 112417682A CN 202011308932 A CN202011308932 A CN 202011308932A CN 112417682 A CN112417682 A CN 112417682A
Authority
CN
China
Prior art keywords
far
antenna
field radiation
parameter
power data
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.)
Granted
Application number
CN202011308932.4A
Other languages
English (en)
Other versions
CN112417682B (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.)
63921 Troops of PLA
Original Assignee
63921 Troops of PLA
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 63921 Troops of PLA filed Critical 63921 Troops of PLA
Priority to CN202011308932.4A priority Critical patent/CN112417682B/zh
Publication of CN112417682A publication Critical patent/CN112417682A/zh
Application granted granted Critical
Publication of CN112417682B publication Critical patent/CN112417682B/zh
Active 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)

Abstract

本发明涉及天线测量技术领域,提供一种天线的远场辐射功率数据的参数拟合方法和装置,包括:获取天线的远场辐射功率数据集,所述远场辐射功率与所述天线的扫描角度一一对应;根据所述远场辐射功率数据集,合理选定待拟合参数,建立远场辐射参数模型;确定待拟合参数初值,并根据间接平差理论计算参数初值的改正量;在所述改正量满足预设阈值条件时,根据改正量和参数初值确定参数拟合结果。本发明的拟合方法适应性强、估计准确、收敛迅速,尤其是圆口径天线,可以准确地刻画圆口径天线远场辐射特性。

Description

天线的远场辐射功率数据的参数拟合方法和装置
技术领域
本发明涉及天线测量技术领域,尤其涉及一种天线的远场辐射功率数据的参数拟合方法和装置。
背景技术
圆形口径天线在许多领域都有广泛的应用。在航天测控领域,几乎所有地面站天线和航天器高增益天线都会采用圆形口径设计。圆形口径天线的远场辐射特性具有旋转对称性,理论辐射方向图也有形式简单的解析表达。
对于已经设计、制造、装配完成的天线,通常会采用实测方式获取其远场辐射功率,以确认天线的增益性能、波束宽度、电轴偏差等关键特性。有时还需要根据天线的辐射特性,进一步进行其他复杂的积分计算,例如计算天线指向接近某些天体(月球、太阳、火星等)时引起的接收噪声温度增量。然而,由于实测功率数据存在测量误差,且分辨率有限,不便于准确、直接地量化其辐射特性;此外,含有误差的离散测量数据也不利于实现高精度积分计算。为此,需要选定合适参数构建天线辐射模型,基于实测数据采用最优估计方法对模型参数进行估计,用以刻画天线远场辐射特性,亦方便其他复杂计算使用。
工程中较为关心的是天线半功率波束宽度这一量级范围内的辐射特性。对圆形口径天线,许多文献和实践中(在精度要求不高的情况下)经常采用二次函数模型拟合上述范围内的远场辐射功率数据集,其优点是模型简单、计算方便;但由于圆形口径天线的理论辐射特性并非抛物线形状,采用二次函数模型无法准确拟合,导致拟合后残差较大,求得的辐射模型与实测方向图存在差异。换言之,采用二次函数模型进行参数拟合所得的增益特性、电轴偏差等并不准确,若利用得到的辐射模型进行其他积分计算,计算结果也将存在一定误差。
发明内容
基于此,本发明实施例提供一种天线的远场辐射功率数据的参数拟合方法和装置,以解决传统方法精度低的问题。
本发明实施例的第一方面,提供一种天线的远场辐射功率数据的参数拟合方法,包括:
确定天线的远场辐射功率数据集,所述远场辐射功率与所述天线的扫描角度一一对应;
根据所述远场辐射功率数据集确定待拟合参数,并建立远场辐射参数模型;
确定所述待拟合参数的初值,并根据间接平差理论计算所述初值的改正量;
在所述改正量满足预设阈值条件时,根据所述改正量和所述初值确定参数拟合结果。
可选的,所述确定天线的远场辐射功率数据集,包括:
获取天线在所述扫描角度θ0n处的远场辐射功率的测量值P0n,n=1,2,…,N;根据
{(θi,Pi)|(θi,Pi)∈{(θ0n,P0n)},Pi≥max(P0n)-10}
确定所述远场辐射功率数据集;其中,i=1,2,…,I,I≤N。
可选的,所述根据所述远场辐射功率数据集确定待拟合参数,并建立远场辐射参数模型,包括:
确定待拟合参数k、α和C,所述远场辐射功率数据集中的功率值Pi与对应的扫描角度θi满足下式:
Figure BDA0002789165550000021
其中,vi为测量和模型误差,J1为一阶贝塞尔函数,k为表征天线波束宽度的因子,α为所述天线的视轴方向对应的角度,C为所述天线的远场辐射功率的最大值;所述待拟合参数构成参数向量X=[k,α,C]T
可选的,所述确定所述待拟合参数的初值,包括:
通过
Figure BDA0002789165550000031
确定所述参数向量X的初值,得到向量初值X0;其中,D为所述天线的口径,λ为所述天线的工作波长,Pi为所述远场辐射功率数据集中的功率值。
可选的,所述根据间接平差理论计算所述初值的改正量,包括:
计算所述远场辐射参数模型对所述参数向量在所述向量初值X0处的Jacobi矩阵B:
Figure BDA0002789165550000032
其中,
Figure BDA0002789165550000033
Figure BDA0002789165550000034
其中,J2为二阶贝塞尔函数;
根据间接平差理论,通过
Figure BDA0002789165550000035
确定所述向量初值的改正量
Figure BDA0002789165550000036
其中,l为所述远场辐射功率数据集中的功率值与基于所述远场辐射参数模型计算的近似值之差,公式如下:
Figure BDA0002789165550000037
其中,P1,P2,...,PI均为所述远场辐射功率数据集中的功率值,θ12,...,θI分别为P1,P2,...,PI对应的所述扫描角度。
可选的,所述天线的远场辐射功率数据的参数拟合方法还包括:
在所述改正量不满足预设阈值条件时,取
Figure BDA0002789165550000041
作为新的初值X'0,X0为所述待拟合参数的初值;
根据间接平差理论计算所述新的初值X'0的改正量
Figure BDA0002789165550000042
在改正量
Figure BDA0002789165550000043
满足所述预设阈值条件时,根据
Figure BDA0002789165550000044
确定所述参数拟合结果
Figure BDA0002789165550000045
否则,迭代计算直至改正量
Figure BDA0002789165550000046
满足所述预设阈值条件。
可选的,所述天线的远场辐射功率数据的参数拟合方法还包括:
根据
Figure BDA0002789165550000047
确定所述参数拟合结果的中误差,其中,
Figure BDA0002789165550000048
为所述向量初值的改正量,I为所述远场辐射功率数据集的元素总数,B为所述远场辐射参数模型对所述参数向量X在向量初值X0处的Jacobi矩阵,l为所述远场辐射功率数据集中的功率值与基于所述远场辐射参数模型计算的近似值之差。
本发明实施例的第二方面,提供一种天线的远场辐射功率数据的参数拟合装置,包括:
数据集确定模块,用于确定天线的远场辐射功率数据集,所述远场辐射功率与所述天线的扫描角度一一对应;
辐射模型建立模块,用于根据所述远场辐射功率数据集确定待拟合参数,并建立远场辐射参数模型;
改正量计算模块,用于确定所述待拟合参数的初值,并根据间接平差理论计算所述初值的改正量;
拟合结果确定模块,用于在所述改正量满足预设阈值条件时,根据所述改正量和所述初值确定参数拟合结果。
本发明实施例的第三方面,提供一种天线的远场辐射功率数据的参数拟合装置,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现如上述实施例的第一方面提供的任一项所述的天线的远场辐射功率数据的参数拟合方法的步骤。
本发明实施例的第四方面提供了一种计算机可读存储介质,所述计算机可读存储介质存储计算机程序,所述计算机程序被处理器执行时实现如实施例的第一方面提供的任一项所述的天线的远场辐射功率数据的参数拟合方法的步骤。
本发明实施例的天线的远场辐射功率数据的参数拟合方法和装置与现有技术相比存在的有益效果是:
本发明主要针对圆口径天线的远场辐射理论特性,合理选定待拟合参数,并建立其远场辐射功率的参数模型,并采用间接平差算法进行参数估计,得到参数拟合结果,本发明的拟合方法适应性强、估计准确、收敛迅速,可以较为准确地刻画圆口径天线远场辐射特性,适应于天线指向偏差标定或其他相关复杂计算场景对远场高精度辐射模型建立和求解的需求。
附图说明
图1是本发明实施例提供的一种天线的远场辐射功率数据的参数拟合方法的实现流程示意图;
图2是本发明实施例提供的一种圆口径天线在不同扫描角度下的远场增益数据图示;
图3是本发明实施例提供的采用本实施例方法进行圆口径天线的远场辐射数据的参数拟合的结果图示和残差图示;
图4是本发明实施例提供的采用传统方法进行圆口径天线的远场辐射数据的参数拟合的结果图示和残差图示;
图5是本发明实施例提供的一种天线的远场辐射功率数据的参数拟合装置的结构示意图;
图6是本发明实施例提供的另一种天线的远场辐射功率数据的参数拟合装置的结构示意图。
具体实施方式
以下描述中,为了说明而不是为了限定,提出了诸如特定系统结构、技术之类的具体细节,以便透彻理解本发明实施例。然而,本领域的技术人员应当清楚,在没有这些具体细节的其它实施例中也可以实现本发明。
为了说明本发明所述的技术方案,下面通过具体实施例来进行说明。
参见图1,为本实施例提供的天线的远场辐射功率数据的参数拟合方法的一个实施例实现流程示意图,详述如下:
步骤S101,确定天线的远场辐射功率数据集,所述远场辐射功率与所述天线的扫描角度一一对应。
本发明主要对圆口径天线的远场辐射功率数据进行参数拟合,构建天线辐射特性的参数模型,刻画天线辐射场、电轴偏差等关键特性,可用于天线指向偏差标定,所得模型可用于进行与天线辐射特性相关的其他复杂计算。
在实际应用中,天线的辐射特性通常是通过对天线方向图的一个或多个切面测量得到。在圆口径天线辐射特性旋转对称的基本前提下,参数拟合实际上是对单个切面的一维测量数据的拟合。此外,实际应用的圆口径天线通常波束较窄(一般小于1°),其远场理论方向图表达式中的sinθ可以近似为θ,近似误差在工程关心的半功率波束宽度尺度上可以忽略。
具体的,设圆口径天线的口径为D(m),工作波长为λ(m),其远场辐射功率在θ0n角度处的测量值为P0n(n=1,2,…,N),其中θ0n的单位为度(°),P0n的单位为dBm(或换算为增益,单位为dBi)。
在一个实施例中,步骤S101所述的获取天线的远场辐射功率数据集的具体实现流程可以包括:
首先,获取天线在所述扫描角度θ0n处的远场辐射功率的测量值P0n,n=1,2,…,N。
然后,根据
{(θi,Pi)|(θi,Pi)∈{(θ0n,P0n)},Pi≥max(P0n)-10}
确定所述远场辐射功率数据集;其中,i=1,2,…,I,I≤N。
选取满足预设条件的测量值Pi,可以更加精准地建立模型,提高参数拟合精度。
步骤S102,根据所述远场辐射功率数据集确定待拟合参数,并建立远场辐射参数模型。
远场辐射参数模型为天线的远场辐射功率模型,以确认天线的增益特性、波束宽度、电轴偏差等关键特性。应理解,在模型中需要有多个参数约束功率与对应角度的关系,使建立的远场辐射参数模型更加准确,例如天线波束宽度、天线的视轴方向对应的角度、天线的远场辐射功率的最大值等等,本实施例选定上述特征参数,基于圆口径天线的理论方向图和窄波束情况下的角度近似建立参数模型,通过参数拟合估计,确定最佳的参数值,从而更加精准地逼近实际的辐射特性。
在一个实施例中,步骤S102所述的根据所述远场辐射功率数据集确定待拟合参数,并建立远场辐射参数模型的具体实现流程可以包括:
确定待拟合参数k、α和C,所述远场辐射功率数据集中的功率值Pi与对应的扫描角度θi满足下式:
Figure BDA0002789165550000071
其中,vi为测量和模型误差,J1为一阶贝塞尔函数,k为表征天线波束宽度的因子,α为所述天线的视轴方向对应的角度,C为所述天线的远场辐射功率的最大值;所述待拟合参数构成参数向量X=[k,α,C]T,即k、α和C这三个参数是为了拟合天线远场辐射功率数据而选定的模型参数,其中,参数k越大,天线波束越窄。
目前,工程上常用二次函数去近似地拟合圆口径天线的辐射功率数据,计算简单,但拟合效果不好,因为圆口径天线的辐射特性并不严格符合二次函数规律。因此,本实施例结合一阶贝塞尔函数构建圆口径天线的远场辐射参数模型,采用间接平差算法估计模型参数,使之更加精准地逼近实际的辐射特性。
步骤S103,确定所述待拟合参数的初值,并根据间接平差理论计算所述初值的改正量。
具体的,本实施例可以根据天线的基本参数对待拟合参数进行初始化,得到参数初值;基于参数初值进行参数估计,即本实施例采用间接平差算法确定参数初值的改正量,对参数初值进行更新,保证参数的精准性。
可选的,步骤S103中所述的确定所述待拟合参数的初值,包括:
通过
Figure BDA0002789165550000081
确定所述参数向量X的初值,得到向量初值X0;其中,D为所述天线的口径,λ为所述天线的工作波长,Pi为所述远场辐射功率数据集中的功率值。
本实施例根据天线的基本参数,例如根据口径和工作波长对待拟合参数进行初始化,使得初始参数更加接近最终拟合结果,降低计算量,提高收敛速度,保证参数拟合的准确、高效。
可选的,步骤S103中所述的根据间接平差理论计算所述初值的改正量的具体实现流程可以包括:
先,计算所述远场辐射参数模型对所述参数向量在所述向量初值X0处的Jacobi矩阵B:
Figure BDA0002789165550000082
其中,
Figure BDA0002789165550000083
Figure BDA0002789165550000084
其中,J2为二阶贝塞尔函数,I为所述远场辐射功率数据集的元素总数,α0为所述天线的视轴方向对应的角度的初值,k0表征天线波束宽度的因子的初值,C0为所述天线的远场辐射功率的最大值的初值。
然后,根据间接平差理论,通过
Figure BDA0002789165550000085
确定所述向量初值X0的改正量
Figure BDA0002789165550000086
其中,l为所述远场辐射功率数据集中的功率值与基于所述远场辐射参数模型计算的近似值之差,公式如下:
Figure BDA0002789165550000091
其中,P1,P2,...,PI均为所述远场辐射功率数据集中的功率值,θ12,...,θI分别为P1,P2,...,PI对应的所述扫描角度。
本实施例通过上述间接平差理论确定的参数初值的改正量更加精准,保证了参数拟合的准确性和精准性。
步骤S104,在所述改正量满足预设阈值条件时,根据所述改正量和所述初值确定参数拟合结果。
本实施例在初始化拟合参数之后,通过间接平差理论计算参数初值的改正量,并迭代求解直至收敛。具体的,本实施例确定改正量
Figure BDA0002789165550000092
的每个元素是否满足预设阈值条件,例如判断改正量
Figure BDA0002789165550000093
中每个元素的绝对值与门限Tol的大小,若
Figure BDA0002789165550000094
退出迭代,根据所述改正量参数初值确定参数拟合结果。例如,将
Figure BDA0002789165550000095
确定为参数的最终值,即参数拟合结果。可选的,本实施例可以根据收敛精度需求设定预设阈值条件,例如Tol=1×10-5
可选的,所述天线的远场辐射功率数据的参数拟合方法还可以包括:
在所述改正量不满足预设阈值条件时,取
Figure BDA0002789165550000096
作为新的初值X'0,X0为所述待拟合参数的初值。
根据间接平差理论计算所述新的初值X'0的改正量
Figure BDA0002789165550000097
在改正量
Figure BDA0002789165550000098
满足所述预设阈值条件时,根据
Figure BDA0002789165550000099
确定所述参数拟合结果
Figure BDA00027891655500000910
否则,迭代计算直至改正量
Figure BDA00027891655500000911
满足所述预设阈值条件。
具体的,取
Figure BDA00027891655500000912
作为新的向量初值X0',更新计算远场辐射参数模型对参数向量X在向量初值X0处的Jacobi矩阵B,以及远场辐射功率数据集中的功率值与基于远场辐射参数模型求得的近似值之差l,从而得到新的改正量
Figure BDA0002789165550000101
若不满足所述预设阈值条件
Figure BDA0002789165550000102
则重复上述步骤,直至满足后停止迭代,得到参数向量X的最终估计值:
Figure BDA0002789165550000103
其中,
Figure BDA0002789165550000104
的第1、2、3元素,即参数k、α和C的最终估计值。
在一个实施例中,所述天线的远场辐射功率数据的参数拟合方法还可以包括:
根据
Figure BDA0002789165550000105
确定所述参数拟合结果的中误差,完成精度评定,该矩阵对角线的第1、2、3元素,即为参数k、α和C估计值的中误差;其中,
Figure BDA0002789165550000106
为向量初值的改正量,I为远场辐射功率数据集的元素总数,B为远场辐射参数模型对参数向量X在向量初值X0处的Jacobi矩阵,l为远场辐射功率数据集中的功率值与基于远场辐射参数模型计算的近似值之差。
本发明针对圆口径天线远场辐射理论特性,结合工程实际,建立其辐射功率数据的参数模型,采用间接平差算法估计得到相关参数,并评定估计精度,拟合方法适应性强、估计准确、收敛迅速,可以较为准确地刻画圆口径天线远场辐射特性,适应于天线指向偏差标定或其他相关复杂计算场景对远场高精度辐射模型建立和求解的需求。
示例性的,以下结合某航天器的X频段(波长0.036m)、4.2m圆口径天线的远场功率(实际为增益)数据为例,对本发明的具体实施方式作进一步说明。
图2所示为该4.2m圆口径天线的远场增益数据,横轴为角度,纵轴为天线增益(dBi),共33点数据。图示中,增益数据对应角度的范围为-0.8°至+0.8°,最大增益为46.0511(dBi),位置大约在-0.1°(电轴偏差约-0.1°)。基于该数据,按照如下流程进行参数拟合。
Step1:从原测量数据中筛选符合条件的数据,记为{(θi,Pi)}(i=1,2,…,I),筛选条件如下:
Pi≥max(P0n)-10=36.0511。
经筛选所得的数据已在图1中用圆圈示出,共计24个,即I=24。
Step2:建立圆口径天线远场辐射参数模型,
Figure BDA0002789165550000111
远场辐射参数模型的三个参数k、α和C构成参数向量X:
X=[k α C]T
Step3:确定参数向量X的初值X0
Figure BDA0002789165550000112
Step4:计算测量模型对参数向量X在X0处的Jacobi矩阵B,即:
Figure BDA0002789165550000113
其中,
Figure BDA0002789165550000114
Figure BDA0002789165550000115
Step5:根据间接平差理论,利用如下公式,求解参数向量X的改正量
Figure BDA0002789165550000116
Figure BDA0002789165550000117
其中,l的计算公式如下:
Figure BDA0002789165550000118
Step6:参数估计,迭代求解直至收敛。根据改正量
Figure BDA0002789165550000121
的计算公式,第一次求解得到
Figure BDA0002789165550000122
设定收敛门限Tol为1×10-5,因
Figure BDA0002789165550000123
Figure BDA0002789165550000124
作为新的初值X0,更新计算Step4、Step5中的B、l,计算对应新的初值的改正量。重复上述判断和计算,表1所示为逐次迭代求解得到的改正量及其分量最大值。第4次迭代计算达到收敛条件,迭代结束,得到参数k、α和C的最终估计值:
k=4.4492
α=-0.0814。
C=46.0332
表1逐次迭代求解得到的改正量及其分量最大值
Figure BDA0002789165550000125
Step7:完成精度评定。按照下式计算所估参数的中误差:
Figure BDA0002789165550000126
该矩阵对角线的第1、2、3元素,为参数k、α和C估计值的中误差,如下:
1.0σk=0.0036
1.0σα=0.0003。
1.0σC=0.0074
图3为本实施例拟合方法对该天线远场增益数据的拟合效果的图示,黑色数据点(左轴)所示为用于参数拟合的远场辐射功率数据集,虚线(左轴)所示为拟合得到的曲线,“+”数据点(右轴)所示为拟合后残差情况。作为对比,图4所示为常用的二次函数拟合结果,可以看出:本实施例的拟合结果更加贴合实测数据,拟合后残差(中误差约0.03dBi)较二次函数拟合后残差(中误差约0.13dBi)显著减小。
上述天线的远场辐射功率数据的参数拟合方法中,主要针对圆口径天线的远场辐射理论特性,合理选定待拟合参数,建立其远场辐射功率的参数模型,采用间接平差算法进行参数估计,得到参数拟合结果。本发明的拟合方法适应性强、估计准确、收敛迅速,可以较为准确地刻画圆口径天线远场辐射特性,适应于天线指向偏差标定或其他相关复杂计算场景对远场高精度辐射模型建立和求解的需求。
本领域技术人员可以理解,上述实施例中各步骤的序号的大小并不意味着执行顺序的先后,各过程的执行顺序应以其功能和内在逻辑确定,而不应对本发明实施例的实施过程构成任何限定。
对应于上文实施例所述的天线的远场辐射功率数据的参数拟合方法,本实施例提供了一种天线的远场辐射功率数据的参数拟合装置。具体参见图5,为本实施例中天线的远场辐射功率数据的参数拟合装置的结构示意图。为了便于说明,仅示出了与本实施例相关的部分。
所述远场辐射功率数据的参数拟合装置主要包括:数据集确定模块110、辐射模型建立模块120、改正量计算模块130和拟合结果确定模块140。
数据集确定模块110用于确定天线的远场辐射功率数据集,所述远场辐射功率与所述天线的扫描角度一一对。
辐射模型建立模块120用于根据所述远场辐射功率数据集确定待拟合参数,并建立远场辐射参数模型。
改正量计算模块130用于确定所述待拟合参数的初值,并根据间接平差理论计算所述初值的改正量。
拟合结果确定模块140用于在所述改正量满足预设阈值条件时,根据所述改正量和所述初值确定参数拟合结果。
上述天线的远场辐射功率数据的参数拟合装置,主要针对圆口径天线的远场辐射理论特性,合理选定待拟合参数,并建立其远场辐射功率的参数模型,并采用间接平差算法进行参数估计,得到参数拟合结果,本发明的拟合方法适应性强、估计准确、收敛迅速,可以较为准确地刻画圆口径天线远场辐射特性,适应于天线指向偏差标定或其他相关复杂计算场景对远场高精度辐射模型建立和求解的需求。
本实施例还提供了一种天线的远场辐射功率数据的参数拟合装置100的示意图。如图6所示,该实施例的天线的远场辐射功率数据的参数拟合装置100包括:处理器150、存储器160以及存储在所述存储器160中并可在所述处理器150上运行的计算机程序161,例如天线的远场辐射功率数据的参数拟合方法的程序。
其中,处理器150在执行存储器160上所述计算机程序161时实现上述天线的远场辐射功率数据的参数拟合方法实施例中的步骤,例如图1所示的步骤101至103。或者,所述处理器150执行所述计算机程序161时实现上述各装置实施例中各模块/单元的功能,例如图5所示模块110至140的功能。
示例性的,所述计算机程序161可以被分割成一个或多个模块/单元,所述一个或者多个模块/单元被存储在所述存储器160中,并由所述处理器150执行,以完成本发明。所述一个或多个模块/单元可以是能够完成特定功能的一系列计算机程序指令段,该指令段用于描述所述计算机程序161在所述远场辐射功率数据的参数拟合装置100中的执行过程。例如,所述计算机程序161可以被分割成数据集确定模块110、辐射模型建立模块120、改正量计算模块130和拟合结果确定模块140,各模块具体功能如下:
数据集确定模块110用于确定天线的远场辐射功率数据集,所述远场辐射功率与所述天线的扫描角度一一对。
辐射模型建立模块120用于根据所述远场辐射功率数据集确定待拟合参数,并建立远场辐射参数模型。
改正量计算模块130用于确定所述待拟合参数的初值,并根据间接平差理论计算所述初值的改正量。
拟合结果确定模块140用于在所述改正量满足预设阈值条件时,根据所述改正量和所述初值确定参数拟合结果。
所述远场辐射功率数据的参数拟合装置100可包括,但不仅限于处理器150、存储器160。本领域技术人员可以理解,图6仅仅是天线的远场辐射功率数据的参数拟合装置100的示例,并不构成对天线的远场辐射功率数据的参数拟合装置100的限定,可以包括比图示更多或更少的部件,或者组合某些部件,或者不同的部件,例如所述远场辐射功率数据的参数拟合装置100还可以包括输入输出设备、网络接入设备、总线等。
所述处理器150可以是中央处理单元(Central Processing Unit,CPU),还可以是其他通用处理器、数字信号处理器(Digital Signal Processor,DSP)、专用集成电路(Application Specific Integrated Circuit,ASIC)、现成可编程门阵列(Field-Programmable Gate Array,FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。
所述存储器160可以是所述远场辐射功率数据的参数拟合装置100的内部存储单元,例如天线的远场辐射功率数据的参数拟合装置100的硬盘或内存。所述存储器160也可以是所述远场辐射功率数据的参数拟合装置100的外部存储设备,例如所述远场辐射功率数据的参数拟合装置100上配备的插接式硬盘,智能存储卡(Smart Media Card,SMC),安全数字(Secure Digital,SD)卡,闪存卡(Flash Card)等。进一步地,所述存储器160还可以既包括所述远场辐射功率数据的参数拟合装置100的内部存储单元也包括外部存储设备。所述存储器160用于存储所述计算机程序以及所述远场辐射功率数据的参数拟合装置100所需的其他程序和数据。所述存储器160还可以用于暂时地存储已经输出或者将要输出的数据。
所属领域的技术人员可以清楚地了解到,为了描述的方便和简洁,仅以上述各功能单元、模型的划分进行举例说明,实际应用中,可以根据需要而将上述功能分配由不同的功能单元、模块完成,即将所述装置的内部结构划分成不同的功能单元或模块,以完成以上描述的全部或者部分功能。实施例中的各功能单元、模块可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中,上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。另外,各功能单元、模块的具体名称也只是为了便于相互区分,并不用于限制本申请的保护范围。上述系统中单元、模块的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。
另外,在本发明各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。
所述集成的模块/单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明实现上述实施例方法中的全部或部分流程,也可以通过计算机程序来指令相关的硬件来完成,所述的计算机程序可存储于一计算机可读存储介质中,该计算机程序在被处理器执行时,可实现上述各个方法实施例的步骤。其中,所述计算机程序包括计算机程序代码,所述计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。所述计算机可读介质可以包括:能够携带所述计算机程序代码的任何实体或装置、记录介质、U盘、移动硬盘、磁碟、光盘、计算机存储器、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、电载波信号、电信信号以及软件分发介质等。需要说明的是,所述计算机可读介质包含的内容可以根据司法管辖区内立法和专利实践的要求进行适当的增减,例如在某些司法管辖区,根据立法和专利实践,计算机可读介质不包括是电载波信号和电信信号。
以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围,均应包含在本发明的保护范围之内。

Claims (10)

1.一种天线的远场辐射功率数据的参数拟合方法,其特征在于,包括:
确定天线的远场辐射功率数据集,所述远场辐射功率与所述天线的扫描角度一一对应;
根据所述远场辐射功率数据集确定待拟合参数,并建立远场辐射参数模型;
确定所述待拟合参数的初值,并根据间接平差理论计算所述初值的改正量;
在所述改正量满足预设阈值条件时,根据所述改正量和所述初值确定参数拟合结果。
2.如权利要求1所述的天线的远场辐射功率数据的参数拟合方法,其特征在于,所述确定天线的远场辐射功率数据集,包括:
获取天线在所述扫描角度θ0n处的远场辐射功率的测量值P0n,n=1,2,…,N;
根据
{(θi,Pi)|(θi,Pi)∈{(θ0n,P0n)},Pi≥max(P0n)-10}
确定所述远场辐射功率数据集;其中,i=1,2,…,I,I≤N。
3.如权利要求1所述的天线的远场辐射功率数据的参数拟合方法,其特征在于,所述根据所述远场辐射功率数据集确定待拟合参数,并建立远场辐射参数模型,包括:
确定待拟合参数k、α和C,所述远场辐射功率数据集中的功率值Pi与对应的扫描角度θi满足下式:
Figure FDA0002789165540000011
其中,vi为测量和模型误差,J1为一阶贝塞尔函数,k为表征天线波束宽度的因子,α为所述天线的视轴方向对应的角度,C为所述天线的远场辐射功率的最大值;所述待拟合参数构成参数向量X=[k,α,C]T
4.如权利要求3所述的天线的远场辐射功率数据的参数拟合方法,其特征在于,所述确定所述待拟合参数的初值,包括:
通过
Figure FDA0002789165540000021
确定所述参数向量X的初值,得到向量初值X0;其中,D为所述天线的口径,λ为所述天线的工作波长,Pi为所述远场辐射功率数据集中的功率值。
5.如权利要求4所述的天线的远场辐射功率数据的参数拟合方法,其特征在于,所述根据间接平差理论计算所述初值的改正量,包括:
计算所述远场辐射参数模型对所述参数向量在所述向量初值X0处的Jacobi矩阵B:
Figure FDA0002789165540000022
其中,
Figure FDA0002789165540000023
Figure FDA0002789165540000024
其中,J2为二阶贝塞尔函数;
根据间接平差理论,通过
Figure FDA0002789165540000025
确定所述向量初值的改正量
Figure FDA0002789165540000026
其中,l为所述远场辐射功率数据集中的功率值与基于所述远场辐射参数模型计算的近似值之差,公式如下:
Figure FDA0002789165540000027
其中,P1,P2,...,PI均为所述远场辐射功率数据集中的功率值,θ12,...,θI分别为P1,P2,...,PI对应的所述扫描角度。
6.如权利要求1所述的天线的远场辐射功率数据的参数拟合方法,其特征在于,所述天线的远场辐射功率数据的参数拟合方法还包括:
在所述改正量不满足预设阈值条件时,取
Figure FDA0002789165540000031
作为新的初值X′0,X0为所述待拟合参数的初值;
根据间接平差理论计算所述新的初值X′0的改正量
Figure FDA0002789165540000032
在改正量
Figure FDA0002789165540000033
满足所述预设阈值条件时,根据
Figure FDA0002789165540000034
确定所述参数拟合结果
Figure FDA0002789165540000035
否则,迭代计算直至改正量
Figure FDA0002789165540000036
满足所述预设阈值条件。
7.如权利要求5所述的天线的远场辐射功率数据的参数拟合方法,其特征在于,所述天线的远场辐射功率数据的参数拟合方法还包括:
根据
Figure FDA0002789165540000037
确定所述参数拟合结果的中误差,其中,
Figure FDA0002789165540000038
为所述向量初值的改正量,I为所述远场辐射功率数据集的元素总数,B为所述远场辐射参数模型对所述参数向量X在向量初值X0处的Jacobi矩阵,l为所述远场辐射功率数据集中的功率值与基于所述远场辐射参数模型计算的近似值之差。
8.一种天线的远场辐射功率数据的参数拟合装置,其特征在于,包括:
数据集确定模块,用于确定天线的远场辐射功率数据集,所述远场辐射功率与所述天线的扫描角度一一对应;
辐射模型建立模块,用于根据所述远场辐射功率数据集确定待拟合参数,并建立远场辐射参数模型;
改正量计算模块,用于确定所述待拟合参数的初值,并根据间接平差理论计算所述初值的改正量;
拟合结果确定模块,用于在所述改正量满足预设阈值条件时,根据所述改正量和所述初值确定参数拟合结果。
9.一种天线的远场辐射功率数据的参数拟合装置,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现如权利要求1至7任一项所述的天线的远场辐射功率数据的参数拟合方法的步骤。
10.一种计算机可读存储介质,所述计算机可读存储介质存储计算机程序,其特征在于,所述计算机程序被处理器执行时实现如权利要求1至7任一项所述的天线的远场辐射功率数据的参数拟合方法的步骤。
CN202011308932.4A 2020-11-20 2020-11-20 天线的远场辐射功率数据的参数拟合方法和装置 Active CN112417682B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011308932.4A CN112417682B (zh) 2020-11-20 2020-11-20 天线的远场辐射功率数据的参数拟合方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011308932.4A CN112417682B (zh) 2020-11-20 2020-11-20 天线的远场辐射功率数据的参数拟合方法和装置

Publications (2)

Publication Number Publication Date
CN112417682A true CN112417682A (zh) 2021-02-26
CN112417682B CN112417682B (zh) 2022-10-28

Family

ID=74773246

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011308932.4A Active CN112417682B (zh) 2020-11-20 2020-11-20 天线的远场辐射功率数据的参数拟合方法和装置

Country Status (1)

Country Link
CN (1) CN112417682B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113533867A (zh) * 2021-07-14 2021-10-22 西安电子科技大学 基于Fourier插值的远场方向图快速测量方法
CN116106644A (zh) * 2023-04-10 2023-05-12 中国人民解放军63921部队 计及角度误差的天线远场辐射功率拟合方法、系统及介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107239602A (zh) * 2017-05-22 2017-10-10 中国电子科技集团公司第四十研究所 一种基于曲线拟合的探头天线模型快速计算方法
CN109631945A (zh) * 2018-12-14 2019-04-16 北京遥测技术研究所 一种遥测地面站指向标校方法
US20190179005A1 (en) * 2017-12-08 2019-06-13 Rohde & Schwarz Gmbh & Co. Kg Near-field measurement system and method for obtaining far-field characteristics
CN110222299A (zh) * 2019-05-05 2019-09-10 中国人民解放军63921部队 针对双变量含有误差的直线拟合问题的方法和设备

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107239602A (zh) * 2017-05-22 2017-10-10 中国电子科技集团公司第四十研究所 一种基于曲线拟合的探头天线模型快速计算方法
US20190179005A1 (en) * 2017-12-08 2019-06-13 Rohde & Schwarz Gmbh & Co. Kg Near-field measurement system and method for obtaining far-field characteristics
CN109631945A (zh) * 2018-12-14 2019-04-16 北京遥测技术研究所 一种遥测地面站指向标校方法
CN110222299A (zh) * 2019-05-05 2019-09-10 中国人民解放军63921部队 针对双变量含有误差的直线拟合问题的方法和设备

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
PEIYUAN LIAN ET AL.: "Feed adjustment method of reflector antenna based on far field", 《IET MICROWAVES, ANTENNAS & PROPAGATION》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113533867A (zh) * 2021-07-14 2021-10-22 西安电子科技大学 基于Fourier插值的远场方向图快速测量方法
CN113533867B (zh) * 2021-07-14 2022-09-06 西安电子科技大学 基于Fourier插值的远场方向图快速测量方法
CN116106644A (zh) * 2023-04-10 2023-05-12 中国人民解放军63921部队 计及角度误差的天线远场辐射功率拟合方法、系统及介质
CN116106644B (zh) * 2023-04-10 2023-06-27 中国人民解放军63921部队 计及角度误差的天线远场辐射功率拟合方法、系统及介质

Also Published As

Publication number Publication date
CN112417682B (zh) 2022-10-28

Similar Documents

Publication Publication Date Title
CN110335297B (zh) 一种基于特征提取的点云配准方法
CN112417682B (zh) 天线的远场辐射功率数据的参数拟合方法和装置
CN109472059A (zh) 基于测量应变的相控阵天线幅相补偿方法
CN111174753B (zh) 基于有理函数模型的光学影像与激光测高数据平差方法
CN110907904B (zh) 一种太阳法方位波束宽度修正方法
CN111046591B (zh) 传感器幅相误差与目标到达角度的联合估计方法
CN112362080B (zh) 航天器天线在轨指向标定的器地数据同步偏差确定方法
CN112444783A (zh) 基于自然均匀目标的合成孔径雷达在轨快速定标方法
CN110806578A (zh) 一种波束控制方法及装置、可读存储介质
CN109444786A (zh) 提高在片负载牵引测量准确度的方法、系统及终端设备
JP6739126B1 (ja) 計測データ処理装置
CN107239602B (zh) 一种基于曲线拟合的探头天线模型快速计算方法
CN112986701A (zh) 一种基于射电源宽带信号的全息测量方法及系统
CN110445557A (zh) 深空测控干涉测量大口径天线指向标校方法及装置
CN111505394A (zh) 基于探头天线位置误差修正的天线球面近场测量方法
CN101873436A (zh) 使用多线内插的镜头亮度补偿方法
CN114485731B (zh) 航天器双轴抛物面天线在轨指向标定方法及设备
CN107015065B (zh) 窄波束天线电轴、相位中心和时延的远场联合标定方法
CN111123250B (zh) 基于模式搜索算法的脉冲多普勒雷达和波束形成方法
CN110907705A (zh) 一种天线增益的测量方法、装置、测量系统、计算机设备以及存储介质
CN116106644B (zh) 计及角度误差的天线远场辐射功率拟合方法、系统及介质
CN107367235B (zh) 一种红外面阵扫描系统的轴系误差标定方法
CN114234904B (zh) 一种基于透镜多波束天线的测角方法、系统及介质
CN117572326B (zh) 微波遥感仪器的辐射定标方法、装置及计算机存储介质
CN112199814B (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