拉格朗日插值法理论及误差分析

合集下载

《拉格朗日插值法》课件

《拉格朗日插值法》课件

在数值分析中的应用
数值积分
01
拉格朗日插值法可用于数值积分,通过插值多项式逼近被积函
数,进而求得积分的近似值。
数值微分
02
利用拉格朗日插值法可以近似求得函数的导数值,用于数值微
分计算。
求解常微分方程
03
通过构造插值多项式,可以将常微分方程转化为代数方程组,
进而求解微分方程的近似解。
在数据拟合中的应用
重要性
拉格朗日插值法是数值分析中的基础方法之一,它为解决各种实际问题提供了重要的数学工具。通过 拉格朗日插值法,我们可以更好地理解和逼近数据,从而为进一步的数值分析和科学计算提供基础。
拉格朗日插值法的历史和发展
历史
拉格朗日插值法由意大利数学家约瑟夫·拉格朗日于18世纪提出。在此之前,人们已经意识到可以通过已知的数 据点来逼近未知的函数值,但缺乏系统的数学方法。拉格朗日的插值法为这个问题提供了一个完整的解决方案, 并在随后的几个世纪中得到了广泛的应用和发展。
深入研究拉格朗日多项式的性质
拉格朗日多项式是拉格朗日插值法的基础,但其性质仍有许多未知之处。未来的研究可以深入探索拉格 朗日多项式的性质,以期为插值法的发展提供新的思路和方法。
THANK YOU
多项式形式
插值多项式的一般形式为 (L(x) = sum_{i=0}^{n} y_i l_i(x)),其中 (l_i(x)) 是拉格朗日插值基函数。
求解插值多项式的系数
系数求解
通过已知的数据点和构造的插值多项 式,求解出多项式的系数。
求解方法
常用的求解方法是高斯消元法或追赶 法,通过求解线性方程组得到插值多 项式的系数。
《拉格朗日插值法》ppt课件
• 引言 • 拉格朗日插值法的基本概念 • 拉格朗日插值法的实现步骤 • 拉格朗日插值法的优缺点分析 • 拉格朗日插值法的应用实例 • 总结与展望

拉格朗日插值法理论及误差分析

拉格朗日插值法理论及误差分析

拉格朗日插值法理论及误差分析拉格朗日插值法理论及误差分析浅析拉格朗日插值法目录:一、引言二、插值及多项式插值的介绍三、拉格朗日插值的理论及实验四、拉格朗日插值多项式的截断误差及实用估计式五、1、截断误差在[a,b]区间上用Ln(x)近似未知或复杂函数f(x),其截断误差是指Rn?x??f?x??Ln?x?通常称Rn?x?为拉格朗日插值余额。

注意到利用公式估计截断误差实际上非常困难。

一是因为它要计算函数f(x)的高阶导数,当f(x)很复杂时,计算量很大,而当f(x)没有可用来计算的表达式时,导数无法准确计算;二是因为即使能得到高阶导数的解析式,但于?的具体位置不知道,所以要估计高阶导数在插值区间上的界一般是非常困难的事情。

因此,公式并不实用。

2、截断误差的实用估计式既然公式估计误差时不实用,那么实际中如何估计截断误差呢?假设插值条件中包含n+2组数据?,n,n?, 1f(xi)?yi , i?0,1那么利用n+1组数据我们可以构造一个n 次拉格朗日插值多项式Ln(x),利用后n+1组数据我们可以构造另一个n次拉格朗日插值多项式L*n(x)。

利用公式知,他们各自的插值余项为f(x)?Ln(x)?1f(n?1)(?)(x?x0)(x?x1)?(x?xn),(n?1)!1f(n?1)(?*)(x?x1)(x?x2)?(x?xn?1), (n?1)!f(x)?L*n(x)?两式相减得L*n(x)?Ln(x)?并可写成1fn?1(?)(x?x1)?(x?xn)(xn?1?x0),(n?1)!L*(x)?Ln(x)1(n?1)f(?)(x?x1)?(x?xn)?n.(n?1)!xn?1?x0注意到上式中利用fn?1(?)?fn?1(?*).该条件在很多情况下是成立的。

利用式可得?Ln(x)?L*n(x)R(x)?f(x)?L(x)?,n? nx0?xn?1? ? *?R*(x)?f(x)?L*(x)?Ln(x)?Ln(x),nn?xn?1x0?式给出了用Ln(x)或L*n(x)作近似计算时的实用误差估计式,它不需要计算高阶导数,也不用估计插值区间上高阶导数的界。

拉格朗日 插值 区间误差限

拉格朗日 插值 区间误差限

拉格朗日插值区间误差限拉格朗日插值方法是一种常用的数值插值方法,用于在给定一组已知数据点的情况下,通过构造一个多项式函数来拟合这些数据点,并在插值区间内求得未知值。

然而,由于插值方法的近似性质,插值结果与真实值之间总会存在一定的误差。

本文将介绍拉格朗日插值法以及其误差限的计算方法。

一、拉格朗日插值法简介拉格朗日插值法是一种基于多项式的插值方法,其基本思想是通过构造一个满足给定数据点的插值多项式来逼近真实的函数曲线。

具体而言,对于给定的n个数据点(xi, yi),拉格朗日插值法的插值多项式可以表示为:P(x) = Σ[ yi * Li(x) ],i=0 to n其中,Li(x)是拉格朗日基函数,定义为:Li(x) = Π[ (x - xj) / (xi - xj) ],j=0 to n,i ≠ j这样,通过求解插值多项式P(x),我们可以在插值区间内求得未知值。

二、插值误差限的计算尽管拉格朗日插值法可以通过构造插值多项式来逼近真实函数曲线,但由于插值方法本质上是一种近似方法,插值结果与真实值之间总会存在一定的误差。

我们可以通过计算插值误差限来评估插值的可靠性。

在拉格朗日插值法中,插值误差限可通过以下等式进行估计:| f(x) - P(x) | ≤ M / (n + 1)! * | x - x0 | * | x - x1 | * ... * | x - xn |其中,f(x)是真实函数的值,P(x)是插值多项式的值,M是插值区间上函数f(x)的最大导数的上界,n是插值多项式的次数。

三、拉格朗日插值法的应用示例为了更好地理解拉格朗日插值法及其误差限的计算方法,我们来看一个具体的示例。

假设我们要通过拉格朗日插值法来估计函数f(x) = sin(x)在区间[0, π]内的某个未知值。

已知在该区间内取了n+1个等间距的数据点(xi, yi),其中i=0, 1, 2, ..., n。

首先,我们可以根据已知数据点构造拉格朗日插值多项式P(x),并计算出未知值的近似值。

拉格朗日插值多项式与泰勒多项式的误差研究详全文

拉格朗日插值多项式与泰勒多项式的误差研究详全文

拉格朗日插值多項i. 式與泰勒多項式的誤差分析朱亮儒★ 曾政清☆ 陳昭地★★國立臺灣師範大學數學系教授☆臺北市立建國高級中學數學教師摘要:本文旨於提供拉格朗日插值多項式與泰勒多項式誤差項估計值的初等簡易證明,並探討其應用價值。

關鍵字:拉格朗日插值多項式、泰勒多項式、誤差項一引言有鑑於教育部99普通高級中學數學課綱在第一冊多項式的運算為迴避解三元一次方程組,首次出現插值多項式及其應用(以不超過三次插值多項式為限><[1][2][3]),99數學課綱包含插值多項式部分如下:求中的.除以的餘式為通過的插值多項式。

若有兩實根,則可寫成的型式。

透過因式定理證明插值多項式的唯一性。

設通過的多項式為,求及.插值多項式:通過的多項式可表示為,求的值。

此處暫不處理下面的題型:「設通過的多項式為,求。

」此類題型將在數學的IV的聯立方程組章節中處理。

此處自然而然讓人想到拉格朗日(Lagrange, J. L.,1736-1816>其人奇事,羅列如下:他出生於義大利西北部的杜林(Turin>,從小就極有數學天分,於18歲開始撰寫數學論文,在數論上曾提出一個著名的定理:「任意正整數都可以表成四個平方數的和」。

他是第一位證明均值定理(The Mean Value Theorem>的大數學家。

(均值定理在高三選修甲微分的單元中會學到<[4]),它是僅次於微積分基本定理的極重要的存在定理>他在30歲時,應腓特烈二世的邀請到柏林作為其宮廷數學大師長達20年之久。

之後接受法國的邀請,到巴黎擔任法國科學院院士,拿破崙<1769-1821, 1804-1815擔任法皇)讚譽他為「數學科學的巍峨金字塔」泰勒定理有拉格朗日誤差的公式<存在性)。

拉格朗日恆等式:,,.具有附加條件的多變數實函數極值拉格朗日乘子定理。

最得意的巨著《分析力學》。

拉格朗日差值誤差公式<[5]):若為區間中相異實數,且,則對每一個,存在,使得,其中為函數在的階拉格朗日插值多項式,而為其插值誤差式。

拉格朗日插值法 (2)

拉格朗日插值法 (2)

在数值分析中,拉格朗日插值法是以法国十八世纪数学家约瑟夫·拉格朗日命名的一种多项式插值方法。

许多实际问题中都用函数来表示某种内在联系或规律,而不少函数都只能通过实验和观测来了解。

如对实践中的某个物理量进行观测,在若干个不同的地方得到相应的观测值,拉格朗日插值法可以找到一个多项式,其恰好在各个观测的点取到观测到的值。

这样的多项式称为拉格朗日(插值)多项式。

数学上来说,拉格朗日插值法可以给出一个恰好穿过二维平面上若干个已知点的多项式函数。

拉格朗日插值法最早被英国数学家爱德华·华林于1779年发现[1],不久后(1783年)由莱昂哈德·欧拉再次发现。

1795年,拉格朗日在其著作《师范学校数学基础教程》中发表了这个插值方法,从此他的名字就和这个方法联系在一起[2]。

对于给定的若n+1个点,对应于它们的次数不超过n 的拉格朗日多项式只有一个。

如果计入次数更高的多项式,则有无穷个,因为所有与相差的多项式都满足条件。

定义对某个多项式函数,已知有给定的k + 1个取值点:其中对应着自变量的位置,而对应着函数在这个位置的取值。

假设任意两个不同的x j都互不相同,那么应用拉格朗日插值公式所得到的拉格朗日插值多项式为:其中每个为拉格朗日基本多项式(或称插值基函数),其表达式为:[3]拉格朗日基本多项式的特点是在上取值为1,在其它的点上取值为0。

存在性对于给定的k+1个点:,拉格朗日插值法的思路是找到一个在一点取值为1,而在其他点取值都是0的多项式。

这样,多项式在点取值为,而在其他点取值都是0。

而多项式就可以满足在其它点取值为0的多项式容易找到,例如:它在点取值为:。

由于已经假定两两互不相同,因此上面的取值不等于0。

于是,将多项式除以这个取值,就得到一个满足“在取值为1,而在其他点取值都是0的多项式”:这就是拉格朗日基本多项式。

唯一性次数不超过k的拉格朗日多项式至多只有一个,因为对任意两个次数不超过k的拉格朗日多项式:和,它们的差在所有k+1个点上取值都是0,因此必然是多项式的倍数。

插值数值实验报告(3篇)

插值数值实验报告(3篇)

第1篇一、实验目的1. 理解并掌握插值法的基本原理和常用方法。

2. 学习使用拉格朗日插值法、牛顿插值法等数值插值方法进行函数逼近。

3. 分析不同插值方法的优缺点,并比较其精度和效率。

4. 通过实验加深对数值分析理论的理解和应用。

二、实验原理插值法是一种通过已知数据点来构造近似函数的方法。

它广泛应用于科学计算、工程设计和数据分析等领域。

常用的插值方法包括拉格朗日插值法、牛顿插值法、样条插值法等。

1. 拉格朗日插值法拉格朗日插值法是一种基于多项式的插值方法。

其基本思想是:给定一组数据点,构造一个次数不超过n的多项式,使得该多项式在这些数据点上的函数值与已知数据点的函数值相等。

2. 牛顿插值法牛顿插值法是一种基于插值多项式的差商的插值方法。

其基本思想是:给定一组数据点,构造一个次数不超过n的多项式,使得该多项式在这些数据点上的函数值与已知数据点的函数值相等,并且满足一定的差商条件。

三、实验内容1. 拉格朗日插值法(1)给定一组数据点,如:$$\begin{align}x_0 &= 0, & y_0 &= 1, \\x_1 &= 1, & y_1 &= 4, \\x_2 &= 2, & y_2 &= 9, \\x_3 &= 3, & y_3 &= 16.\end{align}$$(2)根据拉格朗日插值公式,构造插值多项式:$$P(x) = \frac{(x-x_1)(x-x_2)(x-x_3)}{(x_0-x_1)(x_0-x_2)(x_0-x_3)}y_0 + \frac{(x-x_0)(x-x_2)(x-x_3)}{(x_1-x_0)(x_1-x_2)(x_1-x_3)}y_1 + \frac{(x-x_0)(x-x_1)(x-x_3)}{(x_2-x_0)(x_2-x_1)(x_2-x_3)}y_2 + \frac{(x-x_0)(x-x_1)(x-x_2)}{(x_3-x_0)(x_3-x_1)(x_3-x_2)}y_3.$$(3)计算插值多项式在不同点的函数值,并与实际值进行比较。

拉格朗日插值法知识讲解

拉格朗日插值法知识讲解

拉格朗日插值法5.2 拉格朗日(Lagrange)插值可对插值函数选择多种不同的函数类型,由于代数多项式具有简单和一些良好的特性,例如,多项式是无穷光滑的,容易计算它的导数和积分,故常选用代数多项式作为插值函数。

5.2.1 线性插值问题5.1给定两个插值点其中,怎样做通过这两点的一次插值函数?过两点作一条直线,这条直线就是通过这两点的一次多项式插值函数,简称线性插值。

如图5.1所示。

图5.1 线性插值函数在初等数学中,可用两点式、点斜式或截距式构造通过两点的一条直线。

下面先用待定系数法构造插值直线。

设直线方程为,将分别代入直线方程得:当时,因,所以方程组有解,而且解是唯一的。

这也表明,平面上两个点,有且仅有一条直线通过。

用待定系数法构造插值多项式的方法简单直观,容易看到解的存在性和惟一性,但要解一个方程组才能得到插值函数的系数,因工作量较大和不便向高阶推广,故这种构造方法通常不宜采用。

当时,若用两点式表示这条直线,则有:(5.1)这种形式称为拉格朗日插值多项式。

,,称为插值基函数,计算,的值,易见(5.2)在拉格朗日插值多项式中可将看做两条直线,的叠加,并可看到两个插值点的作用和地位都是平等的。

拉格朗日插值多项式型式免除了解方程组的计算,易于向高次插值多项式型式推广。

线性插值误差定理5.1记为以为插值点的插值函数,。

这里,设一阶连续可导,在上存在,则对任意给定的,至少存在一点,使(5.3)证明令,因是的根,所以可设对任何一个固定的点,引进辅助函数:则。

由定义可得,这样至少有3个零点,不失一般性,假定,分别在和上应用洛尔定理,可知在每个区间至少存在一个零点,不妨记为和,即和,对在上应用洛尔定理,得到在上至少有一个零点,。

现在对求二次导数,其中的线性函数),故有代入,得所以即5.2.2 二次插值问题5.2给定三个插值点,,其中互不相等,怎样构造函数的二次的(抛物线)插值多项式?平面上的三个点能确定一条次曲线,如图5.2所示。

拉格朗日插值性质

拉格朗日插值性质

拉格朗日插值性质
拉格朗日插值是一种多项式插值,它是在已知数据点的插值方法,它的特点是可以在给定的多个数据点上构造一个多项式,使该多项式在这些数据点上取得极值。

拉格朗日插值是基于拉格朗日插值多项式的概念,它是一个多项式,它满足在给定的n个数据点上,其在这些数据点上的值都是相同的,而在其他任意点上的值都是不同的。

拉格朗日插值的优点是拟合效果好,可以很好的拟合离散的数据点,而且拉格朗日插值多项式的曲线也比较平滑,不易出现抖动现象。

缺点是计算量大,拉格朗日插值多项式的构造需要计算大量的系数,所以计算量比较大,而且容易出现过拟合现象。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

浅析拉格朗日插值法目录:一、 引言二、 插值及多项式插值的介绍 三、 拉格朗日插值的理论及实验四、 拉格朗日插值多项式的截断误差及实用估计式 五、 参考文献一、引言插值在数学发展史上是个古老问题。

插值是和拉格朗日(Lagrange )、牛顿(Newton )、高斯(Gauss )等著名数学家的名字连在一起的。

在科学研究和日常生活中,常常会遇到计算函数值等一类问题。

插值法有很丰富的历史渊源,它最初来源人们对天体研究——有若干观测点(我们称为节点)计算任意时刻星球的位置(插值点和插值)。

现在,人们在诸如机械加工等工程技术和数据处理等科研都有很好的应用,最常见的应用就是气象预报。

插值理论和方法能解决在实际中当许多函数表达式未知或形式复杂,如何去构造近似表达式及求得在其他节点处的值的问题。

二、插值及多项式插值1、插值问题的描述设已知某函数关系()y f x =在某些离散点上的函数值:插值问题:根据这些已知数据来构造函数()y f x =的一种简单的近似表达式,以便于计算点,0,1,,i x x i n ≠=的函数值()f x ,或计算函数的一阶、二阶导数值。

xx 0y y1y 1n y -ny 1x 1n x -nx2、插值的几何意义插值的几何意义如图1所示:图1 3、多项式插值 3.1 基本概念假设()y f x =是定义在区间,a b ⎡⎤⎣⎦上的未知或复杂函数,但一直该函数在点01n a x x x b ≤<<<≤处的函数值01,,n y y y 。

找一个简单的函数,例如函数()P x ,使之满足条件(),0,1,2,,,i P x y i n == (3.1)通常把上述01n x x x <<< 称为插值节点,把()P x 称为()f x 的插值多项式,条件(3.1)称为插值条件,并把求()P x 的过程称为插值法。

3.2 插值多项式的存在性和唯一性 如果插值函数是如下m 次的多项式:1011()m m m m m P x a x a x a x a --=+++那么插值函数的构造就是要确定()m P x 表达式中的m+1个系数011,,,m ma a a a -。

由于插值条件包含n+1独立式,只要m=n 就可证明插值函数多项式是唯一存在。

实际上,由n+1个插值条件可得10010100101111111011n n n n n n n n n n nn n n n na x a x a x a y a x a x a x a y a x a x a x a y ------⎧+++=⎪+++=⎪⎨⎪⎪+++=⎩这是一个关于01,,n a a a 的n+1阶线性方程组,且其系数矩阵对应的行列式是线性代数中著名的范德蒙(Vandemonde )行列式。

该行列式得值为 0110(,,)()n in n i j i j V x x x x x ===-∏∏因为i j ≠时,i j x x ≠,所以01(,,)0n n V x x x ≠。

从而证明了上述线性方程组的阶是唯一存在的。

既满足插值条件的多项式唯一存在。

三、 拉格朗日插值的理论及实验 1、拉格朗日插值的理论拉格朗日(Lagrange )插值公式的基本思想是把()n P x 的构造问题转化为n+1个插值基函数()(0,1,,)i l x i n =。

首先我们利用节点直接构造如下多项式:1'1()()()()in n i n x l x x x x ππ++=- 其中101()()()()n n x x x x x x x π+=---,'1011()()()()()n i i i i i i n x x x x x x x x x π+-+=----容易验证该多项式具有性质0,1,i j i l j i ≠⎧=⎨=⎩因此,n 次多项式00110()()()()()nn n n k k k L x l x y l x y l x y l x y ==++=∑一定具有性质()()(),0,1,,,nn i k k i i i k L x l x y l x y i n ====∑既满足插值条件。

我们称()n L x 为拉格朗日插值多项式,()i l x 称为拉格朗日插值及函数。

一次拉格朗日插值多项式又叫做线性插值多项式。

二次拉格朗日插值多项式又叫做抛物线插值多项式。

2、拉格朗日插值实验经过学习掌握拉格朗日插值的理论,学以致用,使学到的知识运用到现实生活中,并运用计算机来解决我们在学习中遇到的一些问题。

以下为运用MATLAB 软件平台上计行拉格朗日插值问题:例:已知在[0,30]内对应的节点x 以及函数值y 如表所示,利用拉格朗日插值多项式求在区间x=2.035,x=9.771,x=17.815,x=26.907所对应的函数值。

在已知数表函数的条件下,拉格朗日插值多项式可用来计算复杂函数或未知函数的函数值,为此我们首先编写如下利用拉格朗日插值多项式方法计算函数值的程序:function y=lagrange(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~=kp=p*(z-x0(j))/(x0(k)-x0(j)); end ends=p*y0(k)+s; end y(i)=s; end上述三重循环给出了拉格朗日插值计算多项式计算任何点x 处的函数值的过程,我们把它标记为lagrange.m文件,接下来我们在MATLAB平台上进行上述例子中的数值试验。

在Command Window中输入的命令及结果如下所示:>> x=[0:2:30];>> y=[0.0 0.45 0.75 1.52 2.12 3.40 3.72 4.17 5.12 5.45 5.67 6.74 7.31 7.85 8.45 8.97];>> lagrange(x,y,2.035)ans =0.3290>> lagrange(x,y,9.771)ans =3.2975>> lagrange(x,y,17.815)ans =5.4483>> lagrange(x,y,26.907)ans =8.6519最后,我们根据拉格朗日插值结果,利用plot命令画出未知函数的图像,命令程序如下:>> x0=[0:2:30];>> y0=lagrange(x,y,x0);>> plot(x0,y0)得到的未知函数图像为:四、拉格朗日插值多项式的截断误差及实用估计式1、截断误差在[a,b]区间上用()n L x 近似未知或复杂函数()f x ,其截断误差是指()()()n n R x f x L x =- (4.1) 通常称()n R x 为拉格朗日插值余额。

注意到利用公式(4.1)估计截断误差实际上非常困难。

一是因为它要计算函数()f x 的高阶导数,当()f x 很复杂时,计算量很大,而当()f x 没有可用来计算的表达式时,导数无法准确计算;二是因为即使能得到高阶导数的解析式,但由于ξ的具体位置不知道,所以要估计高阶导数在插值区间上的界一般是非常困难的事情。

因此,公式(4.1)并不实用。

2、截断误差的实用估计式既然公式(4.1)估计误差时不实用,那么实际中如何估计截断误差呢? 假设插值条件中包含n+2组数据(),i i f x y = 0,1,,,1,i n n =+那么利用n+1组数据我们可以构造一个n 次拉格朗日插值多项式()nL x ,利用后n+1组数据我们可以构造另一个n 次拉格朗日插值多项式*()n L x 。

利用公式(4.1)知,他们各自的插值余项为(1)011()()()()()(),(1)!n n n f x L x f x x x x x x n ξ+-=---+*(1)*1211()()()()()(),(1)!n n n f x L x f x x x x x x n ξ++-=---+两式相减得*11101()()()()()(),(1)!n n n n n L x L x f x x x x x x n ξ++-≈---+并可写成*(1)110()()1()()().(1)!n n n n n L x L x f x x x x n x x ξ++---≈+- (4.2)注意到上式中利用11*()().n n f f ξξ++≈该条件在很多情况下是成立的。

利用式(4.2)可得*01***10()()()()(),()()()()(),n n nn n n n n n n L x L x R x f x L x x x L x L x R x f x L x x x ++⎧-=-≈⎪-⎪⎨-⎪=-≈⎪-⎩ (4.2) 式(4.3)给出了用()n L x 或*()n L x 作近似计算时的实用误差估计式,它不需要计算高阶导数,也不用估计插值区间上高阶导数的界。

总之,拉格朗日插值法的公式结构紧凑,在理论分析中十分方便,然而在计算中,但插值点增加或减少时,所对应的基本多项式就得重新计算而且图像发生很大变化。

像逐次线性插值法、牛顿插值法等都是在拉格朗日插值多项式的基础上延伸出来的。

我们根据实际中的具体问题,为减少插值误差来选取相应的插值法来快速的解决问题。

五、参考文献[1] 数值计算原理 李庆扬,关冶,白峰杉 清华大学出版社 2000.9 [2] 数值分析及其MATLAB 实现 任玉杰 高等教育出版社 7007.3 [3] 数值分析与实验 韩旭里,万中 科学出版社 2006.7[4] 数值分析(第三版)颜庆津 北京航空航天大学出版社 2011.8.15。

相关文档
最新文档