CN111539106A - 基于参数化建模数据的电子侦察装备输出数据仿真方法 - Google Patents

基于参数化建模数据的电子侦察装备输出数据仿真方法 Download PDF

Info

Publication number
CN111539106A
CN111539106A CN202010328443.9A CN202010328443A CN111539106A CN 111539106 A CN111539106 A CN 111539106A CN 202010328443 A CN202010328443 A CN 202010328443A CN 111539106 A CN111539106 A CN 111539106A
Authority
CN
China
Prior art keywords
scanning
frequency
time
reconnaissance equipment
electronic reconnaissance
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
CN202010328443.9A
Other languages
English (en)
Other versions
CN111539106B (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.)
Chengdu Zhongxiangtiandi Network Technology Co Ltd
Original Assignee
Chengdu Zhongxiangtiandi Network Technology 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 Chengdu Zhongxiangtiandi Network Technology Co Ltd filed Critical Chengdu Zhongxiangtiandi Network Technology Co Ltd
Priority to CN202010328443.9A priority Critical patent/CN111539106B/zh
Publication of CN111539106A publication Critical patent/CN111539106A/zh
Application granted granted Critical
Publication of CN111539106B publication Critical patent/CN111539106B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • 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
    • 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
    • G01S7/4004Means for monitoring or calibrating of parts of a radar system
    • G01S7/4021Means for monitoring or calibrating of parts of a radar system of receivers

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了基于参数化建模数据的电子侦察装备输出数据仿真方法,通过计算电子侦察装备输出全脉冲数据,同时加入误差模型和异常模型,让计算机软件仿真输出的全脉冲数据更符合真实装备输出的全脉冲数据;通过此种方法,对电子侦察装备的输出数据进行仿真,可以支持进行电子侦察装备效能仿真与评估、电子侦察装备软件的测试与验证、电子侦察装备的模拟训练,节约实装系统的费用,降低研制周期,解决了装备输出的反馈各不相同,很难达到数据统一标准的效果,在进行仿真输出时,可能因为多个系统不兼容,延迟高,速度缓慢,现有的仿真结果与真实电子侦察装备侦察输出结果存在较大差异的问题。

Description

基于参数化建模数据的电子侦察装备输出数据仿真方法
技术领域
本发明涉及一种输出数据仿真方法,具体涉及基于参数化建模数据的电子侦察装备输出数据仿真方法。
背景技术
随着信息技术的迅猛发展和广泛应用,电子侦察装备作为装备体系作战能力的“倍增器”,已经成为装备体系中最具活力的军事装备。基于信息系统的体系对抗已成为信息化条件下一种主要的战争形态,电子侦察装备的信息化、智能化特征及其对装备体系的融合特性,对电子侦察装备的试验与鉴定也提出了新的要求,世界各国越来越重视电子侦察装备实战环境中的作战试验鉴定,特别要求在逼真作战环境下充分、全面考核电子侦察装备作战效能、作战适用性以及生存能力等,以期在部署前发现电子侦察装备实际作战要求方面存在的问题和一些重要缺陷。
只有在满足战场电磁环境模拟逼真性要求条件下得出的电子装备作战效能、作战适用性等评估结论才会有一定的指导意义,不够逼真的战场电磁环境设计可能会毁掉整个作战试验。国内对武器装备战场复杂电磁环境适应性试验和作战试验的研究与实践起步较晚,相应地对战场电磁环境的构建方法与逼真性研究成为最近几年的研究热点。
现有的电子侦察装备的模型比较简单,主要对工作频率范围、灵敏度进行建模,能支持效能仿真,目前所采用的数据仿真方法,通常是军方用的软件系统,这些软件系统会接入现役的电子侦察装备,对该装备的截获数据和输出数据进行仿真模拟,但是接入每一个装备需要对应开发相应配套的软件,开发难度很大,并且装备输出的反馈各不相同,很难达到数据统一标准的效果,在进行仿真输出时,其因为多个系统不兼容的问题,延迟也高,速度很缓慢,现有的仿真结果与真实电子侦察装备侦察输出结果存在较大差异。
发明内容
本发明所要解决的技术问题是接入不同装备需要对应开发相应配套的软件,开发难度很大,未考察电子侦察设备的测量误差情况,与真实电子侦察装备侦察输出结果存在较大差异,并且电子侦察装备模型不通用,输出延迟高,速度慢,本申请文件提供基于参数化建模数据的电子侦察装备输出数据仿真方法,解决上述的问题。
本发明通过下述技术方案实现:
基于参数化建模数据的电子侦察装备输出数据仿真方法,其特征在于,所述方法包括如下步骤:S1:读取电子侦察设备建模数据和雷达建模数据后对电子侦察波束扫描仿真,实时计算电子侦察设备的波束指向;S2:随后对电子侦察设备的工作频率进行仿真,通过侦察设备频率范围工作带宽、扫描时间、跟踪时间、重点频率集和实时仿真时间参数,实时计算当前电子侦察设备的工作频率参数;S3:在对电子侦察截获关系进行仿真,通过能量截获计算,空域截获关系计算,频率截获关系计算获得侦察设备实时的接收灵敏度、雷达扫描波束、频率,然后判断雷达接收灵敏度、雷达扫描波束、频率是否在侦察设备实时扫描频率范围内;S4:通过导入误差计算量进行全脉冲数据仿真,根据截获的目标雷达参数,对目标雷达的全脉冲波形参数进行仿真,包括全脉冲脉宽、脉冲间隔、脉冲数量、脉幅、脉冲串间隔、脉内调制类型进行仿真;S5:导入异常介入参数后重新进行全脉冲数据仿真,将电子侦察设备建模数据中读取电子侦察设备异常频率、异常脉宽、异常重频、异常方位、异常时间参数以及异常类型导入,对仿真数据进行比对。
在现有技术中,大多都是仅进行截获,而未对截获后的数据进行误差和异常数据的排除,为了避免出现较大的误差,在截获数据后,需要首先对电子侦察设备的测量误差进行对比,让其余截获的数据在一定的匹配范围内才会进行输出,这样避免输出结果与实际截获数据存在误差较大的情况。
进一步地,所述步骤S1中实时计算电子侦察设备的波束指向的方式如下:A1)断扫描方式,如果是扇扫则进步骤b,如果是圆扫则进步骤c,如果是固定则进步骤d;A2)判断扫描方向,如果是顺时针扫描,则进步骤e,否则进步骤f;A3)判断扫描方向,如果是顺时针扫描,则进步骤g,否则进步骤h;A4)天线波束扫描方位,天线波束扫描俯仰;A5)根据扫描时间计算扫描过的角度,angle=speed*time,计算扫描宽度width_deg=e_angle-s_angle,计算扫描过多少个扫描宽度cnt=floor(angle/width_deg),若cnt对2取余为0,则天线波束扫描方位=s_angle+(angle-cnt*width_deg),否则天线波束扫描方位=e_angle-(angle-cnt*width_deg);A6)根据扫描时间计算扫描过的角度,angle=speed*time,计算扫描宽度width_deg=s_angle-e_angle,计算扫描过多少个扫描宽度cnt=floor(angle/width_deg),如果cnt对2取余为0,则天线波束扫描方位=s_angle-(angle-cnt*width_deg),否则天线波束扫描方位=e_angle+(angle-cnt*width_deg);A7)天线波束扫描方位=s_angle+angle,再对扫描方位校验,换算到0-360°范围内;A8)天线波束扫描方位=s_angle-angle,再对扫描方位校验,换算到0-360°范围内;A9)天线波束扫描俯仰=c_rl+el。
进一步地,所述步骤S2中实时计算当前电子侦察设备的工作频率参数的方式如下:B1)判断输入的重点频率集中的频率是否超出扫描范围,超出的则抛弃;B2)然后计算基于工作带宽和最大最小频率的扫描次数cnt=floor((maxFreq-minFreq)/bandWidth);B3)判断fomd((maxFreq-minFreq),bandWidth)是否大于零,如果大于零则cnt加1;B4)计算顺序扫描需要时间normalTime=cnt*scanTime,计算跟踪需要时间followTime=size(freqs)*followTime;B5)计算当前仿真时间处于顺序扫描时间段还是跟踪扫描时间段,realTime=fmod(time,normalTime+followTime),如果realTime小于normalTime,则处于顺序扫描阶段,否则处于跟踪扫描阶段;B6)在正常扫描阶段,计算实时中心频率midFreq=minFreq+bandWidth/2+floor(realTime/scanTime)*bandwidth;B7)在跟踪扫描阶段,计算实时中心频率midFreq=freqs[floor((realTime-normalTime)/followTime)];B8)实时频率最小值realMinFreq=midFreq-bandWidth/2,如果小于最小频率,则使用最小频率;B9)实时频率最大值realMaxFreq=midFreq+bandwidth/2,如果大于最大频率,则使用最大频率。
进一步地,所述步骤S3中能量截获计算根据场景仿真中自身平台姿态信息和其他平台的姿态信息,参考雷达的模型数据、电子侦察设备的模型数据、传播模型,计算电子侦察设备的接收功率,并与接收灵敏度进行比较。
进一步地,所述步骤S3中空域截获关系计算根据作战场景中目标的位置、电子侦察设备的位置,目标雷达扫描波束范围,电子侦察设备波束扫描范围,计算电子侦察设备空域截获关系,通过计算相对天线方位俯仰算法,再基于天线波束宽度以及天线波束当前方位、俯仰,计算目标相对载体的方位、俯仰是否在波束扫描范围内,在范围内,则截获成功,否则截获失败。
进一步地,所述步骤S3中频率截获关系计算根据电子侦察设备工作模式,参考电子侦察设备模型参数,结合雷达模型的数据,计算电子侦察设备空域截获关系,计算出侦察设备实时的频率,然后判断雷达频率是否在侦察设备实时扫描频率范围内。
本发明与现有技术相比,具有如下的优点和有益效果:
1、本发明基于参数化建模数据的电子侦察装备输出数据仿真方法,本发明考虑了在计算电子侦察装备输出全脉冲数据,同时加入误差模型和异常模型,让计算机软件仿真输出的全脉冲数据更符合真实装备输出的全脉冲数据;通过此种方法,对电子侦察装备的输出数据进行仿真,可以支持进行电子侦察装备效能仿真与评估、电子侦察装备软件的测试与验证、电子侦察装备的模拟训练,节约实装系统的费用,降低研制周期;
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施例,对本发明作进一步的详细说明,本发明的示意性实施方式及其说明仅用于解释本发明,并不作为对本发明的限定。
实施例
本发明基于参数化建模数据的电子侦察装备输出数据仿真方法,其特征在于,所述方法包括如下步骤:S1:读取电子侦察设备建模数据和雷达建模数据后对电子侦察波束扫描仿真,实时计算电子侦察设备的波束指向;S2:随后对电子侦察设备的工作频率进行仿真,通过侦察设备频率范围工作带宽、扫描时间、跟踪时间、重点频率集和实时仿真时间参数,实时计算当前电子侦察设备的工作频率参数;S3:在对电子侦察截获关系进行仿真,通过能量截获计算,空域截获关系计算,频率截获关系计算获得侦察设备实时的接收灵敏度、雷达扫描波束、频率,然后判断雷达接收灵敏度、雷达扫描波束、频率是否在侦察设备实时扫描频率范围内;S4:通过导入误差计算量进行全脉冲数据仿真,根据截获的目标雷达参数,对目标雷达的全脉冲波形参数进行仿真,包括全脉冲脉宽、脉冲间隔、脉冲数量、脉幅、脉冲串间隔、脉内调制类型进行仿真;S5:导入异常介入参数后重新进行全脉冲数据仿真,将电子侦察设备建模数据中读取电子侦察设备异常频率、异常脉宽、异常重频、异常方位、异常时间参数以及异常类型导入,对仿真数据进行比对。
在步骤S5中,考虑测量误差情况下的全脉冲数据仿真:根据截获的目标雷达参数,对目标雷达的全脉冲波形参数进行仿真,包括全脉冲脉宽、脉冲间隔、脉冲数量、脉幅、脉冲串间隔、脉内调制类型进行仿真。
从模型建模数据中读出电子侦察设备频率测量误差ERR_R_RF、脉宽测量误差ERR_R_PW、重频测量误差ERR_R_PRI、方位测量误差ERR_R_DOA、时间测量误差ERR_R_TIME;从雷达建模数据中得到雷达脉冲参数,包括工作频率PULSE_RF、脉宽PULSE_PW、重频PULSE_PRI、脉内类型PULSE_MN;从雷达实时仿真中得到雷达的位置LonT、LatT,HT;从电子支援仿真中得到电子支援的位置LonR、LatR,HR从仿真控制可获得当前仿真时间T_SIM;
仿真输出全脉冲的参数如下:输出的全脉冲频率=PULSE_RF+random(0,ERR_R_RF);
输出的全脉冲脉宽=PULSE_PW+random(0,ERR_R_PW);
输出的全脉冲重频=PULSE_PRI+random(0,ERR_R_PRI);
输出的全脉冲方位=方位(根据方位计算公式)+random(0,ERR_R_DOA);
输出的全脉冲时间=T_SIM(精确到ns)+10e9*(雷达与电子支援设备的距离,公里)/C(公里/秒)+random(0,ERR_R_TIME);
如果需要考虑异常数据考虑异常数据情况下的全脉冲数据仿真:从模型建模数据中读取电子侦察设备异常频率、异常脉宽、异常重频、异常方位、异常时间参数以及异常类型(包括单个异常以及多异常类型)。
从模型建模数据中读出电子侦察设备异常频率测量误差ERR_RF、异常脉宽测量误差ERR_PW、异常重频测量误差ERR_PRI、异常方位测量误差ERR_DOA、异常时间测量误差ERR_TIME;读取异常脉冲数量num(包括单一异常和多个异常),异常脉冲数肯定小于全脉冲数;随机生成num个小于全脉冲数量的异常下标[n1,n2…,n];
indexs=rand([1,全脉冲数量],1,num);根据indexs矩阵修改对应下标的脉冲数据,输出异常脉冲的参数如下:
输出的异常脉冲频率=异常频率+random(0,ERR_RF);
输出的异常脉冲脉宽=异常脉宽+random(0,ERR_PW);
输出的异常脉冲重频=异常重频+random(0,ERR_PRI);
输出的异常脉冲方位=异常方位(根据方位计算公式)+random(0,ERR_DOA);输出的异常脉冲时间=T_SIM(精确到ns)+10e9*(雷达与电子支援设备的异常距离,
公里)/C(公里/秒)+random(0,ERR_TIME);
组合异常脉冲,异常脉冲可以单个属性异常或者多个属性异常,如可以异常脉冲中只有频率异常,也可以是频率异常加脉宽异常,以及其他各种情况,可任意组合。
所述步骤S3中能量截获计算根据场景仿真中自身平台姿态信息和其他平台的姿态信息,参考雷达的模型数据、电子侦察设备的模型数据、传播模型,计算电子侦察设备的接收功率,并与接收灵敏度进行比较。在进行能量截获计算时,其具体步骤为基于WGS84将侦察设备和目标雷达的经纬高(lon,lat,height)转换为空间坐标。
赤道半径radiusEquator=6378137.0,极地半径radiusPolar=6356752.3142;
Flattening=(radiusEquator-radiusPolar)/radiusEquator;
eccentricitySquared=2*Flattening–Flattening^2;
N=radiusEquator/(sqrt(1.0-eccentricitySquared*sind(lat)*sind(lat)));
X=(N+height)*cosd(lat)*cosd(lon);
Y=(N+height)*cosd(lat)*sind(lon);
Z=(N*(1.0-eccentricitySquared)+height)*sind(lat);
根据空间两点坐标计算两点距离,侦察设备坐标(x1,y1,z1)和雷达坐标(x2,y2,z2):
Distance=sqrt(pow(x1-x2,2)+pow(y1-y2,2)+pow(z1-z2,2));分别计算侦察设备和雷达设备的天线实时方位俯仰;根据侦察设备和雷达设备所属平台的姿态信息(方位az、俯仰el、横滚rl),分别计算坐标转换矩阵;
Matrix=[cos(rl),0,-sin(rl);sin(el)*sin(rl),cos(el),sin(el)*sin(rl);cos(el)*sin(rl),-sin(el),cos(el)*cos(rl)];
根据侦察设备和雷达设备所属平台的位置信息(经度、纬度、高度),分别将雷达转换到侦察设备载体系中,和把侦察设备转到雷达载体系中;先将目标转换到载体水平系中
Xh=(-X)*cos(az)+(-Y)*(-sin(az));
Yh=(-X)*sin(az)+(-Y)*cos(az);
Zh=Z;
再将其转换到载体系中,[Xa,Ya,Za]=Matrix*[Xh,Yh,Zh],展开矩阵,可得:
Xa=cos(rl)*Xh-sin(rl)*Zh;
Ya=sin(el)*sin(rl)*Xh+cos(el)*Yh+sin(el)*sin(rl)*Zh;
Za=cos(el)*sin(rl)*Xh-sin(el)*Yh+cos(el)*cos(rl)*Zh;
根据相对坐标分别计算雷达相对于侦察设备的方位、俯仰和计算侦察设备相对雷达的方位、俯仰;
根据通用函数表示出相对方位俯仰
AZ=f1(Xa,Ya);EL=f2(Ya,Za,AZ);
theta=f1(X,Y)={theta1=acrtg(|X|/|Y|),当X>0,Y>0时;theta2=π-theta1,当X>0,Y<0时;theta3=theta1-π,当X<0,Y<0时;theta4=-theta1,当X<0,Y>0时};
式中,0<theta1<=π/2,0<|theta|<=π,方位角向右为正,向左为负。
Beta=f2(Y,Z,theta)={beta0=arctg(|Z*cos(theta)|/|Y|),当Z>0;-beta0,当Z<0时};
式中,0<beta0<=π/2,0<|beta|<=π/2,俯仰角向上为正,向下为负。
把上一步计算的方位、俯仰转换到相对其天线主瓣的方位俯仰;
目标相对装备天线的方位=天线相对载体的方位+上一步计算的相对方位;
目标相对装备天线的俯仰=天线相对载体的俯仰+上一步计算的相对俯仰;
计算位于雷达方位、俯仰上的侦察设备天线增益和位于侦察设备方位、俯仰上的雷达天线增益;
计算天线增益方式有多种:
其一:数据形式的天线方向图,提供[方位(0~360)、俯仰(-90~90)、增益(361*181)]形式的天线方向图数据,然后根据上一步计算的相对方位、俯仰进行插值计算天线增益;
其二:采用高斯方向图函数,计算增益因子,公式如下:
F(az,el)=exp(-((az/thetah3dB)^2+(el/thetav3dB)^2))+Fs
其中,az为上一步计算的相对天线主瓣的方位,el为上一步计算的相对天线主瓣的俯仰,thetah3dB为天线主瓣波束水平方向的3dB波束宽度,thetav3dB为天线主瓣波束垂直方向的3dB波束宽度,Fs为天线平均旁瓣电平。
其三:采用单向余弦方向图函数,计算增益因子,公式如下:
F(az,el)=cos(π*(az)/(2*thetah3dB))*cos(π*(el)/(2*thetav3dB))
其中,参数定义同高斯方向图函数。
其四:采用单向辛克形方向图函数,计算增益因子,公式如下:
F(az,el)=(sin(2*π*az/thetah3dB)*sin(2*π*el/thetav3dB))/((*π*az/thetah3dB)*(2*π*el/thetav3dB))。
其参数定义同高斯方向图函数。
根据侦察设备和雷达之间的距离、其他电磁环境参数和选择合适的传播模型来计算传播损耗loss,上步计算的天线增益gain1(雷达),gain2(侦察)、雷达的发射功率power以及,计算侦察设备的接收功率recvPower;
recvPower=power+gain1+gain2–loss;
判断侦察设备的接收功率是否超过其接收灵敏度,如果超过则能量截获成功,否则截获失败。
所述步骤S3中空域截获关系计算根据作战场景中目标的位置、电子侦察设备的位置,目标雷达扫描波束范围,电子侦察设备波束扫描范围,计算电子侦察设备空域截获关系,通过计算相对天线方位俯仰算法,再基于天线波束宽度以及天线波束当前方位、俯仰,计算目标相对载体的方位、俯仰是否在波束扫描范围内,在范围内,则截获成功,否则截获失败。
所述步骤S3中频率截获关系计算根据电子侦察设备工作模式,参考电子侦察设备模型参数,结合雷达模型的数据,计算电子侦察设备空域截获关系,计算出侦察设备实时的频率,然后判断雷达频率是否在侦察设备实时扫描频率范围内。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (6)

1.基于参数化建模数据的电子侦察装备输出数据仿真方法,其特征在于,所述方法包括如下步骤:
S1:读取电子侦察设备建模数据和雷达建模数据后对电子侦察波束扫描仿真,实时计算电子侦察设备的波束指向;
S2:随后对电子侦察设备的工作频率进行仿真,通过侦察设备频率范围工作带宽、扫描时间、跟踪时间、重点频率集和实时仿真时间参数,实时计算当前电子侦察设备的工作频率参数;
S3:在对电子侦察截获关系进行仿真,通过能量截获计算,空域截获关系计算,频率截获关系计算获得侦察设备实时的接收灵敏度、雷达扫描波束、频率,然后判断雷达接收灵敏度、雷达扫描波束、频率是否在侦察设备实时扫描频率范围内;
S4:通过导入误差计算量进行全脉冲数据仿真,根据截获的目标雷达参数,对目标雷达的全脉冲波形参数进行仿真,包括全脉冲脉宽、脉冲间隔、脉冲数量、脉幅、脉冲串间隔、脉内调制类型进行仿真;
S5:导入异常介入参数后重新进行全脉冲数据仿真,将电子侦察设备建模数据中读取电子侦察设备异常频率、异常脉宽、异常重频、异常方位、异常时间参数以及异常类型导入,对仿真数据进行比对。
2.根据权利要求1所述的基于参数化建模数据的电子侦察装备输出数据仿真方法,其特征在于,所述步骤S1中实时计算电子侦察设备的波束指向的方式如下:
A1)断扫描方式,如果是扇扫则进步骤b,如果是圆扫则进步骤c,如果是固定则进步骤d;
A2)判断扫描方向,如果是顺时针扫描,则进步骤e,否则进步骤f;
A3)判断扫描方向,如果是顺时针扫描,则进步骤g,否则进步骤h;
A4)天线波束扫描方位,天线波束扫描俯仰;
A5)根据扫描时间计算扫描过的角度,angle=speed*time,计算扫描宽度width_deg=e_angle-s_angle,计算扫描过多少个扫描宽度cnt=floor(angle/width_deg),若cnt对2取余为0,则天线波束扫描方位=s_angle+(angle-cnt*width_deg),否则天线波束扫描方位=e_angle-(angle-cnt*width_deg);
A6)根据扫描时间计算扫描过的角度,angle=speed*time,计算扫描宽度width_deg=s_angle-e_angle,计算扫描过多少个扫描宽度cnt=floor(angle/width_deg),如果cnt对2取余为0,则天线波束扫描方位=s_angle-(angle-cnt*width_deg),否则天线波束扫描方位=e_angle+(angle-cnt*width_deg);
A7)天线波束扫描方位=s_angle+angle,再对扫描方位校验,换算到0-360°范围内;
A8)天线波束扫描方位=s_angle-angle,再对扫描方位校验,换算到0-360°范围内;
A9)天线波束扫描俯仰=c_rl+el。
3.根据权利要求1所述的基于参数化建模数据的电子侦察装备输出数据仿真方法,其特征在于,所述步骤S2中实时计算当前电子侦察设备的工作频率参数的方式如下:
B1)判断输入的重点频率集中的频率是否超出扫描范围,超出的则抛弃;
B2)然后计算基于工作带宽和最大最小频率的扫描次数cnt=floor((maxFreq-minFreq)/bandWidth);
B3)判断fomd((maxFreq-minFreq),bandWidth)是否大于零,如果大于零则cnt加1;
B4)计算顺序扫描需要时间normalTime=cnt*scanTime,计算跟踪需要时间followTime=size(freqs)*followTime;
B5)计算当前仿真时间处于顺序扫描时间段还是跟踪扫描时间段,realTime=fmod(time,normalTime+followTime),如果realTime小于normalTime,则处于顺序扫描阶段,否则处于跟踪扫描阶段;
B6)在正常扫描阶段,计算实时中心频率midFreq=minFreq+bandWidth/2+floor(realTime/scanTime)*bandwidth;
B7)在跟踪扫描阶段,计算实时中心频率midFreq=freqs[floor((realTime-normalTime)/followTime)];
B8)实时频率最小值realMinFreq=midFreq-bandWidth/2,如果小于最小频率,则使用最小频率;
B9)实时频率最大值realMaxFreq=midFreq+bandwidth/2,如果大于最大频率,则使用最大频率。
4.根据权利要求1所述的基于参数化建模数据的电子侦察装备输出数据仿真方法,其特征在于,所述步骤S3中能量截获计算根据场景仿真中自身平台姿态信息和其他平台的姿态信息,参考雷达的模型数据、电子侦察设备的模型数据、传播模型,计算电子侦察设备的接收功率,并与接收灵敏度进行比较。
5.根据权利要求1所述的基于参数化建模数据的电子侦察装备输出数据仿真方法,其特征在于,所述步骤S3中空域截获关系计算根据作战场景中目标的位置、电子侦察设备的位置,目标雷达扫描波束范围,电子侦察设备波束扫描范围,计算电子侦察设备空域截获关系,通过计算相对天线方位俯仰算法,再基于天线波束宽度以及天线波束当前方位、俯仰,计算目标相对载体的方位、俯仰是否在波束扫描范围内,在范围内,则截获成功,否则截获失败。
6.根据权利要求1所述的基于参数化建模数据的电子侦察装备输出数据仿真方法,其特征在于,所述步骤S3中频率截获关系计算根据电子侦察设备工作模式,参考电子侦察设备模型参数,结合雷达模型的数据,计算电子侦察设备空域截获关系,计算出侦察设备实时的频率,然后判断雷达频率是否在侦察设备实时扫描频率范围内。
CN202010328443.9A 2020-04-23 2020-04-23 基于参数化建模数据的电子侦察装备输出数据仿真方法 Active CN111539106B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010328443.9A CN111539106B (zh) 2020-04-23 2020-04-23 基于参数化建模数据的电子侦察装备输出数据仿真方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010328443.9A CN111539106B (zh) 2020-04-23 2020-04-23 基于参数化建模数据的电子侦察装备输出数据仿真方法

Publications (2)

Publication Number Publication Date
CN111539106A true CN111539106A (zh) 2020-08-14
CN111539106B CN111539106B (zh) 2021-02-19

Family

ID=71975520

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010328443.9A Active CN111539106B (zh) 2020-04-23 2020-04-23 基于参数化建模数据的电子侦察装备输出数据仿真方法

Country Status (1)

Country Link
CN (1) CN111539106B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113065808A (zh) * 2021-05-06 2021-07-02 中国电子科技集团公司第二十九研究所 模拟电子侦察数据上报率指标的方法、设备及存储介质

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1601298A (zh) * 2003-09-26 2005-03-30 清华大学 一种机载雷达的模型化杂波多普勒参数估计方法
CN104021243A (zh) * 2014-05-22 2014-09-03 万建岗 有源相控阵雷达参数化建模及综合优化设计平台
CN104298812A (zh) * 2014-09-11 2015-01-21 西安电子科技大学 一种针对雷达天线伺服系统的快速仿真分析方法
WO2016022255A2 (en) * 2014-07-11 2016-02-11 Advanced Testing Technologies, Inc. Phase noise simulation model for pulse doppler radar target detection
US20170139038A1 (en) * 2015-11-13 2017-05-18 Rohde & Schwarz Gmbh & Co. Kg Method and signal generator for simulation of sea clutter
CN109829251A (zh) * 2019-03-12 2019-05-31 中国人民解放军空军预警学院 基于算子的无源雷达目标截获概率柔性仿真建模方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1601298A (zh) * 2003-09-26 2005-03-30 清华大学 一种机载雷达的模型化杂波多普勒参数估计方法
CN104021243A (zh) * 2014-05-22 2014-09-03 万建岗 有源相控阵雷达参数化建模及综合优化设计平台
WO2016022255A2 (en) * 2014-07-11 2016-02-11 Advanced Testing Technologies, Inc. Phase noise simulation model for pulse doppler radar target detection
CN104298812A (zh) * 2014-09-11 2015-01-21 西安电子科技大学 一种针对雷达天线伺服系统的快速仿真分析方法
US20170139038A1 (en) * 2015-11-13 2017-05-18 Rohde & Schwarz Gmbh & Co. Kg Method and signal generator for simulation of sea clutter
CN109829251A (zh) * 2019-03-12 2019-05-31 中国人民解放军空军预警学院 基于算子的无源雷达目标截获概率柔性仿真建模方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
JUNGMIN PARK等: "Parameter Estimation of Radar Noise Model for Terrain Referenced Navigation Using a New EM Initialization Method", 《IEEE TRANSACTIONS ON AEROSPACE AND ELECTRONIC SYSTEMS》 *
梁广德等: "相控阵雷达信号截获与识别的仿真分析", 《航天电子对抗》 *
陶瞾杨等: "基于数字化样机的雷达目标跟踪技术仿真与验证", 《现代雷达》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113065808A (zh) * 2021-05-06 2021-07-02 中国电子科技集团公司第二十九研究所 模拟电子侦察数据上报率指标的方法、设备及存储介质
CN113065808B (zh) * 2021-05-06 2022-04-12 中国电子科技集团公司第二十九研究所 模拟电子侦察数据上报率指标的方法、设备及存储介质

Also Published As

Publication number Publication date
CN111539106B (zh) 2021-02-19

Similar Documents

Publication Publication Date Title
CN103616671B (zh) 一种相控阵雷达数字仿真系统及其仿真方法
CN111487597B (zh) 一种基于时空频能数据的通用电子侦察截获仿真方法
CN111812599A (zh) 博弈条件下基于低截获性能的组网雷达最优波形设计方法
CN106707255A (zh) 一种相控阵雷达仿真系统及方法
CN111539106B (zh) 基于参数化建模数据的电子侦察装备输出数据仿真方法
CN113225128A (zh) 一种面向无线激光通信的单层相位屏仿真方法
CN114442051A (zh) 一种高逼真度弹载雷达回波模拟方法
JP5199935B2 (ja) 時空間チャネルシミュレータ
CN109085552A (zh) 一种基于试飞数据的杂波半实物化仿真测试方法及系统
CN117111117A (zh) 全空域球顶相控阵天线多目标跟踪测控方法
Bi et al. Optimization method of passive omnidirectional buoy array in on-call anti-submarine search based on improved NSGA-II
Fuchs et al. Model order estimation using a multi-layer perceptron for direction-of-arrival estimation in automotive radar sensors
CN109342829B (zh) 电磁辐射源运动特性等效模拟方法
CN113109797B (zh) 调频连续波凝视雷达的目标探测方法、装置和计算机设备
Tong Modeling and realization of real time electronic countermeasure simulation system based on SystemVue
US6950057B1 (en) Complex radar target simulator
CN106405515B (zh) 一种天基雷达杂波仿真方法
CN109326172B (zh) 一种基于模型架构的雷达训练模拟器
CN115580364A (zh) 一种智能超表面技术辅助无人机信道建模方法
CN114239232A (zh) 基于规划评估双反馈技术的复杂电磁环境方案生成方法
CN112688724A (zh) 基于mimo技术的集群目标协同应答方法及询问应答系统
CN113391275A (zh) 一种采用电子角反的大型复杂目标等效模拟方法和装置
CN118209934B (zh) 空基分布式阵列波束合成方法及装置、电子设备、介质
CN117870872B (zh) 基于星载多波束微波辐射计的亮温图像数据集生成方法
CN117914364B (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
PE01 Entry into force of the registration of the contract for pledge of patent right

Denomination of invention: Simulation method of electronic reconnaissance equipment output data based on parametric modeling data

Effective date of registration: 20210722

Granted publication date: 20210219

Pledgee: Bank of Chengdu science and technology branch of Limited by Share Ltd.

Pledgor: Chengdu zhongxiangtiandi Network Technology Co.,Ltd.

Registration number: Y2021510000161

PE01 Entry into force of the registration of the contract for pledge of patent right
PC01 Cancellation of the registration of the contract for pledge of patent right

Date of cancellation: 20220615

Granted publication date: 20210219

Pledgee: Bank of Chengdu science and technology branch of Limited by Share Ltd.

Pledgor: Chengdu zhongxiangtiandi Network Technology Co.,Ltd.

Registration number: Y2021510000161

PC01 Cancellation of the registration of the contract for pledge of patent right
PE01 Entry into force of the registration of the contract for pledge of patent right

Denomination of invention: Simulation method of electronic reconnaissance equipment output data based on parametric modeling data

Effective date of registration: 20220628

Granted publication date: 20210219

Pledgee: Bank of Chengdu science and technology branch of Limited by Share Ltd.

Pledgor: Chengdu zhongxiangtiandi Network Technology Co.,Ltd.

Registration number: Y2022980009091

PE01 Entry into force of the registration of the contract for pledge of patent right