郑州大学数值分析重点考察内容及各章习题

合集下载

数值分析课后习题及答案

数值分析课后习题及答案

第一章 绪论(12) 第二章 插值法(40-42)2、当2,1,1-=x 时,4,3,0)(-=x f ,求)(x f 的二次插值多项式。

[解]372365)1(34)23(21)12)(12()1)(1(4)21)(11()2)(1()3()21)(11()2)(1(0))(())(())(())(())(())(()(2221202102210120120102102-+=-++--=+-+-⨯+------⨯-+-+-+⨯=----+----+----=x x x x x x x x x x x x x x x x x x x y x x x x x x x x y x x x x x x x x y x L 。

3、给出x x f ln )(=的数值表用线性插值及二次插值计算54.0ln 的近似值。

X 0.4 0.5 0.6 0.7 0.8 x ln -0.916291 -0.693147 -0.510826 -0.357765 -0.223144[解]若取5.00=x ,6.01=x ,则693147.0)5.0()(00-===f x f y ,510826.0)6.0()(11-===f x f y ,则604752.182321.1)5.0(10826.5)6.0(93147.65.06.05.0510826.06.05.06.0693147.0)(010110101-=---=--⨯---⨯-=--+--=x x x x x x x x x y x x x x y x L ,从而6202186.0604752.19845334.0604752.154.082321.1)54.0(1-=-=-⨯=L 。

若取4.00=x ,5.01=x ,6.02=x ,则916291.0)4.0()(00-===f x f y ,693147.0)5.0()(11-===f x f y ,510826.0)6.0()(22-===f x f y ,则 217097.2068475.404115.2)2.09.0(5413.25)24.0(3147.69)3.01.1(81455.45)5.06.0)(4.06.0()5.0)(4.0()510826.0()6.05.0)(4.05.0()6.0)(4.0()693147.0()6.04.0)(5.04.0()6.0)(5.0(916291.0))(())(())(())(())(())(()(22221202102210120120102102-+-=+--+-⨯++-⨯-=----⨯-+----⨯-+----⨯-=----+----+----=x x x x x x x x x x x x x x x x x x x x x x y x x x x x x x x y x x x x x x x x y x L ,从而61531984.0217097.21969765.259519934.0217097.254.0068475.454.004115.2)54.0(22-=-+-=-⨯+⨯-=L补充题:1、令00=x ,11=x ,写出x e x y -=)(的一次插值多项式)(1x L ,并估计插值余项。

数值分析第7章答案

数值分析第7章答案

第七章非线性方程求根一、重点内容提要 (一)问题简介 求单变量函数方程(7.1) 的根是指求(实数或复数),使得.称为方程(7.1)的根,也称为函数的零点.若可以分解为其中m 为正整数,满足,则是方程(7.1)的根.当m=1时,称为单根;当m>1时,称为m 重根.若充分光滑,是方程(7.1)的m 重根,则有(1)()(*)'(*)...(*)0,(*)0m m f x f x f x f x -====≠ 若在[a,b]上连续且,则方程(7.1)在(a,b)内至少有一个实根,称[a,b]为方程(7.1)的有根区间.有根区间可通过函数作图法或逐次搜索法求得. (二)方程求根的几种常用方法 1.二分法设在[a,b]上连续,,则在(a,b)内有根.再设在(a,b)内仅有一个根.令,计算和.若则,结束计算;若,则令,得新的有根区间;若,则令,得新的有根区间.,.再令计算,同上法得出新的有根区间,如此反复进行,可得一有根区间套且110011*,0,1,2,...,()...()22n n n n n n a x b n b a b a b a --<<=-=-==-.故因此,可作为的近似根,且有误差估计 (7.2) 2.迭代法将方程式(7.1)等价变形为 (7.3)若要求满足则;反之亦然.称为函数的一个不动点.求方程(7.1)的根等价于求的不动点由式(7.3)产生的不动点迭代关系式(也称简单迭代法)为 (7.4)函数称为迭代函数.如果对任意,由式(7.4)产生的序列有极限 则称不动点迭代法(7.4)收敛.定理7.1(不动点存在性定理)设满足以下两个条件: 1.对任意有2.存在正常数,使对任意,都有 (7.5) 则在上存在惟一的不动点.定理7.2(不动点迭代法的全局收敛性定理)设满足定理7.1中的两个条件,则对任意,由(7.4)式得到的迭代序列收敛.到的不动点,并有误差估计式 (7.6) 和 (7.7)定理7.3(不动点迭代法的局部收敛性定理)设为的不动点,在的某个邻域连续,且,则迭代法(7.4)局部收敛.收敛阶的概念 设迭代过程(7.4)收敛于方程的根,如果迭代误差当时成产下列渐近关系式(7.8)则称该迭代过程是p 阶收敛的.特别地,p=1时称线性收敛,p>1时称超线性收敛,p=2时称平方收敛.定理7.4(收敛阶定理)对于迭代过程(7.4),如果在所求根的邻近连续,并且 (7.9)则该迭代过程在点的邻近是收敛的,并有(7.10)斯蒂芬森(Steffensen)迭代法 当不动点迭代法(7.4)只有线性收敛阶,甚至于不收敛时,可用斯蒂芬森迭代法进行加速.具体公式为 (7.11) 此法也可写成如下不动点迭代式(7.12)定理7.5(斯蒂芬森迭代收敛定理) 设为式(7.12)中的不动点,则是的不动点;设存在,,则是的不动点,则斯蒂芬森迭代法(7.11)是2阶收敛的. 3.牛顿迭代法牛顿迭代法是一种特殊的不动点迭代法,其计算公式为 其迭代函数为 (7.13)牛顿迭代法的收敛速度 当时,容易证明,,,由定理7.4知,牛顿迭代法是平方收敛的,且(7.14)重根情形的牛顿迭代法 当是的m 重根时,迭代函数在处的导数,且.所以牛顿迭代法求重根只是线性收敛.若的重数m 知道,则迭代式 (7.15)求重根二阶收敛.当m 未知时,一定是函数的单重零点,此时迭代式1()()'()'()['()]()''()0,1,2,...k k k k k k k k k k x f x f x x x x x f x f x f x k μμ+=-=--= (7.16)也是二阶收敛的.简化牛顿法 如下迭代法 称为简化牛顿法或平行弦法.牛顿下山法 为防止迭代不收敛,可采用牛顿下山法.具体方法见教材. 4.弦截法将牛顿迭代法(7.13)中的用在,处的一阶差商来代替,即可得弦截法 (7.17)定理7.6假设在其零点的邻域内具有二阶连续导数,且对任意有,又初值,,则当邻域充分小时,弦截法(7.17)将按阶收敛到.这里p 是方程的正根. 5.抛物线法弦截法可以理解为用过两点的直线方程的根近似替的根.若已知的三个近似根,,用过的抛物线方程的根近似代替的根,所得的迭代法称为抛物线法,也称密勒(Muller)法.当在的邻近有三阶连续导数,,则抛物线法局部收敛,且收敛阶为.二、知识结构图10[1,2]1x x --=≤≤--∈3-3-6k k 32三、常考题型及典型题精解例7-1 证明方程x 在上有一个实根x*,并用二分法求这个根,要求|x -x*|10.若要求|x -x*|10,需二分区间[1,2]多少次?解 设f(x)=x ,则f(1)=-1<0,f(2)=5>0,故方程f(x)=0在[1,2]上有根x*.又因f'(x)=3x -1,所以当x [1,2]时,f'(x)>0,即f (x)=0在[1,2]上有惟一实根x*.用二分法计算结果如表7-1所示.k 0 1 2 3 4 5 6 7 8 9 1 1 1.25 1.25 1.3125 1.3125 1.3125 1.3204 1.3243 1.3243 2 1.5 1.5 1.375 1.375 1.13438 1.3282 1.3282 1.3282 1.32631.5 1.25 1.375 1.3125 1.3438 1.3282 1.3204 1.3243 1.3263 1.3253+ - + - + + - - + +610x e -≤≤⨯≤≤≤≤≥∈-3-39910-6k k k+101此时x =1.3253满足|x -x*|0.9771010,可以作为x*的近2似值.1若要求|x -x*|,只需|x -x*|10即可,解得k+119.932,2即只需把[1,2]二分20次就能满足精度要求.例7-2 已知函数方程(x-2)=1,(1)确定有根区间[a,b];(2)构造不动点迭代公式使之对任意初始近似x [a,b],31|10.k x ---<k 迭代方法均收敛;(3)用所构造的公式计算根的近似值,要求|x1lim lim x x x x x e e e e →+∞→-∞∞∞∞∈解 (1)令f(x)=(x-2)-1,由于f(2)=-1<0,f(3)=-1>0,因此区间[2,3]是方程f(x)=0的一个有根区间.又因f'(x)=(x-1),f(x)=+,f(x)=-1,f'(1)=--1<0,当x>1时f(x)单增,x<1时f(x)单减,故f(x)=0在(-,+)内有且仅有一根x*,即x*[2,3].2'k k x x x x x x e e e e e e e ϕϕϕ-----∈∈≤≤≤∀∈k+100k+1(2)将(x-2)=1等价变形为x=2+,x [2,3].则(x)=2+.由于当x [2,3]时2(x)3,|(x)|=|-|<1故不动点迭代法x =2+,k=0,1,2,...,对x [2,3]均收敛.(3)取x =2.5,利用x =2+进行迭代计算,结果如表7-2所示.473cos 3120cos c k x x x ϕ--+=∈≤4k+10-30k+1k+1k 例 考虑求解方程2的迭代公式2x =4+,k=0,1,2,...3(1)试证:对任意初始近似x R,该方法收敛;(2)取x =4,求根的近似值x ,要求|x -x |10;(3)所给方法的收敛阶是多少?2解 (1)由迭代公式知,迭代函数(x)=4+3{}os ,(,).|'sin |1(,)x x x ϕϕϕ∈-∞+∞≤<-∞+∞∀∈0k 022由于(x)的值域介于(4-)与(4+)之间,且3322(x)|=|-33故根据定理7.1,7.2知,(x)在内存在惟一的不动点x*,且对x R,迭代公式得到的序列x 收敛于x*.(2) 取x =4,迭代计算结果如表7-3所示.此时已满足误差要求,即(3)由于,故根据定理7 .4知方法是线性收敛的,并且有。

数值分析课程重点知识点以及课后重点习题

数值分析课程重点知识点以及课后重点习题

第一章:
知识点:绝对误差,相对误差,有效数字
重点习题:英文版:1、2
中文版:1、2、3
第二章:
知识点:向量、矩阵的范数,条件数
重点习题:英文版:4
中文版:7
第三章:
知识点:Jacobi迭代,G-S迭代(要会写出迭代的形式),收敛性,||M||<1的使用
重点习题:英文版:6
中文版:5
第四章:
知识点:Newton法,弦截法
重点习题:英文版:5、8
中文版:11
第五章:
知识点:Legrange插值公式,(),Newton插值公式(构建
差商表),三次样条,
重点习题:英文版:1、3、5
中文版:1、5、10
第六章:
知识点:最小二乘法(只考一次方程,使用其原始公式或者内积形式公式)
重点习题:英文版:1
中文版:1
第七章:
知识点:数值微分(三点公式——一阶导数、二阶导数情况下的几种不同的公式),数值积分,代数精度,梯形公式,Simpson公式,Guass公式,
重点习题:英文版:1、2
中文版:1、2、6
第八章:
知识点:Eular方法、改进的Eular方法(改进方法为重点,并要求知道证明Eular方法是一阶方法)
重点习题:英文版:1、2
中文版:1、4。

数值分析课后习题与解答

数值分析课后习题与解答

课后习题解答第一章绪论习题一1.设x>0,x*的相对误差为δ,求f(x)=ln x的误差限。

解:求lnx的误差极限就是求f(x)=lnx的误差限,由公式(已知x*的相对误差满足,而,故即2.下列各数都是经过四舍五入得到的近似值,试指出它们有几位有效数字,并给出其误差限与相对误差限。

解:直接根据定义和式(有5位有效数字,其误差限,相对误差限有2位有效数字,有5位有效数字,3.下列公式如何才比较准确?(1)(2)解:要使计算较准确,主要是避免两相近数相减,故应变换所给公式。

(1)(2)4.近似数x*=0.0310,是 3 位有数数字。

5.计算取,利用:式计算误差最小。

四个选项:第二、三章插值与函数逼近习题二、三1. 给定的数值表用线性插值与二次插值计算ln0.54的近似值并估计误差限. 解:仍可使用n=1及n=2的Lagrange插值或Newton插值,并应用误差估计(5.8)。

线性插值时,用0.5及0.6两点,用Newton插值误差限,因,故二次插值时,用0.5,0.6,0.7三点,作二次Newton插值误差限,故2. 在-4≤x≤4上给出的等距节点函数表,若用二次插值法求的近似值,要使误差不超过,函数表的步长h 应取多少?解:用误差估计式(5.8),令因得3. 若,求和.解:由均差与导数关系于是4. 若互异,求的值,这里p≤n+1.解:,由均差对称性可知当有而当P=n+1时于是得5. 求证.解:解:只要按差分定义直接展开得6. 已知的函数表求出三次Newton均差插值多项式,计算f(0.23)的近似值并用均差的余项表达式估计误差.解:根据给定函数表构造均差表由式(5.14)当n=3时得Newton均差插值多项式N3(x)=1.0067x+0.08367x(x-0.2)+0.17400x(x-0.2)(x-0.3) 由此可得f(0.23) N3(0.23)=0.23203由余项表达式(5.15)可得由于7. 给定f(x)=cosx的函数表用Newton等距插值公式计算cos 0.048及cos 0.566的近似值并估计误差解:先构造差分表计算,用n=4得Newton前插公式误差估计由公式(5.17)得其中计算时用Newton后插公式(5.18)误差估计由公式(5.19)得这里仍为0.5658.求一个次数不高于四次的多项式p(x),使它满足解:这种题目可以有很多方法去做,但应以简单为宜。

数值分析 第七章 习题

数值分析 第七章 习题

郑州大学研究生课程(2010-2011学年上学期) 数值分析Numerical Analysis- 1 -第七章 习 题1. 应用对分法求三次方程324100x x +−=在区间[1,2]上的实根,并做误差估计. 2. 若应用对分法求方程sin02xxe π−=在区间[0,1]上误差不超过512的近似根,应对分几次?3. 求方程324100x x +−=在区间[1,2]上的根,并构造如下2个迭代格式:(1) 121110()()4k k kx x x ϕ+==+(2) 3212()410k k k k k x x x x x ϕ+==−−+试判定迭代格式的收敛性,对收敛的格式求近似根,并给出误差估计. 4. 试用迭代的方法证明函数()xf x e −=满足积分方程0()1()xf x f t dt =−∫.5.用不动点迭代法计算s =的近似值.6. 确定求230xx e −=正根的不动点迭代的收敛区间[],a b ,并求出满足4110k k x x −+−<的近似值1k x x +∗=,若要求近似值的误差410ε−≤,应迭代几步?7.6位有效数字.8. 用双点割线法求方程32210200x x x ++−=的根,要求6110k k x x −+−<. 9. 设x a =是()x a ϕ=在区间[],a b 上的不动点,试证:(1)当[]1,C a b ϕ∈且'()0a ϕ≠时,不动点迭代1()n n x x ϕ+=有线性收敛速度;(2)当[]2,C a b ϕ∈且'()0a ϕ=而''()0a ϕ≠,不动点迭代1()n n x x ϕ+=至少有平方收敛速度.10. 设方程1232cos 0x x −+=有迭代式124cos 3k k x x +=+(1)试证:对任意初值0x ,迭代序列{}n x 收敛;(2)取迭代初值04x =,求该方程在误差不超过310−的近似根; (3)此迭代法的收敛阶是多少?证明所得结论.。

郑州大学研究生课程数值分析复习---第四章数值微分与数值积分

郑州大学研究生课程数值分析复习---第四章数值微分与数值积分

郑州大学研究生课程数值分析复习---第四章数值微分与数值积分郑州大学研究生课程(2012-2013学年第一学期)数值分析Numerical Analysis习题课第四章数值微分与数值积分4.2.1 显格式导数与差商的关系数值微分取为导数的近似值,即差商。

00()() ()()()()()2lim lim lim h h h f x h f x h f x f x h f x h f x h f x h h →→→+??′=??+[,]()/i a b n x a ih h b a n =+=?将区间等分为份,是等距节点,是步长。

一、要点回顾()()'()i i i f x h f x f x h+?≈向前差商x ix i +h一阶数值微分由Taylor 展开2()()'()''(),2!f x h f x hf x f x x hξξ+=++≤≤+因此,有误差()()()'()''()()2!i i i f x h f x hR x f x f O h hξ+?=?=?=一阶数值微分()()'()i i i f x f x h f x h≈向后差商x i -hx i一阶数值微分由Taylor 展开2()()'()''(),2!i i i i i hf x h f x hf x f x x hξξ?=?+≤≤+因此,有误差()()()'()''()()2!i i i f x f x h hR x f x f O h h ξ??=?==一阶数值微分()()'()2i i i f x h f x h f x h+??≈中心差商x i -hx i一阶数值微分由Taylor 展开23112322()()'()''()'''(),2!3!()()'()''()'''(),2!3!i i i i i i i i i i i ih hf x h f x hf x f x f x x hh hf x h f x hf x f x f x h x ξξξξ+=+++≤≤+?=?+??≤≤因此,有误差22212()()()'()2 ['''()'''()]'''()()126i i i f x h f x h R x f x hh hf f f O h ξξξ+??=?=+==一阶数值微分()()231000232000120021()()'()''(),2!()()2'()2''(),24()()3()'(),2()4()3()'().2n n n n h f x f x hf x f x O h f x f x hf x h f x O h h f x f x f x f x hf x f x f x f x h=+++=+++??≈?+≈将第一式乘4并减去第二式,除以可得类似可得一阶数值微分数值积分公式求积系数求积节点()()()()nbk k n ak I f f x dx A f x I f =≈∑∫数值积分0()()nn k k k I f A f x ==∑0()()()()(),nbn k k ak R f I f I f f x dx A f x ==?=?∑∫分别称为为数值求积公式和求积公式余项数值积分求积公式的代数精度定义1称求积公式具有m 次代数精度,如果它满足如下两个条件:(i )对所有次数≤m 次的多项式,有(ii )存在m+1次多项式,使得)(x P m 0)()()(=?=m n m m P I P I P R )(1x P m +0)()()(111≠?=+++m n m m P I P I P R 代数精度定义1中的条件(i),(ii)等价于: )()()0(,0)()()()(1≠≤≤== + mknkkxR iimkxxIxR i代数精度012,,,,,()1,,,,.n nx x x n f x x x x =L L 对于给定的一组结点要构造至少有次代数精确度的求积公式则它对于精确成立基本目标代数精度012,,,,n A A A A L 即求积公式的系数满足线性方程组+?=+++?=+++?=+++++1211110022110010n a b x A x A x A a b x A x A x A a b A A A n n nn n n n n n n L M L L 代数精度Newton-Cotes数值积分插值型求积公式上取一组节点在积分区间],[b a bx x x a n ≤<<<≤L 10插值多项式次的作Lagrange n x f )(0()()()nn k k k x f x l x ?==∑为插值基函数),,1,0)((n k x l k L =()(),n x f x ?用作为被积函数的近似有badx x f )(()bn ax dx ?≈∫∫∑==b ank kkdxx l x f 0)()(∑∫==nk bak k dxx l x f 0)()(则,记∫=bak k dx x l A )(∫badx x f )(∑=≈nk k k x f A 0)(Newton-Cotes数值积分()(0,1,,)bk k aA l x dx k n ==∫L 定义:系数由式所确定的求积公式称为插值型求积公式.4.5.11.n n +定理:利用个结点的求积公式至少具有次代数精确度的充分必要条件是它是插值型的Newton-Cotes数值积分等距节点的Newton-Cotes 求积公式],[)(b a C x f ∈设函数插值多项式为的Lagrange x f )([,]a b n 将积分区间分割为等份,nk kh a x k ,,1,0,L =+=为步长其中nab h ?=各节点为0()()()nn k k k x f x l x ?==∑Newton-Cotes数值积分。

数值分析习题(含标准答案)

数值分析习题(含标准答案)

]第一章 绪论姓名 学号 班级习题主要考察点:有效数字的计算、计算方法的比较选择、误差和误差限的计算。

1若误差限为5105.0-⨯,那么近似数有几位有效数字(有效数字的计算) 解:2*103400.0-⨯=x ,325*10211021---⨯=⨯≤-x x 故具有3位有效数字。

2 14159.3=π具有4位有效数字的近似值是多少(有效数字的计算) 解:10314159.0⨯= π,欲使其近似值*π具有4位有效数字,必需!41*1021-⨯≤-ππ,3*310211021--⨯+≤≤⨯-πππ,即14209.314109.3*≤≤π即取( , )之间的任意数,都具有4位有效数字。

3已知2031.1=a ,978.0=b 是经过四舍五入后得到的近似值,问b a +,b a ⨯有几位有效数字(有效数字的计算)解:3*1021-⨯≤-aa ,2*1021-⨯≤-b b ,而1811.2=+b a ,1766.1=⨯b a 2123****102110211021)()(---⨯≤⨯+⨯≤-+-≤+-+b b a a b a b a故b a +至少具有2位有效数字。

2123*****10210065.01022031.1102978.0)()(---⨯≤=⨯+⨯≤-+-≤-b b a a a b b a ab 故b a ⨯至少具有2位有效数字。

4设0>x ,x 的相对误差为δ,求x ln 的误差和相对误差(误差的计算)~解:已知δ=-**xx x ,则误差为 δ=-=-***ln ln xx x x x则相对误差为******ln ln 1ln ln ln xxx x xxx x δ=-=-5测得某圆柱体高度h 的值为cm h 20*=,底面半径r 的值为cm r 5*=,已知cm h h 2.0||*≤-,cm r r 1.0||*≤-,求圆柱体体积h r v2π=的绝对误差限与相对误差限。

(误差限的计算)解:*2******2),(),(h h r r r h r r h v r h v -+-≤-ππ绝对误差限为πππ252.051.02052)5,20(),(2=⨯⋅+⨯⋅⋅⋅≤-v r h v相对误差限为%420120525)5,20()5,20(),(2==⋅⋅≤-ππv v r h v 6设x 的相对误差为%a ,求nx y =的相对误差。

郑州大学研究生课程数值分析复习---第五章 解线性方程组的直接法

郑州大学研究生课程数值分析复习---第五章 解线性方程组的直接法

郑州大学研究生课程(2012-2013学年第一学期)数值分析Numerical Analysis习题课第五章解线性代数方程组的直接法第k 步:消去第k 列依此类推,直到第n-1 步,原方程化为()()(1,...,)k k ik ikkkm ak ai n +=−=设,计算()0k kka≠(1)(1)(1)(1)1112111(2)(2)(2)22222()()n n n n n n nn a a a b x x a a b x b a ⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦L L M M O M 计算(1)()()(1)()()k k k ij ij ik kjk k k i i ik ka a m ab b m b ++=+=+( i ,j = k +1, …, n )一、要点回顾§5.2 Gauss消去法回代过程算法abx )n (nn)n (nn =()()()1() i n 1 , n 2 ,, 1.ni i i iiijjiij i x ba x a=+=−=−−∑L (1)(1)(1)(1)111111()()()(1)(1)(1)11111 iinni i i ii i in n i n n n n n n n n n n a x a x a x ba x a xb a x a x b −−−−−−−−++++=+++=+=L L L L L L L L()() n n nn n n a x b ⎧⎪⎪⎪⎪⎨⎪⎪⎪⎪=⎪⎩乘除运算量:由于计算机中做乘除运算的时间远远超过做加减运算时间,故估计运算量时,往往只估计乘除的次数。

第k 步:消去第k 列()()(1,...,)k k ik ik kkm a k a i n +=−=设,计算()0k kk a ≠计算(1)()()(1)()()k k k ij ij ik kj k k k i i ik k a a m a b b m b ++=+=+( i = k +1, …, n )回代求解:()()n n n n nnx b a =()()()1()ni i i i i ij jiij i x b a x a =+=−∑( i = k +1, …, n )n –k 次(n –k )2次n –k 次n (n+1)/2 次高斯消去法总的乘除运算量为:3233n n n +−§5.2 高斯消去法定理5.2.1是矩阵的顺序主子式A ,0111≠=a D 11110(1,2,,).ii i iia a D i k a a =≠=L MM L L高斯约化的主元素的充要条件()0(1,2,,)i ii a i k ≠=L 0(1,2,,)i D i k ≠=L定理5.2.2若矩阵A 对称正定,则()()01,2,,k kk a k n ≠=L 推论5.2.1如果的顺序主子式A ),1,,2,1(0−=≠n k D k L ⎪⎩⎪⎨⎧=,1)1(11D a 则).,,3,2(/1)(n k D D a k k k kk L ==−记上三角矩阵为,)(n A U (1)111121,n A MMMU L U −−−−==L 为单位下三角矩阵.2131321231111n n n m m m mm m ⎛⎞⎜⎟−⎜⎟⎜⎟=−−⎜⎟⎜⎟⎜⎟−−−⎝⎠M M M O L111121n L M M M −−−−=L 其中为上三角阵.()n AU =三角(LU)分解§5.4 三角分解法§5.4 三角分解法多利特尔(Doolittle)分解定理5.4.1n阶非奇方阵A有唯一的Doolittle分解的充要条件是A的前n-1个顺序主子式0(1,2,..., 1.)k D k n ≠=−§5.4 三角分解法n , 2,i , n , 1,j , 111i111L L ====u a l a u i j j n, , 1k i )/ ( n , ,k j, ,3 ,2 kk1-k 1m imik1-k 1m km kj ∑∑L L L +=−==−====uul alu l a u mkikmj kj n k 计算对例5.4.1用多利特尔分解求解方程组:解设A=LU ,即12312312324211326132.225122x x x x x x A x x x ++=⎧⎛⎞⎪⎜⎟++==⎨⎜⎟⎪++=⎩⎝⎠111213212223313233211100132100,122100u u u l u u l l u ⎛⎞⎛⎞⎛⎞⎜⎟⎜⎟⎜⎟=⎜⎟⎜⎟⎜⎟⎝⎠⎝⎠⎝⎠100211 0.5100 2.51.5.0.50.61000.6LU ⎛⎞⎛⎞⎜⎟⎜⎟∴=⎜⎟⎜⎟⎝⎠⎝⎠§5.4 三角分解法n , 2,i , n, 1,j , 111i111L L ====u a l a u i j j 计算量与Gauss 消去法同.解下三角方程组Ly = b ,即解上三角方程组Ux = y ,即112233100440.5106 4 0.50.6150.6y y y y y y =⎛⎞⎛⎞⎛⎞⎧⎪⎜⎟⎜⎟⎜⎟=⇒=⎨⎜⎟⎜⎟⎜⎟⎪=⎝⎠⎝⎠⎝⎠⎩132231211410 2.51.54 1000.60.61x x x x x x =⎛⎞⎛⎞⎛⎞⎧⎪⎜⎟⎜⎟⎜⎟=⇒=⎨⎜⎟⎜⎟⎜⎟⎪=⎝⎠⎝⎠⎝⎠⎩§5.5 对称正定矩阵的平方根法计算量与Gauss 消去法同.Cholesky分解, A L T A LL L =定理5.5.1:设是对称正定矩阵,则存在唯一的非奇异下三角阵使得且的对角元素皆为正。

《数值分析》练习题及答案解析

《数值分析》练习题及答案解析

《数值分析》练习题及答案解析第一章 绪论主要考查点:有效数字,相对误差、绝对误差定义及关系;误差分类;误差控制的基本原则;。

1. 3.142和3.141分别作为π的近似数具有( )和( )位有效数字.A .4和3B .3和2C .3和4D .4和4 答案:A2. 设 2.3149541...x *=,取5位有效数字,则所得的近似值x=___________ .答案:2.31503.若近似数2*103400.0-⨯=x 的绝对误差限为5105.0-⨯,那么近似数有几位有效数字 解:2*103400.0-⨯=x ,325*10211021---⨯=⨯≤-x x 故具有3位有效数字。

4 . 14159.3=π具有4位有效数字的近似值是多少?解:10314159.0⨯= π,欲使其近似值*π具有4位有效数字,必需!41*1021-⨯≤-ππ,3*310211021--⨯+≤≤⨯-πππ,即14209.314109.3*≤≤π即取( , )之间的任意数,都具有4位有效数字。

第二章 非线性方程求根 主要考查点:二分法N 步后根所在的区间,及给定精度下二分的次数计算;非线性方程一般迭代格式的构造,(局部)收敛性的判断,迭代次数计算; 牛顿迭代格式构造;求收敛阶;1.用二分法求方程012=--x x 的正根,要求误差小于0.05。

(二分法)解:1)(2--=x x x f ,01)0(<-=f ,01)2(>=f ,)(x f 在[0,2]连续,故[0,2]为函数的有根区间。

"(1)计算01)1(<-=f ,故有根区间为[1,2]。

(2)计算041123)23()23(2<-=--=f ,故有根区间为]2,23[。

(3)计算0165147)47()47(2>=--=f ,故有根区间为]47,23[。

(4)计算06411813)813()813(2>=--=f ,故有根区间为]813,23[。

郑州大学数值分析重点考察内容及各章习题

郑州大学数值分析重点考察内容及各章习题

《数值分析》重点考察内容及各章作业答案学院:学号:姓名:重点考察内容基本概念(收敛阶,收敛条件,收敛区域等), 简单欧拉法。

第一章基础掌握:误差的种类,截断误差,舍入误差的来源,有效数字的判断。

了解:误差限,算法及要注意的问题。

第二章插值掌握:Hermite插值,牛顿插值,差商计算,插值误差估计。

了解:Lagrange插值第三章数据拟合掌握:给出几个点求线性拟合曲线。

了解:最小二乘原理第四章数值积分微分掌握:梯形公式,Simpson公式,代数精度,Gauss积分,带权Gauss积分公式推导,复化梯形公式推导及算法。

了解:数值微分,积分余项第五章直接法掌握:LU分解求线性方程组,运算量了解:Gauss消去法,LDL,追赶法第六章迭代法掌握:Jacobi,Gauss-Seidel迭代格式构造,敛散性分析,向量、矩阵的范数、谱半径了解:SOR迭代第七章Nolinear迭代法掌握:牛顿迭代格式构造,简单迭代法构造、敛散性分析,收敛阶。

了解:二分法,弦截法第八章ODE解法掌握:Euler公式构造、收敛阶。

了解:梯形Euler公式、收敛阶,改进Euler公式题目类型:填空,计算,证明综合题第一章 误差1. 科学计算中的误差来源有4个,分别是________,________,________,________。

2. 用Taylor 展开近似计算函数000()()'()()f x f x f x x x ≈+-,这里产生是什么误差?3. 0.7499作34的近似值,是______位有效数字,65.380是舍入得到的近似值,有____几位有效数字,相对误差限为_______. 0.0032581是四舍五入得到的近似值,有_______位有效数字.4. 改变下列表达式,使计算结果比较精确: (1) 11,||1121xx x x --++ (2)||1x(3)1cos ,0,|| 1.xx x x-≠ (4) sin sin ,αβαβ-≈5.采用下列各式计算61)时,哪个计算效果最好?并说明理由。

数值分析复习大纲及课后答案

数值分析复习大纲及课后答案

数值分析复习大纲编者:向穗华时间:2010.5教材:《数值分析(第5版)》. 李庆扬,王能超等编著. 2008年12月第5版.第1章 数值分析与科学计算引论1.1 知识要点总结1. x :准确值2. *x :近似值3. *e :绝对误差 x x e -=**4. *ε:误差限 **ε≤e5. *r e :相对误差 ***xe e r = 6. *r ε:相对误差 ***x r εε=7. *x 具有n 位有效数字,则:)101010(10)1(23121*----⨯++⨯+⨯+±=n n m a a a a x1*1021+-⨯≤-n m x x 1*1021+-⨯=n m ε )1(1*1021--⨯≤n r a ε 8. 误差))(()()(***x x x f x f x f -'≈-误差限 )()())((***'=x x f x f εε9. 误差10≤E E n ,则数值稳定10. 计算函数值问题的条件数10)()(***≥'=x f x f x C p ,则问题是变态的。

11. 避免误差危害,防止有效数字损失,通常要避免两相近数相减和用绝对值很小的数做除数,还要注意运算次序和减少运算次数。

12. 秦九韶n n n n a x a x a x a x p ++++=--1110)( ,求)(*x p 和)(*x p ' 由⎩⎨⎧+==-ii i a x b b a b *100⇒n b x p =)(* 由⎩⎨⎧+==-ii i b x c c b c *100⇒1*)(-='n c x p 1.2 课后习题参考答案1.设0x >,x 的相对误差为δ,求ln x 的误差。

解:令)ln()(x x f =x 的相对误差为δ ∴δ=-=***xx x e r )(x f 的误差为δ=-=-'≈-*****))(()()(x x x x x x f x f x f 所以,ln x 的误差为δ2.设x 的相对误差为2%,求nx 的相对误差。

郑州大学数值分析考试知识点整理

郑州大学数值分析考试知识点整理

郑州⼤学数值分析考试知识点整理题⽬类型:填空(15个共30分),计算,证明综合题第⼀章基础重点考察基本概念(填空题)掌握:误差的种类,截断误差,舍⼊误差的来源,有效数字的判断和定义。

了解:误差限,算法及要注意的问题。

习题:1.误差分类:模型误差 :从实际问题中抽象出数学模型时产⽣的误差。

观测误差 :通过测量得到模型中参数的值导致输⼊数据的误差。

舍⼊误差 :由于计算机字长有限⽽在数值运算的每⼀步所产⽣的误差。

截断误差/⽅法误差 :近似求解时产⽣的误差。

Taylor 展开式近似表达函数产⽣的误差是截断误差/⽅法误差。

2.有效数字定义:设实数x 的近似值为m ka a a x 21.010~?±=其中每个i a 是0,1,2,3....9中的⼀个数字,且01≠a 如果x ~的绝对误差满⾜n k x x x e -?≤-=1021|~||)~(|就称⽤x ~近似x 时具有n 位有效数字,或者说x ~精确到了n 位。

例题:0.2499作为1/4的近似值,有⼏位有效数字?00.24990.249910,0m =?=即,031|0.2499|0.00010.5100.510,34m n n ---=例题:设.......7320508.13==x x 1=1.73, x 2=1.7321, x 3=1.7320是其近似值,问它们分别有⼏位有效数字? 解:31111021 0021.0|3|-?≤≈-=x e 3位 5122102100005.0|3|-?=≤-=x e 5位41331021000051.0|3|-?≤≈-=x e 4位1. 衡量算法优劣的指标有_时间复杂度_,__空间复杂度_.空间复杂度是指:算法需占⽤存储空间的量度。

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

2.两个n 阶矩阵相乘的乘法次数是 3n , 则称两个n 阶矩阵相乘这⼀问题的时间复杂度为3()O n第⼆章插值掌握:差商计算(填空题),插值误差估计(填空题),Hermite 插值(计算题),⽜顿插值。

数值分析计算方法试题集及答案

数值分析计算方法试题集及答案

数值分析复习试题第一章 绪论 一. 填空题 1.*x为精确值x 的近似值;()**x f y =为一元函数()x f y =1的近似值;()**,*y x f y =为二元函数()y x f y ,2=的近似值,请写出下面的公式:**e x x =-:***r x xe x -=()()()*'1**y f x x εε≈⋅ ()()()()'***1**r r x f x y x f x εε≈⋅()()()()()**,**,*2**f x y f x y y x y x yεεε∂∂≈⋅+⋅∂∂()()()()()****,***,**222r f x y e x f x y e y y x y y y ε∂∂≈⋅+⋅∂∂ 2、 计算方法实际计算时,对数据只能取有限位表示,这时所产生的误差叫 舍入误差 。

3、 分别用2.718281,2.718282作数e 的近似值,则其有效数字分别有6 位和 71.73≈(三位有效数字)-211.73 10 2≤⨯。

4、 设121.216, 3.654x x ==均具有3位有效数字,则12x x 的相对误差限为 0.0055 。

5、 设121.216, 3.654x x ==均具有3位有效数字,则12x x +的误差限为 0.01 。

6、 已知近似值 2.4560A x =是由真值T x 经四舍五入得到,则相对误差限为 0.0000204 .7、 递推公式,⎧⎪⎨⎪⎩0n n-1y =y =10y -1,n =1,2,如果取0 1.41y =≈作计算,则计算到10y 时,误差为8110 2⨯;这个计算公式数值稳定不稳定 不稳定 . 8、 精确值 14159265.3*=π,则近似值141.3*1=π和1415.3*2=π分别有 3 位和 4 位有效数字。

9、 若*2.71828x e x =≈=,则x 有 6 位有效数字,其绝对误差限为1/2*10-5 。

数值分析复习---第二章 代数插值

数值分析复习---第二章 代数插值
而满足插值条件(2.5)的n次代数插值多项式Pn(x)也是存 在唯一的.
3/29
郑州大学研究生2012-2013学年课程 数值分析 Numerical Analysis
§2.3 拉格朗日插值方法
n次插值
由抛物插值中构造性方法启发,解决一般的n次代数插值问题.
n次插值基函数
分别构造x0 , x1, …, xn 上的 n 次插值基函数 l0(x), l1(x), …, ln(x),满足
三、典型题目解析 题NAc2-03:
提示: 掌握拉格朗日插值的余项,并会进行不等式的放缩。
24/29
郑州大学研究生2012-2013学年课程 数值分析 Numerical Analysis
三、典型题目解析 题NAc2-04:
提示: 考察差商的性质和牛顿插值多项式的余项。
25/29
郑州大学研究生2012-2013学年课程 数值分析 Numerical Analysis
三、典型题目解析 题NAc2-05:
12/29
郑州大学研究生2012-2013学年课程 数值分析 Numerical Analysis
§2.4 牛顿插值
n次插值
n次Newton插值公式的插值误差为
Rn (x) = (x − x0 )(x − x1)L(x − xn ) f [x, x0 , x1,L, xn ]
n
∏ = f [x, x0 , x1,L, xn ] (x − xi ). i=0
b).
(2.19)
10/29
郑州大学研究生2012-2013学年课程 数值分析 Numerical Analysis
§2.4 牛顿插值
N阶差商表
i xi f (xi) 一阶差商 二阶差商

数值分析 第一章 基础知识.

数值分析  第一章 基础知识.

参考资料
1.关治,陈景良. 数值计算方法. 北京:清华大学出版社,1990. 2.周铁,徐树方等. 计算方法. 北京:清华大学出版社,2006. 3.徐翠微,孙绳武. 计算方法引论. 北京:高等教育出版社,2005. 4.John H.Mathews, Kurtis D.Fink. 数值方法(MATLAB版). 北京:电子 工业出版社,2005. 5.徐士良. 数值分析与算法. 北京:机械工业出版社,2007. 6.葛哲学.精通Matlab.北京:电子工业出版社,2008. 7.任玉杰.数值分析及其MATLAB实现.北京:高等教育出版社,2007.
预备知识
微积分和常微分方程; 线性代数; 数值计算程序设计 (C/Matlab和Mathematica)
5/76
郑州大学2014-2015学年硕士研究生课程 数值分析 Numerical Analysis
参考教材
教材
李庆扬,王能超,易大义.数值分析(第五版).北京:清华大学出版社,2008 李清善,宋士仓. 数值方法. 郑州:郑州大学出版社,2007.
周次 12 13 14 15 16 17 18 19 课次 周二 周二 周二 周二 周二 周二 周二 周二 课时 3课时 3课时 3课时 3课时 3课时 3课时 3课时 3课时
1.8.1 教学内容时间安排
教学内容 备注
第六章 解线性代数方程组的迭代法[1] 第六章 解线性代数方程组的迭代法[2] 第七章 非线性方程的数值解法[1] 第七章 非线性方程的数值解法[2] 第八章 常微分方程初值问题数值解法[1] 第八章 常微分方程初值问题数值解法[2] 习题课 总复习
密码:zzumoe2014(请不要修改)
下载电子教案、算法实现代码及课后习题。

郑州大学研究生课程数值分析 第三章习题

郑州大学研究生课程数值分析 第三章习题
给出数据习题xyxy0100022090750329502525645205088260500753706102514392100428362000331334用一次二次三次多项式及最小二乘原理拟合这些数据并写出正规方程组
郑州大学研究生课程( 2012-2013 学年上学期 )
数值分析 Numerical Analysis
-0.25 1.4392 1.00 4.2836
2.0003
用一次、二次、三次多项式及最小二乘原理拟合这些数据,并写出正规方程组. 3. 用最小二乘原理求一个形如 y = a + bx 的经验公式,使与下列数据相拟合.
2
x y
19
19.0
25
32.3
31
49.0
38
73.3
44 97.8
4. 用最小二乘原理求一个形如 y = Ae 的经验公式,使与下列数据相拟合.
第三章
⎧ x1 − x2 = 5 ⎪ 1. 求 ⎨ x1 + x2 = −4 的最小二乘解. ⎪2 x − 2 x = 10 2 ⎩ 1
2. 给出数据


x y x y
0
-1.00 -0.2209
-0.75 0.3295 0.25 2.5645 0.50 3.1334
-0.5 0.8826 0.75 3.7061
Bx
x y5. 给出实验数据Fra bibliotek1 60
2 30
2 20
4 15
x y
6
4.6
8
4.8
10
4.6
12
4.9
14
5.0
16
5.4
18
5.1
20

数值分析课后习题及答案

数值分析课后习题及答案

数值分析课后习题及答案第一章绪论(12)第二章插值法(40-42)2、当时,,求的二次插值多项式。

[解]。

3、给出的数值表用线性插值及二次插值计算的近似值。

X 0.4 0.5 0.6 0.7 0.8 -0.916291 -0.693147 -0.510826 -0.357765 -0.223144 [解]若取,,则,,则,从而。

若取,,,则,,,则,从而补充题:1、令,,写出的一次插值多项式,并估计插值余项。

[解]由,可知,,余项为,故。

2、设,试利用拉格朗日插值余项定理写出以为插值节点的三次插值多项式。

[解]由插值余项定理,有,从而。

5、给定数据表:,1 2 4 6 7 4 1 0 1 1 求4次牛顿插值多项式,并写出插值余项。

[解]一阶差商二阶差商三阶差商四阶差商 1 42 1 -34 0 6 17 1 0 由差商表可得4次牛顿插值多项式为:,插值余项为。

第三章函数逼近与计算(80-82)26、用最小二乘法求一个形如的经验公式,使它与下列数据相拟合,并求均方误差。

19 25 31 38 44 19.0 32.3 49.0 73.3 97.8[解]由。

又,,,故法方程为,解得。

均方误差为。

27、观测物体的直线运动,得出以下数据:时间t(秒)0 0.9 1.9 3.0 3.9 5.0 距离s(米)0 10 30 5080 110 [解]设直线运动为二次多项式,则由。

,。

又,,,故法方程为,解得。

故直线运动为。

补充题:1、现测得通过某电阻R的电流I及其两端的电压U如下表:I ……U ……试用最小二乘原理确定电阻R的大小。

[解]电流、电阻与电压之间满足如下关系:。

应用最小二乘原理,求R使得达到最小。

对求导得到:。

令,得到电阻R为。

2、对于某个长度测量了n次,得到n个近似值,通常取平均值作为所求长度,请说明理由。

[解]令,求x使得达到最小。

对求导得到:,令,得到,这说明取平均值在最小二乘意义下误差达到最小。

郑州大学研究生课程数值分析复习---第八章 常微分方程数值解法

郑州大学研究生课程数值分析复习---第八章 常微分方程数值解法

郑州大学研究生课程(2012-2013学年第一学期)数值分析Numerical Analysis习题课第八章常微分方程数值解法待求解的问题:一阶常微分方程的初值问题/* Initial-Value Problem */:⎪⎩⎪⎨⎧=∈=0)(],[),(y a y b a x y x f dx dy 解的存在唯一性(“常微分方程”理论):只要f (x , y ) 在[a , b ] ×R 1 上连续,且关于y 满足Lipschitz 条件,即存在与x , y 无关的常数L 使对任意定义在[a , b ] 上的y 1(x ) 和y 2(x ) 都成立,则上述IVP 存在唯一解。

1212|(,)(,)|||f x y f x y L y y −≤−一、要点回顾§8.2 欧拉(Euler)法通常取(常数),则Euler 法的计算格式h h x x i i i ==−+1⎩⎨⎧=+=+)(),(001x y y y x hf y y i i i i i =0,1,…,n ( 8.2 )§8.2 欧拉(Euler)法(1) 用差商近似导数))(,()()()()(1n n n n n n x y x hf x y x y h x y x y +=′+≈+⎩⎨⎧=+=+)(),(01a y y y x hf y y n n n n 差分方程初值问题向前Euler 方法h x y x y x y n n n )()()(1−≈′+))(,()()(1n n n n x y x f hx y x y ≈−+))(,()(n n n x y x f x y =′§8.2 欧拉(Euler)法若用向后差商近似导数,即))(,()()(111++++≈n n n n x y x hf x y x y ⎩⎨⎧=+=+++)(),(0111a y y y x hf y y n n n n hx y x y x y n n n )()()(11−≈′++向后Euler 方法))(,()()(111+++≈−n n n n x y x f h x y x y ))(,()(111+++=′n n n x y x f x y§8.2 欧拉(Euler)法(2)用数值积分方法∫+=−+1))(,()()(1n n x x n n dxx y x f x y x y ∫∫++=′11))(,()(n n n n x x x x dx x y x f dx x y 1111(,)(,()), (,)(,()),()()n n n n n n n n f x y f x y x f x y f x y x y x y y x y ++++≈≈≈≈分别用左矩形和右矩形公式,即代替上式右端的积分,并注意 ,分别得到1111(,)(,)n n n n n n n n y y h f x y y y h f x y ++++=+=+,。

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

《数值分析》重点考察内容及各章作业答案学院:学号:姓名:重点考察内容基本概念(收敛阶,收敛条件,收敛区域等), 简单欧拉法。

第一章基础掌握:误差的种类,截断误差,舍入误差的来源,有效数字的判断。

了解:误差限,算法及要注意的问题。

第二章插值掌握:Hermite插值,牛顿插值,差商计算,插值误差估计。

了解:Lagrange插值第三章数据拟合掌握:给出几个点求线性拟合曲线。

了解:最小二乘原理第四章数值积分微分掌握:梯形公式,Simpson公式,代数精度,Gauss积分,带权Gauss积分公式推导,复化梯形公式推导及算法。

了解:数值微分,积分余项第五章直接法掌握:LU分解求线性方程组,运算量了解:Gauss消去法,LDL,追赶法第六章迭代法掌握:Jacobi,Gauss-Seidel迭代格式构造,敛散性分析,向量、矩阵的范数、谱半径了解:SOR迭代第七章Nolinear迭代法掌握:牛顿迭代格式构造,简单迭代法构造、敛散性分析,收敛阶。

了解:二分法,弦截法第八章ODE解法掌握:Euler公式构造、收敛阶。

了解:梯形Euler公式、收敛阶,改进Euler公式题目类型:填空,计算,证明综合题第一章 误差1. 科学计算中的误差来源有4个,分别是________,________,________,________。

2. 用Taylor 展开近似计算函数000()()'()()f x f x f x x x ≈+-,这里产生是什么误差?3. 0.7499作34的近似值,是______位有效数字,65.380是舍入得到的近似值,有____几位有效数字,相对误差限为_______. 0.0032581是四舍五入得到的近似值,有_______位有效数字.4. 改变下列表达式,使计算结果比较精确: (1) 11,||1121xx x x --++ (2)||1x(3)1cos ,0,|| 1.xx x x-≠ (4) sin sin ,αβαβ-≈5.采用下列各式计算61)时,哪个计算效果最好?并说明理由。

(1)(2)99-(3)6(3- (46. 已知近似数*x 有4位有效数字,求其相对误差限。

上机实验题:1、利用Taylor 展开公式计算 0!kxk x e k ∞==∑,编一段小程序,上机用单精度计算xe 的函数值. 分别取 x =1,5,10,20,-1,-5,-10,-15,-20,观察所得结果是否合理,如不合理请分析原因并给出解决方法. 2、已知定积分1,0,1,2,,206nn x I dx n x ==+⎰,有如下的递推关系1111100(6)61666n n n n n x x x x I dx dx I x x n ---+-===++-⎰⎰可建立两种等价的计算公式 (1) 1016,0.154n n I I I n -=-=取; (2) 12011),0.6n n I nI I n-=-=(取来计算123419,,,,,I I I I I ,编程比较哪种计算的数值结果好,并给出理论分析。

第二章 插值法1. 已知(0)2,(1)1f f ==-,那么差商[1,0]f =_________.2. n 阶差商与导数的关系是01[,,,]n f x x x =__________________.3. 由导数和差商的关系知,[,]i i f x x =__________________。

4. 已知函数()f x 在3,1,4x =的值分别是4,6,9,试构造Lagrange 插值多项式。

5.取节点0120,1,2x x x ===, 对应的函数值和导数值分别为0()1,f x = 11()2,'()2f x f x ==,试建立不超过二次的插值多项式。

(如果将最后一个条件改为2'()2f x =,插值多项式如何计算?)6.已知(0)1,(1)2,'(1)3,(2)9f f f f ====,试建立不超过3次的插值多项式,并写出插值余项.7. 设4()[,]f x C a b ∈,求三次多项式3()p x ,使之满足插值条件11()(),0,1,2'()'()i i p x f x i p x f x ==⎧⎨=⎩8. 设1()P x 是过01,x x 的一次插值多项式,2()[,],f x C a b ∈其中[,]a b 是包含01,x x 的任一区间。

试证明:对任一给定的[,]x a b ∈,在(a,b )上总存在一点ξ,使得101()()()()()()2!f R x f x P x x x x x ξ''=-=--。

9.证明关于互异节点0{}n i i x =的Lagrange 插值基函数0{()}n i i l x =满足恒等式01()()()1n l x l x l x +++≡上机习题:1. 绘制4题的Lagrange 的插值函数的图像。

第三章 数据拟合1. 数据拟合与插值的区别是什么?2. 最小二乘原理是使偏差iδ的___________达到最小3. 求过点(2,3),(0,1),(3,5)的线性拟合函数。

4. 用最小二乘法求一形如2y a bx=+的多项式,使与下列数据相拟合第四章线性方程组的直接解法1. 线性方程组的解法大致可分为_____________,________________。

2. 平方根法和LDL T分解法要求系数矩阵A满足______________。

3. 上三角和下三角方程组的解法分别称为___________,____________。

4. 严格对角占优矩阵的定义是什么?5. 试求下面矩阵的杜利特尔分解(1)6234⎡⎤⎢⎥-⎣⎦。

(2)213 457 285⎡⎤⎢⎥⎢⎥⎢⎥-⎣⎦。

6. 用列主元高斯消去法求解方程组1231521 04313 2063xxx-⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦⎣⎦。

7. 用LU分解法解方程组1232111 6161 10272xxx⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥-=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦。

上机实验题:1.编程实现列主元的高斯消去法2.编程实现LU分解法第五章线性方程组的迭代解法1. 向量(3,2,1,7)T x =--,计算1||||x ,2||||x ,||||x ∞.2. A=312010126-⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦,计算1||||A ,2||||A ,||||A ∞. 3. 2003A ⎡⎤=⎢⎥⎣⎦, 分别计算A 的谱半径()A ρ, 条件数cond ()A ∞,1||||A 4. 矩阵A 的范数与谱半径的关系为__________________________。

5. 求解AX =b 的迭代格式(1)()k k x Bx g +=+收敛的充分必要条件____________________。

6. SOR 迭代法收敛的一个必要条件是松驰因子______________。

7. 写出下面方程的Jacobi 迭代格式1231231231027102854x x x x x x x x x --=⎧⎪-+-=⎨⎪--+=⎩ 8. 给定下列方程组,判断对它们构造的Jacobi 迭代公式和Gauss-Seidel 迭代公式是否收敛(1) 131252728x x x x +=⎧⎨+=⎩ (2) 123121355251285x x x x x x x -+=⎧⎪-+=⎨⎪+=⎩9. 对下列方程组建立收敛的简单迭代公式(提示:先调整方程组)123162132624114x x x -⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥-=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦⎣⎦ 10. 给定方程组123122*********x x x -⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦, (1)分别写出Jacobi 迭代公式和Gauss-Seidel 迭代公式。

(2)证明Jacobi 迭代法收敛,而Gauss-Seidel 迭代法发散。

上机实验题:1. 求解方程组:1231231231027102854x x x x x x x x x --=⎧⎪-+-=⎨⎪--+=⎩以(0)(1,1,1)T x =为初值,当(1)()4||||10k k x x +-∞-<时迭代终止。

(1) 编写Jacobi 迭代法程序 (2) 编写Gauss-Seidel 迭代法程序第六章 数值积分与数值微分1.()ba f x dx ⎰的梯形求积公式是________,Simpson 公式是_______,其代数精度分别为_____,____。

2. n 点Gauss 求积公式的代数精度为___________.3. 确定下列求积公式中的待定系统,使得求积公式的代数精度尽量的高,并指明代数精度 (1)101()()(0)()hh f x dx A f h A f A f h -≈-++⎰(2) 11211()[(1)2()3()]3f x dx f f x f x -≈-++⎰(3) 10100()(0)(1)'(0)f x dx A f A f B f ≈++⎰4. 分别用梯形公式、Simpson 公式、Cotes 公式、Gauss 求积公式计算积分11x e dx -⎰,并估计各种方法的误差。

5. 写出11()f x dx -⎰二点和三点的Gauss-Legendre 求积公式.6. 分别用复化梯形公式和复化Simpson 公式计算下列积分.10,(8)4xdx n x =+⎰7. 确定求积公式10101()()(1)2f x dx A f A f ≈+⎰的求积公式,并求其代数精度。

8. 构造如下形式的Gauss 求积公式:00110()()()x dx A f x A f x ≈+⎰.9. 构造如下形式的Gauss 求积公式: 100111()()()f x dx A f x A f x -≈+⎰.上机实验题:1. 编程实现五点Gauss 积分算法。

第七章 非线性方程与非线性方程组的解法1. 求解非线性方程的根,牛顿法的收敛阶是________,割线法的收敛阶是____________.2. 确定下列方程的有根区间 (1) 32720x x -+= (2) 20x e x -+-=3. 试用牛顿法和弦截法建立计算1,(0)c c≠的迭代格式。

4(0)a>的两种收敛的迭代格式。

5(0)a >4位有效数字。

(迭代求解3次即可)6. 2+++的近似值.7. 设初值00x ≠, 计算1,(0)a a≠的迭代格式 1(2),0,1,2,k k k x x ax k +=-=。

试证:(1)此迭代格式二阶收敛.(2)此迭代格式收敛的充分必要条件为0|1|1ax -<. 上机实验题:1. 用割线法求方程32210200x x x ++-=的根,要求61||10k k x x -+-<第八章 常微分方程初值问题的数值解法1. 求解常微分方程'(0)1y x yy =+⎧⎨=⎩的Euler 公式为______________________, 其局部截断误差的阶数为_________,整体截断误差的阶数为__________.(设步长为h) 2. 应用向前欧拉格式求解初值问题'1,01(0)1y x y x y =-+≤≤⎧⎨=⎩取步长h = 0.1,将计算结果与精确解x y x e -=+对照.。

相关文档
最新文档