CN105510365A - 地面反射系数的现场测量方法 - Google Patents

地面反射系数的现场测量方法 Download PDF

Info

Publication number
CN105510365A
CN105510365A CN201610030989.XA CN201610030989A CN105510365A CN 105510365 A CN105510365 A CN 105510365A CN 201610030989 A CN201610030989 A CN 201610030989A CN 105510365 A CN105510365 A CN 105510365A
Authority
CN
China
Prior art keywords
antenna
ground
ground reflection
value
estimated value
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
CN201610030989.XA
Other languages
English (en)
Other versions
CN105510365B (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.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
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 Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN201610030989.XA priority Critical patent/CN105510365B/zh
Publication of CN105510365A publication Critical patent/CN105510365A/zh
Application granted granted Critical
Publication of CN105510365B publication Critical patent/CN105510365B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
    • G01N23/20Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by using diffraction of the radiation by the materials, e.g. for investigating crystal structure; by using scattering of the radiation by the materials, e.g. for investigating non-crystalline materials; by using reflection of the radiation by the materials
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00Investigating materials by wave or particle radiation
    • G01N2223/05Investigating materials by wave or particle radiation by diffraction, scatter or reflection
    • G01N2223/052Investigating materials by wave or particle radiation by diffraction, scatter or reflection reflection
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00Investigating materials by wave or particle radiation
    • G01N2223/60Specific applications or type of materials
    • G01N2223/614Specific applications or type of materials road surface
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00Investigating materials by wave or particle radiation
    • G01N2223/60Specific applications or type of materials
    • G01N2223/616Specific applications or type of materials earth materials

Landscapes

  • Chemical & Material Sciences (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明涉及一种地面反射系数的现场测量方法,属于雷达信号处理领域。由于地面反射的影响,雷达天线俯仰向方向图会发生严重的波束分裂,进而影响到目标的探测和仰角的测量,为了更精确的得到不同地表情况下的地面反射系数,提出地面反射系数的现场测量方法,该方法包括如下步骤:(1)通道校正,并得到归一化所需各个参量;(2)步进测量,记录每次的上下天线单元比幅值;(3)计算地面反射系数的估计值;(4)根据实际情况设置扩散因子D和地表面高度标准偏差σh,计算地面复介电常数εr的估计值;(5)选取误差最小的εr。该方法可以更好的计算不同地表的地面反射系数,更准确的绘制天线俯仰方向图。

Description

地面反射系数的现场测量方法
技术领域
本发明涉及一种地面反射系数的现场测量方法,用于天线俯仰方向图的估算,属于雷达信号处理领域。
背景技术
雷达天线所接收到的回波信号除了目标反射的直达信号还有回波经地面反射后到达接收天线的信号,如附图1所示,由于地面反射的影响,俯仰向方向图会发生严重的波束分裂,进而影响到目标的探测和仰角的测量,不同仰角目标探测威力和自由空间相比会有较大的偏差,仰角的测量也会受影响,若能对其进行有效的测量,不仅可以更好的估计系统性能,而且可以据此分析误差产生的来源并采取避免措施,从而改善系统精度。
现有的地面反射系数的测量手段主要有:1)数学建模计算法。不同地面的反射系数可根据模型得到,地面反射系数是雷达参数(入射角、频率、极化)和地面参数(粗糙度、介电常数)的复杂函数,对于各个参数都只是粗略估计,由模型得到的系数存在比较大的误差;
2)查表法。由于地面反射的复杂性、组成成分的多样和湿度变化都会影响反射系数的变化,所以查表法得到的地面反射系数也只是差强人意。
发明内容
本发明的目的是为了解决现有技术存在系数误差较大的问题,提供一种地面反射系数的现场测量方法。
本发明的目的是通过如下技术方案实现的。
地面反射系数的现场测量方法,具体步骤如下:
步骤一、将测试天线与接收天线放置在待查场地中间位置,将测试天线连接信号源,并设置信号源频率,需保证该频率避开周围环境电磁干扰;
所述接收天线为带有上下两个天线单元的天线;
步骤二、记录步骤一测试天线信号源功率,记录两个天线单元的信号幅度值,记录两天线之间的距离;
当步骤一所述的接收天线的两个天线单元增益误差大于等于0.5dB时,对接收天线的两个天线单元进行通道校正;
步骤三、接收天线接收测试天线输出的点频信号,然后计算得到两个天线单元的信号比幅值;
在计算比幅值前,将步骤三所述的点频信号进行快速傅里叶变换(FFT)以提高信噪比;并通过多次测量检验接收的数据是否稳定;
步骤四、根据步骤三测得的比幅值,计算地面反射系数的估值:
对雷达上天线单元,测试单元波束仰角为α1,反射波束仰角为θ1;对雷达下天线单元,测试单元波束仰角为α2,反射波束仰角为θ2,上下两天线单元的比幅值为反射路径和直达路径间传播相位差为
其中dbase为测试天线距离接收天线的距离,h1为天线下单元相位中心距地面高度,h2为天线上单元相位中心距地面高度,h3为测试天线相位中心距地面高度,λ为波长,发射天线方向图为FT(θ)=sinc(d/λ*sinθ)(d为天线口径),波束仰角为θd
步骤五、通过测量计算出地表面高度的标准偏差σh;设置扩散因子D小于等于1;然后结合步骤四得到的地面反射系数的估值计算地面水平极化的复介电常数的估计值或者垂直极化的复介电常数的估计值;
水平极化的复介电常数的估计值为:
ϵ r = { s i n θ [ D exp ( - ( 4 πσ h s i n θ λ ) 2 / 2 ) - ρ ^ ] ρ ^ + D exp ( - ( 4 πσ h s i n θ λ ) 2 / 2 ) } 2 + cos 2 θ - - - ( 2 )
垂直极化的复介电常数的估计值为:
ϵ r = - b ± b 2 - 4 a c 2 a - - - ( 3 )
其中:
a = sin 2 θ · [ ρ ^ - D exp ( - ( 4 πσ h s i n θ λ ) 2 / 2 ) ] 2 - - - ( 4 )
b = - [ ρ ^ + D exp ( - ( 4 πσ h s i n θ λ ) 2 / 2 ) ] 2 - - - ( 5 )
c = [ ρ ^ + D exp ( - ( 4 πσ h s i n θ λ ) 2 / 2 ) ] 2 - - - ( 6 )
式中θ为测试单元相对于接收天线上天线单元或者下天线单元的波束仰角。
步骤六、步骤五所得的εr作为初值进行搜索,取使得估计值与步骤二测得的比幅值误差最小的εr作为复介电常数的值,并以此值计算地面反射系数:
ρ ( θ ) = ( ϵ r ) m sin θ - ϵ r - cos 2 θ ( ϵ r ) m sin θ + ϵ r - cos 2 θ D exp ( - ( 4 πσ h sin θ λ ) 2 / 2 ) - - - ( 7 )
其中,水平极化时m=0,垂直极化时m=1。
对步骤六所得的地面反射系数进行验证,方法为:
1、移动测试天线,进行多点测量,能够得到多个两个天线单元的信号幅度值,以及两天线之间的距离;
将多个幅度值进行距离归一化、功率归一化、幅度归一化,并乘以系统误差系数,得到天线俯仰向方向图的实测值。
2、通过步骤六所得的地面反射系数的值计算天线俯仰方向图的估计值:
式中为天线上单元或者下单元地面反射波束仰角。
3、当天线俯仰向方向图的实测值与估计值相吻合时,说明所测得的地面反射系数为有效值。
有益效果
本发明提出的地面反射系数的现场测量方法,配置简单,操作容易,可以方便计算出相应环境的地面反射系数,且相对于传统的模型法和查表法,误差更小,可以更准确的估计天线俯仰方向图。
附图说明
图1是地面反射影响示意图;
图2是地面反射系数的现场测量配置示意图;
图3是上天线单元俯仰方向图;
图4是下天线单元俯仰方向图。
具体实施方式
下面将结合附图及具体实施方式对本发明加以详细说明,需要指出的是,所描述的实施例仅旨在便于对本发明的理解,而不起任何限定作用。
实施例1
地面反射系数的现场测量方法,具体步骤如下:
步骤一、将测试天线与接收天线都为水平极化,且放置在柏油路面中间位置,具体配置图如图2所示,接收天线下单元相位中心距地面高度为1.95m,接收天线上单元相位中心距地面高度为2.47m,测试天线单元相位中心距地面高度为2.21m,将测试天线连接信号源,并设置信号源频率为656MHz,发射功率0dB,在测试天线距离接收天线2m处进行接收通道校正;
步骤二、记录步骤一测试天线信号源功率0dB,记录两个天线单元的信号幅度值,记录两天线之间的距离2m;
步骤三、将测试天线移动到10m处,接收天线接收测试天线输出的点频信号,然后计算得到两个天线单元的信号比幅值;
在计算比幅值前,将步骤三所述的点频信号进行快速傅里叶变换(FFT)以提高信噪比;
测量3次,看数据是否稳定,若稳定,记录比幅值作为测试数据。
步骤四、根据步骤三测得的比幅值,计算地面反射系数的估值:
步骤五、通过测量计算出地表面高度的标准偏差σh=0.01;设置扩散因子D=1;然后结合步骤四得到的地面反射系数的估值计算地面水平极化的复介电常数的估计值或者垂直极化的复介电常数的估计值;
水平极化的复介电常数的估计值为:
ϵ r = { s i n θ [ D exp ( - ( 4 πσ h s i n θ λ ) 2 / 2 ) - ρ ^ ] ρ ^ + D exp ( - ( 4 πσ h s i n θ λ ) 2 / 2 ) } 2 + cos 2 θ - - - ( 10 )
垂直极化的复介电常数的估计值为:
ϵ r = - b ± b 2 - 4 a c 2 a - - - ( 11 )
其中:
a = sin 2 θ · [ ρ ^ - D exp ( - ( 4 πσ h s i n θ λ ) 2 / 2 ) ] 2 - - - ( 12 )
b = - [ ρ ^ + D exp ( - ( 4 πσ h s i n θ λ ) 2 / 2 ) ] 2 - - - ( 13 )
c = [ ρ ^ + D exp ( - ( 4 πσ h s i n θ λ ) 2 / 2 ) ] 2 - - - ( 14 )
步骤六、采用步骤五中的介电常数,并将实部和虚部各取其本身0.5倍~1.5倍作为变化范围,步进选取能使实测值和理论值总误差最小的介电常数。最后得到εr=2.4451+1.6111i,并据此计算地面反射系数:
ρ ( θ ) = ( ϵ r ) d s i n θ - ϵ r - cos 2 θ ( ϵ r ) d s i n θ + ϵ r - cos 2 θ D exp ( - ( 4 πσ h s i n θ λ ) 2 / 2 ) - - - ( 15 )
对步骤六所得的地面反射系数进行验证,方法为:
1、移动测试天线,进行多点测量,测试天线按11m~20m之间每1m步进一次,20m~30m之间每2m步进一次,30m~60m每5m步进一次(也可按照所需要测量的角度范围反推距离步进的范围),共进行21次步进测量,能够得到多个两个天线单元的信号幅度值,以及两天线之间的距离;
将多个幅度值进行距离归一化、功率归一化、幅度归一化,并乘以系统误差系数,得到天线俯仰向方向图的实测值。
其中,上天线单元的校正幅度值为6.8376×108-1.0256×108i,下天线单元的校正幅度值为6.3109×108-9.9383×107i。因系统误差影响,归一化后的天线俯仰方向图需要乘以系统误差系数(本次实验为0.49)。
2、通过步骤六所得的地面反射系数的值计算天线俯仰方向图的估计值:
3、当天线俯仰向方向图的实测值与估计值相吻合时,说明所测得的地面反射系数为有效值。
附图3和附图4分别为根据本文所提的地面反射系数的现场测量方法进行实验得出的天线俯仰方向图。方向一致所以实验测得的地面反射系数是准确的。
以上所述的具体描述,对发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施实例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (5)

1.地面反射系数的现场测量方法,其特征在于:具体步骤如下:
步骤一、将测试天线与接收天线放置在待查场地中间位置,将测试天线连接信号源,并设置信号源频率,需保证该频率避开周围环境电磁干扰;
步骤二、记录步骤一测试天线信号源功率,记录两个天线单元的信号幅度值,记录两天线之间的距离;
步骤三、接收天线接收测试天线输出的点频信号,然后计算得到两个天线单元的信号比幅值;
步骤四、根据步骤三测得的比幅值,计算地面反射系数的估值:
对雷达上天线单元,测试单元波束仰角为α1,反射波束仰角为θ1;对雷达下天线单元,测试单元波束仰角为α2,反射波束仰角为θ2,上下两天线单元的比幅值为反射路径和直达路径间传播相位差为其中dbase为测试天线距离接收天线的距离,h1为天线下单元相位中心距地面高度,h2为天线上单元相位中心距地面高度,h3为测试天线相位中心距地面高度,λ为波长,发射天线方向图为FT(θ)=sinc(d/λ*sinθ)(d为天线口径),波束仰角为θd
步骤五、通过测量计算出地表面高度的标准偏差σh;设置扩散因子D小于等于1;然后结合步骤四得到的地面反射系数的估值计算地面水平极化的复介电常数的估计值或者垂直极化的复介电常数的估计值;
水平极化的复介电常数的估计值为:
垂直极化的复介电常数的估计值为:
其中:
式中θ为测试单元相对于接收天线上天线单元或者下天线单元的波束仰角;
步骤六、步骤五所得的εr作为初值进行搜索,取使得估计值与步骤二测得的比幅值误差最小的εr作为复介电常数的值,并以此值计算地面反射系数:
其中,水平极化时m=0,垂直极化时m=1。
2.如权利要求1所述的地面反射系数的现场测量方法,其特征在于:对所述步骤六所得的地面反射系数进行验证,方法为:
(1)、移动测试天线,进行多点测量,能够得到多个两个天线单元的信号幅度值,以及两天线之间的距离;
将多个幅度值进行距离归一化、功率归一化、幅度归一化,并乘以系统误差系数,得到天线俯仰向方向图的实测值;
(2)、通过步骤六所得的地面反射系数的值计算天线俯仰方向图的估计值:
式中为天线上单元或者下单元地面反射波束仰角;
(3)、当天线俯仰向方向图的实测值与估计值相吻合时,说明所测得的地面反射系数为有效值。
3.如权利要求1或2所述的地面反射系数的现场测量方法,其特征在于:所述接收天线为带有上下两个天线单元的天线。
4.如权利要求1所述的地面反射系数的现场测量方法,其特征在于:当步骤一所述的接收天线的两个天线单元增益误差大于等于0.5dB时,需对接收天线的两个天线单元进行通道校正。
5.如权利要求1所述的地面反射系数的现场测量方法,其特征在于:在计算比幅值前,将步骤三所述的点频信号进行快速傅里叶变换(FFT)以提高信噪比;并通过多次测量检验接收的数据是否稳定。
CN201610030989.XA 2016-01-18 2016-01-18 地面反射系数的现场测量方法 Expired - Fee Related CN105510365B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610030989.XA CN105510365B (zh) 2016-01-18 2016-01-18 地面反射系数的现场测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610030989.XA CN105510365B (zh) 2016-01-18 2016-01-18 地面反射系数的现场测量方法

Publications (2)

Publication Number Publication Date
CN105510365A true CN105510365A (zh) 2016-04-20
CN105510365B CN105510365B (zh) 2018-05-15

Family

ID=55718497

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610030989.XA Expired - Fee Related CN105510365B (zh) 2016-01-18 2016-01-18 地面反射系数的现场测量方法

Country Status (1)

Country Link
CN (1) CN105510365B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111399006A (zh) * 2020-04-07 2020-07-10 北京理工大学 一种高灵敏度gnss载波跟踪环路优化方法
CN113850999A (zh) * 2021-09-03 2021-12-28 杭州海康威视数字技术股份有限公司 车位检测装置和用于监控车位的具有雷达的摄像机装置

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102647693B1 (ko) * 2016-11-28 2024-03-15 주식회사 에이치엘클레무브 레이더 장치 및 그의 오차 보정 방법

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU1601590A1 (ru) * 1988-06-20 1990-10-23 Предприятие П/Я М-5566 Способ измерени коэффициента отражени листовых материалов
CN103308909A (zh) * 2013-03-06 2013-09-18 西安电子科技大学 基于多参数交替搜索的米波雷达低仰角测高方法
CN103885054A (zh) * 2014-02-28 2014-06-25 西安电子科技大学 一种基于分布源反射模型的米波雷达低仰角测高方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU1601590A1 (ru) * 1988-06-20 1990-10-23 Предприятие П/Я М-5566 Способ измерени коэффициента отражени листовых материалов
CN103308909A (zh) * 2013-03-06 2013-09-18 西安电子科技大学 基于多参数交替搜索的米波雷达低仰角测高方法
CN103885054A (zh) * 2014-02-28 2014-06-25 西安电子科技大学 一种基于分布源反射模型的米波雷达低仰角测高方法

Non-Patent Citations (9)

* Cited by examiner, † Cited by third party
Title
李宗谦等: "从雷达后向散射系数反演土壤湿度与复介电常数", 《中国科学(E辑)》 *
李相平等: "地杂波后向散射系数经验模型研究", 《舰船电子工程》 *
罗续成等: "平直地面对一类圆极化波的等效反射系数", 《宇航学报》 *
郦能敬等: "《对空情报雷达总体论证理论与实践》", 31 December 2008, 国防工业出版社 *
陈伯孝等: "《现代雷达系统分析与设计》", 30 September 2012, 国防工业出版社 *
马保科等: "地面电磁(光)波散射特性研究", 《西安工业大学学报》 *
黄南晖: "地质雷达探测的波场分析", 《地球科学—中国地质大学学报》 *
黄忠来等: "利用探地雷达频谱反演层状介质几何与电性参数", 《地球物理学报》 *
齐国雷等: "FDTD 方法分析高功率微波粗糙地面散射特性", 《强激光与粒子束》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111399006A (zh) * 2020-04-07 2020-07-10 北京理工大学 一种高灵敏度gnss载波跟踪环路优化方法
CN111399006B (zh) * 2020-04-07 2020-11-24 北京理工大学 一种高灵敏度gnss载波跟踪环路优化方法
CN113850999A (zh) * 2021-09-03 2021-12-28 杭州海康威视数字技术股份有限公司 车位检测装置和用于监控车位的具有雷达的摄像机装置

Also Published As

Publication number Publication date
CN105510365B (zh) 2018-05-15

Similar Documents

Publication Publication Date Title
CN103064073B (zh) 一种基于频率捷变改变雷达目标特性的方法
CN103674898B (zh) 一种基于激光测距的反射率测试自动校正方法
CN101900692B (zh) 大面积土壤湿度测量方法
CN103064098B (zh) 卫星导航数字多波束发射阵列天线指向偏差修正方法
CN109669075B (zh) 基于开口矩形波导的介质复介电常数无损反射测量方法
TWI540792B (zh) 天線系統的遠場調校系統
CN102854401A (zh) 一种阵列天线方向图时域测量方法
CN103728321B (zh) 多功能材料电磁参数测试系统及测试方法
CN106291502B (zh) 目标rcs测量中背景提取与抵消的最大概率时域处理方法
CN103412227A (zh) 一种基于绕射抑制的频选雷达天线罩传输性能测试系统及其方法
CN108037374A (zh) 一种阵列天线近场标定方法
CN105136073B (zh) 一种在边坡形变监测中的气象校正模型
CN102664690B (zh) 一种智能天线doa估计性能的有线测试方法
CN105510365A (zh) 地面反射系数的现场测量方法
CN102967853B (zh) 雷达散射截面测试场背景信号获取方法
CN103217211A (zh) 一种基于合成孔径原理的变电站噪声源分布测量方法
CN106199220B (zh) 基于光程差校正的阵列天线相位一致性测量方法
CN102798778A (zh) 内场天线测量系统信号传递环节建模方法
CN103257340A (zh) 一种利用雷达卫星标定多台地面接收机幅度一致性的方法
CN104237876B (zh) 目标支架到定标支架的距离的确定方法及装置
CN112859028A (zh) 一种外场地物散射时域回波采集及谱分析方法及系统
CN111487474B (zh) 数值孪生电磁测量系统
CN113917241B (zh) 一种天线方向图快速测量和预估方法、系统、设备及终端
CN115524540A (zh) 场地散射不确定度仿真模型构建方法、装置、设备及介质
Zheng et al. A simple method for on-wafer antenna gain measurement

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

Granted publication date: 20180515

Termination date: 20200118