陈磊10203多项式插值的震荡现象Word格式文档下载.docx
- 文档编号:20831195
- 上传时间:2023-01-25
- 格式:DOCX
- 页数:12
- 大小:80.28KB
陈磊10203多项式插值的震荡现象Word格式文档下载.docx
《陈磊10203多项式插值的震荡现象Word格式文档下载.docx》由会员分享,可在线阅读,更多相关《陈磊10203多项式插值的震荡现象Word格式文档下载.docx(12页珍藏版)》请在冰豆网上搜索。
,对
进行多项式插值,观察比较
的图形与插值多项式
图形.
(2)取
进行三次样条插值(取自然边界条件),观察比较
的图形与样条插值函数
(3)对函数
重复上述第
(1)项.
三、实验结果的分析
图1:
当n=2时函数
的拉格朗日多项式插值图
图2:
当n=5时函数
图3:
当n=10时函数
图4:
当n=20时函数
结果分析:
多项式插值逼近结果如图所示,图1为插值结点为2,图2为插值结点为5,图3为插值结点为10,图4为插值结点为20。
当输入插值多项式的次数由2个结点增加为5个结点时,函数图形越接近,当插值结点增加到10个结点时,曲线在中间段比较光滑,但是在图3两端时会出现与原函数
偏差会很大。
当插值次数增加到20个结点时,可以看出,适当提高插值多项式次数,可以提高逼近的精度,同时也出现偏差的地方偏差也会很大。
当输入插值多项式的次数由2个结点时原函数图与插值图的偏差还是很大的,随着插值结点由5个到10个图像变得越光滑,只会在两端出现偏差较大现象。
当插值次数增加到20时,可以看出,适当提高插值多项式次数,可以提高逼近的精度,但同时也出现偏差的地方偏差也会增大。
第二问程序图:
当n=10时,观察比较
从上图可以看出图像拟合程度比较好,当n=10时,图像上在两端结点分布较多,在中间结点数较少;
当n=5时,图像可以看出结点分布较均匀,每段曲线间长度大致相等。
4、实验程序
第一问和第三问程序:
functiont_charpt2
result=inputdlg({'
选择试验题目,若选2.1,2.3,请输入1:
'
},'
charpt_2'
1,{'
1'
});
Nb=str2num(char(result));
if(Nb~=1)&
&
(Nb~=2),errordlg('
选择试验题目错误!
);
return;
end
promps={'
请选择试验函数,若选f(x),请输入f,若选H(x),请输入H'
};
result=inputdlg(promps,'
charpt2'
f'
Nb_f=char(result);
if(Nb_f~='
Nb_f~='
H'
)
errordlg('
试验函数选择错误!
请输入插值多项式的次数N:
10'
Nd=str2num(char(result));
if(Nd<
1)
插值多项式的次数输入错误!
switchNb_f
case'
f=inline('
1./(1+x.^2)'
a=-5;
b=5;
x./(1+x.^4)'
if(Nb==1)
x0=linspace(a,b,Nd+1);
y0=feval(f,x0);
x=a:
0.1:
b;
y=Lagrange(x0,y0,x);
clf;
fplot(f,[ab],'
co'
holdon;
plot(x,y,'
b--'
xlabel('
x'
ylabel('
y=f(x)oandy=Ln(x)--'
%--------------------------------------------------
functiony=Lagrange(x0,y0,x)
n=length(x0);
m=length(x);
fori=1:
m
z=x(i);
s=0.0;
fork=1:
n
p=1.0;
forj=1:
if(j~=k)
p=p*(z-x0(j))/(x0(k)-x0(j));
end
s=s+p*y0(k);
y(i)=s;
第二问程序:
当n=10时;
x=-5:
0.01:
5;
y=1./(1+x.^2);
x1=-5:
1:
y1=interp1(x,y,x1,'
spline'
plot(x1,y1,'
o'
x,y)
holdon
y'
y1;
当n=5时;
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 陈磊 10203 多项式 震荡 现象