RU2509319C1 - Способ радиолокационного определения времени окончания активного участка баллистической траектории - Google Patents

Abstract

Изобретение относится к устройствам траекторной обработки радиолокационной информации. Достигаемый технический результат изобретения - повышение вероятности определения времени окончания активного участка (АУТ) баллистической траектории за счет исключения измерений угла места и азимута из обрабатываемых выборок. Для этого на вход устройства определения времени окончания АУТ подают данные измерений дальности и радиальной скорости ракеты через одинаковые интервалы времени, равные периоду обзора радиолокационной станции, вычисляют произведения дальности на радиальную скорость, формируют фиксированную выборку типа «скользящего окна» значений произведений дальности на радиальную скорость, находят оценку скорости изменения произведения дальности на радиальную скорость путем оптимального взвешенного суммирования выборки значений произведений дальности на радиальную скорость, вычисляют среднеквадратическую ошибку (СКО) оценки, вычисляют отношение оценки скорости к СКО этой оценки, в каждом новом положении «скользящего окна» сравнивают отношение оценки скорости к СКО этой оценки с заданным порогом, решение об окончании АУТ принимают в момент времени, когда величина полученного отношения превышает заданный порог, величину которого выбирают в соответствии с требуемой вероятностью определения времени окончания АУТ. 3 ил., 6 табл.

Description

Изобретение относится к устройствам траекторной обработки радиолокационной информации и может быть использовано в РЛС и в автоматизированных системах управления радиолокационных подразделений.

Задачу определения времени окончания активного участка траектории (АУТ) необходимо решать для того, чтобы не допустить появления методических ошибок определения параметров баллистической траектории. В частности, координаты точки падения ракет малой и средней дальности полета могут определяться с недолетом или перелетом от нескольких десятков до нескольких сотен километров.

В качестве примера рассмотрим траекторию ракеты типа «Скад» с дальностью стрельбы 300 км. Параметры этой траектории приведены в таблице 1:

время от момента пуска t;

измеряемые полярные координаты: дальность r, угол места ε, азимут β, радиальная скорость $$\dot{r}$$

;

преобразованные координаты: высота z=rsinε, декартовы горизонтальные координаты х=rcosεsinβ и y=rcosεcosβ, произведение дальности на радиальную скорость $$r\dot{r}$$

;

модуль скорости V.

Активный участок траектории заканчивается на 65-ой секунде полета ракеты. Координаты ракеты измеряются в РЛС метрового диапазона (РЛС МДВ) «Резонанс-НЭ». Среднеквадратические ошибки (СКО) измерения координат в РЛС «Резонанс-НЭ»: дальности σ_r=300 м, радиальной скорости $${\sigma }_{\dot{r}}=\mathrm{1,5}$$

м/с, угла места и азимута σ_β=σ_ε=90 мин. РЛС находится в точке падения ракеты. (Вооружение ПВО и РЭС России. Альманах. М.: Издательство НО «Лига содействия оборонным предприятиям», 2011. - С.356-361).

Таблица 1
t, c	r, км	$$\dot{r}$$ , м/с	$$r\dot{r}$$ , м/с	ε, град	β, град	z, км	x, y, км	V, м/c	АУТ
35	297,83	-350,5	-104,39	0,16	45	0.52	210,566	598
40	295,84	-446,3	-132,03	0,67	45	3,46	209,145	729
45	293,35	-497,1	-145,82	1,28	45	6,55	207,349	877
50	290,28	-674,2	-195,71	2,0	45	10,13	205,101	1044
55	286,58	-808,7	-231,76	2,86	45	14,3	202,343	1234
60	282,17	-957,7	-270,23	3,87	45	19,04	199,039	1452
65	277,15	-1011	-280,14	5,04	45	24,35	195,189	1528
70	272,13	-994,1	-270,52	6,23	45	29,53	191,258	1496	ПУТ
75	267,2	-979,2	-261,64	7,4	45	34,41	187,337	1465
80	262,34	-966	-253,42	8.56	45	39,05	183,409	1436
85	257,54	-954,4	-245,8	9,72	45	43,48	179,465	1407

Для анализа величины методических ошибок прогноза точки падения ракеты вычислим оценки модуля скорости $$\widehat{V}$$

в различных точках на активном и на пассивном участках баллистической траектории по выборкам типа «скользящего окна» из пяти измерений декартовых координат x_i, y_i, z_i:

$$\widehat{V}=\sqrt{{\widehat{V}}_z^2+{\widehat{V}}_x^2+{\widehat{V}}_y^2}$$

где оценки скорости изменения высоты $${\widehat{V}}_z$$

, и горизонтальных координат $${\widehat{V}}_x$$

, $${\widehat{V}}_y$$

вычислялись по формулам:

$${\widehat{V}}_z=\frac{1}{T_0}{\displaystyle \sum _{i=1}^n{z}_i{\eta }_V\left(i\right)},$$

$${\widehat{V}}_x=\frac{1}{T_0}{\displaystyle \sum _{i=1}^n{x}_i{\eta }_V\left(i\right)},$$

$${\widehat{V}}_y=\frac{1}{T_0}{\displaystyle \sum _{i=1}^n{y}_i{\eta }_V\left(i\right)};$$

где

$${\eta }_V\left(i\right)=\frac{6}{n\left(n^2-1\right)\left(n^2-4\right)}\left[\left(n+1\right)\left(n+2\right)\left(6n-7\right)-2i\left(16n^2-19\right)+30i^2\left(n-1\right)\right]$$

- весовой коэффициент оценки скорости маневрирующей цели в реальном режиме времени;

Т₀ - период обзора РЛС;

n - количество измерений в выборке. (Кузьмин С.З. Цифровая обработка радиолокационной информации. - М: «Радио и связь», 1967, С.305-306).

Результаты расчетов, приведенные в таблице 2, показали, что модуль скорости и составляющие вектора скорости определяются с методическими ошибками, если в выборке имеются измерения, произведенные на АУТ, то есть при включенном ракетном двигателе. Ошибки устраняются, если в выборке отсутствуют измерения, произведенные на АУТ, и она состоит только из измерений, произведенных после выключения двигателя, то есть на пассивном участке траектории (ПУТ). В приведенном примере методические ошибки устраняются через 15 секунд после окончания АУТ. При увеличении длительности «скользящего окна» (объема выборки) будет соответственно увеличиваться время появления методических ошибок после окончания АУТ. При этом значения оценок скорости ракеты в точках, находящихся на АУТ, меньше истинных значений, а в точках, находящихся на ПУТ, больше истинных значений скорости. По этой причине дальность до точки падения, вычисленная по формуле $$L=\frac{V^2\sin 2\vartheta }{g}+zctg\vartheta$$

, будет определяться с недолетом, либо с перелетом. В приведенной формуле ϑ - угол наклона траектории, g - ускорение силы тяжести. (Жаков А.М., Пигулевский Ф.А. Управление баллистическими ракетами. - М: «Военное издательство», 1965, С.15).

Таблица 2
Участок	АУТ			ПУТ
Время полета, м/с	55	60	65	70	75	80	85	90
Оценка скорости, м/с	1216	1422	1614	1655	1559	1444	1407,4	1379,5
Истинная скорость, м/с	1234	1452	1528	1496	1465	1436	1407	1380
Методич. ошибка $$\delta {\widehat{V}}_n={\widehat{V}}_n-V_{ист}$$ , м/с	-18(-312)	-30(-106)	+94	+159	+88	+7,3	+0,4	-0,5

Примечание: в скобках приведена разность между оценкой модуля скорости на АУТ и значением максимальной скорости (М_макс=1528 м/с) в конце АУТ.

В приведенном примере при измерении скорости ракеты за 5 секунд до окончания АУТ, то есть на 60-й секунде полета, координаты точки падения будут определяться с недолетом около 47-ми км. Если измерять скорость через 5 секунд после выключения двигателя, то перелет будет равен 45-ти км.

Таким образом, для исключения методических ошибок прогноза точки падения определение (оценивание) начальных параметров движения ракеты должно производиться после выключения ракетного двигателя, то есть в точке, находящейся на пассивном участке траектории, а измерения координат ракеты, произведенные до выключения ракетного двигателя, то есть на АУТ, не должны использоваться. Поэтому выборки значений измеренных координат, по которым вычисляются параметры движения ракеты (модуль скорости, угол наклона траектории к горизонту, курс и три координаты), должны формироваться после определения времени окончания АУТ.

Известны способы определения времени окончания АУТ средствами разведки инфракрасного и оптического диапазона по факелу ракетного двигателя. (Колгашкин Ю.Г. Комплексы самолетного базирования для обнаружения стартующих БР средней и малой дальности. / Международная конференция по проблемам глобальной защиты от баллистических ракет. М.: МАК «Вымпел» 1993. С.126-128).

Из радиолокационных способов аналогами заявляемому способу являются способы определения момента окончания маневра. Это связано с тем, что на АУТ ракета совершает маневр большой интенсивности. Вначале ракета движется вертикально вверх. Затем, на участке выведения, ракета движется по дуге и угол наклона траектории изменяется от 90° до расчетного значения для заданной дальности полета. Двигатель выключается, когда скорость достигнет заданной величины. Поэтому на АУТ скорости изменения высоты и горизонтальных декартовых координат увеличиваются с переменными ускорениями, величины которых зависят от силы тяги двигателя и в несколько раз больше ускорения силы тяжести. На ПУТ ракета становится неманеврирующей целью, так как летит по баллистической кривой с постоянным вертикальным ускорением, равным ускорению силы тяжести. Величина скорости изменения горизонтальных декартовых координат практически не изменяется, так как горизонтальные составляющие ускорения примерно равны нулю. Таким образом, устройства, реализующие способы определения времени окончания АУТ, и устройства, реализующие способы обнаружения маневра, должны решать одинаковую задачу - обнаружение факта изменения характера движения ракеты. (Жаков А.М., Пигулевский Ф.А. Управление баллистическими ракетами. - М: «Военное издательство», 1965, С.10-11).

Известны способы обнаружения времени окончания маневра путем сравнения оценок ускорения декартовых координат со среднеквадратической ошибкой (СКО) этих оценок. (Кузьмин С.З. Цифровая обработка радиолокационной информации. - М: «Радио и связь», 1967, С.310-311). Решение об окончании маневра и, следовательно, об окончании АУТ, принимается в момент времени, когда оценки ускорения декартовых координат становятся меньше СКО этих оценок. Основным недостатком этих способов являются высокие требования к точности измерения угла места и азимута. Поэтому в РЛС, размеры антенн которых соизмеримы с длиной волны и величины ошибок измерения угла места и азимута достигают значений до нескольких градусов, оценки ускорений меньше СКО как на активном, так и на пассивном участке и определить границу между ними практически невозможно.

Наиболее близким по своей сущности к заявляемому способу, то есть прототипом, является способ определения времени окончания АУТ по абсолютной величине приращения скорости изменения декартовой координаты. (Кузьмин С.З. Цифровая обработка радиолокационной информации. - М: «Радио и связь», 1967, С.347). Для этого находят оценки скорости изменения декартовой координаты, например координаты х, в текущем обзоре $${\widehat{V}}_x\left(N\right)$$

и в предыдущем обзоре $${\widehat{V}}_x\left(N-1\right)$$

по одинаковым выборкам типа «скользящего окна» значений координаты и вычисляют их разность, то есть абсолютную величину приращения скорости:

$$\left|\Delta V_x\right|=\left|{\widehat{V}}_x\left(N\right)-{\widehat{V}}_x\left(N-1\right)\right|$$

.

Затем делят эту абсолютную величину приращения скорости на среднеквадратическую ошибку оценки скорости $${\sigma }_{V_x}$$

и сравнивают с некоторым наперед заданным порогом или числом a в каждом новом положении «скользящего окна».

Таким образом, устройство определения времени окончания АУТ способом-прототипом представляет собой пороговое устройство и работает по следующему алгоритму:

если $$\frac{\left|\Delta V_x\right|}{{\sigma }_{V_x}}>a$$

- ракета находится на активном участке траектории, то есть имеется маневр;

если $$\frac{\left|\Delta V_x\right|}{{\sigma }_{V_x}}\le a$$

- ракета находится на пассивном участке траектории, то есть маневр отсутствует.

По аналогичному правилу определяется время окончания АУТ при использовании абсолютного приращения оценки скорости изменения второй горизонтальной декартовой координаты $$\left|\Delta V_y\right|=\left|{\widehat{V}}_y\left(N\right)-{\widehat{V}}_y\left(N-1\right)\right|$$

:

если $$\frac{\left|\Delta V_y\right|}{{\sigma }_{V_y}}>a$$

- ракета находится на активном участке траектории, то есть имеется маневр;

если $$\frac{\left|\Delta V_y\right|}{{\sigma }_{V_y}}\le a$$

- ракета находится на пассивном участке траектории, то есть маневр отсутствует.

По своей сущности данный алгоритм является алгоритмом выбора гипотезы модели движения ракеты. Гипотеза 1 - ракета на АУТ, то есть маневрирующая цель. Вероятность правильного выбора гипотезы 1 равна вероятности того, что абсолютное приращение скорости изменения декартовой координаты больше СКО оценки этой скорости $$\left|\Delta V_y\right|\ge {\sigma }_{V_y}$$

. Данная вероятность рассчитывается по следующей формуле:

$$p_{АУТ}\left(\left|\Delta V\right|\ge {\sigma }_{\widehat{V}}\right)=p_{АУТ}\left(\frac{\left|\Delta V\right|}{\sigma \widehat{V}}\right)={\displaystyle \underset{-\frac{\left|\Delta V\right|}{{\sigma }_{\widehat{V}}}}{\overset{\frac{\left|\Delta V\right|}{{\sigma }_{\widehat{V}}}}{\int }}{e}^{-t^2/2}}dt$$

Вероятность правильного выбора гипотезы 2 равна вероятности того, что $$\left|\Delta V_y\right|\le {\sigma }_{V_y}$$

, рассчитывается по следующей формуле:

$$p_{ПУТ}\left(\left|\Delta V\right|\le {\sigma }_{\widehat{V}}\right)=p_{ПУТ}\left(\frac{{\sigma }_{\widehat{V}}}{\left|\Delta V\right|}\right)={\displaystyle \underset{-\frac{{\sigma }_{{}_{\widehat{V}}}}{\left|\Delta V\right|}}{\overset{\frac{{\sigma }_{{}_{\widehat{V}}}}{\left|\Delta V\right|}}{\int }}{e}^{-t^2/2}}dt$$

Вероятность определения времени окончания АУТ будет равна минимальной величине вышеприведенных вероятностей, то есть $$p_{АУ{Т}_{МИН}}$$

либо $$p_{ПУ{Т}_{МИН}}$$

. (Бронштейн И.Н., Семендяев К.А. Справочник по математике для инженеров и учащихся ВТУЗОВ. М:. «Наука», 1980, С.92-93).

Число значений декартовых координат в обрабатываемой выборке (в «скользящем окне») не должно превышать 4-5, чтобы не увеличивалось запаздывание выявления факта окончания АУТ.

Так как скорости изменения горизонтальных декартовых координат практически не изменяются на пассивном участке, то их оценки и СКО этих оценок вычисляются по формулам для линейной траектории:

,

.

Весовые коэффициенты оценки скорости на линейной траектории вычисляются по формуле: $${\eta }_V\left(i\right)=\frac{12i-6-6i}{n\left(n^2-1\right)}$$

(Кузьмин С.З. Цифровая обработка радиолокационной информации. - М.: «Радио и связь», 1967, С.301).

Значения СКО оценок вычисляются по формулам:

, $${\sigma }_{\widehat{V}}{}_y=\frac{{\sigma }_y}{T_0}\sqrt{\frac{12}{n\left(n^2-1\right)}}$$

(Там же, С.308).

Формулы вычисления СКО измерения координат имеют вид:

$${\sigma }_x=\sqrt{{\left(r\cos \epsilon \cos \beta {\sigma }_{\beta }\right)}^2+{\left(r\sin \epsilon \sin \beta {\sigma }_{\epsilon }\right)}^2+{\left(\cos \epsilon \sin \beta {\sigma }_r\right)}^2}$$

,

$${\sigma }_y=\sqrt{{\left(r\cos \epsilon \sin \beta {\sigma }_{\beta }\right)}^2+{\left(r\sin \epsilon \cos \beta {\sigma }_{\epsilon }\right)}^2+{\left(\cos \epsilon \cos \beta {\sigma }_r\right)}^2}$$

.

Ошибки измерения дальности σ_r можно не учитывать, так как их влияние несущественно. Ошибки измерения угла места σ_ε оказывают существенное влияние на СКО оценок, если высота БЦ z_ср=r_срsinε_ср соизмерима с величиной горизонтальных декартовых координат. В остальных случаях доминирующее влияние оказывают ошибки измерения азимута σ_β. При этом СКО σ_x, а следовательно, и $${\sigma }_{\widehat{V}}{}_x$$

, принимают максимальные значения при азимуте, равном нулю, и минимальные значения - при азимуте, равном 90°. СКО $${\sigma }_{\widehat{V}}{}_y$$

, наоборот, максимальны при азимуте 90° и минимальны при азимуте 0°. При азимуте, равном 45°, СКО одинаковы, то есть $${\sigma }_{\widehat{V}}{}_x={\sigma }_{\widehat{V}}{}_y$$

. Поэтому для определения времени окончания АУТ при азимутах БЦ меньше 45° целесообразно использовать СКО $${\sigma }_{\widehat{V}}{}_y$$

, а при азимутах больше 45° - СКО $${\sigma }_{\widehat{V}}{}_x$$

.

В качестве примера в таблицах 3-5 приведены значения абсолютного приращения оценок скорости изменения декартовой координаты y|ΔV_y| за период обзора РЛС, СКО $${\sigma }_{\widehat{V}}{}_y$$

и величины отношений $$\frac{\left|\Delta V_y\right|}{{\sigma }_{\widehat{V}}{}_y}$$

, рассчитанные по выборкам типа «скользящего окна» из 3-х значений декартовой координаты y при ошибках измерения угловых координат в РЛС МДВ «Резонанс-НЭ», равным σ_ε=σ_β=90 мин. В таблице 3 приведены результаты расчетов при азимуте 10°, в таблице 4 - при азимуте 45°, в таблице 5 - при азимуте 0°. В качестве исходных данных использовались данные таблицы 1. За начало отсчета времени принят момент окончания АУТ на 65-й секунде полета ракеты.

Как видно из данных таблицы 3 и 4, при азимутах ракеты от 10° до 45° величина приращений оценок скорости изменения декартовой координаты y за период обзора РЛС |ΔV_y| меньше СКО оценок скорости $${\sigma }_{\widehat{V}}{}_y$$

как на активном, так и на пассивном участках, поэтому граница между ними, а значит и время окончания АУТ, определяется с низкой вероятностью. Вероятность $$p_{АУТ}\left(\left|\Delta V\right|\ge {\sigma }_{\widehat{V}}\right)$$

и, следовательно, вероятность правильного определения времени окончания АУТ при азимуте 10° не превышает 0,5, а при больших азимутах приближается к нулю. Аналогичное соотношение |ΔV_x| и $${\sigma }_{\widehat{V}}{}_x$$

наблюдается при азимутах от 45° до 80°.

Повышать точность оценок скорости за счет увеличения объема выборок нецелесообразно, так как длительность «скользящего окна» становится соизмеримой с продолжительностью АУТ ракет малой и средней дальности. Кроме того, увеличивается запаздывание выявления факта окончания АУТ.

Таблица 3
Участок	Активный участок			Пассивный участок
Время, сек	-10	0		+5	+10	+15	+20
|ΔVy|, м/с	134	153		87	10	1	3
$${\sigma }_{\widehat{V}}{}_y$$ при β=10°, м/с	187	195		201	209	215	225
$$\frac{\left|\Delta V_y\right|}{{\sigma }_{\widehat{V}}{}_y}$$	0,7	0,8		0,4	0,05	0,005	0,013
Таблица 4
Участок	Активный участок			Пассивный участок
Время, сек	-10		0	+5	+10	+15	+20
|ΔVy|, м/с	96		110	63	7	0,7	1,5
$${\sigma }_{\widehat{V}}{}_y$$ при β=45°, м/с	750		726	707	700	686	672
$$\frac{\left|\Delta V_y\right|}{{\sigma }_{\widehat{V}}{}_y}$$	0,13		0,15	0,09	0,01	0,001	0,02

Задача решается только при азимутах близких 0° или 90. В этом случае на величину СКО оценки скорости изменения горизонтальных декартовых координат доминирующее влияние оказывают ошибки измерения угла места.

Таблица 5
Участок	Активный участок		Пассивный участок
Время, сек.	-10	0	+5	+10	+15	+20
|ΔVy|, м/с	136	155	116	17	1	3
$${\sigma }_{\widehat{V}}{}_y$$ при β=0°, м/с	67	99	116	126	143	159
$$\frac{\left|\Delta V_y\right|}{{\sigma }_{\widehat{V}}{}_y}$$	2	1,5	1	0,13	0,007	0,02

На графиках фиг.1 показаны значения |ΔV_y| и $${\sigma }_{\widehat{V}}{}_y$$

, вычисленные по данным измерений в РЛС «Резонанс - НЭ», в различных точках траектории ракеты «Скад» при азимутах 0° и 10°.

Таким образом, основным недостатком прототипа, являются низкие вероятности определения времени окончания АУТ при грубых измерениях азимута и угла места. Поэтому в РЛС МДВ, либо в других РЛС, размеры антенны которых соизмеримы с длиной волны, использовать данный способ практически невозможно.

Техническим результатом настоящего изобретения является разработка нового способа, при использовании которого повышается вероятность определения времени окончания АУТ за счет исключения измерений азимута и угла места из обрабатываемых выборок.

В предлагаемом изобретении вычисляют оценки скорости изменения произведения дальности r_i на радиальную скорость $${\stackrel{\cdot }{r}}_i$$

путем оптимального

взвешенного суммирования фиксированной выборки типа «скользящего окна» значений произведений дальности на радиальную скорость $$r_i{\dot{r}}_i$$

:

$$\widehat{\dot{U}}\left(r_i{\dot{r}}_i\right)=\frac{1}{T_0}{\displaystyle \sum _{i=1}^n{r}_i{\dot{r}}_i{\eta }_V\left(i\right)}.$$

При анализе различных типов баллистических траекторий была выявлена следующая закономерность: оценки $$\stackrel{\hat{\cdot }}{U}\left(r_i{\stackrel{\cdot }{r}}_i\right)$$

отрицательны на активном участке траектории и положительны на пассивном участке. Данная закономерность показана на графике фиг.2. Поэтому, так же, как в прототипе, вычисляются СКО оценок:

$${\sigma }_{\widehat{\dot{U}}\left(r_i{\dot{r}}_i\right)}=\frac{\sigma \left(r_i{\dot{r}}_i\right)}{T_0}\sqrt{\frac{12}{n\left(n^2-1\right)}},$$

$$\sigma \left(r_i{\dot{r}}_i\right)=\sqrt{{\left(r_i{\sigma }_{\dot{r}}\right)}^2+{\left({\dot{r}}_i{\sigma }_r\right)}^2},$$

где r_i, $${\dot{r}}_i$$

- измеренные значения дальности и радиальной скорости;

σ_r, $${\sigma }_{\dot{r}}$$

- среднеквадратические ошибки измерения дальности и радиальной скорости.

Полученные оценки $$\widehat{\dot{U}}\left(r_i{\dot{r}}_i\right)$$

делят на величину среднеквадратической ошибки этих оценок $${\sigma }_{\widehat{\dot{U}}\left(r_i{\dot{r}}_i\right)}$$

. В каждом новом положении «скользящего окна» сравнивают полученные отношения $$\frac{\widehat{\dot{U}}\left(r_i{\dot{r}}_i\right)}{{\sigma }_{\widehat{\dot{U}}\left(r_i{\dot{r}}_i\right)}}$$

с порогом а. Решение об окончании АУТ и о начале ПУТ принимают в момент времени, когда отношение $$\frac{\widehat{\dot{U}}\left(r_i{\dot{r}}_i\right)}{{\sigma }_{\widehat{\dot{U}}\left(r_i{\dot{r}}_i\right)}}$$

становится больше числа а:

если $$\frac{\widehat{\dot{U}}\left(r_i{\dot{r}}_i\right)}{{\sigma }_{\widehat{\dot{U}}\left(r_i{\dot{r}}_i\right)}}<a$$

- ракета на активном участке траектории;

если $$\frac{\widehat{\ddot{U}}\left(r_i^2\right)}{{\sigma }_{\widehat{\ddot{U}}\left(r_i^2\right)}}\ge a$$

- ракета на пассивном участке траектории.

В отличие от прототипа, величина СКО оценки зависит только от ошибок измерения радиальной скорости и дальности. При этом доминирующее влияние оказывают ошибки измерения радиальной скорости. Ошибки измерения радиальной скорости и дальности не зависят от размеров антенны и могут быть уменьшены до нескольких метров в секунду и до нескольких десятков или сотен метров.

В качестве примера в таблице 6 приведены значения оценок скорости изменения произведения дальности на радиальную скорость $$\widehat{\dot{U}}\left(r_i{\dot{r}}_i\right)$$

, СКО их определения $${\sigma }_{\widehat{\dot{U}}\left(r_i{\dot{r}}_i\right)}$$

по выборкам из 3-х измерений радиальной скорости и дальности, произведенных в РЛС метрового диапазона «Резонанс-НЭ» ( $${\sigma }_{\dot{r}}=\mathrm{1,5}$$

м/с, σ_r=300 м), а также величины отношений $$\frac{\left|\widehat{\dot{U}}\left(r_i{\dot{r}}_i\right)\right|}{{\sigma }_{\widehat{\dot{U}}\left(r_i{\dot{r}}_i\right)}}.$$

Таблица 6
Участок	Активный участок			Пассивный участок
Время, сек	-10	-5	0	+5	+10	+15	+20	+25
$$\widehat{\dot{U}}\left(r_i{\dot{r}}_i\right),$$ км2/с2	-8,59	-7,45	-4,84	-0,03	+1,83	+1,71	+1,58	+1,47
$${\sigma }_{\widehat{\dot{U}}\left(r_i{\dot{r}}_i\right)},$$ км2/с2	0,07	0,07	0,07	0,07	0,07	0,07	0,07	0,07
$$\frac{\left|\widehat{\dot{U}}\left(r_i{\dot{r}}_i\right)\right|}{{\sigma }_{\widehat{\dot{U}}\left(r_i{\dot{r}}_i\right)}}$$	123	106	70	0,4	26	24	22	21

Как видно из данных таблицы 6, значения оценок скорости изменения произведения дальности на радиальную скорость $$\widehat{\dot{U}}\left(r_i{\dot{r}}_i\right)$$

меньше нуля (отрицательны) на активном участке и больше нуля (положительны) на пассивном участке траектории. При этом оценки $$\widehat{\dot{U}}\left(r_i{\dot{r}}_i\right)$$

на пассивном участке траектории превышают значения СКО $${\sigma }_{\widehat{\dot{U}}\left(r_i{\dot{r}}_i\right)}$$

более чем в двадцать раз. На АУТ абсолютная величина оценок более чем в сто раз превышает величину СКО этих оценок. Поэтому, даже при увеличении ошибок измерения радиальной скорости в 10 раз, вероятности правильного выбора гипотезы 1, то есть $$p_{АУТ}\left(\frac{\left|\widehat{\dot{U}}\left(r_i{\dot{r}}_i\right)\right|}{{\sigma }_{\widehat{\dot{U}}\left(r_i{\dot{r}}_i\right)}}\right)$$

, и гипотезы 2, то есть $$p_{ПУТ}\left(\frac{\left|\widehat{\dot{U}}\left(r_i{\dot{r}}_i\right)\right|}{{\sigma }_{\widehat{\dot{U}}\left(r_i{\dot{r}}_i\right)}}\right)$$ , практически равны единице. Поэтому величину порога а можно выбирать в интервале от 2-х до 3-х и время окончания АУТ будет определяться также с вероятностью, близкой единице.

Принцип работы устройства определения времени окончания АУТ заявляемым способом по выборке из трех значений произведений дальности на радиальную скорость поясняется схемой, приведенной на фиг.3. В состав устройства входят блок 1 преобразования измеренных полярных координат, то есть вычисления произведений дальности на радиальную скорость, блок 2 оценивания первого приращения произведений дальности на радиальную скорость за период обзора, блок 3 деления оценки первого приращения на период обзора и вычисления оценки скорости изменения произведений дальности на радиальную скорость, блок 4 вычисления среднеквадратической ошибки оценки скорости изменения произведения дальности на радиальную скорость, блок 5 деления оценки скорости на среднеквадратическую ошибку и пороговое устройство 6.

На вход блока 1 устройства подаются через равные промежутки времени, равные периоду обзора РЛС, данные измерений радиальной скорости и дальности, и вычисляются произведения дальности на радиальную скорость $$r_i{\dot{r}}_i$$

. Значение произведения дальности на радиальную скорость в текущем обзоре $$r_1{\dot{r}}_1$$

подается на вход первой линии задержки блока 2 оценивания первого приращения произведений дальности на радиальную скорость, а в предыдущем обзоре $$r_2{\dot{r}}_2$$

не используются, так как весовой коэффициент равен нулю. Значение произведения дальности на радиальную скорость $$r_3{\dot{r}}_3$$

, задержанное на время, равное двум периодам обзора (2Т₀), умножается на весовой коэффициент, равный (+0,5) и одновременно с текущим взвешенным значением произведения дальности на радиальную скорость поступает на вход сумматора. Подобная схема для оптимального оценивания параметров приведена на 303-й странице упомянутой монографии Кузьмина С.З. «Цифровая обработка радиолокационной информации». При увеличении количества измерений в «скользящем окне» до числа n необходимо использовать (n-1) линий задержки, умножителей, а также заранее вычислить соответствующие числу n весовые коэффициенты. При реализации способа на цифровой вычислительной машине (ЦВМ) вычисление и запоминание значений произведений дальности на радиальную скорость производится в оперативном запоминающем устройстве, заранее вычисленные весовые коэффициенты хранятся в долговременном запоминающем устройстве, а операции суммирования и умножения производятся в арифметическом устройстве. Сигнал с выхода сумматора поступает на вход делителя 3, где производится деление на период обзора, в результате чего получается оценка скорости изменения произведения дальности на радиальную скорость. Полученная оценка делится на величину среднеквадратической ошибки определения этой оценки, вычисленной в блоке 4 по данным измерений дальности и радиальной скорости. Полученный результат подается на пороговое устройство 6, где сравнивается с порогом, величина которого зависит от заданной вероятности определения времени окончания АУТ.

При превышении величины отношения $$\frac{\widehat{\dot{U}}\left(r_i{\dot{r}}_i\right)}{{\sigma }_{\widehat{\dot{U}}\left(r_i{\dot{r}}_i\right)}}$$

заданного порога принимается решение об окончании активного участка и о начале пассивного участка баллистической траектории. После принятия такого решения информация о времени окончания АУТ выдается потребителем на устройства (системы) определения параметров движения ракеты по выборкам радиолокационных измерений, произведенных на пассивном участке траектории, то есть после выключения ракетного двигателя. В результате этого параметры движения ракеты и ее траектория рассчитываются без методических ошибок, появляющихся при наличии в обрабатываемых выборках измерений, произведенных на АУТ.

Заявляемое изобретение соответствует условиям новизны и изобретательского уровня. Признаками изобретения, совпадающими с признаками прототипа, являются операции преобразования измеренных полярных координат, вычисление оценок скорости изменения преобразованных координат и значений СКО этих оценок, деление оценок на СКО и принятие решения об окончании активного участка траектории по результатам сравнения величины отношения оценки скорости к СКО этой оценки с величиной заданного порога.

Новизна заключается в следующем:

формируются фиксированные выборки произведений дальности на радиальную скорость, а не выборки значений декартовых координат;

вычисляются оценки скорости изменения произведения дальности на радиальную скорость, а не оценки скорости изменения горизонтальных декартовых координат;

с порогом сравниваются значения отношения оценки скорости к СКО этой оценки, а не отношения абсолютной величины приращения оценок скоростей изменения декартовых координат в соседних обзорах к СКО оценки скорости;

решение об окончании АУТ принимают при превышении заданного порога величиной отношения оценки скорости к СКО этой оценки, а не в момент времени, когда отношение абсолютной величины приращения оценок скоростей изменения декартовых в соседних обзорах к СКО оценки скорости становится меньше заданного порога.

Промышленная применимость заявляемого изобретения подтверждается возможностью определения времени окончания АУТ с помощью РЛС типа «Резонанс-НЭ», находящихся на вооружении ПВО, и других станций, измеряющих радиальную скорость ракеты с ошибками до 10-20 м/с.

Использование предлагаемого способа радиолокационного определения времени окончания активного участка позволит повысить вероятность определения времени окончания АУТ и устранить методические ошибки определения параметров движения ракеты, экстраполяции баллистической траектории и прогноза точки падения ракеты за счет исключения измерений дальности, азимута, угла места, радиальной скорости, произведенных на активном участке, из состава обрабатываемых выборок.

Claims (1)

1. Способ радиолокационного определения времени окончания активного участка баллистической траектории, заключающийся в том, что производят преобразование измеренных полярных координат баллистической ракеты, формируют фиксированную выборку типа «скользящего окна» значений преобразованных координат, вычисляют оценку скорости изменения преобразованных координат и ее среднеквадратическую ошибку, вычисляют отношение оценки скорости к среднеквадратической ошибке этой оценки, решение о времени окончания активного участка принимают по результатам сравнения полученного отношения оценки скорости к среднеквадратической ошибке этой оценки в каждом новом положении «скользящего окна» с порогом, соответствующим заданной вероятности, отличающийся тем, что при преобразовании полярных координат вычисляют произведения измеренных значений дальности на измеренные значения радиальной скорости, формируют фиксированную выборку полученных произведений, оценку скорости изменения произведения дальности на радиальную скорость находят путем оптимального взвешенного суммирования выборки значений произведений дальности на радиальную скорость, а решение об окончании активного участка траектории принимают в момент времени, когда значение отношения оценки скорости изменения произведения дальности на радиальную скорость к среднеквадратической ошибке этой оценки становится больше заданного порога.

Images