特殊函数的计算机仿真应用23.pdfVIP

第二十三章 特殊函数的计算机仿真应用 23.1 连带勒让德函数 勒让德函数 球函数 23.1.1 连带勒让德函数 MATLAB计算连带勒让德函数用下面语句： p=legendre(n, x) 图 23.1 三阶连带勒让德函数图形分布 例23.1.1 画出所有三阶连带勒让德函数的图形。 仿真程序如下，如图23.1所示。 x=0:0.01:1; y=legendre(3,x); plot(x,y(1,:),-,x,y(2,:),-.,x,y(3,: ),:,x,y(4,:),--) legend(P_3^0,P_3^1,P_3^2,P_3^3) 图 23.1 三阶连带勒让德函数图形分布 23.1.2 勒让德多项式 例23.1.2 画出勒让德多项式的图形 x=0:0.01:1; y1=legendre(1,x); 图23.2 勒让德多项式曲线 y2=legendre(2,x); y3=legendre(3,x); y4=legendre(4,x); y5=legendre(5,x); plot(x,y1(1,:),x,y2(1,:),x,y3(1,:),x,y4(1,:),x,y5(1,:)) title(Legendre) 图23.2 勒让德多项式曲线 23.1.3.球谐函数 根据公式 m (l−| m |)!(2l +1) |m| imϕ θϕ θ Y ( , ) P (cos )e l l (l+| m |)!4π 可以利用连带勒让德函数绘出球谐函数的图形分布. 23.1.4. 勒让德多项式的母函数图形 根据勒让德多项式的母函数公式，有 ⎧ ∞ n ∑r P (cosθ) (r 1) ⎪ n 1 ⎪ n 0 ⎨ 2 ∞ 1−2 cos + r θ r − −1 ⎪∑r n P (cosθ) (r 1) ⎪ n n 0 ⎩ 其中满足r 1 ，且在第一象限的母函数（即上式左端） 图形分布如图23.3 所示 其仿真程序为 [X,Z]=meshgrid([0:0.1:3],[0:0.1:2]); [Q,R]=cart2pol(X,Z); R(find(R 1))=NaN; u=1./sqrt(1-2.*R.*cos(Q)+R.^2); meshc(X,Z,u) Rin=R; Rin(find(Rin1))=NaN; Rout=R; Rout(find(Rout1))=NaN; Uin=1; 图23.3 勒让德母函数在第一象限图形 Uout=1./Rout; for k=1:20