CN106053968B - 使用时频分析测量的单道天线远场天线因子估计方法 - Google Patents

使用时频分析测量的单道天线远场天线因子估计方法 Download PDF

Info

Publication number
CN106053968B
CN106053968B CN201610388610.2A CN201610388610A CN106053968B CN 106053968 B CN106053968 B CN 106053968B CN 201610388610 A CN201610388610 A CN 201610388610A CN 106053968 B CN106053968 B CN 106053968B
Authority
CN
China
Prior art keywords
antenna
formula
ground
far field
ref
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
CN201610388610.2A
Other languages
English (en)
Other versions
CN106053968A (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.)
CASIC Space Engineering Development Co Ltd
Original Assignee
CASIC Space Engineering Development Co Ltd
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 CASIC Space Engineering Development Co Ltd filed Critical CASIC Space Engineering Development Co Ltd
Priority to CN201610388610.2A priority Critical patent/CN106053968B/zh
Priority to PCT/CN2016/093510 priority patent/WO2017206320A1/zh
Publication of CN106053968A publication Critical patent/CN106053968A/zh
Application granted granted Critical
Publication of CN106053968B publication Critical patent/CN106053968B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R29/00Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
    • G01R29/08Measuring electromagnetic field characteristics
    • G01R29/10Radiation diagrams of antennas

Landscapes

  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Variable-Direction Aerials And Aerial Arrays (AREA)

Abstract

本发明提出一种使用时频分析测量的单通道天线远场天线参数估计方法,该方法应用于对数周期偶极数组天线,天线之间的间距为10米,该方法应用时间频率分析方法从而有效抑制一些该单一天线方法中出现的非期望反射,并且利用由时间频率方法所估计的天线距离估计远场天线因子,另外,使用天线辐射中心距离修正共识来估计对数周其偶极数组天线的远场天线因子,与其他方法所估计的远场天线因子相比,所估计的天线因子具有更好的远场增益相容性,天线距离为5米的地方,在频率范围从300兆赫兹到1000兆赫兹的范围内频率差异低于0.3dB。

Description

使用时频分析测量的单道天线远场天线因子估计方法
技术领域
本发明涉及一种天线因子估计的方法,特别是采用时间-频率方法测量的针对单道远场天线的因子估计方法。
背景技术
目前的远场天线因子测量方法存在如下的缺陷,首先不能使用时间域的减法将地面反射波和待测天线的固有自反射波分开测量;另外,之前所提出的测量方法是用固定的天线具有,由此,远场天线因子估计的主要的测量不确定度大部分和天线的距离有关,因此不能降低测量的不确定度;第三,不能确定每个频率下的天线距离,从而不能有效地确定估计方法来估计远场天线因子;第四,不能有效利用时域分析方法和时频分析方法,并使用其中的短时傅立叶变换技术进行天线因子的确定,从而估计精度不高。
由此提出一种使用时频分析测量的单通道天线远场天线参数估计方法,该方法应用于对数周期偶极数组天线,天线之间的间距为10米,该方法应用时间频率分析方法从而有效抑制一些该单一天线方法中出现的非期望反射,并且利用由时间频率方法所估计的天线距离估计远场天线因子,另外,使用天线辐射中心距离修正共识来估计对数周其偶极数组天线的远场天线因子,与其他方法所估计的远场天线因子相比,所估计的天线因子具有更好的远场增益相容性,天线距离为5米的地方,在频率范围从300兆赫兹到1000兆赫兹的范围内频率差异低于0.3dB。
发明内容
本发明的目的通过如下方案实现:
一种使用时频分析测量的单通道天线远场天线因子估计方法,该方法用于标定对数周期偶极数组天线,包括如下步骤:
(1)在消声室内的地面上设置所有的测量对象,将对数周期偶极距天线阵列作为被测对象,天线垂直于地面并位于地面之上的间距为1米到5米之间的范围;
(2)通过矢量网络分析仪测量频域的反射系数s11(ω),并且将其转换为相应的时域反射系数s11(t);
(3)通过公式(1)计算s11(t),其中公式使用傅立叶逆变换F-1(ω)以及汉宁频域窗W(ω),其中公式(1)为:
s11(t)=F-1(W(ω)s11(ω)) (1)
其中s11(t)包括s11intrinsic(t),s11ground_ref(t)以及其他非期望的波形,其中s11intrinsic(t)为天线内部结构的固有反射,s11ground_ref(t)为地面的垂直入射反射反射波,其他非期望反射波形具有非期望的环境反射波,如天线罩等;
(4)拾取垂直地面入射的反射波从而估计自由空间天线因子,采用将s11(t)平均的方式估计s11intrinsic(t),随着垂直对天线进行扫描将s11(t)平均,保留不变的部分s11intrinsic(t),并将其从s11(t)提取出;
(5)通过公式(2)将s11intrinsic(t)从时域的s11(t)中减去可以得到直接反射波s11ground_ref(t),公式(2)如下:
s11ground_ref(t)=s11(t)-s11intrinsic(t) (2);
(6)使用短时傅立叶变换以及汉宁时间窗的方法计算s11ground_ref(t)的时间-频率响应,从而确定每个频率的天线距离D=z(ω),其中,汉宁窗时域宽度为20ns,傅立叶变换的次数为8196;
(7)使用公式(3)由估计的天线距离估计天线因子:
af2(ω,z)=η0·[Z0·λ0·D(ω)·s11ground_ref(ω,z)]-1 (3)
其中,η0=120πΩ为自由空间特征阻抗,Z0=50Ω为同轴线缆的特征阻抗,并且λ0为相应的自由空间波长,D(ω)为待测对象幅射点到地面的双向距离,在对数周期偶极数组天线为测量对象的情况下,该双向距离的公式如公式(4)所示:
D(ω)=z+2×d1-f-λ (4)
其中,z为从待测对象的顶部到地面的双向距离,d1-f为天线顶部到对数周期偶极数组天线辐射点的距离;
(8)将天线的幅值中心位置应用于公式(4)获得远场天线因子估计值的公式(5):
其中,affar(ω)为远场天线因子。
附图说明
附图1为根据本发明实施例使用天线子标定方法建立天线因子测量值的示意图;
附图2为根据本发明实施例的对数周期偶极数组天线s11(t)的标定结果示意图;
附图3为根据本发明实施例的对数周期偶极数组天线的s11_ground_ref(t)估计值示意图;
附图4为采用本发明实施例方法估计的远场天线因子与使用其它方法估计的天线因子之间的差异值示意图。
附图标记:1-待测天线 2-地面 3-网络分析仪
具体实施方式
一种使用时频分析测量的单通道天线远场天线因子估计方法,首先该方法用于标定对数周期偶极数组天线,其实验建立方法示意性的如图1所示。所有的测量都是建立在消声室内的地面上,将对数周期偶极距天线阵列作为被测对象,天线垂直于地面并位于地面之上的间距为1米到5米之间的范围;通过矢量网络分析仪测量频域的反射系数s11(ω),并且将其转换为相应的时域反射系数s11(t);通过公式(1)计算s11(t),其中公式使用傅立叶逆变换F-1(ω)以及汉宁频域窗W(ω),公式(1)如下所示:
s11(t)=F-1(W(ω)s11(ω)) (1)
其中s11(t)包括s11intrinsic(t),s11ground_ref(t)以及其他非期望的波形,其中s11intrinsic(t)为天线内部结构的固有反射,s11ground_ref(t)为地面的垂直入射反射反射波,其他非期望反射波形具有非期望的环境反射波,如天线罩等,所估计的时域波形如图2所示。为了估计自由空间天线因子,我们必须拾取垂直地面入射的反射波,首先采用将s11(t)平均的方式估计s11intrinsic(t),这是因为s11intrinsic(t))的相位不会变化,即使距离地面的位置高度变化了,而相反的是,随着相对地面的双向传播,s11ground_ref(t)的相位会发生变化,随着垂直对天线进行扫描将s11(t)平均,保留不变的部分s11intrinsic(t),并将其从s11(t)提取出,也就是说,在时间轴上以固定的延迟时间出现不变的s11intrinsic(t)。
通过公式(2)将s11intrinsic(t)从时域的s11(t)中减去可以得到直接反射波s11ground_ref(t),公式(2)如下所示:
s11ground_ref(t)=s11(t)-s11intrinsic(t) (2);
所估计的s11ground_ref(t)如图3所示。
使用短时傅立叶变换以及汉宁时间窗的方法计算s11ground_ref(t)的时间-频率响应,从而确定每个频率的天线距离D=z(ω),其中,汉宁窗时域宽度为20ns,傅立叶变换的次数为8196;使用公式(3)由估计的天线距离估计天线因子,其中公式(3)如下所示:
af2(ω,z)=η0·[Z0·λ0·D(ω)·s11ground_ref(ω,z)]-1 (3)
其中,η0=120πΩ为自由空间特征阻抗,Z0=50Ω为同轴线缆的特征阻抗,并且λ0为相应的自由空间波长,D(ω)为待测对象幅射点到地面的双向距离,在对数周期偶极数组天线为测量对象的情况下,该双向距离的公式如公式(4)所示:
D(ω)=z+2×d1-f-λ (4)
其中,z为从待测对象的顶部到地面的双向距离,d1-f为天线顶部到对数周期偶极数组天线辐射点的距离;
将天线的幅值中心位置应用于公式(4)获得远场天线因子估计值的计算公式(5):
其中,affar(ω)为远场天线因子。
附图4表示使用本发明所述方法所估计的天线因子和相应的其他方法估计的远场天线因子之间的差值。三种测量情况分别是固定长度,短时傅立叶变换,以及短时傅立叶变换-单元长度。结果发现,采用本发明的方法,天线距离为5米的地方,在频率范围从300兆赫兹到1000兆赫兹的范围内频率差异低于0.3dB。

Claims (1)

1.使用时频分析测量的单通道天线远场天线因子估计方法,该方法用于标定对数周期偶极数组天线,其特征在于包括如下步骤:
(1)在消声室内的地面上设置所有的测量对象,将对数周期偶极数组天线作为被测对象,天线垂直于地面并位于地面之上的间距为1米到5米之间的范围;
(2)通过矢量网络分析仪测量频域的反射系数S11(ω),并且将其转换为相应的时域反射系数S11(t);
(3)通过公式(1)计算S11(t),其中公式使用傅立叶逆变换F-1(ω)以及汉宁频域窗W(ω),其中公式(1)为:
S11(t)=F-1(W(ω)s11(ω))(1)
其中S11(t)包括S11intrinsic(t),S11ground_ref(t)以及其他非期望的反射波形,其中S11intrinsic(t)为天线内部结构的固有反射,S11ground_ref(t)为垂直地面入射的反射波,其他非期望反射波形具有非期望的环境反射波;
(4)拾取垂直地面入射的反射波S11ground_ref(t)从而估计自由空间天线因子,采用将S11(t)平均的方式估计S11intrinsic(t);S11intrinsic(t)的相位不会变化,即使距离地面的位置高度变化了,而相反的是,随着相对地面的双向传播,S11ground_ref(t)的相位会发生变化,随着对天线进行垂直扫描从而将S11(t)平均,保留不变的部分为S11intrinsic(t),并将其从S11(t)提取出;即在时间轴上以固定的延迟时间出现不变的S11intrinsic(t);
(5)通过公式(2)将S11intrinsic(t)从时域的S11(t)中减去可以得到垂直地面入射的反射波S11ground_ref(t),公式(2)如下:
S11ground_ref(t)=S11(t)-S11intrinsic(t)(2);
(6)使用短时傅立叶变换以及汉宁时间窗的方法计算S11ground_ref(t)的时间-频率响应,从而确定每个频率的天线距离D(ω)=f(z),其中,汉宁窗时域宽度为20ns,傅立叶变换的次数为8196;
(7)使用公式(3)由估计的天线距离估计天线因子:
af2(ω,z)=η0·[Z0·λ0·D(ω)·S11ground_ref(ω,z)]-1(3)
其中,η0=120πΩ为自由空间特征阻抗,Z0=50Ω为同轴线缆的特征阻抗,并且λ0为相应的自由空间波长,D(ω)为待测对象幅射点到地面的双向距离,在对数周期偶极数组天线为测量对象的情况下,该双向距离的公式如公式(4)所示:
D(ω)=z+2×d1-f-λ(4)
其中,z为从待测对象的顶部到地面的双向距离,d1-f为天线顶部到对数周期偶极数组天线辐射点的距离;
(8)将天线的幅值中心位置应用于公式(4)获得远场天线因子估计值的公式(5):
其中,affar(ω)为远场天线因子。
CN201610388610.2A 2016-06-03 2016-06-03 使用时频分析测量的单道天线远场天线因子估计方法 Active CN106053968B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201610388610.2A CN106053968B (zh) 2016-06-03 2016-06-03 使用时频分析测量的单道天线远场天线因子估计方法
PCT/CN2016/093510 WO2017206320A1 (zh) 2016-06-03 2016-08-05 使用时频分析测量的单道天线远场天线因子估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610388610.2A CN106053968B (zh) 2016-06-03 2016-06-03 使用时频分析测量的单道天线远场天线因子估计方法

Publications (2)

Publication Number Publication Date
CN106053968A CN106053968A (zh) 2016-10-26
CN106053968B true CN106053968B (zh) 2019-09-10

Family

ID=57170269

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610388610.2A Active CN106053968B (zh) 2016-06-03 2016-06-03 使用时频分析测量的单道天线远场天线因子估计方法

Country Status (2)

Country Link
CN (1) CN106053968B (zh)
WO (1) WO2017206320A1 (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112595902B (zh) * 2020-12-24 2022-07-15 中国科学技术大学 一种阵列天线系统故障诊断方法及装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5204685A (en) * 1992-07-01 1993-04-20 The United States Of America As Represented By The Secretary Of The Air Force ARC range test facility
CN1744379A (zh) * 2005-09-28 2006-03-08 中山大学 多端口控制频率固定波束可调的微带漏波天线
CN103245841A (zh) * 2013-04-28 2013-08-14 西北工业大学 一种基于球面近场扫描外推的天线方向图测试方法
CN104931799A (zh) * 2015-04-30 2015-09-23 中国电子科技集团公司第四十一研究所 一种在片天线的电性能测试系统及方法
CN105445566A (zh) * 2015-11-13 2016-03-30 西北工业大学 球面多探头天线测试数据处理方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1466183A1 (en) * 2002-01-18 2004-10-13 Her Majesty in Right of Canada as Represented by the Minister of Industry Antenna array for the measurement of complex electromagnetic fields
JP2004354362A (ja) * 2003-05-29 2004-12-16 Taiyo Musen Co Ltd アンテナの遠方界放射パターン測定法
CN101308177B (zh) * 2008-07-11 2010-08-25 西安电子科技大学 主动反射面天线的电性能预测方法
CN102818942B (zh) * 2012-08-24 2014-08-20 湖北航天技术研究院计量测试技术研究所 天线远场参数校准装置及校准方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5204685A (en) * 1992-07-01 1993-04-20 The United States Of America As Represented By The Secretary Of The Air Force ARC range test facility
CN1744379A (zh) * 2005-09-28 2006-03-08 中山大学 多端口控制频率固定波束可调的微带漏波天线
CN103245841A (zh) * 2013-04-28 2013-08-14 西北工业大学 一种基于球面近场扫描外推的天线方向图测试方法
CN104931799A (zh) * 2015-04-30 2015-09-23 中国电子科技集团公司第四十一研究所 一种在片天线的电性能测试系统及方法
CN105445566A (zh) * 2015-11-13 2016-03-30 西北工业大学 球面多探头天线测试数据处理方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
FDTD中近远场变换计算天线方向图的新方法;吴锋涛、尹家贤、袁乃昌;《微波学报》;20050430;第21卷(第02期);全文

Also Published As

Publication number Publication date
WO2017206320A1 (zh) 2017-12-07
CN106053968A (zh) 2016-10-26

Similar Documents

Publication Publication Date Title
US8102298B2 (en) Ultra-wideband radar waveform calibration for measurements of a heterogeneous material
Martorella et al. On the fractal dimension of sea surface backscattered signal at low grazing angle
JP2010535344A (ja) パラメトリック時間領域手法を用いた地上反射波軽減
CN109799495A (zh) 一种用于高保真阵列处理的宽带时延估计方法
CN109471095A (zh) 基于快速迭代插值的fmcw雷达距离估计算法
Stock et al. Multiple baseline lightning interferometry-Improving the detection of low amplitude VHF sources
CN106053968B (zh) 使用时频分析测量的单道天线远场天线因子估计方法
Cosoli et al. A real-time and offline quality control methodology for SeaSonde high-frequency radar currents
Shvets et al. Estimating the lower ionosphere height and lightning location using multimode “tweek” atmospherics
Wyatt et al. Measuring high and low waves with HF radar
Costa et al. Adjustment of state space models in view of area rainfall estimation
Bialer et al. A time-of-arrival estimation algorithm for OFDM signals in indoor multipath environments
Lustmann et al. Indoor channel spectral statistics, K-factor and reverberation distance
CN107315169B (zh) 基于二阶统计量相似度的杂波协方差矩阵估计方法
de Paolo et al. Improving SeaSonde radial velocity accuracy and variance using radial metrics
Eisen et al. Identifying isochrones in GPR profiles from DEP-based forward modeling
Sinatra et al. Calibration of spatial rain scanner using rainfall depth of rain gauges
Alonso-Arroyo et al. Two dedicated soil moisture experiments using the scatterometric properties of GNSS-reflectometry
Adhyapak et al. Free space antenna factor computation using time domain gating and deconvolution filter for site validation of fully anechoic rooms
Kurokawa et al. Antenna gain and phase center estimation for broadband antenna using metal ground plane
Chang et al. Polarization measurements in a UWB multipath channel
Santos et al. Scatterer detection by successive cancellation for UWB-method and experimental verification
Wagenaars et al. Algorithms for arrival time estimation of partial discharge pulses in cable systems
Mertens et al. Determination of the stability of a pulse GPR system and quantification of the drift effect on soil material characterization by full-wave inversion
Shvets et al. A Technique for Automatic Monitoring the Lower Ionosphere and Lightning Location by Tweekatmospherics

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20190821

Address after: 430416 Xuelinlute No. 1, Shuangliu Street, Xinzhou District, Wuhan City, Hubei Province

Applicant after: Aerospace Science and Engineering Space Engineering Development Co., Ltd.

Address before: 518000 Guangdong city of Shenzhen province Nanshan District two road, Taoyuan Street dragon dragons Mingyuan 4 Building 405

Applicant before: Shenzhen Fanxi Electronics Co., Ltd.