Функция Матье

редактировать
Не путать с функцией Массьё.

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

d 2 у d Икс 2 + ( а - 2 q потому что 2 Икс ) у знак равно 0 , {\ displaystyle {\ frac {d ^ {2} y} {dx ^ {2}}} + (a-2q \ cos 2x) y = 0,}

где и - параметры. Впервые они были представлены Эмилем Леонаром Матье, который столкнулся с ними во время изучения вибрирующих эллиптических пластиков. У них есть приложения во многих областях физических наук, таких как оптика, квантовая механика и общая теория относительности. Они обычно возникают в задачах, связанных с периодическим движением, или при анализе краевых задач уравнений в частных производных, обладающих эллиптической симметрией. а {\ displaystyle a} q {\ displaystyle q}

СОДЕРЖАНИЕ
  • 1 Определение
    • 1.1 Функции Матье
    • 1.2 Модифицированные функции Матье
    • 1.3 Нормализация
  • 2 Теория Флоке
  • 3 Другие типы функций Матье
    • 3.1 Второй вид
    • 3.2 Дробный порядок
  • 4 Явное представление и вычисление
    • 4.1 Первый вид
    • 4.2 Второй вид
    • 4.3 Измененные функции
  • 5 Недвижимость
    • 5.1 Качественное поведение
    • 5.2 Размышления и переводы
    • 5.3 Ортогональность и полнота
    • 5.4 Интегральные тождества
    • 5.5 Асимптотические разложения
  • 6 приложений
    • 6.1 Уравнения с частными производными
    • 6.2 Динамические проблемы
    • 6.3 Квантовая механика
  • 7 См. Также
  • 8 Примечания
  • 9 ссылки
  • 10 Внешние ссылки
Определение

Функции Матье

В некоторых случаях функция Матье относится к решениям дифференциального уравнения Матье для произвольных значений и. Когда не может возникнуть путаницы, другие авторы используют этот термин для обозначения - или -периодических решений, которые существуют только для особых значений и. Точнее, для заданных (реальных) таких периодических решений существует бесконечное число значений, называемых характеристическими числами, условно индексируемых как две отдельные последовательности и, для. Соответствующие функции обозначены и соответственно. Иногда их также называют косинус-эллиптическими и синус-эллиптическими, или функциями Матье первого рода. а {\ displaystyle a} q {\ displaystyle q} π {\ displaystyle \ pi} 2 π {\ displaystyle 2 \ pi} а {\ displaystyle a} q {\ displaystyle q} q {\ displaystyle q} а {\ displaystyle a} а п ( q ) {\ displaystyle a_ {n} (q)} б п ( q ) {\ displaystyle b_ {n} (q)} п знак равно 1 , 2 , 3 , {\ Displaystyle п = 1,2,3, \ ldots} ce п ( Икс , q ) {\ displaystyle {\ text {ce}} _ {n} (x, q)} se п ( Икс , q ) {\ displaystyle {\ text {se}} _ {n} (x, q)}

В результате предположения, что это реально, как характеристические числа, так и связанные функции являются действительными. q {\ displaystyle q}

ce п ( Икс , q ) {\ displaystyle {\ text {ce}} _ {n} (x, q)}и могут быть дополнительно классифицированы по четности и периодичности (как в отношении) следующим образом: se п ( Икс , q ) {\ displaystyle {\ text {se}} _ {n} (x, q)} Икс {\ displaystyle x}

Функция Паритет Период
ce п ,   п  четный {\ displaystyle {\ text {ce}} _ {n}, \ n {\ text {even}}} четный π {\ displaystyle \ pi}
ce п ,   п  странный {\ displaystyle {\ text {ce}} _ {n}, \ n {\ text {odd}}} четный 2 π {\ displaystyle 2 \ pi}
se п ,   п  четный {\ displaystyle {\ text {se}} _ {n}, \ n {\ text {even}}} странный π {\ displaystyle \ pi}
se п ,   п  странный {\ displaystyle {\ text {se}} _ {n}, \ n {\ text {odd}}} странный 2 π {\ displaystyle 2 \ pi}

Индексация с целым числом, кроме того, служит для организации характеристических чисел в порядке возрастания, удобна в том, что и станет пропорционально и, как. Поскольку это целое число, это приводит к классификации и как функций Матье (первого рода) целого порядка. Для общих и, кроме них, могут быть определены решения, включая функции Матье дробного порядка, а также непериодические решения. п {\ displaystyle n} ce п ( Икс , q ) {\ displaystyle {\ text {ce}} _ {n} (x, q)} se п ( Икс , q ) {\ displaystyle {\ text {se}} _ {n} (x, q)} потому что п Икс {\ displaystyle \ cos nx} грех п Икс {\ Displaystyle \ грех nx} q 0 {\ displaystyle q \ rightarrow 0} п {\ displaystyle n} ce п {\ displaystyle {\ text {ce}} _ {n}} se п {\ displaystyle {\ text {se}} _ {n}} а {\ displaystyle a} q {\ displaystyle q}

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

Тесно связаны модифицированные функции Матье, также известные как радиальные функции Матье, которые являются решениями модифицированного дифференциального уравнения Матье.

d 2 у d Икс 2 - ( а - 2 q шиш 2 Икс ) у знак равно 0 , {\ displaystyle {\ frac {d ^ {2} y} {dx ^ {2}}} - (a-2q \ cosh 2x) y = 0,}

которое можно связать с исходным уравнением Матье, взяв. Соответственно, модифицированные функции Матье первого рода интегрального порядка, обозначаемые и, определяются из Икс ± я Икс {\ Displaystyle х \ к \ pm {\ rm {я}} х} Ce п ( Икс , q ) {\ displaystyle {\ text {Ce}} _ {n} (x, q)} Se п ( Икс , q ) {\ displaystyle {\ text {Se}} _ {n} (x, q)}

Ce п ( Икс , q ) знак равно ce п ( я Икс , q ) . Se п ( Икс , q ) знак равно - я se п ( я Икс , q ) . {\ displaystyle {\ begin {align} {\ text {Ce}} _ {n} (x, q) amp; = {\ text {ce}} _ {n} ({\ rm {i}} x, q). \\ {\ text {Se}} _ {n} (x, q) amp; = - {\ rm {i}} \, {\ text {se}} _ {n} ({\ rm {i}} x, q). \ end {выравнивается}}}

Эти функции являются действительными, когда они реальны. Икс {\ displaystyle x}

Нормализация

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

0 2 π ce п ( Икс , q ) 2 d Икс знак равно 0 2 π se п ( Икс , q ) 2 d Икс знак равно π {\ displaystyle \ int _ {0} ^ {2 \ pi} {\ text {ce}} _ {n} (x, q) ^ {2} dx = \ int _ {0} ^ {2 \ pi} { \ text {se}} _ {n} (x, q) ^ {2} dx = \ pi}

а также требовать и как. ce п ( Икс , q ) + потому что п Икс {\ displaystyle {\ text {ce}} _ {n} (x, q) \ rightarrow + \ cos nx} se п ( Икс , q ) + грех п Икс {\ displaystyle {\ text {se}} _ {n} (x, q) \ rightarrow + \ sin nx} q 0 {\ displaystyle q \ rightarrow 0}

Теория Флоке
Основная статья: теория Флоке

Многие свойства дифференциального уравнения Матье можно вывести из общей теории обыкновенных дифференциальных уравнений с периодическими коэффициентами, называемой теорией Флоке. Центральный результат - теорема Флоке:

Теорема Флоке - уравнение Матье всегда имеет хотя бы одно решение, такое что, где - постоянная, которая зависит от параметров уравнения и может быть действительной или комплексной. у ( Икс ) {\ Displaystyle у (х)} у ( Икс + π ) знак равно σ у ( Икс ) {\ Displaystyle у (х + \ пи) = \ сигма у (х)} σ {\ displaystyle \ sigma}

Характерные числа естественно связать с теми значениями, которые в результате. Однако следует отметить, что теорема гарантирует только существование по крайней мере, одно решение, удовлетворяющее, когда уравнение Матье на самом деле имеет два независимых решения для любого заданного,. Действительно, оказалось, что с равным одной из характеристических чисел, уравнение Матье имеет только одно периодическое решение (то есть, с периодом или), и это решение является одним из,. Другое решение является непериодическим, обозначается и, соответственно, и называется функцией Матье второго рода. Этот результат формально можно сформулировать как теорему Инса: а ( q ) {\ Displaystyle а (д)} а {\ displaystyle a} σ знак равно ± 1 {\ displaystyle \ sigma = \ pm 1} у ( Икс + π ) знак равно σ у ( Икс ) {\ Displaystyle у (х + \ пи) = \ сигма у (х)} а {\ displaystyle a} q {\ displaystyle q} а {\ displaystyle a} π {\ displaystyle \ pi} 2 π {\ displaystyle 2 \ pi} ce п ( Икс , q ) {\ displaystyle {\ text {ce}} _ {n} (x, q)} se п ( Икс , q ) {\ displaystyle {\ text {se}} _ {n} (x, q)} fe п ( Икс , q ) {\ displaystyle {\ text {fe}} _ {n} (x, q)} ge п ( Икс , q ) {\ displaystyle {\ text {ge}} _ {n} (x, q)}

Теорема Инса - Определите в основном периодическую функцию как удовлетворяющую. Тогда, за исключением тривиального случая, уравнение Матье никогда не имеет двух (независимых) принципиально периодических решений при одинаковых значениях и. у ( Икс + π ) знак равно ± у ( Икс ) {\ Displaystyle у (х + \ пи) = \ пм у (х)} q знак равно 0 {\ displaystyle q = 0} а {\ displaystyle a} q {\ displaystyle q} Пример из теоремы Флоке, с,, (реальная часть, красный, мнимая часть, зеленый) п ( а , q , Икс ) {\ Displaystyle Р (а, д, х)} а знак равно 1 {\ displaystyle a = 1} q знак равно 1 / 5 {\ displaystyle q = 1/5} μ 1 + 0,0995 я {\ displaystyle \ mu \ приблизительно 1 + 0,0995i}

Эквивалентное утверждение теоремы Флоке состоит в том, что уравнение Матье допускает комплексное решение вида

F ( а , q , Икс ) знак равно exp ( я μ Икс ) п ( а , q , Икс ) , {\ Displaystyle F (а, д, х) = \ ехр (я \ му \, х) \, Р (а, д, х),}

где - комплексное число, показатель Флоке (или иногда показатель Матье), а - комплексная функция, периодическая по с периодом. Пример показан справа. μ {\ displaystyle \ mu} п {\ displaystyle P} Икс {\ displaystyle x} π {\ displaystyle \ pi} п ( а , q , Икс ) {\ Displaystyle Р (а, д, х)}

Другие типы функций Матье

Второй вид

Поскольку уравнение Матье является дифференциальным уравнением второго порядка, можно построить два линейно независимых решения. Теория Флоке гласит, что если равно характеристическому числу, одно из этих решений можно считать периодическим, а другое - непериодическим. Периодическое решение является одной из функций и называется функцией Матье первого рода целого порядка. Непериодическая обозначается либо и, соответственно, и называется функцией Матье второго рода (целого порядка). Непериодические решения неустойчивы, т. Е. Расходятся как. а {\ displaystyle a} ce п ( Икс , q ) {\ displaystyle {\ text {ce}} _ {n} (x, q)} se п ( Икс , q ) {\ displaystyle {\ text {se}} _ {n} (x, q)} fe п ( Икс , q ) {\ displaystyle {\ text {fe}} _ {n} (x, q)} ge п ( Икс , q ) {\ displaystyle {\ text {ge}} _ {n} (x, q)} z ± {\ Displaystyle г \ rightarrow \ pm \ infty}

Вторые решения соответствуют модифицированным функциям Матье и естественно определяются как и. Ce п ( Икс , q ) {\ displaystyle {\ text {Ce}} _ {n} (x, q)} Se п ( Икс , q ) {\ displaystyle {\ text {Se}} _ {n} (x, q)} Fe п ( Икс , q ) знак равно - я fe п ( Икс я , q ) {\ displaystyle {\ text {Fe}} _ {n} (x, q) = - я {\ text {fe}} _ {n} (xi, q)} Ge п ( Икс , q ) знак равно ge п ( Икс я , q ) {\ displaystyle {\ text {Ge}} _ {n} (x, q) = {\ text {ge}} _ {n} (xi, q)}

Дробный порядок

Функции Матье дробного порядка можно определить как те решения и, нецелое число, которые превращаются в и как. Если иррационально, они непериодичны; однако они остаются ограниченными как. ce п ( Икс , q ) {\ displaystyle {\ text {ce}} _ {p} (x, q)} se п ( Икс , q ) {\ displaystyle {\ text {se}} _ {p} (x, q)} п {\ displaystyle p} потому что п Икс {\ displaystyle \ cos px} грех п Икс {\ displaystyle \ sin px} q 0 {\ displaystyle q \ rightarrow 0} п {\ displaystyle p} Икс {\ Displaystyle х \ rightarrow \ infty}

Важное свойство решений и, для нецелых, является то, что они существуют для того же значения. Напротив, when является целым числом и никогда не встречается для одного и того же значения. (См. Теорему Инса выше.) ce п ( Икс , q ) {\ displaystyle {\ text {ce}} _ {p} (x, q)} se п ( Икс , q ) {\ displaystyle {\ text {se}} _ {p} (x, q)} п {\ displaystyle p} а {\ displaystyle a} п {\ displaystyle p} ce п ( Икс , q ) {\ displaystyle {\ text {ce}} _ {p} (x, q)} se п ( Икс , q ) {\ displaystyle {\ text {se}} _ {p} (x, q)} а {\ displaystyle a}

Эти классификации кратко изложены в таблице ниже. Аналогично определяются модифицированные аналоги функции Матье.

Классификация функций Матье
Заказ Первый вид Второй вид
интеграл ce п ( Икс , q ) {\ displaystyle {\ text {ce}} _ {n} (x, q)} fe п ( Икс , q ) {\ displaystyle {\ text {fe}} _ {n} (x, q)}
интеграл se п ( Икс , q ) {\ displaystyle {\ text {se}} _ {n} (x, q)} ge п ( Икс , q ) {\ displaystyle {\ text {ge}} _ {n} (x, q)}
Дробное

( нецелое) п {\ displaystyle p}

ce п ( Икс , q ) {\ displaystyle {\ text {ce}} _ {p} (x, q)} se п ( Икс , q ) {\ displaystyle {\ text {se}} _ {p} (x, q)}
Явное представление и вычисление

Первый вид

Функции Матье первого рода можно представить в виде ряда Фурье :

ce 2 п ( Икс , q ) знак равно р знак равно 0 А 2 р ( 2 п ) ( q ) потому что ( 2 р Икс ) ce 2 п + 1 ( Икс , q ) знак равно р знак равно 0 А 2 р + 1 ( 2 п + 1 ) ( q ) потому что [ ( 2 р + 1 ) Икс ] se 2 п + 1 ( Икс , q ) знак равно р знак равно 0 B 2 р + 1 ( 2 п + 1 ) ( q ) грех [ ( 2 р + 1 ) Икс ] se 2 п + 2 ( Икс , q ) знак равно р знак равно 0 B 2 р + 2 ( 2 п + 2 ) ( q ) грех [ ( 2 р + 2 ) Икс ] {\ displaystyle {\ begin {align} {\ text {ce}} _ {2n} (x, q) amp; = \ sum _ {r = 0} ^ {\ infty} A_ {2r} ^ {(2n)} (q) \ cos (2rx) \\ {\ text {ce}} _ {2n + 1} (x, q) amp; = \ sum _ {r = 0} ^ {\ infty} A_ {2r + 1} ^ {(2n + 1)} (q) \ cos \ left [(2r + 1) x \ right] \\ {\ text {se}} _ {2n + 1} (x, q) amp; = \ sum _ { r = 0} ^ {\ infty} B_ {2r + 1} ^ {(2n + 1)} (q) \ sin \ left [(2r + 1) x \ right] \\ {\ text {se}} _ {2n + 2} (x, q) amp; = \ sum _ {r = 0} ^ {\ infty} B_ {2r + 2} ^ {(2n + 2)} (q) \ sin \ left [(2r + 2) x \ right] \\\ конец {выровнено}}}

Коэффициенты разложения и являются функциями, но не зависят от. Подстановкой в ​​уравнение Матье можно показать, что они подчиняются трехчленным рекуррентным соотношениям в нижнем индексе. Например, для каждого находит А j ( я ) ( q ) {\ Displaystyle A_ {j} ^ {(я)} (q)} B j ( я ) ( q ) {\ Displaystyle B_ {j} ^ {(я)} (q)} q {\ displaystyle q} Икс {\ displaystyle x} ce 2 п {\ displaystyle {\ text {ce}} _ {2n}}

а А 0 - q А 2 знак равно 0 ( а - 4 ) А 2 - q ( А 4 + 2 А 0 ) знак равно 0 ( а - 4 р 2 ) А 2 р - q ( А 2 р + 2 + А 2 р - 2 ) знак равно 0 , р 2 {\ displaystyle {\ begin {align} aA_ {0} -qA_ {2} amp; = 0 \\ (a-4) A_ {2} -q (A_ {4} + 2A_ {0}) amp; = 0 \\ (a-4r ^ {2}) A_ {2r} -q (A_ {2r + 2} + A_ {2r-2}) amp; = 0, \ quad r \ geq 2 \ end {align}}}

Будучи повторения второго порядка в индексе, всегда можно найти два независимых решения, и таким образом, что общее решение может быть выражено в виде линейной комбинации двух:. Более того, в этом частном случае асимптотический анализ показывает, что один из возможных вариантов выбора фундаментальных решений обладает свойством 2 р {\ displaystyle 2r} Икс 2 р {\ displaystyle X_ {2r}} Y 2 р {\ displaystyle Y_ {2r}} А 2 р знак равно c 1 Икс 2 р + c 2 Y 2 р {\ displaystyle A_ {2r} = c_ {1} X_ {2r} + c_ {2} Y_ {2r}}

Икс 2 р знак равно р - 2 р - 1 ( - е 2 q 4 ) р [ 1 + О ( р - 1 ) ] Y 2 р знак равно р 2 р - 1 ( - 4 е 2 q ) р [ 1 + О ( р - 1 ) ] {\ displaystyle {\ begin {align} X_ {2r} amp; = r ^ {- 2r-1} \ left (- {\ frac {e ^ {2} q} {4}} \ right) ^ {r} \ left [1 + {\ mathcal {O}} (r ^ {- 1}) \ right] \\ Y_ {2r} amp; = r ^ {2r-1} \ left (- {\ frac {4} {e ^ {2} q}} \ right) ^ {r} \ left [1 + {\ mathcal {O}} (r ^ {- 1}) \ right] \ end {выравнивается}}}

В частности, конечно, тогда как расходится. Таким образом, записывая, мы видим, что для сходимости представления в виде ряда Фурье необходимо выбрать такое, что. Эти варианты соответствуют характеристическим номерам. Икс 2 р {\ displaystyle X_ {2r}} Y 2 р {\ displaystyle Y_ {2r}} А 2 р знак равно c 1 Икс 2 р + c 2 Y 2 р {\ displaystyle A_ {2r} = c_ {1} X_ {2r} + c_ {2} Y_ {2r}} ce 2 п {\ displaystyle {\ text {ce}} _ {2n}} а {\ displaystyle a} c 2 знак равно 0 {\ displaystyle c_ {2} = 0} а {\ displaystyle a}

В общем, однако, решение трехчленной рекуррентности с переменными коэффициентами не может быть представлено простым способом, и, следовательно, нет простого способа определить из условия. Более того, даже если приблизительное значение характеристического числа известно, его нельзя использовать для получения коэффициентов путем численного повторения повторяемости в сторону увеличения. Причина в том, что до тех пор, пока характеристическое число приближается только к характеристическому числу, оно не является идентичным, и расходящееся решение в конечном итоге доминирует для достаточно больших. а {\ displaystyle a} c 2 знак равно 0 {\ displaystyle c_ {2} = 0} А 2 р {\ displaystyle A_ {2r}} р {\ displaystyle r} а {\ displaystyle a} c 2 {\ displaystyle c_ {2}} 0 {\ displaystyle 0} Y 2 р {\ displaystyle Y_ {2r}} р {\ displaystyle r}

Чтобы преодолеть эти проблемы, требуются более сложные полуаналитические / численные подходы, например, с использованием расширения непрерывной дроби, преобразованием рекуррентности в проблему собственных значений матрицы или реализацией алгоритма обратной рекуррентности. Сложность трехчленного рекуррентного отношения - одна из причин, по которой существует мало простых формул и тождеств, включающих функции Матье.

На практике функции Матье и соответствующие характеристические числа могут быть вычислены с использованием предварительно упакованного программного обеспечения, такого как Mathematica, Maple, MATLAB и SciPy. Для малых значений и низкого порядка они также могут быть выражены пертурбативно в виде степенных рядов, что может быть полезно в физических приложениях. q {\ displaystyle q} п {\ displaystyle n} q {\ displaystyle q}

Второй вид

Есть несколько способов представления функций Матье второго рода. Одно представление в терминах функций Бесселя :

fe 2 п ( Икс , q ) знак равно - π γ п 2 р знак равно 0 ( - 1 ) р + п А 2 р ( 2 п ) ( - q )   Я [ J р ( q е я Икс ) Y р ( q е - я Икс ) ] , где  γ п знак равно { 2 ,  если  п знак равно 0 2 п ,  если  п 1 fe 2 п + 1 ( Икс , q ) знак равно π q 2 р знак равно 0 ( - 1 ) р + п А 2 р + 1 ( 2 п + 1 ) ( - q )   Я [ J р ( q е я Икс ) Y р + 1 ( q е - я Икс ) + J р + 1 ( q е я Икс ) Y р ( q е - я Икс ) ] ge 2 п + 1 ( Икс , q ) знак равно - π q 2 р знак равно 0 ( - 1 ) р + п B 2 р + 1 ( 2 п + 1 ) ( - q )   Re [ J р ( q е я Икс ) Y р + 1 ( q е - я Икс ) - J р + 1 ( q е я Икс ) Y р ( q е - я Икс ) ] ge 2 п + 2 ( Икс , q ) знак равно - π q 4 ( п + 1 ) р знак равно 0 ( - 1 ) р + п B 2 р + 2 ( 2 п + 2 ) ( - q )   Re [ J р ( q е я Икс ) Y р + 2 ( q е - я Икс ) - J р + 2 ( q е я Икс ) Y р ( q е - я Икс ) ] {\ displaystyle {\ begin {align} {\ text {fe}} _ {2n} (x, q) amp; = - {\ frac {\ pi \ gamma _ {n}} {2}} \ sum _ {r = 0} ^ {\ infty} (- 1) ^ {r + n} A_ {2r} ^ {(2n)} (- q) \ {\ text {Im}} [J_ {r} ({\ sqrt { q}} e ^ {ix}) Y_ {r} ({\ sqrt {q}} e ^ {- ix})], \ quad {\ text {where}} \ gamma _ {n} = \ left \ { {\ begin {array} {cc} {\ sqrt {2}}, amp; {\ text {if}} n = 0 \\ 2n, amp; {\ text {if}} n \ geq 1 \ end {array}} \ right. \\ {\ text {fe}} _ {2n + 1} (x, q) amp; = {\ frac {\ pi {\ sqrt {q}}} {2}} \ sum _ {r = 0 } ^ {\ infty} (- 1) ^ {r + n} A_ {2r + 1} ^ {(2n + 1)} (- q) \ {\ text {Im}} [J_ {r} ({\ sqrt {q}} e ^ {ix}) Y_ {r + 1} ({\ sqrt {q}} e ^ {- ix}) + J_ {r + 1} ({\ sqrt {q}} e ^ { ix}) Y_ {r} ({\ sqrt {q}} e ^ {- ix})] \\ {\ text {ge}} _ {2n + 1} (x, q) amp; = - {\ frac { \ pi {\ sqrt {q}}} {2}} \ sum _ {r = 0} ^ {\ infty} (- 1) ^ {r + n} B_ {2r + 1} ^ {(2n + 1) } (- q) \ {\ text {Re}} [J_ {r} ({\ sqrt {q}} e ^ {ix}) Y_ {r + 1} ({\ sqrt {q}} e ^ {- ix}) - J_ {r + 1} ({\ sqrt {q}} e ^ {ix}) Y_ {r} ({\ sqrt {q}} e ^ {- ix})] \\ {\ text { ge}} _ {2n + 2} (x, q) amp; = - {\ frac {\ pi q} {4 (n + 1)}} \ sum _ {r = 0} ^ {\ infty} (- 1) ^ {r + n} B_ {2r + 2} ^ {(2n + 2)} (- q) \ {\ text {Re}} [J_ {r} ({\ sqrt {q}} e ^ {ix }) Y_ {r + 2} ({\ sqrt {q}} e ^ {- ix}) - J_ {r + 2} ({\ sqrt {q}} e ^ {ix}) Y_ {r} ({ \ sqrt {q}} e ^ {- ix})] \ конец {выровнено}}}

где и и функция Бесселя первого и второго рода. п , q gt; 0 {\ displaystyle n, qgt; 0} J р ( Икс ) {\ Displaystyle J_ {r} (х)} Y р ( Икс ) {\ Displaystyle Y_ {r} (х)}

Измененные функции

Традиционный подход к числовому вычислению модифицированных функций Матье заключается в использовании ряда произведений функций Бесселя. Для больших и необходимо тщательно выбирать форму ряда, чтобы избежать ошибок вычитания. п {\ displaystyle n} q {\ displaystyle q}

Характеристики

Существует относительно немного аналитических выражений и тождеств, включающих функции Матье. Более того, в отличие от многих других специальных функций, решения уравнения Матье, вообще говоря, не могут быть выражены через гипергеометрические функции. В этом можно убедиться, преобразовав уравнение Матье в алгебраическую форму с помощью замены переменной: т знак равно потому что ( Икс ) {\ Displaystyle т = \ соз (х)}

( 1 - т 2 ) d 2 у d т 2 - т d у d т + ( а + 2 q ( 1 - 2 т 2 ) ) у знак равно 0. {\ displaystyle (1-t ^ {2}) {\ frac {d ^ {2} y} {dt ^ {2}}} - t \, {\ frac {dy} {dt}} + (a + 2q (1-2t ^ {2})) \, y = 0.}

Поскольку это уравнение имеет нерегулярную особую точку на бесконечности, его нельзя преобразовать в уравнение гипергеометрического типа.

Качественное поведение

Примерные графики функций Матье первого рода Участок под разную ce 1 ( Икс , q ) {\ displaystyle {\ text {ce}} _ {1} (x, q)} q {\ displaystyle q}

Для малых, и ведут себя так же, как и. Для произвольных они могут значительно отличаться от своих тригонометрических аналогов; однако в целом они остаются периодическими. Кроме того, для какой - либо реальной, и есть именно простые нули в, и, как нули кластера о. q {\ displaystyle q} ce п {\ displaystyle {\ text {ce}} _ {n}} se п {\ displaystyle {\ text {se}} _ {n}} потому что п Икс {\ displaystyle \ cos nx} грех п Икс {\ Displaystyle \ грех nx} q {\ displaystyle q} q {\ displaystyle q} ce м ( Икс , q ) {\ displaystyle {\ text {ce}} _ {m} (x, q)} se м + 1 ( Икс , q ) {\ displaystyle {\ text {se}} _ {m + 1} (x, q)} м {\ displaystyle m} 0 lt; Икс lt; π {\ Displaystyle 0 lt;х lt;\ pi} q {\ displaystyle q \ rightarrow \ infty} Икс знак равно π / 2 {\ Displaystyle х = \ пи / 2}

Для и поскольку модифицированные функции Матье имеют тенденцию вести себя как периодические функции с затуханием. q gt; 0 {\ displaystyle qgt; 0} Икс {\ Displaystyle х \ rightarrow \ infty}

В последующем, и факторы, от разложения Фурье для и можно ссылаться (см Явного представления и вычисление). Они зависят от и не зависят от. А {\ displaystyle A} B {\ displaystyle B} ce п {\ displaystyle {\ text {ce}} _ {n}} se п {\ displaystyle {\ text {se}} _ {n}} q {\ displaystyle q} п {\ displaystyle n} Икс {\ displaystyle x}

Размышления и переводы

Благодаря их четности и периодичности, а также имеют простые свойства при отражениях и переводах, кратных: ce п {\ displaystyle {\ text {ce}} _ {n}} se п {\ displaystyle {\ text {se}} _ {n}} π {\ displaystyle \ pi}

ce п ( Икс + π ) знак равно ( - 1 ) п ce п ( Икс ) se п ( Икс + π ) знак равно ( - 1 ) п se п ( Икс ) ce п ( Икс + π / 2 ) знак равно ( - 1 ) п ce п ( - Икс + π / 2 ) se п + 1 ( Икс + π / 2 ) знак равно ( - 1 ) п se п + 1 ( - Икс + π / 2 ) {\ displaystyle {\ begin {align} amp; {\ text {ce}} _ {n} (x + \ pi) = (- 1) ^ {n} {\ text {ce}} _ {n} (x) \ \ amp; {\ text {se}} _ {n} (x + \ pi) = (- 1) ^ {n} {\ text {se}} _ {n} (x) \\ amp; {\ text {ce} } _ {n} (x + \ pi / 2) = (- 1) ^ {n} {\ text {ce}} _ {n} (- x + \ pi / 2) \\ amp; {\ text {se}} _ {n + 1} (x + \ pi / 2) = (- 1) ^ {n} {\ text {se}} _ {n + 1} (- x + \ pi / 2) \ end {выровнено}}}

Также можно писать функции с отрицательными значениями через функции с положительными значениями: q {\ displaystyle q} q {\ displaystyle q}

ce 2 п + 1 ( Икс , - q ) знак равно ( - 1 ) п se 2 п + 1 ( - Икс + π / 2 , q ) ce 2 п + 2 ( Икс , - q ) знак равно ( - 1 ) п ce 2 п + 2 ( - Икс + π / 2 , q ) se 2 п + 1 ( Икс , - q ) знак равно ( - 1 ) п ce 2 п + 1 ( - Икс + π / 2 , q ) se 2 п + 2 ( Икс , - q ) знак равно ( - 1 ) п se 2 п + 2 ( - Икс + π / 2 , q ) {\ displaystyle {\ begin {align} amp; {\ text {ce}} _ {2n + 1} (x, -q) = (- 1) ^ {n} {\ text {se}} _ {2n + 1 } (- x + \ pi / 2, q) \\ amp; {\ text {ce}} _ {2n + 2} (x, -q) = (- 1) ^ {n} {\ text {ce}} _ {2n + 2} (- x + \ pi / 2, q) \\ amp; {\ text {se}} _ {2n + 1} (x, -q) = (- 1) ^ {n} {\ text { ce}} _ {2n + 1} (- x + \ pi / 2, q) \\ amp; {\ text {se}} _ {2n + 2} (x, -q) = (- 1) ^ {n} {\ text {se}} _ {2n + 2} (- x + \ pi / 2, q) \ end {align}}}

Более того,

а 2 п + 1 ( q ) знак равно б 2 п + 1 ( - q ) б 2 п + 2 ( q ) знак равно б 2 п + 2 ( - q ) {\ displaystyle {\ begin {align} amp; a_ {2n + 1} (q) = b_ {2n + 1} (- q) \\ amp; b_ {2n + 2} (q) = b_ {2n + 2} (- q) \ конец {выровнено}}}

Ортогональность и полнота

Как и их тригонометрические аналоги и, периодические функции Матье и удовлетворяют соотношениям ортогональности потому что п Икс {\ displaystyle \ cos nx} грех п Икс {\ Displaystyle \ грех nx} ce п ( Икс , q ) {\ displaystyle {\ text {ce}} _ {n} (x, q)} se п ( Икс , q ) {\ displaystyle {\ text {se}} _ {n} (x, q)}

0 2 π ce п ce м d Икс знак равно 0 2 π se п se м d Икс знак равно δ п м π 0 2 π ce п se м d Икс знак равно 0 {\ displaystyle {\ begin {align} amp; \ int _ {0} ^ {2 \ pi} {\ text {ce}} _ {n} {\ text {ce}} _ {m} dx = \ int _ { 0} ^ {2 \ pi} {\ text {se}} _ {n} {\ text {se}} _ {m} dx = \ delta _ {nm} \ pi \\ amp; \ int _ {0} ^ {2 \ pi} {\ text {ce}} _ {n} {\ text {se}} _ {m} dx = 0 \ end {align}}}

Более того, при фиксированном и рассматриваемом как собственное значение уравнение Матье имеет форму Штурма-Лиувилля. Это означает, что собственные функции и образуют полный набор, т.е. любую - или -периодическую функцию можно разложить в ряд по и. q {\ displaystyle q} а {\ displaystyle a} ce п ( Икс , q ) {\ displaystyle {\ text {ce}} _ {n} (x, q)} se п ( Икс , q ) {\ displaystyle {\ text {se}} _ {n} (x, q)} π {\ displaystyle \ pi} 2 π {\ displaystyle 2 \ pi} Икс {\ displaystyle x} ce п ( Икс , q ) {\ displaystyle {\ text {ce}} _ {n} (x, q)} se п ( Икс , q ) {\ displaystyle {\ text {se}} _ {n} (x, q)}

Интегральные тождества

Решения уравнения Матье удовлетворяют классу интегральных тождеств относительно ядер, которые являются решениями χ ( Икс , Икс ) {\ Displaystyle \ чи (х, х ')}

2 χ Икс 2 - 2 χ Икс 2 знак равно 2 q ( потому что 2 Икс - потому что 2 Икс ) χ {\ displaystyle {\ frac {\ partial ^ {2} \ chi} {\ partial x ^ {2}}} - {\ frac {\ partial ^ {2} \ chi} {\ partial x '^ {2}} } = 2q \ left (\ cos 2x- \ cos 2x '\ right) \ chi}

Точнее, если решает уравнение Матье с заданными и, то интеграл ϕ ( Икс ) {\ Displaystyle \ фи (х)} а {\ displaystyle a} q {\ displaystyle q}

ψ ( Икс ) C χ ( Икс , Икс ) ϕ ( Икс ) d Икс {\ Displaystyle \ пси (х) \ эквив \ int _ {C} \ чи (х, х ') \ фи (х') dx '}

где - путь в комплексной плоскости, также решает уравнение Матье с тем же и при соблюдении следующих условий: C {\ displaystyle C} а {\ displaystyle a} q {\ displaystyle q}

  • χ ( Икс , Икс ) {\ Displaystyle \ чи (х, х ')} решает 2 χ Икс 2 - 2 χ Икс 2 знак равно 2 q ( потому что 2 Икс - потому что 2 Икс ) χ {\ displaystyle {\ frac {\ partial ^ {2} \ chi} {\ partial x ^ {2}}} - {\ frac {\ partial ^ {2} \ chi} {\ partial x '^ {2}} } = 2q \ left (\ cos 2x- \ cos 2x '\ right) \ chi}
  • В рассматриваемых регионах существует и является аналитическим ψ ( Икс ) {\ Displaystyle \ psi (х)} χ ( Икс , Икс ) {\ Displaystyle \ чи (х, х ')}
  • ( ϕ χ Икс - ϕ Икс χ ) {\ displaystyle \ left (\ phi {\ frac {\ partial \ chi} {\ partial x '}} - {\ frac {\ partial \ phi} {\ partial x'}} \ chi \ right)} имеет такое же значение в конечных точках C {\ displaystyle C}

Используя соответствующую замену переменных, уравнение для может быть преобразовано в волновое уравнение и решено. Например, одно решение. Примеры полученных таким образом тождеств: χ {\ displaystyle \ chi} χ ( Икс , Икс ) знак равно грех ( 2 q 1 / 2 грех Икс грех Икс ) {\ Displaystyle \ чи (х, х ') = \ зп (2q ^ {1/2} \ грех х \ грех х')}

se 2 п + 1 ( Икс , q ) знак равно se 2 п + 1 ( 0 , q ) π q 1 / 2 B 1 ( 2 п + 1 ) 0 π грех ( 2 q 1 / 2 грех Икс грех Икс ) se 2 п + 1 ( Икс , q ) d Икс ( q gt; 0 ) Ce 2 п ( Икс , q ) знак равно ce 2 п ( π / 2 , q ) π А 0 ( 2 п ) 0 π потому что ( 2 q 1 / 2 шиш Икс потому что Икс ) ce 2 п ( Икс , q ) d Икс       ( q gt; 0 ) {\ displaystyle {\ begin {align} {\ text {se}} _ {2n + 1} (x, q) amp; = {\ frac {{\ text {se}} '_ {2n + 1} (0, q)} {\ pi q ^ {1/2} B_ {1} ^ {(2n + 1)}}} \ int _ {0} ^ {\ pi} \ sinh (2q ^ {1/2} \ sin x \ sin x ') {\ text {se}} _ {2n + 1} (x', q) dx '\ qquad (qgt; 0) \\ {\ text {Ce}} _ {2n} (x, q) amp; = {\ frac {{\ text {ce}} _ {2n} (\ pi / 2, q)} {\ pi A_ {0} ^ {(2n)}}} \ int _ {0} ^ {\ pi} \ cos (2q ^ {1/2} \ cosh x \ cos x ') {\ text {ce}} _ {2n} (x', q) dx '\ qquad \ \ \ (qgt; 0) \ конец {выровнено}}}

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

Также существуют интегральные отношения между функциями первого и второго рода, например:

fe 2 п ( Икс , q ) знак равно 2 п 0 Икс ce 2 п ( τ , - q )   J 0 ( 2 q ( потому что 2 Икс - потому что 2 τ ) ) d τ , п 1 {\ displaystyle {\ text {fe}} _ {2n} (x, q) = 2n \ int _ {0} ^ {x} {\ text {ce}} _ {2n} (\ tau, -q) \ J_ {0} \ left ({\ sqrt {2q (\ cos 2x- \ cos 2 \ tau)}} \ right) d \ tau, \ qquad n \ geq 1}

актуально для любых сложных и реальных. Икс {\ displaystyle x} q {\ displaystyle q}

Асимптотические разложения

Следующие асимптотические разложения справедливы для,, и: q gt; 0 {\ displaystyle qgt; 0} Я ( Икс ) знак равно 0 {\ displaystyle {\ text {Im}} (х) = 0} Re ( Икс ) {\ Displaystyle {\ текст {Re}} (х) \ rightarrow \ infty} 2 q 1 / 2 шиш Икс q 1 / 2 е Икс {\ displaystyle 2q ^ {1/2} \ cosh x \ simeq q ^ {1/2} e ^ {x}}

Ce 2 п ( Икс , q ) ( 2 π q 1 / 2 ) 1 / 2 ce 2 п ( 0 , q ) ce 2 п ( π / 2 , q ) А 0 ( 2 п ) е - Икс / 2 грех ( q 1 / 2 е Икс + π 4 ) Ce 2 п + 1 ( Икс , q ) ( 2 π q 3 / 2 ) 1 / 2 ce 2 п + 1 ( 0 , q ) ce 2 п + 1 ( π / 2 , q ) А 1 ( 2 п + 1 ) е - Икс / 2 потому что ( q 1 / 2 е Икс + π 4 ) Se 2 п + 1 ( Икс , q ) - ( 2 π q 3 / 2 ) 1 / 2 se 2 п + 1 ( 0 , q ) se 2 п + 1 ( π / 2 , q ) B 1 ( 2 п + 1 ) е - Икс / 2 потому что ( q 1 / 2 е Икс + π 4 ) Se 2 п + 2 ( Икс , q ) ( 2 π q 5 / 2 ) 1 / 2 se 2 п + 2 ( 0 , q ) se 2 п + 2 ( π / 2 , q ) B 2 ( 2 п + 2 ) е - Икс / 2 грех ( q 1 / 2 е Икс + π 4 ) {\ displaystyle {\ begin {align} {\ text {Ce}} _ {2n} (x, q) amp; \ sim \ left ({\ frac {2} {\ pi q ^ {1/2}}}} \ справа) ^ {1/2} {\ frac {{\ text {ce}} _ {2n} (0, q) {\ text {ce}} _ {2n} (\ pi / 2, q)} {A_ {0} ^ {(2n)}}} \ cdot e ^ {- x / 2} \ sin \ left (q ^ {1/2} e ^ {x} + {\ frac {\ pi} {4}} \ right) \\ {\ text {Ce}} _ {2n + 1} (x, q) amp; \ sim \ left ({\ frac {2} {\ pi q ^ {3/2}}} \ right) ^ {1/2} {\ frac {{\ text {ce}} _ {2n + 1} (0, q) {\ text {ce}} '_ {2n + 1} (\ pi / 2, q) } {A_ {1} ^ {(2n + 1)}}} \ cdot e ^ {- x / 2} \ cos \ left (q ^ {1/2} e ^ {x} + {\ frac {\ pi } {4}} \ right) \\ {\ text {Se}} _ {2n + 1} (x, q) amp; \ sim - \ left ({\ frac {2} {\ pi q ^ {3/2 }}} \ right) ^ {1/2} {\ frac {{\ text {se}} '_ {2n + 1} (0, q) {\ text {se}} _ {2n + 1} (\ pi / 2, q)} {B_ {1} ^ {(2n + 1)}}} \ cdot e ^ {- x / 2} \ cos \ left (q ^ {1/2} e ^ {x} + {\ frac {\ pi} {4}} \ right) \\ {\ text {Se}} _ {2n + 2} (x, q) amp; \ sim \ left ({\ frac {2} {\ pi q ^ {5/2}}} \ right) ^ {1/2} {\ frac {{\ text {se}} '_ {2n + 2} (0, q) {\ text {se}}' _ { 2n + 2} (\ pi / 2, q)} {B_ {2} ^ {(2n + 2)}}} \ cdot e ^ {- x / 2} \ sin \ left (q ^ {1/2} e ^ {x} + {\ frac {\ pi} {4}} \ right) \ end {align}}}

Таким образом, модифицированные функции Матье экспоненциально убывают при большом действительном аргументе. Подобные асимптотические разложения можно записать для и ; они также экспоненциально затухают при большом действительном аргументе. Fe п {\ displaystyle {\ text {Fe}} _ {n}} Ge п {\ displaystyle {\ text {Ge}} _ {n}}

Для четных и нечетных периодических функций Матье и связанных с ними характеристических чисел можно также получить асимптотические разложения для больших. В частности, для характеристических чисел примерно с нечетным целым числом, т. Е. c е , s е {\ displaystyle ce, se} а {\ displaystyle a} q {\ displaystyle q} N {\ displaystyle N} N N 0 знак равно 2 п + 1 , п знак равно 1 , 2 , 3 , . . . , {\ Displaystyle N \ приблизительно N_ {0} = 2n + 1, n = 1,2,3,...,}

а ( N ) знак равно - 2 q + 2 q 1 / 2 N - 1 2 3 ( N 2 + 1 ) - 1 2 7 q 1 / 2 N ( N 2 + 3 ) - 1 2 12 q ( 5 N 4 + 34 N 2 + 9 ) {\ displaystyle a (N) = - 2q + 2q ^ {1/2} N - {\ frac {1} {2 ^ {3}}} (N ^ {2} +1) - {\ frac {1} {2 ^ {7} q ^ {1/2}}} N (N ^ {2} +3) - {\ frac {1} {2 ^ {12} q}} (5N ^ {4} + 34N ^ {2} +9)}
- 1 2 17 q 3 / 2 N ( 33 N 4 + 410 N 2 + 405 ) - 1 2 20 q 2 ( 63 N 6 + 1260 N 4 + 2943 N 2 + 41807 ) + О ( q - 5 / 2 ) {\ displaystyle \; \; \; \; \; \; \; \; - {\ frac {1} {2 ^ {17} q ^ {3/2}}} N (33N ^ {4} + 410N ^ {2} +405) - {\ frac {1} {2 ^ {20} q ^ {2}}} (63N ^ {6} + 1260N ^ {4} + 2943N ^ {2} +41807) + { \ mathcal {O}} (д ^ {- 5/2})}

Обратите внимание на симметрию при замене и на и, что является важной особенностью расширения. Условия этого расширения были получены явно до срока заказа включительно. Здесь только приблизительно нечетное целое число, потому что в пределе всех минимальных сегментов периодического потенциала становятся фактически независимыми гармоническими осцилляторами (следовательно, нечетным целым). При уменьшении становится возможным туннелирование через барьеры (на физическом языке), что приводит к расщеплению характеристических чисел (в квантовой механике, называемых собственными значениями), соответствующих четным и нечетным периодическим функциям Матье. Это расщепление достигается с помощью граничных условий (в квантовой механике это обеспечивает разбиение собственных значений на энергетические зоны). Граничные условия: q 1 / 2 {\ displaystyle q ^ {1/2}} N {\ displaystyle N} - q 1 / 2 {\ displaystyle -q ^ {1/2}} - N {\ displaystyle -N} | q | - 7 / 2 {\ displaystyle | q | ^ {- 7/2}} N {\ displaystyle N} q {\ displaystyle q \ rightarrow \ infty} потому что 2 Икс {\ displaystyle \ cos 2x} N 0 {\ displaystyle N_ {0}} q {\ displaystyle q} а а {\ displaystyle a \ rightarrow a _ {\ mp}}

( d c е N 0 - 1 d Икс ) π / 2 знак равно 0 , c е N 0 ( π / 2 ) знак равно 0 , ( d s е N 0 d Икс ) π / 2 знак равно 0 , s е N 0 + 1 ( π / 2 ) знак равно 0. {\ displaystyle {\ bigg (} {\ frac {dce_ {N_ {0} -1}} {dx}} {\ bigg)} _ {\ pi / 2} = 0, \; \; ce_ {N_ {0 }} (\ pi / 2) = 0, \; \; {\ bigg (} {\ frac {dse_ {N_ {0}}} {dx}} {\ bigg)} _ {\ pi / 2} = 0, \; \; se_ {N_ {0} +1} (\ pi / 2) = 0.}

Наложив эти граничные условия на асимптотические периодические функции Матье, связанные с указанным выше разложением для, получаем а {\ displaystyle a}

N - N 0 знак равно 2 ( 2 π ) 1 / 2 ( 16 q 1 / 2 ) N 0 / 2 е - 4 q 1 / 2 [ 1 2 ( N 0 - 1 ) ] ! [ 1 - 3 ( N 0 2 + 1 ) 2 6 q 1 / 2 + 1 2 13 q ( 9 N 0 4 - 40 N 0 3 + 18 N 0 2 - 136 N 0 + 9 ) + . . . ] . {\ displaystyle N-N_ {0} = \ mp 2 {\ bigg (} {\ frac {2} {\ pi}} {\ bigg)} ^ {1/2} {\ frac {(16q ^ {1 / 2}) ^ {N_ {0} / 2} e ^ {- 4q ^ {1/2}}} {[{\ frac {1} {2}} (N_ {0} -1)]!}} { \ bigg [} 1 - {\ frac {3 (N_ {0} ^ {2} +1)} {2 ^ {6} q ^ {1/2}}} + {\ frac {1} {2 ^ { 13} q}} (9N_ {0} ^ {4} -40N_ {0} ^ {3} + 18N_ {0} ^ {2} -136N_ {0} +9) +... {\ bigg]}. }

Соответствующие характеристические числа или собственные значения затем следует разложением, т. Е.

а ( N ) знак равно а ( N 0 ) + ( N - N 0 ) ( а N ) N 0 + . . . . {\ Displaystyle а (N) = а (N_ {0}) + (N-N_ {0}) {\ bigg (} {\ frac {\ partial a} {\ partial N}} {\ bigg)} _ { N_ {0}} +....}

Вставка соответствующих выражений выше дает результат

а ( N ) а ( N 0 ) знак равно - 2 q + 2 q 1 / 2 N 0 - 1 2 3 ( N 0 2 + 1 ) - 1 2 7 q 1 / 2 N 0 ( N 0 2 + 3 ) - 1 2 12 q ( 5 N 0 4 + 34 N 0 2 + 9 ) - . . . {\ displaystyle a (N) \ rightarrow a _ {\ mp} (N_ {0}) = - 2q + 2q ^ {1/2} N_ {0} - {\ frac {1} {2 ^ {3}}} (N_ {0} ^ {2} +1) - {\ frac {1} {2 ^ {7} q ^ {1/2}}} N_ {0} (N_ {0} ^ {2} +3) - {\ frac {1} {2 ^ {12} q}} (5N_ {0} ^ {4} + 34N_ {0} ^ {2} +9) -...}
( 16 q 1 / 2 ) N 0 / 2 + 1 е - 4 q 1 / 2 ( 8 π ) 1 / 2 [ 1 2 ( N 0 - 1 ) ] ! [ 1 - N 0 2 6 q 1 / 2 ( 3 N 0 2 + 8 N 0 + 3 ) + . . . ] . {\ Displaystyle \; \; \; \; \; \; \; \; \; \; \; \; \; \; \ mp {\ frac {(16q ^ {1/2}) ^ {N_ { 0} / 2 + 1} e ^ {- 4q ^ {1/2}}} {(8 \ pi) ^ {1/2} [{\ frac {1} {2}} (N_ {0} -1)]!}} {\ bigg [} 1 - {\ frac {N_ {0}} {2 ^ {6} q ^ {1/2}}} (3N_ {0} ^ {2} + 8N_ {0} +3) +... {\ bigg]}.}

Ведь это собственные значения, связанные с четными собственными функциями Матье или (т.е. со знаком минус сверху) и нечетными собственными функциями Матье или (т.е. со знаком плюс). Явные и нормированные разложения собственных функций можно найти в или. N 0 знак равно 1 , 3 , 5 , . . . {\ displaystyle N_ {0} = 1,3,5,...} c е N 0 {\ displaystyle ce_ {N_ {0}}} c е N 0 - 1 {\ displaystyle ce_ {N_ {0} -1}} s е N 0 + 1 {\ displaystyle se_ {N_ {0} +1}} s е N 0 {\ displaystyle se_ {N_ {0}}}

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

Приложения

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

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

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

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

В оставшейся части этого раздела подробно рассматривается анализ двумерного уравнения Гельмгольца. В прямоугольных координатах уравнение Гельмгольца имеет вид

( 2 Икс 2 + 2 у 2 ) ψ + k 2 ψ знак равно 0 , {\ displaystyle \ left ({\ frac {\ partial ^ {2}} {\ partial x ^ {2}}} + {\ frac {\ partial ^ {2}} {\ partial y ^ {2}}} \ вправо) \ psi + k ^ {2} \ psi = 0,}

Эллиптические координаты определяются как

Икс знак равно c шиш μ потому что ν у знак равно c грех μ грех ν {\ Displaystyle {\ begin {align} x amp; = c \ cosh \ mu \ cos \ nu \\ y amp; = c \ sinh \ mu \ sin \ nu \ end {align}}}

где, и положительная постоянная. Уравнение Гельмгольца в этих координатах имеет вид 0 μ lt; {\ Displaystyle 0 \ Leq \ му lt;\ infty} 0 ν lt; 2 π {\ Displaystyle 0 \ Leq \ Nu lt;2 \ pi} c {\ displaystyle c}

1 c 2 ( грех 2 μ + грех 2 ν ) ( 2 μ 2 + 2 ν 2 ) ψ + k 2 ψ знак равно 0 {\ displaystyle {\ frac {1} {c ^ {2} (\ sinh ^ {2} \ mu + \ sin ^ {2} \ nu)}} \ left ({\ frac {\ partial ^ {2}} {\ partial \ mu ^ {2}}} + {\ frac {\ partial ^ {2}} {\ partial \ nu ^ {2}}} \ right) \ psi + k ^ {2} \ psi = 0}

Постоянные кривые - конфокальные эллипсы с фокусным расстоянием ; следовательно, эти координаты удобны для решения уравнения Гельмгольца на областях с эллиптическими границами. Разделение переменных через дает уравнения Матье μ {\ displaystyle \ mu} c {\ displaystyle c} ψ ( μ , ν ) знак равно F ( μ ) грамм ( ν ) {\ Displaystyle \ пси (\ му, \ ню) = F (\ му) G (\ ню)}

d 2 F d μ 2 - ( а - c 2 k 2 2 шиш 2 μ ) F знак равно 0 d 2 грамм d ν 2 + ( а - c 2 k 2 2 потому что 2 ν ) грамм знак равно 0 {\ displaystyle {\ begin {align} amp; {\ frac {d ^ {2} F} {d \ mu ^ {2}}} - \ left (a - {\ frac {c ^ {2} k ^ {2 }} {2}} \ ch 2 \ mu \ right) F = 0 \\ amp; {\ frac {d ^ {2} G} {d \ nu ^ {2}}} + \ left (a - {\ frac {c ^ {2} k ^ {2}} {2}} \ cos 2 \ nu \ right) G = 0 \\\ конец {выровнено}}}

где - постоянная разделения. а {\ displaystyle a}

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

  • Периодичность по, т.е. ν {\ displaystyle \ nu} ψ ( μ , ν ) знак равно ψ ( μ , ν + 2 π ) {\ Displaystyle \ пси (\ му, \ ню) = \ пси (\ му, \ ню +2 \ пи)}
  • Непрерывность смещения по межфокальной линии: ψ ( 0 , ν ) знак равно ψ ( 0 , - ν ) {\ Displaystyle \ psi (0, \ nu) = \ psi (0, - \ nu)}
  • Непрерывность производной по межфокальной линии: ψ μ ( 0 , ν ) знак равно - ψ μ ( 0 , - ν ) {\ Displaystyle \ psi _ {\ mu} (0, \ nu) = - \ psi _ {\ mu} (0, - \ nu)}

Для данного, это ограничивает решения теми, которые имеют форму и, где. Это то же самое, что и ограничение допустимых значений для данного. Ограничения возникают из-за наложения физических условий на некоторую ограничивающую поверхность, такую ​​как эллиптическая граница, определяемая с помощью. Например, зажим мембраны при наложении, что, в свою очередь, требует k {\ displaystyle k} Ce п ( μ , q ) ce п ( ν , q ) {\ displaystyle {\ text {Ce}} _ {n} (\ mu, q) {\ text {ce}} _ {n} (\ nu, q)} Se п ( μ , q ) se п ( ν , q ) {\ displaystyle {\ text {Se}} _ {n} (\ mu, q) {\ text {se}} _ {n} (\ nu, q)} q знак равно c 2 k 2 / 2 {\ Displaystyle д = с ^ {2} к ^ {2} / 2} а {\ displaystyle a} k {\ displaystyle k} k {\ displaystyle k} μ знак равно μ 0 gt; 0 {\ displaystyle \ mu = \ mu _ {0}gt; 0} μ знак равно μ 0 {\ displaystyle \ mu = \ mu _ {0}} ψ ( μ 0 , ν ) знак равно 0 {\ displaystyle \ psi (\ mu _ {0}, \ nu) = 0}

Ce п ( μ 0 , q ) знак равно 0 Se п ( μ 0 , q ) знак равно 0 {\ displaystyle {\ begin {align} {\ text {Ce}} _ {n} (\ mu _ {0}, q) = 0 \\ {\ text {Se}} _ {n} (\ mu _ { 0}, q) = 0 \ end {выровнено}}}

Эти условия определяют нормальные режимы работы системы.

Динамические проблемы

В динамических задачах с периодически меняющимися силами уравнение движения иногда принимает форму уравнения Матье. В таких случаях знание общих свойств уравнения Матье - особенно в отношении устойчивости решений - может быть важным для понимания качественных особенностей физической динамики. Классическим примером в этом отношении является перевернутый маятник. Другие примеры:

Квантовая механика

Функции Матье играют роль в некоторых квантово-механических системах, особенно с пространственно-периодическими потенциалами, такими как квантовый маятник и кристаллические решетки.

Модифицированное уравнение Матье возникает также при описании квантовой механики сингулярных потенциалов. Для частного сингулярного потенциала радиальное уравнение Шредингера V ( р ) знак равно грамм 2 / р 4 {\ Displaystyle V (г) = г ^ {2} / г ^ {4}}

d 2 у d р 2 + [ k 2 - ( + 1 ) р 2 - грамм 2 р 4 ] у знак равно 0 {\ displaystyle {\ frac {d ^ {2} y} {dr ^ {2}}} + \ left [k ^ {2} - {\ frac {\ ell (\ ell +1)} {r ^ {2 }}} - {\ frac {g ^ {2}} {r ^ {4}}} \ right] y = 0}

можно преобразовать в уравнение

d 2 φ d z 2 + [ 2 час 2 шиш 2 z - ( + 1 2 ) 2 ] φ знак равно 0. {\ displaystyle {\ frac {d ^ {2} \ varphi} {dz ^ {2}}} + \ left [2h ^ {2} \ cosh 2z- \ left (\ ell + {\ frac {1} {2 }} \ right) ^ {2} \ right] \ varphi = 0.}

Преобразование достигается следующими заменами

у знак равно р 1 / 2 φ , р знак равно γ е z , γ знак равно я грамм час , час 2 знак равно я k грамм , час знак равно е я π / 4 ( k грамм ) 1 / 2 . {\ displaystyle y = r ^ {1/2} \ varphi, r = \ gamma e ^ {z}, \ gamma = {\ frac {ig} {h}}, h ^ {2} = ikg, h = e ^ {I \ pi / 4} (кг) ^ {1/2}.}

Решая уравнение Шредингера (для этого конкретного потенциала) в терминах решений модифицированного уравнения Матье, можно получить свойства рассеяния, такие как S-матрица и поглощательная способность.

Смотрите также
Заметки
Рекомендации
Внешние ссылки
Последняя правка сделана 2024-01-01 11:41:23
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте