Просмотр исходного текста страницы Теория вероятностей в артиллерии. Часть 1. Разрушение одиночной ДЗОТ

== Введение ==
Артиллерийское дело всегда было тесно связано со многими разделами математики. Движение снаряда в канале [[ствол]]а [[Артиллерийское орудие|орудия]], его [[Полёт снаряда|полёт в воздушной среде]], а также его действие на различного рода преграды являются сложными комплексами физических явлений. Для точного ведения огня необходимо рассмотреть действие очень многих причин на процесс стрельбы, выделить из них главные и отбросить несущественные, чтобы создать его физическую модель. Как только такая модель создана, физика и химия отходят в сторону и на первую роль выходит математика — нужно аналитически или численно решить уравнения и неравенства, требуемые для определения искомых величин — например максимальной дальности огня, дульной энергии [[снаряд]]а или износа канала ствола. По большому счёту все эти задачи являются детерминированными — зная конструкцию орудия и боеприпасов к нему, используемые в их изготовлении материалы, разработчики в лице конструкторов и учёных получат необходимые баллистические сведения о новой или модернизированной артиллерийской системе ещё до того, как она будет отстреляна на полигоне.

Но в погоне за точностью, дальнобойностью и пробивным действием не следует забывать и об экономности: от трёх первых факторов не будет никакого толка, если снаряды требуют экзотичных и редких материалов для своего изготовления, а орудия можно будет сделать только с привлечением самых что ни на есть высококвалифицированных рабочих в количестве пары штук в месяц. Здесь в игру вступают инженеры-технологи, которые по мере возможностей стараются объединить идеи учёных и конструкторов с суровой реальностью производственных возможностей. В итоге применительно к артиллерийской системе и боеприпасам для неё получается некоторый компромисс между изначальным проектом, и теми орудиями и снарядами, которые поставляются в войска. Причём далеко не всегда получается так, чтобы имелось изобилие боеприпасов и можно было бы не оглядываться на их расход на учениях и в бою. Здесь скорее имеет место обратная ситуация, когда этот расход нужно тщательно планировать, дабы сберечь не такой уж бесконечный запас снарядов и не беспредельный ресурс орудий. Не случайно именно это весьма настойчиво рекомендуется в наставлениях, руководствах службы и таблицах стрельбы. Поскольку на движение снаряда влияет целый сонм случайных изменений многих параметров, детерминистические методы расчёта дают только усреднённый результат, а на практике реальные значения величин имеют некоторый разброс около полученных из теории значений. Эти отклонения применительно к координатам точки падения снаряда называются рассеиванием. Из-за него даже в случае максимально точной наводки орудия поражение цели вовсе не является гарантированным результатом. Как следствие, возникает необходимость оценки количества снарядов, требуемых для выполнения учебной или боевой задачи. Тут на помощь артиллеристу вновь приходит математика, точнее её раздел, называемый теорией вероятности — учение о случайных событиях и закономерностях в их мире. В данной статье автор как раз и предлагает рассмотреть пример её использования на практическом примере, который является частью профессиональной подготовки артиллеристов.

== Постановка задачи ==
Рассмотрим ту же мысленную учебно-боевую задачу, которая обсуждалась в предыдущей статье «[[Что скрывается за внешней простотой]]» — [[152-мм гаубица образца 1943 года (Д-1)|152-мм гаубица обр. 1943 г. (Д-1)]] обстреливает дерево-земляную огневую точку (ДЗОТ) условного противника, удалённую от огневой позиции на 5000 метров. Будем считать, что наводка орудия с учётом угла места, деривации, метеорологических условий, состояния артиллерийской системы и боеприпасов уже точно выполнена — то есть центр эллипса рассеивания координат точек разрывов совмещён с целью. Для ведения огня назначена [[Осколочно-фугасный снаряд|осколочно-фугасная граната]] ОФ-530 с установкой взрывателя на замедленное (фугасное) действие, [[метательный заряд]] можно использовать любой. Необходимо оценить расход гранат, требуемых для разрушения ДЗОТ.

Очевидно, что для решения поставленной задачи нужна и некоторая информация о цели и разрушительных свойствах используемых боеприпасов. Пусть обитаемое пространство вражеской ДЗОТ в плане представляет собой прямоугольник со сторонами 3×4 м. Также известно, что при взрыве гранаты ОФ-530 в грунте средней плотности образуется воронка диаметром 3,5 м и глубиной 1,2 м. Поэтому для разрушения ДЗОТ в общем случае прямое попадание в него не требуется. Будем считать, что если точка разрыва удалена не далее чем на 1 метр по каждой из своих координат в горизонтальной плоскости от любой точки обитаемого пространства ДЗОТ, последняя гарантированно разрушается совместным действием механического воздействия, ударной волны и осколков. Иными словами, граната должна попасть тоже в прямоугольную область размером 5×6 м, центр которой совпадает с центром обитаемого пространства ДЗОТ.

== Определения и допущения ==
Строгие математические определения нужных нам для решения поставленной задачи понятий и величин можно найти в любом учебнике по теории вероятностей. На субъективный взгляд автора, в этом вопросе ничто не заменит замечательную книгу Елены Сергеевны Вентцель «Теория вероятностей», выдержавшую целый ряд переизданий, да ещё и написанную по материалам лекций по теории стрельбы в одном из военных высших учебных заведений. Здесь же поневоле для краткости изложения придётся воспользоваться некоторыми упрощениями и приближениями, которые, тем не менее, дадут достаточно точный результат.

=== Математические величины ===
Введём следующие обозначения:

:'''P<sub>1</sub>''' — вероятность разрушения ДЗОТ одним выстрелом. Её можно определить как отношение числа попавших в определённую выше окрестность ДЗОТ гранат к общему числу израсходованных боеприпасов '''N''', если условия каждого выстрела и полёта снаряда совершенно идентичны, а N стремится к бесконечности. В реальности, конечно, количество выстрелов ограничено и это отношение является так называемой относительной частотой требуемого нам события (попадания гранаты в заданную область и разрушения ДЗОТ как следствие). Но при N в диапазоне 10—20 отличие вероятности от относительной частоты уже не столь велико; а при N порядка 100 этой разницей уже можно пренебречь для любых практических расчётов.

:'''P<sub>i</sub>''' — вероятность того, что ДЗОТ окажется разрушенным после i = 2, 3, 4, … выстрелов; при условии, что все выпущенные ранее гранаты — 1-я, … , (i-1)-я — не попали в заданную окрестность цели. Это является следствием того факта, что для разрушения ДЗОТ достаточно только одного близкого попадания и после него продолжать стрельбу просто незачем. Соответственно математическое выражение для P<sub>i</sub> для этого условия согласно правилам сложения и умножения вероятностей независимых друг от друга событий будет следующим:

:: P<sub>i</sub> = P<sub>1</sub> + (1 - P<sub>1</sub>) × P<sub>1</sub> + (1 - P<sub>1</sub>)<sup>2</sup> × P<sub>1</sub> + … + (1 - P<sub>1</sub>)<sup>i - 1</sup> × P<sub>1</sub>

::Для других боевых задач, когда цель выводится из строя не гарантированно, а с некоторой вероятностью, условие и формула для P<sub>i</sub> могут быть гораздо более сложными.

: '''F(x, y)''' — двумерная плотность вероятности разрушения ДЗОТ одним выстрелом. Это действительная функция от координат x и y в горизонтальной плоскости, обладающая тем свойством, что при достаточно малых размерах прямоугольной области Δx × Δy в окрестности точки с координатами (x<sub>0</sub>, y<sub>0</sub>) вероятность P<sub>1</sub> может быть вычислена по формуле:

:: P<sub>1</sub> = F(x<sub>0</sub>, y<sub>0</sub>) × Δx × Δy

:: Критерием малости Δx и Δy можно считать небольшое изменение функции F (например менее 1%) при любом изменении её аргументов в диапазоне
:: (x<sub>0</sub> - ½Δx … x<sub>0</sub> + ½Δx, y<sub>0</sub> - ½Δy … y<sub>0</sub> + ½Δy)

::Более строгое определение плотности вероятности уже требует знания высшей математики, поскольку P<sub>1</sub> будет ничем иным, как двумерным интегралом от F(x, y) по требуемой нам окрестности. Однако для наших целей вполне пригодится и сильно упрощённый его вариант.

=== Алгоритм решения и используемые приближения ===
В артиллерийском деле принято считать функцию F(x, y) подчиняющейся так называемому нормальному закону распределения случайной величины — когда отклонение от её среднего значения вызывается совместным действием многих и независимых друг от друга факторов. Ещё это распределение называют гауссовым, в честь подробно исследовавшего его знаменитого немецкого математика Карла Гаусса. В подавляющем большинстве случаев это допущение хорошо соответствует действительности — небольшие случайные вариации в массе снаряда, сообщаемой ему начальной скорости (как по модулю, так и по направлению) и атмосферных условий друг с другом ни статистически, ни причинно между собой не связаны. Бывают однако и исключения — разброс точек попадания шрапнельных пуль или стреловидных поражающих элементов на горизонтальную плоскость при разрыве боеприпаса на нисходящей ветви своей траектории нормальному закону не подчиняется. Также закон распределения точек попадания снарядов на наклонную плоскую поверхность при небольших углах встречи отличен от нормального.

Математически функция F записывается следующим образом:

: F(x, y) = exp(-x² / 2S<sub>x</sub>²) × exp(-y² / 2S<sub>y</sub>²) / (2 × π × S<sub>x</sub> × S<sub>y</sub>) ,

: где начало координат совмещено с целью; ось X считается направленной вдоль проекции линии цели на горизонтальную плоскость уровня цели по высоте, а ось Y перпендикулярна оси X и лежит в той же плоскости. Величины S<sub>Д</sub> и S<sub>б</sub> называются стандартными или среднеквадратичными отклонениями по осям X и Y соответственно. Под exp(z) понимается показательная функция с возводимым в степень z основанием натурального логарифма e ≈ 2,7182… .

Однако в артиллерийском деле принято использовать в расчётах не стандартные, а так называемые срединные отклонения, которые составляют приблизительно 2/3 (более точно 0,674) от стандартных отклонений для гауссова распределения:

: [[Срединное отклонение]] по дальности — В<sub>Д</sub> = 0,674 × S<sub>x</sub>
: Боковое срединное отклонение (или срединное отклонение по фронту) — В<sub>б</sub> = 0,674 × S<sub>y</sub>

Для практического применения задание F(x, y) в аналитическом виде подходит плохо. Только ЭВМ способна быстро рассчитать нам вероятность по известным Δx и Δy, В<sub>Д</sub> и В<sub>б</sub> (последние приведены в таблицах стрельбы); даже с помощью инженерного калькулятора это займёт достаточно много времени с высокой вероятностью ошибки при выполнении операций под действием стресса. Однако нормальное распределение обладает рядом замечательных свойств, строго математический вывод которых превосходно изложен в учебнике Е. С. Вентцель. Нам будет полезно следующее из них: при большом числе выстрелов 50% X-координат попаданий попадут в отрезок от -В<sub>Д</sub> до В<sub>Д</sub>. Если увеличить этот отрезок в два раза, то внутри него очутится 82% всех попаданий, а если в три — то 96%. Оставшимися 4% на практике обычно пренебрегают, причём львиная доля из них уложится в ещё раз расширенный интервал от -4В<sub>Д</sub> до 4В<sub>Д</sub>. Абсолютно также дело обстоит и с аналогичной величиной по оси Y, только срединное отклонение здесь равно  В<sub>б</sub>. Такой характер распределения координат точек попаданий довольно удобен при расчётах на практике, именно этим (с определённой долей исторической традиции) и объясняется использование срединного отклонения в артиллерийском деле вместо среднеквадратичного, которое более привычно для учёных и статистиков. Для сравнения, при нормальном законе распределения случайной величины с нулевым средним в интервал от минус среднеквадратичного отклонения до плюс среднеквадратичного отклонения попадает около 69% всех выпавших её значений, вместо 50% для срединного отклонения.

В нашем случае отклонения по дальности и фронту не зависят друг от друга, что математически выражается в возможности представить функцию F(x, y) в виде произведения F’(x) × F”(y). Без особого ущерба для точности в каждом из интервалов с шагом в одно срединное отклонение значения функций F’(x) и F”(y) можно считать постоянными. Эти значения должны при их перемножении на длину интервала давать тот же результат, что и определённый интеграл от каждой из этих функций на этом интервале. Таким образом их плавные колоколообразные формы подменяются чем-то вроде ступенчатого пьедестала: 

[[Файл:Плотность_вероятности_1.png|thumb|center|800px|
Иллюстрация принятой в артиллерийском деле аппроксимации нормального распределения ступенчатой функцией. Стандартное (среднеквадратичное) отклонение равно 1, соответствующее ему срединное отклонение — 0,67. По оси абсцисс нижняя шкала чёрного цвета размечена в единицах, равных одному стандартному (среднеквадратичному) отклонению, верхняя шкала зелёного цвета размечена в единицах, равных одному срединному отклонению.]]
<br clear="all" />

Это существеннейшим образом упрощает счёт, позволяя обойтись только тривиальными арифметическими действиями. В самом деле, вероятность попадания в прямоугольник со сторонами 2×В<sub>Д</sub> и 2×В<sub>б</sub> равна 25% — поскольку приблизительно в половине случаев при большом числе выстрелов произойдёт попадание в нужный интервал по оси X и, как следствие независимости отклонений по дальности и фронту, только половина от этой половины количества разрывов ляжет требуемым образом ещё и по оси Y. В силу сделанного выше допущения шансы попасть в любые точки внутри рассматриваемого прямоугольника одинаковы, поэтому вероятность попасть в цель площадью S внутри него есть не что иное, как одна четвёртая от отношения S к его площади:

: P<sub>1</sub> = ¼ × S / (4 × В<sub>Д</sub> × В<sub>б</sub>)

Теперь у нас есть практически пригодный алгоритм решения поставленной задачи. Осталось лишь подставить в него S — площадь ДЗОТ и его близкой окрестности в плане, равную 5 × 6 = 30 м², и значения В<sub>Д</sub> и В<sub>б</sub> для различных метательных зарядов из таблиц стрельбы для 152-мм гаубицы обр. 1943 г. (Д-1).

== Результаты расчётов ==
{| cellpadding="2" style="background:#fff; border-collapse:collapse; border:1px solid #996; margin:1ex 0; text-align:center;  width:98%;"
|- style="background-color:#996; color:#fff;"
| colspan="5" style="border:1px solid #996"| '''Расчёт вероятности разрушения удалённого на 5000 м ДЗОТ с S = 30 м² для различных зарядов одним выстрелом'''
|-
| style="border:1px solid #996; width:auto; background-color:#ddb;" rowspan="2"| Заряд
| style="border:1px solid #996; width:20%; background-color:#ddb;"| '''В<sub>Д</sub>'''
| style="border:1px solid #996; width:20%; background-color:#ddb;"| '''В<sub>б</sub>'''
| style="border:1px solid #996; width:23%; background-color:#ddb;"| '''P<sub>1</sub>''' = ¼×'''S''' / ('''4'''×'''В<sub>Д</sub>'''×'''В<sub>б</sub>''')
| style="border:1px solid #996; width:23%; background-color:#ddb;"| '''P<sub>1</sub>''' теор.
|-
| style="border:1px solid #996; width:20%; background-color:#ddb;"| м
| style="border:1px solid #996; width:20%; background-color:#ddb;"| м
| style="border:1px solid #996; width:20%; background-color:#ddb;"| безразмерно
| style="border:1px solid #996; width:20%; background-color:#ddb;"| безразмерно
|-
| style="border:1px solid #996; text-align:left"| Полный
| style="border:1px solid #996"| 18
| style="border:1px solid #996"| 2,6
| style="border:1px solid #996"| 0,040
| style="border:1px solid #996"| 0,043
|-
| style="border:1px solid #996; text-align:left"| Первый
| style="border:1px solid #996"| 18
| style="border:1px solid #996"| 3,0
| style="border:1px solid #996"| 0,034
| style="border:1px solid #996"| 0,038
|-
| style="border:1px solid #996; text-align:left"| Второй
| style="border:1px solid #996"| 18
| style="border:1px solid #996"| 3,2
| style="border:1px solid #996"| 0,033
| style="border:1px solid #996"| 0,036
|-
| style="border:1px solid #996; text-align:left"| Третий
| style="border:1px solid #996"| 17
| style="border:1px solid #996"| 3,4
| style="border:1px solid #996"| 0,032
| style="border:1px solid #996"| 0,036
|-
| style="border:1px solid #996; text-align:left"| Четвёртый
| style="border:1px solid #996"| 24
| style="border:1px solid #996"| 3,4
| style="border:1px solid #996"| 0,023
| style="border:1px solid #996"| 0,025
|-
| style="border:1px solid #996; text-align:left"| Пятый
| style="border:1px solid #996"| 26
| style="border:1px solid #996"| 3,5
| style="border:1px solid #996"| 0,021
| style="border:1px solid #996"| 0,023
|-
| style="border:1px solid #996; text-align:left"| Шестой
| style="border:1px solid #996"| 28
| style="border:1px solid #996"| 3,7
| style="border:1px solid #996"| 0,018
| style="border:1px solid #996"| 0,020
|}

Приведённая в самом правом столбце величина '''P<sub>1</sub>'''&nbsp;теор. является рассчитанной путём интегрирования по точной формуле вероятностью разрушения ДЗОТ одним выстрелом. Как видно из данных таблицы, упрощённый алгоритм вносит некоторую погрешность, которая тем больше, чем выше кучность боя орудия — с возрастанием срединных отклонений В<sub>Д</sub> и В<sub>б</sub> приближённый результат стремится к точному. Это связано с тем, что в нужном диапазоне изменения аргументов x и y при больши́х В<sub>Д</sub> и В<sub>б</sub> функция F меняется мало и вполне может быть заменена неким постоянным значением. По дальности это условие заведомо выполняется (наименьшее В<sub>Д</sub>, равное 17 метрам, в разы превосходит половину продольного размера окрестности ДЗОТ), а вот по фронту — не совсем (наименьшее В<sub>б</sub>, равное 2,6 метра, лишь чуть больше, чем половина поперечного размера цели), см. следующий рисунок:

[[Файл:Эллипс рассеивания у Д-1.gif|thumb|center|800px|Эллипс рассеивания с полуосями В<sub>Д</sub> = 28 м и В<sub>б</sub> = 3,7 м у 152-мм гаубицы обр. 1943 г. (Д-1) для шестого заряда и дальности 5000 метров. В закрашенную красным цветом площадь приходится около 20% всех разрывов при ударной стрельбе. Бордовым обозначена окрестность ДЗОТ, попадание в которую ведёт к выполнению поставленной боевой задачи.]]
<br clear="all" />

Впрочем максимальное отличие в 4% между точными и приближёнными величинами на практике сколь нибудь значимой роли не играет, в чём можно убедиться, рассмотрев вероятности разрушения цели серией из 10, 20, 30 и 40 выстрелов:

{| cellpadding="2" style="background:#fff; border-collapse:collapse; border:1px solid #996; margin:1ex 0; text-align:center;  width:98%;"
|- style="background-color:#996; color:#fff;"
| colspan="5" style="border:1px solid #996"| '''Расчёт вероятности разрушения удалённого на 5000 м ДЗОТ с S = 30 м² для различных зарядов серией из N выстрелов'''
|-
| style="border:1px solid #996; width:auto; background-color:#ddb;"| Заряд
| style="border:1px solid #996; width:20%; background-color:#ddb;"| '''P<sub>10</sub>'''
| style="border:1px solid #996; width:20%; background-color:#ddb;"| '''P<sub>20</sub>'''
| style="border:1px solid #996; width:20%; background-color:#ddb;"| '''P<sub>30</sub>'''
| style="border:1px solid #996; width:20%; background-color:#ddb;"| '''P<sub>40</sub>'''
|-
| style="border:1px solid #996; text-align:left"| Полный
| style="border:1px solid #996"| 0,34 / 0,36
| style="border:1px solid #996"| 0,56 / 0,58
| style="border:1px solid #996"| 0,71 / 0,73
| style="border:1px solid #996"| 0,80 / 0,83
|-
| style="border:1px solid #996; text-align:left"| Первый
| style="border:1px solid #996"| 0,29 / 0,32
| style="border:1px solid #996"| 0,50 / 0,54
| style="border:1px solid #996"| 0,65 / 0,69
| style="border:1px solid #996"| 0,75 / 0,79
|-
| style="border:1px solid #996; text-align:left"| Второй
| style="border:1px solid #996"| 0,29 / 0,31
| style="border:1px solid #996"| 0,49 / 0,52
| style="border:1px solid #996"| 0,63 / 0,67
| style="border:1px solid #996"| 0,74 / 0,77
|-
| style="border:1px solid #996; text-align:left"| Третий
| style="border:1px solid #996"| 0,28 / 0,31
| style="border:1px solid #996"| 0,48 / 0,52
| style="border:1px solid #996"| 0,62 / 0,67
| style="border:1px solid #996"| 0,73 / 0,77
|-
| style="border:1px solid #996; text-align:left"| Четвёртый
| style="border:1px solid #996"| 0,21 / 0,22
| style="border:1px solid #996"| 0,37 / 0,40
| style="border:1px solid #996"| 0,50 / 0,53
| style="border:1px solid #996"| 0,61 / 0,64
|-
| style="border:1px solid #996; text-align:left"| Пятый
| style="border:1px solid #996"| 0,19 / 0,21
| style="border:1px solid #996"| 0,35 / 0,37
| style="border:1px solid #996"| 0,47 / 0,50
| style="border:1px solid #996"| 0,57 / 0,61
|-
| style="border:1px solid #996; text-align:left"| Шестой
| style="border:1px solid #996"| 0,17 / 0,18
| style="border:1px solid #996"| 0,30 / 0,33
| style="border:1px solid #996"| 0,42 / 0,45
| style="border:1px solid #996"| 0,52 / 0,55
|-
| style="border:1px solid #996; text-align:left; font-size:80%; padding-left:10px;" colspan="5"| Значение до наклонной черты соответствует вероятности, рассчитанной по приближённому алгоритму, после — точному.
|}

Стоит заметить, что точная формула даёт чуть больший шанс на успех, поскольку внутри интервала от -В<sub>Д</sub> до В<sub>Д</sub> попадания всё же распределены не равномерно, как предполагается в используемом приближении, а тяготеют к нулю. Но, как видно из следующей таблицы, этот выигрыш не является сколь-нибудь существенно значимым в пересчёте на количество снарядов, требуемое для поражения цели с заданной вероятностью, например 50%:

{| cellpadding="2" style="background:#fff; border-collapse:collapse; border:1px solid #996; margin:1ex 0; text-align:center;  width:98%;"
|- style="background-color:#996; color:#fff;"
| colspan="6" style="border:1px solid #996"| '''Расчёт количества выстрелов, нужного для разрушения удалённого на 5000 м ДЗОТ с S = 30 м² с заданной вероятностью'''
|-
| style="border:1px solid #996; width:auto; background-color:#ddb;"| Заряд
| style="border:1px solid #996; width:16%; background-color:#ddb;"| '''50%'''
| style="border:1px solid #996; width:16%; background-color:#ddb;"| '''75%'''
| style="border:1px solid #996; width:16%; background-color:#ddb;"| '''90%'''
| style="border:1px solid #996; width:16%; background-color:#ddb;"| '''99%'''
| style="border:1px solid #996; width:16%; background-color:#ddb;"| '''99,9%'''
|-
| style="border:1px solid #996; text-align:left"| Полный
| style="border:1px solid #996"| 17 / 16
| style="border:1px solid #996"| 34 / 32
| style="border:1px solid #996"| 57 / 53
| style="border:1px solid #996"| 113 / 105
| style="border:1px solid #996"| 170 / 158
|-
| style="border:1px solid #996; text-align:left"| Первый
| style="border:1px solid #996"| 21 / 18
| style="border:1px solid #996"| 41 / 36
| style="border:1px solid #996"| 67 / 60
| style="border:1px solid #996"| 134 / 119
| style="border:1px solid #996"| 200 / 179
|-
| style="border:1px solid #996; text-align:left"| Второй
| style="border:1px solid #996"| 21 / 19
| style="border:1px solid #996"| 42 / 38
| style="border:1px solid #996"| 69 / 63
| style="border:1px solid #996"| 138 / 126
| style="border:1px solid #996"| 206 / 189
|-
| style="border:1px solid #996; text-align:left"| Третий
| style="border:1px solid #996"| 22 / 19
| style="border:1px solid #996"| 43 / 38
| style="border:1px solid #996"| 71 / 63
| style="border:1px solid #996"| 142 / 126
| style="border:1px solid #996"| 213 / 189
|-
| style="border:1px solid #996; text-align:left"| Четвёртый
| style="border:1px solid #996"| 30 / 28
| style="border:1px solid #996"| 60 / 55
| style="border:1px solid #996"| 99 / 91
| style="border:1px solid #996"| 198 / 182
| style="border:1px solid #996"| 297 / 273
|-
| style="border:1px solid #996; text-align:left"| Пятый
| style="border:1px solid #996"| 33 / 30
| style="border:1px solid #996"| 66 / 60
| style="border:1px solid #996"| 109 / 99
| style="border:1px solid #996"| 217 / 198
| style="border:1px solid #996"| 326 / 297
|-
| style="border:1px solid #996; text-align:left"| Шестой
| style="border:1px solid #996"| 39 / 35
| style="border:1px solid #996"| 77 / 69
| style="border:1px solid #996"| 127 / 114
| style="border:1px solid #996"| 254 / 228
| style="border:1px solid #996"| 381 / 342
|-
| style="border:1px solid #996; text-align:left; font-size:80%; padding-left:10px;" colspan="5"| Значение до наклонной черты соответствует количеству выстрелов, рассчитанному по приближённому алгоритму, после — точному.
|}

Однако третья таблица неинформативна в практическом смысле: из неё трудно найти ответ на вопрос, а сколько же в среднем нужно снарядов для поражения цели и каков будет разброс значений возле этой средней величины? Естественно, что теория вероятностей знает, как вычислить эти характеристики. Но для большей наглядности и меньшего числа математических выкладок автор предлагает решить этот вопрос путём имитационного моделирования. Конечно, было бы неплохо взять настоящую гаубицу Д-1, выкатить её на полигон и по результатам достаточно многочисленной серии выстрелов получить данные из натурного эксперимента. Но для более-менее приемлемой точности нужно произвести около одной тысячи выстрелов, что является крайне дорогостоящим мероприятием во всех отношениях. К счастью, вычислительные мощности в наше время стали настолько дешёвыми, что интересующий нас результат вполне можно получить за пару минут. Для этого нужно только составить несложную программу, которая выдаёт пару псевдослучайных величин с описанным выше нормальным законом распределения до тех пор, пока они не удовлетворят критерию попадания в заданную область; при этом считается число таких «виртуальных выстрелов» до достижения «попадания». На языке PHP реализация этого алгоритма может иметь такой вид:

<source lang="php">
  define(SX, 17 / 0.674);      // Пересчёт срединного отклонения по дальности в среднеквадратичное
  define(SY, 3.4 / 0.674);     // То же самое для бокового срединного отклонения
  define(DX, 6 / 2);           // Вычисляем "полудлину" ДЗОТ и её окрестности
  define(DY, 5 / 2);           // Вычисляем "полуширину" ДЗОТ и её окрестности

// Инициализируем генератор псевдослучайных чисел случайным образом по машинному таймеру
  srand ((double) microtime() * 1000000);

  $n = 0;   // Число "выстрелов"
  do
  {
    $n++;
    $coordinates = Normal_distributed(SX, SY);                        // Получаем координаты "разрыва"
  } while (abs($coordinates[0]) > DX || abs($coordinates[1]) > DY);   // Если не попали, "стреляем" ещё раз

  echo "Hit after $n shots.";  // Выводим сообщение о количестве "выстрелов"
  exit;

// Функция, выдающая пару независимых нормально распределённых случайных величин с заданными среднеквадратичными
// отклонениями. Алгоритм взят из книги:
// Дж. Форсайт, М. Малькольм, К. Моулер "Машинные методы математических вычислений"

  function Normal_distributed($std1 = 1.0, $std2 = 1.0)
  {
    do
    {
      $v1 = 1 - 2.0 * rand(0, 999999999) / 1000000000;
      $v2 = 1 - 2.0 * rand(0, 999999999) / 1000000000;
      $s  = $v1 * $v1 + $v2 * $v2;
    } while ($s > 1);

    $coeff = sqrt(-2.0 * log($s) / $s);
    return array(0 => $std1 * $v1 * $coeff, $std2 * $v2 * $coeff);
  }
</source>

Путём несложной модификации этой программы можно получить много интересных результатов, из которых автору показались наиболее интересными следующие: распределение точек «разрывов» на горизонтальной плоскости и функция распределения числа выстрелов в серии до разрушения ДЗОТ. Графическая иллюстрация первого результата моделирования для случая, когда попадание было достигнуто 39-ым «выстрелом на третьем заряде», приведена ниже:

[[Файл:Рассеивание у Д-1 (численное моделирование).png|thumb|center|800px|Зелёным пунктиром отмечены удаления от точки прицеливания от одного до четырёх срединных отклонений по дальности по оси X и от одного до четырёх срединных боковых отклонений по оси Y. Площадь, попадание в которую приводит к разрушению ДЗОТ, выделена серым цветом с чёрной окантовкой.]]
<br clear="all" />

Вообще же в «обстрелах на третьем заряде нескольких ДЗОТ, удалённых на 5 км» наблюдались как случаи попадания в них первой гранатой, так и серии за полторы сотни выстрелов до достижения успеха. Но среднее число лежало где-то в окрестности 30. Для максимального его уточнения автором был проведён численный эксперимент по «разрушению» 100000 ДЗОТ с фиксацией количества целей, «приведённых к молчанию» серией из N выстрелов. Его результаты изображены на следующем графике (который есть не что иное, как обобщённая функция плотности вероятности для дискретной случайной величины N, умноженная на некий коэффициент):

[[Файл:Распределение по числу выстрелов в серии.png|thumb|center|800px|Распределение числа «разрушенных ДЗОТ» по количеству «выстрелов в серии до их разрушения».]]
<br clear="all" />

Для знающих математику здесь нет абсолютно ничего нового — чистой воды экспоненциальное распределение со средним, равным обратной вероятности попадания в ДЗОТ (28 выстрелов), и равным ей же среднеквадратичным отклонением. Эти величины можно уже с полным правом использовать в практических расчётах, но с одной очень важной оговоркой. Впрочем об этом чуть позже.

== Комментарии ==
Как мог бы заметить сам Капитан Очевидность, вероятность поразить цель зависит от номера метательного заряда и при заданном удалении и отсутствии препятствий на пути снаряда имеется определённая возможность повысить шансы на успех, выбрав полный заряд с наименьшей крутизной траектории. Это рекомендовалось ещё в другой замечательной книге «Артиллерия», изданной почти 60 лет назад, в 1953 году. Но при этом надо учитывать тот факт, что стрельба на полном заряде ведёт к быстрому износу орудия, а на малых зарядах при приблизительно вдвое меньшей вероятности попадания ресурс ствола возрастает в 4-5 (а иногда и больше) раз. Поэтому надо найти определённый компромисс между шансами на попадание (т. е. расходом боеприпасов) и сбережением артиллерийской системы. В нашем случае третий заряд будет неплохим выбором с этой точки зрения.

Стоит также обратить внимание на то обстоятельство, что третья таблица в пух и прах развеивает популярное поверье о том, что если шансы на успех в отдельно взятой и независимой от исхода других попытке составляют 4%, то двадцати пяти попыток будет достаточно для гарантированного выпадения благоприятного исхода. Как видно из данных стрельбы на полном заряде, поразить цель, имея даже не 25, а 35 выстрелов, удастся только в 75% случаев. Если нам нужно обстрелять, к примеру, 10 ДЗОТ на удалении в 5000 метров и для каждой из них предусмотреть 35 осколочно-фугасных гранат, то весьма высоки шансы на то, что от одной до трёх целей сохранят боеспособность. Это много! В связи с этим можно сделать некоторое «лирическое отступление» в область … компьютерных игр, а именно фэнтезийной походовой стратегии Heroes of Might and Magic V. Среди разнообразных существ, которыми игрок может укомплектовать своё войско (или встретить в составе вражеских сил), есть привидения. Волей разработчиков этой бесплотной нежити дарована способность в трети случаев не получать урон от противника. И по результатам тестов серии уклонений от ударов могли оказываться куда более длинными, чем три подряд таких события, при том, что сами привидения врагу урон наносят. Это существенно стало влиять на исход боя, поэтому создателям игры пришлось ограничить эту способность, добавив к ней условие «не больше трёх уклонений подряд». Естественно, что в реальной жизни такого «божественного вмешательства» ожидать не приходится и надо быть готовым, что в ряде случаев число попыток будет куда больше, чем подсказывает интуиция на базе житейского опыта. На том, собственно говоря, игорные дома с незапамятных времён и делают свой бизнес. А в артиллерии, если мы хотим гарантированно выполнить поставленную задачу, приходится полагаться на работу тыла и снабженческих служб, чтобы снарядов всегда хватало — поскольку даже при самой точной наводке и при самом лучшем уходе за орудием рассеивание неизбежно.

== Использованная литература ==
* Вентцель Е. С. ''Теория вероятностей'' — М.:Наука, 1964
* Никифоров Н. Н., Туркин П. И., Жеребцов А. А., Галиенко С. Г. ''Артиллерия'' / Под общ. ред. Чистякова М. Н. — М.: Воениздат МО СССР, 1953.
* Форсайт Дж., Малькольм М., Моулер К. ''Машинные методы математических вычислений'' / Пер. с англ. Икрамова Х. Д. — М.: Мир, 1980.
* ''152-мм гаубица обр. 1943 г. Руководство службы'' —  М.: Воениздат МО СССР, 1958
* ''Таблицы стрельбы 152-мм гаубицы обр. 1943 г.'' — М.: Воениздат МО СССР, 1968

[[Категория:Статьи по теоретическим основам артиллерии]]