Файл: Олянюк, П. В. Оптимальный прием сигналов и оценка потенциальной точности космических измерительных комплексов.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 23.10.2024

Просмотров: 77

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

1.Основная форма представления матриц перехода

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

делении

 

составляющих

исходной

g

и

конечной

систем

параметров q t

в общем

случае

может быть представлено

следующей

зависимостью:

 

 

PT,Ug P.,

 

PTjUgP.2

 

 

 

 

PT,UB P . m

 

 

 

 

 

 

P U L P . ,

P U L P . 2

 

 

 

 

PfolJg2 u e P./H

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

(VI.4.1)

 

 

P T „,U f f P . ,

 

P T m U e P . 2

 

 

 

 

P T m U g P . m

 

 

 

 

которая

показывает,

что при не нулевой матрице

U g

 

матрица U q

будет

особенной

(clet (J q =0)

тогда,

когда

один

или

несколько

 

столбцов Р.у

будут нулевыми

или между

ними

будет

наблюдаться

пропорциональ­

ность.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так,

 

например,

если у'-й столбец

 

матрицы

Р

представляет

собой

нуль-вектор, то для

всех

п и

k (п, & = 1 , 2, . . . , ш)

имеем

тождество

 

 

 

 

Р Т , . и г Р . п ^ Р ^ и е Р . ,

 

о,

 

 

 

 

 

(VI.4.2)

т. е. у'-я строка

и у'-й столбец

матрицы

U q

вырождаются в

нулевые. Если

между у'-м и t-м столбцами

матрицы

Р

наблюдается

пропорциональ­

ность (пусть P.J =

аР . ; ) , то для

всех

элементов

матрицы

U q j содержа­

щих указанные

столбцы, будет справедливо

соотношение

 

 

 

 

 

 

PT . UgP . y/P^P . , - ^

Р т у и г Р . „ / Р ! г и г Р . „ s

a.

 

(VI.4.3)

Это

свидетельствует

о том, что между

у'-м и £-м столбцами

(строка­

ми) матрицы \J q имеется

пропорциональность,

т. е. матрица

U q

 

особен­

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

2. Критерии анализа матриц перехода

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

дания

кеплеровых параметров

орбиты,

для которых определитель

мат­

рицы

Р имеет минимальное

значение

или равен нулю. Это будут

об-

11"

 

 

 

163-


ласти пространства кеплеровых элементов орбиты, в которых матрица

Рблизка к сингулярной (особенной) или сингулярна. Прежде всего по­

пытаемся

найти

те области

задания

 

кеплеровых

элементов

орбиты, для

которых

определитель

матрицы

стремится

к нулю. При этом

будем по­

лагать,

что

определитель

матрицы

равен

нулю

только

в тех случаях,

когда

столбцы

матрицы

Р

вырождаются

в нуль-векторы

 

или

между

ни лги

наблюдается

пропорциональность

(линейная

зависимость).

 

Исследуем

матрицу

перехода,

записанную

в

форме

(VI.2.23), т.' е.

исследуем

столбцы

Р, 7

матрицы

перехода

на

вырождение

их в

нуле­

вые векторы и на пропорциональность между ними.

 

 

 

 

 

 

Из

предыдущего материала

видно, что ни одна

из матриц,

входящих

в качестве

 

сомножителей

 

в матрицы-столбцы

Р.,-,

не вырождается в

нуль-матрицы

ни при каких

реальных

значениях

кеплеровых

 

параметров

орбиты и текущего аргумента Е. Поэтому

элементы Р.у

матрицы

пере­

хода

могут

принимать

нулевые

значения только

в

тех случаях,

когда

линейная комбинация элементов матриц, зависящих от аргумента Е (эле­

ментов

матрицы g ,

и

производных

от ее

составляющих

по

парамет­

рам а, е, М0),

и элементов

других матриц-сомножителей

(элементов

мат­

риц G,

Н, S, и их производных

соответственно по углу

Q,

i , м) будет

равна

нулю.

 

 

 

 

 

 

 

 

 

 

 

 

В

качестве критерия вырождения любого из векторов

Р . ,

в

нуль-

вектор

можно

принять

критерий

равенства

нулю

модуля

вектора

Р.

или его квадрата. В

самом

деле,

Р.у

может

быть

нуль-вектором

тогда ;

ц только тогда, когда

имеется равенство

 

 

 

 

 

 

 

 

 

 

 

 

Р у . Р . у ^ О .

 

 

 

 

(VI.4.4)

 

В

качестве

критерия

линейной

зависимости

(пропорциональности)

двух

векторов Р.у и Р.£

можно принять критерий

равенства

нулю мо­

дуля разности двух векторов, один из которых есть

исследуемый

вектор,

например р ^

а

второй — произведение некоторой вещественной

скаляр­

ной

величины

а

на

второй исследуемый

вектор Р.,-. Если

составить

произведение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( Р . / -

аР.,-)т ( Р . / -

<*Р./) = е,

 

 

(VI.4.5)

то можно показать,

что это выражение

будет равно

нулю только в том

случае,

если между

составляющими

векторов

Р.( - и Р . ; имеется

пропор­

циональность, а а есть коэффициент пропорциональности. Следовательно, задача исследования на пропорциональность двух векторов есть задача

отыскания

такого

вещественного значения коэффициента

а,

не равного

нулю, при котором

с учетом области задания параметров

q

будет вы :

полниться

равенство

 

 

 

( Р . / -

аР.,-)т

( Р . ; - а Р . , ) =

0.

(VI.4.6)

"Если это равенство

выполняется

при а = 0,

то

фактически

выполняется

условие

(VI.4.4), т. е. Р.у

есть

нуль-вектор.

 

 

 

Можно показать,

что если

Р.^ и Р.,-

ортогональны, то

выражение

((VI.4.5)

приводится к соотношению

 

 

 

 

 

P:.P./ + a2P:,P., = p

 

(VI.4.7)

1G4


н может быть равно нулю в

той

области задания

параметров q ; в ко­

торой оба

вектора .нулевые,

так

как

ортогональные

векторы линейно

независимы.

 

 

 

 

 

 

Необходимо отметить,

что

условие (VI.4.6) выполняется тогда, когда

неравенство

Буняковского

Коши,

используемое в работе [4] для ха­

рактеристики общих свойств навигационных методов, обусловленных осо­ бенностями фундаментальных матриц с усредненными элементами, пре­

вращается

в

равенство.

 

 

 

 

 

На основе проведенных исследований по определению и анализу

квадратов

модулей

векторов

матрицы

Р

(VI.2.23) можно

прийти к

сле­

дующим выводам.

 

 

 

 

 

 

Квадрат

модуля

первого

вектора

Ру,

Р-1 при эллиптическом движе­

нии не равен

нулю

ни при

каких значениях кеплеровых

элементов.

Од­

нако, когда время t0 соответствует моменту прохода КА над экватором, производные от координатных составляющих вектора состояния g по углу наклонения орбиты равны нулю. Кроме того, квадрат модуля ско­ ростных составляющих этого вектора для параболических орбит равен бесконечности, для гиперболических орбит его значение становится отри­

цательным.

Это говорит о

том, что для характеристики параболических

и гиперболических орбит необходимо использовать

систему других пара­

метров. В

частности, для

параболических орбит

достаточно сохранить

только пять параметров, так как шестой может быть определен из огра­

ничений,

накладываемых

на

существование параболической орбиты

[28].

Для гиперболических орбит большая полуось

теряет то

значение,

ко­

торое она имела при эллиптическом движении.

 

 

 

Для

эллиптических

орбит

квадрат модуля

третьего

вектора РУ3

Р-з

не имеет нулевых значений. Однако он показывает, что производные от

координатных составляющих вектора состояния gr

 

по

параметру

угло­

вого

расстояния

восходящего узла при i, а>=90°

и

£ = 0 (момент

вре­

мени / 0 соответствует проходу КА точки вертекса

при

полярной

орбите)

равны

•нулю. При

этих значениях наклонения орбиты

и

углового

расстоя­

ния перигея для случая, когда момент ^0 соответствует проходу КА над

экватором,

производные

от скоростных

составляющих

по

параметру

углового расстояния восходящего узла равны нулю.

 

 

 

Определяя квадрат

модуля

второго,

четвертого

и последующих

столбцов матрицы Р . приходим к выводу,

что эти столбцы не вырож­

даются в

нулевые, так как квадраты их

модулей определяются

суммой

квадратов

элементов последнего

сомножителя, который

ни

при

каких

значениях аргумента Е и других внутриплоскостных параметров не вы­ рождается в нуль-вектор.

 

При исследовании матрицы

Р

на

пропорциональность

между

столб­

цами получены данные, которые показывают,

что

между

1-м

и

2-,

4-,

5-,

6-м векторами (столбцами) выражения (VI.2.23)

пропорциональности

не

наблюдается. Эти векторы взаимно ортогональны, а так

как

ни один

из

них не

является ненулевым,

то

условие (VI.4.6)

не выполняется.

 

 

Следует

подчеркнуть,

что

при

исследовании матрицы

перехода

на

пропорциональность между

столбцами

особое

внимание следует

обратить

на

те векторы, которые содержат производные от составляющих

 

вектора

состояния

по

параметрам

одинаковой

размерности

(в частности,

по

уг­

ловым элементам),' и на те области задания параметров, в которых по

физическим соображениям тот или иной элемент

орбиты

теряет

смысл.

Так, исследования

выражения (VI.4.5)

для

второго и

шестого

столбцов

матрицы показали,

что оно становится

равным

нулю

при

а =

—1

для

е = 0 . Равенство

нулю указанного выражения

говорит о

том, что

 

в

этом

случае матрица

перехода становится особенной

и

определить все

шесть

165


кеплеровых

параметров не представляется

возможным. Это

н понятно,

так как в

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

ружности

параметр углового расстояния

перигея ю теряет

смысл. При

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

Выражение

(VI.4.5) для

второго

и

третьего

столбцов

матрицы пере­

хода

Р при условии,

если

долгота

восходящего

узла

теряет

физический

смысл

(/=0), удовлетворяет

требованию

(VI.4.6). Это говорит

о том, что

между

вторым

и третьим столбцами

матрицы

(производными

от

векто­

ра состояния JT по скалярным величинам Q и

со] наблюдается

пропор­

циональность.

Поэтому

кеплеровы

элементы орбиты

при

небольших уг­

лах наклонения малоэффективны. Для характеристики этих орбит необ­

ходимо применять

другую систему параметров. В частности, в

работе

[3] рекомендуется

вместо угла наклонения использовать косинус

этого

угла, однако эффективность такой замены требует дальнейшего иссле­ дования.

 

Для почти

круговых орбит

при

небольших углах наклонения наряду

с пропорциональностью между 2-м

и 3-м

наблюдается

также

пропорцио­

нальность

между 3-м и 6-м столбцами матрицы.

 

 

 

 

 

 

 

Для определения пропорциональности между столбцами

Р.,- и Р.,-

могут быть использованы и другие

критерии. В частности,

по

аналогии

со

скалярным

произведением

векторов

в трехмерном

 

евклидовом

про­

странстве можно ввести «угол» 0 между

векторами

P.j

и Р . ( ,

опреде­

лив его из соотношения

 

 

 

 

 

 

 

 

 

 

 

^ ^ ^ ( P j . P . ^ P y . P . / J C P ^ P . , ) .

 

 

 

 

(VI.4.8)

Значения

параметров, удовлетворяющих

условию

cos2 0=I,

определяют

ту

область

их возможного задания,

в которой матрица

Р

становится

син­

гулярной из-за пропорциональности ее столбцов. Исследования с исполь­

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

полу­

ченных

выводов

и

показали,

что

косинус угла

между

столбцами Р , 2

1 1

Р.з,

Р.2

и

Р.6.

Р-з и

Р-с

связан

соответственно

с

параметрами

i,

е, i

и е

такой

функциональной

зависимостью,

которая

при i,

е -+0

приводит

косинус угла к своему максимальному значению.

 

 

 

 

При

i,

е = 0 в

системе

кеплеровых

параметров,

как

в

системе

отсче­

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

КА,

между не­

которыми

ее

составляющими

имеется

линейная зависимость.

Этими со­

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

между столбцами

матрицы

перехода, но и столбцами

и

строками

-матри­

цы

вторых производных

АКФ

сигнала

по

определяемым

параметрам

q . Дл я

орбит с

эксцентриситетом или наклонением,

близким к

нулю,

матрица

вторых

производных АКФ сигнала

независимо

от

расположе­

ния

используемых

измерительных

средств

и их комплексирования

являет­

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

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

166


ства

кеплеровых параметров орбиты (последние преломляются

в свой­

ствах

матрицы

перехода

р), наблюдаются

в тех областях их

задания,

в которых один или несколько элементов

теряют

физический

смысл.

При этом для однозначного определения движения КА можно

использо­

вать

меньшее

количество

независимых обобщенных

параметров

(случай

вырождения эллипса в окружность или параболу) или взамен некоторых кеплеровых элементов орбиты ввести новые параметры (случай эллип­ тических экваториальных и гиперболических орбит).

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

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

VI.5. Особенности определения

параметров

в прямоугольных

и криволинейных системах

отсчета

 

 

/. Прямоугольные системы

 

В качестве исходной

системы параметров g

выберем систему эле­

ментов, характеризующих начальные условия движения КА в некоторой

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

связан­

ной с наблюдателем. Тогда, выбирая в качестве определяемых

конечных

параметров движения координаты и скорость КА в любой другой прямо­

угольной координатной

системе q, например,

в

различных

геоцентриче­

ских

системах

отсчета,

мы

обнаруживаем, что

вследствие

ортогональ­

ности матрицы-перехода, характеризующей преобразования

дифференци­

алов

(VI.2.5),

модуль

ее

определителя равен

единице. А

это означает,

что возможности уточнения параметров движения в произвольной пря­

моугольной

системе отсчета

не зависят от области задания координат

и скорости

и одинаковы для

всех прямоугольных координатных систем.

2. Криволинейные системы

Для удобства в качестве исходной будем рассматривать такую пря­ моугольную геоцентрическую систему отсчета, плоскость OXY которой совмещена с основной плоскостью криволинейных систем координат, а ось ОХ совпадает с некоторым характерным направлением, относительно

которого

в криволинейных системах

производится отсчет углов. При этом

под

матрицей

 

перехода Р

будем

подразумевать одну из матриц пре­

образования,

определяемого

соотношением

(VI.i2.-15). Матрицы указанного

преобразования

квазидиагональные,

поэтому их определители

выра­

жаются

через

определители

матриц

непосредственного

.преобразования

дифференциалов

W u , Wc ( |,

1 1 Wr -

Последние легко могут быть опреде­

лены, так как матрицы непосредственного

преобразования

дифференциа­

лов

являются

диагональными

(VI.2.12).

 

 

 

167