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

shimbirev_b_p_teoriya_figury_zemli

.pdf
Скачиваний:
100
Добавлен:
04.06.2015
Размер:
14.99 Mб
Скачать

Формула (111.16) аналогична

формуле

 

(III.7),

разница

лишь

в

замене

соз 0 на соз ф. Поэтому,

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

ранее выводы,

получим

 

 

 

 

 

Д "

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(111.18)

 

 

рП+1 Р „ ( с о з ф ) ( р > Д ) ,

 

 

 

 

 

 

 

п=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

п=О2 Дп+1

Л ( с о з ф ) ( Р < я ) .

 

 

 

 

 

 

 

(III.19)

Выше было доказано, что

рпР (соз 0) — целый однородный

гармонический

многочлен в переменных

х, у,

г;

это

же

свойство

 

должно

принадлежать

 

Р(х,у,г)

и

рпРп

 

(соз

 

г(;),

 

так

как

эта

функция

 

легко

получается

из

рпРп (соз

0)

путем

 

 

 

 

 

 

простого преобразования координат точ-

 

 

 

ки

М,

которое

сводится

лишь

к пово-

 

 

 

роту

осей,

а

это

не

может

 

повлиять

 

 

 

на

алгебраические

свойства

многочлена

 

 

 

относительно

 

координат

 

х,

у,

г

точ-

 

 

 

ки

Р.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом рпРп

(соз ф)

можно

 

 

 

рассматривать как шаровую функцию сте-

 

 

 

пени

п

 

в

переменных

х, у, г,

а

Р п

 

 

 

(соз ф) — как

 

сферическую

 

функцию

 

 

 

координат 0 и Я. Выражение Рп

(соз ф)/р"+1

 

 

 

будет являться шаровой функцией сте-

 

 

 

пени —п — 1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Определим теперь явный вид шаро-

 

 

 

вых и сферических функций в общем

 

 

 

случае. Возьмем произвольную гармо-

 

 

 

ническую

функцию

Уг

(р,

0,

К)

внутри

 

 

 

сферы радиуса В и разложим

ее

в

бес-

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

(11.44),

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

(11.42)

в

виде

 

 

 

 

 

 

 

 

 

 

 

Разложим ядро интеграла по степеням р. Для этого используем выражение (111.19)

ЯП+1 Рп( С081])).

Дифференцируя его по В, получим

70

Следовательно,

00 СО

п=О

п=О

 

СО

п=О

Подставляя полученное значение ядра интеграла в (III.20) и интегрируя почленно, найдем

М р ,

е,

=

 

а

 

 

 

п= о

 

 

Введя элемент поверхности

сферы

единичного радиуса йсо =

и обо-

значив

 

 

 

 

 

 

Упф,

=

ш

Я,')Л, (соз ф) Л»,

(111.21)

получим

 

 

 

 

 

 

оэ

 

 

 

 

 

 

 

 

 

VI (Р, е, я) =

^ Л 0 '

(П 1 -22)

 

 

 

п=0

 

Функцию Уп

(0, А) можно

рассматривать как произвольную

сферическую

функцию степени п, так как Рп

(соз ф) является сферической функцией коорди-

нат 0 и Я и это свойство не может измениться при интегрировании ее по переменным 0' и Я'. Поэтому члены ряда (111.22) —рп Уп (0, к) будут являться шаровыми функциями степени п.

Решим аналогичную задачу для гармонической функции Уе (р, 0, к) вне сферы радиуса В. Для этого в интеграл Пуассона (11.41) подставим значение его ядра, используя (11.44),

о

Мы имели (111.18)

п= О

после дифференцирования получим

- й г ( Т Ь 2 п ^ » < с о в * > '

п=О

следовательно,

с»

п=0

71

Подставляя это выражение в (111.23) и почленно интегрируя, найдем

 

со

 

УДР, 0,

=

А,')Л,(С081>)ДГ,

 

п=О

ст

вводя опять элемент поверхности сферы единичного радиуса и используя соотношение (111.21), получим

оо

п=О

0П+1 (11124. )

Каждая шаровая функция этого ряда имеет отрицательную степень

 

 

 

(х,

у, г) =

 

 

.

 

 

Пределы рядов (111.22) и (111.24) при р -V Я должны быть равны значению

заданной на сфере функции

 

 

 

 

 

 

 

 

 

ПшУДр,

0, А) = /(0, X),

 

 

 

р->-В

 

 

 

 

 

 

 

 

П т У Д р , 0, А ) = / ( 0 , я).

 

 

 

р-*-В

 

 

 

 

 

 

 

В то же время, полагая

(111.22) и (111.24) р

Я,

получим

 

 

 

 

 

 

 

 

00

 

 

Н т V,

(р,

0,

X) =

V, (Я,

0,

X) =

2

Уп (9-

*•)

Р-УВ

 

 

 

 

 

 

п=О

 

 

и

 

 

 

 

 

 

 

 

 

П т У Д р ,

0,

Х) = Уе(Я,

0,

X) =

2

У„(в,

Я).

Р-+К

 

 

 

 

 

 

П=1

 

 

Следовательно, разрыв, который имеют интегралы Пуассона (11.41) и (11.42)

при р = Я, устраняется в рядах (111.22) и (111.24). На

поверхности сферы о

оба ряда (111.22) и (111.24) дают

 

00

 

/(9, А) = 2 Уп (0, Л).

(111.25)

п=О

 

§ 15. РАЗЛОЖЕНИЕ ПРОИЗВОЛЬНОЙ ФУНКЦИИ В РЯД ПО СФЕРИЧЕСКИМ ФУНКЦИЯМ

Полученную в предыдущем разделе формулу (111.25) можно рассматривать как разложение произвольной функции / (0, X) в ряд по сферическим функциям У (называемым игреками Лапласа).

Условия, которым должна удовлетворять функция / (0, X) для того, чтобы разложение ее в бесконечный ряд сферических функций было возможно, для всех функций, с которыми приходится иметь дело на практике, обычно выполняются. Поэтому полученный результат (111.25) можно сформулировать так:

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

72

Выразим функцию Уп (0, А) в явном виде. Для этого следует сферическую функцию Р„ (соз ф), входящую в выражение (111.21), представить как функцию координат 0 и А,

Рп (соз ф) = Рп (соз 0) Рп (соз 0') + 2

(га — к) !

[соз кХ соз кХ' + 31п &Азт кХ"] X

 

 

(п+к)

!

 

 

 

 

 

 

к=1

 

 

 

 

 

 

 

 

х Р „ Д 0 )

 

РпЛП

 

 

 

(III.26)

Доказательства этого соотношения мы не приводим. Формула (111.26)

показывает, что сферическая

функция

степени

п

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

в виде линейной комбинации основных функций

 

 

 

 

Функции вида Рп (СОЗ 0), Рпк Ф)

СОЗ кх

и!

Р пк

(0) 81П

кХ.

 

Рпк(В)созкХ

и

 

 

Рпкф)ёткХ

 

 

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

сферическими функциями.

Функция Рпк (0) назы-

вается присоединенной функцией Лежандра и вычисляется по формуле

Рпк (9) =

8111* 0

й **

п (С°8 ^

.

 

(Ш.27)

 

 

 

 

(йсозб)К

 

 

 

При к = 0 присоединенная функция превращается в полином Лежандра

 

Рп о (0) = Р„(созО),

 

 

 

(111.28)

который называют главной сферической

функцией.

 

 

 

Главная сферическая функция является функцией лишь одной перемен-

ной — полярного расстояния

0 (или широты <р).

Поскольку 0 =

90 —- ф, то

Рп

(соз 0) =

Рп

(81П ф).

 

 

 

Присоединенные сферические

функции

есть

функции двух

переменных:

0 (или ф) и X. Таким образом,

на основании

(111.26) можно сделать вывод, что

сферическая функция степени п содержит одну главную сферическую функцию — полином Лежандра Рп (соз 0) и 2п присоединенных сферических функций. Подставляя (111.26) в (111.21), найдем явное выражение для произвольной

сферической функции Уп

(0, X), как функции координат 0 и А

 

У п ф , а,) = 2

пксо8кХ

+

Впк8ткХ)Рпк(Щ,

(111.29)

где

 

к=о

 

 

 

 

 

 

2ге + 1

 

 

 

 

 

 

А„

Д / ( 0 ' ,

А') Рп

(соз0')Ло

 

 

 

 

 

 

 

 

 

2ге+1

 

 

В т =

0

 

 

А пк '-

(п — к) !

X") соз кХ'Рпкф') (1а

(III. 30)

(га +

А)!

 

 

2п+1

(га — к) !

 

 

 

 

В,пк :

{п + к){

- ^ / ( 0 ' ,

Х')зткХ'РПкФ')с1а

 

Подставив (111.29) в (111.25), найдем разложение заданной на поверхности сферы функции в бесконечный ряд сферических функций

/(0, Я) = 2 2 пксовкХ +

ВпкзткХ)Рпкф).

п=о к-О

§16. КЛАССИФИКАЦИЯ СФЕРИЧЕСКИХ ФУНКЦИЙ

Втеории специальных функций доказывается, что на основном интервале переменной соз 0 от —1 до + 1 полином Лежандра Рп (соз 0) имеет п неравных вещественных корней, расположенных симметрично по отношению к соз 0 = 0.

 

 

 

|!м

Этим корням соответствуют п зна-

р г(ж в )

 

 

чений полярного

расстояния

0,

1,0

 

 

симметричных

 

по

отношению

 

 

 

к экватору сферы; таким образом,

0,5

 

 

полином

Лежандра

Рп

(соз 0)

 

 

как бы делит всю сферу на (п +

1)

• Г "

 

зону; он обращается в нуль

на

ш/\торА/

параллелях,

разделяющих

эти

 

 

зоны, а внутри зон принимает

 

1

попеременно положительные и от-

-0,5

рицательные

значения

(рис.

22,

 

90°

125''го'

т и

23, 24). По этой причине поли-

 

Рис. 22

 

 

номы Лежандра

носят

название

 

 

 

зональных

сферических

функций

 

 

 

 

 

 

 

 

 

 

(гармоник).

 

 

 

 

 

Следует

заметить,

что каждая

зональная гармоника симметрична относи-

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

Вдоль произвольной параллели зональные гармоники в среднем не равны нулю.

Обратимся теперь к присоединенным функциям; общее выражение при-

соединенной функции имеет

вид

 

 

 

 

 

 

 

> 9

акРп

(соз 6)

соз к

= 1 2 , 3 ,

. . . ,

п).

 

 

 

 

(й соз 0)

з т

У

 

'

 

 

 

При к — п многочлен

апРп

(соз 9)

обращается

в

величину

постоянную

 

 

(ЙСО8 0)П

 

 

 

 

 

 

и поэтому сферические функции принимают вид

 

 

 

 

 

 

зт"©созпХ и 8т"0зтгаА,.

 

 

 

 

 

Первый множитель

з т "

0 в этих выражениях обращается

в нуль

только

на полюсах сферы (0 =

0° и 0 =

180°); второй — на 2п меридианах, которыми

ограничиваются сферические секторы

(двухугольники),

где

з т

пХ и

соз пХ

принимают попеременно положительные и отрицательные значения (рис. 25, 26). Отметим, что вдоль меридиана гармоника сохраняет постоянный знак, но величина ее меняется с изменением широты; эти гармоники называются секториальными.

74

Ри(СОЗ 0)

В случаях, когда (0 < к < га) многочлен

(с0о8,^ имеет п — к веще-

ственных корней, которым соответствуют п — к значений полярного расстояния 0, ими определяются га — к параллелей, которыми вся сфера делится на га — к + 1 зону, кроме того, каждый из множителей вш кХ и соз кХ обращается в нуль для 2к значений долготы X, т. е. на 2к меридианах, отстоящих друг от друга на 180°/к. Этой сеткой меридианов и параллелей вся сфера делится на сферические четырехугольники (кроме полярных областей, где образуются треугольники); в каждых двух прилежащих четырехугольниках данные

Р/,4 005 М

функции попеременно положительны и отрицательны (рис. 27, 28, 29). В пределах каждого отдельного сферического четырехугольника (1еззега) величина функции меняется и с широтой, и с долготой; эти гармоники получили наименование тессералъных.

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

Таким образом, в системе 2га + 1 основных сферических функций степени га имеется одна зональная, две секториальные и 2га — 2 тессеральные. Все они являются функциями, осциллирующими на сфере. Повышая степень сферических функций, мы как бы облекаем сферу правильной системой постепенно уменьшающихся участков, в которых происходит перемена знака сферических функций. Таким образом, задача разложения данной функции (0, X) в ряд сферических функций есть вопрос наилучшего приближения к заданной совокупности ее значений на сфере путем комбинации осциллирующих функций.

Проиллюстрируем все вышесказанное рядом примеров. На рис. 22, 23, 24, 25, 26, 27, 28 и 29 показан ход изменения некоторых сферических функций

76

Рис. 28 -5"

115"

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

мы видим, что

зональная гармоника второй степени — Р2

(соз 0) на северном

и

южном

полюсах

принимает

максимальные

значения

2 (сов 0) =

+ 1 ] ,

а

на экваторе

имеет

минимум

2 (соз 0) = —х/2]. Всю

поверхность

сферы

эта гармоника

делит на три области: первая область — от северного полюса

до широты ф^- =

35° 20', в которой функция положительная, вторая область —

от широты

фК

=

35° 20' до широты ф8

= 35° 20', где гармоника отрицательна,

и

третья

область — от широты

ф8 =

35° 20' до

южного

полюса, в которой

гармоника

вновь имеет положительные

значения.

 

 

 

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

Секториальные гармоники Р22 соз (см. рис. 25) и Р4 4 соз (см. рис. 26) делят всю поверхность сферы на сферические секторы, в которых знаки поочередно меняются, первая функция делит сферу на четыре сектора, а вторая — на восемь.

Наконец, тессеральная гармоника Р41 соз X (см. рис. 27) делит сферу на два сектора и четыре зоны, функция Р42 соз на четыре сектора и три зоны (см. рис. 28) и, наконец, функция Р43 соз ЗА — на шесть секторов и две зоны (см. рис. 29).

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

78

Полагая п = 2, получим

Р3 (СОЗ 0) = СОЗ вР2 (СОЗ 0)— у Р х (соз 0) = - | соз 0 ( - | соз2 0 -

-

— соз 0 = у соз3 0 — соз 6,

при п — 3 будем иметь

Р4 (соз 0) = ^ соз 0Р3 (соз 0) - - | Р2 (соз 0) = 1 соз 0 ( у соз3 0 — соз в ) —

- т ( Т с о з 2 0 - У ) = - Х с о з 4 0 - Х с ° 8 2 0 + 4 -

Продолжая этот процесс далее, можно получить полином Лежандра любой

степени

п.

 

 

 

 

 

 

 

 

 

 

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

ненные функции Лежандра в следующем порядке:

 

 

 

 

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

аргументу соз 0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а

1

п

 

 

Рп0

 

 

 

к

рпк

 

р пп

 

0

 

 

1

 

 

 

 

 

 

 

 

1

 

СОЗ 0

 

 

 

 

 

3111 0

 

2

 

3 С032 0

1

 

 

1

3 81П0СО8 0

 

3 ЗЩ2 0

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

3

 

-Ц- С083

0

СОЗ 0

 

1

- | - 8 1 П 0 ( 5 С О 8 2 0

— 1 )

15 81П 3 0

 

 

 

 

 

 

 

2

15 з т 2 0 С08 0

 

 

4

35

СОЗ4 Э _ ^ С О 3 2 0 + А

 

1

- | - З Ш 0 ( 7 СОЗ® 0 — 3

с о з 9)

105 81114 0

— -

 

 

о

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

4^-81112 0 (7С032 9 — 1)

 

 

 

 

 

 

 

 

 

3

105 зш3 0 соз 0

 

 

5 ^

 

1 0

ч а

,

5

г Л

1

У Я П Е ^СОЗ4 0 - | - - С О 3 2 0 + ^ - )

945 з т 5

0

^СО8 5 0 —

СО830 +

"2^- СОЗ 0 )

 

 

 

 

 

 

 

2

31П2 0 ^ СОЗ3 0

• СОЗ 0 ^

 

 

 

 

 

 

 

 

 

3

^ В 1 П З Е ^ 0 0 8 2 0 — 1 . )

 

 

 

 

 

 

 

 

 

4

945 81Ц4 0 соз 0

 

 

79

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]