CN108662955A - 一种基于光子探测的激光引信回波仿真方法 - Google Patents
一种基于光子探测的激光引信回波仿真方法 Download PDFInfo
- Publication number
- CN108662955A CN108662955A CN201810581202.8A CN201810581202A CN108662955A CN 108662955 A CN108662955 A CN 108662955A CN 201810581202 A CN201810581202 A CN 201810581202A CN 108662955 A CN108662955 A CN 108662955A
- Authority
- CN
- China
- Prior art keywords
- photon
- laser
- target
- intersection point
- fuze
- 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
Links
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F42—AMMUNITION; BLASTING
- F42C—AMMUNITION FUZES; ARMING OR SAFETY MEANS THEREFOR
- F42C13/00—Proximity fuzes; Fuzes for remote detonation
- F42C13/02—Proximity fuzes; Fuzes for remote detonation operated by intensity of light or similar radiation
-
- 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 Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Optical Radar Systems And Details Thereof (AREA)
Abstract
本发明公开了一种基于光子探测的激光引信回波仿真方法,首先建立激光引信坐标和目标模型,其次初始化发射激光并确定光子发射时的位置和方向,依次发射光子;然后判断光子能否照射到目标,计算被目标反射后的光子移动方向,再判断被目标反射后的光子能否被激光引信接收,记录光子被接收的时刻和能量;最后判断全部光子仿真是否结束,如果是统计每个时刻回波光子的总能量,就得到了激光引信的目标回波。本发明的仿真方法能够解决激光引信回波仿真精度低的问题。
Description
技术领域
本发明属于数值仿真技术领域,具体涉及一种激光引信回波仿真的方法。
背景技术
激光引信回波仿真是评估激光引信探测性能、测试目标散射特性重要且有效的手段。现有方案均是以激光束作为探测单元计算被探测目标回波,存在以下不足:(1)未考虑激光束截面光场强度分布的不均匀性;(2)未考虑激光引信接收窗口对回波角度的要求。上述不足导致基于激光束探测的激光引信回波仿真存在较大的误差。
发明内容
有鉴于此,本发明提供了一种基于光子探测的激光引信回波仿真方法,能够解决激光引信回波仿真精度低的问题。
一种基于光子探测的激光引信回波仿真方法,该方法的实现步骤如下:
步骤一、建立激光引信坐标系:以激光引信发射窗口为坐标原点,发射激光轴向为z轴,接收窗口相对发射窗口所在方向为x轴,按左手法则建立坐标系;
步骤二、建立目标模型:采用三角形网格将目标表面划分为多个面元,提取每个三角形面元的顶点位置、法线向量、反射率和双向反射分布函数;
步骤三、初始化发射激光:将激光引信发射的激光信号离散化,并将发射激光功率转换为发射的光子数量,得到光子的初始能量;
步骤四:确定光子发射时的位置和方向,依次发射光子;
步骤五、判断光子能否照射到目标:依次判断光子移动轨迹所在直线与目标每个面元是否有交点;若与所有面元均没有交点,则光子照射不到目标;若只与一个面元有交点,则光子照射到目标,且唯一的交点即为光子在目标上的照射位置;若与多个面元有交点,则光子照射到目标,距离光子发射位置最近的交点为光子在目标上的照射位置;
如果照射到目标进入下一步,如果未照射到目标则进入步骤九;
步骤六、计算被目标反射后的光子移动方向:根据目标双向反射分布函数计算得到目标散射方向概率分布,利用散射方向概率分布抽样确定反射后的光子移动方向;
步骤七、判断被目标反射后的光子能否被激光引信接收:若被目标反射后光子移动方向朝向激光引信一侧,计算光子到达激光引信接收窗口平面的位置,若光子到达位置在接收窗口内,且入射角度符合接收视场角,则光子被激光引信接收,成为回波光子;
如果是则进入下一步,如果否则进入步骤九;
步骤八:记录光子被接收的时刻和能量;
步骤九:判断全部光子仿真是否结束,如果是统计每个时刻回波光子的总能量,就得到了激光引信的目标回波;如果否则返回步骤四。
进一步地,所述步骤四中所述光子发射位置为:
其中ω0为激光束腰半径,ξ1,ξ2为标准正态分布随机数;
所述光子发射方向为:
其中θ0=|(θ'/2)·ξ3|为光子发射方向天顶角,θ'为激光光束发散角,ξ3为标准正态分布随机数,为光子发射方向方位角,ξ4为[0,1]区间上的均匀分布随机数。
进一步地,所述步骤五中判断光子移动轨迹所在直线与目标面元是否有交点的方法是:首先计算光子移动轨迹所在直线与面元所在平面之间的交点,再判断交点是否在面元内部;若交点在面元内部,则光子移动轨迹所在直线与面元有交点,否则没有交点。
进一步地,所述步骤八中,光子被接收时的能量为:
E'=ηeηrηtE
其中E为光子发射时的初始能量,ηe为激光引信发射光学系统的透过率,ηr为激光引信接收光学系统的透过率,ηt为目标的反射率。
有益效果:
本发明的方法采用光子作为探测单元,利用光子在空间中的分布来模拟激光束在空间中的光场强度分布,同时加入了激光回波角度判别,更加真实的模拟激光引信的回波接收过程,从而提高了激光引信目标回波的模拟精度。
附图说明
图1为基于光子探测的激光引信回波仿真流程图;
图2为激光引信坐标系示意图;
图3为目标模型图;
图4为激光引信发射激光的波形图;
图5为激光引信目标回波的波形图。
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
本发明提供了一种基于光子探测的激光引信回波仿真方法,以脉冲激光引信对距离5m处直升机的回波仿真为例,其仿真步骤如附图1所示,
步骤一、建立激光引信坐标系:以激光引信发射窗口为坐标原点,发射激光轴向为z轴,接收窗口相对发射窗口所在方向为x轴,按左手法则建立,如图2所示。
步骤二、建立目标模型:采用三角形网格将直升机表面划分为多个面元,如图3所示;提取每个三角形面元的顶点位置、法线向量、反射率和双向反射分布函数,格式如下:
第1列:面元编号;
第2~4列:顶点1坐标(tx1,ty1,tz1);
第5~7列:顶点2坐标(tx2,ty2,tz2);
第8~10列:顶点3坐标(tx3,ty3,tz3);
第11列:面元法线(VT1,VT2,VT3);
第12列:面元反射率ηt;
第13列:双向反射分布函数。
步骤三、初始化发射激光:激光引信发射激光为脉冲信号,峰值功率75W,如图4所示。设定激光引信峰值功率对应的光子数量为7500万个,每个光子的初始能量为:
其中PM为激光引信发射激光峰值功率,NM为峰值功率对应的光子数量。
以0.1ns间隔对发射激光脉冲信号进行离散化,激光引信在每个时刻发射的光子数量为:
其中P(i)为第i个时刻的激光功率。
步骤四、激光引信按照时间顺序依次发射每个时刻的光子,光子的初始位置为:
其中ω0为激光束腰半径,ξ1,ξ2为标准正态分布随机数。
光子发射方向为:
其中θ0=|(θ'/2)·ξ3|为光子发射方向天顶角,θ'为激光光束发散角,ξ3为标准正态分布随机数,为光子发射方向方位角,ξ4为[0,1]区间上的均匀分布随机数。
步骤五、判断光子能否照射到目标:依次判断光子移动轨迹所在直线与目标每个面元是否有交点,判断方法是:首先计算光子沿移动轨迹到达面元所在平面需要的距离:
然后计算光子移动轨迹所在直线与面元所在平面的交点:
最后判断交点是否在面元内部,若:
S1+S2+S3=S
则光子移动轨迹所在直线与面元有交点,否则没有交点;其中S为面元面积,S1、S2、S3为交点与面元两个顶点所围成的三角形区域面积。
若光子移动轨迹所在直线与目标的所有面元均没有交点,则光子照射不到目标;若只与一个面元有交点,则光子照射到目标,且唯一的交点即为光子在目标上的照射位置;若与多个面元有交点,则光子照射到目标,距离光子发射位置最近的交点为光子在目标上的照射位置。
如果照射到目标进入下一步,如果未照射到目标则进入步骤九;
步骤六、计算被目标反射后的光子移动方向:被朗伯体目标反射后的光子移动方向为:
其中为反射方位角,为[0,2π]区间上的均匀分布随机数,θ1为反射天顶角,根据双向反射分布函数得出的散射方向概率分布抽样确定,朗伯体目标的双向反射分布函数为:
P(θ1)=cosθ1
步骤七、判断被目标反射后的光子能否被激光引信接收:若被目标散射后的光子移动方向满足:
u′z<0
说明光子移动方向朝向激光引信一侧,计算光子到达激光引信接收窗口平面的位置:
其中(x2,y2,z2)为光子在目标上的照射位置,L2为光子从目标到达激光引信接收窗口平面的移动距离。
若光子到达位置在接收窗口内:
(x3-dtr)2+y3 2≤Rr 2
且入射角度符合接收视场角:
则光子被激光引信接收,成为回波光子;其中dtr为激光引信收发间距,Rr为接收窗口半径,θview为接收视场角。
如果是则进入下一步,如果否则进入步骤九;
步骤八、计算并保存光子被接收的时刻:
其中t0为光子发射时刻,L为光子在探测过程中的总移动距离,c为光速。
计算并保存光子被接收时的能量:
E'=ηeηrηtE
其中ηe为激光引信发射光学系统的透过率,ηr为激光引信接收光学系统的透过率,ηt为目标的反射率。
步骤九、判断全部光子仿真是否结束,如果是统计每个时刻回波光子的总能量:
Pr(t)=∑E'(t)
其中E'(t)为t时刻被接收光子的能量。
至此就得到了激光引信的目标回波,如图5所示,发射出的探测激光照射到直升机的两个距离不同的部位,因此得到两个回波峰值。如果否则返回步骤四。
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (4)
1.一种基于光子探测的激光引信回波仿真方法,其特征在于,该方法的实现步骤如下:
步骤一、建立激光引信坐标系:以激光引信发射窗口为坐标原点,发射激光轴向为z轴,接收窗口相对发射窗口所在方向为x轴,按左手法则建立坐标系;
步骤二、建立目标模型:采用三角形网格将目标表面划分为多个面元,提取每个三角形面元的顶点位置、法线向量、反射率和双向反射分布函数;
步骤三、初始化发射激光:将激光引信发射的激光信号离散化,并将发射激光功率转换为发射的光子数量,得到光子的初始能量;
步骤四:确定光子发射时的位置和方向,依次发射光子;
步骤五、判断光子能否照射到目标:依次判断光子移动轨迹所在直线与目标每个面元是否有交点;若与所有面元均没有交点,则光子照射不到目标;若只与一个面元有交点,则光子照射到目标,且唯一的交点即为光子在目标上的照射位置;若与多个面元有交点,则光子照射到目标,距离光子发射位置最近的交点为光子在目标上的照射位置;
如果照射到目标进入下一步,如果未照射到目标则进入步骤九;
步骤六、计算被目标反射后的光子移动方向:根据目标双向反射分布函数计算得到目标散射方向概率分布,利用散射方向概率分布抽样确定反射后的光子移动方向;
步骤七、判断被目标反射后的光子能否被激光引信接收:若被目标反射后光子移动方向朝向激光引信一侧,计算光子到达激光引信接收窗口平面的位置,若光子到达位置在接收窗口内,且入射角度符合接收视场角,则光子被激光引信接收,成为回波光子;
如果是则进入下一步,如果否则进入步骤九;
步骤八:记录光子被接收的时刻和能量;
步骤九:判断全部光子仿真是否结束,如果是统计每个时刻回波光子的总能量,就得到了激光引信的目标回波;如果否则返回步骤四。
2.如权利要求1所述的基于光子探测的激光引信回波仿真方法,其特征在于,所述步骤四中所述光子发射位置为:
其中ω0为激光束腰半径,ξ1,ξ2为标准正态分布随机数;
所述光子发射方向为:
其中θ0=|(θ'/2)·ξ3|为光子发射方向天顶角,θ'为激光光束发散角,ξ3为标准正态分布随机数,为光子发射方向方位角,ξ4为[0,1]区间上的均匀分布随机数。
3.如权利要求1所述的基于光子探测的激光引信回波仿真方法,其特征在于,所述步骤五中判断光子移动轨迹所在直线与目标面元是否有交点的方法是:首先计算光子移动轨迹所在直线与面元所在平面之间的交点,再判断交点是否在面元内部;若交点在面元内部,则光子移动轨迹所在直线与面元有交点,否则没有交点。
4.如权利要求1所述的基于光子探测的激光引信回波仿真方法,其特征在于,所述步骤八中,光子被接收时的能量为:
E'=ηeηrηtE
其中E为光子发射时的初始能量,ηe为激光引信发射光学系统的透过率,ηr为激光引信接收光学系统的透过率,ηt为目标的反射率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810581202.8A CN108662955B (zh) | 2018-06-07 | 2018-06-07 | 一种基于光子探测的激光引信回波仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810581202.8A CN108662955B (zh) | 2018-06-07 | 2018-06-07 | 一种基于光子探测的激光引信回波仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108662955A true CN108662955A (zh) | 2018-10-16 |
CN108662955B CN108662955B (zh) | 2019-08-16 |
Family
ID=63775619
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810581202.8A Expired - Fee Related CN108662955B (zh) | 2018-06-07 | 2018-06-07 | 一种基于光子探测的激光引信回波仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108662955B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110095035A (zh) * | 2019-04-22 | 2019-08-06 | 南京理工大学 | 导弹引信探测仿真方法 |
CN110986702A (zh) * | 2019-11-26 | 2020-04-10 | 天津津航技术物理研究所 | 一种用于激光引信回波采集的窗口自动遮挡装置 |
CN114372348A (zh) * | 2021-12-13 | 2022-04-19 | 北京理工大学 | 一种弹载线阵激光成像引信快速仿真方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5315506A (en) * | 1987-05-22 | 1994-05-24 | University Of Michigan | Correction for compton scattering by analysis of spatially dependent energy spectra employing regularization |
CN102594440A (zh) * | 2012-02-22 | 2012-07-18 | 大连大学 | 一种光子传输性能仿真方法 |
CN102706447A (zh) * | 2012-06-07 | 2012-10-03 | 中国航天科工集团第三研究院第八三五八研究所 | 一种激光探测系统会聚光斑质量测试方法 |
CN104197794A (zh) * | 2014-08-13 | 2014-12-10 | 上海无线电设备研究所 | 一种大视场目标探测激光引信收发光学系统 |
CN105486180A (zh) * | 2015-12-15 | 2016-04-13 | 西安电子科技大学 | 基于波束分解和局部照射的激光引信近场回波功率计算 |
-
2018
- 2018-06-07 CN CN201810581202.8A patent/CN108662955B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5315506A (en) * | 1987-05-22 | 1994-05-24 | University Of Michigan | Correction for compton scattering by analysis of spatially dependent energy spectra employing regularization |
CN102594440A (zh) * | 2012-02-22 | 2012-07-18 | 大连大学 | 一种光子传输性能仿真方法 |
CN102706447A (zh) * | 2012-06-07 | 2012-10-03 | 中国航天科工集团第三研究院第八三五八研究所 | 一种激光探测系统会聚光斑质量测试方法 |
CN104197794A (zh) * | 2014-08-13 | 2014-12-10 | 上海无线电设备研究所 | 一种大视场目标探测激光引信收发光学系统 |
CN105486180A (zh) * | 2015-12-15 | 2016-04-13 | 西安电子科技大学 | 基于波束分解和局部照射的激光引信近场回波功率计算 |
Non-Patent Citations (2)
Title |
---|
王凤杰,陈慧敏: "脉冲激光引信云雾回波特性仿真", 《光学 精密工程》 * |
王凤杰,陈慧敏等: "云雾后向散射激光回波特性研究", 《红外与激光工程》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110095035A (zh) * | 2019-04-22 | 2019-08-06 | 南京理工大学 | 导弹引信探测仿真方法 |
CN110095035B (zh) * | 2019-04-22 | 2021-06-29 | 南京理工大学 | 导弹引信探测仿真方法 |
CN110986702A (zh) * | 2019-11-26 | 2020-04-10 | 天津津航技术物理研究所 | 一种用于激光引信回波采集的窗口自动遮挡装置 |
CN110986702B (zh) * | 2019-11-26 | 2021-10-15 | 天津津航技术物理研究所 | 一种用于激光引信回波采集的窗口自动遮挡装置 |
CN114372348A (zh) * | 2021-12-13 | 2022-04-19 | 北京理工大学 | 一种弹载线阵激光成像引信快速仿真方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108662955B (zh) | 2019-08-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103278808B (zh) | 一种多线扫描式激光雷达装置 | |
CN105137415B (zh) | 激光测距仪接收视场标定及光轴平行度测量的装置和方法 | |
CN108445471B (zh) | 一种单光子激光雷达多探测器条件下的测距精度评估方法 | |
CN108662955B (zh) | 一种基于光子探测的激光引信回波仿真方法 | |
CN204044359U (zh) | 一种二维扫描式激光测距装置 | |
CN102243301B (zh) | 激光测距机检测装置 | |
CN110045341B (zh) | 雷达高度表低截获性能测试方法 | |
CN105486180B (zh) | 基于波束分解和局部照射的激光引信近场回波功率计算方法 | |
CN111060899A (zh) | 星地一体化激光雷达回波波形仿真方法及系统 | |
CN108519589A (zh) | 基于无源目标的星载激光测高仪足印定位方法 | |
CN110068803A (zh) | 一种雷达设备空中定标试验装置及方法 | |
CN114509744B (zh) | 一种激光雷达的测距检测率的评估方法、装置及设备 | |
CN107976686A (zh) | 一种多视场角海洋激光雷达及其视场角择优方法 | |
CN115935719A (zh) | 一种空中目标光子计数激光雷达探测仿真的方法和系统 | |
CN102323592A (zh) | 一种目标回波信号的归一化方法 | |
CN207937597U (zh) | 一种激光探测雷达 | |
CN107504862A (zh) | 一种全方位高精度激光定位方法 | |
CN102262230A (zh) | 卫星的测距方法和装置 | |
CN111880168A (zh) | 一种基于无源数字阵列雷达的目标定位方法 | |
CN107102310B (zh) | 一种多路激光雷达探测方法 | |
CN105445744A (zh) | 一种线激光目标探测系统和目标探测方法 | |
CN102200630A (zh) | 基于目标反射信号的光束瞄准系统中三角扫描偏差校准方法 | |
CN109211223A (zh) | 一种集成化接近测量敏感器及系统 | |
CN115033991A (zh) | 一种用于智能驾驶汽车仿真的毫米波雷达模拟方法及装置 | |
CN104865564A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190816 |
|
CF01 | Termination of patent right due to non-payment of annual fee |