CN103049660A - 一种逐发瞄准火炮的击发频率的计算方法 - Google Patents
一种逐发瞄准火炮的击发频率的计算方法 Download PDFInfo
- Publication number
- CN103049660A CN103049660A CN2012105665549A CN201210566554A CN103049660A CN 103049660 A CN103049660 A CN 103049660A CN 2012105665549 A CN2012105665549 A CN 2012105665549A CN 201210566554 A CN201210566554 A CN 201210566554A CN 103049660 A CN103049660 A CN 103049660A
- Authority
- CN
- China
- Prior art keywords
- alpha
- sigma
- cannon
- error
- delta
- 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
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种逐发瞄准火炮的击发频率的计算方法。该方法运用随机穿越理论,考虑逐发瞄准火炮射击过程中存在射击准备时间,根据火炮的射击诸元误差穿越射击门的随机穿越频率,以及射击诸元误差在射击门内和射击门外的概率,从而得到射击诸元误差的随机穿越周期概率密度函数;然后根据火炮的随机击发周期的密度函数,计算出火炮的平均击发周期和平均射频。本发明可以为逐发瞄准火炮的射击频率的设计和检测提供一种有效的计算方法。
Description
技术领域
本发明属于武器系统效能分析领域,用于对逐发瞄准火炮的射击频率这一战技指标在进行武器系统设计和验证时进行计算。
背景技术
传统的逐发瞄准火炮,如地炮、突击炮等,其射击方式是:每一次击发完成后,射击冲击载荷使得火炮的身管产生随机振动,火炮在高低角和方位角上的射击诸元误差的大小和方向都是随机的,需要等待振动的过渡过程基本结束,并且射击准备(后座、复进、装填、瞄准)工作完成,火炮才可以进行下一次击发。对现代的逐发瞄准自动火炮,特别是行进间射击的自动火炮,不论是车载还是舰载,虽都配置了稳定装置,也难使其射击过程中的随机振动达到需求的稳定程度。为解决这一问题,行进间射击的直瞄火炮均以其火控解算诸元为中心、而间瞄火炮则以其陀螺装置构建的复瞄器所保持的火控解算诸元为中心,设置一个射击门。如果射击门为椭圆,射击诸元误差对其随机穿越过程如附图1所示。图中,坐标原点为火控解算射击诸元;以a,b为长短半轴的椭圆为预置的射击门,记为Ω;Z(t)=[x(t),y(t)]T∈R2中的x(t),y(t)分别为火炮方位角、高低角对应的射击诸元误差;称t1,t3,t5…为Z(t)对于射击门Ω的正(由外向内)穿越点,t2,t4,t6…为Z(t)对于射击门Ω的负穿越点。如果对Z(t)采用连续的检测装置(如自整角机或正余弦变压器),再考虑到火炮击发后,还需要射击准备时间Tδ,为能尽快的自行生成击发指令,逐发瞄准体制下的火炮击发瞬时的选择准则应是:第一个正穿越点t1是首发击发瞬时;后一击发瞬时与前一击发瞬时所对应的两个正穿越点间的时间间隔是大于Tδ中的最小者。很显然,受到上述随机振动过程的影响,逐发瞄准火炮的射击频率是随机的。与之不同的是对于速射武器,如高炮、机枪等,为了保证较高的射频,其射击过程中不需要等待每一次冲击振动结束就进行连续射击,其射频往往采用固定值。
由于逐发瞄准的火炮射频的均值是火力反应速度的关键性指标,它已被纳入战技指标体系。因此必须对火炮射频的随机特性加以分析和研究,以获得最佳的火炮射频。当前,有关逐发瞄准的火炮射频是通过大量的试验数字仿真或者试凑的办法来获得,需要花费大量的时间,并且缺乏理论依据。目前,对武器系统的射频进行研究的为文献也有所报道,但是都没有考虑逐发瞄准火炮在实际射击过程中存在需要射击准备时间这一不可或缺的约束,因而其研究结果不符合武器系统射击过程的实际情况。本发明针对上述逐发瞄准火炮的射击特点,特别是考虑逐发瞄准火炮存在需要射击准备时间这一不可或缺的约束,在研究火炮的射击诸元误差对射击门的随机穿越特性的基础上,给出了逐发瞄准火炮的随机射频的分布函数及其均值的数学计算方法。该方法适用于所有逐发瞄准的自动击发武器。
对于由炮手直接操瞄的非自动火炮,它的瞄准镜的不灵敏区的边界等效于射击门;对坦克炮而言,它已有了射击门,但其主要任务是提高首发命中率,如果在后续的射击过程中要求它们 进行多发射击,为提高它们的射速,本方法同样适用。
发明内容
本发明的目的在于提供一种逐发瞄准火炮的击发频率的计算方法。
实现方案:一种逐发瞄准火炮的击发频率的计算方法,根据测量射击准备时间、火炮的射击诸元误差穿越射击门的随机穿越频率以及射击诸元误差在射击门内和射击门外的概率,计算射击诸元误差的随机穿越周期概率的分布,然后得出火炮的随机击发周期的密度函数,最后计算出火炮的平均击发周期和平均射频。具体步骤如下:
1)根据战技指标要求或实际测量获得射击准备时间Tδ;
2)获得射击诸元误差穿越射击门的随机穿越平均频率λ;
3)获得射击诸元误差在射击门内概率α0和射击门外的概率α1;
4)根据前三步获得的参数带入公式
式中
得到射击诸元误差的n个随机穿越周期概率密度函数fch(t,n);
5)将fch(t,n)带入公式
步骤2中所述的λ可以根据火炮高低角和方位角的射击诸元误差的均方差和误差的导数的均方差,通过理论计算求的;也可以根据实际数据的检测直接获得。针对不同的射击门,理论计算方法不同,当射击门是椭圆时,λ的计算公式是:
步骤3中所述的α0和α1可以通过实测的射击诸元误差的均方差在射击门内外的占比统计获得或通过理论计算获得。当射击门为椭圆时,通过公式
式中σx,σy为射击诸元误差的均方差。
本发明所述的计算方法符合逐发瞄准火炮的实际射击规则,可以为逐发瞄准自动火炮的平均射频的论证、设计和检测提供一种有效的计算方法。为该类武器的设计方和验证方对平均射频这一新的战技指标提供了严谨的数学分析工具。
附图说明
图1为射击诸元误差穿越椭圆射击门示意图。
图2为射击诸元误差穿越椭圆射击门展开示意图。
图3处理后的火炮射击诸元误差数据小比例尺图。
具体实施方式
具体计算步骤如下:
1)根据战技指标要求或实际测量获得射击准备时间Tδ;
2)获得射击诸元误差穿越射击门的随机穿越平均频率λ;
λ可以根据火炮高低角和方位角的射击诸元误差的均方差和误差的导数的均方差,通过理论计算求的;也可以根据实际数据的检测直接获得。针对不同的射击门,理论计算方法不同,例如射击门是椭圆的λ的计算公式是:
如果射击门是矩形或其他不规则形状,可以通过相关文献的方法进行计算。
3)获得射击诸元误差在射击门内概率α0和射击门外的概率α1;
α0和α1可以通过实测的射击诸元误差的均方差在射击门内外的占比统计获得或通过理论计算获得。当射击门为椭圆时,通过公式
式中σx,σy为射击诸元误差的均方差,
式中
为虚自变量第一类Bessel函数。其收敛非常快,少许二、三项就有很高的精度,是I0(z)的常用近似表达。
4)根据前三步获得的参数带入公式
式中
得到射击诸元误差的n个随机穿越周期概率密度函数fch(t,n);
fch(t,n)证明过程如下:
若随机穿越频率λ为时不变常数时,单个滞留时间Tin=t服从指数分布为
单个待机时间Tout=t服从指数分布为
单个随机周期Tch=t服从卷积指数分布为
由于相互独立的n个连续随机变量的和的密度函数等于它们的密度函数的n-1次卷积,记*为卷积符号,考虑到滞留时间与待机时间都是相互独立的,故有:
式中
式中
5)将fch(t,n)带入公式
得到火炮的随机击发周期的密度函数;
将图1中的所有穿越点依时间轴展开,如图2所示,图2给出了从t>0点出发的三个可能的火炮射击诸元误差:它们分别表示在射击准备区间(0,Tδ]内出现n=0,1,2,3个正穿越点的四种情形,实际上,它们会依概率密度函数fch(t,n)给出的概率出现更多的正穿越点。依前述的逐发瞄准射击准则,如果击发周期发生于第n个随机周期的终端,那么,这n个随机周期之和,即击发周期Tj,它的分布密度为fch(t,n),Tj=t≥Tδ。又,上述事件等价于区间t∈(0,Tδ]内,必有n-1个随机周期,它发生的概率为fch(t,n-1)。显然,如果击发周期Tj=t≥Tδ发生在第n个随机周期的终端,那么,它出现的概率密度应为fch(Tδ,n-1)·f(t,n),t≥Tδ。根据全概率公式,对fch(Tδ,n-1)·f(t,n),t中n≥0的所有整数求和,可得火炮的随机击发周期的密度函数。
下面以2个具体的示例说明本发明的实施方式。
实施例1:
本例应用于逐发瞄准火炮射频的论证和设计。
(1)设椭圆射击门半径ρ=limil,以不相关的两个正态白噪声分别通过传递函数为
的成型滤波器,以其输出的不相关的各态历经函数x(t),y(t)作为在两个方向上的射击诸元误差的模拟量(设单位为mil),显然有
(2)计算出射击诸元误差对射击门的随机穿越频率,根据公式
λ=0.2067[次/秒]
(3)再计算出射击诸元误差在射击门内和射击门外的概率,根据公式
α0=0.1353,α1=0.8647
(4)若取射击准备时间Tδ=0.8s,再将α0,α1,λ带入公式
可计算出平均击发周期和平均射频:
实施例2
本例应用于武器系统的检验。
现有以40ms采样周期的实测火炮射击诸元误差数据三组,分别整合成具有2998个采样点、历时119.88秒的方位角x(k)、高低角y(k)曲线,如图3所示。
若取射击准备时间Tδ=1.9s,椭圆形射击门:短轴a=lmil,长轴b=1.768mil,现分析其平均射频γj。统计出相应的
λ=1.0331次/秒
α0=0.0832,α1=0.9168
将Tδ,α0,α1,λ带入公式
可计算出平均射频:
虽然结合了具体示例描述了本发明的具体实施方式,但是对于本领域技术人员,在不脱离本发明原理的前提下,还可以做出若干变形和改进,这些也应该视为本发明的保护范围。
Claims (5)
1.一种逐发瞄准火炮的击发频率的计算方法,其特征在于根据测量射击准备时间、火炮的射击诸元误差穿越射击门的随机穿越频率以及射击诸元误差在射击门内和射击门外的概率,计算射击诸元误差的随机穿越周期概率的分布,然后得出火炮的随机击发周期的密度函数,最后计算出火炮的平均击发周期和平均射频,具体步骤如下:
1)测量射击准备时间Tδ;
2) 获得射击诸元误差穿越射击门的随机穿越频率λ;
3)获得射击诸元误差在射击门内概率α0和射击门外的概率α1;
4)根据前三步获得的参数带入公式 ,式中
2.根据权利要求1所述的逐发瞄准火炮的击发频率的计算方法,其特征在于步骤2中所述的随机穿越频率λ根据火炮高低角和方位角的射击诸元误差的均方差和误差的导数的均方差理论计算求的,或根据实际数据的检测直接获得。
4.根据权利要求1所述的逐发瞄准火炮的击发频率的计算方法,其特征在于步骤2中所述的α0和α1通过实测的射击诸元误差的均方差在射击门内外的占比统计获得或通过理论计算获得。
5.根据权利要求1或4所述的逐发瞄准火炮的击发频率的计算方法,其特征在于当射击门为椭圆时,α0和α1通过公式得到:
式中σx,σy为射击诸元误差的均方差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210566554.9A CN103049660B (zh) | 2012-12-21 | 2012-12-21 | 一种逐发瞄准火炮的击发频率的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210566554.9A CN103049660B (zh) | 2012-12-21 | 2012-12-21 | 一种逐发瞄准火炮的击发频率的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103049660A true CN103049660A (zh) | 2013-04-17 |
CN103049660B CN103049660B (zh) | 2015-07-22 |
Family
ID=48062294
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210566554.9A Expired - Fee Related CN103049660B (zh) | 2012-12-21 | 2012-12-21 | 一种逐发瞄准火炮的击发频率的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103049660B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103134382A (zh) * | 2013-02-02 | 2013-06-05 | 中北大学 | 一种基于炮口脉冲噪声的火炮射速测试方法 |
CN106815378A (zh) * | 2015-11-27 | 2017-06-09 | 南京理工大学 | 连续波强激光武器动态毁伤概率的非毁检测方法 |
CN107908926A (zh) * | 2017-11-29 | 2018-04-13 | 中国人民解放军63850部队 | 一种射弹散布具有相关性的高炮毁伤概率确定方法 |
-
2012
- 2012-12-21 CN CN201210566554.9A patent/CN103049660B/zh not_active Expired - Fee Related
Non-Patent Citations (3)
Title |
---|
李文彬: "超高射频火炮点火控制装置设计及内弹道过程仿真", 《中国优秀硕士学位论文全文数据库》 * |
王军 等: "具有射击门体制的武器系统射击延时分析", 《火力与指挥控制》 * |
陶德进 等: "基于共有分量分解的速射火炮毁伤概率计算模型", 《兵工学报》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103134382A (zh) * | 2013-02-02 | 2013-06-05 | 中北大学 | 一种基于炮口脉冲噪声的火炮射速测试方法 |
CN106815378A (zh) * | 2015-11-27 | 2017-06-09 | 南京理工大学 | 连续波强激光武器动态毁伤概率的非毁检测方法 |
CN106815378B (zh) * | 2015-11-27 | 2020-04-28 | 南京理工大学 | 连续波强激光武器动态毁伤概率的非毁检测方法 |
CN107908926A (zh) * | 2017-11-29 | 2018-04-13 | 中国人民解放军63850部队 | 一种射弹散布具有相关性的高炮毁伤概率确定方法 |
CN107908926B (zh) * | 2017-11-29 | 2021-09-03 | 中国人民解放军63850部队 | 一种射弹散布具有相关性的高炮毁伤概率确定方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103049660B (zh) | 2015-07-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP2623921B1 (en) | Low-altitude low-speed small target intercepting method | |
KR101119882B1 (ko) | 사격제원 성능 분석 장치 및 그 동작 방법 | |
CN103049660B (zh) | 一种逐发瞄准火炮的击发频率的计算方法 | |
CN107908926A (zh) | 一种射弹散布具有相关性的高炮毁伤概率确定方法 | |
EP1604167B1 (en) | Method of making a projectile in a trajectory act at a desired point at a calculated point of time | |
CN206594195U (zh) | 远程枪弹测速鉴定装置 | |
CN105512405A (zh) | 一种mpc引射喷管直径的优化设计方法 | |
RU2290594C1 (ru) | Способ стрельбы боевой машины по скоростной цели (варианты) и система для его реализации | |
Donzier et al. | Small arm fire acoustic detection and localization systems: gunfire detection system | |
CN103186714B (zh) | 一种基于羽流场的固体推进剂二次烟雾特性虚拟评价方法 | |
KR101380790B1 (ko) | 3차원 교전 분석 장치 및 그 방법 | |
CN112629325A (zh) | 智能化城市低空反无人机系统 | |
KR102534405B1 (ko) | 공간정보기반 공역과 곡사화기 탄도 충돌 위협 분석 장치 및 방법 | |
RU2138757C1 (ru) | Способ стрельбы боевой машины по высокоскоростной цели и система для его осуществления | |
EP2120002B1 (de) | Verfahren zur Klassifikation von anfliegenden RAM-Geschossen | |
Tan et al. | Comparative study on the combat effectiveness of shipborne weapon against air based on ADC method | |
CN116793625A (zh) | 一种冲击波与破片复合冲击试验装置及方法 | |
RU2280836C1 (ru) | Способ защиты летательных аппаратов от управляемых ракет и система для его реализации | |
Jia-Lin et al. | Research on Requirements Analysis Method of Space Information Support Equipment System for Anti-missile Combat | |
CN117291008A (zh) | 一种试验剖面的制定和应力施加方法 | |
Massey et al. | A System-of-System Design of a Guided Projectile Mortar Defense System | |
Brown | Electric energy weapon implementation on an airborne platform | |
Kui et al. | General Design of Research-oriented Rehearsal System under Data Rules | |
Hongwei et al. | Measurement Method of the Air Miss Distance and ExplosionDistance of Proximity Fuse | |
CN118565273A (zh) | 一种液气相变发射装置、超高速侵彻试验装置及其测试方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
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: 20150722 Termination date: 20161221 |