数值计算课后复习资料

合集下载

数值计算方法 第3章复习

数值计算方法 第3章复习

1 第3章 插值法与数据拟合一、考核知识点拉格朗日插值法及其余项、牛顿插值、最小二乘法、超定方程组。

二、考核要求:1.熟练掌握拉格朗日插值法及其余项。

2.掌握牛顿插值。

3.了解最小二乘法的基本思想,熟练掌握求最小二乘多项式与超定方程组最小二乘解的方法。

三、重、难点分析例1 已知,3)9(,2)4(==f f 用线性插值计算)5(f ,并估计误差。

解 取插值节点x 0= 4,x 1= 9,两个插值基函数分别为)9(51)(1010--=--=x x x x x x l )4(51)(0101-=--=x x x x x x l 故有 565)4(53)9(52)()()(11001+=-+--=+=x x x y x l y x l x L 2.25655)5()5(1=+=≈L f 误差为 )(2)95)(45(!2)()5(2ξξf f R ''-=--''=例2 求过点(0,1)、(1,2)、(2,3)的三点插值多项式。

解:由Lagrange 插值公式又0120120,1,2;1,2,3x x x y y y ======故例3已知f(0)=8, f(1)= -7.5, f(2)= -18;用牛顿插值法求f(x)在[0,2]之间的近似零点。

0201122012010210122021()()()()()()()()()()()()()x x x x x x x x x x x x L x y y y x x x x x x x x x x x x ------=++------2(1)(2)(0)(2)(0)(1)()123(01)(02)(10)(12)(20)(21)1x x x x x x L xx ------=⨯+⨯+⨯------=+2例4求下列超定方程组的最小二乘解。

⎪⎩⎪⎨⎧=-=+=+2724212121x x x x x x1 解 令 ⎪⎩⎪⎨⎧--=-+=-+=2724213212211x x u x x u x x u23222121u u u x x ++=)(ϕ221221221)2()72()4(--+-++-+=x x x x x x 由 ⎪⎪⎩⎪⎪⎨⎧=-+=∂∂=-+=∂∂0)1662(20)1323(2212211x x x x x x ϕϕ得法方程组 ⎩⎨⎧=+=+166213232121x x x x解得 7231=x 7112=x所以最小二乘解为 7231=x 7112=x2 解 方程组写成矩阵形式为 正规方程组为即解得12114127112x x ⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎣⎦⎢⎥⎢⎥-⎣⎦⎣⎦12114111111127121121112x x ⎡⎤⎡⎤⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥--⎣⎦⎣⎦⎣⎦⎢⎥⎢⎥-⎣⎦⎣⎦1232132616x x ⎡⎤⎡⎤⎡⎤=⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦122311,77x x ==。

数值分析复习资料

数值分析复习资料

数值分析复习资料一、重点公式第一章 非线性方程和方程组的数值解法 1)二分法的基本原理,误差:~12k b ax α+--<2)迭代法收敛阶:1lim0i pi ic εε+→∞=≠,若1p =则要求01c <<3)单点迭代收敛定理:定理一:若当[],x a b ∈时,[](),x a b ϕ∈且'()1x l ϕ≤<,[],x a b ∀∈,则迭代格式收敛于唯一的根;定理二:设()x ϕ满足:①[],x a b ∈时,[](),x a b ϕ∈, ②[]121212,,, ()(),01x x a b x x l x x l ϕϕ∀∈-≤-<<有 则对任意初值[]0,x a b ∈迭代收敛,且:110111i i iii x x x llx x x lαα+-≤---≤-- 定理三:设()x ϕ在α的邻域内具有连续的一阶导数,且'()1ϕα<,则迭代格式具有局部收敛性;定理四:假设()x ϕ在根α的邻域内充分可导,则迭代格式1()i i x x ϕ+=是P 阶收敛的 ()()()0,1,,1,()0j P j P ϕαϕα==-≠ (Taylor 展开证明)4)Newton 迭代法:1'()()i i i i f x x x f x +=-,平方收敛 5)Newton 迭代法收敛定理:设()f x 在有根区间[],a b 上有二阶导数,且满足: ①:()()0f a f b <; ②:[]'()0,,f x x a b ≠∈;③:[]'',,f x a b ∈不变号④:初值[]0,x a b ∈使得''()()0f x f x <;则Newton 迭代法收敛于根α。

6)多点迭代法:1111111()()()()()()()()()i i i i i i i i i i i i i i i f x f x f x x x x x f x f x f x f x f x f x x x -+-----=-=+----收敛阶:P =7)Newton 迭代法求重根(收敛仍为线性收敛),对Newton 法进行修改 ①:已知根的重数r ,1'()()i i i i f x x x rf x +=-(平方收敛) ②:未知根的重数:1''()(),()()()i i i i u x f x x x u x u x f x +=-=,α为()f x 的重根,则α为()u x 的单根。

数值分析期末复习资料

数值分析期末复习资料

数值分析期末复习资料数值分析期末复习题型:一、填空 二、判断 三、解答(计算) 四、证明第一章误差与有效数字一、有效数字1、定义:若近似值X*的误差限是某一位的半个单位,该位到x*的第一位非零数字共有n 位,就说x*有n 位有效数字。

2、两点理解:(1) 四舍五入的一定是有效数字(2) 绝对误差不会超过末位数字的半个单位eg. ・§丄% 3、 定理1 (P6):若x*具有n 位有效数字,则其相对误差虧疗茲T 4、考点:(1)计算有效数字位数:一个根据定义理解,一个根据定理1 (P7例题3) 二、避免误差危害原则 1、原则:(1) 避免大数吃小数(方法:从小到大相加;利用韦达定理:xl*x2= c / a ) 避免相近数相减(方法:有理化)eg. V777-77 =c ・2 X2sin7 或 减少运算次数(方法:秦九韶算法)eg.P20习题14 三. 数值运算的误差估计 1、公式:(1) 一元函数:I £*( f 3))1 Q |「(於)1・| £*(力|或其变形公式求相对误差(两边同时除以f (卅))eg. P19习题1、2、5(2) (3) ln(x + £)- In x = In 1;1 — cos X =(2)多元函数(P8) eg. P8例4, P19习题4第二章插值法一、插值条件1、定义:在区间[a, b]上,给定n+1个点,aWxoVx[V・・・VxWb的函数值yi=f(xi),求次数不超过n的多项式P(x),饋兀)=儿 i =0,1,2,…,力2、定理:满足插值条件、n+1个点、点互异、多项式次数Wn的P(x)存在且唯一二、拉格朗日插值及其余项1、n次插值基函数表达式(P26 (2.8))2、插值多项式表达式(P26 (2.9))3、插值余项(P26 (2.12)):用于误差估计4、插值基函数性质(P27 (2. 17及2. 18)) eg. P28例1三、差商(均差)及牛顿插值多项式1、差商性质(P30):(1)可表示为函数值的线性组合(2)差商的对称性:差商与节点的排列次序无关(3)均差与导数的关系(P31 (3.5))2、均差表计算及牛顿插值多项式例:已知X=1,4,9的平方根为1,2,3,利用牛顿基本差商公式求"的近似值。

数值计算方法复习要点

数值计算方法复习要点

第一章引论计算方法解决问题的主要思想计算方法的精髓:以直代曲、化繁为简1、采用“构造性”方法构造性方法是指具体地把问题的计算公式构造出来。

这种方法不但证明了问题的存在性,而且有了具体的计算公式,就便于编制程序上机计算。

2、采用“离散化”方法把连续变量问题转为求离散变量问题。

例:把定积分离散成求和,把微分方程离散成差分方程。

3、采用“递推化”方法将复杂的计算过程归结为简单过程的多次重复。

由于递推算法便于编写程序,所以数值计算中常采用“递推化”方法。

4、采用“近似代替”方法计算机运算必须在有限次停止,所以数值方法常表现为一个无穷过程的截断,把一个无限过程的数学问题,转化为满足一定误差要求的有限步来近似替代。

算法的可行性分析时间复杂度、空间复杂度分析算法的复杂性(包含时间复杂性和空间复杂性)。

时间复杂度是算法耗费时间的度量。

算法的空间复杂度是指算法需占用存储空间的量度算法的可靠性分析良态算法、病态算法一个算法若运算过程中舍入误差的积累对最后计算结果影响很大,则称该算法是不稳定的或病态算法,反之称为稳定算法或良态算法。

误差的来源1、模型误差我们所建立的数学模型是对实际问题进行抽象简化而得到的。

因而总是近似的,这就产生了误差。

这种数学模型解与实际问题的解之间出现的误差,称为模型误差。

2、观测误差观测到的数据与实际数据之差。

3、截断误差数学模型的准确解与计算方法的准确解之间的误差。

4、舍入误差由于计算机字长有限,原始数据在计算机上表示会产生误差,每次计算又会产生新的误差,这种误差称为舍入误差。

绝对误差、相对误差定义2 记x*为x的近似数,称E(x)=x-x*为近似数x*的绝对误差,|E(x)|为绝对误差限。

定义3 称Er(x)=(x-x*)/x为近似数x*的相对误差。

实际运算时也将Er*(x)=(x-x*)/x*称为近似数x*的相对误差。

“四舍五入”:即尾数是4或以下则舍去,尾数是6或以上则进1,如果尾数是5,则规定:前面一位数字是偶数则舍去,奇数则进1。

(完整)数值计算方法复习

(完整)数值计算方法复习

2016计算方法复习务必通过本提纲例子和书上例子掌握如下书本内容:1. 会高斯消去法;会矩阵三角分解法;会Cholesky 分解的平方根法求解方程组2. 会用插值基函数;会求Lagrange, 会计算差商和Newton 插值多项式和余项3. 会Jacobi 迭代、Gauss —Seidel 迭代的分量形式,迭代矩阵,谱半径,收敛性4. 会写非线性方程根的Newton 迭代格式;斯蒂芬森加速5. 会用欧拉预报-校正法和经典四阶龙格—库塔法求解初值问题6. 会最小二乘法多项式拟合7. 会计算求积公式的代数精度;(复化)梯形公式和(复化)辛普生公式求积分;高斯-勒让德求积公式第1章、数值计算引论(一)考核知识点误差的来源类型;绝对误差和绝对误差限,相对误差和相对误差限,有效数字;误差的传播。

(二) 复习要求1。

了解数值分析的研究对象与特点。

2。

了解误差来源与分类,会求有效数字; 会简单误差估计. 3.了解误差的定性分析及避免误差危害。

(三)例题例1. 设x =0.231是精确值x *=0。

229的近似值,则x 有2位有效数字。

例2. 为了提高数值计算精度, 当正数x 充分大时, 应将)1ln(2--x x 改写为)1ln(2++-x x .例3. 3*x 的相对误差约是*x 的相对误差的1/3 倍.第2章、非线性方程的数值解法(一)考核知识点对分法;不动点迭代法及其收敛性;收敛速度; 迭代收敛的加速方法;埃特金加速收敛方法;Steffensen 斯特芬森迭代法;牛顿法;弦截法. (二) 复习要求1.了解求根问题和二分法.2。

了解不动点迭代法和迭代收敛性;了解收敛阶的概念和有关结论。

3。

理解掌握加速迭代收敛的埃特金方法和斯蒂芬森方法。

4。

掌握牛顿法及其收敛性、下山法, 了解重根情形. 5.了解弦截法. (三)例题1。

为求方程x 3―x 2―1=0在区间[1.3,1.6]内的一个根,把方程改写成下列形式,并建立相应的迭代公式,迭代公式不收敛的是( )(A )11,1112-=-=+k k x x x x 迭代公式 (B )21211,11kk x x x x +=+=+迭代公式(C ) 3/12123)1(,1k k x x x x +=+=+迭代公式 (D )231x x =-迭代公式11221+++=+k k kk x x x x 解:在(A)中,2/32)1(21)(,11)(,11--='-=-=x x x x x x ϕϕ2/3)16.1(21->=1.076故迭代发散。

数值计算方法复习

数值计算方法复习

数值计算方法复习1.数值求解数值求解是通过数值计算方法来寻找一个给定方程的数值解。

常见的数值求解方法包括二分法、牛顿法、割线法和迭代法等。

-二分法是一种用于求解单调函数方程的数值方法。

它将函数方程的解限定在一个区间内,然后通过缩小区间的方式来逼近解。

二分法的思想是通过不断将区间一分为二,并判断解是否在其中一半区间内,从而缩小解的范围。

-牛顿法是一种用于求解非线性方程的数值方法。

它利用函数方程的切线来逼近解。

牛顿法的核心思想是通过不断迭代逼近解的位置,使得迭代序列逐渐收敛到解。

-割线法是一种求解非线性方程的数值方法,类似于牛顿法。

它通过连结两个近似解点,得到一个割线,然后以割线和x轴的交点作为下一次迭代的近似解点。

-迭代法是一种通过迭代计算来逼近解的数值方法。

迭代法的核心思想是通过不断更新迭代序列的值,使得序列逐渐收敛到解。

2.插值与拟合插值与拟合是通过已知数据点来推断出未知数据点的数值计算方法。

-插值是通过已知数据点在这些点之间进行推断。

常见的插值方法包括拉格朗日插值和分段线性插值。

拉格朗日插值通过构造一个n次多项式函数来拟合已知数据点,从而推断出未知数据点的值。

分段线性插值是指将数据点之间的区间划分为若干段,然后在每段区间内使用线性插值来推断未知数据点的值。

-拟合是通过已知数据点在这些点之间进行逼近。

常见的拟合方法包括最小二乘拟合和多项式拟合。

最小二乘拟合通过使得残差的平方和最小来找到最优拟合函数。

多项式拟合是指通过构造一个n次多项式函数来拟合已知数据点,从而得到一个逼近函数。

3.数值积分数值积分是通过数值计算方法来近似计算函数的定积分。

常见的数值积分方法包括矩形法、梯形法、辛普森法和龙贝格法等。

-矩形法是一种通过将积分区间划分为若干个小区间,然后在每个小区间上通过函数的平均值来近似计算定积分的方法。

-梯形法是一种通过将积分区间划分为若干个小区间,然后在每个小区间上通过线性插值来近似计算定积分的方法。

数值计算方法 第5章复习

数值计算方法 第5章复习

第五章 常微分方程数值解法一、考核知识点:欧拉法,改进欧拉法,龙格-库塔法。

二、考核要求:1.熟练掌握用欧拉法,改进欧拉法求微分方程近似解的方法。

2.了解龙格-库塔法的基本思想;掌握用龙格-库塔法求微分方程近似解的方法。

3.了解稳定性。

三、重、难点分析例1 用欧拉法,预估——校正法求一阶微分方程初值问题⎩⎨⎧=-='1)0(y y x y ,在0=x (0.1)0,1近似解 解 (1)用1.0=h 欧拉法计算公式n n n n n n x y y x y y 1.09.0)(1.01+=-+=+,01n =,计算得 9.01=y 82.01.01.09.09.02=⨯+⨯=y(2)用预估—校正法计算公式1,0)(05.01.09.0)0(111)0(1=⎩⎨⎧-+-+=+=++++n y x y x y y x y y n n n n n n n n n计算得91.01=y ,83805.02=y例2、取0.1h =, 用改进欧拉法预测-校正公式求初值问题⎩⎨⎧=++='1)0(12y y x y 在0.10.2x =,处的近似值. 计算过程保留3位小数.解:改进欧拉法预测-校正公式为2n 12211111(,)(1)[(,)(,)](2)22n n n n n n n n n n n n n n n n n y y hf x y y h x y h h y y f x y f x y y x y x y ++++++⎧=+=+++⎪⎨=++=+++++⎪⎩,由h =0.1,x 0=0,y 0=1,x 1=0.1,有 ⎪⎩⎪⎨⎧=+++++==+++=227.1)2.11.0102(21.012.1)101(1.0122121y y , 由h =0.1,x 1=0.1,y 1=1.227,x 2=0.2,有 ⎪⎩⎪⎨⎧=+++++==+++=528.1)488.12.0227.11.02(21.0227.1488.1)227.11.01(1.0227.122222y y , 故,所求y (0.1)≈y 1=1.227 y (0.2)≈y 2=1.528。

数值计算方法 第6章复习

数值计算方法 第6章复习

第6章 逐次逼近法一、考核知识点:向量范数与矩阵范数及其性质,谱半径,对角占优矩阵,迭代法的收敛性,雅可比迭代法与高斯-塞德尔迭代法及其收敛性。

简单迭代法,牛顿迭代法,割线法,收敛性。

二、考核要求:1.了解向量范数的定义、性质;了解矩阵范数的定义、性质,知道谱半径的定义。

2.了解对角占优矩阵;了解迭代法的收敛性。

3.熟练掌握雅可比迭代法其收敛性的判断。

4.熟练掌握高斯-塞德尔迭代法其收敛性的判断。

5.熟练掌握用牛顿迭代法、割线法求方程近似根的方法。

了解其收敛性。

6.掌握用简单迭代法求方程的方法近似根的方法。

了解其收敛性。

三、重、难点分析例1 已知向量X=(1,-2,3),求向量X 的三种常用范数。

解 3max ==∞i ix X,14,612211====∑∑==ni ini i xXx X例2证明 ,1∞∞≤≤Xn X X证明 因为 11max X x x x Xn i i p i i=<==∑=∞∞=≤=≤∑Xn x n x n x i ip ni i max 1所以,1∞∞≤≤Xn X X例3 已知矩阵⎥⎦⎤⎢⎣⎡-=2212A ,求矩阵A 的三种常用范数。

解 4max 31==∑=∞j ij i a A,∑===ni ij ja A 114max , 39)9)(4(36135228522822122122122===--=+-=--=-⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡-⎥⎦⎤⎢⎣⎡-=λλλλλλλλA I A A A A T T例4 已知方程组⎪⎪⎪⎭⎫⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎭⎫ ⎝⎛121212212321x x x a a a (1)写出解此方程组的雅可比法迭代公式 (2)证明当4>a 时,雅可比迭代法收敛 (3)取5=a ,TX )101,51,101()0(=,求出)2(X 。

解 (1)对3,2,1=i ,从第i 个方程解出i x ,得雅可比法迭代公式为:⎪⎪⎪⎩⎪⎪⎪⎨⎧=--=--=--=+++ ,1,0,)21(1)222(1)21(1)(2)(1)1(3)(3)(1)1(2)(3)(2)1(1m x x a x x x a x x x a x m m m m m m m m m (2)当4>a 时,A 为严格对角占优矩阵,所以雅可比迭代法收敛。

数值计算方法复习题

数值计算方法复习题

fuxiti例1证明方程1-x-sin x=0在区间[0,1]内有一个根,使用二分法求误差不超过0.5×10-4的根要迭代多少次?证明令f(x)=1-x-sin x,∵f(0)=1>0,f(1)=-sin1<0∴f(x)=1-x-sin x=0在[0,1]有根.又f'(x)=1-c os x>0(x∈[0.1]),故f(x)=0在区间[0,1]内有唯一实根.给定误差限ε=0.5×10-4,有只要取n=14.例4选择填空题1. 设函数f(x)在区间[a,b]上连续,若满足,则方程f(x)=0在区间[a,b]一定有实根.答案:f(a)f(b)<0解答:因为f(x)在区间[a,b]上连续,在两端点函数值异号,由连续函数的介值定理,必存在c,使得f(c)=0,故f(x)=0一定有根.2. 用简单迭代法求方程f(x)=0的实根,把方程(x)=0表成x=ϕ(x),则f(x)=0的根是( )(A)y=x与y=ϕ(x)的交点(B) y=x与y=ϕ(x)交点的横坐标(C) y=x与x轴的交点的横坐标(D) y=ϕ(x)与x轴交点的横坐标答案:(B)解答:把f(x)=0表成x=ϕ(x), 满足x=ϕ(x)的x是方程的解,它正是y=x与y=ϕ(x)的交点的横坐标.3.为求方程x3―x2―1=0在区间[1.3,1.6]内的一个根,把方程改写成下列形式,并建立相应的迭代公式,迭代公式不收敛的是( )(A)(B)(C)(D)答案:(A)解答:在(A)中故迭代发散.在(B)中,故迭代收敛.在(C)中,,故迭代收敛.在(D)中,类似证明,迭代收敛.例3填空选择题:1. 用高斯列主元消去法解线性方程组作第1次消元后的第2,3个方程分别为。

解答1. 选a21=2为主元,作行互换,第1个方程变为:2x1+2x2+3x3=3,消元得到是应填写的内容。

一、解答下列问题:1) 数值计算中,最基础的五个误差概念(术语)是 , , , , .2) 分别用 2.718281, 2.718282 作数e 的近似值 ,它们的有效位数分别有位, 位; 又取73.13≈ (三位有效数字),则≤-73.13 .3)为减少乘除法运算次数,应将算式32)1(7)1(51318---+-+=x x x y 改写成4)为减少舍入误差的影响,应将算式 9910- 改写成 5)递推公式 ⎪⎩⎪⎨⎧=-==-,2,1,110210n y y y n n如果取41.120≈=y 作计算,则计算到10y 时,误差有这个计算公式数值稳定不稳定 ?1) 绝对误差 , 相对误差 , 有效数字 , 截断误差 , 舍入误差 。

数值计算与最优化复习1答案

数值计算与最优化复习1答案

《数算与最化》复值计优习1一.填空题1、Matlab中,清除屏幕的命令是( clc )。

2、y=100:-4:1,length(y)=( 25 )。

3、设y有5位有效数字,绝对误差限为( 0.5ⅹ10-4 )。

4、设pi的相对误差限≤0.000125,pi的有效数字位为( 4 )。

5、x i =i(i=1,2,3),f(x i )=x i ^3,f(x 1,x 2,x 3)=( 6 )。

6、用抛物型求积公式计算⎰-++114)12(dx x x 的值为( 8/3 )。

7、f(x)=x 3-6x 2+5x+1=0的Newton迭代公式为( x k+1=x k -(x k ^3-6x k ^2+5x k +1)/(3x k ^2-12x k +5) )。

8、A-qI的最大特征值为n,A的最大特征值为( n+q )。

9、Euler方法的局部离散误差为( O(h 2) )。

10、B的最小特征值为m,B -1-zI的最大特征值为( 1/m-z )。

11、用松驰迭代法解方程组,要求迭代收敛,松驰因子应满足( 0<w<2 )。

12、4次Newton-cotes求积公式的代数精度为( 5 )。

二.1、四阶Runge-Kutta方法的局部离散误差为O(h 5)。

( √ )2、非线性规划中的目标函数只能是非线性函数,约束函数可以为非线性函数。

( ⅹ )3、clear是Matlab中清内存变量命令。

( √ )4、Gauss型求积公式的代数精度不一定比Newton-cotes求积公式代数精度高。

(√ )三.计算题1、已测得直角三角形的斜边c和一直角边a的近似值为c *=5 , a *=3 ,若最大可能的误差分别是±0.15和±0.12,试求直角边所对应的角A可能的绝对误差?解:A=arcsin(a/c)2、设x 0=-1 , x 1=0 , x 2=1 , x 3=2 , Lagrange 插值,并求出x=1.5时的值,估计出误差的范围。

数值计算方法复习知识点

数值计算方法复习知识点

数值计算方法复习知识点数值计算方法是研究计算数值解的方法和数值计算的理论。

它是计算数学的一个分支,主要用于解决无法用解析方法求解的数学模型问题。

本文将综述数值计算方法的一些重要知识点,包括插值与逼近、数值微分与数值积分、线性方程组的直接解法与迭代解法以及常微分方程的数值解法。

一、插值与逼近1.插值:插值是利用已知数据点构造一个函数,使得该函数在给定的数据点上与已知函数完全相等。

常见的插值方法有拉格朗日插值和牛顿插值。

2. 逼近:逼近是从已知数据点构造一个函数,使得该函数在给定的数据点附近与已知函数近似相等。

逼近常用的方法有最小二乘逼近和Chebyshev逼近。

二、数值微分与数值积分1.数值微分:数值微分是通过计算差分商来近似计算函数的导数。

常见的数值微分方法有前向差分、后向差分和中心差分。

2.数值积分:数值积分是通过近似计算定积分的值。

常见的数值积分方法有中矩形法、梯形法和辛普森法。

三、线性方程组的直接解法与迭代解法1.直接解法:直接解法是通过一系列数学运算直接计算线性方程组的解。

常见的直接解法有高斯消元法和LU分解法。

2. 迭代解法:迭代解法是通过迭代计算逼近线性方程组的解的方法。

常见的迭代解法有Jacobi迭代法和Gauss-Seidel迭代法。

四、常微分方程的数值解法1.常微分方程:常微分方程是描述动力系统的数学模型,常用来描述物理系统、生物系统等。

常微分方程的数值解法主要包括初始值问题的一阶常微分方程和常微分方程组的数值解法。

2.常微分方程的数值解法:常微分方程的数值解法有欧拉方法、改进的欧拉方法、龙格-库塔方法等。

这些方法都是将微分方程转化为递推方程,通过迭代计算逼近微分方程的解。

总结:数值计算方法是求解数学模型的重要工具,在科学计算、工程设计和经济管理等领域有广泛的应用。

本文回顾了数值计算方法的一些重要知识点,包括插值与逼近、数值微分与数值积分、线性方程组的直接解法与迭代解法以及常微分方程的数值解法。

数值计算方法 第4章复习

数值计算方法 第4章复习

第4章 数值积分与微分一、考核知识点:梯形公式及其余项,辛卜生公式及其余项,复合梯形公式及其余项,复合辛卜生公式,高斯求积公式,代数精度,。

二、考核要求:1.了解代数精度概念,掌握插值型求积公式代数精度的判别方法。

2.熟练掌握梯形、复化梯形公式及其余项;熟练掌握辛卜生,熟练掌握运用它们计算定积分的近似值。

三、重、难点分析例1 在区间]1,1[-上,求以1,0,1321==-=x x x 为节点的插值求积公式。

解:由系数计算公式得 ⎰⎰⎰---=++==-+-+==----=11211111031)11()1(,34)10)(10()1)(1(,31)11()1(dx x x A dx x x A dx x x A 所以求积公式为)1(31)0(34)1(31)(11f f f dx x f ++-≈⎰- 例2求积公式)2(31)1(34)0(31)(20f f f dx x f ++≈⎰的代数精确度为( )。

解 由于此公式为3个节点的插值求积公式,代数精度至少为2。

令3)(x x f =,代入插值求积公式得左边=441204203==⎰x dx x ,右边4231134)0(31333=++=, 所以 左边=右边 再令4)(x x f =,代入插值求积公式得左边=532204=⎰dx x ,右边=320231134031444=++ 所以 左边≠右边 所以此公式具有3次代数精度。

例3 用梯形公式和4=n 的复合梯形公式求积分⎰+101x dx ,并估计误差。

解 (1) 梯形公式因为 1,0==b a ,11)(+=x x f ,代入梯形公式得 则75.0]111101[21)]1()0([211110=+++=+≈+⎰f f dx x (2) 复合梯形公式 因为 414=-=a b h 和复合梯形公式得 )]1())43()21()41((2)0([811110f f f f f dx x ++++≈+⎰ 697.0]21)746454(21[81≈+++⨯+= 因为 11)(+=x x f , 3)1(2)(x x f +='' , 2)(max 102=''=≤≤x f M x 所以 96116122)3(12)()(23=⨯≤''-=f n a b f R 注意:在用复合梯形公式计算 积分时注意系数的排列。

数值计算方法复习要点

数值计算方法复习要点

数值计算方法复习要点1.近似方法的概念和意义:近似方法是指通过一系列逼近计算步骤来得到问题的数值解。

在实际问题中,很多问题无法通过解析方法来求解,数值计算方法提供了一种有效的途径。

近似方法的正确性和稳定性对于数值计算方法的可靠性至关重要。

2.插值方法:插值方法是指通过已知数据点构造一个函数来逼近未知数据点的数值方法。

常见的插值方法有拉格朗日插值和牛顿插值。

在复习插值方法时,需要掌握插值多项式的构造方法和插值误差估计的技巧。

3.数值微分与数值积分:数值微分与数值积分是数值计算方法中的核心内容。

数值微分用于求取函数的导数近似值,常见的数值微分方法有差分法和微分方程法。

数值积分则是用于求取函数的积分近似值,常见的数值积分方法有梯形法则、辛普森法则和高斯积分法则。

4.非线性方程求解:非线性方程求解是数值计算方法中的重要问题之一、常见的非线性方程求解方法有二分法、牛顿法、割线法和试位法等。

在复习非线性方程求解时,要理解这些方法的基本原理和收敛性条件,并学会分析其收敛速度和稳定性。

5.线性方程组求解:线性方程组求解是数值计算方法中的另一个重要问题。

常见的线性方程组求解方法有高斯消元法、LU分解法和迭代法等。

在复习线性方程组求解时,需要理解这些方法的基本原理和收敛性条件,并学会分析其计算复杂度和稳定性。

6.数值解常微分方程:数值解常微分方程是数值计算方法的一个重要应用领域。

常见的数值解常微分方程的方法有欧拉法、改进欧拉法、龙格-库塔法等。

在复习数值解常微分方程时,需要掌握这些方法的基本原理和实现技巧,并学会分析其精度和稳定性。

8.线性插值和非线性插值:线性插值是插值方法的一种简单形式,即通过已知的两个数据点之间的线性关系来逼近未知数据点的值。

非线性插值则是通过已知的多个数据点之间的非线性关系来逼近未知数据点的值。

理解线性插值和非线性插值的原理和应用场景对于选择合适的插值方法具有重要意义。

以上是数值计算方法复习的一些重点要点,通过理解和掌握这些要点,可以为进一步深入学习和应用数值计算方法奠定基础。

数值计算方法复习题9

数值计算方法复习题9

数值计算方法复习题9习题九1. 取步长h = 0.1,分别用欧拉法与改进的欧拉法解下列初值问题(1);(2)准确解:(1);(2);欧拉法:,,,改进的欧拉法:,,,2. 用四阶标准龙格—库塔法解第1题中的初值问题,比较各法解的精度。

,,,3. 用欧拉法计算下列积分在点处的近似值。

0.5000,1.1420,2.5011,7.24504. 求下列差分格式局部截断误差的首项,并指出其阶数。

(1),2(2),3;(3),4(4),45.用Euler法解初值问题取步长h=0.1,计算到x=0.3(保留到小数点后4位).解: 直接将Eulerr法应用于本题,得到由于,直接代入计算,得到6.用改进Euler法和梯形法解初值问题取步长h=0.1,计算到x=0.5,并与准确解相比较.解:用改进Euler法求解公式,得计算结果见下表用梯形法求解公式,得解得精确解为7.证明中点公式(7.3.9)是二阶的,并求其局部截断误差主项.证明根据局部截断误差定义,得将右端Taylor展开,得故方法是二阶的,且局部截断误差主项是上式右端含h3的项。

8.用四阶R-K方法求解初值问题取步长h=0.2.解直接用四阶R-K方法其中计算结果如表所示:9.对于初值问题解因f'(y)=-100,故由绝对稳定区间要求(1)用Euler法解时,(2)用梯形法解时,绝对稳定区间为,由因f对y是线性的,故不用迭代,对h仍无限制。

(3)用四阶R-K方法时,10. (1) 用Euler法求解,步长h应取在什么范围内计算才稳定?(2) 若用梯形法求解,对步长h有无限制? (3) 若用四阶R-K方法求解,步长h如何选取?解:用四阶显式Adams公式先要算出,而,其余3点可用四阶R-K方法计算。

由,得由计算得再由四步四阶Adams显式方法得11.用四步四阶的Adams显式方法求解初值问题取h=0.1.(1)用形如的线性二步法解(2)试确定参数,使方法具有尽可能高的阶数,并求出局部截断误差主项.解本题仍利用局部截断误差的Taylor展开,要确定参数,可令解得而方法得局部截断。

《数值计算方法》复习资料

《数值计算方法》复习资料

实用文档《数值计算方法》复习资料第一章数值计算方法与误差分析第二章非线性方程的数值解法第三章线性方程组的数值解法第四章插值与曲线拟合第五章数值积分与数值微分第六章常微分方程的数值解法自测题课程的性质与任务数值计算方法是一门应用性很强的基础课,在学习高等数学,线性代数和算法语言的基础上,通过本课程的学习及上机实习、使学生正确理解有关的基本概念和理论,掌握常用的基本数值方法,培养应用计算机从事科学与工程计算的能力,为以后的学习及应用打下良好基础。

第一章数值计算方法与误差分析一考核知识点误差的来源类型;绝对误差和绝对误差限,相对误差和相对误差限,有效数字;绝对误差的传播。

二复习要求1.知道产生误差的主要来源。

2.了解绝对误差和绝对误差限、相对误差和相对误差限和有效数字等概念以及它们之间的关系。

3.知道四则运算中的误差传播公式。

实用文档三例题例 1 设x*= =3.1415926⋯近似值 x=3.14 = 0.314× 101,即 m=1,它的绝对误差是- 0.001 592 6 ,⋯有即 n=3,故 x=3.14 有 3 位有效数字 .x=3.14准确到小数点后第 2 位 .又近似值 x=3.1416,它的绝对误差是0.0000074 ⋯,有即 m=1,n= 5, x=3.1416 有 5 位有效数字 .而近似值x=3.1415,它的绝对误差是0.0000926 ⋯,有即 m=1,n= 4, x=3.1415 有 4 位有效数字 .这就是说某数有s 位数,若末位数字是四舍五入得到的,那么该数有s 位有效数字;例 2指出下列各数具有几位有效数字,及其绝对误差限和相对误差限:2.000 4-0.002 009 0009 000.00解因为 x1=2.000 4= 0.200 04× 101, 它的绝对误差限 0.000 05=0.5 × 10 1―5,即m=1,n=5, 故 x=2.000 4 有 5 位有效数字 . a1=2,相对误差限x2=- 0.002 00,绝对误差限0.000 005,因为 m=-2,n=3 ,x2=- 0.002 00 有 3 位有效数字 . a1=2 ,相对误差限r ==0.002 5实用文档x3=9 000 ,绝对误差限为0.5× 100,因为 m=4, n=4, x3=9 000 有 4 位有效数字, a=9 ,相对误差限r== 0.000 056x4=9 000.00 ,绝对误差限0.005,因为 m=4, n=6, x4=9 000.00 有 6 位有效数字,相对误差限为r== 0.000 000 56由 x3与 x4可以看到小数点之后的0,不是可有可无的,它是有实际意义的.例 3 ln2=0.69314718⋯,精确到10-3的近似值是多少?解精确到 10-3= 0.001,意旨两个近似值x1,x2满足,由于近似值都是四舍五入得到的,要求满足,近似值的绝对误差限应是=0.0005,故至少要保留小数点后三位才可以。

数值计算方法 第2章复习

数值计算方法 第2章复习

第2章 线性方程组直接解法一、考核知识点:高斯消元法,主元消元法(列主元消元法),追赶法,矩阵的三角分解。

二、考核要求:1.了解高斯消元法、主元消元法、追赶法的基本思想和使用条件2.掌握矩阵的三角分解(LU 分解,LDU 分解)、追赶法。

3.熟练掌握用列主元消元法求解线性方程组的方法。

三、重、难点分析例1 用列主元消元法解方程组。

⎪⎩⎪⎨⎧=++=++=++53368435532321321321x x x xx x x x x注意:每次消元时主元的选取是各列中系数最大的。

解 第1列主元为3,交换第1、2方程位置后消元得, ⎪⎪⎪⎩⎪⎪⎪⎨⎧=+=-=++331351313168433232321x x x x x x x第2列主35,元为交换第2、3方程位置后消元得 ⎪⎪⎪⎩⎪⎪⎪⎨⎧=-=+=++5252331356843332321x x x x x x回代解得 2,2,1123==-=x x x例2.将矩阵A 进行三角分解(LU 分解,LDU 分解)其中⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=1332222224A 说明:一般进行矩阵的三角分解采用紧凑格式。

即应用矩阵乘法和矩阵相等原则进行矩阵的三角分解(或代入公式求得相应元素)。

在分解时注意矩阵乘法、矩阵求逆等代数运算。

解: 9,2;1,121,21;2,2,43322123132321321232312212222113131112121131312121111=-=-=-=-==-=-====-======r r r l a l r l a r r l a r a a l a a l a r a r a r则矩阵的LU 分解为 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡--=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----911224122112111332222224 因为对角阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=914D ,则⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡--==-111212111R D U 所以矩阵的LDU 分解为 ⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡--⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡--=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----11121211914122112111332222224。

数值计算复习资料

数值计算复习资料

第一章 绪论§1 绪论:数值分析的研究内容 §2 误差的来源和分类 §3 误差的表示 §4 误差的传播 §5 算法设计的若干原则一、误差的分类(绝对误差,相对误差)例1-1 设 x *=是由精确值x 经过四舍五入得到的近似值。

问 x 的绝对误差限ε和相对误差限η各是多少解:因为 x =x * ± ,所以绝对误差限为ε= 相对误差限为二、有效数字定义 设数 x 的近似值可以表示为其中 m 是整数,αi (i=1,2, …, n ) 是0到9 中的一个数字,而α1 ≠ 0. 如果其绝对误差限为则称近似数 x* 具有 n 位有效数字。

结论:通过四舍五入原则求得的近似数,其有效数字就是从末尾到第一位非零数字之间的所有数字。

例1-2 下列近似数是通过四舍五入的方法得到的,试判定它们各有几位有效数字:x 1* =87540,x 2*=8754×10, x 3*=, x 4*= ×10-2%23.018.2005.0*≈==x εηm n x 10.021*⨯±=ααα *1102m nx x --≤⨯111x x *-≤5511101-*⨯≤-x x 所以已知有5位有效数字。

同理可以写出可以得出 x 2 , x 3 , x 4 各具有4、3、4 位有效数字。

例1-3 已知 e =……, 试判断下面两个近似数各有几位有效数字解:由于而e 1有7位有效数字。

同理:e 2 只有6位有效数字。

三、算法设计的若干原则• 1:两个很接近的数字不做减法:• 2: 不用很小得数做分母(不用很大的数做分子)练习: 求方程 x 2-56x +1=0 的两个根,使它们至少具有四位有效数字•第二章 插值与拟合1、Lagrange 插值多项式,Newton 插值多项式的构造与插值余项估计,及证明过程。

2、 Hermite 插值多项式的构造与插值余项估计,带导数条件的插值多项式的构造方法,基于承袭性的算法,基函数法, 重节点差商而所以 7161102110210000005.00000001.0--⨯=⨯=≤=- e e 510.8754010x *=⨯而1221102x x *-≤⨯520.875410x *=⨯54221102x x *--≤⨯5331102x x *--≤⨯230.34510x *-=⨯-23331102x x *--≤⨯6441102x x *--≤⨯240.345010x *-=⨯24441102x x *---≤⨯718281.2,718282.221==e e 6110210000005.00000001.0-⨯=≤=- e e 11102718282.0718282.2⨯==e 615210211021000005.00000008.0--⨯=⨯=<=- e e表的构造;3、分段插值及三次样条插值的构造4、最小二乘拟合• 掌握Lagrange 插值多项式的构造方法及具体结构 • 掌握Lagrange 插值多项式误差分析方法和证明方法 • 掌握Newton 插值多项式的形式及误差 • 掌握差商表的构造过程 关于离散数据:Newton 插值多项式:例1-3 已知f (x ) 的五组数据(1,0)、(2,2)、(3,12)、(4,42)、(5,116),求 N 4 (x )。

数值分析计算方法复习

数值分析计算方法复习
矩阵的三角分解定理:设A为n阶矩阵,若A的顺序主 子式Di ≠0(i=1,2,…,n-1),则A可分解为一个单位下三 角矩阵L和一个上三角矩阵U的乘积,且这种分解是唯 一的。
Gauss列主元素消元法
基本思想:在系数矩阵A的第一列中选绝对值最大的 元素作主元,设该元素所在行为第i1行,交换第一 行与第i1行,进行第一次消元;再在第2-n行的第 二列中选绝对值最大的元素作主元,设该元素所在 行为第i2行,交换第二行与第i2行,进行第二次消 元,……直到消元过程完成为止。
则对任意x0 [a, b],由xn+1=(xn )得到的迭代序列{xn } 收敛到(x)的不动点x *,并有误差估计:
|
xn
x*
|
1
1
L
|
xn1
xn
|
Ln xn x * 1 L x1 x0
局部收敛性定理 设x*为(x)的不动点, `(x)在x*的 某邻域连续,且|`(x*)|<1,则不动点迭代法xk+1=(xk )
若对任意 x0[a, b],由上述迭代得序列{xk},有极限
lim
k
xk
x*
则称迭代过程收敛,且x*=(x*)为(x)的不动点。
(二)不动点的存在性和迭代法的收敛性
▪ 不动点的存在性定理 设(x)C[a, b]且满足以下两
个条件:
(1)对于任意x [a, b],有a≤(x)≤b; (2)若(x)在[a, b]一阶连续,且存在常数0<L<1,
第二章 非线性方程的求根方法
▪ *引言 ▪ 方程求根的二分法 ▪ 迭代法及其收敛性
不动点迭代法 不动点的存在性与迭代法的收敛性 迭代收敛的加速方法
▪ Newton迭代法:
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

习 题 一 解 答1.取3.14,3.15,227,355113作为π的近似值,求各自的绝对误差,相对误差和有效数字的位数。

分析:求绝对误差的方法是按定义直接计算。

求相对误差的一般方法是先求出绝对误差再按定义式计算。

注意,不应先求相对误差再求绝对误差。

有效数字位数可以根据定义来求,即先由绝对误差确定近似数的绝对误差不超过那一位的半个单位,再确定有效数的末位是哪一位,进一步确定有效数字和有效数位。

有了定理2后,可以根据定理2更规范地解答。

根据定理2,首先要将数值转化为科学记数形式,然后解答。

解:(1)绝对误差:e(x)=π-3.14=3.14159265…-3.14=0.00159…≈0.0016。

相对误差:3()0.0016()0.51103.14r e x e x x -==≈⨯ 有效数字:因为π=3.14159265…=0.314159265…×10,3.14=0.314×10,m=1。

而π-3.14=3.14159265…-3.14=0.00159…所以│π-3.14│=0.00159…≤0.005=0.5×10-2=21311101022--⨯=⨯ 所以,3.14作为π的近似值有3个有效数字。

(2)绝对误差:e(x)=π-3.15=3.14159265…-3.14=-0.008407…≈-0.0085。

相对误差:2()0.0085()0.27103.15r e x e x x --==≈-⨯ 有效数字:因为π=3.14159265…=0.314159265…×10,3.15=0.315×10,m=1。

而π-3.15=3.14159265…-3.15=-0.008407…所以│π-3.15│=0.008407……≤0.05=0.5×10-1=11211101022--⨯=⨯ 所以,3.15作为π的近似值有2个有效数字。

(3)绝对误差:22() 3.14159265 3.1428571430.0012644930.00137e x π=-=-=-≈-L L 相对误差:3()0.0013()0.4110227r e x e x x--==≈-⨯ 有效数字:因为π=3.14159265...=0.314159265 (10)22 3.1428571430.3142857143107==⨯,m=1。

而22 3.14159265 3.1428571430.0012644937π-=-=-L L 所以221322 3.14159265 3.1428571430.0012644930.0057110.510101022π----=-=≤=⨯=⨯=⨯L L 所以,227作为π的近似值有3个有效数字。

(4)绝对误差:355() 3.14159265 3.141592920.00000027050.000000271113e x π=-=-=-≈-L L 相对误差:7()0.000000271()0.86310355113r e x e x x--==≈-⨯ 有效数字:因为π=3.14159265...=0.314159265 (10)355 3.141592920.31415929210113==⨯,m=1。

而355 3.14159265 3.141592920.0000002705113π-=-=-L L 所以6617355 3.14159265 3.141592920.00000027050.0000005113110.510101022π----=-=≤=⨯=⨯=⨯L L 所以,355113作为π的近似值有7个有效数字。

指出:①实际上,本题所求得只能是绝对误差限和相对误差限,而不是绝对误差和相对误差。

2、用四舍五入原则写出下列各数的具有五位有效数字的近似数。

346.7854,7.000009,0.0001324580,0.600300解:346.7854≈346.79,7.000009≈7.0000,0.0001324580≈0.00013246,0.600300≈0.60030。

指出:注意0。

只要求写出不要求变形。

3、下列各数都是对准确数进行四舍五入后得到的近似数,试分别指出他们的绝对误差限和相对误差限和有效数字的位数。

12340.0315,0.3015,31.50,5000x x x x ====。

分析:首先,本题的准确数未知,因此绝对误差限根据四舍五入规则确定。

其次,应当先求绝对误差限,再求相对误差限,最后确定有效数字个数。

有效数字由定义可以直接得出。

解:由四舍五入的概念,上述各数的绝对误差限分别是1234()0.00005,()0.00005,()0.005,()0.5x x x x εεεε====由绝对误差和相对误差的关系,相对误差限分别是111222333444()0.00005()0.16%,0.0315()0.00005()0.02%,0.3015()0.005()0.002%,31.5()0.5()0.01%.5000x x x x x x x x x x x x εδεδεδεδ==≈==≈==≈==≈ 有效数字分别有3位、4位、4位、4位。

指出:本题显然是直接指出有效数位、直接写出绝对误差,用定义求出相对误差。

4.0.1%。

解:设取n 个有效数字可使相对误差小于0.1%,则 111100.1%2n a -⨯<,而34≤≤,显然13a =,此时, 1111110100.1%223n n a --⨯=⨯<⨯, 即13110106n --⨯<, 也即461010n ⨯>所以,n=4。

3.162≈。

5、在计算机数系F(10,4,-77,77)中,对31120.14281100.31415910x x =⨯=-⨯与,试求它们的机器浮点数()(1,2)i fl x i =及其相对误差。

解:3333111111112222()0.142810,(())()0.14281100.1428100.0000110,()0.314210,(())()0.31415910(0.314210)0.0004110fl x e fl x x fl x fl x e fl x x fl x =⨯=-=⨯-⨯=⨯=-⨯=-=-⨯--⨯=⨯其相对误差分别是3112310.00001100.000041100.007%,0.013%0.1428100.314210e e ⨯⨯=≈=≈-⨯-⨯。

6、在机器数系F(10,8,L,U)中,取三个数4220.2337125810,0.3367842910,0.3367781110x y z -=⨯=⨯=-⨯,试按(),()x y z x y z ++++两种算法计算x y z ++的值,并将结果与精确结果比较。

解:422222222(())(0.23371258100.3367842910)0.3367781110(0.00000023100.3367842910)0.33677811100.33678452100.33677811100.0000064110fl x y z -++=⨯+⨯-⨯=⨯+⨯-⨯=⨯-⨯=⨯ 42242222(())0.2337125810(0.33678429100.3367781110)0.23371258100.00000618100.00000023100.00000618100.0000064110fl x y z --++=⨯+⨯-⨯=⨯+⨯=⨯+⨯=⨯精确计算得:4222222220.23371258100.33678429100.3367781110(0.00000023371258100.3367842910)0.33677811100.33678452371258100.33677811100.000064137125810x y z -++=⨯+⨯-⨯=⨯+⨯-⨯=⨯-⨯=⨯ 第一种算法按从小到大计算,但出现了两个数量级相差较大的数相加,容易出现大数吃小数.而第二种算法则出现了两个相近的数相减,容易导致有效数位的减少。

计算结果证明,两者精度水平是相同的。

***在机器数系F(10,8,L,U)中,取三个数4220.2337125810,0.3367842910,0.3367781110x y z --=⨯=⨯=-⨯,试按(),()x y z x y z ++++两种算法计算x y z ++的值,并将结果与精确结果比较。

解:42222222222(())(0.23371258100.3367842910)0.3367781110(0.00233713100.3367842910)0.33677811100.33912142100.33677811100.00003391100.33677811100.336744210fl x y z -----++=⨯+⨯-⨯=⨯+⨯-⨯=⨯-⨯=⨯-⨯=-⨯ 42242242222(())0.2337125810(0.33678429100.3367781110)0.2337125810(0.00003368100.3367781110)0.23371258100.33674742100.00000023100.33674742100.3367471910fl x y z ----++=⨯+⨯-⨯=⨯+⨯-⨯=⨯-⨯=⨯-⨯=-⨯第一种算法是按从小到大的顺序计算的,防止了大数吃小数,计算更精确。

精确计算得:42220.23371258100.33678429100.33677811100.0000233712580.003367842933.6778110.00339121415833.67781133.6744197858420.3367441978584210x y z --++=⨯+⨯-⨯=+-=-=-=-⨯显然,也是第一种算法求出的结果和精确结果更接近。

7、某计算机的机器数系为F(10,2,L,U),用浮点运算分别从左到右计算及从右到左计算10.40.30.20.040.030.020.01+++++++试比较所得结果。

解:从左到右计算得10.40.30.20.040.030.020.010.1100.04100.03100.02100.00100.00100.00100.00100.19101.9+++++++=⨯+⨯+⨯+⨯+⨯+⨯+⨯+⨯=⨯=从右到左计算得111110.40.30.20.040.030.020.010.010.020.030.040.20.30.410.1100.2100.3100.4100.20.30.410.10.20.30.410.11010.1100.1100.2102----+++++++=+++++++=⨯+⨯+⨯+⨯++++=++++=⨯+=⨯+⨯=⨯= 从右到左计算避免了大数吃小数,比从左到右计算精确。

相关文档
最新文档