RU2592747C1 - Способ стрельбы по движущейся цели - Google Patents
Способ стрельбы по движущейся цели Download PDFInfo
- Publication number
- RU2592747C1 RU2592747C1 RU2015118108/28A RU2015118108A RU2592747C1 RU 2592747 C1 RU2592747 C1 RU 2592747C1 RU 2015118108/28 A RU2015118108/28 A RU 2015118108/28A RU 2015118108 A RU2015118108 A RU 2015118108A RU 2592747 C1 RU2592747 C1 RU 2592747C1
- Authority
- RU
- Russia
- Prior art keywords
- projectile
- vector
- target
- trajectory
- miss
- Prior art date
Links
Images
Landscapes
- Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)
Abstract
Изобретение относится к технике стрельбы по двигающимся целям и может использоваться в системах обнаружения и определения траектории полета поражающих целей. Технический результат - повышение точности. Для этого определяют точку встречи и момента встречи цели и снаряда, на основании оценок координат положения цели формируют начальную опорную траекторию движения снаряда и вектор промаха снаряда, осуществляют линеаризацию траектории снаряда в окрестности опорной траектории по углам прицеливания, на основании которой формируют матрицу частных производных вектора положения снаряда по углам прицеливания, которая удовлетворяет системе линейных однородных дифференциальных уравнений; при превышении длины вектора промаха своего максимально возможного допустимого порогового значения с использованием матрицы частных производных осуществляют коррекцию углов прицеливания артиллерийского орудия и повторяют формирование опорной траектории снаряда и вектора промаха, а при отсутствии превышения длины вектора промаха своего максимально возможного допустимого порогового значения осуществляют стрельбу по цели, используя последние величины углов прицеливания. 1 ил.
Description
Изобретение относится к технике стрельбы по движущимся целям и может использоваться в системах обнаружения, определения траектории и поражения целей.
Известен способ прицеливания при атаке скоростных целей истребителем по спрямленной траектории, заключающийся в том, что цель сопровождают с помощью бортовых радиолокационных или/и оптических систем и при этом определяют текущие значения векторов ее относительной дальности, скорости и линейного ускорения, прогнозируют вектор дальности упрежденного положения цели на суммарное время выхода истребителя на границу области эффективного применения оружия и полета снаряда, пущенного с границы этой области до встречи с целью, направляют полет истребителя в процессе атаки в точку выхода истребителя на границу области эффективного применения оружия и полета снаряда (RU 2170907 C2).
Известен способ стрельбы управляемым снарядом с лазерной полуактивной головкой самонаведения по движущейся цели, включающий топографическую привязку целеуказателя и огневой позиции к местности, обнаружение цели целеуказателем. В этом способе задают набор вероятных маршрутов движения цели и, после обнаружения цели, выбирают ближайший к измеренным координатам цели маршрут, с учетом которого осуществляют расчет прогнозируемой точки встречи снаряда с целью (RU 2347999 C2).
Известен способ определения точности наведения и сближения снаряда с целью по наблюдаемым параметрам их траекторного движения, в котором с помощью станций сопровождения авиационно-космических объектов в едином базисе измеряют текущие значения векторов дальности до снаряда и цели, вычитают из первого второй и определяют текущее значение вектора дальности до снаряда относительно цели, по измеренным в момент начала наблюдения векторам дальности цели и снаряда определяют значения векторов их скорости, на основе полученных данных прогнозируют траектории их полета до расчетной точки встречи (RU 2355990 C2).
Однако данные известные способы не являются универсальными. Каждый из этих способов имеет свою ограниченную область применения, при выходе за рамки которой расчетная точка встречи или не будет определяться, или будет определяться некорректно.
Технический результат заявленного способа заключается в повышении точности стрельбы артиллерийской установки за счет обеспечения возможности формирования траектории движения снаряда, учитывающей функциональную зависимость движения от начальных условий полета снаряда.
Для этого в способе стрельбы по цели, сопровождаемой радиолокационной системой, и двигающейся по произвольной траектории, описываемой уравнением y=g(t), у∈ℝ3, поражающим снарядом, двигающимся по траектории, описываемой системой обыкновенных дифференциальных уравнений
, x(0)=х0, x∈ℝn, заключающемся в определении точки встречи и момента встречи цели и снаряда, когда координаты положения цели и снаряда совпадают xi(Т)-yi(Т)=0, i=1, 2, 3, при этом на основании оценок координат положения цели, путем приближенного задания начальных углов прицеливания артиллерийского орудия, формируют начальную опорную траекторию движения снаряда и вектор промаха снаряда, осуществляют линеаризацию траектории снаряда в окрестности опорной траектории по углам прицеливания, на основании которой формируют матрицу частных производных вектора положения снаряда по углам прицеливания, которая удовлетворяет системе линейных однородных дифференциальных уравнений; при превышении длины вектора промаха своего максимально возможного допустимого порогового значения с использованием матрицы частных производных осуществляют коррекцию углов прицеливания артиллерийского орудия и повторяют формирование опорной траектории снаряда и вектора промаха, а при отсутствии превышения длины вектора промаха своего максимально возможного допустимого порогового значения осуществляют стрельбу по цели, используя последние величины углов прицеливания.
Здесь y - вектор, состоящий из 3-х компонент y1, y2, y3, соответствующих положению цели по каждой оси прямоугольной системы координат;
t - время;
g - вектор-функция, состоящая из 3-х компонент gi(t), g2(t), g3(t), описывающих движение цели вдоль каждой оси прямоугольной системы координат,
ℝ3 - 3-мерное Евклидово пространство;
x - вектор фазовых координат снаряда, первые 3 координаты которого описывают положение снаряда по осям прямоугольной системы координат,
х0 - вектор начальных условий движения снаряда,
f - вектор-функция, описывающая скорость изменения фазовых координат снаряда,
n - количество координат (длина) вектора фазовых координат снаряда.
Существо заявленного способа состоит в следующем.
Необходимо определить точку встречи цели и поражающего снаряда. Рассматриваемый способ разработан для снарядов, динамика движения которых в процессе полета может явным образом зависеть от начальных условий, например начального угла бросания. Так, начальный угол бросания может оказывать влияние на траекторию штатного снаряда как непосредственно в момент выстрела, так и на всем протяжении полета.
Считаем, что движение снаряда описывается системой обыкновенных нелинейных дифференциальных уравнений, а движение цели известно и описывается гладкими функциями времени:
Фазовый вектор x(t)∈ℝn, описывающий движение снаряда, удовлетворяет системе обыкновенных дифференциальных уравнений:
Часть компонент вектора начальных условий движения снаряда задана, а часть не задана и подлежит определению (например, углы прицеливания):
a∈ℝn-m - заданная часть начальных условий, a z∈ℝm - незаданная часть начальных условий фазового вектора, которая может варьироваться. Знак Т в верхней индексации вектора и матрицы обозначает операцию транспонирования.
Задача состоит в нахождении такого вектора начальных условий z, при котором в некоторый момент времени T траектории x(t) и y(t) пересекутся:
где i=1, 2, 3.
Предположим, что первые 3 координаты фазовых векторов снаряда и цели соответствуют векторам положения этих точек в 3-мерном евклидовом пространстве ℝ3. На вектор z могут быть наложены дополнительные ограничения в форме равенств:
Задача (2)-(5) является классической нелинейной краевой задачей.
При фиксированном векторе Z фазовый вектор-функция x является решением задачи Коши и в каждый момент времени является функцией своих начальных условий. Поэтому, если k+3=m+1, то эта задача сводится к решению системы k+3 нелинейных уравнений (4), (5) относительно m+1 неизвестных (m неизвестных начальных условий плюс время встречи).
Метод Ньютона, применяемый к решению этой системы уравнений, предполагает на каждом шаге вычисление матрицы частных производных фазового вектора решения системы дифференциальных уравнений (2) (матрицы Якоби) по начальным условиям z. Матрица частных производных Н(t) вектор-функции решения этой системы по начальным условиям z
удовлетворяет системе линейных однородных дифференциальных уравнений в вариациях:
Численное интегрирование системы дифференциальных уравнений (6), (7) одновременно с системой дифференциальных уравнений (2) позволяет организовать итерационный вычислительный процесс метода Ньютона.
Если правая часть дифференциального уравнения, описывающего движение точки 1, явно зависит от вектора начальных условий Z, то задача о точке встречи имеет некоторые специфические особенности. В этом случае дифференциальное уравнение (2) принимает вид:
где x(t)∈ℝn, z∈ℝm, и уравнение в вариациях для матрицы Якоби решения этой системы отличается от уравнения (6).
В любом случае краевая задача о точке встречи сводится к системе нелинейных уравнений относительно вектора z и момента встречи t.
где Ф(z,t)∈ℝk+3, z∈ℝm.
x1(z,t), x2(z,t), x3(z,t) - первые 3 компоненты вектора x(t) - решения системы дифференциальных уравнений (8) с начальными условиями (3).
Далее для решения системы уравнений применяется метод Ньютона (9), при этом (q+1)-й шаг итерационного вычисления имеет вид:
где J(z(q), t(q)) - матрица Якоби вектор-функции Ф:
Jφ(z) - матрица Якоби вектор-функции φ(z), (Jφ(z))i,j=∂φi/∂zj,
0k×1 - нулевой вектор длиной k,
Для вычисления матрицы J(z, t) необходимо иметь выражение для градиента функции xi(z, t) по переменным z, t. Частная производная функции xi(z, t) определяется следующим образом:
где x(t) - решение системы дифференциальных уравнений (8).
Далее вычисляют векторы частных производных вектор-функции x(z, t) по начальным условиям z1, …, zm. Для определенности вычисляют вектор и формулируют обобщение системы дифференциальных уравнений (6), (7).
При фиксированных начальных условиях z вектор hi(t) является решением следующей линейной неоднородной системы дифференциальных уравнений:
где х(t) - решение системы дифференциальных уравнений (8) с начальными условиями (2). Через δj обозначен j-й орт.
Для того чтобы найти градиент ∇zxi(z,t) по переменным z1, …, zm надо проинтегрировать m систем линейных неоднородных дифференциальных уравнений на интервале времени [0; t], получив при этом m векторов столбцов hi(t) - решений этих систем в момент времени t. Тогда:
Такие вычисления дают возможность организации итерационного вычислительного процесса метода Ньютона.
На прилагаемой фигуре показана принципиальная блок-схема алгоритма расчета точки встречи.
В блоке 1 осуществляется формирование начального приближения, в блоке 2 - расчет вектор-функции Ф и ее матрицы Якоби J, в блоке 3 - формирование вектора поправки начальных условий траектории снаряда, в блоке 4 осуществляется поправка траектории снаряда, в блоке 5 - расчет невязки, в блоке 6 - проверка того, что траектория снаряда найдена с заданной точностью.
Из физических соображений выбираются начальные условия z(0), удовлетворяющие системе уравнений (5).
Блок 1. Формирование начального приближения.
Интегрируется система дифференциальных уравнений (8) с достаточно малым шагом h по времени:
Одновременно интегрируется m систем дифференциальных уравнений:
На каждом шаге интегрирования вычисляется расстояние между снарядом и целью:
Если расстояние на текущем шаге увеличилось по сравнению с предыдущим, то фиксируется конечное время T(0)=t.
Блок 2. Расчет вектор-функции Ф и матрицы Якоби J.
В рассматриваемом блоке производится численное интегрирование систем дифференциальных уравнений на шаге с номером q итерационного процесса, а затем, по результатам интегрирования, формируется вектор-функция Ф(z(q), Т(q)) и ее матрица Якоби J(z(q), T(q)).
Начальные условия (3) системы дифференциальных уравнений (2) заданы и являются входными данными этого блока. Надо численно проинтегрировать на интервале [0; Т(q)] следующие системы дифференциальных уравнений:
Интегрирование можно производить любым известным численным методом, например методом Эйлера.
Блок 3. Формирование вектора поправки начальных условий траектории снаряда.
d=J-1Ф, где d - вектор поправки начальных условий, a J-1 - матрица, обратная матрице J.
Блок 4. Коррекция начальных условий расчетной траектории, проводимая в этом блоке, осуществляется по формуле:
Блок 5. Расчет невязки, осуществляемый по формуле:
Блок 6. Проверка точности.
Проверяется выполнение неравенства:
Δ<ε, где ε - некоторое малое наперед заданное число.
В предложенном способе стрельбы по цели поражающим снарядом определяют точку встречи и момент встречи цели и снаряда, для чего проводится линеаризация траектории снаряда по начальным условиям движения (углам наведения), определяется зависимость движения снаряда от его начальных условий. При этом матрица частных производных вектора положения снаряда по начальным условиям его движения удовлетворяет системе обыкновенных линейных неоднородных дифференциальных уравнений. Для решения возникающей краевой задачи применяется метод Ньютона. Все полученные результаты для системы дифференциальных уравнений описывают движение штатного снаряда.
Claims (1)
- Способ стрельбы по цели, сопровождаемой радиолокационной системой, двигающейся по произвольной траектории, описываемой уравнением y=g(t), y∈R3, поражающим снарядом, двигающимся по траектории, описываемой системой обыкновенных дифференциальных уравнений
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2015118108/28A RU2592747C1 (ru) | 2015-05-14 | 2015-05-14 | Способ стрельбы по движущейся цели |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2015118108/28A RU2592747C1 (ru) | 2015-05-14 | 2015-05-14 | Способ стрельбы по движущейся цели |
Publications (1)
Publication Number | Publication Date |
---|---|
RU2592747C1 true RU2592747C1 (ru) | 2016-07-27 |
Family
ID=56557016
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU2015118108/28A RU2592747C1 (ru) | 2015-05-14 | 2015-05-14 | Способ стрельбы по движущейся цели |
Country Status (1)
Country | Link |
---|---|
RU (1) | RU2592747C1 (ru) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU16211U1 (ru) * | 2000-06-09 | 2000-12-10 | Ростовский военный институт ракетных войск | Устройство адаптивного управления |
RU16959U1 (ru) * | 2000-08-21 | 2001-02-27 | Ростовский военный институт ракетных войск | Устройство адаптивного управления |
RU2170907C1 (ru) * | 2000-01-17 | 2001-07-20 | Мамошин Владимир Романович | Способ прицеливания при атаке скоростных целей истребителем по спрямленной траектории и устройство для его реализации |
RU2347999C2 (ru) * | 2007-03-20 | 2009-02-27 | Государственное унитарное предприятие "Конструкторское бюро приборостроения" | Способ стрельбы управляемым снарядом с лазерной полуактивной головкой самонаведения по движущейся цели (варианты) |
-
2015
- 2015-05-14 RU RU2015118108/28A patent/RU2592747C1/ru active IP Right Revival
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2170907C1 (ru) * | 2000-01-17 | 2001-07-20 | Мамошин Владимир Романович | Способ прицеливания при атаке скоростных целей истребителем по спрямленной траектории и устройство для его реализации |
RU16211U1 (ru) * | 2000-06-09 | 2000-12-10 | Ростовский военный институт ракетных войск | Устройство адаптивного управления |
RU16959U1 (ru) * | 2000-08-21 | 2001-02-27 | Ростовский военный институт ракетных войск | Устройство адаптивного управления |
RU2347999C2 (ru) * | 2007-03-20 | 2009-02-27 | Государственное унитарное предприятие "Конструкторское бюро приборостроения" | Способ стрельбы управляемым снарядом с лазерной полуактивной головкой самонаведения по движущейся цели (варианты) |
Non-Patent Citations (1)
Title |
---|
Сейдж Э.П., Уайт Ч.С. Оптимальное управление системами. -М.: Радио и связь, 1982, стр.81-91. Красовский Н.Н Оптимальное управление в обыкновенных динамических системах. УМН, 1965, том 20, выпуск 3(123), стр.153-174 . * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR102140097B1 (ko) | 총포 기반 대공 방어용 사격 통제 방법 | |
US6421116B1 (en) | Method for determining the relative movement between missile and target | |
Fonod et al. | Estimation enhancement by cooperatively imposing relative intercept angles | |
Fairfax et al. | Position estimation for projectiles using low-cost sensors and flight dynamics | |
US7400289B1 (en) | Plume-to-hardbody offset compensation in boosting missiles | |
Annabattula et al. | Advanced submarine integrated weapon control system | |
RU2592747C1 (ru) | Способ стрельбы по движущейся цели | |
Roh et al. | Trajectory optimization using Cramér-Rao lower bound for bearings-only target tracking | |
US11061107B1 (en) | System and method for intercepting an exo-atmospheric target using passive ranging estimation | |
KR102031929B1 (ko) | 연속적 시변 편향을 이용한 종말 선도각 제어 장치 및 방법 | |
Ullah et al. | Active vehicle protection using angle and time-to-go information from high-resolution infrared sensors | |
Raj et al. | Estimation of line-of-sight rate in a homing Missile Guidance loop using optimal filters | |
Ramezani et al. | Simulation‐Based Early Prediction of Rocket, Artillery, and Mortar Trajectories and Real‐Time Optimization for Counter‐RAM Systems | |
Zheng et al. | Terminal efficiency of fragment air-to-air missile using Monte Carlo method | |
Kale Güvenç | Range to-go estimation for a tactical missile with a passive seeker | |
Barbaric et al. | Passive ranging using image intensity and contrast measurements | |
JP4046626B2 (ja) | 飛翔体誘導装置 | |
Jahan et al. | Cannon fired ball with relative velocity | |
RU2499218C1 (ru) | Способ защиты объекта от средств воздушного нападения и система для его осуществления | |
Kim et al. | Launch point prediction employing the smoothing IPDA algorithm in 3-D cluttered environments | |
US11841211B2 (en) | Method for fire control of an anti-aircraft gun | |
AU785510B1 (en) | A method of a weapon guidance by target state estimation | |
KR102489644B1 (ko) | 30 mm 개틀링 함포의 실시간 사격 통제 명령 산출 장치 및 방법 | |
US11940249B2 (en) | Method, computer program and weapons system for calculating a bursting point of a projectile | |
RU2355990C2 (ru) | Способ определения точности наведения и сближения снаряда с целью по наблюдаемым параметрам их траекторного движения |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
MM4A | The patent is invalid due to non-payment of fees |
Effective date: 20170515 |
|
NF4A | Reinstatement of patent |
Effective date: 20190605 |