插值方程求根拟合(2009建模课).docVIP

插 值 问题：已知y是以x为自变量的一元函数，当x分别取下列数据时 ， 对应的y的值分别是 . 求任意一点x处的函数值y . 名词：称数据为插值节点，简称节点；称数据为节点处的函数值。 插值法：根据，构造函数f(x)使得 . Lagrange插值 定义1 称为Lagrange插值基函数，简称基函数。 定义2 称 为Lagrange插值多项式。 注：L(x)就可以作为原插值问题的答案，这种解法称为Lagrange插值法。 用Matlab做Lagrange插值计算： 将定义1、定义2编写成Matlab程序，程序清单如下： function y=lagr1(x0,y0,x) n=length(x0);m=length(x); for i=1:m z=x(i); s=0.0; for k=1:n p=1.0; for j=1:n if j~=k p=p*(z-x0(j))/(x0(k)-x0(j)); end end s=s+y0(k)*p; end y(i)=s; end 该程序文件名记为“lagr1.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.8, ……, 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=lagr1(x0,y0,x) (注释：输出所要求的那76个函数值) plot(x,y) (注释：以x为横坐标、y为纵坐标画出图象) 执行结果： 2. 分段线性插值 3. 三次样条插值 定义3 若有一个分段函数S(x)，满足下面三条： (1). 在每个小区间 上是一个三次多项式； (2). 在每个分段点 处二阶连续可导； (3). . 则称S(x)是插值问题的三次样条插值函数。 用Matlab做三次样条插值计算： 输入数据： 算函数值：y=spline(x0,y0,x) 画出图象：plot(x,y) 如：上面例2中，用三次样条插值法确定函数y=y(x)，该函数图象为 对函数， 令节点为 -5、-3、-1、1、3、5，令， 分别用Lagrange法、分段线性法、三次样条法作插值计算，并画图象 解 程序清单为： hold on x0=[-5,-3,-1,1,3,5];y0=1./(1+x0.^2); x=-5:0.1:5; y1=1./(1+x.^2); y2=lagr1(x0,y0,x); y3=spline(x0,y0,x); plot(x,y1),plot(x,y2),plot(x0,y0),plot(x,y3) hold off 运行结果： 方程求跟 有少数方程求根很容易，如：方程有两个根，是1、-2；方程有一个根，是1. 大多数方程求根不容易，如 就求不出准确根（即：一点误差都没有的根），只能用数值解法求近似根. f(x)连续，f(a)与f(b)异号，ab,则方程f(x)=0在区间(a,b)内至少有一个根，称(a,b)是该方程的一个有根区间. 若已知(a,b)内有且仅有一个根，则称(a,b)是一个单根区间. 确定了单根区间(a,b)后，常用的数值求根的方法是：二分法；牛顿法. (1).二分法 例 用二分法求在(1,2)内的一个根，误差0.05. 解 设，则 有根区间是(1,2) ，有根区间(1,1.5) ， 有根区间(1,1.25) , 有根区间(1.125,1.25) 有根区间(1.125,1.1875) 取为近似根，可使误差小于 . (2). 牛顿法 方程f(x)=0. 估计一个近似根(误差较小)，记为，在曲线上的点处作切线，切线方程为，该切线与OX轴交点记为，则， . 在曲线上的点处作切线，该切线与OX轴交点记为，则， . 一般地有 称为牛顿迭代公式，这样得到的数列的极限就是准确根. 例 ，用牛顿法求根。 解 ，根据牛顿迭代公式得 得到方程的一个近似根1.1640，误差小于0.0001. 注意：若f(x)的导数不易求，则，估计两个近似根(误差较小)，分别记为，用替代牛顿迭代公式中的，得到新的迭代公式 . 通常求方程的根