Что делает калькулятор таблицы присоединённых полиномов Лежандра?
Этот инструмент рассчитывает таблицу значений присоединённой функции Лежандра \(P_n^m(x)\) (степень \(n\), порядок \(m\)) на выбранном интервале \(x\) и строит соответствующую кривую. Это чистая математика: результат одинаков в любой стране, без единиц измерения и каких-либо национальных особенностей. Присоединённые полиномы Лежандра встречаются повсюду в физике и прикладной математике — в сферических гармониках, при решении уравнения Лапласа в сферических координатах, в мультипольных разложениях и в квантовой механике углового момента.
Как пользоваться
Введите целую степень \(n\) (0, 1, 2, ...) и целый порядок \(m\), для которого \(-n \le m \le n\). Задайте начальное значение \(x\) (от -1 до 1), шаг и количество строк. Значения по умолчанию \(n = 2\), \(m = 1\), начало = -1, шаг = 0,02, 101 строка проходят \(x\) от -1 до +1 включительно. Выберите тип A (соглашение Wolfram) или тип B (соглашение Maple): для вещественных \(x\) на интервале (-1, 1) они совпадают по модулю и различаются лишь знаком или фазой множителя.
Разбор формулы
Для целого \(n\) и \(0 \le m \le n\) применяется выражение $$P_n^m(x) = (-1)^m\,(1-x^2)^{m/2}\,\frac{d^m}{dx^m}P_n(x)$$ которое вычисляется через численно устойчивую рекуррентную формулу: \(P_m^m = (-1)^m(2m-1)!!(1-x^2)^{m/2}\), \(P_{m+1}^m = x(2m+1)P_m^m\), далее \((l-m)P_l^m = (2l-1)x P_{l-1}^m - (l+m-1)P_{l-2}^m\). Для отрицательного \(m\): \(P_n^{-m} = (-1)^m\left[\frac{(n-m)!}{(n+m)!}\right]P_n^m\). Такая замкнутая рекуррентность позволяет избежать переполнения из-за гамма-функции, характерного для прямой записи через 2F1 при положительном целом \(m\).
Разбор примера
При \(n = 2\), \(m = 1\) функция имеет вид $$P_2^1(x) = -3x\sqrt{1-x^2}$$ При \(x = 0\) значение равно 0; при \(x = 0{,}5\) оно составляет \(-3(0{,}5)(0{,}866025) = -1{,}299038\); при \(x = -0{,}5\) получаем \(+1{,}299038\). Кривая стартует с 0 (\(x = -1\)), поднимается примерно до \(+1{,}1547\) вблизи \(x = -0{,}577\), пересекает ноль при \(x = 0\), опускается примерно до \(-1{,}1547\) вблизи \(x = +0{,}577\) и возвращается к 0 при \(x = +1\).
Частые вопросы
Почему \(n\) и \(m\) должны быть целыми? Полином обрывается (становится конечным) только при неотрицательном целом \(n\); рекуррентная формула и множители \((n\pm m)!\) требуют целого \(m\) с условием \(-n \le m \le n\).
Что за значение показано в карточке результата? В верхнем блоке выводятся \(x\) и \(P_n^m(x)\) для средней строки таблицы (медианный индекс) — быстрая контрольная точка для проверки кривой.
Что получается при \(m = 0\)? \(P_n^0(x)\) — это обычный полином Лежандра \(P_n(x)\).
Замкнутые формы присоединённых функций Лежандра P_n^m(x)
Присоединённые функции Лежандра \(P_n^m(x)\) для целого степени \(n\) и порядка \(0\le m\le n\) следуют из \(P_n^m(x)=(-1)^m(1-x^2)^{m/2}\dfrac{d^m}{dx^m}P_n(x)\). Множитель \((-1)^m\) — это фаза Кондона–Шортли, включённая в соглашение Тип A (совпадает с Wolfram); соглашение Тип B (Maple) его опускает, поэтому его элементы с нечётным \(m\) отличаются только знаком. Таблица ниже приводит явные формы в соответствии с Типом A.
| \(n\) | \(m\) | \(P_n^m(x)\) (Тип A, со знаком) |
|---|---|---|
| 0 | 0 | \(1\) |
| 1 | 0 | \(x\) |
| 1 | 1 | \(-\sqrt{1-x^2}\) |
| 2 | 0 | \(\tfrac{1}{2}(3x^2-1)\) |
| 2 | 1 | \(-3x\sqrt{1-x^2}\) |
| 2 | 2 | \(3(1-x^2)\) |
| 3 | 0 | \(\tfrac{1}{2}(5x^3-3x)\) |
| 3 | 1 | \(-\tfrac{3}{2}(5x^2-1)\sqrt{1-x^2}\) |
| 3 | 2 | \(15x(1-x^2)\) |
| 3 | 3 | \(-15(1-x^2)^{3/2}\) |
В качестве проверки, при \(x=0.5\) элемент \(P_2^1\) даёт \(-3(0.5)\sqrt{1-0.25}=-1.5\sqrt{0.75}=\) -1.299038. Столбец \(m=0\) воспроизводит обычные полиномы Лежандра \(P_n(x)\), например \(P_3^0(x)=\tfrac12(5x^3-3x)\), которые можно табулировать с помощью калькулятора таблицы полиномов Лежандра.
Ключевые термины и переменные
- Степень \(n\)
-
Неотрицательное целое число (
degreeN), определяющее порядок подлежащего полинома Лежандра \(P_n(x)\), который является полиномом степени \(n\). - Порядок \(m\)
-
Целое число (
orderM), управляющее числом берущихся производных. Для действительных результатов на \((-1,1)\) обычно используют \(0\le m\le n\); когда \(m>n\) функция тождественно равна нулю, потому что \(m\)-я производная полинома степени \(n\) обращается в нуль. - Аргумент \(x\)
-
Точка вычисления (
initialXплюс \(i\cdot\)stepX). Функции действительны при \(-1\le x\le 1\); в физике \(x=\cos\theta\). - Тип A (Wolfram / Кондон–Шортли)
-
Включает множитель фазы \((-1)^m\). Это соглашение, используемое
LegendrePWolfram и стандартными учебниками квантовой механики. - Тип B (Maple)
- Опускает фазовый множитель \((-1)^m\), поэтому \(P_n^m\) (Тип B) \(=(-1)^m\,P_n^m\) (Тип A). Величины идентичны; только знак элементов с нечётным \(m\) отличается.
- Двойной факториал \((2m-1)!!\)
- Произведение нечётных целых чисел \((2m-1)(2m-3)\cdots 3\cdot 1\), где \((-1)!!=1\). Он появляется в старшем коэффициенте \(P_m^m(x)=(-1)^m(2m-1)!!\,(1-x^2)^{m/2}\); например, \(P_3^3\) использует \(5!!=15\). Смотрите калькулятор двойного факториала для этих значений.
- Соотношение для отрицательного порядка
- При \(m>0\), \(P_n^{-m}(x)=(-1)^m\dfrac{(n-m)!}{(n+m)!}\,P_n^{m}(x)\), связывающее положительные и отрицательные порядки через факториалы.
Интерпретация таблицы и графика
Несколько структурных свойств позволяют вам проверить табулированные значения и построенную кривую:
- Чётность. \(P_n^m(-x)=(-1)^{n+m}P_n^m(x)\). Когда \(n+m\) чётно, график симметричен относительно \(x=0\); когда \(n+m\) нечётно, он антисимметричен (и, следовательно, проходит через начало координат).
- Нули во внутренней части. На открытом интервале \((-1,1)\), \(P_n^m(x)\) имеет ровно \(n-m\) простых нулей. Например \(P_3^1\) имеет два нуля во внутренней части, в то время как \(P_n^n\) не имеет ни одного.
- Поведение на концах. Из-за множителя \((1-x^2)^{m/2}\) каждая функция с \(m>0\) обращается в нуль при \(x=\pm 1\). При \(m=0\) значения равны \(P_n(1)=1\) и \(P_n(-1)=(-1)^n\).
- Величина вблизи краёв. При большем \(m\) множитель \((1-x^2)^{m/2}\) резко подавляет кривую при \(x\to\pm1\), поэтому наибольшие отклонения происходят в направлении средней части диапазона.
Эти функции — \(\theta\)-зависимая часть сферических гармоник \(Y_n^m(\theta,\phi)\): записав \(x=\cos\theta\), получаем \(Y_n^m\propto P_n^m(\cos\theta)\,e^{im\phi}\). Внутренние нули становятся нодальными кругами широты, а обращение в нуль при \(m>0\) на концах соответствует тому, что гармоники стремятся к нулю на полюсах. Те же самые значения \(P_n^m\) таким образом непосредственно поступают в вычисление сферических гармоник при выбранных \(\theta\) и \(\phi\).
