CN113341238B - 一种利用太阳辐射测量天线方向图的方法 - Google Patents
一种利用太阳辐射测量天线方向图的方法 Download PDFInfo
- Publication number
- CN113341238B CN113341238B CN202110635512.5A CN202110635512A CN113341238B CN 113341238 B CN113341238 B CN 113341238B CN 202110635512 A CN202110635512 A CN 202110635512A CN 113341238 B CN113341238 B CN 113341238B
- Authority
- CN
- China
- Prior art keywords
- antenna
- formula
- sun
- angle
- solar radiation
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R29/00—Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
- G01R29/08—Measuring electromagnetic field characteristics
- G01R29/10—Radiation diagrams of antennas
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- 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)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Evolutionary Computation (AREA)
- Computing Systems (AREA)
- Geometry (AREA)
- Algebra (AREA)
- Computer Hardware Design (AREA)
- Electromagnetism (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明公开了一种利用太阳辐射测量天线方向图的方法,其利用太阳作为辐射源即可实现天线方向图测试,解决了在微波暗室测试天线方向图具有局限性的问题;且只需在晴天条件下利用系统自身观测太阳辐射强度,无需其他测试仪器仪表和特殊测试环境就可实现天线方向图测量。因此本发明利用太阳辐射法测量微波天线方向图的技术避免了测试过程对微波暗室特殊环境的依赖,通过自动跟踪扫描太阳,就可直接测量出天线方向性函数,能够在外场环境实现,降低了天线测试成本和周期,提高了测试效率。
Description
【技术领域】
本发明属于微波天线测试领域,涉及一种测量天线方向图的方法,具体涉及一种利用太阳辐射测量天线方向图的方法。
【背景技术】
传统测试天线方向图的方法是在微波暗室内利用标准增益天线和辐射源等实现天线方向图的测试,该测试过程复杂,需要精密的测试仪器,费用高,可重复性差,且需要特殊的测量环境,无法在野外等环境下实现。因此,需要考虑一种测量方法较为简单,对环境要求不高,易于实现且成本低的天线测试方法。
【发明内容】
针对上述问题,本发明提供一种利用太阳辐射测量天线方向图的方法,不仅避免了对微波暗室的依赖,而且可在外场环境易于实现,能够实现自动测量,大幅度降低了产品的测试成本,具有测量方法简单、适用性广、成本低等优点。
本发明是通过以下技术方案实现的,提供一种利用太阳辐射测量天线方向图的方法,该方法利用太阳作为辐射源进行天线方向图的测量,实现在野外环境即可实时测量天线方向图而不依赖传统的微波暗室天线方向图的测量方法,具体包括以下步骤:
S1实时计算太阳方位角和高度角;
S2控制转台,使天线指向太阳,所述转台为方位和俯仰转台;
S3根据计算获得天线在太阳中心周围扫描的亮温散点图;
S4移动天线方位,使得太阳不在天线波束内,控制天线俯仰对天空进行扫描,计算各个仰角的大气衰减;
S5校准大气衰减,计算没有大气衰减时观测到的太阳辐射到达天线的亮温;
S6利用高斯辐射模型拟合得到天线方向性函数及天线波束宽度;
S7校准计算天线增益;
S8计算天线有效接收面积和孔径效率。
特别的,所述S1根据日地轨道关系,计算太阳方位角Az和高度角El,具体按照如下公式进行计算:
T0=(ts-12)·15° (1),
El=arcsin(sinθlatsinδ+cosδcosT0) (2),
于公式(1)中,ts为太阳时角,于公式(2)和公式(3)中,δ为太阳赤纬,θlat为微波天线所在纬度。
特别的,所述S3具体通过以下方法实现:
S31固定转台的天线仰角,通过步进方式控制转台天线指向太阳中心位置,控制方位转台在角度-10°~10°转动天线扫描太阳,观测各频点太阳辐射强度,所述太阳辐射强度为经过大气衰减到达天线的信号和大气自身辐射信号;
S32改变转台的天线仰角,重复S31,天线的方位和仰角步进范围在太阳中心-10°~10°范围内,实现对太阳中心周围的栅格扫描;
S33按如下公式计算天线观测的亮温:
于公式(4)中,所述为天线方位角,θ为俯仰角,Ωs为太阳波束立体角,ΩA为天线波束立体角,Tm为大气平均辐射温度,Tbg=2.75K为宇宙背景辐射亮温,Tsun为太阳平均辐射亮温,τ(θ)为天线指定角度的大气厚度。
特别的,所述S4中观测天空时的亮温按如下公式计算:
太阳辐射到达天线被接收的亮温为公式(4)和公式(5)相减,整理后得到如下公式:
按如下公式计算大气衰减,即天线指定角度的大气厚度τ(θ):
将公式(7)带入公式(6)中进行大气衰减校准,整理后得到如下公式,经计算即可得到没有大气衰减的太阳辐射到达天线的亮温,
于公式(8)中,x、y为天线波束中心与太阳中心在方位和俯仰两个方向上的角距。
特别的,所述S5中天线方向性函数及天线波束宽度分别按照如下公式计算:
于公式(9)和公式(10)中,式中θH为天线H面的天线波束宽度,θE天线E面的天线波束宽度,G为天线增益模式,经过公式(9)和公式(10) 即可求得天线三维方向性函数和天线波束宽度。
特别的,所述S6中天线最大增益按照如下公式计算:
于公式(11)中,ΩD为太阳立体角,R为比例因子,天线模型的校准因子Cr按如下公式计算:
公式(12)中,θD为太阳角径,θA为天线波束宽度;
比例因子R按照如下公式积分得到:
特别的,所述S7中天线有效接收面积Ag和孔径效率η按照如下公式计算:
η=Ae/Ag*100 (15),
于公式(14)中,Ae为天线有效孔径。
本发明提供一种利用太阳辐射测量天线方向图的方法,具有以下有益效果:
1)同传统的微波暗室测量方法相比,本发明对环境要求低,无需其他精密测试仪表,且不需要微波暗室的特殊测试环境就可以实现微波天线方向图测试,不但具有与微波暗室相同的测量功能特点,而且减少了测试的复杂度,降低了测试成本;
2)本发明可以利用系统自身实现天线方向图的自动化测量,计算过程简单结果准确并且迅速,大幅度降低了系统测试周期和调试难度,可用于系统性能的监测和评估。
【附图说明】
图1为利用微波辐射计观测太阳示意图;
图2为利用微波辐射计扫描太阳辐射亮温结果图;
图3为利用太阳扫描数据拟合的天线方向图。
图4为在方位和俯仰过太阳中心扫描的数据拟合得到的方向图;
图5在微波暗室测量的天线方向图与利用太阳测量的天线方向图对比结果。
【具体实施方式】
于本发明中,所述的方位扫描指的是指天线仰角固定,只转动天线方位转台进行观测,俯仰扫描指的是固定天线方位转台,转动天线俯仰转台进行扫描观测。为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进一步详细说明。
以下实施例以地基多通道微波辐射计为例,利用位于西安市秦岭大气科学试验基地的MWP967KV型地基多通道微波辐射计进行了太阳扫描观测实验。
本发明提供一种利用太阳辐射测量天线方向图的方法,具体包括以下步骤:
S1根据日地轨道关系,实时计算太阳方位和高度角,具体按照如下公式进行计算:
T0=(ts-12)·15° (1),
El=arcsin(sinθlatsinδ+cosδcosT0) (2),
于公式(1)中,ts为太阳时角,于公式(2)和公式(3)中,δ为太阳赤纬,θlat为微波天线所在纬度。
S2控制方位和俯仰转台使得天线指向太阳,微波辐射计观测太阳示意图如图1所示。
S3根据计算获得天线在太阳中心周围扫描的亮温散点图,具体按照如下方法进行实施:
S31固定转台的天线仰角,通过步进方式控制转台天线指向太阳中心位置,控制方位转台在角度-10°~10°转动天线扫描太阳,观测各频点太阳辐射强度,所述太阳辐射强度为经过大气衰减到达天线的信号和大气自身辐射信号;
S32改变天线仰角,重复S31,天线方位和仰角步进范围在太阳中心 -10°~10°范围(根据实际情况设定范围),距离太阳近的时候步进角度小扫描点密集,距离太阳中心远时步进角度大一点扫描点稀疏一些。最终实现对太阳中心周围的栅格扫描;
S33扫描完成后,按如下公式计算天线观测的亮温:
于公式(4)中,所述为天线方位角,θ为俯仰角,Ωs为太阳波束立体角,ΩA为天线波束立体角,Tm为大气平均辐射温度,Tbg=2.75K为宇宙背景辐射亮温,Tsun为太阳平均辐射亮温,τ(θ)为天线指定角度的大气厚度。
S4移动天线方位使得太阳不在天线波束内,控制天线俯仰对天空进行扫描,由于在S4中扫描时间短,忽略扫描过程中的大气变化,利用俯仰天空扫描数据计算各个仰角的大气衰减,并校准大气衰减,计算没有大气衰减时观测到的太阳辐射到达天线的亮温,所述观测天空时的亮温按如下公式进行计算:
太阳辐射到达天线被接收的亮温为公式(4)和公式(5)相减,整理后得到如下公式:
按如下公式计算大气衰减,即天线指定角度的大气厚度τ(θ):
经过S3得到天线在太阳周围扫描的亮温散点图,亮温散点图经过S4 大气衰减校准后得到图2,该实施例图2给出了四个频点(但不限于)的扫描结果,由图可以看出天线正对太阳中心,观测到的数据亮温最大,天线波束偏离太阳中心时,亮温逐渐减小。
S5利用高斯模型对图2中的观测数据进行最小二乘拟合,得到天线波束宽度和天线方向性函数,如图3所示给出四个频点的三维天线方向图。由于测试天线方向图通常只测量两个方向,因此可对扫描过程进行简化,在方位和俯仰扫描时天线通过太阳中心进行正交扫描,拟合的天线方向图如图4所示,所述天线方向性函数及天线波束宽度分别按照如下公式计算:
于公式(9)和公式(10)中,式中θH为天线H面的天线波束宽度,θE天线E面的天线波束宽度,G为天线增益模式,经过公式(9)和公式(10) 即可求得天线三维方向性函数和天线波束宽度。
S6根据S5得到的天线波束宽度和天线方向性函数校准计算天线最大增益,具体按如下公式计算:
于公式(11)中,ΩD为太阳立体角,R为比例因子,天线模型的校准因子Cr按如下公式计算:
公式(12)中,θD为天线角径,θA为天线波束宽度。
比例因子R按照如下公式积分得到:
S7计算天线有效接收面积和孔径效率,具体按如下公式进行计算:
η=Ae/Ag*100 (15),
于公式(14)中,Ae为天线有效孔径。
经过上述操作,最终获得的测试结果如表1所示。
表1天线方向图测试结果
最终获得的天线方向图如图5所示,并结合传统微波暗室测量进行对比,由图5可知由两者测量结果在主瓣上完全吻合,因此本发明在达到传统测量结果精准度的同时,具有以下优点:
1)完全不依赖于微波暗室等传统的测量方法和技术,即可实现在任意时间任意外场环境下实时进行天线方向图的监测和测量,具有测量精度高、适用性更广的优点;
2)操作简单实用性强,不需要专业的测量设备和工具,降低了测试成本,提高了测试效率。
Claims (3)
1.一种利用太阳辐射测量天线方向图的方法,其特征在于,包括以下步骤:
S1实时计算太阳方位角和高度角;
S2控制转台,使天线指向太阳,所述转台为方位和俯仰转台;
S3根据计算获得天线在太阳中心周围扫描的亮温散点图,所述S3具体通过以下方法实现:
S31固定转台的天线仰角,通过步进方式控制转台天线指向太阳中心位置,控制方位转台在角度-10°~10°转动天线扫描太阳,观测各频点太阳辐射强度,所述太阳辐射强度为经过大气衰减到达天线的信号和大气自身辐射信号;
S32改变转台的天线仰角,重复S31,天线的方位和仰角步进范围在太阳中心-10°~10°范围内,实现对太阳中心周围的栅格扫描;
S33按如下公式计算天线观测的亮温:
于公式(4)中,所述为天线方位角,θ为俯仰角,Ωs为太阳波束立体角,ΩA为天线波束立体角,Tm为大气平均辐射温度,Tbg=2.75K为宇宙背景辐射亮温,Tsun为太阳平均辐射亮温,τ(θ)为天线指定角度的大气厚度;
S4移动天线方位,使得太阳不在天线波束内,控制天线俯仰对天空进行扫描,计算各个仰角的大气衰减,并校准大气衰减,计算没有大气衰减时观测到的太阳辐射到达天线的亮温,所述S4中观测天空时的亮温按如下公式计算:
太阳辐射到达天线被接收的亮温为公式(4)和公式(5)相减后得到,整理后得到如下公式:
按如下公式计算大气衰减,即天线指定角度的大气厚度τ(θ):
将公式(7)计算的大气衰减代入公式(6)中进行大气衰减校准,整理后得到如下公式,经计算即可得到没有大气衰减的太阳辐射到达天线的亮温,
于公式(8)中,x、y为天线波束中心与太阳中心在方位和俯仰两个方向上的角距;
S5利用高斯辐射模型拟合得到天线方向性函数及天线波束宽度,所述S5中天线方向性函数及天线波束宽度分别按照如下公式计算:
于公式(9)和公式(10)中,θH为天线H面的天线波束宽度,θE为天线E面的天线波束宽度,G为天线增益模型,经过公式(9)和公式(10)通过拟合即可求得天线三维方向性函数和天线波束宽度;
S6校准计算天线增益,所述S6中天线最大增益按照如下公式计算:
于公式(11)中,ΩD为太阳立体角,R为比例因子,天线模型的校准因子Cr按如下公式计算:
公式(12)中,θD为太阳角径,θA为天线波束宽度;
比例因子R按照如下公式积分得到:
S7计算天线有效接收面积和孔径效率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110635512.5A CN113341238B (zh) | 2021-06-08 | 2021-06-08 | 一种利用太阳辐射测量天线方向图的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110635512.5A CN113341238B (zh) | 2021-06-08 | 2021-06-08 | 一种利用太阳辐射测量天线方向图的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113341238A CN113341238A (zh) | 2021-09-03 |
CN113341238B true CN113341238B (zh) | 2022-12-16 |
Family
ID=77475093
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110635512.5A Active CN113341238B (zh) | 2021-06-08 | 2021-06-08 | 一种利用太阳辐射测量天线方向图的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113341238B (zh) |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105842259B (zh) * | 2016-03-18 | 2018-07-24 | 国家海洋技术中心 | 星载盐度计冷空外定标方法和装置 |
FR3071057B1 (fr) * | 2017-09-14 | 2021-02-12 | Centre Nat Rech Scient | Systeme de mesure des composantes du rayonnement solaire |
CN107942146A (zh) * | 2017-10-31 | 2018-04-20 | 西安空间无线电技术研究所 | 一种微波辐射计天线辐射效率测试方法及系统 |
CN109581078B (zh) * | 2018-11-30 | 2021-09-14 | 电子科技大学 | 一种适用于半空间环境中天线的方向图测量系统及方法 |
CN112526227B (zh) * | 2020-11-20 | 2022-09-20 | 国家卫星气象中心(国家空间天气监测预警中心) | 星载微波辐射计天线方向特性的测量方法及测量装置 |
-
2021
- 2021-06-08 CN CN202110635512.5A patent/CN113341238B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN113341238A (zh) | 2021-09-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108037374B (zh) | 一种阵列天线近场标定方法 | |
CN110618408B (zh) | 一种精密测距系统天线相位中心的系统标定方法 | |
Kummer et al. | Antenna measurements—1978 | |
CN109031359B (zh) | 一种基于微波暗室的gnss接收天线绝对相位中心校正方法 | |
CN109283511B (zh) | 一种宽覆盖多波束接收基阵校准方法 | |
CN110308746A (zh) | 一种适用于三轴转台式测控天线的星体标校方法 | |
CN109581279A (zh) | 一种超宽带多波束系统测向精度的校正方法及装置 | |
CN103616569A (zh) | 一种毫米波平面近场测试相位修正方法 | |
CN105606906B (zh) | 一种毫米波相控阵测试标定方法 | |
CN109712201B (zh) | 一种广域相机的定位能力标定装置及标定方法 | |
CN110887568B (zh) | 一种月球观测系统 | |
CN111366793B (zh) | 一种减小截断误差的平面天线近场测量方法及系统 | |
CN110058204A (zh) | 一种基于方向图匹配的星载天线波束中心定标方法 | |
CN107167807B (zh) | 三维天线阵列综合孔径辐射计分段图像反演方法 | |
CN104375129B (zh) | 一种分布式阵列相参合成雷达发射相参参数校准方法 | |
WO2019071917A1 (zh) | 卫星跟踪方法 | |
CN110907904A (zh) | 一种太阳法方位波束宽度修正方法 | |
CN111929651A (zh) | 一种综合孔径辐射计系统相位残差标定方法 | |
CN116520035A (zh) | 一种采用平面近场快速反演方法的二维方向图测试方法 | |
CN113341238B (zh) | 一种利用太阳辐射测量天线方向图的方法 | |
CN115865142B (zh) | 一种雷达天线波束指向校准方法、系统、装置及存储介质 | |
Fernández et al. | Advances in antenna measurement and characterization using unmanned aerial vehicles | |
CN105929361A (zh) | 单天线最优化比幅无线电测向系统及方法 | |
CN113702402B (zh) | 一种水汽辐射计方位仰角自动校准方法 | |
CN114265025A (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 |