shimbirev_b_p_teoriya_figury_zemli
.pdfФормула (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')Ло |
|
|||
|
|
4я |
|
|
|
|
|
|
|
2ге+1 |
|
|
В т = |
0 |
|
|
|
А пк '- |
(п — к) ! |
X") соз кХ'Рпкф') (1а |
(III. 30) |
|||||
2я |
(га + |
А)! |
||||||
|
||||||||
|
2п+1 |
(га — к) ! |
|
|
|
|
||
В,пк : |
2я |
{п + к){ |
- ^ / ( 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 соз 2Х (см. рис. 25) и Р4 4 соз 4Х (см. рис. 26) делят всю поверхность сферы на сферические секторы, в которых знаки поочередно меняются, первая функция делит сферу на четыре сектора, а вторая — на восемь.
Наконец, тессеральная гармоника Р41 соз X (см. рис. 27) делит сферу на два сектора и четыре зоны, функция Р42 соз 2Х на четыре сектора и три зоны (см. рис. 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