插值方程求根拟合.docxVIP

插值 问题：已知y是以x为自变量的一元函数，当x 分别取下列数据时 xQxix2--xn , 对应的y的值分别是儿，X，『2,…，儿? 求任意一点x处的函数值y ? 名词：称数据兀为插值节点，简称节点；称数据X 为节点兀?处的函数值。 插值法：根据（无,必），i = 0,l,2,…,兀,构造函数f（x） 使得/(兀)=必，i = 0,l,2,…,n . 1. Lagrange 插值 定义1 lj(x)=-兀0）…（勺一勺.1）（?_?+l) lj(x)= -兀0）…（勺一勺.1）（? _?+l)…(勺_兀) （其中，丿分别取0,1, 2, ?..,/!） 称为Lagrange插值基函数，简称基函数。 定义 2 称 ￡(x) = yjlj(x)为 Lagrange 插值 7=0 多项式。 注：L(x)就可以作为原插值问题的答案，这种解 法称为Lagrange插值法。 用Matlab做Lagrange插值计算： 将定义1、定义2编写成Matlab程序，程序清单 如下： function y=lagrl(xO,yO,x) n=leng th(x0);m=le ngth(x); for i=l:m z=x(i); s=0.0; for k=l:n P=1.0； for j=l:n ifj ?二k p=p*(z-xO(j))/(xO(k)-xO(j)); end s=s+y0(k)*p; end y(i)=s； end 该程序文件名记为“lagrl.m”，存放在Matlab工作区。 例 X 0 3 5 7 9 11 12 13 14 15 y 0 1.2 1.7 2.0 2.1 2.0 1.8 1.2 1.0 1.6  求当X分别取0, 0.2, 0.4, 0.6, 0.……,15时，对应的 函数值y,并画出函数y=y(x)的图象. 用Matlab作计算： x0 二[0,3,5,7,9,11,12,13,14,15]; y0=[0,1.2,1.7,2,2.1,2,1.8,1.2,1,1.6]; x=0:0.2:15; y=lagrl(xO,yO,x) (注释：输111所要求的那76个函数值) plot(x,y) (注释：以x为横坐标、y为纵坐标画出图象) 执行结果： 2.分段线性插值 2.分段线性插值 3.三次样条插值 定义3若有一个分段函数S(x),满足下面三条: .在每个小区间[g?+J (心0,1,2,???也-1) 上是一个三次多项式； .在每个分段点兀(=1,2,???/-1)处二阶连 续可导； . S(xJ = % (z = 0,1,2,???,/!). 则称S(x)是插值问题广 X°州…耳的三次样 U：儿X…儿 络桶值函数- 用Matlab做三次样条插值计算： 输入数据:x0 = [x0,xJ;yO = [yo,--]; 输入数据: x = x0:（小步长）：￡； 算函数值：y=spline(xO,yO,x) 画出图象：plot(x,y) 如：上面例2中，用三次样条插值法确定函数 y=y(x),该函数图象为 例3 对函数y =—厶，xg [-5,5], ? 1 + JC 令节点兀为?5、?3、?1、1、3、5,令y.=—^—f 1 1 1 + # 分别用Lagrange法、分段线性法、三次样条法作插值 计算，并画图象 解程序清单为： hold on x0=[-5,-3,-l ,1,3,5] ;y0= l./(l +xO.A2); x=-5:0」:5; yl=l./(l+x.A2); y2=lagrl(x0,y0,x); y3=spline(x0,y0,x); plot(x,y 1 )，plot(x,y2),plot(x0,y0),plot(x,y3) hold off 运行结果： 方程求跟有少数方程求根很容易，女山方程”+ — 2 = 0有 两个根,是1、-2；方程ln(x) = 0有一个根，是］. 大多数方程求根不容易，如 x4+x-3 = 0, xln(x) = 1, x = e~x 就求不出准确根(即：一点误差都没有的根)，只能用 数值解法求近似根. f(x)连续，班小与班1))异号，ab,则方程f (x) =0 在区间(a, b)内至少有一个根，称(a, b)是该方程的一 个有根区间. 若已知(a, b)内有且仅有一个根，则称(a, b)是一 个单根区间. 确定了单根区间(a, b)后，常用的数值求根的方法 是：二分法；牛顿法. .-分法 例用二分法求〃+兀一3二0在(1,2)内的一个根， 误差50. 05. 解 设 /(x) = x4+x-3，贝!J /(1)0,/(2)0,有根区间是(1,2) /(1)0,/(2)0, 有根区间是(1,2) /(1.5) = 1.54+1.5-30, 有根区间(1,1.