数值分析论文-插值方法在数学分析中的应用

V数值分析中插值方法的分析与应用数值分析中插值方法的分析与应用摘要：数值分析是高等学校理工科一门重要的基础课程，主要研究数学方法的数值求解。数值分析是各种计算性科学的联系纽带和共性基础，是一门兼有基础性、应用性和边缘性的交叉学科，数值分析中插值法包括拉格朗日插值法、牛顿插值法、埃尔米特插值法等。本文主要介绍了各种插值方法的计算分析和推导，通过简单的例题进行算法分析并编程得出计算结果。关键字：数值分析；数值求解；插值法1绪论在最近的几十年中，随着计算机的发展，计算数学和应用数学中的各种方法也相应发展起来，特别是应用数学，它已经越来越渗透到其它非理工学科和各行各业中，尤其表现在生命科学、政治、军事、经济等非传统数学应用领域.同时许多教师在实践中也认识到，现有的大学数学教学内容与实际要求相去甚远.比如，几位大学计算机系毕业的学生，在面对工作中所遇见的一个非线性方程求根的问题时，他们既不知道该如何利用计算机编程求解，也不知道该如何利用计算机软件求解.某单位在LAMOS望远镜设计中，有一个复杂的概率计算问题，这个概率涉及到一个重积分，而且重积分的区问不能解析给出，负责计算的学生面对此问题感到不知所措.兴起于80年代末90年代初的数学建模比赛在一定程度上弥补了这个缺憾，参赛选手们通过参加比赛，激发了他们对数学的兴趣，也培养了他们应用数学工具解决实际问题的能力.虽然数学建模活动对学生的创造能力、应用能力有所帮助，但参加这个活动的学生毕竟是少数，这些做法并没有真正使广大学生掌握应用数学对实际问题的分析处理能力.那么,有没有这样一门课程，它既是必修课程，又具有像数学建模那样培养学生分析问题、解决问题能力的课程呢?事实上，现有的数学课程中，数值分析课程本身就具有一定的理论教学与实践的意义.数值分析是一门介绍适合于在计算机上使用的数值分析方法的课程，有时也称为计算方法课程，与其它相关数学课程相比，数值分析方法是偏重于应用的一门课程，其中的理论和方法不仅在其他专业课程中常常运用，而且在解决实际问题中也常常会用到.数值分析方法课程的基础是数学分析、线性代数、微分方程等数学理论，这些理论都为普通工科高等数学教育所覆盖，它的内容大体包括三个部分：数值逼近、数值代数、微分方程数值求解。2多种插值方法的分析比较数值分析插值法是一种古老的数学方法，它来自生产实践。利用计算机解决工程问题与常规手工计算的差异就在于它特别的计算方法.电机设计中常常需要通过查曲线、表格或通过作图来确定某一参量，如查磁化曲线、查异步电动机饱和系数曲线等.手工设计时，设计者是通过寻找坐标的方法来实现.用计算机来完成上述工作时，采用数值插值法来完成。因此学好数值分析的插值法很重2.1插值方法的定义表1插值数据表F(x₁)F(x₂)F(x₃)F(x₄)·求函数f(x)的近似表达式p(x)的方法.插值方法的必要条件是误差函数或余项R(x)=f(x)-p(x)满足关系式R(xi)=0(i=0,2,…n)当插值函数p(x)是多项式时，称为代数插值方法.代数插值方法有Lagrange插值方法，逐次线性插值法Newton插值方法，Hermite插值方法，分段插值方法和样条插值方法等.其基本思想都是用高次代数多项式或分段的低次多项式作被插函数p(x)的近似表达式。2.2常用的几种插值方法n+1个n次多项式k=0,1,……,n称为Lagrange插值基函数件L,(xk)=f(xk),k=0,1,……,n的n次多项式为Lagrange插值多项式，称为插值余项，其中ξ=ξ(x)∈(a,b)1.差商的定义f(x)关于x;的零阶差商2.Newton插值多项式设给定的n+1个互异点(x,f(x)),k=0,1,……,n,x₁≠x,,i≠j,N₇(xk)=f(xk),k=0,1,……,n的n次多项式N,(x)=f[x₀]+f[x₀,x₁](x-x₀)+……+f[x₀,x₁,…,x,](x-x₀)……(x-x,-)为Newton插值多项式，称H₂+(x;)=f,H₂n+(x;)=f,,i=0,1,…,n其中a,(x)=[1-2(x-x,)l(x2.2.4分段插值设在区间[a,b]上给定n+1个插值节点2.2.5样条插值给定这些点的函数值y,=f(x,)。②在每个区间[x;,x₄](i=0,1,2,……,n)上是3次多项式；称s(x)为3次样条插值函数。上的函数值y,=f(x,),i=0,1,2,……,n+m,构造一个有理插值R(x;)=f(x;),i=0,1,2,……,n+m则称R(x)为点集{x₀,x₁,x₂,……,x,+m-,X+m}上的有理插值函数3插值算法的编程1.Lagrange插值方法的编程以0.32,0.34,0.36为节点，分别用线性插值和抛物插值求正弦函数在0.3367处的近似值。拉格朗日插值程序清单：#include<stdio.h>doubleLagrange(doublex,intn,doublea[],doubleb[]);yoidmain()doublea[20],b[20],x,printf("请按顺序输入x对应的y的数值：\n");{/*for(i=0;i<n;i++)printf("%f",a[i]);y=Lagrange(x,n,a,b):}doubleLagrange(doublex,intn,doublea[],{L=L*(x-a[j])/(a[i]-a[j]);}}y+=L*b[i];returny;2.Newton插值方法的编程的近似值。牛顿插值程序清单：#include<stdio.h>#include<stdlib.h>#include<iostream.h>doubleNewton(doublex,intn,doublea[],doubleb[]);voidmain()printf("请输入%d个不同的x的数值：\n",n);printf("请按顺序输入x对应的y的数值：\n");{)y=Newton(x,n,a,b);}doubleNewton(doublex,intn,doublea[],doubleb[])doublet=1.0,zhi=b[0],p=0.{t=(x-a[i-1])*t;p=Chashang(0,i,a,b)*t;}ʃc=(Chashang(a1+1,b1,a,b)-Chashang(al,b1-1,a,b))/(a[b1]-a[a1]);4结论通过以上