книги из ГПНТБ / Селезнев В.П. Инерциальная навигация летательных аппаратов
.pdfВ соответствии со структурной схемой на фиг. 3.4,6, получаем следующее дифференциальное уравнение:
(Р2 + ki М з Р : k\ k 2 g + |
-= pV |
ki - |
+ |
|
кх кг До ' |
(k4 i |
/е, к2 кяр) Д,3. |
(3.25) |
|
При соблюдении условия инвариантности k\k2 ---- —*= 0 и вве- |
||||
дении обозначений |
|
|
|
R |
|
|
|
|
|
к \ кг g к* — т « у-2- |
<чо = |
|
kyk2k3 = 2dm0х |
получим уравнение, совпадающее с (3.23). В зависимости от знака можно увеличить (ki <4 0) или уменьшить (k4 > 0) период астроинерциальной системы по сравнению с периодом М. Шулера, не нарушая условия инвариантности (3.22) с точностью до погреш
ности навигационного корректора.
Подбором коэффициента k3 может быть получен любой декре мент затухания d. Для быстрого затухания колебаний необходимо иметь d = 0,7 : 1,0. Внешняя информация может подаваться в астроинерциальную систему не непрерывно, а только в короткие промежутки времени: например, при начальной выставке (для лик видации начальных погрешностей), для ликвидации накопивших ся погрешностей в режиме «памяти», когда видимость небесных светил отсутствует. Для того чтобы процесс коррекции был по воз
можности |
коротким, необходимо сокращать |
время |
переходного |
процесса. |
С этой целью необходимо подбирать |
х2 > |
1 за счет из |
менения передаточного коэффициента k4. Но значительное умень шение периода инерциальной системы приводит к тому, что она начинает возмущаться помехами от внешней информации.
Погреш ности комплексной астроинерциальной системы
Погрешности возникают от ошибок внешней коррекции Ар и
Ай, поворотных и других ускорений Аа. Рассмотрим два наиболее характерных случая:
1) погрешности |
Д(3 и Ар |
и ускорение До |
не зависят от вре |
мени и являются постоянными величинами; |
|
||
2) погрешности |
Д|3 и Д,3 |
являются случайными функциями |
|
времени. |
|
|
|
В первом случае установившаяся погрешность инерциальной |
|||
вертикали будет: |
|
|
|
Руст- |
- ^ 1 д |
и - — ДР4-—2. |
(3.26) |
|
х2 |
а>0 х |
|
Отсюда следует, что при больших значениях коэффициента измене ния частоты колебаний х(х > i) величины погрешностей, вызван-, ных ошибками в информации о скорости летательного аппарата
144
(*Р> и поворотными ускорениями (Да), существенно уменье шаются. Преобладающее значение приобретает погрешность, вы званная ошибками навигационного корректора в измерении коор
динат местонахождения |
(руст Дз). |
О влиянии изменения частоты |
колебаний инерциальной |
вертикали |
на погрешности коррекции |
можно судить из следующего примера. Современные радионавига ционные системы определяют местонахождение летательного аппа рата с точностью до ДS x = 1 км на дальностях в несколько сот километров, а допплеровские радиолокационные системы измеряют путевую скорость с относительной ошибкой до 0,3 — 0,4%. Если коррекция астроинерциальной системы ведется одновременно: по координатам — от радионавигационной системы, а по скоро сти — от допплеровской радиолокационной системы, то
др = Д5, |
1 |
= 0',54; |
Др |
0,0041/, |
R |
6,37-10G |
|
|
R |
Положив скорость летательного |
аппарата |
V = |
1 км/сек, получим |
|
Др = |
0,22-10-2 угл. мин/сек. |
|
При таких погрешностях навигационной информации получим сле дующие погрешности астроинерциальной системы:
для х |
1 |
2d |
0,004 V |
|
|
|
|
■3'.7: |
|||
|
|
I Р у с т | |
R |
||
для х = |
10 |
>uo |
|
||
2d .0,004Ух = 0',54 + 0',37. |
|||||
: Р у ст I = А ? |
|||||
|
|
Ю»0 |
R |
|
Из этого примера следует, что десятикратное увеличение частоты колебаний инерциальной вертикали позволяет существенно умень шить ее погрешности до уровня, не превышающего одной угловой
минуты.
Произведем оценку погрешностей системы для случая, когда
погрешности Др и Др являются случайными функциями време ни. Оценивая эти погрешности их спектральными плотностями Д и Д5 р , можем написать выражение для среднеквадратической по грешности вертикали
—2 |
1 |
оо |
|
|
|
|
|
|
! (гю) |2 |
(ш) du |
те I' |
Ф2 (lot) |2 ДSfi (и>) d и), |
(3.27> |
||||
е* = |
7Г |
|||||||
|
оо |
|
|
|
|
|||
где |
|
Ф, (1ш) ]2 = |
|
Ц)р(х2 — |
I)2 |
. |
(3.28) |
|
|
|
|
|
|
||||
|
|
|
(<оо X2 — |
(О2)2 + |
4 d 2 0)0 * 2 “ >2 |
|
|
|
|
|
Ф2(/«))р = |
|
Ad 2 а>о х2___________ |
|
(3.29) |
||
|
|
|
|
|
|
(о>0 X2--<02)2 -f- 4d2 Мр X2 0)2
10. В. Г1. Селезнев |
145 |
Спектральные плотности ошибок Д53(ш) и \Sa («4 определяются н основном свойствами измерительных систем, дающих величиныS x и V и колебаниями объектов, на которых установлены системы Влияние колебаний объектов на величины погрешностей особенно значительно в тех случаях, когда результаты измерения зависят от положения объекта по отношению к вертикали места. В таких слу чаях колебания, являющиеся случайными функциями времени, определяют в значительной мере спектральные характеристики по грешностей Д^е и Д . Допплеровский измеритель путевой ско рости и является такой системой, погрешности которой определя
ются колебаниями объекта. При из мерении координат местонахожде ния объекта при помощи радиона вигационных систем влияние коле баний менее значительно. Учитывая
сказанное, можно взять |
для спект |
||
ральных |
плотностей |
Д5р(ш) |
и |
Д5р(<») |
следующие выражения |
[8]: |
Фиг. 3.5. Графики спектральных плотностей
4х, а»
(3.30)
4/., 32 (г2 + |
- 2) |
Д.^э (и>) |
(3.31) |
( г 2 ю2)2— |
« 2 |
где
,2 ^ >4 + Й2,
)-i — параметр затухания; й — резонансная частота.
На фиг. 3.5 даны графики функции ДS$ и ДS? .
Подставляя выражения (3.28) — (3.31) в (3.27), получим после интегрирования и преобразований
|
|
|
~2 |
“ 2 —2 |
5 |
(3.32) |
|
|
|
|
*р 'Г |
||
где еЛ.р и |
— |
среднеквадратические значения |
погрешностей |
|||
инерциальной |
системы, проистекающие от неточной |
информации |
||||
о координатах и скорости объекта. При этом имеем: |
|
|||||
72 |
‘"о |
(*2 — 1 )2 (2й?% -л + |
),,) |
(3.33) |
||
Едгр — |
|
a»o*2 + 2</Q0xX1) |
|
|||
|
Хз(х? + |
|
||||
“ 2. |
16rf2o2 [«ш2'/.2 Дд + Kxz 2(Дс 4- «>2 Х 2 Д0) |
|
||||
|
|
|
пг2 ш0 у.Д |
|
||
где |
|
|
|
|
||
|
|
у-4 ) 2 + |
(2А — v2 ) (2/гг2 — v 2 io 2 -a 2) ; |
|||
д = |
( 24 |
— |
Да = 22[—\z*— u^x4) J- (2Л — v*)(z2-— v2)];
146
Лс - - |
{z* - |
*«) -J- (2 h - v-) -f (2hz2 - |
®« xM; ' |
|
|||
= |
Z i — |
<■>< 7* - I - z2 (2h |
V2 ); |
|
|
||
/г == V 2 d2 — 1— 2 d Y d 2— 1 + \' 2 d — 1 I - 2d Y d 2 - T ; |
|||||||
2A -- |
2<o2x2 (2a!2 - 1), |
|
|
|
|||
v2- 2 z 2 - |
4ш*. |
|
|
|
|
|
|
Проведем |
анализ |
погрешностей гх? и гх%. |
Из выражения |
||||
(3.33) следует, |
что при |
/•= |
1, т. е. в случае системы с периодом |
||||
Шулера, позиционная погрешность |
еХд обращается в нуль. |
С воз |
|||||
растанием х погрешность |
гх? |
возрастает и при большом |
сокра |
щении периода инерциальной системы стремится к своему макси
мальному значению s ^ max -- 2з|з. |
|
|
||
Скоростная погрешность также зависит от -/.. |
Если предполо |
|||
жить, что d |
1 и X, < (о0, то при изменении / от |
х = 1 до •/. -> оо |
||
погрешность |
|
|
|
|
|
|
« Ь - 4 . |
|
|
|
|
°>о 2:2 |
|
|
при х = 1 начинает возрастать, достигая максимума |
||||
|
- 2 |
при |
х —— |
|
|
|
|
шо |
|
При дальнейшем возрастании коэффициента погрешность на |
||||
чинает уменьшаться, достигая |
нуля при х-^со. |
Такой характер |
||
зависимости погрешности ех^ |
от х |
объясняется тем, что с ростом |
частоты инерциальной системы амплитудные значения частотной характеристики убывают обратно пропорционально х2 при общем падении спектральной плотности на частотах выше резонансной. Отсюда следует сделать вывод о необходимости увеличения часто ты инерциальной системы до значений, превышающих резонанс
ную частоту 12 . |
При таком увеличении частоты инерциальной си |
||||
стемы появятся |
неучитываемые погрешности от угловых колеба |
||||
ний самолета. |
|
|
|
|
|
Выбор |
коэффициента х |
можно |
обусловить |
требованием |
|
«>0 у. < 12. |
12 •= |
0,02 11сек, то |
х < 16,3. По-видимому, десятикрат |
||
Если |
|||||
ное уменьшение периода будет приемлемым. |
|
||||
Если параметр затухания |
> 12, |
то функция зд., |
не имеет |
максимума по х и монотонно убывает при возрастании х.
Ю* |
14/ |
§3.5. КОМПЛЕКСНЫЕ АСТРОИНЕРЦИАЛЬНЫЕ СИСТЕМЫ НАВИГАЦИИ
ВТРЕХМЕРНОМ ПРОСТРАНСТВЕ
Для навигации летательного аппарата в трехмерном простран стве необходимо знать три координаты местонахождения (X, Y, Z) и составляющие скорости (X, Y, Z) относительно избранной системы координат.
Применение гироинерциальных и астрономических систем на вигации имеет следующие особенности. Гироинерциальная система достаточно точно измеряет скорости полета, а измерение координат местонахождения производится с накоплением значительных по грешностей. Астрономические системы отличаются высокой точ ностью изменения координат местонахождения, но измерение ско рости,получаемое дифференцированием координат, дает чрезмер но большие погрешности. Объединение этих систем в общую комп лексную систему навигации позволяет взаимно скомпенсировать их недостатки, а также значительно улучшить динамические характе ристики астроинерциальной системы и повысить ее точность.
Структурные особенности комплексной астроинерциальной системы
Астроинерциальная система включает фотоследящие системы за тремя небесными телами или звездами, гироскопическую плат форму, три акселерометра, оси чувствительности которых ориенти рованы по осям выбранной системы координат, интеграторы для вычисления скоростей и координат и три вычислителя. На фиг. 3.6 представлена возможная структурная схема комплексной астрои нерциальной системы1
Однократное интегрирование сигналов акселерометров ах, ау, аг дает скорости полета, а двухкратное интегрирование — мгно венные значения координат. Полученные значения координат по
даются в |
вычислитель Bi |
для |
формирования |
компенсационных |
сигналов |
g XK, g yK, g z к |
и в |
вычислитель В2, |
дающий угловые |
координаты для управления телескопами Ti, 'Гг, Т3._Сигналы рас согласования AS, Avj, АС между задаваемыми S, т], С и фактиче
скими S, т), С угловыми координатами пеленгуемых небесных тел или звезд поступают в вычислитель В3 поправок АА\ АУ, AZ к ко ординатам местоположения. Эти поправки передаются через уси лительные звенья цепи обратной связи в схему измерения и инте грирования ускорений.
Особенностью схемы является то, что сигналы обратной связи^ суммируются с выходными сигналами акселерометров и первых' интеграторов, вследствие чего, как будет показано ниже, удается обеспечить устранение неустойчивости, сокращение периода и демп фирование инерциальной системы.
1) Б о д н е р В. А., С е л е з н е в В. П. К теории невозмущаемых систем с тремя каналами автокомпенсации ускорений от сил тяготения, корректируемых внешней информацией. Известия АН СССР, ОТН, «Энергетика и автоматика», 1960, № 3.
148
|
|
|
|
|
Кгх |
Кгх |
IX |
0 |
В г |
|
|
|
|
|
|
|
|
||
|
|
|
9 ^ |
Р |
Р |
|
SX |
|
|
|
d Y |
|
|
|
|
|
|
||
|
|
Нзу |
|
|
% |
|
|
||
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
л |
|
|
5, |
|
MZ |
|
*4; |
»Ук. |
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
х?я |
I Kzi |
|
|
|
|
|
|
|
3z ZL |
tE |
|
|
|
|
|
|
|
|
|
|
|
|
||
Ф ii г. |
3.6. Структурная схема комплексной аетроинерциальной систе |
||||||||
мы, осуществляющей измерения в трехмерном пространстве: |
|||||||||
Т ,, Т 2, Т3 — телескопы: |
Bj, В2, |
Вя — вычислители; k\x , k<iK f |
ф t |
k^z — пере |
|||||
|
|
|
|
даточные коэффициенты |
|
|
|
||
Вычислитель Bt должен формировать сигналы компенсации |
|||||||||
ускорений от сил тяготения по формулам: |
|
|
|
||||||
|
П |
|
|
|
П |
|
|
||
Я х к |
Ц |
/ |
' т |
( S , |
• - * / ) . |
V |
/ - ^ < - S |
' y - К.), |
|
|
. |
, |
Ri |
|
|
. |
, |
Ri |
|
|
|
|
|
|
|
(■=. 1 |
|
(3.34) |
|
|
|
|
|
■я, к |
У / |
^ A S z - Z i ) , |
|
||
|
|
|
|
|
|
||||
|
|
|
|
|
i-1 |
R i |
|
|
|
|
|
|
|
|
|
|
|
|
149
где т„ Х (, Yt, Z — |
соответственно |
массы- и координаты |
центров |
|||
/?, — |
небесных тел; |
летательного |
аппарата до |
не |
||
расстояние от |
||||||
S x, S y , S z — |
бесных тел; |
местоположения, |
измеренные |
на |
||
координаты |
||||||
|
вигационной |
системой; |
|
|
|
|
f - - постоянная тяготения. |
|
S 4, S z |
и |
|||
Вычислитель В2 на основе входных координат 5,., |
координат пеленгуемых небесных тел определяет угловые коорди наты небесных тел относительно инерциальной системы координат. Связь между координатами летательного аппарата и угловыми ко ординатами небесных тел определяется уравнениями поверхностей положения:
? - |
F x [X , |
У , |
Z, |
Л',, |
У/, |
Zf); |
|
rt = F 2 (X, |
У, |
Z, |
X i% УР Z ()\ |
|
|||
C= |
f'8(AC |
У, |
Z, |
X it |
Kf, |
Z,). |
(3.35) |
Примеры этих уравнений можно найти в § 2.11. Поскольку вместо
координат X, Y, |
Z вычислитель В2 получает измеренные координа |
||
ты S x, S y, S z, |
то вычисленные |
им |
координаты $, т(, " отлича |
ются от измеренных телескопами координат £, rt, С |
|||
На основе рассогласования |
Д$ = |
f — Дт) = ^’—v(, Л^==^_Л-* |
в вычислителе В3 производится вычисление поправок ДА', ДК, \Z .
при этом используются уравнения (3.35), из которых для прира щений получаем
_ d F x\ Х + |
dFJ Д У + dFx ДА; |
|
||
дХ |
дУ |
l z |
|
|
_ д Г 2 ДА'4- |
ДК + |
dF, ДА; |
|
|
~ дХ |
дУ |
dZ |
|
|
dFs |
дУ ДГ-(- |
dF, |
ДА. |
|
дХ |
dZ |
|
||
Отсюда |
Dv |
|
|
Г) |
|
д7 |
|||
D |
У, |
-JL |
||
D ’ |
^ |
D |
||
где |
|
|
|
|
dF\ |
dX l |
dP\ |
|
|
dX |
|
|
||
дУ |
dZ |
|
|
|
dF2 dF, |
dF2 |
|
|
|
dX |
дУ |
Tz |
|
|
dFt |
dF, |
dF, |
|
|
aX |
dY |
dZ |
|
|
(3.36)
(3.37)
150
а определители Dx, Dy, D z получаются из D заменой соответствую щих столбцов правыми частями уравнений (3.36). Таким образом, уравнения (3.37) определяют программу работы вычислителя Вз.
Уравнения комплексной астроинерциальной системы
Предположим, что система координат xyz, по осям которой на правлены оси чувствительности акселерометров, повернута отно сительно инерциальной системы координат x0y0Zo на углы а, {3, у (фиг; 3.7). Измеряемые акселерометрами ускорения будут:
ах = |
+ Аа*г : |
|
|
«у= (hu + |
А%г : |
’ (3.38) |
|
а* - |
+ |
дя*г • |
|
где |
|
|
|
ax,= X + g x\ |
|
||
= |
Г + g,r |
|
|
a2o = |
’Z - ' - g z. |
(3.39) |
Ф и г. 3.7. Отклонение системы координат (xyz) акселе рометров относительно инерциальной (.vof/nzo) системы * координат
Для малых углов я, |3, у получаем
к а л-г — 8а — |
ааг„; |
|
Аа>у = |
Т«-’„ — Ря»-.,; |
|
Да*г = |
вл.,-0- |
( 3. 40) |
151
Из структурной схемы (фиг. 3.6) получаем выражения для из
меряемых скоростей: |
|
|
|
|
|
|
$х = (ах — ёх К — Кх АХ) — - Кх Кх АХ + Х 0-> |
|
|||||
|
|
|
Р |
|
|
|
Sy = (а, - |
g y к - |
Л3у ДК) ^ |
— Лг3у kiy ДК + К0; |
(3.41) |
||
|
|
|
Р |
|
|
|
S z = [ a z - |
g 2 к - |
А3г AZ) |
- |
k3z kiz AZ + А |
|
|
и координат: |
|
|
Р |
|
|
|
|
|
|
|
|
|
|
|
5 , = |
^ 5 , + |
^ |
; |
|
|
|
|
|
D |
|
|
|
|
5 у = |
^ - 5 у + |
Г0; |
(3.42) |
||
|
|
|
Р |
|
|
|
|
S z = - ^ S z + Z Q, |
|
||||
|
|
|
р |
|
|
|
где Хо, Уо, Z0 и j 0, |
Y0, Z n— начальные значения координат и ско |
|||||
ростей. |
|
|
|
|
|
|
Выражения (3.42) являются уравнениями астроинерциальной |
||||||
системы. Запишем эти уравнения в приращениях. Полагая |
|
|||||
АХ = Sx — X; |
ДГ = 5 У- К ; |
A Z = S z - Z ; |
(3.43) |
|||
Аёх== ё х ~ |
ёхк, |
Аёу = ё у ~ ё У«> Аёг = ё г ~ ё г |
(3.44) |
и Кх Кх — КуКу — |
— |
получаем вместо (3.42) |
|
||
(р2 + |
2dx о хр + |
22) ДА"= |
Д g x + |
АаХг; |
|
(р2 + |
2dy ilyр + |
22) ДК = |
Agy + |
ДаУг ; |
(3.45) |
(p2 + 2 dz Qzp + U2z) b Z = b g z + baZT,. |
|
||||
где |
|
|
|
|
|
i l l ^ k s x , |
il2y |
|
ill |
- К г |
|
2 dx 2 ^ = Кх kix, |
2 dy 2 у = /?3у kiy, |
2 d ZQ г = k%z kiz. |
|
Погрешности системы автокомпенсации ускорений
'от сил тяготения
Погрешности Agz, Д^у, Д^г, формируемые на основе сигналов S x, S v, S z, отличны от нуля и равны:
A g x ^ — ^ A X + AAY+BAZ-,
Аёу — ААХ — Т)2 ДГ + CAZ; |
(3.46) |
Д g z = ЯАХ + C A Y - y]2 AZ,
152
где |
■-Е/5 1 - 3 |
'X - X t 2 |
|
|
/=■1 |
R l |
Ri |
|
|
lTli 1 - 3 |
Y - Y, 's |
|
/«1 |
R l |
Ri |
|
П |
|
z —zt 2 |
|
Е Л 1 - 3 |
||
|
|
я / |
Д, |
П« , ( X - X t) ( Y - Y t) .
л= з £ /
ы\ |
Ri |
R l |
|
|
|
||
n |
m, (АГ-Л-,) (£ -■ ? ,). |
||
B = |
|||
Л |
|
||
- » E |
R l |
||
r-1 |
Ri |
||
|
|
||
П |
|
{Y - Уi)(Z ~ Z , ) |
|
c-«S/ |
|||
я? |
|||
i-i |
|
||
|
|
||
Подставляя (3.46) в (3.45), получим |
|||
{p2 + 2d, Й,/> + <4 + -^) ДА' - |
ЛАК - /^ Z = Aa,r ; |
(3.47)
(3.48)
— ЛДА" Ч- (/>•’ -!- 2dу 9.ур |
|
| (->2т)р ДГ — CAZ = |
Д%г ; |
(3.49) |
- 5ДЛ" - СДГ + (р 2 f 2 |
^ |
‘2 , р 4- <Z2 -г,2) AZ - |
Aa*r . |
|
Параметры системы могут быть выбраны из условия устойчивости системы, хорошего демпфирования и сокращения периода.
Погрешности в измерении ускорений Ай*,., АаУг , ДaZ{ , вызван
ные отклонением телескопов от направлений на звезды, в режиме слежения могут быть сделаны пренебрежимо малыми.
Рассмотрим устойчивость системы автокомпенсации ускоре ний от сил тяготения в предположении, что отклонения телескопов отсутствуют (До*г = ДаУг = а2г = 0), а коэффициенты системы
(3.49) меняются медленно по сравнению с временем протекания процессов в системе автокомпенсации. В этих предположениях устойчивость системы определяется свойствами характеристическо го уравнения
F A P ) |
- |
В |
|
А (Р) =■ — А |
- |
С - F x (p)Fy(P)Fz (p)~ |
|
~ В |
F Z(P) |
|
|
~ C ' F X{P) |
&Fy(p) - A2 Fz { р ) ~ 2 А В С = 0 , |
(3.50) |
153