（92941部隊(duì)94分隊(duì) 葫蘆島 125001）

1 引言

火炮初速是火控系統(tǒng)解算或射表計(jì)算的重要參數(shù)之一，直接影響火炮系統(tǒng)的命中精度［1］。而產(chǎn)生火炮初速誤差的因素很多，如裝藥量微小差異、身管磨損變化、藥溫變化、彈丸質(zhì)量差異等，這些因素都很難在工藝或使用中加以控制，因而造成了初速誤差。隨著雷達(dá)測量技術(shù)與信號處理技術(shù)的發(fā)展，一些自動化程度較高的新型中、大口徑艦炮和自行榴彈炮上都裝有測速雷達(dá)［2］，用于實(shí)時完成火炮初速的測量，并將測量結(jié)果預(yù)測模型轉(zhuǎn)換，再通過數(shù)據(jù)接口或網(wǎng)絡(luò)接口傳到火控系統(tǒng)進(jìn)行下一次射擊的命中解算。由于這種初速測量的滯后性，因而測量的結(jié)果只能作為后續(xù)射擊預(yù)測值使用，勢必造成一定的修正誤差。

目前，火炮所采用的初速預(yù)測模型的統(tǒng)計(jì)樣本量不盡相同，主要是因?yàn)椴煌吞柕幕鹋谠谶B續(xù)射擊時初速分布規(guī)律有所差異，即使同一型號初速分布規(guī)律也可能不同。本文重點(diǎn)以實(shí)際測量數(shù)據(jù)統(tǒng)計(jì)分析不同初速分布規(guī)律下，采用不同樣本量的初速預(yù)測模型的精度。

2 初速分布規(guī)律

把一門火炮的初速V0變化看成是一個隨機(jī)過程，那么可以把這種過程分成兩類：一是其變化過程有確定形式或者說有必然規(guī)律，也就是說可以用一個（或幾個）時間t確定的函數(shù)來表示，即V（t）；二是其變化沒有必然的變化規(guī)律，即其變化過程不能用一個（或幾個）時間t的確定函數(shù)來表示［3］。彈種、裝藥號、彈重級、裝藥溫度、批次、裝藥存貯時間、火炮身管磨損與燒蝕等這些隨機(jī)因素，造成初速散布的隨機(jī)性［4］。

火炮作為利用密集度對敵目標(biāo)進(jìn)行覆蓋打擊的武器，一般作戰(zhàn)要求以最大射速進(jìn)行長連射或短連射進(jìn)行效力射擊［5］。本文主要研究短時間內(nèi)、射擊條件無太大變化的效力射擊時的初速分布與預(yù)測問題；對于歷史的初速數(shù)據(jù)，由于射擊條件不同，可能帶來較大統(tǒng)計(jì)誤差，因此不作為目前效力射擊的初速預(yù)測。

效力射擊過程中，火炮、彈藥的技術(shù)狀態(tài)以及氣象參數(shù)等影響因素基本保持不變，只是身管、藥室溫度會有些影響。試驗(yàn)表明，效力射時的初速分布規(guī)律會遇到以下兩種情況：

一種是初速分布存在一定弱相關(guān)關(guān)系（見圖1）：

如果首發(fā)彈的初速與第二發(fā)彈初速存在明顯的弱相關(guān)關(guān)系，對于該型炮為何出現(xiàn)這種現(xiàn)象目前還沒有辦法解釋。如圖1初速分布可以用這種弱相關(guān)的數(shù)學(xué)模型可表示為

式中：v（i）為第i發(fā)初速值；系數(shù)A為修正系數(shù)；B為隨機(jī)誤差；i＝1，2，…，N。

圖1 某型火炮七組七連發(fā)射擊的各發(fā)炮彈初速分布圖

另一種是一組初速分布為正態(tài)隨機(jī)分布，數(shù)學(xué)模型可表示為

式中：v為初速值為組平均初速；Δv為隨機(jī)誤差（滿足正態(tài)分布）。

3 初速的預(yù)測

火炮系統(tǒng)對初速數(shù)據(jù)的要求［6］：

1）各種條件下，精確預(yù)測初速；

2）快速初速預(yù)測與快速修正。

在上述條件下，一方面，要求測速雷達(dá)精度要高，預(yù)測模型要精確；另一方面，要求能夠向火炮系統(tǒng)快速提供預(yù)測結(jié)果，但這兩個方面卻成為了矛盾。根據(jù)統(tǒng)計(jì)學(xué)理論，當(dāng)統(tǒng)計(jì)射擊發(fā)數(shù)N達(dá)不到參數(shù)估計(jì)的要求，會降低了預(yù)測精度［7］，而反過來，如果統(tǒng)計(jì)射擊發(fā)數(shù)N過大，不僅對提高預(yù)測精度不利，還要浪費(fèi)火炮系統(tǒng)修正時間，那么如何既保證預(yù)測精度又能滿足一定的修正實(shí)時性？無疑在測速雷達(dá)精度一定的條件下，合理的預(yù)測方法是解決這一矛盾的關(guān)鍵。

3．1 樣本量的確定方法［8］

對于某型火炮估計(jì)值＾σV0、＾μV0在不同的試驗(yàn)條件下是不同的，而σV0一般對于某型火炮及彈藥是確定的，因此，我們利用估計(jì)精度值來確定樣本的大小。即：

即子樣大小滿足：

對于一般火炮試驗(yàn)確定值為σV0≈2m／s～3m／s，測速雷達(dá)的精度（相對精度）為1‰，綜合考慮，γ＝50%較為合適。通過上式計(jì)算，統(tǒng)計(jì)射擊發(fā)數(shù)N≥4。也就是說，要想滿足統(tǒng)計(jì)的預(yù)測精度，樣本數(shù)要大于等4發(fā)，無疑系統(tǒng)將在4×1／r（r為火炮射速）后進(jìn)行初速修正。

3．2 預(yù)測模型

效力射擊時，為了減小標(biāo)準(zhǔn)化轉(zhuǎn)換誤差，測速雷達(dá)所測初速數(shù)據(jù)一般不經(jīng)標(biāo)準(zhǔn)化轉(zhuǎn)換而直接按照一定的預(yù)測方法將數(shù)據(jù)直接通過網(wǎng)絡(luò)送至火控進(jìn)行使用［9］。

初速散布滿足式（1）的預(yù)測模型：

對于A、B值的估算，由于首發(fā)彈與第二發(fā)彈初速波動較大，并且先驗(yàn)數(shù)據(jù)較少，對于同一類型的其它火炮也不一定適用，很難進(jìn)行準(zhǔn)確統(tǒng)計(jì)估算，所以采取下面簡單方法進(jìn)行預(yù)測：

Δw是測得ν1后預(yù)測第二發(fā)的修正量，采用先前試驗(yàn)數(shù)據(jù)求取首發(fā)與第二發(fā)初速的平均差值作為Δw。

第三發(fā)彈的預(yù)測：由于第三發(fā)彈的初速不同于前兩發(fā)彈的初速，為了減少ν1的權(quán)系數(shù)而不用ν1，ν2和的一半作為預(yù)測值，公式為

第四發(fā)以后的彈丸初速預(yù)測均采用了前三發(fā)的滑動平均，即：

前三發(fā)彈初速的預(yù)測是依賴試驗(yàn)經(jīng)驗(yàn)值，其預(yù)測效果好壞，對后面的預(yù)測結(jié)果并沒有多大影響。

初速散布滿足式（2）的預(yù)測模型：

采用滑動平均值法。分別采用鄰近前3、4、5發(fā)初速與求平均值進(jìn)行初速預(yù)測。

累計(jì)平均值法。將之前測速雷達(dá)所測初速數(shù)據(jù)進(jìn)行平均，作為下一發(fā)的初速預(yù)測值［10］。

4 結(jié)果計(jì)算

由于試驗(yàn)條件限制，我們采用了幾組初速數(shù)據(jù)進(jìn)行分析，其中射擊條件、氣象等因素的影響這里不做分析，只對預(yù)測效果進(jìn)行分析。對于式（1）需要事先得到初速初值與Δw，這在實(shí)際作戰(zhàn)使用中并無太大意義，屬于初速預(yù)測的特例，該方法應(yīng)該慎用［11］。因此，我們只對具有普遍意義的滑動平均值法進(jìn)行分析。由于預(yù)測值的滯后性，火控系統(tǒng)中解算使用的是當(dāng)前射擊之前的初速統(tǒng)計(jì)值，為了更加直觀觀察分析預(yù)測效果優(yōu)劣，我們使用當(dāng)前的初速實(shí)測值Vi＋1實(shí)與之前初速預(yù)測值ˉVi＋1預(yù)一次差D來表示。

即：

用D3、D4、D5、D平均分別表示當(dāng)前實(shí)測值與鄰近前3發(fā)、前4發(fā)、前5發(fā)、前i發(fā)預(yù)測值的一次差，結(jié)果如圖2所示。將實(shí)測值與幾種預(yù)測模型預(yù)測值一次差結(jié)果進(jìn)行分段統(tǒng)計(jì)，統(tǒng)計(jì)結(jié)果如表2所示。

表1 某型火炮連續(xù)效力射的實(shí)測第一組初速值

利用同樣的方法，對另外一組數(shù)據(jù)進(jìn)行了分析，這一組的射擊條件、彈藥型號與第一組有所不同。一次差的分布圖如圖3所示，一次差統(tǒng)計(jì)結(jié)果如表4所示。

圖2 第一組實(shí)測值與預(yù)測值一次差分布圖

表2 第一組初速實(shí)測值與預(yù)測值一次差統(tǒng)計(jì)結(jié)果

表3 某型火炮連續(xù)效力射的實(shí)測第二組初速值

圖3 第二組實(shí)測值與預(yù)測值一次差分布圖

表4 第二組初速實(shí)測值與預(yù)測值一次差統(tǒng)計(jì)結(jié)果

根據(jù)上述分析，可以得到以下結(jié)論：

1）初速分布為隨機(jī)散布，不同射擊條件下初速散布會有所不同，第一組比第二組散布??；

2）初速隨機(jī)散布越大，初速預(yù)測偏差越大。第一組預(yù)測精度比第二組高很多；

3）采用的3 發(fā)、4 發(fā)、5 發(fā)的初速滑動平均預(yù)測這三種模型，預(yù)測精度提高并不明顯，而N發(fā)初速累計(jì)平均預(yù)測值精度會明顯高些，顯然符合統(tǒng)計(jì)學(xué)規(guī)律，因此為了使系統(tǒng)快速得到預(yù)測值，N取當(dāng)前射擊的前3發(fā)作為統(tǒng)計(jì)值可滿足要求；

4）N發(fā)累計(jì)平均值預(yù)測方法雖然預(yù)測精度有所提高，但也只是在效力射后期預(yù)測精度的提高，從作戰(zhàn)使用上看，對高射速長連射較為有利。

5 結(jié)語

對于某種火炮初速的實(shí)時預(yù)測模型，如果有條件應(yīng)該在大量試驗(yàn)的基礎(chǔ)上掌握其初速分布規(guī)律的基礎(chǔ)上進(jìn)行建立，能提高火炮首發(fā)命中精度。而一般情況下火炮初速分布屬正態(tài)隨機(jī)散布，采用N發(fā)初速平均值的預(yù)測模型是合適的?？傊?，測速雷達(dá)實(shí)際作用是提高了當(dāng)時射擊條件下火控系統(tǒng)的初速修正精度，從而提高火炮系統(tǒng)的射擊精度。

［1］李立緯，劉志偉，陳亮．提高艦炮射擊精度的若干途徑［J］．四川兵工學(xué)報，2011，33（8）：48-50．

［2］胡江，黃景德，解維河．基于測速雷達(dá)的艦炮初速測量技術(shù)研究［J］．艦船電子工程，2011，31（6）：94-96．

［3］馬百雙，劉昌錦．炮口初速測量方法綜述［J］．四川兵工學(xué)報，2011，32（8）：48-50．

［4］朱玉仙，崔曉光．概率論與數(shù)理統(tǒng)計(jì)［M］．沈陽：東北大學(xué)出版社，1998．

［5］張喜發(fā)，盧興華．火炮燒蝕內(nèi)彈道學(xué)［M］．北京：國防工業(yè)出版社，1998．

［6］李笑娟，倪晉平．高射速火炮連發(fā)測量技術(shù)［J］．兵工自動化，30（6）：68-70．

［7］金達(dá)根．實(shí)驗(yàn)外彈道學(xué)［M］．兵器工業(yè)出版社，1991．

［8］郭錫福．遠(yuǎn)程火炮武器系統(tǒng)射擊精度分析［M］．北京：國防工業(yè)出版社，2004．

［9］邱志明．艦炮武器系統(tǒng)分析［M］．北京：兵器工業(yè)出版社，1999．

［10］丁毅，陳曉明，王國博，等．炮口初速測量校準(zhǔn)與修正方法［J］．四川兵工學(xué)報，2011，32（10）：67-69．

［11］盧金柱，魯小強(qiáng)．火炮初速灰色預(yù)測方法研究［J］．艦船電子工程，2010，30（11）：35-37．