原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
数值分析实验指导书数值分析实验指导书
《数值分析》
实 验 指 导 书
潍坊学院数学与信息科学学院
2012年04月
目 录
目录 I
实验一 插值与曲线拟合的最小二乘法 1
实验二 数值积分 4
实验三 解线性方程组的直接法 9
实验四 解线性方程组的迭代法 11
实验五 非线性方程的数值解法 13
实验六 常微分方程数值解法 17
实验一 插值与曲线拟合的最小二乘法
一、实验目的：
1．了解拉格朗日插值法、牛顿插值法、曲线拟合最小二乘法的基本原理和方法；
2．掌握拉格朗日插值多项式牛顿插值多项式的用法；
3．掌握最小二乘原理，会求拟合函数及超定方程组的最小二乘解。
二、实验内容：
1．用拉格朗日插值公式和牛顿插值公式确定函数值；
2．对函数f (x)进行拉格朗日插值和牛顿插值；
3．利用Polyfit拟合幂函数，利用Polyfit拟合多项式。
三、实验过程：
1．给定函数四个点的数据如下:
，
试用插值公式确定函数在处的函数值。
MATLAB程序如下：
X=[1.1,2.3,3.9,5.1]; Y =[3.877,4.726,4.651 ,2.117];
p1=poly(X(1)); p2=poly(X(2));
p3=poly(X(3)); p4=poly(X(4));
l01= conv ( conv (p2, p3), p4)/(( X(1)- X(2))* ( X(1)- X(3)) * ( X(1)- X(4))),
l11= conv ( conv (p1, p3), p4)/(( X(2)- X(1))* ( X(2)- X(3)) * ( X(2)- X(4))),
l21= conv ( conv (p1, p2), p4)/(( X(3)- X(1))* ( X(3)- X(2)) * ( X(3)- X(4))),
l31= conv ( conv (p1, p2), p3)/(( X(4)- X(1))* ( X(4)- X(2)) * ( X(4)- X(3))),
l0=poly2sym (l01),
l1=poly2sym (l11),l2=poly2sym (l21), l3=poly2sym (l31),
P = l01* Y(1)+ l11* Y(2) + l21* Y(3) + l31* Y(4),
运行后输出的基函数l0，l1，l2和l3为
l0 =-1/24*x^3+1/8*x^2-1/12*x，l1 =1/4*x^3-1/4*x^2-x+1
l2 =-1/3*x^3+4/3*x，l3 =1/8*x^3+1/8*x^2-1/4*x
输入程序
L=poly2sym (P),x=2.101; Y = polyval(P,x)
运行后输出插值多项式和插值为
L=-629/5376*x^3+31433/53760*x^2-63029765850741/281474976710656*x+2010616283501353/562949953421312
Y =4.5969
输入程序
L=poly2sym (P),x=4.234; Y = polyval(P,x)
运行后输出插值多项式和插值为
L=-629/5376*x^3+31433/53760*x^2-63029765850741/281474976710656*x+2010616283501353/562949953421312
Y =4.2244
L=145616387951645/9007199254740992*x^3-2517512191700115/4503599627370496*x^2+14477/6000*x+2007/1000
Y = 3.4290
输入程序
syms M; x=2.101;
R3=M*abs((x-X(1))*(x-X(2)) *(x-X(3)) *(x-X(4)))/24
运行后输出误差限为
R3 =435065974692861/36028797018963968*M
2．在区间上取结点数，等距间隔的节点为插值点，对于函数进行拉格朗日插值。
MATLAB程序如下
t=-5: 1:5;
ft=(1+t.*t).\5;
t1=-5:1:5;
ft1=(1+t1.*t1).\5;
y1=Lagran(t1,ft1,t);
plot(t,ft,b:,t,y1,g+);
xlabel(x);ylabel(y);
对一组数据做拉格朗日的M文件如下
Lagran．m
function fi =Lagran(x,f,xi)
fi=zeros(size(xi));
npl=length(f);
for i=1:npl
z=ones(siz
文档评论（0）