RU2539966C2 - Способ определения координат неподвижного радиомаяка с активным ответчиком - Google Patents

Abstract

Изобретение относятся к радиотехнике и может быть использовано, например, для определения местоположения неподвижного аварийного радиобуя, оборудованного радиомаяком с активным ответчиком. Достигаемый технический результат - повышение точности определения местоположения неподвижного радиомаяка с активным ответчиком и снижение количества измерителей дальности. Указанный результат достигается тем, что перемещающимся запросчиком последовательно излучается в N точках пространства сигнал запроса, траектория перемещения запросчика выбирается так, чтобы угол между направлениями от возможного местоположения радиомаяка на крайние точки траектории перемещения был больше, либо равен 180°, ответчик на радиомаяке принимает и переизлучает сигнал запроса, а координаты радиомаяка вычисляются по совокупности измеренных дальностей до радиомаяка в каждой из N точек измерения на траектории перемещения передатчика запросчика. 6 ил.

Description

Изобретение используется в радиотехнике.

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

1. Способ определения координат [1, 2], реализующий разностно-дальномерный метод. Такой способ заключаются в приеме и измерении разностей задержек сигнала источника радиоизлучения (ИРИ) в пунктах приема с известным местоположением, решении гиперболических уравнений на центральном приемном пункте и вычислении на этой основе координат ИРИ.

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

2. Способ определения координат разностно-доплеровским методом [3, с.246…248]. Сущность способа заключается в приеме сигнала ИРИ на двух перемещающихся измерителях, ретрансляции сигналов на центральный пункт обработки, вычислении разности доплеровских частот и разности радиальных скоростей измерителей с последующим вычислением координат ИРИ по разности доплеровских частот и разности радиальных скоростей.

Недостаток данного способа - необходимость использования не менее двух измерительных пунктов, кроме того, данные пункты должны перемещаться с достаточно большой скоростью, и данный способ не позволяет измерять координаты ИРИ, излучающих непрерывные или квазинепрерывные сигналы, т.к. в данном случае нельзя точно измерить разность частот сигнала [3, с.251].

Из известных способов определения координат неподвижного радиомаяка с активным ответчиком по технической сущности наиболее близким аналогом (прототипом) предлагаемого способа является способ, реализующий дальномерный метод [3, с.220…228], который включает посылку сигналов запроса передатчиками не менее чем с двух наземных пунктов, переизлучение сигналов запроса передатчиком ответчика на радиомаяке, координаты которого надо определить, измерение дальностей от наземных пунктов до радиомаяка R₁ и R₂ по величине задержки сигнала ответчика относительно сигналов запросчика и вычисление координат ответчика, как точки пересечения окружностей с радиусами R₁ и R₂, центры которых находятся в местах расположения запросчиков.

Недостатки способа прототипа:

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

2. При неизменной точности измерения дальностей для снижения ошибки определения координат радиомаяка с активным ответчиком необходимо увеличивать число стационарных приемо-передающих пунктов (запросчиков).

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

Поставленная цель достигается тем, что в известном способе определения координат неподвижного радиомаяка с активным ответчиком, включающем посылку сигнала запроса передатчиком запросчика, прием и переизлучение сигнала запроса ответчиком на радиомаяке, согласно изобретению излучение сигнала запроса осуществляют перемещающимся передатчиком запросчика в N точках с координатами [x_i,y_i], $$\text{ i}=\overline{\text{1,N}}$$

, измерение совокупности дальностей $${\tilde{R}}_i$$

, $$\text{i}=\overline{\text{1,N}}$$

по величине задержки сигнала ответчика относительно сигнала запросчика выполняют в каждой из N точек измерения на траектории перемещения запросчика, причем траектория перемещения запросчика выбирается так, чтобы угол между направлениями от возможного местоположения радиомаяка на крайние точки траектории перемещения был больше 180°, а вычисление координат [x,y] неподвижного радиомаяка с активным ответчиком осуществляется путем определения аргументов минимизации функции F(x,y):

$$\left[x,y\right]=\underset{\langle x,\text{ y}\rangle }{\arg \min }\left\{\text{F}\left(\text{x, y}\right)\right\}=\underset{\langle x,\text{ y}\rangle }{\arg \min }\left\{{\displaystyle \sum _{i=1}^N{\left(\sqrt{{\left(x_i-x\right)}^2+{\left(y_i-y\right)}^2}-{\tilde{R}}_i\right)}^2}\right\}$$

,

где $$F\left(x,\text{ y}\right)={\displaystyle \sum _{i=1}^N{\left(\sqrt{{\left(x_i-x\right)}^2+{\left(y_i-y\right)}^2}-{\tilde{R}}_i\right)}^2}$$

Сопоставительный анализ с прототипом показывает, что заявляемый способ отличается введением новой операции - перемещения запросчика по заданной траектории (так, чтобы угол между направлениями от возможного местоположения радиомаяка на крайние точки траектории перемещения был больше 180°). Таким образом, заявляемый способ соответствует критерию изобретения «новизна».

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

Поясним подробнее сущность заявляемого способа.

Предположим, что неподвижный радиомаяк с активным ответчиком находится в заданной области двухмерного пространства. Обозначим через [x,y] - местоположение радиомаяка, а координатами [x_i,y_i], $$\text{ i}=\overline{\text{1,N}}$$

- множество из N позиций на траектории перемещения запросчика, в которых измеряется дальность до радиомаяка. Аддитивные ошибки измерения дальностей будем считать независимыми нормальными случайными величинами с нулевым математическим ожиданием и одинаковыми дисперсиями $${\sigma }_R^2$$

. Тогда функция правдоподобия для независимых измерений представляет собой произведение условных гауссовских функций плотности вероятности:

где $${\tilde{R}}_i$$

, $$\text{ i}=\overline{\text{1,N}}$$

- вектор измеренных расстояний.

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

где $$F\left(x,\text{ y}\right)={\displaystyle \sum _{i=1}^N{\left(\sqrt{{\left(x_i-x\right)}^2+{\left(y_i-y\right)}^2}-{\tilde{R}}_i\right)}^2.}\left(\text{3}\right)$$

Определим на основании полученной функции правдоподобия (1) потенциальную точность оценки местоположения источника радиоизлучения дальномерным методом. Для этого воспользуемся выражением логарифма этой функции правдоподобия, применяя которое можно построить эталонную ковариационную матрицу, являющейся нижней границей дисперсии оценки координат [x,y]. Такая граница определяется неравенством Рао-Крамера [4, с.533…538]:

$$COV\left(x,\text{ y}\right)\ge {Ф}^{-1}=\left[\begin{array}{cc}{D}_x& K_{xy}\\ K_{xy}& D_y\end{array}\right],\text{ (4)}$$

где D_x - дисперсия оценки по оси X; K_xy - корреляционный момент оценки координат; D_y - дисперсия оценки по оси Y;

Ф - информационная матрица Фишера; E[•] - символ математического ожидания. Найдем значения частных производных в (5).

С этой целью представим логарифм функции правдоподобия (1) в следующем виде

где $$S_i\left[x,y\right]=\sqrt{{\left(x_i-x\right)}^2+{\left(y_i-y\right)}^2}.\text{ (7)}$$

Вычисление в общем виде производных в (5) с учетом (6) и (7) дает следующие результаты:

Так как по условию постановки задачи случайная величина $$\left(S_i\left[x,\text{ y}\right]-{\tilde{R}}_i\right)$$

распределена по нормальному закону с нулевым математическим ожиданием, то

$$\begin{array}{l}E\left[{\displaystyle \sum _{i=1}^N\left\{\left(S_i\left[x,\text{ y}\right]-{\tilde{R}}_i\right)\cdot \frac{{\partial }^2{S}_i\left[x,\text{ y}\right]}{\partial x^2}\right\}}\right]=\\ =E\left[{\displaystyle \sum _{i=1}^N\left\{\left(S_i\left[x,\text{ y}\right]-{\tilde{R}}_i\right)\cdot \frac{{\partial }^2{S}_i\left[x,\text{ y}\right]}{\partial y^2}\right\}}\right]=E\left[{\displaystyle \sum _{i=1}^N\left\{\left(S_i\left[x,\text{ y}\right]-{\tilde{R}}_i\right)\cdot \frac{{\partial }^2{S}_i\left[x,\text{ y}\right]}{\partial x\partial y}\right\}}\right]=0\text{. (9)}\end{array}$$

Тогда на основании (6)-(9) и с учетом того, что математическое ожидание функции, не имеющей случайных параметров, равно значению этой функции, получаем общие формулы для вычисления матрицы Фишера:

где

$$\{\begin{array}{l}\frac{\partial S_i\left[x,\text{ y}\right]}{\partial x}=\frac{x-x_i}{\sqrt{{\left(x_i-x\right)}^2+{\left(y_i-y\right)}^2}};\\ \frac{\partial S_i\left[x,\text{ y}\right]}{\partial y}=\frac{y-y_i}{\sqrt{{\left(x_i-x\right)}^2+{\left(y_i-y\right)}^2}}\end{array}.\text{ (11)}$$

Таким образом, матрица Фишера приобретает окончательный вид:

$$Ф=\left[\begin{array}{l}\frac{\text{1}}{{\sigma }_{\text{R}}^{\text{2}}}{\displaystyle \sum _{i=1}^N{\left(\frac{\partial S_i\left[x,\text{ y}\right]}{\partial x}\right)}^2}\frac{\text{1}}{{\sigma }_{\text{R}}^{\text{2}}}{\displaystyle \sum _{i=1}^N\frac{\partial S_i\left[x,\text{ y}\right]}{\partial x}\cdot \frac{\partial S_i\left[x,\text{ y}\right]}{\partial y}}\\ \frac{\text{1}}{{\sigma }_{\text{R}}^{\text{2}}}{\displaystyle \sum _{i=1}^N\frac{\partial S_i\left[x,\text{ y}\right]}{\partial y}\cdot \frac{\partial S_i\left[x,\text{ y}\right]}{\partial x}\frac{\text{1}}{{\sigma }_{\text{R}}^{\text{2}}}{\displaystyle \sum _{i=1}^N{\left(\frac{\partial S_i\left[x,\text{ y}\right]}{\partial y}\right)}^2}}\end{array}\right]=-\frac{\text{1}}{{\sigma }_{\text{R}}^{\text{2}}}\cdot \left[\begin{array}{cc}A& -B\\ -B& C\end{array}\right],\text{ (12)}$$

где

$$A={\displaystyle \sum _{i=1}^N\frac{{(x-x_i)}^2}{{\left(x_i-x\right)}^2+{\left(y_i-y\right)}^2}};\text{ (13)}$$

$$B={\displaystyle \sum _{i=1}^N\frac{(y-y_i)\cdot (y-y_i)}{{\left(x_i-x\right)}^2+{\left(y_i-y\right)}^2}};\text{ (14)}$$

$$C={\displaystyle \sum _{i=1}^N\frac{{(y-y_i)}^2}{{\left(x_i-x\right)}^2+{\left(y_i-y\right)}^2}};\text{ (15)}$$

Полученное выражение (12) для расчета матрицы Фишера свидетельствуют о том, что ковариационная матрица оценки координат [x,y], определяемая неравенством Рао-Крамера (4), не содержит в себе вектор измерения дальностей $${\tilde{R}}_i$$

, $$\text{ i}=\overline{\text{1,N}}$$

, не зависит от алгоритма поиска аргументов минимизации функции правдоподобия и определяется только конкретной совокупностью координат точек измерения дальностей [x_i,y_i], $$\text{ i}=\overline{\text{1,N}}$$

задающей вид траектории перемещения запросчика.

Поскольку расчетные оценки координат радиомаяка [x^*,y^∗], получаемые методом максимального правдоподобия в рамках сделанных выше допущений, асимптотически нормальны с математическим ожиданием [x,y] и ковариационной матрицей Ф^-1 [4, с.546…547], то функция плотности вероятности [x^*,y^∗] может быть описана двумерным гауссовским распределением [5, с.231]:

$$\begin{array}{l}W\left[{\text{x}}^{\text{*}}{\text{, y}}^{\ast }\right]=\frac{1}{2\pi \cdot {\sigma }_x\cdot {\sigma }_y\sqrt{1-r_{xy}^2}}\exp \cdot \\ \cdot \left\{-\frac{1}{2\sqrt{1-r_{xy}^2}}\left[\frac{{\left(x^{*}-x\right)}^2}{{\sigma }_x^2}+\frac{2\cdot r_{xy}\cdot \left(x^{*}-x\right)\cdot \left(y^{*}-y\right)}{{\sigma }_x\cdot {\sigma }_y}+\frac{{\left(y^{*}-y\right)}^2}{{\sigma }_y^2}\right]\right\},\text{ (16)}\end{array}$$

где $${\sigma }_x^2$$

, $${\sigma }_y^2$$

, r_xy - определяются из ковариационной матрицы (5): $${\sigma }_x^2=D_x$$

; $${\sigma }_y^2=D_y$$

; $$r_{xy}=K_{xy}/\sqrt{D_x\cdot D_y}$$

.

В геометрической интерпретации совместная двумерная плотность вероятности (16) представляет собой холмообразную поверхность, вершина которой находится над точкой [x,y] плоскости x0y. Сечение поверхности W[x^*,y^*] плоскостью, параллельной x0y, представляет собой эллипс рассеяния [5, с.235], уравнение проекции которого на плоскость x0y имеет вид:

$$\left[\frac{{\left(x^{*}-x\right)}^2}{{\sigma }_x^2}+\frac{2\cdot r_{xy}\cdot \left(x^{*}-x\right)\cdot \left(y^{*}-y\right)}{{\sigma }_x\cdot {\sigma }_y}+\frac{{\left(y^{*}-y\right)}^2}{{\sigma }_y^2}\right]=p_e^2\text{. (17)}$$

где p_e - параметр, определяющий вероятность попадания расчетной оценки координат в эллипс рассеяния (например, вероятность попадания в эллипс рассеяния с параметром p_e - 1,73 равна 0.95 [6, c.26]).

Размеры большой (a) и малой (b) полуосей эллипса рассеяния, а также угол (β) его ориентации (угол между большой осью эллипса и положительным направлением оси 0x, отсчитываемый против часовой стрелки) рассчитываются по формулам [7, с.580]:

$$a\left(b\right)=\sqrt{2}\cdot {\sigma }_R\cdot p_e/\sqrt{A+C-(+)\sqrt{{\left(A-C\right)}^2+4B^2}},\text{ (18)}$$

$$\beta =arctg\left[2B/\left(C-A\right)\right]/2,\text{ -}\pi \text{/2}<\beta \le \pi \text{/2 (19)}$$

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

Без потери общности рассуждений радиомаяк можно поместить в начало системы координат: x=0, y=0, а i-ю позицию измерения дальности представить относительно местоположения радиомаяка в полярных координатах [R_i, α_i], $$\text{ i}=\overline{\text{1,N}}$$

, где как показано на фиг.1а: α_i - угол между направлением от радиомаяка на i-ю позицию измерителя и осью 0x: α_i=arctg[y_i/x_i]; R_i - истинное расстояние от радиомаяка до i-й позиции измерения дальности. Тогда в предположении малых некоррелированных ошибок измерений дальностей (R_i>>σ_R, $$\text{ i}=\overline{\text{1,N}}$$

) вектор измеренных расстояний $${\tilde{R}}_i$$

, $$\text{ i}=\overline{\text{1,N}}$$

будет нормально распределенным, а выражения для расчета элементов матрицы Фишера (13)-(15) значительно упростятся:

$$A={\displaystyle \sum _{i=1}^N\frac{{\left(x-x_i\right)}^2}{{\left(x_i-x\right)}^2+{\left(y_i-y\right)}^2}}={\displaystyle \sum _{i=1}^N\frac{{\left(R_i\cdot \cos \left[{\alpha }_i\right]\right)}^2}{{\left(R_i\cdot \cos \left[{\alpha }_i\right]\right)}^2+{\left(R_i\cdot \sin \left[{\alpha }_i\right]\right)}^2}={\displaystyle \sum _{i=1}^N{\left(\cos \left[{\alpha }_i\right]\right)}^2}};\text{ (20)}$$

$$B=-{\displaystyle \sum _{i=1}^N\frac{\left(y-y_i\right)\cdot \left(x-x_i\right)}{{\left(x_i-x\right)}^2+{\left(y_i-y\right)}^2}=-{\displaystyle \sum _{i=1}^N\frac{\left(R_i\cdot \sin \left[{\alpha }_i\right]\right)\times \left(R_i\cdot \cos \left[{\alpha }_i\right]\right)}{{\left(R_i\cdot \cos \left[{\alpha }_i\right]\right)}^2+{\left(R_i\cdot \sin \left[{\alpha }_i\right]\right)}^2}}}=-{\displaystyle \sum _{i=1}^N\sin \left[{\alpha }_i\right]\times }\cos \left[{\alpha }_i\right];\text{ (21)}$$

$$C={\displaystyle \sum _{i=1}^N\frac{{\left(y-y_i\right)}^2}{{\left(x_i-x\right)}^2+{\left(y_i-y\right)}^2}}={\displaystyle \sum _{i=1}^N\frac{{\left(R_i\cdot \sin \left[{\alpha }_i\right]\right)}^2}{{\left(R_i\cdot \cos \left[{\alpha }_i\right]\right)}^2+{\left(R_i\cdot \sin \left[{\alpha }_i\right]\right)}^2}}={\displaystyle \sum _{i=1}^N{\left(R_i\cdot \sin \left[{\alpha }_i\right]\right)}^2.\text{ (22)}}$$

Так как по условиям постановки задачи радиомаяк неподвижен, позиции измерения на траектории движения запросчика могут быть перенумерованы произвольно, например, так, чтобы угол между направлением от радиомаяка на i-ю позицию запросчика α_i всегда увеличивался (фиг.1б).

Следовательно, с учетом (20)-(22) можно утверждать, что в рамках принятых ограничений и допущений (R_i>>σ_R, $$\text{ i}=\overline{\text{1,N}}$$

) любая траектория перемещения запросчика относительно радиомаяка может быть однозначно задана только возрастающей последовательностью углов α_i и углом визирования базы измерений¹ Θ (¹ Угол визирования базы измерений задается между направлениями от радиомаяка на крайние точки траектории движения эапросчика (на фиг.1б - это угол между направлениями наточку «1» и «8»)). В связи с этим задача обоснования траектории перемещения запросчика, обеспечивающей при фиксированном количестве позиций измерения - дальности минимальную ошибку оценки координат радиомаяка, сводится фактически к задаче обоснования угла визирования базы измерений, минимизирующего площадь эллипса ошибок: S=π·a·b.

Пусть угол между направлением от радиомаяка на i-ю позицию запросчика и осью 0x возрастает по линейному закону: $${\alpha }_i=i\times \frac{\Theta }{N}$$

. Тогда суммы в (20)-(22) могут быть представлены в замкнутой форме [7]:

$$A={{\displaystyle \sum _{i=1}^N\left(\cos \left[{\alpha }_i\right]\right)}}^2=\frac{N}{2}+\frac{1}{2}\cdot \sin \left[\Theta \right]\cdot \cos \left[\left(N+1\right)\cdot \frac{\Theta }{N}\right]/\sin \left[\frac{\Theta }{N}\right];\text{ (23)}$$

$$B=-{\displaystyle \sum _{i=1}^N\sin \left[{\alpha }_i\right]\times \cos \left[{\alpha }_i\right]=-\frac{1}{2}\cdot \sin \left[\Theta \right]\cdot \sin \left[\left(N+1\right)\cdot \frac{\Theta }{N}\right]/\sin \left[\frac{\Theta }{N}\right]};\text{ (24)}$$

$$C={\displaystyle \sum _{i=1}^N{\left(\sin \left[{\alpha }_i\right]\right)}^2=\frac{N}{2}}-\frac{1}{2}\cdot \sin \left[\Theta \right]\cdot \cos \left[\left(N+1\right)\cdot \frac{\Theta }{N}\right]/\sin \left[\frac{\Theta }{N}\right].\text{ (25)}$$

С учетом (23)-(25) выражения (18) для расчета большой и малой полуосей эллипса ошибок приобретают вид:

$$a\left(b\right)=\sqrt{2}\cdot {\sigma }_R/\sqrt{A+C-(+)\sqrt{{\left(A-C\right)}^2+4B^2}}=\sqrt{2}\cdot {\sigma }_R/\sqrt{N-(+)\sqrt{{\left(\sin \left[\Theta \right]/\sin \left[\frac{\Theta }{N}\right]\right)}^2}},\text{ (26)}$$

а площадь эллипса ошибок вычисляется по формуле:

$$S=\pi \cdot a\cdot b=\frac{2\cdot \pi \cdot {\sigma }_R^2\cdot p_e^2}{\sqrt{N^2-{\left(\sin \left[\Theta \right]/\sin \frac{\Theta }{N}\right)}^2}}.\text{ (27)}$$

Очевидно, что площадь эллипса ошибок будет минимальной при максимальном значении знаменателя в (27). Это достигается для $$\sin \left[\Theta \right]/\sin \frac{\Theta }{N}=0$$

, то есть когда:

$$\Theta =k\cdot \pi ,\text{ (28)}$$

где k - целое число, k=1, 2, 3….

Проиллюстрируем справедливость (28) следующим примером. Пусть вероятность попадания расчетной оценки координат радиомаяка в эллипс рассеяния задана равной 0.95, что соответствует параметру p_e=1.73 в выражении (17), а траектория движения запросчика представляет собой дугу с центром в точке нахождения ИРИ [x,y] и неизменным радиусом R.

В общем виде оценка координат радиомаяка как аргумент минимизации функции правдоподобия (1) может быть вычислена через определение глобального экстремума целевой функции (3) F(x,y) путем решения системы уравнений:

$$\{\begin{array}{l}\frac{\partial F\left(x,y\right)}{\partial x}=0,\\ \frac{\partial F\left(x,y\right)}{\partial y}=0.\end{array}\text{ (29)}$$

Вследствие нелинейности целевой функции (3), характерный вид которой иллюстрируется на фиг.2, поиск корней уравнений (29) может производиться численными методами на ЭВМ (например, с использованием функций пакета Matlab).

Результаты решения уравнений (29), полученные для некоторых частных случаев с использованием имитационного моделирования средствами Matlab R2012a, приведены на фиг.3-5.

Здесь радиомаяк располагается в центре осей координат[x=0, y=0], а результаты абсолютного большинства численных оценок координат радиомаяка группируются в виде точек вокруг них в пределах, очерчиваемых эллипсом ошибок. Bee расчеты проводились в предположении того, что среднеквадратическое отклонение (СКО) ошибки измерения дальности составляет σ_R=100 метров, а величина выборки измерений одинакова для всех траекторий движения и равна N=100. Графики подтверждают вывод о минимальности площади эллипса ошибок при угле визирования базы измерений, кратном 180° (на фиг.4 и 5 эллипсы вырождаются в окружности одинаковой площади).

На фиг.6 приведен построенный с использованием формулы (27) график нормированной площади эллипса ошибок $$S/\left({\sigma }_R^2\cdot p_e^2\right)$$

при N=100, который иллюстрирует зависимость точности оценки координат радиомаяка от угла визирования базы измерений.

Видно, что при фиксированном количестве пунктов измерения дальности N угол визирования базы измерений должен быть кратен 180°. Причем, как следует из графика на фиг.6, увеличение угла визирования базы измерений больше 180° не приводит к сколь-нибудь заметному изменению точности оценки координат радиомаяка. Поэтому для достижения максимальной точности определения координат радиомаяка целесообразно выбирать траекторию перемещения запросчика так, чтобы угол между направлениями от возможного местоположения радиомаяка на крайние точки траектории перемещения запросчика был заведомо больше 180°, но не превышал его значительно.

Найдем радиус максимального отклонения оценки координат радиомаяка при угле визирования 180°, то есть когда эллипс вырождается в окружность. С этой целью приравняем формулу площади эллипса ошибок (27) к выражению расчета площади окружности:

где r - радиус окружности, равный максимальному отклонению оценки координат при фиксированной вероятности, задаваемой параметром p_e.

Из (30) получаем формулу расчета радиуса максимального отклонения оценки координат радиомаяка для угла визирования 180° и вероятности попадания этой оценки внутрь окружности 0,99 (p_e=2,14):

$$r=\frac{\sqrt{2}\cdot {\sigma }_R\cdot p_e}{\sqrt{N}}=\frac{\sqrt{2}\cdot {\sigma }_R\cdot 2,14}{\sqrt{N}}\approx \frac{3\cdot {\sigma }_R}{\sqrt{N}}.\text{ (31)}$$

Для сравнения полученной оценки (31) с ошибкой местоопределения радиомаяка, достигаемой способом прототипом, воспользуемся приведенной в [3, с.221] формулой:

$$r_{П}=\frac{3\cdot {\sigma }_R}{\sqrt{2}}\cdot \cos ec\gamma \approx 2,121\cdot {\sigma }_R\cdot \cos ec\gamma \text{, (32)}$$

где γ - угол пересечения окружностей с радиусами R₁ и R₂, центры которых находятся в местах расположения стационарных запросчиков.

Примечание: В формуле (32), в отличие от приведенной в [3, с.221] формулы, дополнительно включен множитель 3 (по правилу «трех сигм») для получения радиуса попадания оценки внутрь окружности с вероятностью 0,99.

Поскольку для прототипа формула (32) позволяет проводить расчет точности оценки координат при двух измерениях дальностей, а в предлагаемом способе в общем случае реализуется N измерений, то корректнее предположить, что в прототипе также осуществляется N измерений дальностей (по N/2 измерений каждым запросчиком). Так как при многократных измерениях СКО ошибки измерения дальностей уменьшается пропорционально корню от числа измерений [5, с.459, ф.(11.7.4)], выражение (32) можно представить в виде:

$$r_{П}\approx \frac{2,121\cdot {\sigma }_R}{\sqrt{N/2}}\cdot \cos ec\gamma \approx \frac{\text{3}\cdot {\sigma }_{\text{R}}}{\sqrt{N}}\cdot \cos ec\gamma \text{. (33)}$$

Так как неравенство cosec γ≥1 всегда выполняется, из сравнения (31) и (33) можно утверждать, что при одинаковом числе измерений также всегда выполняется неравенство: r_П > r. Это обусловлено тем, что вероятность выполнения равенства cosec γ=1, достигаемого при γ=90°, для непрерывной случайной величины практически равна нулю.

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

Источники информации

1. Патент RU №2000129837, опубл. 20.10.2002 г.

2. Патент RU №2204145, опубл. 05.10.2003 г.

3. Кондратьев B.C. и др. Многопозиционные радиотехнические системы. - М.: Радио и связь, 1986. - 264 с.

4. Крамер Г. Математические методы статистики / Перевод с англ. под ред. А.Н. Колмогорова. - М.: МИР, 1975.

5. Вентцель Е.С., Овчаров Л.А. Теория вероятностей и ее инженерные приложения. - М.: Наука, 1988.

6. Сайбель А.Г. Основы теории точности радиотехнических методов местоопределения. - М.: ОБОРОНГИЗ, 1958.

7. Кукес И.С., Старик М.Е. Основы радиопеленгации. - М.: Сов. радио, 1964.

8. Прудников А.П., Брычков Ю.А., Маричев О.И. Интегралы и ряды. Элементарные функции. - М.: Наука, 1981.

Claims (1)

1. Способ определения координат неподвижного радиомаяка с активным ответчиком, включающий посылку сигнала запроса передатчиком запросчика, прием и переизлучение сигнала запроса ответчиком па радиомаяке, отличающийся тем, что излучение сигнала запроса осуществляют перемещающимся передатчиком запросчика в N точках с координатами [x_i, y_i],

   

   , измерение совокупности дальностей

   

   ,

   

   по величине задержки сигнала ответчика относительно сигнала запросчика выполняют в каждой из N точек измерения на траектории перемещения передатчика запросчика, причем траектория перемещения запросчика выбирается так, чтобы угол между направлениями от возможного местоположения радиомаяка на крайние точки траектории перемещения был больше 180°, а вычисление координат [x, y] неподвижного радиомаяка с активным ответчиком осуществляется путем определения аргументов минимизации (функции F(x, y):
   

   

   
   где

   

Images