СПОСОБ НАВИГАЦИИ ЛЕТАТЕЛЬНЫХ АППАРАТОВ Российский патент 2014 года по МПК G01S5/02 

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

Для обеспечения навигации летательных аппаратов (ЛА) определяют местоположение ЛА путем измерения их текущих координат. Местоположение ЛА определяют с помощью радиоволн в радиолокационных системах навигации.

Известен способ спутниковой радионавигации ГЛОНАСС, обеспечивающий непрерывность навигационно-временного обеспечения [1]. Это означает, что существует возможность определения координат местоположения, например, ЛА в любое время суток на базе спутниковых радионавигационных систем, включающих навигационные спутники (НС), наземный комплекс (НК), в составе наземной станции траекторных измерений и наземного вычислительного центра.

Реализация способа [1] заключается в следующем.

Используют несколько навигационных спутников (НС).

Формируют каждым НС сложный радиосигнал, который представляет собой гармоническое колебание с цифровой фазовой модуляцией.

Излучают сформированные НС сложные радиосигналы.

Принимают излученные сложные радиосигналы.

Измеряют дальность до ЛА с помощью принятых излученных сложных сигналов и доплеровское смещение частоты излученного сложного сигнала используя НК (наземную станцию траекторных измерений).

Определяют (вычисляют) три пространственные координаты ЛА, составляющие вектора скорости ЛА и текущее время с помощью НК (наземного вычислительного центра).

Передают вычисленную информацию с НК на борт НС.

Принимают переданную вычисленную информацию на борту НС.

Запоминают принятую информацию от НК на борту НС.

Передают на ЛА запомненную принятую информацию с борта НС в навигационном сообщении.

Принимают на ЛА навигационное сообщение с борта НС и используют полученную информацию о местоположении ЛА для управления движением ЛА (коррекции движения ЛА).

Недостатками способа [1] являются низкая помехоустойчивость работы при реализации способа при действии искусственных (преднамеренных) помех, которые делают способ [1] неработоспособным, и отсутствие автономности работы способа, поскольку необходим информационный контакт с несколькими НС. В результате ЛА также теряет свою работоспособность, поскольку на время действия помех информация о местоположении ЛА не поступает на борт ЛА, который движется по некорректируемой траектории.

Известен способ навигации летательных аппаратов [2], выбранный за прототип.

Способ [2] заключается в следующем.

Используют картосличительные корреляционно-экстремальные способы навигации (КЭСН), обеспечивающие отыскание и слежение за оптимальным режимом движения ЛА [3]. КЭСН основаны на сравнении текущих карт местности, полученных с использованием радиоволн, с эталонными картами той же местности, априори находящимися на ЛА, в основе которых лежит определение местоположения ЛА с последующим управлением движением ЛА путем коррекции их местоположения. Эталонные карты устанавливают на ЛА до момента начала движения над заданной поверхностью местности, а текущие карты получают во время движения ЛА. По отклонениям эталонных карт местности от текущих в заданной точке траектории движения ЛА определяют отклонение фактической траектории от заданной. В результате сравнения эталонных и текущих карт местности вырабатывают поправки в координаты ЛА для управления движением путем коррекции местоположения ЛА.

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

Способ [2] реализуется следующим образом.

Выбирают участок местности, размер которого в плановых координатах должен быть достаточен для определения местоположения ЛА.

Для мерного участка составляют эталонную карту местности в виде цифровой матрицы, размеры которой определяются размерами мерного участка.

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

Эталонная карта местности до начала движения находится на ЛА.

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

Составляют текущую карту путем измерения параметров мерного участка с помощью радиоволн. Лучи радиоволн излучают следующим образом. Первым излучают луч, направление распространения которого расположено в плоскости, ортогональной плоскости горизонта мерного участка, или в плоскости, которая находится под углом к плоскости горизонта мерного участка (первый луч). Далее излучают лучи, направления распространения которых не совпадают с направлением распространения первого луча, причем направления распространения одной части лучей находятся слева (по направлению движения ЛА) от первого луча, а другой части - справа (по направлению движения ЛА) от первого луча. Количество лучей справа и слева от первого луча одинаковое.

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

Вычисляют местные координаты проекции точки траектории ЛА на плоскость плановых координат.

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

Вычисляют для всех гипотез внутри квадрата неопределенности слагаемые показателя близости.

Проводят по завершении всех измерений поиск экстремума показателя близости.

Определяют поправки к координатам местоположения ЛА в плановых координатах мерного участка на основе анализа взаимных смещений эталонной и текущей карт местности мерного участка.

После определения местоположения ЛА в плановых координатах мерного участка вычисляют высоту ЛА над поверхностью мерного участка в координатах мерного участка (в точке определения местоположения ЛА в плановых координатах мерного участка).

Выдают поправки к координатам местоположения ЛА в плановых координатах мерного участка по трем координатам.

Управляют движением ЛА путем коррекции его местоположения по трем координатам по мере прохождения мерного участка.

Недостатком способа [2] является то, что при составлении текущей карты местности отсутствует информация о направлении движения летательного аппарата, скорости и ускорении движения ЛА.

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

Технический результат достигается тем, что способ навигации летательных аппаратов заключается в использовании эталонной карты местности, выборе участка местности (мерный участок), находящегося в пределах эталонной карты, составлении текущей карты мерного участка путем измерений наклонных дальностей с помощью многолучевого режима измерения при помощи радиоволн, излучаемых в виде лучей, определении разности результатов многолучевых измерений по текущей карте, сравнении текущей и эталонных карт в пределах квадрата неопределенности для определения показателя близости путем сравнения разности результатов многолучевых измерений текущей и эталонной карт, определении плановых координат и высоты летательного аппарата, вычислении сигнала коррекции траектории движения и управлении движением летательного аппарата по трем координатам (плановые координаты и высота) путем коррекции его местоположения. Перед вычислением сигнала коррекции через равные интервалы времени Δt после получения первой текущей карты мерного участка получают вторую текущую карту и третью текущую карту мерного участка, определяют координаты (плановые координаты и высота) второго и третьего местоположений летательного аппарата в плановых координатах эталонной карты путем последовательного сравнения первой текущей и эталонной карт мерного участка в пределах первого квадрата неопределенности, представляющего собой массив данных M₁, для определения показателя близости, второй текущей и эталонной карт мерного участка в пределах второго квадрата неопределенности, представляющего собой массив данных M₂, для определения показателя близости, третьей текущей и эталонной карт мерного участка в пределах третьего квадрата неопределенности, представляющего собой массив данных M₃, для определения показателя близости. Размеры второго и третьего квадратов неопределенности значительно меньше размеров первого квадрата неопределенности M₂<<M₁ и M₃<<M₁. Определяют направление движения, скорость и ускорение летательного аппарата путем сравнения координат (плановые координаты и высота) первого, второго и третьего местоположений летательного аппарата в плановых координатах эталонной карты, вычисляют сигнал коррекции траектории движения. Управляют движением летательного аппарата путем коррекции его местоположения используя полученные значения координат (плановые координаты и высота) и трехмерного вектора движения летательного аппарата (направление, скорость и ускорение движения по трем координатам).

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

В результате способ навигации позволяет повысить точность определения местоположения ЛА за счет использования информации о параметрах траектории движения ЛА - о трехмерном векторе движения ЛА (направлении, скорости и ускорении движения).

Способ навигации ЛА поясняют следующие чертежи:

- на фигуре 1 представлено взаимное расположение систем координат при движении ЛА;

- на фигуре 2 показано измерение наклонной дальности по отдельным лучам радиоволн;

- на фигуре 3 показан вектор движения ЛА, причем координаты ЛА (x₁, y₁, z₁), где (x₁, y₁) - плановые координаты и z₁ - высота, получены с помощью первой текущей карты, координаты ЛА (x₂, y₂, z₂) получены с помощью второй текущей карты, а координаты ЛА (x₃, y₃, z₃) получены с помощью третьей текущей карты, стрелка - вектор движения ЛА, указывающий направление движения, $$\overrightarrow{V}$$ - вектор скорости движения ЛА, представленный в виде векторов скоростей $${\overrightarrow{V}}_1$$ и $${\overrightarrow{V}}_2$$ , S - расстояние между местоположениями ЛА, представленное в виде расстояний S₁ и S₂.

Способ навигации ЛА осуществляется следующим образом.

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

При анализе будем применять следующие системы координат, показанные на фиг.1, для:

- эталонной карты рельефа местности (большой прямоугольник на фиг.1);

- квадрата неопределенности (квадрат на фиг.1);

- возможных положений ЛА внутри квадрата неопределенности в момент начала измерений (точки на фиг.1) и точек траектории ЛА, в которых производят измерения (ромбы на фиг.1).

Местная система координат является левой прямоугольной декартовой системой координат O_rx_ry_rz_r с началом O_r. При этом оси O_rx_r и O_ry_r лежат в плоскости местного горизонта, то есть x_r и y_r являются плановыми координатами ЛА. Относительно плоскости O_rx_ry_r определяют высоты рельефа местности и ЛА. Таким образом, указанные высоты являются соответствующими значениями координаты z_r. Начало координат O_r выбрано так, что ось O_rx_r направлена в расчетную точку появления ЛА над участком местности, соответствующим эталонной карте. Ось O_rx_r считается коллинеарной по отношению к горизонтальной составляющей расчетного вектора скорости ЛА (фиг.1). Плановая система координат неподвижна и связана с эталонной картой.

Для того чтобы привязать плановые координаты к эталонной карте, используют дискретную систему координат. При этом оси N_x и N_y сонаправлены с осями местной системы O_rx_ry_r (фиг.1). Нулевым индексам в дискретной системе координат соответствует левый нижний угол эталонной карты. Предполагаемое значение дискретных координат ЛА в момент начала сбора данных O_r обозначим как (n_xrЭ, n_yrЭ).

Связь между координатами дискретной и местной систем координат определяется выражениями вида

$$\begin{array}{c}{n}_{xr}=n_{xЭ}+x_r/{\Delta }_{Э},\\ n_{yr}=n_{yЭ}+y_r/{\Delta }_{Э},\end{array}(1)$$

где Δ_Э - размер ячейки сетки эталонной карты.

При определении координат точек измерений в виде точек пересечения подстилающей поверхности и лучей радиоволн используют систему координат, связанную с текущим положением ЛА. Начало координат этой системы (точка O_a) поместим в текущую точку проекции траектории ЛА на плоскость плановых координат. Ось O_ax_a по направлению совпадает с текущим направлением вектора скорости ЛА.

Связь между координатами системы, связанной с ЛА, и местной системы координат определяется выражениями

$$\begin{array}{c}{x}_r=x_a\cos {\alpha }_x+y_a\sin {\alpha }_x+x_{rc},\\ y_r=y_a\cos {\alpha }_x-x_a\sin {\alpha }_x+y_{rc},\end{array}(2)$$

где x_rc и y_rc - координаты проекции текущего положения ЛА на плоскость плановых координат; α_x - курс ЛА в текущей точке измерений (угол между вектором скорости ЛА и осью O_rx_r, на фиг.1 - это угол между осями x_r и x_a).

Исходными данными для вычислений в КЭСН являются:

- эталонная карта, представляющая собой массив данных о рельефе местности, элементы которого представляют собой высоты рельефа в узлах координатной сетки эталонной карты на плоскости O_rx_ry_r (фиг.1);

- данные, определяемые углом между лучами α_R радиоволн;

- данные измерений, отличных от полученных с помощью лучей радиоволн, поступающие с каждым результатом измерений: значения угла α_z крена, угла α_x курса, угла α_t тангажа, скорости τ ЛА (фиг.2);

- текущая карта, представляющая собой массив данных М в виде набора значений наклонных дальностей по всем трем лучам, полученных в каждом измерении.

Текущая карта местности определяется матрицей H(n_x, n_y) размером N_x×N_y, а эталонная карта - матрицей H_Э(n_x, n_y) размером N_xЭ×N_yЭ. Дискретные значения n_sx и n_sy соответствуют горизонтальному и вертикальному смещению текущей карты относительно эталонной и отсчитываются от левого нижнего угла эталонной карты, для которого n_sx=n_sy ~ ® •

В предлагаемом способе используют дифференциально-разностный алгоритм обработки многолучевых измерений, который рассмотрим на примере трехлучевого КЭСН. Лучи радиоволн излучают последовательно в следующем порядке: первым - центральный луч (по направлению движения ЛА), вторым - левый луч и третьим - правый луч относительно центрального луча.

Используемый алгоритм состоит в определении разности измерений левого и центрального лучей, правого и центрального лучей текущего измерения, а также в вычислении разности измерений центрального луча в текущем измерении и в предыдущем.

Показатель близости данных получим в виде

$$\begin{array}{l}P\left(n_{sx},n_{sy}\right)={\displaystyle \sum _{n_x=1}^{N_{xR}{}^{-1}}[\left|H^{12}\left(n_x\right)-H_{Э}^{12}\left(nx,nsx,nsy\right)\right|+}\\ +\left|H^{32}\left(n_x\right)-H_{Э}^{32}\left(n_x,n_{sx},n_{sy}\right)\right|+\\ +|\left[H^{22}\left(n_x\right)-H_{Э}^{22}\left(n_x,n_{sx},n_{sy}\right)\right]-\\ -\left[H^{22}\left(n_x-1\right)-H_{Э}^{22}\left(n_x,n_{sx},n_{sy}\right)\right]|,(3)\end{array}$$

где для обозначения разностей измерений применена следующая индексация: верхний индекс 12 означает вычитание результата измерения по первому лучу из результатов измерения по левому, а верхний индекс 32 - вычитание результата по первому лучу из результата по правому, верхний индекс 22 означает результат измерения значения высоты по первому лучу в данном и в предыдущем измерении.

Обработку данных каждого измерения производят следующим образом.

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

Вычисляют местные координаты проекции точки траектории ЛА на плоскость плановых координат.

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

Проводят по завершении всех измерений поиск экстремума показателя близости.

Определяют поправки к координатам местоположения ЛА в плановых координатах мерного участка на основе анализа взаимных смещений эталонной и текущей карт местности мерного участка.

Вычисляют местоположение ЛА в плановых координатах мерного участка.

Вычисляют высоту ЛА в координатах мерного участка.

Получают последовательно через равные промежутки времени Δt вторую текущую карту мерного участка и третью текущую карту мерного участка.

Используют при определении второго и третьего местоположений ЛА в плановых координатах эталонной карты второй и третий квадраты неопределенности, представляющие собой массивы данных M₂ и M₃ соответственно, размеры которых значительно меньше размеров первого квадрата неопределенности, представляющего собой массив данных M_x, M₂<<M₁ и M₃<<M₃.

Вычисляют последовательно второе и третье местоположения ЛА в плановых координатах мерного участка на основе измерений с помощью лучей радиоволн наклонных дальностей от ЛА до поверхности мерного участка используя последовательно вторую и третью текущие карты.

Получают значения вектора движения летательного аппарата: координаты трех местоположений ЛА (фиг.3), а также скорость и ускорение движения ЛА в трехмерном пространстве.

Выдают поправки к координатам местоположения ЛА (по трем координатам: плановым координатам и по высотам) используя значения трехмерного вектора движения ЛА.

Рассмотрим данный алгоритм подробно.

По измеренным значениям наклонных дальностей до элементов подстилающей поверхности D_j (j∈[1, 3] - номер луча радиоволн) определяют координаты точек измерений в системе координат, связанной с ЛА, а также измеряемые значения высоты рельефа H_j.

Согласно фиг.2 получим систему уравнений для каждого луча j в следующем виде

$$\begin{array}{l}{x}_j=D_j\cos \left({\alpha }_z^j-{\alpha }_R^j\right)\cdot \cos \left({\alpha }_t^j\right),\\ y_j=D_j\sin \left({\alpha }_z^j-{\alpha }_R^j\right),\\ H_j=D_j\cos \left({\alpha }_z^j-{\alpha }_R^j\right)\cdot \sin \left({\alpha }_t^j\right).(4)\end{array}$$

По полученным данным о скорости и углах эволюции определяют координаты проекции точки траектории ЛА на плоскость плановых координат.

Для первого измерения эти координаты считаются нулевыми (то есть положение ЛА во время первого измерения совпадает с началом координат системы O_rx_ry_r - фиг.1)

$$x_{rc}^1=\mathrm{0,}{y}_{rc}^1=0.(5)$$

Местные координаты ЛА во время выполнения следующих измерений вычисляют по рекуррентным соотношениям

$$\begin{array}{l}{x}_{rc}^i=x_{rc}^{i-1}+{\nu }^{i-1}{T}_c\cos {\alpha }_x^{i-1},\\ y_{rc}^i=y_{rc}^{i-1}+{\nu }^{i-1}{T}_c\sin {\alpha }_x^{i-1},(6)\end{array}$$

где T_c - период поступления измеренных данных, i=1, 2,… - номер текущего измерения.

Зная координаты точек измерений в системе координат, связанной с ЛА (O_ax_ay_a), и местные координаты начала O_a, определяют координаты точек измерений в местной системе координат

$$\begin{array}{l}{x}_r^i=x_j\cos {\alpha }_x^j+y_j\sin {\alpha }_x^j+x_{rc},\\ y_r^i=y_j\cos {\alpha }_x^j+y_j\sin {\alpha }_x^j+y_{rc}.(7)\end{array}$$

Перебор гипотез заключается в следующем. Для каждого из возможных положений ЛА в пределах квадрата неопределенности вычисляют координаты (в дискретной системе координат) точек измерений, определенных ранее в местной системе координат. Для этого точку O_r помещают поочередно в узлы координатной сетки внутри квадрата неопределенности и выполняют вычисление показателя близости данных. Координаты точек измерений в дискретной системе координат определяют как

$$\begin{array}{l}{n}_x=x_r/{\Delta }_{Э}+n_{xrЭ}+n_{hx},\\ n_y=y_r/{\Delta }_{Э}+n_{yrЭ}+n_{hy},(8)\end{array}$$

где n_hx и n_hy - смещения точки O_r для различных гипотез, определенные в количестве узлов сетки эталонной карты.

Показатель близости данных P(n_hx, n_hy) для данного случая имеет вид

$$\begin{array}{l}P\left(n_{hx},n_{hy}\right)={\displaystyle \sum _{k+1}^K[\left|H_k^{12}-H_{Э,k}^{12}\left(n_{hx},n_{hy}\right)\right|+\left|H_k^{32}-H_{Э,k}^{32}\left(n_{hx},n_{hy}\right)\right|+}\\ +\left|H_k^{22}-H_{Э,k}^{22}\left(n_{hx},n_{hy}\right)\right|].(9)\end{array}$$

Здесь n_hx и n_hy - смещения точки O_r для различных гипотез; К - количество измерений; $$H_k^{12}$$ - разность измеренных в k-м измерении значений высоты по первому и левому лучам; $$H_k^{32}$$ - разность измеренных в k-м измерении значений высоты по первому и правому лучам; $$H_k^{22}$$ - разность значений высоты измеренных в k-м и (k-1)-м измерениях; $$H_{Э,k}^{12}(\cdot )$$ - разность определенных для некоторой гипотезы (для определенного значения n_hx и n_hy) значений высоты по данным эталонной карты на k-м измерении по первому и левому лучам; $$H_{Э,k}^{32}(\cdot )$$ - разность определенных для некоторой гипотезы значений высоты по данным трехмерной эталонной карты на k-м измерении по первому и правому лучам; $$H_{Э,k}^{22}(\cdot )$$ - разность определенных для некоторой гипотезы значений высоты по данным эталонной карты на k-м и (k-1)-м измерениях по первому лучу.

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

По завершении всех измерений производят поиск глобального экстремума показателя близости данных P(n_hx, n_hy) и определяют поправки к координатам местоположения ЛА по плановым координатам.

Определяют поправки к координатам местоположения ЛА по высоте в момент определения поправок местоположения ЛА в плановых координатах мерного участка.

Описанную выше последовательность операций применяют к каждой из трех текущих карт мерного участка.

Для этого получают последовательно через равные интервалы времени Δt первую, вторую и третью текущие карты мерного участка.

Размеры первого квадрата неопределенности при определении первого местоположения ЛА в плановых координатах эталонной карты определяются величиной допустимых отклонений фактической траектории ЛА от заданной.

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

При этом не происходит снижение точности измерений, поскольку при переборе гипотез уменьшилось количество переборов за счет исключения гипотез, не влияющих на определение местоположения ЛА, а время измерения существенно сокращается. Поэтому обработка информации по второй и третьей текущим картам увеличивает время обработки информации по текущим картам не в 3 раза, а значительно меньше - не более чем на 20%.

Определяют координаты первого (x₁, y₁, z₁), второго (x₂, y₂, z₂) и третьего (x₃, y₃, z₃) местоположений ЛА в плановых координатах эталонной карты, сравнивая последовательно первую текущую и эталонную карты мерного участка, вторую текущую и эталонную карты мерного участка и третью текущую и эталонную карты мерного участка.

Определяют по координатам трех местоположений ЛА (с использованием первой, второй и третьей текущих картах мерного участка) вектор направления движения ЛА (фиг.3) - скорость движения, направление движения и ускорение ЛА.

Определяют направление и скорость V₁ движения ЛА после нахождения второго местоположения ЛА в плановых координатах мерного участка эталонной карты (фиг.3) используя вторую текущую карту мерного участка.

Скорость V₁ движения ЛА определим в виде

$$V_1=\Delta S_1/\Delta t_1,(10)$$

где ΔS₁ - расстояние между точкой с координатами (x₁, y₁, z₁) на эталонной карте (первое местоположение ЛА в плановых координатах мерного участка эталонной карты) и точкой с координатами (x₂, y₂, z₂) на эталонной карте (второе местоположение ЛА в плановых координатах мерного участка эталонной карты) (фиг.3); Δt₁=Δt - временной интервал между моментом времени t₁, в который была получена первая текущая карта, и моментом времени t₂=t₁+Δt, в который была получена вторая текущая карта.

Направление движения ЛА определим как направление вектора скорости $$\overrightarrow{V}={\overrightarrow{V}}_1$$ из точки с координатами (x₁, y₁, z₁) в точку с координатами (x₂, y₂, z₂) в плановых координатах мерного участка эталонной карты (фиг.3).

Уточняют направление и скорость движения ЛА после нахождения третьего местоположения ЛА в плановых координатах мерного участка эталонной карты и определяют ускорение движения ЛА (фиг.3).

Скорость V₂ движения ЛА определим в виде

$$V_2=\Delta S_2/\Delta t_2,(11)$$

где ΔS₂ - расстояние между точкой с координатами (x₂, y₂, z₂) на эталонной карте (второе местоположение ЛА в плановых координатах мерного участка эталонной карты) и точкой с координатами (x₂, y₂, z₂) на эталонной карте (третье местоположение ЛА в плановых координатах мерного участка эталонной карты) (фиг.3); Δt₂=Δt - временной интервал между моментом времени t₂, в который была получена вторая текущая карта, и моментом времени t₃=t₂+Δt, в который была получена третья текущая карта.

Направление движения ЛА определим как направление вектора скорости $$\overrightarrow{V}={\overrightarrow{V}}_2$$ из точки с координатами (x₂, y₂, z₂) в точку с координатами (x₃, y₃, z₃) в плановых координатах мерного участка эталонной карты (фиг.3).

Уточненную (среднюю) V_cp скорость движения ЛА определим как

$$V_{cp}=\left(V_1+V_2\right)/2.(12)$$

Уточненное направление движения ЛА определим как направление вектора скорости V_cp из точки с координатами (x₁, y₁, z₁) в точку с координатами (x₃, y₃, z₃) в плановых координатах мерного участка эталонной карты.

Определяют ускорение ЛА.

Ускорение W движения ЛА определим как

$$W=\Delta V/\Delta t,(13)$$

где ΔV=V₂-V₁ - изменение скорости движения ЛА за время Δt.

Вычисляют сигнал коррекции траектории движения и управляют движением ЛА путем коррекции его местоположения дополнительно используя значения вектора движения ЛА (направление движения, скорость и ускорение движения в трехмерном пространстве).

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

Реализуется повышенная точность управления движением ЛА, поскольку определяют поправки по трем координатам, скорости и ускорению движения и выдают сигналы коррекции местоположения ЛА, соответственно, по трем координатам, скорости и ускорению движения, чего не было в прототипе.

Обработка информации по второй и третьей текущим картам увеличивает время обработки информации по текущим картам не в 3 раза, а значительно меньше - не более чем на 20%.

Важно отметить, что рассмотренный способ навигации ЛА с предложенным алгоритмом сохраняет свои положительные свойства при различном количестве лучей радиоволн и обеспечивает работоспособность и для трех, и для большего количества лучей радиоволн. Число используемых лучей радиоволн определяется только временем, в течение которого обеспечивается измерение местоположения ЛА при движении над мерным участком поверхности.

ЛИТЕРАТУРА

1. ГЛОНАСС. Принципы построения и функционирования / Под ред. А.И.Перова, В.Н.Харисова. - М.: Радиотехника, 2005. - С.688 (с.12-15).

2. Патент РФ №2338158, МПК G01C 21/00 (2006.01). Способ навигации летательных аппаратов / Хрусталев А.А., Кольцов Ю.В., Егоров С.Н. // Изобретения. Полезные модели. - 2008. - Опубл. 10.11.2008. - Бюл. №31 (прототип).

3. Белоглазов И.Н., Джанджгава Г.И., Чигин Г.П. Основы навигации по геофизическим полям. - М.: Наука, 1985. - С.328 (с.10-11, 19-22, 25-34)

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

Способ навигации летательных аппаратов, заключающийся в использовании эталонной карты местности, выборе участка местности (мерный участок), находящегося в пределах эталонной карты, составлении текущей карты мерного участка путем измерений наклонных дальностей с помощью многолучевого режима измерения при помощи радиоволн, излучаемых в виде лучей, определении разности результатов многолучевых измерений по текущей карте, сравнении текущей и эталонных карт в пределах квадрата неопределенности для определения показателя близости путем сравнения разности результатов многолучевых измерений текущей и эталонной карт, определении плановых координат и высоты летательного аппарата, вычислении сигнала коррекции траектории движения и управлении движением летательного аппарата по трем координатам (плановые координаты и высота) путем коррекции его местоположения, отличающийся тем, что перед вычислением сигнала коррекции через равные интервалы времени Δt после получения первой текущей карты мерного участка, получают вторую текущую карту и третью текущую карту мерного участка, определяют координаты (плановые координаты и высота) второго и третьего местоположений летательного аппарата в плановых координатах эталонной карты путем последовательного сравнения первой текущей и эталонной карт мерного участка в пределах первого квадрата неопределенности, представляющего собой массив данных M₁, для определения показателя близости, второй текущей и эталонной карт мерного участка в пределах второго квадрата неопределенности, представляющего собой массив данных M₂, для определения показателя близости, третьей текущей и эталонной карт мерного участка в пределах третьего квадрата неопределенности, представляющего собой массив данных M₃, для определения показателя близости, причем размеры второго и третьего квадратов неопределенности значительно меньше размеров первого квадрата неопределенности M₂<<M₁ и M₃<<M₁, определяют направление движения, скорость и ускорение летательного аппарата путем сравнения координат (плановые координаты и высота) первого, второго и третьего местоположений летательного аппарата в плановых координатах эталонной карты, вычисляют сигнал коррекции траектории движения и управляют движением летательного аппарата путем коррекции его местоположения, используя полученные значения координат (плановые координаты и высота) и трехмерного вектора движения летательного аппарата (направление, скорость и ускорение движения по трем координатам).