Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Сферическая астрономия

..pdf
Скачиваний:
5
Добавлен:
12.11.2023
Размер:
15.77 Mб
Скачать

мы Ьс + Ь о - Ь с х Ь а , причем Ьо является определяющей констан­ той, Ьс вычисляется на основе эфемерид БЕ405.

При таком определении разность между ТПВ и ТТ, измеряемая -на поверхности Земли, остается меньше 2 мс в течение нескольких тысячелетий.

Время Терн (аргумент эфемерид ОЕ405) с точки зрения потреби­ телей совпадает с ТВЪ, определяемым в настоящей резолюции.

Постоянный член ТПВо выбран для согласования нового време­ ни ТПВ с широко используемым разложением ТПВ - ТТ Файрхеда и Бретаньона. Однако наличие этого члена означает, что время ТБ В не синхронизовано с ТТ, ТСС и ТСВ на эпоху 1977, январь 1,0 ТА1 в геоцентре.

ОСНОВНЫЕ МАТЕМАТИЧЕСКИЕ ОПРЕДЕЛЕНИЯ

С.1. Матричная алгебра

Матрицей А размера т х п называется таблица скалярных вели­ чин:

( ац

«12

***

«1п ^

«21

«22

***

«2п

 

 

 

—[«#Ь

\«гп1

«т2

 

«ТПП/

п \т — число строк, п — число столбцов матрицы. Элемент ац называется элементом матрицы и расположен в г-ой строке и .7 -ом столбце.

Для матриц определены следующие операции:

1. Матрицы А и В размера т х п равны = В)> если для всех г

и2 равны их элементы: = Ь^.

2.Сумма матриц А и В размера т х п есть матрица С размера

т х п\

С = А + В \ Сц = -}- Ъц.

3.Произведение матрицы А размера т х п на скаляр а есть мат­ рица В размера т х п :

В = а.А\ Ьц = аа^ .

4.Произведение матрицы А размера ш х п н а матрицу В размера

пх к есть матрица С размера га х к:

п

С = АВ\ Сц =

1=1

Матрица АТУтранспонированная по отношению к матрице А раз­ мера га х га, есть матрица Ат = [а^] размера пхтп. Справедливы сле­ дующие соотношения:

 

-I- В )т = Ат + В т ;

(аА)т = аА т ; (А В )Т = В ТАТ .

 

 

 

Матрица А размера п х п называется квадратной матрицей по­

рядка га.

 

 

 

 

Квадратная матрица А называется диагональной, если все недиа­

гональные элементы а ц {г ф з) равны нулю.

 

 

 

Диагональная матрица А размера п х п называется единичной:

А

если все диагональные элементы равны единице: а ц

=

1,

ац

0 при г ф з. Если I

единичная матрица размера п

х

п,

то для любой матрицы В размера п х п справедливы равенства

1В = В1 = В .

Квадратная матрица А называется симметричной, если Ат = А, т. е. если

Квадратная матрица А называется невырожденной, если она име­ ет единственную обратную матрицу А~1Уопределяемую условием:

А - 1А = АА~г = I .

Квадратная матрица А называется ортогональной, если АтА = ААТ = /, т. е. если Ат = Л-1.

Если квадратные матрицы А и В одного порядка невырождены, скаляр а ф 0, то

(.А В У 1 = В ' 1 А ' 1; (а А )'1 = а~1А~1\ (Д "1)’ 1 = А.

Квадратная матрица не вырождена тогда и только тогда, когда ее строки (столбцы) линейно независимы.

С.2. Линейная алгебра

Система линейных уравнений

п

 

^

1) ...) тп

 

 

 

 

3=1

 

 

 

 

 

эквивалентна матричному уравнению

 

 

 

( а п

^12

• •

а\п ^

( х Л

А х \

а21

0.22

*• •

02п

Х 2

*>2

Ах — В или

 

 

 

 

 

\ а7п1

0>т2

• •

о,тп)

\х п/

\ Ь т /

Если матрица А квадратная и невырожденная, то матричное уравне­ ние имеет единственное решение:

х = А~1В.

Если матрица А является ортогональной и ее детерминант равен единице (с1е! А = 1), то линейное преобразование у = Ах называет­ ся вращением.

Система т уравнений

п

/г(жь ...,х„) =

- ы = 0

(*= 1,2

 

 

3=1

 

 

линейно независима, если из условия

^г/г(#ъ • • •, хп) =

0 при

всех значениях х х ,...,х п следует, что Ах

= Л2 = ... = Ат

= 0. В

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

Векторы ах, аг, аз линейно независимы, если из уравнения

Ахах + Лга2 Азаз = 0

следует, что Лх = А2 = А3 = 0. В противном случае векторы ах, аг, аз линейно зависимы и по крайней мере один из них, например, ах мо­ жет быть выражен в виде линейной комбинации ах = + Сз^з остальных векторов.

Любой вектор а в трехмерном пространстве может быть пред­ ставлен в виде разложения:

а = ахех + а2 в2 + азез

относительно линейно-независимых векторов ех,в2 ,ез. Эти векто­ ры называются базисом. В трехмерном пространстве числа ах, а2 , аз являются координатами вектора а в системе координат, определяе­ мой базисными векторами ех, в2 , ез.

С.З. Декартовы прямоугольные и сферические координаты вектора

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

и = их1 + иу$ + иге

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

Базисными векторами сферической системы координат являет­ ся правая тройка единичных векторов 1 г , 1 $ , 1 а , т . е. 1 Г х \$ = 1 а , на­ правленных в сторону увеличения соответствующих координат: г — расстояния, в — кошироты, Л —долготы. Разложение вектора и по базису имеет вид:

и= иг\г -I- ив\в + и \\\

иг, ив, их компоненты вектора и в базисе 1г, I#, 1Д.

Декартовы координаты вектора и = (их, иу, иг) выражаются че­ рез сферические координаты (г, в, Л) следующим образом:

г = |и|.

Единичные векторы 1Г, 1 0 , 1л в декартовых координатах имеют вид:

.

1

ГСОЗ в СОЗ

д г

16 =

 

С08 в 81П Л

г д в

— 81П О

1Л =

Г 81П 0 д \

Преобразование между компонентами вектора в декартовом и сферическом базисе имеет вид:

 

Лиг \

^81П0СО8 А

СОЗ 0СОЗ А

—зтА

П у

= м щ

; м =

81П в81П А

СОЗ в81П А

соз А

и /

\ « л /

\

С08 0

—зт 0

0

Обратное преобразование имеет вид:

^иг\ ( их^

=М Т

СЛ . Элементы дифференциального и интегрального

исчисления

Дифференциал функции у = /(х ) в точке х , если он существует, равен:

йу = ~^.^х ~ /'(х)с1х,

/ '( х ) = сХу/йх производная функции.

Дифференциал функции у = /(а?1 , #2 , ... , хп), если он существу­ ет, равен:

л

, . д / ,

0 / ^

* =

' +

+ •••+

причем частные производные

вычисляются в рассмат­

риваемой точке.

 

 

Если / ( х) — действительная функция, имеющая в интервале а < х < Ь п-ю производную /( п), то

/(х ) = /(а )

+ /\а ) (х - а) + ^ /" (а )(х - а ) 2 + ...+

+ т—

/ ( п-1)(ж - а)"-1 + Яп(х), (а < х < Ь),

( п -

1)!

где Яп(х) называется остаточным членом.

Градиентом скалярной функции Р(г) = Р (х, у , г) называется векторная функция, определяемая формулой:

1 п/

\

— 1 -|-

9 Р л

д Р щ

 

 

 

§гаа Р(г) =

 

+ — к в декартовых координатах,

, ^

ч

дР .

1 Э Р .

1

 

дР

 

§гас1 Р(т) =

-д -1 г +

- яЗ-1 * + — г-

5

дт-1Л в сфер, координатах.

 

 

дг

 

г дО

г

0 д \

 

Полный дифференциал

скалярной функции .Р(г) = Р (х, у, г \

соответствующий перемещению точки на Лх = Лх И - Лу$ + Лг к,

 

 

ч

дР

 

дР

дР

,

,

 

ЛР{г) =

-т^сЬг +

 

-I-

 

= Л* •§гас1 ^(г).

Дифференциал Лг радиуса-вектора г вдоль кривой С, описывае­ мой уравнением

{ х =

2/ = 2/(0

2 = 2(<),

определяется в каждой точке кривой г = (х{1), у(<), г(1)^ формулой:

= <1х\ + йу$ + <1гк = { ^ 1 + ^

Вектор йг направлен по касательной к кривой С в точке с радиусомвектором г.

Квадрат элемента длины:

Аз2

= Лх2 + Лу2 + Лг2

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

Аз2

= Лг2 + г2АО2 -|- г2 з т 2 0 ЛХ2

в сферических координатах.

Преобразование дифференциалов из сферической в декартову систему координат имеет вид:

С.5. Криволинейные координаты

Если в области V трехмерного евклидова пространства заданы соотношения, ставящие в соответствие каждой точке (х, у , г) тройку чисел х \ г = 1,2,3, причем

X1 = Л (ж, 2/, г), х2 = / 2(х, у, г), х 3 = / 3(х, у, г),

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

Условие хг = /г(х,у,г) = сопз!, г = 1,2,3 определяет коорди­ натную поверхность. Две координатных поверхности, соответству­ ющие различным координатам х г, х*, {ъ ф з) пересекаются по коор­ динатной линии, соответствующей третьей координате х к.

Единичные векторы \\(х1 ,х 2,х 3)у \2(х1,х 2ух3)у \ъ(х1 ,х 2,х 3)ука­ сательные к координатным линиям а?1, х2уж3, являются локальными базисными векторами.

В качестве локальных базисных векторов можно выбрать тройку направленных по координатным линиям векторов (не обязательно единичных), которые определяются уравнениями:

е \ { х 1 , х 2 , х 3 ) = у / д п \ 1 , е 2 ( х 1 , х 2 у х 3 ) = у 'р м Ь ,

егС®1,®2,®3) = у^зз^з,

где ргг — компоненты метрического тензора:

Локальные базисные векторы ех, ег, ез могут быть выражены че­ рез орты 1^, к декартовой системы координат по формулам:

д х ,

ду , дг

е< = а +

+

Функции

являются направляющими косинусами орта

вг по отношению к декартовым осям О#, Оу, Ог.

 

В качестве локальных базисных векторов можно выбрать тройку

векторов е ^ е ^ е 3:

 

 

е 1(аг1 ж2, х3)

е2 х е3

ез X в!

з2 (гг1, гс2, ж3)

 

 

V911922933 ?

 

в 1 х е 2

е3^ 1, ж2, х3)

у /9119229^3 ^

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

В базисе е х ,

е

2координаты, е з

вектора а называются ковариант-

ными координатами:

 

 

 

 

 

 

 

а ,г =

а*(ж1, ж2,ж3) = а • е^,

г =

1,2,3,

 

а в базисе е

1 , е 32, контравариантными:е

 

 

 

 

аг = а*(ж1,ж2,ж3) = а • е

гг,=

1,2,3.

 

Справедливы формулы:

 

 

 

 

 

 

 

 

 

 

 

1,

при г = к

 

©г *

9г^1

|®г| — у/ 9ггч

 

{

при г ф к.

 

 

 

 

 

 

 

О,

 

В частном случае прямоугольных декартовых координат х 1 = х ,

х 2 = у> х3 = г имеем:

 

 

 

 

 

е х

=1 =

е

1 х = 2 =1 , 2 е= 2 е =1

 

3 =

3 е=е 3 =1

к .

Система криволинейных координат х 1, х 2, х3 является ортого­ нальной, если

дгк(х1 ,х 2,х 3) = 0 при г ф к

вкаждой точке (ж1, х2>ж3). Координатные линии, а значит, и векто­ ры локального базиса 1х, 12 , 13 будут перпендикулярны друг к другу

вкаждой точке.

Элемент объема в криволинейных координатах:

= у/дЛх1вх2Лх3, у /д = у/9и 922933

Интеграл по объему от функции

/(х,у,г) =/[х(х123),у(х123),г(х12,х3)]

по ограниченной области V

 

Д

/ у( , жг)6У, = Д

Д у,г)йхЛуйг/ ( ж ,

=

=Д /[х(хг,| * х2, х3),у(хг, х2, х3),г(хг,х2, х3)\^/дЛх1Лх2Лх3.

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

В сферических координатах ^ = г2 з т 0.

С.6. Сферические функции

Гравитационный потенциал V во всех точках, находящихся на поверхности и вне Земли, удовлетворяет уравнению Лапласа:

 

 

У2С/ =

д2и

д2и

д2и

0.

 

 

 

 

 

дх2 + ду2

+ ду2

 

 

 

Оператор V называется «набла». В сферических координатах (г, 0, Л)

уравнение Лапласа имеет вид:

 

 

 

 

 

^

= 1 А

/ ^

+ _ 2

д

( . йд и

 

1

д2и

_____

8,11 " г й

+

г2 з т 2 0 9А2

 

г2 9г \

)

Г2 81П 0 90

Решение уравнения Лапласа есть:

 

 

 

 

 

 

 

г1

 

 

 

 

 

 

 

V = \

или

^ х Р{п{соз0)ег т \

 

 

 

 

I

 

 

 

 

 

 

где /, га — целые числа, причем / > 0 и |га| <

г =

\ / —1- Функции

Р/т (/х) называются присоединенными функциями Лежандра степе­

ни I и порядка га. Функции

 

есть решения дифференциально­

го уравнения:

 

 

 

 

Л2Р

ЛР

 

~2

1

1(1 + 1 ) -

 

Р = 0.

(1 _ А* )“ГТ “

— Н

1 - /х 2

 

(1/1