Polynomial sequence
График полиномиальной функции Якоби P n ( α , β ) {\displaystyle P_{n}^{(\alpha ,\beta )}} с n = 10 {\displaystyle n=10} и α = 2 {\displaystyle \alpha =2} и β = 2 {\displaystyle \beta =2} в комплексной плоскости из − 2 − 2 i {\displaystyle -2-2i} к 2 + 2 i {\displaystyle 2+2i} с цветами, созданными с помощью функции Mathematica 13.1 ComplexPlot3D В математике полиномы Якоби (иногда называемые гипергеометрическими полиномами ) P n ( α , β ) ( x ) {\displaystyle P_{n}^{(\alpha ,\beta )}(x)} являются классом классических ортогональных многочленов . Они ортогональны по весу. ( 1 − x ) α ( 1 + x ) β {\displaystyle (1-x)^{\alpha }(1+x)^{\beta }} на интервале [ − 1 , 1 ] {\displaystyle [-1,1]} . Полиномы Гегенбауэра , а, следовательно, и полиномы Лежандра , Цернике и Чебышева , являются частными случаями полиномов Якоби. [1]
Полиномы Якоби были введены Карлом Густавом Якобом Якоби .
Полиномы Якоби определяются через гипергеометрическую функцию следующим образом: [2]
P n ( α , β ) ( z ) = ( α + 1 ) n n ! 2 F 1 ( − n , 1 + α + β + n ; α + 1 ; 1 2 ( 1 − z ) ) , {\displaystyle P_{n}^{(\alpha ,\beta )}(z)={\frac {(\alpha +1)_{n}}{n!}}\,{}_{2}F_{1}\left(-n,1+\alpha +\beta +n;\alpha +1;{\tfrac {1}{2}}(1-z)\right),} где ( α + 1 ) n {\displaystyle (\alpha +1)_{n}} является символом Поххаммера (падающего факториала). В этом случае ряд для гипергеометрической функции конечен, поэтому получается следующее эквивалентное выражение:
P n ( α , β ) ( z ) = Γ ( α + n + 1 ) n ! Γ ( α + β + n + 1 ) ∑ m = 0 n ( n m ) Γ ( α + β + n + m + 1 ) Γ ( α + m + 1 ) ( z − 1 2 ) m . {\displaystyle P_{n}^{(\alpha ,\beta )}(z)={\frac {\Gamma (\alpha +n+1)}{n!\,\Gamma (\alpha +\beta +n+1)}}\sum _{m=0}^{n}{n \choose m}{\frac {\Gamma (\alpha +\beta +n+m+1)}{\Gamma (\alpha +m+1)}}\left({\frac {z-1}{2}}\right)^{m}.} Эквивалентное определение даёт формула Родригеса : [1] [3]
P n ( α , β ) ( z ) = ( − 1 ) n 2 n n ! ( 1 − z ) − α ( 1 + z ) − β d n d z n { ( 1 − z ) α ( 1 + z ) β ( 1 − z 2 ) n } . {\displaystyle P_{n}^{(\alpha ,\beta )}(z)={\frac {(-1)^{n}}{2^{n}n!}}(1-z)^{-\alpha }(1+z)^{-\beta }{\frac {d^{n}}{dz^{n}}}\left\{(1-z)^{\alpha }(1+z)^{\beta }\left(1-z^{2}\right)^{n}\right\}.} Если α = β = 0 {\displaystyle \alpha =\beta =0} , то оно сводится к полиномам Лежандра :
P n ( z ) = 1 2 n n ! d n d z n ( z 2 − 1 ) n . {\displaystyle P_{n}(z)={\frac {1}{2^{n}n!}}{\frac {d^{n}}{dz^{n}}}(z^{2}-1)^{n}\;.} Альтернативное выражение для реального аргумента [ редактировать ] Серьезно x {\displaystyle x} полином Якоби альтернативно можно записать как
P n ( α , β ) ( x ) = ∑ s = 0 n ( n + α n − s ) ( n + β s ) ( x − 1 2 ) s ( x + 1 2 ) n − s {\displaystyle P_{n}^{(\alpha ,\beta )}(x)=\sum _{s=0}^{n}{n+\alpha \choose n-s}{n+\beta \choose s}\left({\frac {x-1}{2}}\right)^{s}\left({\frac {x+1}{2}}\right)^{n-s}} и для целого числа n {\displaystyle n}
( z n ) = { Γ ( z + 1 ) Γ ( n + 1 ) Γ ( z − n + 1 ) n ≥ 0 0 n < 0 {\displaystyle {z \choose n}={\begin{cases}{\frac {\Gamma (z+1)}{\Gamma (n+1)\Gamma (z-n+1)}}&n\geq 0\\0&n<0\end{cases}}} где Γ ( z ) {\displaystyle \Gamma (z)} это гамма-функция .
В частном случае, когда четыре величины n {\displaystyle n} , n + α {\displaystyle n+\alpha } , n + β {\displaystyle n+\beta } , n + α + β {\displaystyle n+\alpha +\beta } являются неотрицательными целыми числами, полином Якоби можно записать как
P n ( α , β ) ( x ) = ( n + α ) ! ( n + β ) ! ∑ s = 0 n 1 s ! ( n + α − s ) ! ( β + s ) ! ( n − s ) ! ( x − 1 2 ) n − s ( x + 1 2 ) s . {\displaystyle P_{n}^{(\alpha ,\beta )}(x)=(n+\alpha )!(n+\beta )!\sum _{s=0}^{n}{\frac {1}{s!(n+\alpha -s)!(\beta +s)!(n-s)!}}\left({\frac {x-1}{2}}\right)^{n-s}\left({\frac {x+1}{2}}\right)^{s}.} ( 1 )
Сумма распространяется на все целые значения s {\displaystyle s} для которого аргументы факториалов неотрицательны.
P 0 ( α , β ) ( z ) = 1 , {\displaystyle P_{0}^{(\alpha ,\beta )}(z)=1,} P 1 ( α , β ) ( z ) = ( α + 1 ) + ( α + β + 2 ) z − 1 2 , {\displaystyle P_{1}^{(\alpha ,\beta )}(z)=(\alpha +1)+(\alpha +\beta +2){\frac {z-1}{2}},} P 2 ( α , β ) ( z ) = ( α + 1 ) ( α + 2 ) 2 + ( α + 2 ) ( α + β + 3 ) z − 1 2 + ( α + β + 3 ) ( α + β + 4 ) 2 ( z − 1 2 ) 2 . {\displaystyle P_{2}^{(\alpha ,\beta )}(z)={\frac {(\alpha +1)(\alpha +2)}{2}}+(\alpha +2)(\alpha +\beta +3){\frac {z-1}{2}}+{\frac {(\alpha +\beta +3)(\alpha +\beta +4)}{2}}\left({\frac {z-1}{2}}\right)^{2}.} Полиномы Якоби удовлетворяют условию ортогональности
∫ − 1 1 ( 1 − x ) α ( 1 + x ) β P m ( α , β ) ( x ) P n ( α , β ) ( x ) d x = 2 α + β + 1 2 n + α + β + 1 Γ ( n + α + 1 ) Γ ( n + β + 1 ) Γ ( n + α + β + 1 ) n ! δ n m , α , β > − 1. {\displaystyle \int _{-1}^{1}(1-x)^{\alpha }(1+x)^{\beta }P_{m}^{(\alpha ,\beta )}(x)P_{n}^{(\alpha ,\beta )}(x)\,dx={\frac {2^{\alpha +\beta +1}}{2n+\alpha +\beta +1}}{\frac {\Gamma (n+\alpha +1)\Gamma (n+\beta +1)}{\Gamma (n+\alpha +\beta +1)n!}}\delta _{nm},\qquad \alpha ,\ \beta >-1.} Согласно определению, они не имеют единичной нормы по весу. Это можно исправить, разделив на квадратный корень правой части приведенного выше уравнения, когда n = m {\displaystyle n=m} .
Хотя это не дает ортонормированного базиса, иногда предпочтительнее альтернативная нормализация из-за ее простоты:
P n ( α , β ) ( 1 ) = ( n + α n ) . {\displaystyle P_{n}^{(\alpha ,\beta )}(1)={n+\alpha \choose n}.} Полиномы имеют соотношение симметрии
P n ( α , β ) ( − z ) = ( − 1 ) n P n ( β , α ) ( z ) ; {\displaystyle P_{n}^{(\alpha ,\beta )}(-z)=(-1)^{n}P_{n}^{(\beta ,\alpha )}(z);} таким образом, другое терминальное значение равно
P n ( α , β ) ( − 1 ) = ( − 1 ) n ( n + β n ) . {\displaystyle P_{n}^{(\alpha ,\beta )}(-1)=(-1)^{n}{n+\beta \choose n}.} The k {\displaystyle k} производная явного выражения приводит к
d k d z k P n ( α , β ) ( z ) = Γ ( α + β + n + 1 + k ) 2 k Γ ( α + β + n + 1 ) P n − k ( α + k , β + k ) ( z ) . {\displaystyle {\frac {d^{k}}{dz^{k}}}P_{n}^{(\alpha ,\beta )}(z)={\frac {\Gamma (\alpha +\beta +n+1+k)}{2^{k}\Gamma (\alpha +\beta +n+1)}}P_{n-k}^{(\alpha +k,\beta +k)}(z).} Полином Якоби P n ( α , β ) {\displaystyle P_{n}^{(\alpha ,\beta )}} является решением линейного однородного дифференциального уравнения второго порядка [1]
( 1 − x 2 ) y ″ + ( β − α − ( α + β + 2 ) x ) y ′ + n ( n + α + β + 1 ) y = 0. {\displaystyle \left(1-x^{2}\right)y''+(\beta -\alpha -(\alpha +\beta +2)x)y'+n(n+\alpha +\beta +1)y=0.} Рекуррентное соотношение для полиномов Якоби фиксированной α {\displaystyle \alpha } , β {\displaystyle \beta } является: [1]
2 n ( n + α + β ) ( 2 n + α + β − 2 ) P n ( α , β ) ( z ) = ( 2 n + α + β − 1 ) { ( 2 n + α + β ) ( 2 n + α + β − 2 ) z + α 2 − β 2 } P n − 1 ( α , β ) ( z ) − 2 ( n + α − 1 ) ( n + β − 1 ) ( 2 n + α + β ) P n − 2 ( α , β ) ( z ) , {\displaystyle {\begin{aligned}&2n(n+\alpha +\beta )(2n+\alpha +\beta -2)P_{n}^{(\alpha ,\beta )}(z)\\&\qquad =(2n+\alpha +\beta -1){\Big \{}(2n+\alpha +\beta )(2n+\alpha +\beta -2)z+\alpha ^{2}-\beta ^{2}{\Big \}}P_{n-1}^{(\alpha ,\beta )}(z)-2(n+\alpha -1)(n+\beta -1)(2n+\alpha +\beta )P_{n-2}^{(\alpha ,\beta )}(z),\end{aligned}}} для n = 2 , 3 , … {\displaystyle n=2,3,\ldots } .Пишу для краткости a := n + α {\displaystyle a:=n+\alpha } , b := n + β {\displaystyle b:=n+\beta } и c := a + b = 2 n + α + β {\displaystyle c:=a+b=2n+\alpha +\beta } , это становится с точки зрения a , b , c {\displaystyle a,b,c}
2 n ( c − n ) ( c − 2 ) P n ( α , β ) ( z ) = ( c − 1 ) { c ( c − 2 ) z + ( a − b ) ( c − 2 n ) } P n − 1 ( α , β ) ( z ) − 2 ( a − 1 ) ( b − 1 ) c P n − 2 ( α , β ) ( z ) . {\displaystyle 2n(c-n)(c-2)P_{n}^{(\alpha ,\beta )}(z)=(c-1){\Big \{}c(c-2)z+(a-b)(c-2n){\Big \}}P_{n-1}^{(\alpha ,\beta )}(z)-2(a-1)(b-1)c\;P_{n-2}^{(\alpha ,\beta )}(z).} Поскольку полиномы Якоби можно описать с помощью гипергеометрической функции, рекурренты гипергеометрической функции дают эквивалентные рекурренты полиномов Якоби. В частности, смежные отношения Гаусса соответствуют тождествам
( z − 1 ) d d z P n ( α , β ) ( z ) = 1 2 ( z − 1 ) ( 1 + α + β + n ) P n − 1 ( α + 1 , β + 1 ) = n P n ( α , β ) − ( α + n ) P n − 1 ( α , β + 1 ) = ( 1 + α + β + n ) ( P n ( α , β + 1 ) − P n ( α , β ) ) = ( α + n ) P n ( α − 1 , β + 1 ) − α P n ( α , β ) = 2 ( n + 1 ) P n + 1 ( α , β − 1 ) − ( z ( 1 + α + β + n ) + α + 1 + n − β ) P n ( α , β ) 1 + z = ( 2 β + n + n z ) P n ( α , β ) − 2 ( β + n ) P n ( α , β − 1 ) 1 + z = 1 − z 1 + z ( β P n ( α , β ) − ( β + n ) P n ( α + 1 , β − 1 ) ) . {\displaystyle {\begin{aligned}(z-1){\frac {d}{dz}}P_{n}^{(\alpha ,\beta )}(z)&={\frac {1}{2}}(z-1)(1+\alpha +\beta +n)P_{n-1}^{(\alpha +1,\beta +1)}\\&=nP_{n}^{(\alpha ,\beta )}-(\alpha +n)P_{n-1}^{(\alpha ,\beta +1)}\\&=(1+\alpha +\beta +n)\left(P_{n}^{(\alpha ,\beta +1)}-P_{n}^{(\alpha ,\beta )}\right)\\&=(\alpha +n)P_{n}^{(\alpha -1,\beta +1)}-\alpha P_{n}^{(\alpha ,\beta )}\\&={\frac {2(n+1)P_{n+1}^{(\alpha ,\beta -1)}-\left(z(1+\alpha +\beta +n)+\alpha +1+n-\beta \right)P_{n}^{(\alpha ,\beta )}}{1+z}}\\&={\frac {(2\beta +n+nz)P_{n}^{(\alpha ,\beta )}-2(\beta +n)P_{n}^{(\alpha ,\beta -1)}}{1+z}}\\&={\frac {1-z}{1+z}}\left(\beta P_{n}^{(\alpha ,\beta )}-(\beta +n)P_{n}^{(\alpha +1,\beta -1)}\right)\,.\end{aligned}}} полиномов Производящая функция Якоби определяется выражением
∑ n = 0 ∞ P n ( α , β ) ( z ) t n = 2 α + β R − 1 ( 1 − t + R ) − α ( 1 + t + R ) − β , {\displaystyle \sum _{n=0}^{\infty }P_{n}^{(\alpha ,\beta )}(z)t^{n}=2^{\alpha +\beta }R^{-1}(1-t+R)^{-\alpha }(1+t+R)^{-\beta },} где
R = R ( z , t ) = ( 1 − 2 z t + t 2 ) 1 2 , {\displaystyle R=R(z,t)=\left(1-2zt+t^{2}\right)^{\frac {1}{2}}~,} и ветвь квадратного корня выбирается так, что R ( z , 0 ) = 1 {\displaystyle R(z,0)=1} . [1]
Для x {\displaystyle x} в интерьере [ − 1 , 1 ] {\displaystyle [-1,1]} , асимптотика P n ( α , β ) {\displaystyle P_{n}^{(\alpha ,\beta )}} для больших n {\displaystyle n} определяется формулой Дарбу [1]
P n ( α , β ) ( cos θ ) = n − 1 2 k ( θ ) cos ( N θ + γ ) + O ( n − 3 2 ) , {\displaystyle P_{n}^{(\alpha ,\beta )}(\cos \theta )=n^{-{\frac {1}{2}}}k(\theta )\cos(N\theta +\gamma )+O\left(n^{-{\frac {3}{2}}}\right),} где
k ( θ ) = π − 1 2 sin − α − 1 2 θ 2 cos − β − 1 2 θ 2 , N = n + 1 2 ( α + β + 1 ) , γ = − π 2 ( α + 1 2 ) , 0 < θ < π {\displaystyle {\begin{aligned}k(\theta )&=\pi ^{-{\frac {1}{2}}}\sin ^{-\alpha -{\frac {1}{2}}}{\tfrac {\theta }{2}}\cos ^{-\beta -{\frac {1}{2}}}{\tfrac {\theta }{2}},\\N&=n+{\tfrac {1}{2}}(\alpha +\beta +1),\\\gamma &=-{\tfrac {\pi }{2}}\left(\alpha +{\tfrac {1}{2}}\right),\\0<\theta &<\pi \end{aligned}}} и " O {\displaystyle O} " член однороден на интервале [ ε , π − ε ] {\displaystyle [\varepsilon ,\pi -\varepsilon ]} для каждого ε > 0 {\displaystyle \varepsilon >0} .
Асимптотика полиномов Якоби вблизи точек ± 1 {\displaystyle \pm 1} задается формулой Мелера – Гейне
lim n → ∞ n − α P n ( α , β ) ( cos ( z n ) ) = ( z 2 ) − α J α ( z ) lim n → ∞ n − β P n ( α , β ) ( cos ( π − z n ) ) = ( z 2 ) − β J β ( z ) {\displaystyle {\begin{aligned}\lim _{n\to \infty }n^{-\alpha }P_{n}^{(\alpha ,\beta )}\left(\cos \left({\tfrac {z}{n}}\right)\right)&=\left({\tfrac {z}{2}}\right)^{-\alpha }J_{\alpha }(z)\\\lim _{n\to \infty }n^{-\beta }P_{n}^{(\alpha ,\beta )}\left(\cos \left(\pi -{\tfrac {z}{n}}\right)\right)&=\left({\tfrac {z}{2}}\right)^{-\beta }J_{\beta }(z)\end{aligned}}} где пределы едины для z {\displaystyle z} в ограниченной области .
Асимптотика снаружи [ − 1 , 1 ] {\displaystyle [-1,1]} является менее явным.
Выражение ( 1 ) позволяет выразить d-матрицу Вигнера d m ′ , m j ( ϕ ) {\displaystyle d_{m',m}^{j}(\phi )} (для 0 ≤ ϕ ≤ 4 π {\displaystyle 0\leq \phi \leq 4\pi } )в терминах полиномов Якоби: [4]
d m ′ m j ( ϕ ) = ( − 1 ) m − m ′ − | m − m ′ | 2 [ ( j + M ) ! ( j − M ) ! ( j + N ) ! ( j − N ) ! ] 1 2 ( sin ϕ 2 ) | m − m ′ | ( cos ϕ 2 ) | m + m ′ | P j − m ( | m − m ′ | , | m + m ′ | ) ( cos ϕ ) , {\displaystyle d_{m'm}^{j}(\phi )=(-1)^{\frac {m-m'-|m-m'|}{2}}\left[{\frac {(j+M)!(j-M)!}{(j+N)!(j-N)!}}\right]^{\frac {1}{2}}\left(\sin {\tfrac {\phi }{2}}\right)^{|m-m'|}\left(\cos {\tfrac {\phi }{2}}\right)^{|m+m'|}P_{j-m}^{(|m-m'|,|m+m'|)}(\cos \phi ),}
где M = max ( | m | , | m ′ | ) , N = min ( | m | , | m ′ | ) {\displaystyle M=\max(|m|,|m'|),N=\min(|m|,|m'|)} .
^ Перейти обратно: а б с д и ж Сегё, Габор (1939). «IV. Полиномы Якоби». Ортогональные полиномы . Публикации коллоквиума. Том. XXIII. Американское математическое общество. ISBN 978-0-8218-1023-1 . МР 0372517 . Определение находится в IV.1; дифференциальное уравнение – в IV.2; Формула Родригеса находится в IV.3; производящая функция находится в IV.4; рекуррентное отношение находится в IV.5; асимптотическое поведение находится в VIII.2 ^ Абрамовиц, Милтон ; Стегун, Ирен Энн , ред. (1983) [июнь 1964 г.]. «Глава 22» . Справочник по математическим функциям с формулами, графиками и математическими таблицами . Серия «Прикладная математика». Том. 55 (Девятое переиздание с дополнительными исправлениями десятого оригинального издания с исправлениями (декабрь 1972 г.); первое изд.). Вашингтон, округ Колумбия; Нью-Йорк: Министерство торговли США, Национальное бюро стандартов; Дуврские публикации. п. 561. ИСБН 978-0-486-61272-0 . LCCN 64-60036 . МР 0167642 . LCCN 65-12253 . ^ П. К. Суетин (2001) [1994], «Полиномы Якоби» , Энциклопедия математики , EMS Press ^ Биденхарн, LC; Лук, Джей Ди (1981). Угловой момент в квантовой физике . Чтение: Аддисон-Уэсли. Эндрюс, Джордж Э.; Аски, Ричард; Рой, Ранджан (1999), Специальные функции , Энциклопедия математики и ее приложений, том. 71, Издательство Кембриджского университета , ISBN 978-0-521-62321-6 , МР 1688958 , ISBN 978-0-521-78988-2 Коорнвиндер, Том Х.; Вонг, Родерик СК; Кукук, Рулоф; Свартау, Рене Ф. (2010), «Ортогональные полиномы» , в Олвере, Фрэнке В.Дж .; Лозье, Дэниел М.; Буасверт, Рональд Ф.; Кларк, Чарльз В. (ред.), Справочник NIST по математическим функциям , издательство Кембриджского университета, ISBN 978-0-521-19225-5 , МР 2723248 .
скрывать Международный Национальный Другой