原创力文档- 1、本文档共8页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
6.1.2 特殊函数
TP321MA/124 MATLAB 教程
% notebook -setup;
% ctrl-enter
% alt-d
由特殊函数构成的完备正交基函数序列可分为两类：
在区间 [-1，1] 的两个端点满足自然边界条件且权函数为 1 的正交基，有
计算连带勒让德函数 的语句为
P = legendre(n,x);
y = P(m+1,:)
在区间 [0，1] 的左端点满足自然边界条件,右端点满足齐次第一类边界条件且权函数为 x 的正交基,有
计算贝塞尔函数Jn的语句为
Besselj(n,x);
第二章中的S-L定理已经在理论上保证了上述函数系的正交性。借助数值计算方法来验证基函数的正交性，则是本小节的工作。
编程算例 1 试验证勒让德函数 的正交性。
解
x = linspace(-1,1,100);
P = legendre(3,x);
plot(x,P(1,:));
x = linspace(-1,1,100);
P = legendre(4,x);
plot(x,P(1,:));
%验证正交性
format long;figure;
x = linspace(-1,1,100);
P = legendre(3,x);
plot(x,P(1,:));
y1 = P(1,:);
hold on;
P = legendre(4,x);
plot(x,P(1,:));
y2 = P(1,:);
y1*y2/100
ans =
1.443289932012704e-017
编程算例 2 勒让德函数 的绘图。
x = linspace(-1,1,100);
P = legendre(1,x);
plot(x,P(1,:));
hold on;
x = linspace(-1,1,100);
P = legendre(2,x);
plot(x,P(1,:));
hold on;
x = linspace(-1,1,100);
P = legendre(3,x);
plot(x,P(1,:));
hold on;
x = linspace(-1,1,100);
P = legendre(4,x);
plot(x,P(1,:));
x = linspace(-1,1,1000);
P = legendre(1,x);
y1=P(1,:);
P = legendre(2,x);
y2=P(1,:);
P = legendre(3,x);
y3=P(1,:);
P = legendre(4,x);
y4=P(1,:);
figure;
plot(x,y1,r-,x,y2,g--,x,y3,b:,x,y4,k-.),grid;
title(勒让德函数的图形);
xlabel(x);
ylabel(legendre);
legend(P1,P2,P3,P4);
编程算例 3 连带勒让德函数 的绘图。
x = linspace(-1,1,1000);
P = legendre(1,x);
y1=P(2,:);
P = legendre(2,x);
y2=P(2,:);
P = legendre(3,x);
y3=P(2,:);
P = legendre(4,x);
y4=P(2,:);
figure;
plot(x,-y1,r-,x,-y2,g--,x,-y3,b:,x,-y4,k-.),grid;
%title(连带勒让德函数的图形);
xlabel(x);
ylabel(legendre);
legend(P(1)1,P(1)2,P(1)3,P(1)4);
文档评论（0）