插值法在测井数据处理中的应用_基于Excel的数据处理

插值法在测井数据处理中的应用_基于Excel的数据处理
插值法在测井数据处理中的应用_基于Excel的数据处理

牛顿插值法原理及应用

牛顿插值法 插值法是利用函数f (x)在某区间中若干点的函数值,作出适当的特定函数,在这些点上取已知值,在区间的其他点上用这特定函数的值作为函数f (x)的近似值。如果这特定函数是多项式,就称它为插值多项式。当插值节点增减时全部插值基函数均要随之变化,这在实际计算中很不方便。为了克服这一缺点,提出了牛顿插值。牛顿插值通过求各阶差商,递推得到的一个公式: f(x)=f[x0]+f[x0,x1](x-x0)+f[x0,x1,x2](x-x0)(x-x1)+...f[x0,...xn](x-x0 )...(x-xn-1)+Rn(x)。 插值函数 插值函数的概念及相关性质[1] 定义:设连续函数y-f(x) 在区间[a,b]上有定义,已知在n+1个互异的点 x0,x1,…xn上取值分别为y0,y1,…yn (设a≤ x1≤x2……≤xn≤b)。若在函数类中存在以简单函数P(x) ,使得P(xi)=yi,则称P(x) 为f(x)的插值函数. 称x1,x2,…xn 为插值节点,称[a,b]为插值区间。 定理:n次代数插值问题的解存在且唯一。

牛顿插值法C程序 程序框图#include void main() { float x[11],y[11][11],xx,temp,newton; int i,j,n; printf("Newton插值:\n请输入要运算的值:x="); scanf("%f",&xx); printf("请输入插值的次数(n<11):n="); scanf("%d",&n); printf("请输入%d组值:\n",n+1); for(i=0;i

三次样条插值方法的应用

CENTRAL SOUTH UNIVERSITY 数值分析实验报告

三次样条插值方法的应用 一、问题背景 分段低次插值函数往往具有很好的收敛性,计算过程简单,稳定性好,并且易于在在电子计算机上实现,但其光滑性较差,对于像高速飞机的机翼形线船体放样等型值线往往要求具有二阶光滑度,即有二阶连续导数,早期工程师制图时,把富有弹性的细长木条(即所谓的样条)用压铁固定在样点上,在其他地方让他自由弯曲,然后沿木条画下曲线,称为样条曲线。样条曲线实际上是由分段三次曲线并接而成,在连接点即样点上要求二阶导数连续,从数学上加以概括就得到数学样条这一概念。下面我们讨论最常用的三次样条函数及其应用。 二、数学模型 样条函数可以给出光滑的插值曲线(面),因此在数值逼近、常微分方程和偏微分方程的数值解及科学和工程的计算中起着重要的作用。 设区间[]b ,a 上给定有关划分b x x n =<<<= 10x a ,S 为[]b ,a 上满足下面条件的函数。 ● )(b a C S ,2∈; ● S 在每个子区间[]1,+i i x x 上是三次多项式。 则称S 为关于划分的三次样条函数。常用的三次样条函数的边界条件有三种类型: ● Ⅰ型 ()()n n n f x S f x S ''0'',==。 ● Ⅱ型 ()()n n n f x S f x S ''''0'''',==,其特殊情况为()()0''''==n n x S x S 。 ● Ⅲ型 ()() 3,2,1,0,0==j x S x S n j j ,此条件称为周期样条函数。 鉴于Ⅱ型三次样条插值函数在实际应用中的重要地位,在此主要对它进行详细介绍。 三、算法及流程 按照传统的编程方法,可将公式直接转换为MATLAB 可是别的语言即可;另一种是运用矩阵运算,发挥MATLAB 在矩阵运算上的优势。两种方法都可以方便地得到结果。方法二更直观,但计算系数时要特别注意。这里计算的是方法一的程序,采用的是Ⅱ型边界条件,取名为spline2.m 。 Matlab 代码如下: function s=spline2(x0,y0,y21,y2n,x) %s=spline2(x0,y0,y21,y2n,x) %x0,y0 are existed points,x are insert points,y21,y2n are the second

几种插值法比较与应用

多种插值法比较与应用 (一)Lagrange 插值 1. Lagrange 插值基函数 n+1个n 次多项式 ∏ ≠=--=n k j j j k j k x x x x x l 0)( n k ,,1,0ΛΛ= 称为Lagrange 插值基函数 2. Lagrange 插值多项式 设给定n+1个互异点))(,(k k x f x ,n k ,,1,0ΛΛ=,j i x x ≠,j i ≠,满足插值条件 )()(k k n x f x L =,n k ,,1,0ΛΛ= 的n 次多项式 ∏∏ ∏=≠==--==n k n k j j j k j k k n k k n x x x x x f x l x f x L 0 00 ))(()()()( 为Lagrange 插值多项式,称 ∏=+-+=-=n j j x n n x x n f x L x f x E 0 )1()()!1()()()()(ξ 为插值余项,其中),()(b a x x ∈=ξξ (二)Newton 插值 1.差商的定义 )(x f 关于i x 的零阶差商 )(][i i x f x f = )(x f 关于i x ,j x 的一阶差商

i j i j j i x x x f x f x x f --= ][][],[ 依次类推,)(x f 关于i x ,1+i x ,……,k i x +的k 阶差商 i k i k i i k i i k i i i x x x x f x x f x x x f --= +-+++++] ,,[],,[],,,[111ΛΛΛΛΛ 2. Newton 插值多项式 设给定的n+1个互异点))(,(k k x f x ,n k ,,1,0ΛΛ=,j i x x ≠,j i ≠, 称满足条件 )()(k k n x f x N =,n k ,,1,0ΛΛ= 的n 次多项式 )()](,,,[)](,[][)(10100100---++-+=n n n x x x x x x x f x x x x f x f x N ΛΛΛΛΛ 为Newton 插值多项式,称 ],[,)(],,,[)()()(010b a x x x x x x f x N x f x E n j j n n ∈-=-=∏=ΛΛ 为插值余项。 (三)Hermite 插值 设],[)(1b a C x f ∈,已知互异点0x ,1x ,…,],[b a x n ∈及所对应的函数值为0f ,1f ,…,n f ,导数值为'0f ,'1f ,…,'n f ,则满足条件 n i f x H f x H i i n i i n ,,1,0,)(,)(''1212Λ===++ 的12+n 次Hermite 插值多项式为 )()()(0 '12x f x f x H j n j j j n j i n βα∏∏=++= 其中 )())((,)]()(21[)(2 2'x l x x x l x l x x x j j j j j j j j ---=βα

几种插值法的应用和比较

插值法的应用与比较 信科1302 万贤浩 13271038 1格朗日插值法 在数值分析中,拉格朗日插值法是以法国十八世纪数学家约瑟夫·路易斯·拉格朗日命名的一种多项式插值方法.许多实际问题中都用函数来表示某种内在联系或规律,而不少函数都只能通过实验和观测来了解.如对实践中的某个物理量进行观测,在若干个不同的地方得到相应的观测值,拉格朗日插值法可以找到一个多项式,其恰好在各个观测的点取到观测到的值.这样的多项式称为拉格朗日(插值)多项式.数学上来说,拉格朗日插值法可以给出一个恰好穿过二维平面上若干个已知点的多项式函数.拉格朗日插值法最早被英国数学家爱德华·华林于1779年发现,不久后由莱昂哈德·欧拉再次发现.1795年,拉格朗日在其著作《师范学校数学基础教程》中发表了这个插值方法,从此他的名字就和这个方法联系在一起. 1.1拉格朗日插值多项式 图1 已知平面上四个点:(?9, 5), (?4, 2), (?1, ?2), (7, 9),拉格朗日多项式:)(x L (黑色)穿过所有点.而每个基本多项式:)(00x l y ,)(11x l y , )(22x l y 以及)(x l y ??各穿过对应的一点,并在其它的三个点的x 值上取零. 对于给定的若1+n 个点),(00y x ,),(11y x ,………),(n n y x ,对应于它们的次数不超过n 的拉格朗日多项式L 只有一个.如果计入次数更高的多项式,则有无穷个,因为所有与L 相差 ))((10x x x x --λ……)(n x x -的多项式都满足条件. 对某个多项式函数,已知有给定的1+k 个取值点: ),(00y x ,……,),(k k y x ,

几种常用的插值方法

数学系 信息与计算科学1班 李平 指导老师:唐振先 摘要:插值在诸如机械加工等工程技术和数据处理等科学研究中有许多直接的应用,在很多领域都要用插值的办法找出表格和中间值,插值还是数值积分微分方程数值解等数值计算的基础。本文归纳了几种常用的插值方法,并简单分析了其各自的优缺点。 关键词:任意阶多项式插值,分段多项式插值。 引言:所谓插值,通俗地说就是在若干以知的函数值之间插入一些未知函数值,而插值函数的类型最简单的选取是代数多项式。用多项式建立插值函数的方法主要用两种:一种是任意阶的插值多项式,它主要有三种基本的插值公式:单项式,拉格朗日和牛顿插值;另一种是分段多项式插值,它有Hermite 和spine 插值和分段线性插值。 一.任意阶多项式插值: 1.用单项式基本插值公式进行多项式插值: 多项式插值是求通过几个已知数据点的那个n-1阶多项式,即P n-1(X)=A 1+A 2X+…A n X n-1,它是一个单项式基本函数X 0,X 1…X n-1的集合来定义多项式,由已知n 个点(X,Y )构成的集合,可以使多项式通过没数据点,并为n 个未知系数Ai 写出n 个方程,这n 个方程组成的方程组的系数矩阵为Vandermonde 矩阵。 虽然这个过程直观易懂,但它都不是建立插值多项式最好的办法,因为Vandermonde 方程组有可能是病态的,这样会导致单项式系数不确定。另外,单项式中的各项可能在大小上有很大的差异,这就导致了多项式计算中的舍入误差。 2.拉格朗日基本插值公式进行插值: 先构造一组插值函数L i (x ) =011011()()()() ()()()() i i n i i i i i i n x x x x x x x x x x x x x x x x -+-+--------,其中i=0,… n.容易看出n 次多项式L i (x )满足L i (x )=1,(i=j );L i (x )=0,(i ≠j ),其中i=0,1…n ,令L i (x )=0()n i i i y l x =∑这就是拉格朗日插值多项式。与单项式基本 函数插值多项式相比,拉格朗日插值有2个重要优点:首先,建立插值多项式不需要求解方程组;其次,它的估计值受舍入误差要小得多。拉格朗日插值公式结构

插值法的分类与应用

插值法的方法与应用 武汉科技大学城市建设学院 琚婷婷 结构工程 201108710014 【摘要】文章讨论插值法在数值分析中的中心地位和重要作用,比较插值法间的优缺点,应用以及各种方法之间的相互联系。 【关键词】插值法;应用。 1.插值问题的提出 在许多实际问题及科学研究中,因素之间往往存在着函数关系,但是这些关系的显示表达式不一定都知道,通常只是由观察或测试得到一些离散数值,所以只能从这些数据构造函数的近似表达式,有时虽然给出了解析表达式,但由于解析表达式过于复杂,使用或计算起来十分麻烦。这就需要建立函数的某种近似表达,而插值法就是构造函数的近似表达式的方法。 2.插值法的数学表达 由于代数多项式是最简单而又便于计算的函数,所以经常采用多项式作为插值函数,称为多项式插值。多项式插值法有拉格朗日插值法,牛顿插值法、埃尔米特插值法,分段插值法和样条插值法等。其基本思想都是用高次代数多项式或分段的低次多项式作为被插值函数f (x)的近似解析表达式。 3.常用多项式插值公式构造 (I)拉格朗日插值 n 次拉格朗日插值多项式p n (x)对可表示为 p n (x)= y i l i (x)n i=0= y i ( x ?x j x i ?x j n j ≠0i ≠j n i=0) 其中l i x ,i =0,1,2???,n 称为插值基函数,插值余项为: R n (x)= f (x)- p n (x)=f n +1 (ξ) n+1 ! (x ?x i )n i=0 拉格朗日插值多项式在理论分析中非常方便,因为它的结构紧凑,利用基函

数很容易推导和形象的描述算法,但是也有一些缺点,当插值节点增加、减少或其位置变化时,整个插值多项式的结构都会改变,这就不利于实际计算,增加了算法复杂度,此时我们通常采用牛顿插值多项式算法。 (2)牛顿插值多项式 牛顿插值多项式为 N(x)=f(x0)+f x0,x1(x?x0)++???+f[x0,x1,???,x n](x?x0)(x?x1)???(x?x n?1)用它插值时,首先要计算各阶差商,而各高阶差商可归结为一阶差商的逐次计算。一般情况讨论的插值多项式的节点都是任意分布的,但是在实际应用中,出现了很多等距节点的情形,这时的插值公式可以进一步简化,在牛顿均差插值多项式中各阶均差用相应的差分代替,就得到了各种形式的等距节点插值公式,常用的是牛顿前插与后插公式。 (3)分段插值 在整个插值区间上,随着插值节点的增多,插值多项式的次数必然增高,而高次插值会产生Runge现象,不能有效的逼近被插函数,人们提出用分段的低次多项式分段近似被插函数,这就是分段插值法。构造分段插值多项式的方法仍然是基函数法,即先在每个插值节点上构造分段线性插值基函数,再对基函数作线性组合。它的优点在于只要节点间距充分小,总能获得所要求的精度,即收敛性总能得到保证,另一优点是它的局部性质,即如果修改某个数据,那么插值曲线仅仅在某个局部范围内受到影响。 (4)Hermite插值 分段线性插值的算法简单,计算量小,然而从整体上看,逼近函数不够光滑,在节点处,逼近函数的左右导数不相等,若要求逼近函数与被逼近函数不仅在插值节点上取相同的函数值,而且还要求逼近函数与被逼近函数在插值节点上取相同的若干阶导数值,这类问题称为Hermite插值。 (5)样条插值 通常我们用到的分段三次埃尔米特插值构造的是一个整体上具有一阶光滑性的插值多项式,但在实际中,对光滑性的要求更高。如飞机外形的理论模型,舶体放样等型值线等常要求有二阶的光滑度。工程上常用的是3次样条函数s(x)。其基本思想是将插值区间n等分后,在每一个小区间上,采用分段3次Hermite

拉格朗日插值算法及应用实验报告

数学与统计学院课程设计报告 课程:数值分析 题目:拉格朗日插值的算法设计及应用年级:三年级 专业:数学与应用数学 学号:08063008 姓名:肖天天 指导教师:宁娣 2010年 12 月 8 日 数学与统计学院本科课程设计

拉格朗日插值的算法设计及应用 【摘要】 本文简介拉格朗日插值,它的算法及程序和拉格朗日在实际生活中的运 用。运用了拉格朗日插值的公式,以及它在MA TLAB 中的算法程序,并用具体例子说明。拉格朗日插值在很多方面都可以运用,具有很高的应用价值。 【关键词】 拉格朗日;插值;公式;算法程序;应用;科学。 一、绪论 约瑟夫·拉格朗日(Joseph Louis Lagrange),法国数学家、物理学家。他在数学、力学和天文学三个学科领域中都有历史性的贡献,其中尤以数学方面的成就最为突出。拉格朗日对流体运动的理论也有重要贡献,提出了描述流体运动的拉格朗日方法。数据建模有两大方法:一类是插值方法,另一类是拟合函数一般的说,插值法比较适合数据准确或数据量小的情形。然而Lagrange 插值有很多种,1阶,2阶,…n 阶。我们可以利用拉格朗日插值求方程,根据它的程序求原方程的图像。下面我具体介绍分析一下拉格朗日插值的算法设计及应用。 二、正文 1、基本概念 已知函数y=f(x)在若干点i x 的函数值i y =()i x f (i=0,1,???,n )一个差值问题就是求一“简单”的函数p(x):p(i x )=i y ,i=0,1,???,n, (1) 则p(x)为f(x)的插值函数,而f(x)为被插值函数会插值原函数,0x ,1x ,2x ,...,n x 为插值节点,式(1)为插值条件,如果对固定点- x 求f(- x )数值解,我们称- x 为一个插值节点,f(- x )≈p(- x )称为- x 点的插值,当- x ∈[min(0x ,1x ,2x ,...,n x ),max(0x ,1x ,2x ,...,n x )]时,称为内插,否则称为外插式外推,特别地,当p(x)为不超过n 次多项式时称为n 阶Lagrange 插值。 2、Lagrange 插值公式 (1)线性插值)1(1L 设已知0x ,1x 及0y =f(0x ) ,1y =f(1x ),)(1x L 为不超过一次多项式且满足 )(01x L =0y ,)(11x L =1y ,几何上,) (1x L 为过(0x ,0y ),(1x ,1y )的直线,从而得到 )(1x L =0y + 101x x y y --(x-0x ). (2)

几种常用的插值方法

几种常用的插值方法 数学系信息与计算科学1班平 指导老师:唐振先 摘要:插值在诸如机械加工等工程技术和数据处理等科学研究中有许多直接的应用,在很多领域都要用插值的办法找出表格和中间值,插值还是数值积分微分方程数值解等数值计算的基础。本文归纳了几种常用的插值方法,并简单分析了其各自的优缺点。 关键词:任意阶多项式插值,分段多项式插值。 引言:所谓插值,通俗地说就是在若干以知的函数值之间插入一些未知函数值,而插值函数的类型最简单的选取是代数多项式。用多项式建立插值函数的方法主要用两种:一种是任意阶的插值多项式,它主要有三种基本的插值公式:单项式,拉格朗日和牛顿插值;另一种是分段多项式插值,它有Hermite和spine插值和分段线性插值。 一.任意阶多项式插值: 1.用单项式基本插值公式进行多项式插值: 多项式插值是求通过几个已知数据点的那个n-1阶多项式,即P n-1(X)=A1+A2X+…A n X n-1,它是一个单项式基本函数X0,X1…X n-1的集合来定义多项式,由已知n个点(X,Y)构成的集合,可以使多项式通过没数据点,并为n个未知系数Ai写出n个方程,这n个方程组成的方程组的系数矩阵为Vandermonde 矩阵。 虽然这个过程直观易懂,但它都不是建立插值多项式最好的办法,因为Vandermonde方程组有可能是病态的,这样会导致单项式系数不确定。另外,单项式中的各项可能在大小上有很大的差异,这就导致了多项式计算中的舍入误差。 2.拉格朗日基本插值公式进行插值:

先构造一组插值函数L i (x ) =011011()()()() ()()()() i i n i i i i i i n x x x x x x x x x x x x x x x x -+-+--------,其中i=0,… n.容易看出n 次多项式L i (x )满足L i (x )=1,(i=j );L i (x )=0,(i ≠j ),其中i=0,1…n ,令L i (x )=0()n i i i y l x =∑这就是拉格朗日插值多项式。与单项式基本函 数插值多项式相比,拉格朗日插值有2个重要优点:首先,建立插值多项式不需要求解方程组;其次,它的估计值受舍入误差要小得多。拉格朗日插值公式结构紧凑,在理论分析中很方便,但是,当插值节点增加、减少或其位置变化时全部插值函数均要随之变化,从而整个插值公式的结构也将发生变化,这在实际计算是非常不利的。 3.使用牛顿均差插值公式进行多项式进行插值: 首先,定义均差,f 在xi,xj 上的一阶均差()()[,]j i i j j i f x f x f x x x x -=-,其中(i ≠j)。f 在 x i ,x j ,x k 的二阶均差f[x i ,x j ,x k ]= [,][,] i j j k j k f x x f x x x x --,k 阶均 f[x i …x k ]= 10[][] k i k k f x x f x x x x ---。 由此得出牛顿均值插值多项式的公式为Pn(x)=f[x 0]+f[x 0-x 1](x-x 0)+…+f[x 0,x n ](x-x 0)…(x-x n-1)。实际计算中经常利用下表给出的均差表直接构造牛顿插值公式 , , … …

几种插值法的应用和比较论文(数学类)

几种插值法的应用与比较 作者:*** 指导老师:*** 摘要本文主要介绍了几种常用插值法的应用和比较,针对每个插值法,经过详细的论证和讨论,给出了每个插值法的优点和缺点.通过对数学插值法的研究、比较及应用的讨论及总结,从而得出所讨论插值方法的各自优势,以方便用户选择合适的插值法. 关键词拉格朗日插值重心拉格朗日插值分段线性插值 1 引言 在许多实际问题及科学研究中,因素之间往往存在着函数关系,但是这些关系的显示表达式不一定都知道,通常只是由观察或测试得到一些离散数值,所以只能从这些数据构造函数的近似表达式,有时虽然给出了解析表达式,但由于解析表达式过于复杂,计算起来十分麻烦.这就需要建立函数的某种近似表达,而插值法就是构造函数的近似表达式的方法. 由于代数多项式是最简单而又便于计算的函数,所以经常采用多项式作为插值函数,称为多项式插值.多项式插值法有拉格朗日插值法,牛顿插值法、埃尔米特插值法,分段插值法和样条插值法等.其基本思想都是用高次代数多项式或分段的低次多项式作为被插值函数的近似解析表达式. 2拉格朗日插值法 在数值分析中,拉格朗日插值法是以法国十八世纪数学家约瑟夫·路易斯·拉格朗日命名的一种多项式插值方法.许多实际问题中都用函数来表示某种内在联系或规律,而不少函数都只能通过实验和观测来了解.如对实践中的某个物理量进行观测,在若干个不同的地方得到相应的观测值,拉格朗日插值法可以找到一个多项式,其恰好在各个观测的点取到观测到的值.这样的多项式称为拉格朗日(插值)多项式.数学上来说,拉格朗日插值法可以给出一个恰好穿过二维平面上若干个已知点的多项式函数.拉格朗日插值法最早被英国数学家爱德华·华林于1779年发现,不久后由莱昂哈德·欧拉再次发现.1795年,拉格朗日在其著作《师范学校数学基础教程》中发表了这个插值方法,从此他的名字就和这个方法联系在一起. 2.1 拉格朗日插值多项式

几种插值法的应用和比较教学提纲

几种插值法的应用和 比较

插值法的应用与比较 信科1302 万贤浩 13271038 1格朗日插值法 在数值分析中,拉格朗日插值法是以法国十八世纪数学家约瑟夫·路易斯·拉格朗日命名的一种多项式插值方法.许多实际问题中都用函数来表示某种内在联系或规律,而不少函数都只能通过实验和观测来了解.如对实践中的某个物理量进行观测,在若干个不同的地方得到相应的观测值,拉格朗日插值法可以找到一个多项式,其恰好在各个观测的点取到观测到的值.这样的多项式称为拉格朗日(插值)多项式.数学上来说,拉格朗日插值法可以给出一个恰好穿过二维平面上若干个已知点的多项式函数.拉格朗日插值法最早被英国数学家爱德华·华林于1779年发现,不久后由莱昂哈德·欧拉再次发现.1795年,拉格朗日在其著作《师范学校数学基础教程》中发表了这个插值方法,从此他的名字就和这个方法联系在一起. 1.1拉格朗日插值多项式 图1

已知平面上四个点:(?9, 5), (?4, 2), (?1, ?2), (7, 9),拉格朗日多项式:)(x L (黑色)穿过所有点.而每个基本多项式:)(00x l y ,)(11x l y , )(22x l y 以及)(x l y ??各穿过对应的一点,并在其它的三个点的x 值上取零. 对于给定的若1+n 个点),(00y x ,),(11y x ,………),(n n y x ,对应于它们的次数不超过n 的拉格朗日多项式L 只有一个.如果计入次数更高的多项式,则有无穷个,因为所有与L 相差))((10x x x x --λ……)(n x x -的多项式都满足条件. 对某个多项式函数,已知有给定的1+k 个取值点: ),(00y x ,……,),(k k y x , 其中i x 对应着自变量的位置,而i y 对应着函数在这个位置的取值. 假设任意两个不同的i x 都互不相同,那么应用拉格朗日插值公式所得到的拉格朗日插值多项式为: )()(0x l y x L j k j j ∑==, 其中每个)(x l j 为拉格朗日基本多项式(或称插值基函数),其表达式为: )()()()()()()()()(111100,0k j k j j j j j j j k j i i i j i j x x x x x x x x x x x x x x x x x x x x x l --------=--=++--≠=∏ΛΛ, 拉格朗日基本多项式()x l i 的特点是在j x 上取值为1,在其它的点i x ,j i ≠ 上取值为0. 例:设有某个多项式函数f ,已知它在三个点上的取值为: ? 10)4(=f , ? 25.5)5(=f , ? 1)6(=f , 要求)18(f 的值. 首先写出每个拉格朗日基本多项式:

插值法会计运用

《会计之友(下旬刊)》2008年04期 《财务管理》教学中插值法的快速理解和掌握 [09-03-2614:10:00]作者:田笑丰 摘要在时间价值及内部报酬率计算时常用到插入法,但初学者对该方法并不是很容易理解和掌握。本文根据不同情况分门别类。利用相似三角形原理推导出插入法计算用公式。并将其归纳为两类:加法公式和减法公式,简单易懂、理解准确、便于记忆、推导快捷。 关键词插入法;近似直边三角形;相似三角形 时间价值原理正确地揭示了不同时点上资金之间的换算。是财务决策的基本依据。为此,财务人员必须了解时间价值的概念和计算方法。但在教学过程中。笔者发现大多数教材插值法(也叫插入法)是用下述方法来进行的。如高等教育出版社2000年出版的《财务管理学》P62对贴现期的。 事实上,这样计算的结果是错误的。最直观的判断是:系数与期数成正向关系。而4.000更接近于3.791。那么最后的期数n应该更接近于5,而不是6。正确结果是:n=6-0.6=5.4(年)。由此可见,这种插入法比较麻烦,不小心时还容易出现上述错误。 笔者在教学实践中用公式法来进行插值法演算,效果很好,现分以下几种情况介绍其原理。 一、已知系数F和计息期n。求利息率i 这里的系数F不外乎是现值系数(如:复利现值系数PVIF年金现值系数PVIFA)和终值系数(如:复利终值系数FVIF、年金终值系数FVIFA)。 (一)已知的是现值系数那么系数与利息率(也即贴现率)之间是反向关系:贴现率越大系数反而越小,可用图1表示。 图1中。F表示根据题意计算出来的年金现值系数(复利现值系数的图示略有不同,在于i可以等于0,此时纵轴上的系数F等于1),F为在相应系数表中查到的略大于F的那个系数,F对应的利息率即为i。查表所得的另一个比F略小的系数记作F,其对应的利息率为i。 (二)已知的是终值系数 那么系数与利息率之间是正向关系:利息率越大系数也越大。其关系可用图2表示。 图2中,F表示根据题意计算出来的某种终值系数。F为在相应系数表中查到的略小于F的那个系数。F对应的利息率仍记作i,查表所得的另一个比F略大的系数记作F,其对应的利息率即为i。 上面两图中,二者往往相差1%,最多也不超过5%,故曲边三角形ABC和ADE可近

牛顿插值法原理及应用

牛顿插值法原理及应用

牛顿插值法 插值法是利用函数f (x)在某区间中若干点的函数值,作出适当的特定函数,在这些点上取已知值,在区间的其他点上用这特定函数的值作为函数f (x)的近似值。如果这特定函数是多项式,就称它为插值多项式。当插值节点增减时全部插值基函数均要随之变化,这在实际计算中很不方便。为了克服这一缺点,提出了牛顿插值。牛顿插值通过求各阶差商,递推得到的一个公式: f(x)=f[x0]+f[x0,x1](x-x0)+f[x0,x1,x2](x-x0)(x-x1)+...f[x0,...xn](x-x 0)...(x-xn-1)+Rn(x)。 插值函数 插值函数的概念及相关性质[1] 定义:设连续函数y-f(x) 在区间[a,b]上有定义,已知在n+1个互异的点 x0,x1,…xn上取值分别为y0,y1,…yn (设a≤ x1≤x2……≤xn≤b)。若在函数类中存在以简单函数P(x) ,使得P(xi)=yi,则称P(x) 为f(x)的插值函数. 称x1,x2,…xn 为插值节点,称[a,b]为插值区间。 定理:n次代数插值问题的解存在且唯一。

牛顿插值法C程序 程序框图#include void main() { float x[11],y[11][11],xx,temp,newton; int i,j,n; printf("Newton插值:\n请输入要运算的值:x="); scanf("%f",&xx); printf("请输入插值的次数(n<11):n="); scanf("%d",&n); printf("请输入%d组值:\n",n+1); for(i=0;i

常用数据插值方法的实际应用

常用数据插值方法的实际应用 任务下达书 一、实验目的 科学计算的很多问题和函数有关,但除少数代数多项式和有理分式之外,绝大多数函数都不能在计算机上直接进行求解,因此需要用计算机上可直接求解的函数插值已知函数。数据插值的基本出发点是根据函数f(x)的一组数据{x i,f(x i)}(i=0,1,…,n),按照插值条件构造简单易算的函数P(x)去近似f(x)。但对于实际问题所对应的函数,并不是插值点越多、插值多项式的次数越高,就越能获得高的精度。 通过本实验,试图达到如下目的: (1) 了解插值的基本原理和Lagrange插值、Newton插值、分段线性插值、样条插值的基本思想; (2) 验证数据插值中的Runge现象,掌握避免Runge现象的方法; (3) 掌握将数值方法转换为计算机程序的方法;掌握用高级程序设计语言实现 Lagrange插值、Newton插值、分段线性插值、样条插值的方法;培养学生上机操作和实际动手的能力; (4) 通过范例展现求解实际问题的初步建模过程;提高探索和解决问题的能力。 (5) 通过撰写实验报告,促使自己提炼思想,按逻辑顺序进行整理,并以他人能 领会的方式表达自己思想形成的过程和理由。 二、实验内容 假设两物体的运行轨迹满足如下方程: 物体1:y1 =x / (1+x2),x∈[-5,5] 物体2:y2 =1 / (1+25x2),x∈[-5,5] 其图形如图1所示。

图1 函数y1 =x / (1+x 2)和方程y2 =1 / (1+25x 2)在区间[-5,5]上的图形 但是在实际应用中,一般并不能确知物体的运动轨迹方程,只能通过测量一些x 坐标点的值,然后通过数据插值的方法来求其他未知点的函数值。设已测量到不同x 点的函数值如下表1所示,试采用Lagrange 插值、Newton 插值、分段线性插值、样条插值等方法进行插值,并绘制插值函数的图形,在对不同插值方法结果比较的基础上,确定一种最好的插值方法。 表1 不同x 点处物体的函数值 x -5.0 -4.5 -4.0 -3.5 -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0 y1 -0.1923 -0.2118 -0.2353 -0.2642 -0.3000 -0.3448 -0.4000 -0.4615 -0.5000 -0.4000 0 y2 0.0016 0.0020 0.0025 0.0033 0.0044 0.0064 0.0099 0.0175 0.0385 0.1379 1.0000 不同插值方法的结果比较可以采用如下方式:将计算结果参考表2的形式排列,对照函数的准确值,观察同一方案、不同的插值节点数n 的计算结果的变化状态,不同方案结果的精度比较,有无Runge (龙格)现象发生。取4位小数计算。 x 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 y1 0.4000 0.5000 0.4615 0.4000 0.3448 0.3000 0.2642 0.2353 0.2118 0.1923 y2 0.1379 0.0385 0.0175 0.0099 0.0064 0.0044 0.0033 0.0025 0.0020 0.0016

相关主题
相关文档
最新文档