CN104765029B - 一种基于射线追踪确定表面波导盲区信息的方法 - Google Patents

一种基于射线追踪确定表面波导盲区信息的方法 Download PDF

Info

Publication number
CN104765029B
CN104765029B CN201510178625.1A CN201510178625A CN104765029B CN 104765029 B CN104765029 B CN 104765029B CN 201510178625 A CN201510178625 A CN 201510178625A CN 104765029 B CN104765029 B CN 104765029B
Authority
CN
China
Prior art keywords
mrow
ray
msub
height
blind area
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
CN201510178625.1A
Other languages
English (en)
Other versions
CN104765029A (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.)
Xidian University
Original Assignee
Xidian University
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 Xidian University filed Critical Xidian University
Priority to CN201510178625.1A priority Critical patent/CN104765029B/zh
Publication of CN104765029A publication Critical patent/CN104765029A/zh
Application granted granted Critical
Publication of CN104765029B publication Critical patent/CN104765029B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/40Means for monitoring or calibrating

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

一种基于射线追踪确定表面波导盲区信息的方法,步骤如下:获取修正折射率剖面数据;确定计算范围,包括最大传播距离,最大传播高度,高度步进值,发射天线高度以及发射天线初始发射仰角及射线仰角间隔,判断是否对射线进行追踪;根据Snell定律的二阶泰勒近似模式对射线进行追踪,获得每一根射线的传播轨迹;对于所有符合追踪条件的射线进行陷获与否的判定并分类,提取出射射线的离散数据;获得波导顶部盲区的起始位置参数;计算表面波导顶部盲区面积。本发明引入射线追踪法对电波射线进行判定、归类,提取及数据离散化处理,简单直观并且准确地获得盲区的起始位置、范围,使盲区计算模型的复杂程度降低,从而实现对流层表面波导盲区面积的快速计算。

Description

一种基于射线追踪确定表面波导盲区信息的方法
技术领域
本发明涉及雷达探测及电波传播计算,尤其涉及一种确定表面波导盲区面积的计算方法。
背景技术
雷达的实际探测受制于其工作体制、信号检测能力、天线类型等因素,会存在一定的电磁盲区。与此同时,不同的探测环境也会对其电磁盲区造成影响。这一点在低空对流层海上大气环境中体现的尤为明显。由于海面上空,对流层底部往往存在折射率梯度极端异常的现象,会形成大气波导,使得一部分电波陷获在波导层内,损耗变小,传播到原本无法探测到的区域,形成超视距传播。但是另一方面而言,又会改变电波的正常传播轨迹,形成电波无法到达的电磁盲区。如今,对于如何利用大气波导的普遍存在性来实现海上微波雷达超视距探测已经成为一个研究热点。但对于伴随形成的电磁盲区的细致研究还较少有成果涉及。电磁盲区对于实际海上军事意义是极为重大的——近海面低空突防与反突防问题,随着我国海洋战略的逐步推进,必须获得足够的重视。
目前有关海上对流层大气波导的电磁盲区分析主要是针对蒸发波导而言的,对于表面波导的分析还相对停留在存在性的说明上,缺乏有效的定量的计算与分析。因而本发明利用射线追踪方法(Ray Trace)提出一种针对海上表面波导盲区信息的分析方法。
发明内容
本发明的目的是提供一种可以基于射线追踪的确定对流层表面波导盲区的计算方法。
为了实现上述目的,本发明采取如下的技术解决方案:
一种基于射线追踪确定表面波导盲区信息的方法,包括以下步骤:
步骤1、获取修正折射率剖面数据;
步骤2、确定计算范围,包括最大传播距离Xmax,最大传播高度Hmax,高度步进值h,发射天线高度Ht以及发射天线初始发射仰角及射线仰角间隔,判断是否对射线进行追踪;
步骤3、根据Snell定律的二阶泰勒近似模式对射线进行追踪,获得每一根射线的传播轨迹;
步骤4、对于所有符合追踪条件的射线进行陷获与否的判定并分类,提取出射射线的离散数据;
步骤4a、陷获状态的判断步骤如下:判断射线是否能够传播到达最大传播高度Hmax,如果是则将该射线归为出射射线,否则归为陷获射线,统计出射射线的数目;
步骤4b、对于每一根出射射线,提取表征其位置的三维离散数据[rn,xni,yni],其中,rn表示第n条出射射线,yni表示第n条出射射线上第i个步进点的高度值,xni表示第n条出射射线上对应于yni的水平距离;
步骤5、获得波导顶部盲区的起始位置参数;
步骤5a、找出xni的最大值所在的射线rxmax,获取表征该射线轨迹位置的二维离散数据[xxmaxi,yxmaxi];
步骤5b、获取射线rxmax上最接近波导层顶高的步进点P1的坐标(XP1,YP1),以及射线rxmax上最接近最大传播高度的步进点P2的坐标(XP2,YP2),射线rxmax在P1点到P2点之间的曲线轨迹表征了表面波导顶部盲区的开始位置信息;
步骤6、计算表面波导顶部盲区面积S;
式中的xi、xi+1、yi、yi+1分别表示相邻步进点的坐标值。
更具体的,所述步骤2中判断是否对射线进行追踪的条件为:若发射天线初始发射仰角绝对值大于3°,则不对其进行追踪,否则对其进行追踪。
更具体的,所述步骤2中判断是否对射线进行追踪的条件为:若发射天线高度大于波导顶高,则不对其进行追踪,否则对其进行追踪。
本发明从电波传播机理上,确定表面波导顶部盲区的相关信息,引入射线追踪法对电波射线判定、归类,提取离散化数据,可以简单、直观并且准确地获得盲区的起始位置、范围,并给出盲区面积的快速计算方法,在军用海上低空突防与预警反突防,民用通信等领域具有显著的作用。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中需要使用的附图做简单介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明方法的流程图;
图2为采用本发明的方法计算的对流层表面波导的盲区面积的结果图。
具体实施方式
为了让本发明的上述和其它目的、特征及优点能更明显,下文特举本发明实施例,并配合所附图示,做详细说明如下。
在下面的描述中阐述了很多具体细节以便于充分理解本发明,但是本发明还可以采用其它不同于在此描述的其它方式来实施,本领域技术人员可以在不违背本发明内涵的情况下做类似推广,因此本发明不受下面公开的具体实施例的限制。
结合图1所示,本发明方法包括以下步骤:发射天线发出射线,
步骤1、获取修正折射率剖面数据;
根据表面波导三折线模型获取表面波导四个基本输入参数:基础层高度HB、波导层厚度HD、基础层大气修正折射率随高度变化斜率SE1以及波导层大气修正折射率随高度变化斜率SE2,根据以上四参数获得修正折射率剖面数据;
步骤2、确定计算范围:最大传播距离Xmax(单位:千米),最大传播高度Hmax(单位:米),高度步进值h(单位:米),发射天线高度Ht(单位:米)以及发射天线初始发射仰角及射线仰角间隔,判断是否对射线进行追踪;
追踪的条件为:若发射天线初始发射仰角绝对值大于3°,或发射天线高度大于波导顶高,即Ht>(HB+HD),则不对其进行追踪,若发射天线初始发射仰角绝对值小于等于3°,或发射天线高度小于等于波导顶高,则进行射线追踪;
步骤3、根据Snell定律的二阶泰勒近似模式对射线进行追踪,获得每一根射线的传播轨迹;
将大气视为未知环境,通过射线传播的前一步对后一步进行判断,射线上相邻两点的高度分别为h1、h2,该两点的仰角分别为θ1、θ2,大气修正折射指数m1、m2通过步骤1中修正折射率剖面数据获得,m在h1与h2高度间随距离线性变化,斜率为k,根据如下方程组通过迭代求解获得每一根射线的传播轨迹;
本步骤采用现有技术中的方法来获取射线的传播轨迹,此处不做赘叙;
步骤4、对于所有符合追踪条件的射线进行陷获与否的判定并分类,提取出射射线的离散数据;
步骤4a、陷获状态的判断步骤如下:判断射线是否能够传播到达最大传播高度Hmax,如果是则将该射线归类为出射射线,否则归类为陷获射线,统计出射射线的数目;
步骤4b、对于每一根出射射线,提取表征其位置的三维离散数据[rn,xni,yni],其中,rn表示第n条出射射线,yni表示传播高度范围内第n条出射射线上第i个步进点处的高度值,xni表示第n条出射射线上对应于yni的水平距离;
出射射线的三维离散数据从步骤3获得的射线的传播轨迹获取;
本步骤中,对发射天线发出的射线按照其陷获状态进行分类,并提取每一条射线的位置数据,从而方便研究电波在陷获或非陷获状态时不同的传播特点,为研究波导顶部电磁盲区信息打下基础;
步骤5、获得波导顶部盲区的起始位置参数;
步骤5a、结合图2,图2中x轴为传播距离,y轴为传播高度,对于最大传播高度Hmax处,找出xni的最大值所在的射线rxmax,即找出所有向上出射射线中位于最外部的一根,取出表征该射线轨迹位置的二维离散数据[xxmaxi,yxmaxi];
步骤5b、获取射线rxmax上最接近波导层顶高(HB+HD)的步进点P1的坐标(XP1,YP1),以及射线rxmax上最接近最大传播高度Hmax的步进点P2的坐标(XP2,YP2),射线rxmax在P1点到P2点之间的曲线轨迹表征了表面波导顶部盲区的开始位置信息;
本步骤中,通过对出射射线的位置关系比较,确定出当前大气参数、相关雷达参数计算条件下,电波传播所能到达的最远处轨迹曲线,并给出其离散点的水平距离及高度数据,作为盲区起始位置,更为直观地确定了表面波导顶部盲区的位置信息;
步骤6、计算表面波导顶部盲区面积S;
如图2所示,盲区区域可以看作是由以下几条曲线:y=YP2、y=YP1和x=Xmax围成,将整个盲区以直线x=XP2垂直划分为两部分:射线rxmax在P1点到P2点之间的曲线以下、YP1高度以上的相邻步进点之间所围成的梯形区域和长为(Xmax-XP2)、宽为[Hmax-(HB+HD)]的矩形区域,
当XP1<x<XP2时,运用微分求和思想对射线rxmax在P1点到P2点曲线以下,YP1高度以上的相邻步进点之间所围成的区域的面积S1累加求和;在高度步进值h取值较小的情况下可以近似表示曲线以下盲区部分面积;
当XP2<x<Xmax时,根据矩形的长(Xmax-XP2)和宽[Hmax-(HB+HD)]求出矩形区域的面积S2
则盲区面积S为:
式中的xi、xi+1、yi、yi+1分别表示相邻步进点的坐标值,yi从YP1开始取起,yi+1从YP1以上一个步进点取起,yi+1的最后一个取值步进点为YP2,xi+1,xi的取值与yi+1,yi相对应。
本发明在获取表面波导顶部盲区大小信息时,通过确定盲区的边界位置,根据射线的位置数据微分求和获得面积大小,相较于现有使用传播损耗计算法获得盲区大小,更加直观,边界更加清晰。
本发明提出的基于射线追踪确定表面波导的盲区计算方法,从电波传播机理上,确定了表面波导顶部盲区的相关信息。引入射线追踪法对电波射线进行判定、归类,提取及数据离散化处理,可以简单、直观并且准确地获得盲区的起始位置、范围,可以使得计算结果更具清晰的物理意义,反映出盲区的位置分布信息。并且具有较之于传统的抛物方程方法更为快速方便、运算量小等优势,使得盲区的计算模型的复杂程度降低,从而可以实现针对对流层表面波导盲区面积的快速计算。该问题的研究在军用海上低空突防与预警反突防,以及民用通信等领域具有显著的作用。
下面结合图2,以一个算例实验对本发明方法的效果进行说明:
1)算例计算条件
设定大气表面波导四参数分别为:
HB=50m,HD=50m,SE1=100M/km,SE2=-300M/km;
确定计算参数的范围为:
Xmax=300km,Hmax=500m,h=1m,Ht=60m;
设定初始仰角范围[-1°,1°],射线仰角间隔0.005°;
2)计算结果分析
经过射线追踪计算,取出的最外部射线段为
P1(40.30km,100m)
P2(122.44km,500m);
P1点与P2点之间射线段为顶部盲区起始位置。
经计算,阴影部分面积S1=10.46km2,矩形区域面积S2=71.00km2,则盲区总面积S=81.46km2
以上所述,仅是本发明的较佳实施例而已,并非对本发明做任何形式上的限制,虽然本发明已以较佳实施例揭露如上,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,当可利用上述揭示的技术内容做出些许更动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。

Claims (1)

1.一种基于射线追踪确定表面波导盲区信息的方法,其特征在于,包括以下步骤:
步骤1、获取修正折射率剖面数据;
步骤2、确定计算范围,包括最大传播距离Xmax,最大传播高度Hmax,高度步进值h,发射天线高度Ht以及发射天线初始发射仰角及射线仰角间隔,判断是否对射线进行追踪,追踪的条件为:若发射天线初始发射仰角绝对值大于3°,或发射天线高度大于波导顶高,则不对其进行追踪,若发射天线初始发射仰角绝对值小于等于3°,或发射天线高度小于等于波导顶高,则进行射线追踪;
步骤3、根据Snell定律的二阶泰勒近似模式对射线进行追踪,获得每一根射线的传播轨迹;
步骤4、对于所有符合追踪条件的射线进行陷获与否的判定并分类,提取出射射线的离散数据;
步骤4a、陷获状态的判断步骤如下:判断射线是否能够传播到达最大传播高度Hmax,如果是则将该射线归为出射射线,否则归为陷获射线,统计出射射线的数目;
步骤4b、对于每一根出射射线,提取表征其位置的三维离散数据[rn,xni,yni],其中,rn表示第n条出射射线,yni表示第n条出射射线上第i个步进点的高度值,xni表示第n条出射射线上对应于yni的水平距离;
步骤5、获得波导顶部盲区的起始位置参数;
步骤5a、找出xni的最大值所在的射线rxmax,获取表征该射线轨迹位置的二维离散数据[xxmaxi,yxmaxi];
步骤5b、获取射线rxmax上最接近波导顶高的步进点P1的坐标(XP1,YP1),以及射线rxmax上最接近最大传播高度的步进点P2的坐标(XP2,YP2),射线rxmax在P1点到P2点之间的曲线轨迹表征了表面波导顶部盲区的开始位置信息;
步骤6、计算表面波导顶部盲区面积S;
<mrow> <mi>S</mi> <mo>=</mo> <msub> <mi>S</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>S</mi> <mn>2</mn> </msub> <mo>=</mo> <mo>&amp;lsqb;</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <munderover> <mo>&amp;Sigma;</mo> <msub> <mi>Y</mi> <mrow> <mi>P</mi> <mn>1</mn> </mrow> </msub> <msub> <mi>Y</mi> <mrow> <mi>P</mi> <mn>2</mn> </mrow> </msub> </munderover> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>+</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <mo>-</mo> <msub> <mi>Y</mi> <mrow> <mi>P</mi> <mn>1</mn> </mrow> </msub> <mo>&amp;CenterDot;</mo> <mrow> <mo>(</mo> <msub> <mi>X</mi> <mrow> <mi>P</mi> <mn>2</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>X</mi> <mrow> <mi>P</mi> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mrow> <mo>(</mo> <msub> <mi>X</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>X</mi> <mrow> <mi>P</mi> <mn>2</mn> </mrow> </msub> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <mrow> <mo>(</mo> <msub> <mi>Y</mi> <mrow> <mi>P</mi> <mn>2</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>Y</mi> <mrow> <mi>P</mi> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
式中的xi、xi+1、yi、yi+1分别表示相邻步进点的坐标值。
CN201510178625.1A 2015-04-15 2015-04-15 一种基于射线追踪确定表面波导盲区信息的方法 Expired - Fee Related CN104765029B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510178625.1A CN104765029B (zh) 2015-04-15 2015-04-15 一种基于射线追踪确定表面波导盲区信息的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510178625.1A CN104765029B (zh) 2015-04-15 2015-04-15 一种基于射线追踪确定表面波导盲区信息的方法

Publications (2)

Publication Number Publication Date
CN104765029A CN104765029A (zh) 2015-07-08
CN104765029B true CN104765029B (zh) 2017-09-05

Family

ID=53646979

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510178625.1A Expired - Fee Related CN104765029B (zh) 2015-04-15 2015-04-15 一种基于射线追踪确定表面波导盲区信息的方法

Country Status (1)

Country Link
CN (1) CN104765029B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109900973B (zh) * 2019-02-25 2021-06-01 西安电子科技大学 基于规则矩形折射率单元的带有光强衰减的射线追踪方法
WO2020200454A1 (en) * 2019-04-04 2020-10-08 Siemens Industry Software And Services B.V. A method for computer-implemented simulation of radar raw data

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103135108A (zh) * 2011-11-23 2013-06-05 北京雷音电子技术开发有限公司 大气波导实时探测和诊断方法及终端
CN103278804B (zh) * 2013-06-14 2017-02-08 张永刚 一种波导超视距雷达

Also Published As

Publication number Publication date
CN104765029A (zh) 2015-07-08

Similar Documents

Publication Publication Date Title
US11688196B2 (en) Fish biomass, shape, and size determination
WO2021249071A1 (zh) 一种车道线的检测方法及相关设备
CN104200484B (zh) 基于云团特征分析的分布式光伏系统超短期出力预测方法
CN105241465B (zh) 一种道路更新的方法
CN106093946B (zh) 一种适用于场面监视雷达的目标凝聚方法及装置
CN103208115B (zh) 基于测地线距离的图像显著性区域检测方法
CN102938062B (zh) 一种基于内容的文档图像倾斜角估计方法
CN105974412B (zh) 一种用于合成孔径雷达的目标特征提取方法
CN104484868B (zh) 一种结合模板匹配和图像轮廓的运动目标航拍跟踪方法
CN102252623A (zh) 基于视频差异分析的输电线路导/地线覆冰厚度测量方法
CN110263713A (zh) 车道线检测方法、装置、电子设备及存储介质
CN106156758B (zh) 一种sar海岸图像中海岸线提取方法
CN102722706A (zh) 基于粒子滤波的红外弱小目标检测与跟踪方法及装置
CN102999926A (zh) 一种基于底层特征融合的图像视觉显著性计算方法
CN103500453B (zh) 基于伽玛分布和邻域信息的sar图像显著性区域检测方法
Zheng et al. A Decision‐Making Method for Ship Collision Avoidance Based on Improved Cultural Particle Swarm
CN104765029B (zh) 一种基于射线追踪确定表面波导盲区信息的方法
CN104535960A (zh) 一种基于rfid的室内快速定位方法
CN110852243A (zh) 一种基于改进YOLOv3的道路交叉口检测方法及装置
CN102722702A (zh) 多特征融合的粒子滤波视频目标跟踪方法
CN114781766B (zh) 水文站点的水文信息预测方法、装置、设备及存储介质
CN103268586B (zh) 一种基于热扩散理论的窗融合方法
US8374391B1 (en) Model of visibility of targets
CN112130166A (zh) 基于反光板网络的agv定位方法及装置
CN108629227A (zh) 在图像中确定车辆左右边界的方法及系统

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
EXSB Decision made by sipo to initiate 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: 20170905