Associated legendre polynomials fail after certain degree

7 次查看（过去 30 天）
ailbeildce
ailbeildce 2017-12-16
编辑： Elvis Alexander Agüero Vera 2023-1-16
Hi,
I am using legendre polynomials for an application on spherical harmonics. However the code
legendre(170,0.5)
where 170 is the degree/order fail, giving me Inf or NaN. Is this considered a bug or is there way to aid the issue using higher precision somehow?
Best

请先登录，再回答此问题。

采纳的回答

David Goodmanson
David Goodmanson 2017-12-19
Hi ailbeildce,
Try legendre(n,x,'norm') or legendre(n,x,'sch'). Each of these normalizes the associated legendre function slightly differently, and both leave out a factor in front that gets out of hand in a big way as m gets large [where m is the upper parameter in Pmn, 0<=m<=n, and m=0 corresponds to the usual Pn].
With either of those options, n can go up to at least 2400.
You can see what the factors are in 'doc legendre'. You will have to check, but I think the 'norm' option for Pmn gives you
Int{-1,1} Pmn(x)^2 dx = 1,
appropriate for spherical harmonics.
  3 个评论
显示 1更早的评论
David Goodmanson
David Goodmanson 2018-3-30
You're very welcome. I should probably know, but what is a PM feature?
Elvis Alexander Agüero Vera
Elvis Alexander Agüero Vera 2023-1-16
编辑：Elvis Alexander Agüero Vera 2023-1-16
I guess he refers to a private mesage.
Somewhat related question: I also need to calculate with efficiency the derivatives of the legendre Polynomials. I would appreciate a fast way of computing that.
Also, why is it that
f = matlabFunction(diff(legendreP(50, x), x))
is so unstable for degrees greater than, say, 50?

请先登录，再进行评论。

更多回答（1 个）

Walter Roberson
Walter Roberson 2017-12-16
If you have the symbolic toolbox you can work with it
https://www.mathworks.com/help/symbolic/legendrep.html
  2 个评论
ailbeildce
ailbeildce 2017-12-17
编辑：ailbeildce 2017-12-17
legendre() gives out more information than legendreP. Although I don't know if there's a way to generate Y_l^m where m!=0 with legendreP.
Walter Roberson
Walter Roberson 2017-12-17
For integer m you can see https://en.wikipedia.org/wiki/Associated_Legendre_polynomials#Definition_for_non-negative_integer_parameters_%E2%84%93_and_m which the formula given in terms of derivatives. As the different orders correspond to different numbers of derivatives of the Legendre polynomial, you can find the different orders in a loop.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Polynomials 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by