Вычисление радиальных функций Матье-Ханкеля

Вычисление радиальных функций матье-ханкеля

Н.И. Волвенко, V курс, Институт математики и компьютерных наук ДВГУ, Т.В. Пак тАУ научный руководитель, доцент, к.ф.-м.н., и.о. зав. кафедрой КТ

Функции Матье, в отличие от широко известных специальных функций, таких как полиномы Лежандра, функции Бесселя и Неймана, изучены ещё недостаточно полно. Почти все используемые методы расчёта связаны с разложением в ряды по более простым цилиндрическим и т.п. функциям. Недостаток таких методов в том, что они достаточно громоздки и имеют ограниченную применимость.

Функции Матье возникают при разделении переменных в уравнении Гельмгольца:

, (1)

где - некоторая вещественная положительная константа и - оператор Лапласа.

Эллиптические координаты , допускающие разделение переменных связаны с декартовыми: , .

Полагая Вав методе разделения переменных, получаем уравнения:

, ,

где Ва- константа разделения. Эти уравнения являются вариантами уравнений Матье.

Дифференциальное уравнения Матье имеет вид

, (2)

где обычно переменная Ваимеет вещественное значение, а Ва- заданный вещественный ненулевой параметр.

Собственные значения Ваи граничные условия

Ва(3)

соответствуют чётным функциям Матье , а собственные значения Ваи граничные условия

Ва(4)

нечётным функциям Матье

В силу свойств симметрии уравнение (2) имеет 4 типа периодических решений, называемых функциями Матье 1-ого рода: чётную π-периодическую, чётную 2π-периодическую, нечётную 2π-периодическую, нечётную π-периодическую функции, которые чаще всего обозначаются таким образом: , , , .

Собственные значения , отвечающие функциям , , , , обозначаются через , , , .

Модифицированное уравнение Матье

Ва(5)

получается из уравнения Матье (2) подстановкой . В зависимости от того, будет в (5) или , это уравнение имеет либо решение , либо решение , которые являются соответственно чётной и нечётной функциями от ξ.

Функции, являющиеся решениями уравнения (5), называются радиальными функциями Матье (РФМ).

Различают РФМ 1, 2, 3 и 4 рода: , , , .

Вычисление функций Матье I рода

Радиальные функции Матье первого рода являются решениями ОДУ второго порядка

, Ва(6)

удовлетворяющие в нуле условию

, если Ва(7)

, если

И на бесконечности условию

~, Ва(8)

где - задано, а Ва() - собственные значения задачи (2), (3), (4),

Параметр Ваиспользуются для различия случаев использования чётного или нечётного номера собственного значения для π и 2π периодических собственных функций:

Для решения задачи (6)-(8) используем модификацию метода фазовых функций.

Введём замену переменных:

Ва(9)

Ва(10)

Здесь Ва- "масштабирующая" функция, положительная на , удовлетворяющая условию Вапри , её выбор находится в нашем распоряжении.

Подставляя (9), (10) в исходное уравнение (6) задачи для Ваи :

Ва(11)

Ва(12)

где Ваи .

Для совместного решения задач Коши для Ваи Ваиспользуется следующий приём. Функцию ищем в точках . На каждом из отрезков Вавспомогательные функции Ванаходятся, как решение задач Коши

Ва(13)

где .

Поскольку для любых решений Ваи , уравнений (12) и (13) справедливо соотношение , получаем рекуррентные формулы ВлназадВ» для вычисления , ,

, , (14)

причём .

Итак, краткий алгоритм решения задачи (6)-(8) состоит в следующем:

1. Решаются совместно задачи Коши (11), (12) запоминая в точках разбиения отрезка Вавеличины , , ;

2. Полагая , по формуле (14) вычисляем , ;

3. По формуле (10) вычисляем функции , ;

4. Из (9) и (10) получаем выражение для производной функции

.

В качестве сглаживающей функции предлагается следующая функция

, где .

Вычисление функций Матье III рода

Волновая радиальная функция Матье-Ханкеля третьего рода является решением обыкновенного дифференциального уравнения второго ворядка на полубесконечном интервале:

, . (15)

Условие на бесконечности

~, . (16)

Для уравнения (15) условие (16) эквивалентно условию:

,

и при достаточно больших Валинейному соотношению:

, .

Ва(17)

Решение задачи (17) существует, единственно и при достаточно больших Вапредставимо асимптотическим рядом .

Рассмотрим алгоритм нахождения функций . Для их вычисления нужно перенести граничное условие

,

где , справа налево от точки Вадо точки .

Воспользуемся вариантом ортогональной дифференциальной прогонки.

По всему отрезку Вапереносим соотношение

,

потребовав выполнение условия для всех , , где Ваи Ваудовлетворяют системе дифференциальных уравнений 1-ого порядка

.

Функции Матье 3-его рода ищем по формуле:

,

где .

Функции Матье 2-ого рода вычисляются по формуле:

.

функция матье дифференциальное уравнение

Описанные алгоритмы вычисления радиальных функций эллиптического цилиндра опробованы в широком диапазоне изменения параметров. Точность результатов определяется точностью используемого метода Рунге-Кутта для решения соответствующих задач Коши.

Литература

1. Абрамов А.А., Дышко А.Л., Пак Т.В. и др. Численные методы решения задач на собственные значения для систем обыкновенных дифференциальных уравнений с особенностями. тАУ Третья конференция по дифференциальным уравнениям и приложениям. тАУ Тезисы докладов. Руссе, Болгария, 1985. тАУ с.4.

2. Миллер У. мл. Симметрия и разделение переменных / Пер. с англ. тАУ М.: Мир, 1981. тАУ 342 с.

3. Справочник по специальным функциям с формулами, графиками таблицами. / Под редакцией М. Абрамовица, И. Стигана. тАУ М. тАУ 1979. тАУ 832 с.:ил.

Вместе с этим смотрят:

РЖнварiантнi пiдпростори. Власнi вектори i власнi значення лiнiйного оператора

Актуальные проблемы квантовой механики

Алгебра и алгебраические системы

Анализ эмпирического распределения

Аналитическая теория чисел. L-функция Дирихле