How do I specify m for orthpoly::legendre ?

Question

Billy 2011-6-14

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/9485-how-do-i-specify-m-for-orthpoly-legendre

采纳的回答： Billy

I want to get the associated Legendre polynomials in mupad, but only the Legendre (not assoc.) are avail.

http://www.mathworks.com/help/toolbox/mupad/orthpoly/legendre.html

Is there a way to get this or can I call the MATLAB legendre func from Mupad?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Billy 2011-6-14

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/9485-how-do-i-specify-m-for-orthpoly-legendre#answer_13074

在 MATLAB Online 中打开

Using http://mathworld.wolfram.com/AssociatedLegendrePolynomial.html, you can express the associated legendre polynomials in terms of the unassoc ones:

use(orthpoly,legendre)
P := (l,m) -> piecewise( 
[m>0, (-1)^m * (1-x^2)^(m/2) * diff( legendre( l, x ), x $ m ) ], 
[m=0, legendre( l, x )],
[m<0, (-1)^(-m) * (l+m)!/(l-m)! * P( l,-m )] )

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Answer 2

Billy 2012-5-17

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/9485-how-do-i-specify-m-for-orthpoly-legendre#answer_48151

在 MATLAB Online 中打开

In fact a better way to do this:

Define Legendre polynomials:

P := (n) -> 1/(2^n * n!) * diff( (x^2-1)^n, x $ n )
// and the associated legendre polynomials:
Plm := (l,m) -> ( 1 - x^2 )^(m/2) * diff( P(l), x $ m )
// The normalizing factor:
Klm := (l,m) -> sqrt( ( (2*l+1)/(4*PI) ) * ( ( (l-m)!) / ((l+m)!) ) )
// Finally, the spherical harmonics
Ylm := (l,m,t,p) -> piecewise(
[m>0,   sqrt(2)*Klm(l,m) * (Plm(l,m)|x=cos(t)) * cos(m*p) ],
[m=0,   Klm(l,m) * (Plm(l,m)|x=cos(t)) ],
[m<0,   sqrt(2)*Klm(l,-m) * (Plm(l,-m)|x=cos(t)) *sin(-m*p) ]
)
// Generate a plot:
plot( plot::Spherical( [ abs( Ylm(4,0,t,p) ), p, t], p=0..2*PI, t=0..PI,
UMesh=60, VMesh=60, FillColorType=Flat, Color=[.8,.5,.2] ), Scaling=Constrained )
// truncate values greater than 0
lt0 := (v) -> (1-heaviside(v)) * v 
// truncate values less than 0
gt0 := (v) -> heaviside(v)*v 
// plot
plot(
plot::Spherical( [ abs(gt0(Ylm(5,4,t,p))), p, t], p=0..2*PI, t=0..PI,
UMesh=60, VMesh=60, FillColorType=Flat, Color=[0,.6,0] ),
plot::Spherical( [ abs(lt0(Ylm(5,4,t,p))), p, t], p=0..2*PI, t=0..PI,
UMesh=60, VMesh=60, FillColorType=Flat, Color=[.4,0,.0] ),
Scaling=Constrained
)

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Walter Roberson 2012-5-17

(Note: the above is MuPAD code and must be run from within a MuPad notepad, or coded carefully into an evalin(symengine) call, or written to a file and the file called upon within MuPAD)

请先登录，再进行评论。