CN105548971B - 一种针对海面蒸发波导超视距探测雷达的目标隐身方法 - Google Patents
一种针对海面蒸发波导超视距探测雷达的目标隐身方法 Download PDFInfo
- Publication number
- CN105548971B CN105548971B CN201610015373.5A CN201610015373A CN105548971B CN 105548971 B CN105548971 B CN 105548971B CN 201610015373 A CN201610015373 A CN 201610015373A CN 105548971 B CN105548971 B CN 105548971B
- Authority
- CN
- China
- Prior art keywords
- array
- evaporation
- target
- horizon
- parameter
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/38—Jamming means, e.g. producing false echoes
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
一种针对海面蒸发波导超视距探测雷达的目标隐身方法,步骤如下:测量欲隐身目标周围环境的大气参数,确定蒸发波导参数;根据射线方程计算蒸发波导中电波射线轨迹;计算可以实现蒸发波导超视距传播的边沿射线跳距,确定波束轨迹跳距,在欲隐身目标周围K倍波束轨迹跳距范围内部署相干声源阵列,通过声波干涉产生的声压分布,改变欲隐身目标周围大气的大气折射指数空间分布结构,使得超视距探测雷达的信号不能有效传播至欲探测目标并返回至探测雷达。本发明利用声波干涉原理,在海面蒸发波导环境中产生特定结构的声压分布,改变蒸发波导特征参数或者破坏蒸发波导结构,使得基于蒸发波导环境的超视距探测雷达系统无法有效探测目标,实现有效隐身。
Description
技术领域
本发明涉及雷达目标探测技术,特别涉及一种针对海面蒸发波导超视距探测雷达的目标隐身方法。
背景技术
大气波导是由于大气折射指数随高度迅速下降而形成的一种区别于标准大气的异常大气结构,它能够使电磁波射线向下弯曲的曲率大于地球表面的曲率,从而将电磁能量陷获在波导结构内形成大气波导传播。在贴近海面的几十米高度空间内,大气表面层内海水与大气的双重交互作用使得大气结构变化更加剧烈,蒸发波导就是频繁出现在海洋环境中的一种大气波导形式。
海面蒸发波导是由于海面水汽蒸发,使得在海面上很小高度范围内的大气湿度随高度锐减而形成的陷获层,它的形成与洋面上存在较强的垂直方向的湿度梯度有关。当大气折射指数N(N与折射率n的关系为N=(n-1)×106)在垂直方向的梯度满足dN/dh<-0.157N单位/m,或者当空间某点处修正折射率M(M与大气折射指数N的关系为M=N+0.157h,h为该点处的地面高度值)在垂直方向的梯度满足dM/dh<0时,电磁波在蒸发波导中传播时可能被海面蒸发波导层状结构所约束,出现超折射现象。海面蒸发波导一般发生在海洋大气环境几十米高度以下的海面大气环境中,具有稳定性好、持续时间长以及水平方向上延伸距离远等特点。
当海上雷达系统采用合适的频率和极化状态、且发射仰角满足一定条件时,雷达信号可沿海面蒸发波导传播,实现超视距探测海上目标。蒸发波导环境为雷达探测视距以外的目标提供了一种新途径,蒸发波导的出现大大地增加了电子对抗的灵活性,增加了电子侦察与反侦察、电子干扰与反干扰、电子欺骗与反欺骗、电子隐身与反隐身、电子摧毁与反摧毁等攻防手段。
由于受到地球曲率的影响,传统舰载雷达对海上目标的最远探测距离不超过40公里,利用空中预警机进行远距离目标追踪,往往成本高昂、维护费用高并且数量极为有限,而且由于空中预警机目标过于明显,容易受到威胁。但是,基于蒸发波导的超视距雷达的探测距离往往可以达到300公里以上,远远超过了视距雷达的探测范围。由于其成本低、精度高以及远距离探测等优点,蒸发波导超视距雷达探测应用愈来愈广泛。海面蒸发波导环境中超视距雷达的应用使得岸基和舰载微波雷达系统在海面蒸发波导环境条件下的探测能力明显增强。
躲避雷达探测的隐身技术与雷达探测技术是一对相互博弈的对抗技术,针对特定环境下雷达探测技术的隐身方法,对于增强特殊目标的生存能力具有重要的理论意义和实用价值。隐身技术通过降低特殊目标的可探测信息特征,使对方探测系统难以发现或者发现距离缩短。广义的隐身技术分为有源隐身技术和无源隐身技术两类。其中,有源隐身技术主要是指在被探测目标周围或附近设置假目标、利用地形匹配以及施放烟幕和利用光或电子干扰等手段隐蔽己方的真实目标。无源隐身技术主要是依靠减少目标的可探测信息特征,使敌方各种探测系统不能发现目标或发现概率减低的技术。
在雷达对抗过程中,海面蒸发波导环境中雷达超视距探测的应用越来越广泛,尤其在军事方面显示出强大的功能。同时,针对基于蒸发波导环境超视距雷达探测的目标隐身方法也成为学术界关注的焦点。目前,针对海面蒸发波导环境中超视距雷达探测的隐身方法主要聚焦于舰艇本身的设计,包括对军舰的外形结构、舰体材料、武器系统以及动力系统等进行隐身处理,大大减少了雷达反射/散射面积,使得对方超视距雷达尽可能小的探测到目标,这类隐身方式为无源隐身技术。此外,一些有源隐身技术也得到了很广泛的应用,比如喷雾自卫系统、释放烟幕弹、释放箔条等手段也可以减少对方超视距雷达探测到我方目标的真实位置。目前,尚无文献涉及通过干预蒸发波导环境而实现针对蒸发波导超视距雷达探测的隐身方法。
发明内容
本发明的目的是针对基于海面蒸发波导环境实现的超视距探测目标的雷达探测技术,提供一种通过人为扰动大气而破坏波导环境结构实现目标隐身的方法。
为了实现上述目的,本发明采取如下的技术解决方案:
一种针对海面蒸发波导超视距探测雷达的目标隐身方法,包括以下步骤:
S100、测量欲隐身目标周围环境的大气参数,从而确定蒸发波导参数;
S200、根据射线方程计算蒸发波导中电波射线轨迹;
S300、确定声源部署范围;计算可以实现蒸发波导超视距传播的边沿射线跳距,确定波束轨迹跳距,在欲隐身目标周围K倍波束轨迹跳距范围内部署相干声源阵列,K≥1;
S400、设定相干声源阵列的结构及声波频率,使得在欲隐身目标周围蒸发波导局部空间垂直剖面内出现ΔNC_C_E的极值,ΔNC_C_E为相干声波引起大气折射指数变化的有效值,步骤如下:
根据阵列结构获得声源阵列的阵列因子|S|关于空间位置参量的函数;
计算中央声压p0,在z=d的平面内搜索pi,j=0.35p0对应的位置,得到声压极大值级数坐标;
根据声压极大值级数坐标获取空间位置参量的搜索范围,在空间位置参量的搜索范围内搜索阵元间距、频率的取值,求解声源阵列的阵列因子|S|关于空间位置参量的导数或者偏导数为零的方程;
根据求解得到的声源阵列中阵元间距和声波频率确定阵列结构;
S500、设定声源阵列中阵元的声波功率,使pi,j=0.35p0对应位置声压的声压级SPL≥250dB。
在上述技术方案的基础上,所述波束轨迹跳距为由射线方程计算得到的电磁波波束的边沿射线中跳距中较大者,所述边沿射线的跳距为其相邻两次与海面发生相互作用点之间的距离。
在上述技术方案的基础上,所述步骤S300中的K≥3。
在上述技术方案的基础上,所述步骤S400中出现的ΔNC_C_E的极值的数量大于等于5。
在上述技术方案的基础上,所述相干声源阵列由多个半柱面波扩音器沿直线均匀间隔排列而成,声源阵列的阵列因子|S|表示为:其中,A为扩音器的数量,b为相邻扩音器间的距离,λ为声波的波长,θ为空间某点在x-z平面内的投影相对于坐标原点的位置矢量与z轴的夹角,β为两个相邻扩音器之间的初始相位差。
在上述技术方案的基础上,所述相干声源阵列由间隔布置的Ax×Ay个半球面波扩音器构成,声源阵列的阵列因子|S|表示为:其中,Ax为沿x轴设置的扩音器的数量,Ay为沿y轴设置的扩音器的数量,bx为x方向上相邻扩音器间的距离,by为y方向上相邻扩音器间的距离,为空间某点相对于坐标原点的位置矢量在x-y平面的投影与x轴的夹角,βx和βy分别为沿x方向和y方向的两个相邻扩音器之间的初始相位差。
由以上技术方案可知,本发明针对海面蒸发波导存在的情况下,利用相干声源阵列的声波干涉效应,在大气空间产生稳定的、特定结构的声压分布,改变海面蒸发波导中大气折射率垂直方向的梯度分布特性,从而引起大气折射率垂直方向梯度的变化,破坏海面蒸发波导结构,最终中断超视距雷达信号在海面蒸发波导中的传播,使得基于海面蒸发波导环境实现超视距探测目标的雷达,无法有效利用海面蒸发波导发现欲隐身的海上目标。
附图说明
图1为波导环境大气折射指数梯度结构及超视距传播示意图;
图2为超视距传播参数说明图;
图3为本发明方法的流程图;
图4a为确定边沿射线跳距的示意图;
图4b为确定声源阵列范围的示意图;
图5为本发明方法的示意图;
图6为蒸发波导结构折射指数梯度分布变化的效果示意图;
图7为设定声源阵列的结构及声波频率步骤的流程图;
图8为采用平面线阵的声源阵列的示意图;
图9为针对图8所示声源阵列的结构和频率优化示意图;
图10为采用矩形阵列的声源阵列的示意图;
图11为针对图10所示声源阵列的结构和频率优化示意图。
以下结合附图对本发明的具体实施方式作进一步详细地说明。
具体实施方式
下面结合附图对本发明进行详细描述,在详述本发明实施例时,为便于说明,表示器件结构的附图会不依一般比例做局部放大,而且所述示意图只是示例,其在此不应限制本发明保护的范围。需要说明的是,附图均采用非常简化的形式且均使用非精准的比例,仅用以方便、清晰地辅助说明本发明实施例的目的。
如图1所示:当海面上空一定空间范围内的大气折射指数N随高度变化的梯度满足dN/dh<-0.157N单位/m或者修正折射率M在垂直方向的梯度满足dM/dh<0时,该大气层结称为蒸发波导。蒸发波导的下界面为海面,上界面为大气分层面。如果无线电波的波长小于蒸发波导环境对应的临界波长λc,且发射角度满足特定的条件时,无线电波可以实现类似于金属波导内的传播模式而实现超视距传播。
如图2所示,水平极化波和垂直极化波在蒸发波导中的临界波长λc_h、λc_v与波导厚度d、天线架高d0和修正指数的梯度dM/dh有密切的关系。水平极化波在蒸发波导中的临界波长为垂直极化波在蒸发波导中的临界波长为除了dM/dh和无线电波的波长需要满足上述条件外,实现超视距传播还需要电磁波射线与平行于海面的平面之间的夹角需要小于最大临界角θc,φ表示天线主瓣的半功率波束宽度。
从上述电波传播理论可以看出,雷达信号要实现超视距传播的充要条件为:存在海面蒸发波导环境,即沿雷达信号传播路径的大气参数满足dN/dh<-0.157N单位/m或者dM/dh<0;雷达信号的工作波长小于极限波长;雷达天线主瓣波束轴与平行于海面的平面之间的夹角,即初始发射角小于最大临界角θc,其中的极限波长可根据波导环境特征参数和信号的极化形式得出,最大临界角根据天线的位置和波导环境特性参数得出。由此可见,基于蒸发波导环境实现的海上超视距雷达探测对波导环境参数依赖性非常强。
基于此,本发明方法的基本思路是:利用声波干涉原理,在海面蒸发波导环境中产生特定结构的声压分布,以引起大气折射率空间分布的变化,从而破坏蒸发波导环境中的大气折射指数梯度分布结构,通过改变欲隐身目标周围的大气参数来改变蒸发波导特征参数或者破坏蒸发波导结构,使得基于蒸发波导环境的超视距探测雷达系统无法有效探测目标,进而实现欲隐身目标的有效隐身。
由于声波是弹性纵波,所以声波传播过程中会引起对流层局部大气产生周期运动,使得局部大气的压强发生周期性变化。而大气折射指数与大气压强存在密切的关系,所以声波扰动可以改变对流层局部大气折射指数,如果声波在局部空间发生干涉,那么声压在干涉区域将出现有规律的空间分布,即出现压强起伏最大和压强起伏最小的稳定分布,所以在干涉区域即出现了人为可以控制的折射指数“不均匀体”,大气折射指数“不均匀体”的空间分布特性由声源阵列的阵列因子|S|决定,这些不均匀体的形状和尺寸可以通过声波的频率和相干声源的结构参数来控制,本发明正是利用以上相干声波改变大气折射指数空间分布结构的理论机理来实现目标隐身。
以上是本发明的核心思想,下面将结合本发明实施例中的附图,对本发明实施例的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下获得的所有其他实施例,都属于本发明保护的范围。
在下面的描述中阐述了很多具体细节以便于充分理解本发明,但是本发明还可以采用其它不同于在此描述的其它方式来实施,本领域技术人员可以在不违背本发明内涵的情况下做类似推广,因此本发明不受下面公开的具体实施例的限制。
图3为本发明方法的流程图,下面结合图3,对本发明的目标隐身方法进行说明。本发明目标隐身方法通过在欲隐身的海上目标周围设置相干声源阵列,利用相干声源阵列发出的相干声波的干涉效应,改变大气折射率的空间结构分布,从而破坏形成海面蒸发波导的必要条件,截断海面蒸发波导中超视距雷达信号的传播,最终使得欲隐身海上目标有效隐身。其具体步骤如下:
S100、测量欲隐身目标周围环境的大气参数,大气参数包括大气温度、湿度和压强,从而确定蒸发波导参数,蒸发波导参数包括蒸发波导的高度、强度和厚度;
S200、根据射线方程计算蒸发波导中电波射线轨迹;
根据蒸发波导参数,由射线方程计算电磁波射线在蒸发波导环境中的传播射线轨迹,其中,n为大气折射率,l为射线上某点相对射线起点的轨迹长度,为射线上某点处位置相对于射线起点的位置矢量;可采用传统积分法、泰勒级数近似法和混合法对射线方程进行求解,计算波导环境中射线轨迹时混合法更适用;
S300、计算可以实现蒸发波导超视距传播的边沿射线跳距,确定波束轨迹跳距,从而确定声源部署范围,在欲隐身目标周围K倍波束轨迹跳距范围内部署相干声源阵列,K≥1;优选的,K≥3;
假设由射线方程计算得到的电磁波波束中边沿射线的轨迹如图4a所示,边沿射线1的跳距为其相邻两次与海面发生相互作用点之间的距离,即A、B两点之间的“虚线”长度,边沿射线2的跳距为其相邻两次与海面发生相互作用点之间的距离,即C、D两点之间的“虚线”长度,将边沿射线1和边沿射线2的跳距中较大者定义为波束轨迹跳距s,则在欲隐身目标周围K倍波束轨迹跳距s范围内部署相干声源阵列,如图4b所示,在图4b中4个箭头方向的K·s长度范围内部署相干声源;当边沿射线指射线仰角小于临界仰角θc的最大仰角对应的射线时,可以认为边沿射线是射线仰角等于θc的电波射线;
S400、设定相干声源阵列的结构及声波频率,使得在欲隐身目标周围蒸发波导局部空间垂直剖面内出现w个ΔNC_C_E的极值,ΔNC_C_E为相干声波引起大气折射指数变化的有效值;优选的,w≥5;
如图5和图6所示,本发明通过在欲隐身目标周围部署相干声波源,利用声波干涉产生的声压分布,改变欲隐身目标周围大气的大气折射指数空间分布结构,破坏无线电波实现超视距传播所需充要条件中的“dN/dh<-0.157N单位/m或者dM/dh<0”,使得超视距探测雷达的信号不能有效传播至欲探测目标并且返回至探测雷达,实现目标隐身;
相干声波扰动下,空间某点处相干声波引起的声压有效值p=0.707·ps|S|,与声压有效值对应的大气折射指数起伏的有效值其中,T为环境温度,P为大气压强,e为水汽压,ps为声源阵列中的阵元在空间某点处位置的声压,|S|为声源阵列的阵列因子,ps=pAexp(-α),pA为每个阵元表面处的声压幅度、α为从声源表面至空间某点处位置的声压衰减值,pA和α可以根据声波辐射理论求得,大气折射指数起伏的有效值ΔNC_C_E的极值(极大值或者极小值)由声源阵列的阵列因子|S|空间分布函数决定,ΔNC_C_E的极值出现时,声源阵列的阵列因子|S|关于空间位置参量θ、的导数或者偏导数为零;因此,通过求解方程可以确定声源阵列的结构和声波频率;
如图7所示,设定声源阵列的结构及声波频率包括以下子步骤:
a.根据阵列结构获得声源阵列的阵列因子|S|关于空间位置参量的函数;由声阵列理论可知,根据不同的阵列结构,声源阵列的阵列因子|S|是关于空间位置参量θ或者关于θ和的函数,当声源阵列为线阵时,|S|是关于空间位置参量θ的函数,当声源阵列为除了线阵之外的其它结构时,|S|是关于空间位置参量θ和的函数;
b.计算中央声压p0,在z=d的平面内搜索声压pi,j=0.35p0对应的位置,得到声压极大值级数坐标,其中,d为蒸发波导厚度;|S|取最大值时对应的位置就是产生中央声压p0的位置,声源阵列相对于该位置呈对称分布状态,p0=0.707·ps0|S|,ps0为声源阵列中的阵元在|S|为最大值时位置的声压;当|S|是关于空间位置参量θ的函数时,pi,j表示在z=d的平面内由于θ变化使得|S|出现的第j个极大值,当|S|是关于空间位置参量θ和的函数时,pi,j表示在z=d的平面内由于θ变化使得|S|出现的第j个极大值、由于变化使得|S|出现的第i个极大值;
c.根据声压极大值级数坐标获取空间位置参量的搜索范围,在空间位置参量的搜索范围内搜索阵元间距b、频率f的取值,求解方程或方程组空间位置参量θ的搜索范围是[θ1,π/2),θ1=θi,j,θi,j为pi,j=0.35p0对应位置的θ取值,在同一个θ取值条件下空间位置参量的搜索范围为[0,2π);
d.根据求解得到的声源阵列中阵元间距b和声波频率f确定阵列结构;
S500、设定声源阵列中阵元的声波功率,使pi,j=0.35p0对应位置声压的声压级SPL≥250dB。
声源阵列结构不同,声源阵列的阵列因子|S|也不同,通过调整阵列的阵元间距、阵元声波频率的取值,使得在欲隐身目标周围蒸发波导局部空间垂直剖面内出现w个ΔNC_C_E的极大值或者极小值;下面以两个例子对声源阵列的结构及声波频率的设定步骤进行说明;
当声源阵列采用如图8所示的平面线阵时,该声源阵列由多个半柱面波扩音器沿直线均匀间隔排列而成,由干涉理论可知,声源阵列的阵列因子|S|表示为:此时,|S|是关于空间位置参量θ的函数,其中,A为扩音器的数量,b为相邻扩音器间的距离,λ为声波的波长,θ为空间某点在x-z平面内的投影相对于坐标原点的位置矢量与z轴的夹角(图8),β为两个相邻扩音器之间的初始相位差;
即对于该平面线阵,声源阵列的结构及声波频率的设定通过方程确定,计算中央声压,搜索声压pi,j=0.35p0对应的位置,得到声压极大值级数坐标(j),如图9所示,根据z=d的平面内声压周期分布结果,得到波导环境垂直剖面内空间位置参量θ的搜索范围[θ1,π/2),在该搜索范围内搜索阵元间距b、频率f的取值,使得b/λ的取值满足方程从而得到声源阵列的结构及各阵元的声波频率。优选的,在θ∈[θ1,π/2)条件下解的个数w≥5。
当声源阵列采用如图10所示的平面矩形阵列时,该声源阵列由间隔布置的Ax×Ay个半球面波扩音器构成,声源阵列的阵列因子|S|表示为:此时,|S|是空间位置参量θ、的函数,其中,Ax为沿x轴设置的扩音器的数量,Ay为沿y轴设置的扩音器的数量,bx为x方向上相邻扩音器间的距离,by为y方向上相邻扩音器间的距离,θ为空间某点在x-z平面内的投影相对于坐标原点的位置矢量与z轴的夹角,为空间某点相对于坐标原点的位置矢量在x-y平面的投影与x轴的夹角(图10),βx和βy分别为沿x方向和y方向的两个相邻扩音器之间的初始相位差;
即对于该平面矩形阵列,声源阵列的结构及声波频率的设定通过方程组确定,计算中央声压,搜索声压pi,j=0.35p0对应的位置,得到声压极大值级数坐标(i,j),如图11所示,根据z=d的平面内声压周期分布结果,得到波导环境垂直剖面内空间位置参量θ的搜索范围[θ1,π/2),在同一个θ取值条件下另一空间位置参量的搜索范围为[0,2π),在该两个空间位置参量的搜索范围内搜索阵元间距bx、by和频率f的取值,使得bx/λ、by/λ的取值满足方程组从而得到声源阵列的结构及各阵元的声波频率。优选的,在θ∈[θ1,π/2)和条件下解的个数w≥5;
设定声波功率使得pi,j=0.35p0对应位置声压的声压级SPL≥250dB。
本发明提供的目标隐身方法通过在欲隐身目标周围部署一定结构的相干声波源,通过声波干涉产生的声压分布,改变欲隐身目标周围大气的大气折射指数空间分布结构,破坏无线电波实现超视距传播所需的条件,使得探测雷达的信号不能有效传播至欲探测目标并且返回至探测雷达,是一种全新的隐身方法,相比于已有的隐身方法,本发明操作更灵活、作用范围更广、应用成本更低,尤其是可以配合现有的尚不具备隐身能力的舰艇,大大减小了这类舰艇受到的来自水面的威胁,从而可以使己方舰艇发挥更佳的作战性能,是一种更可靠并且实用的针对海面蒸发波导环境中超视距雷达探测的隐身方法,有望成为隐身技术领域的一项重大突破。
当然,本发明的技术构思并不仅限于上述实施例,还可以依据本发明的构思得到许多不同的具体方案,例如,声源阵列的结构除了采用均匀平面线阵、均匀平面矩形阵以外,还可以采用平面环形阵列、非均匀平面线阵列、非均匀平面矩形阵列、立体阵列,诸如此等改变以及等效变换均应包含在本发明所述的范围之内。
Claims (6)
1.一种针对海面蒸发波导超视距探测雷达的目标隐身方法,其特征在于,包括以下步骤:
S100、测量欲隐身目标周围环境的大气参数,从而确定蒸发波导参数;
S200、根据射线方程计算蒸发波导中电波射线轨迹;
S300、确定声源部署范围;计算可以实现蒸发波导超视距传播的边沿射线跳距,确定波束轨迹跳距,在欲隐身目标周围K倍波束轨迹跳距范围内部署相干声源阵列,K≥1;
S400、设定相干声源阵列的结构及声波频率,使得在欲隐身目标周围蒸发波导局部空间垂直剖面内出现ΔNC_C_E的极值,ΔNC_C_E为相干声波引起大气折射指数变化的有效值,步骤如下:
根据阵列结构获得声源阵列的阵列因子|S|关于空间位置参量的函数;
计算中央声压p0,在z=d的平面内搜索pi,j=0.35p0对应的位置,得到声压极大值级数坐标,d为蒸发波导厚度,z为蒸发波导厚度方向的坐标,pi,j为搜索声压;
根据声压极大值级数坐标获取空间位置参量的搜索范围,在空间位置参量的搜索范围内搜索阵元间距、频率的取值,求解声源阵列的阵列因子|S|关于空间位置参量的导数或者偏导数为零的方程;
根据求解得到的声源阵列中阵元间距和声波频率确定阵列结构;
S500、设定声源阵列中阵元的声波功率,使pi,j=0.35p0对应位置声压的声压级SPL≥250dB。
2.如权利要求1所述的针对海面蒸发波导超视距探测雷达的目标隐身方法,其特征在于:所述波束轨迹跳距为由射线方程计算得到的电磁波波束的边沿射线中跳距中较大者,所述边沿射线的跳距为其相邻两次与海面发生相互作用点之间的距离。
3.如权利要求1所述的针对海面蒸发波导超视距探测雷达的目标隐身方法,其特征在于:所述步骤S300中的K≥3。
4.如权利要求1所述的针对海面蒸发波导超视距探测雷达的目标隐身方法,其特征在于:所述步骤S400中出现的ΔNC_C_E的极值的数量大于等于5。
5.如权利要求1所述的针对海面蒸发波导超视距探测雷达的目标隐身方法,其特征在于:所述相干声源阵列由多个半柱面波扩音器沿直线均匀间隔排列而成,声源阵列的阵列因子|S|表示为:其中,A为扩音器的数量,b为相邻扩音器间的距离,λ为声波的波长,θ为空间某点在x-z平面内的投影相对于坐标原点的位置矢量与z轴的夹角,β为两个相邻扩音器之间的初始相位差。
6.根据权利要求1所述的针对海面蒸发波导超视距探测雷达的目标隐身方法,其特征在于:所述相干声源阵列由间隔布置的Ax×Ay个半球面波扩音器构成,声源阵列的阵列因子|S|表示为:其中,Ax为沿x轴设置的扩音器的数量,Ay为沿y轴设置的扩音器的数量,bx为x方向上相邻扩音器间的距离,by为y方向上相邻扩音器间的距离,λ为声波的波长,θ为空间某点在x-z平面内的投影相对于坐标原点的位置矢量与z轴的夹角,为空间某点相对于坐标原点的位置矢量在x-y平面的投影与x轴的夹角,βx和βy分别为沿x方向和y方向的两个相邻扩音器之间的初始相位差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610015373.5A CN105548971B (zh) | 2016-01-11 | 2016-01-11 | 一种针对海面蒸发波导超视距探测雷达的目标隐身方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610015373.5A CN105548971B (zh) | 2016-01-11 | 2016-01-11 | 一种针对海面蒸发波导超视距探测雷达的目标隐身方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105548971A CN105548971A (zh) | 2016-05-04 |
CN105548971B true CN105548971B (zh) | 2017-12-15 |
Family
ID=55828283
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610015373.5A Active CN105548971B (zh) | 2016-01-11 | 2016-01-11 | 一种针对海面蒸发波导超视距探测雷达的目标隐身方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105548971B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2766934C1 (ru) * | 2020-09-24 | 2022-03-16 | Акционерное общество "Научно-производственное предприятие "Салют" (АО "НПП "Салют") | Мобильная береговая радиолокационная станция загоризонтного обнаружения с повышенной скрытностью излучения |
CN112213707A (zh) * | 2020-10-14 | 2021-01-12 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | 一种蒸发波导环境下舰载脉冲雷达作用范围的预测方法 |
CN113255190B (zh) * | 2021-06-07 | 2023-06-20 | 西安电子科技大学 | 一种不规则非均匀介质体雷达散射截面的计算方法 |
CN118050728B (zh) * | 2024-04-16 | 2024-06-25 | 中国水利水电第十四工程局有限公司 | 一种用于航道安全监控的目标获取方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5040157A (en) * | 1990-10-15 | 1991-08-13 | The United States Of America As Represented By The Secretary Of The Navy | Expendable virtual vertical sensing array |
CN101995574A (zh) * | 2010-11-03 | 2011-03-30 | 中国科学院声学研究所 | 一种近场聚焦波束形成定位法 |
CN102288970A (zh) * | 2011-07-13 | 2011-12-21 | 中国人民解放军海军工程大学 | 一种探测大气波导环境的方法、系统及探测机 |
CN102540162A (zh) * | 2011-12-12 | 2012-07-04 | 中国船舶重工集团公司第七二四研究所 | 基于海杂波的低空电磁波传播特性估计方法 |
CN104678384A (zh) * | 2013-11-28 | 2015-06-03 | 中国科学院声学研究所 | 一种波束域的声压差互相关谱分析水下目标速度估计方法 |
-
2016
- 2016-01-11 CN CN201610015373.5A patent/CN105548971B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5040157A (en) * | 1990-10-15 | 1991-08-13 | The United States Of America As Represented By The Secretary Of The Navy | Expendable virtual vertical sensing array |
CN101995574A (zh) * | 2010-11-03 | 2011-03-30 | 中国科学院声学研究所 | 一种近场聚焦波束形成定位法 |
CN102288970A (zh) * | 2011-07-13 | 2011-12-21 | 中国人民解放军海军工程大学 | 一种探测大气波导环境的方法、系统及探测机 |
CN102540162A (zh) * | 2011-12-12 | 2012-07-04 | 中国船舶重工集团公司第七二四研究所 | 基于海杂波的低空电磁波传播特性估计方法 |
CN104678384A (zh) * | 2013-11-28 | 2015-06-03 | 中国科学院声学研究所 | 一种波束域的声压差互相关谱分析水下目标速度估计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105548971A (zh) | 2016-05-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105548971B (zh) | 一种针对海面蒸发波导超视距探测雷达的目标隐身方法 | |
CN108318032A (zh) | 一种考虑攻防对抗的无人机航迹智能规划方法 | |
Wu et al. | Ray tracing based wireless channel modeling over the sea surface near Diaoyu islands | |
CN108388732B (zh) | 掠海目标多径散射特性仿真计算方法及系统 | |
CN110031801A (zh) | 一种高速临近空间飞行器的探测方法及系统 | |
CN108089158A (zh) | 一种全向雷达角反射器阵列 | |
Hu et al. | Solution and optimization of aircraft swarm cooperating anti-stealth formation configuration | |
Wei et al. | A model for calculating electromagnetic scattering from target in evaporation duct | |
Gaoqian et al. | Radar netting technology & its development | |
CN114089304A (zh) | 一种动态交会状态目标近场散射回波的求解方法 | |
Chen et al. | 3D representation of radar coverage in complicated environment | |
Chu et al. | Acoustic Scattering Characteristics of Underwater One-dimensional Corner Reflector Array | |
Niu et al. | The Array Optimization Model for Distributed Antisubmarine of Non-cluster USV Based on Coverage | |
Papastamatis et al. | Tropospheric Propagation Modeling in the Modern Air Defense Environment | |
Tian et al. | Anti‐multipath interference method of semi‐active seek | |
Zou et al. | Using radar echo to confront anti-radiation missiles | |
Hou et al. | RCS Reduction Based on Optimization Algorithms and Novel Polarization Conversion Meta-surface | |
Matuszewski | Method of radiolocation object shield zone calculation for ground jammer stations | |
Li et al. | Three-dimensional Fitting Method for Scattering Characteristics of Objects | |
Hou et al. | Study on the controlled atmospheric refractivity fluctuation impairing the imaging quality of the heterodyne detection lidar | |
Zou et al. | Modeling and Simulation of Anti-ship Missile Passive Jamming Signal | |
Wang et al. | Reach on Dynamic EME Adaptability in Battlefield Based on System Synergy | |
Zhang et al. | Single station positioning of identification friend or foe signal | |
He et al. | Reflection on the Makeup Protection of the Target of People’s Livelihood Under the Condition of Informatization | |
Li et al. | A Method for Parallel Calculation of Radar Detection Capability Based on 3D Subdivision Grid. |
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 |