原创力文档- 1、本文档共4页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
龙格现象matlab算法
实验报告
课程名称：___计算方法____________指导老师：___程晓良________成绩：__________________
实验名称：___观察龙格现象________________实验类型：________________同组学生姓名：__________
一、实验目的和要求（必填） 二、实验内容和原理（必填）
三、主要仪器设备（必填） 四、操作方法和实验步骤
五、实验数据记录和处理 六、实验结果与分析（必填）
七、讨论、心得
一问题描述
例外发生了：龙格在研究多项式插值的时候，发现有的情况下，并非取节点（日期数）越多多项式就越精确。著名的例子是f（x）=1/(1+25x^2).它的插值函数在两个端点处发生剧烈的波动，造成较大的误差。
二相关公式
MATLAB程序
一、取等距节点，n=5,10,15,20
for n = 5:5:20
subplot(2,2,n/5)
syms x;
f = 1/(1+25*x^2);
x1=sym(zeros(n+1));
W=sym(ones(n+1));
L=sym(0);
for i=0:n
x1(i+1)=-1+2*i/n;
end
for i=0:n
for j=0:n
if j~=i
w=(x-x1(j+1))/(x1(i+1)-x1(j+1));
W(i+1)=W(i+1)*w;
end
end
L=L+W(i+1)*(1/(1+25*x1(i+1)^2));
end
LL(n)=simplify(L);
x=-1:0.01:1;
y1=subs(f,x);
y2=subs(L,x);
plot(x,y1,b);hold on;
plot(x,y2,r);hold off;
title([?-oˉêyf(x)=1/(1+25*x^2)ó?,num2str(n),′?2??μoˉêy]);
xlabel(x);ylabel(y);
legend(?-oˉêy,2??μoˉêy);
grid on
end
二、取节点Xj=cosjπ/n,j=0,1,…,n.n分别取…,50
for n = 5:5:50
subplot(2,5,n/5)
syms x;
f = 1/(1+25*x^2);
x1=sym(zeros(n+1));
W=sym(ones(n+1));
L=sym(0);
for i=0:n
x1(i+1)=cos(i*pi/n);
end
for i=0:n
for j=0:n
if j~=i
w=(x-x1(j+1))/(x1(i+1)-x1(j+1));
W(i+1)=W(i+1)*w;
end
end
L=L+W(i+1)*(1/(1+25*x1(i+1)^2));
end
LL(n)=simplify(L);
x=-1:0.01:1;
y1=subs(f,x);
y2=subs(L,x);
plot(x,y1,b);hold on;
plot(x,y2,r);hold off;
title([?-oˉêyf(x)=1/(1+25*x^2)ó?,num2str(n),′?2??μoˉêy]);
xlabel(x);ylabel(y);
legend(?-oˉêy,2??μoˉêy);
grid on
end
四、实验分析
x=cos为节点时随着n的增大
五心得与体会
3130102509____
日期：________________
地点：_______桌号
装 订 线
文档评论（0）