华南理工大学数值分析试题2014年

华南理工大学数值分析试题2014年
华南理工大学数值分析试题2014年

数值分析学期期末考试试题与答案(A)

期末考试试卷(A 卷) 2007学年第二学期 考试科目: 数值分析 考试时间:120 分钟 学号 姓名 年级专业 一、判断题(每小题2分,共10分) 1. 用计算机求 1000 1000 1 1 n n =∑时,应按照n 从小到大的顺序相加。 ( ) 2. 为了减少误差,进行计算。 ( ) 3. 用数值微分公式中求导数值时,步长越小计算就越精确。 ( ) 4. 采用龙格-库塔法求解常微分方程的初值问题时,公式阶数越高,数值解越精确。( ) 5. 用迭代法解线性方程组时,迭代能否收敛与初始向量的选择、系数矩阵及其演变方式有 关,与常数项无关。 ( ) 二、填空题(每空2分,共36分) 1. 已知数a 的有效数为0.01,则它的绝对误差限为________,相对误差限为_________. 2. 设1010021,5,1301A x -????????=-=-????????-???? 则1A =_____,2x =______,Ax ∞ =_____. 3. 已知5 3 ()245,f x x x x =+-则[1,1,0]f -= ,[3,2,1,1,2,3]f ---= . 4. 为使求积公式 1 1231 ()()(0)33 f x dx A f A f A f -≈- ++? 的代数精度尽量高,应使1A = ,2A = ,3A = ,此时公式具有 次的代数精度。 5. n 阶方阵A 的谱半径()A ρ与它的任意一种范数A 的关系是 . 6. 用迭代法解线性方程组AX B =时,使迭代公式(1) ()(0,1,2,)k k X MX N k +=+=产 生的向量序列{ }() k X 收敛的充分必要条件是 . 7. 使用消元法解线性方程组AX B =时,系数矩阵A 可以分解为下三角矩阵L 和上三角矩

2016华工计算机计算方法(数值分析)考试试卷

考完试了,顺便把记得的题目背下来,应该都齐全了。我印象中也就只有这些题,题目中的数字应该是对的,我也验证过,不过也不一定保证是对的,也有可能我也算错了。还有就是试卷上面的题目可能没有我说的这么短,但是我也不能全把文字背下来,大概意思就是这样吧。每个部分的题目的顺序可能不是这样,但总体就是这四大块。至于每道题目的分值,我记得的就写出来了,有些题目没注意。我题目后面写的结果都是我考试时算出来的,考完了也懒得验证了,可能不一定对,自己把握吧,仅供参考。 华南理工大学2016计算机计算方法(数值分析)考试试卷 一填空题(16分) 1.(6分)X* = 3.14,准确值x = 3.141592,求绝对误差e(x*) =,相对误差e r(x*) =,有效数 位是。 2.(4分)当插值函数的n越大时,会出现龙格现象,为解决这个问题,分段函数不一个 不错的办法,请写出分段线性插值、分段三次Hermite插值和三次样条插值各自的特点。 3.(3分)已知x和y相近,将lgx – lgy变换成可以使其计算结果更准确。 4.(3分)已知2x3– 3x2 +2 = 0,求牛顿迭代法的迭代式子。 解题思路:1. 这里的绝对误差和相对误差是没有加绝对值的,而且要注意是用哪个数减去哪个数得到的值,正负号会不一样;2. 可以从它们函数的连续性方面来说明;3. 只要满足课本所说的那几个要求就可以;这个记得迭代公式就可以直接写,记不住可以自己推导,就是用泰勒展开式来近似求值得到的迭代公式。 我最终的结果是: 1.-0.001592 -0.000507 3 2.分段线性插值保证了插值函数的连续性,但是插值函数的一次导数不一定连续; 分段三次Hermite既保证了插值函数的连续性,也保证了其一次导数的连续性; 三次样条插值保证了插值函数及其一次导数和二次导数的连续性 3.lg(x/y) 4.x k+1 = x k– (2x3– 3x2 +2)/(6x2 -6x) 二计算题(64分) 1.已知f(x) = x3–x -1,用对分法求其在[0 , 2]区间内的根,误差要满小于0.2,需要对分多少 次?请写出最后的根结果。 解题思路:每次求区间的中值并计算其对应的函数值,然后再计算下一个区间中值及函数值,一直到两次区间中值的绝对值小于0.2为止。 我最终算得的对分次数是4,根的结果为11/8. 2. (1)请根据以上数据构造Lagrange三次插值函数; (2)请列出差商表并写出Newton三次插值函数。 解题思路:(1) 直接按照书本的定义把公式列出来就可以了,这个要把公式记住了才行,不然也写不了;(2)差商表就是计算Newton三次插值函数过程中计算到的中间值及结果值,可以先在草稿上按照Newton公式的计算过程把公式写出来,然后把中间用到的值

数值分析课后习题部分参考答案

数值分析课后习题部分参考答案 Chapter 1 (P10) 5.求.2的近似值x *,使其相对误差不超过 0.1%。 解:,2 =1.4…。 、I * * n 设X 有n 位有效数字,则|e(x)|乞0.5 10 10』。 0.5 101 』 从而,|e r (x)匡 。 1 故,若0.5 101』乞0.1%,则满足要求。 解之得,n 丄4。 x =1.414。 (P10) 7.正方形的边长约100cm ,问测量边长时误差应多大,才能保证面积的误差不超 过 1 cm 2。 解:设边长为a ,则a 100 cm 。 设测量边长时的绝对误差为 e ,由误差在数值计算的传播,这时得到的面积的绝对误差有如 下估计::、2 100 e 。按测量要求,|2 100 e|_1 解得,|e|空 0.5 10 °。 Chapter 2 (P47) 5.用三角分解法求下列矩阵的逆矩阵 1 1 -P A = 2 1 0 。 -1 0丿 解:设A 4 '。分别求如下线性方程组: V '0 ' ■0 ' 0 ,A B = 1 ,AY = 0 6

"(1)1 (1)1 (-1) -1 " (2)2 (1) -1 (0)2 .(1)1 (-1)2 (0)-3」 「 1 0 0^ 「1 1 —1 ' 即,L = 2 1 0 ,U = 0 -1 2 。 2 b 1° 0 -3」 经直接三角分解法的回代程,分别求解方程组, 3 1 3 2 (P47)6.分别用平方根法和改进平方根法求解方程组: (1 2 1 -3、 1 5、 2 5 0 -5 1 X 2 2 1 0 14 1 X 3 16 1—3 -5 1 15丿

数值分析试题及答案汇总

数值分析试题 一、 填空题(2 0×2′) 1. ?? ????-=? ?????-=32,1223X A 设x =是精确值x *=的近似值,则x 有 2 位 有效数字。 2. 若f (x )=x 7-x 3+1,则f [20,21,22,23,24,25,26,27]= 1 , f [20,21,22,23,24,25,26,27,28]= 0 。 3. 设,‖A ‖∞=___5 ____,‖X ‖∞=__ 3_____, ‖AX ‖∞≤_15_ __。 4. 非线性方程f (x )=0的迭代函数x =?(x )在有解区间满足 |?’(x )| <1 ,则使用该迭代 函数的迭代解法一定是局部收敛的。 5. 区间[a ,b ]上的三次样条插值函数S (x )在[a ,b ]上具有直到 2 阶的连续导数。 6. 当插值节点为等距分布时,若所求节点靠近首节点,应该选用等距节点下牛顿差商 公式的 前插公式 ,若所求节点靠近尾节点,应该选用等距节点下牛顿差商公式的 后插公式 ;如果要估计结果的舍入误差,应该选用插值公式中的 拉格朗日插值公式 。 7. 拉格朗日插值公式中f (x i )的系数a i (x )的特点是:=∑=n i i x a 0)( 1 ;所以当 系数a i (x )满足 a i (x )>1 ,计算时不会放大f (x i )的误差。 8. 要使 20的近似值的相对误差小于%,至少要取 4 位有效数字。 9. 对任意初始向量X (0)及任意向量g ,线性方程组的迭代公式x (k +1)=Bx (k )+g (k =0,1,…)收 敛于方程组的精确解x *的充分必要条件是 ?(B)<1 。 10. 由下列数据所确定的插值多项式的次数最高是 5 。 11. 牛顿下山法的下山条件为 |f(xn+1)|<|f(xn)| 。 12. 线性方程组的松弛迭代法是通过逐渐减少残差r i (i =0,1,…,n )来实现的,其中的残差 r i = (b i -a i1x 1-a i2x 2-…-a in x n )/a ii ,(i =0,1,…,n )。 13. 在非线性方程f (x )=0使用各种切线法迭代求解时,若在迭代区间存在唯一解,且f (x )

11:数值分析试题2009~2010

中国石油大学(北京)2009--2010学年第一学期 研究生期末考试试题A (闭卷考试) 课程名称:数值分析 注:计算题取小数点后四位 一、填空题(共30分,每空3分) 1、 已知x =0.004532是由准确数a 经四舍五入得到的近似值,则x 的绝对误差 界为_______________。 2、数值微分公式()() '()i i i f x h f x f x h +-≈ 的截断误差为 。 3、已知向量T x =,求Householder 变换阵H ,使(2,0)T Hx =-。 H = 。 4、利用三点高斯求积公式 1 1 ()0.5556(0.7746)0.8889(0)0.5556(0.7746) f x d x f f f -≈-++? 导出求积分 4 0()f x dx ?的三点高斯求积公式 。 5、4 2 ()523,[0.1,0.2,0.3,0.4,0.5]_____.f x x x f =+-=若则 6、以n + 1个互异节点x k ( k =0,1,…,n ),(n >1)为插值节点的 Lagrange 插值基函数为l k (x)( k =0,1,…,n ),则 (0)(1)__________.n k k k l x =+=∑ 7、已知3()P x 是用极小化插值法得到的cos x 在[0,4]上的三次插值多项式,则3()P x 的 截断误差上界为3()cos ()R x x P x =-≤_________. 8、已知向量(3,2,5)T x =-,求Gauss 变换阵L ,使(3,0,0)T Lx =。L =_________. 9、设3 2 ()(7)f x x =-, 给出求方程()0f x =根的二阶收敛的迭代格式_________。

华南理工大学数值分析试题-14年下-A

《数值分析》A 卷 第 1 页 共 2 页 华南理工大学研究生课程考试 《数值分析》试卷A (2015年1月9日) 1. 考前请将密封线内各项信息填写清楚; 所有答案请按要求填写在本试卷上; 课程代码:S0003004; 4. 考试形式:闭卷; 5. 考生类别:硕士研究生; 本试卷共八大题,满分100分,考试时间为150分钟。 一.(12分)解答下列问题 1.欲计算下式: ()13(1)2(1)(2)7(1)(2)(3)6(1)(2)(3)(4),P x x x x x x x x x x x =+-+------+---- 2.设有递推公式 0161,1,2,n n y y y n -?=??=-=?? *001.732y y = 作实际计算,问计算到10y 时误差为初始误差*00y y -的多少 这个计算过程数值稳定吗 ? . (14分)解答下列问题 1. 若2()63f x x =+,则[1,2,3]f 和[1,2,34]f ,的值分别是多少? 2. 1012 . (10分) 设f 在互易节点i x 上的值()()0,1,....i i f f x i n ==。试证明:f 在节点i x 上n 次最小二乘拟合多项式()n p x 与f 在节点i x 上的n 次Lagrange 插值多项式()n L x 一致,()()=n n p x L x 。 . (12分) 按代数精度的定义,构造下列形式的求积公式(即确定参数,A B ,α): ()()()11f x dx Af Bf αα-≈-+? Gauss 型求积公式。

《数值分析》A 卷 第 2 页 共 2 页 五. (14分) 已知线性代数方程组Ax=b 为: (1) 用顺序高斯消去法求解方程组Ax=b ; (2) 先由(1)的消元过程直接写出A 的LU 分解,再利用该LU 分解求解方程组Ax=b 。 六. (12分) 对方程组323,,121Ax b A b ????===????-???? ,拟用迭代法 (1)()()(),0,1,k k k x x Ax b k α+=+-= 求解,试确定实数α的取值范围,使得该迭代公式收敛。 七. (14分) 欲求方程 ln 2 (1)x x x -=> 的根,试 (1)证明 [3, 4] 为方程的一个有根区间; (2)在区间 [3, 4] 上构造一个收敛的不动点迭代公式; (3)求所构造迭代公式的收敛阶。 八. (12分) 对初值问题 ()()00 y f xy y x y '=???=?? (1)试利用Taylor 展开公式推导下列数值求解公式: ()()()212 n n n n n n n n n n h y y hf x y f x y y x f x y +'=+++???? (2)指出上述公式是几阶公式。 ??????? ?????????=????????????????????????????????-----n n n n n n n n b b b b x x x x d u u u v d v d v d 12112112111221100 0000 . 0)/(,0,1 1,,,≠-≠∑-=n i i i i n i i i i i d v u d d b v u d 已知其中

数值分析课后习题部分参考答案

数值分析课后习题部分参考答案 Chapter 1 (P10)5. 求2的近似值* x ,使其相对误差不超过%1.0。 解: 4.12=。 设* x 有n 位有效数字,则n x e -??≤10 105.0|)(|* 。 从而,1 105.0|)(|1* n r x e -?≤。 故,若%1.010 5.01≤?-n ,则满足要求。 解之得,4≥n 。414.1* =x 。 (P10)7. 正方形的边长约cm 100,问测量边长时误差应多大,才能保证面积的误差不超过12 cm 。 解:设边长为a ,则cm a 100≈。 设测量边长时的绝对误差为e ,由误差在数值计算的传播,这时得到的面积的绝对误差有如下估计:e ??≈1002。按测量要求,1|1002|≤??e 解得,2 105.0||-?≤e 。 Chapter 2 (P47)5. 用三角分解法求下列矩阵的逆矩阵: ???? ? ??--=011012111A 。 解:设()γβα=-1 A 。分别求如下线性方程组: ????? ??=001αA ,????? ??=010βA ,???? ? ??=100γA 。 先求A 的LU 分解(利用分解的紧凑格式),

???? ? ??-----3)0(2)1(1)1(2)0(1)1(2)2(1)1(1)1(1 )1(。 即,? ??? ? ??=121012001L ,??? ?? ??---=300210111U 。 经直接三角分解法的回代程,分别求解方程组, ????? ??=001Ly 和y U =α,得,???? ? ??-=100α; ???? ? ??=010Ly 和y U =β,得,???????? ??=3231 31β; ???? ? ??=100Ly 和y U =γ,得,;???????? ??--=3132 31γ。 所以,??????? ? ? ?---=-313 2132310 313101A 。 (P47)6. 分别用平方根法和改进平方根法求解方程组: ????? ? ? ??=??????? ????????? ??----816 2115153114015052 31214321x x x x 解: 平方根法: 先求系数矩阵A 的Cholesky 分解(利用分解的紧凑格式),

数值分析试卷及其答案

1、(本题5分)试确定7 22 作为π的近似值具有几位有效数字,并确定其相对误差限。 解 因为 7 22 =3.142857…=1103142857 .0-? π=3.141592… 所以 312102 11021005.0001264.0722--?=?=<=- π (2分) 这里,3,21,0=-=+-=n n m m 由有效数字的定义可知7 22 作为π的近似值具有3位有效数字。 (1分) 而相对误差限 3102 1 0005.00004138.0001264.07 22-?= <≈= -= π π πε r (2分) 2、(本题6分)用改进平方根法解方程组:??? ?? ??=????? ??????? ??--654131*********x x x ; 解 设???? ? ??????? ? ?????? ??===????? ??--11111 1 131321112323121 32 132 31 21 l l l d d d l l l LDL A T 由矩阵乘法得: 5 7,21,215 27 ,25,2323121321- ==-== -==l l l d d d (3分) 由y D x L b Ly T 1 ,-==解得 T T x y )9 23 ,97,910(,)563, 7,4(== (3分) 3、(本题6分)给定线性方程组???????=++-=+-+=-+-=-+17 7222382311387 510432143213 21431x x x x x x x x x x x x x x 1)写出Jacoib 迭代格式和Gauss-Seidel 迭代格式; 2)考查Jacoib 迭代格式和Gauss-Seidel 迭代格式的敛散性; 解 1)Jacoib 迭代格式为

数值分析试题

华南理工大学研究生课程考试《数值分析》试卷 A (2015年1月9日)注意事项: 1. 考前请将密封线内各项信息填写清楚; 2. 所有答案请按要求填写在本试卷上; 3. 课程代码:S0003004; 4. 考试形式:闭卷; 5. 考生类别:硕士研究生; 6. 本试卷共八大题,满分100分,考试时间为150分钟。一.(12分)解答下列问题1.欲计算下式:()13(1)2(1)(2)7(1)(2)(3)6(1)(2)(3)(4),P x x x x x x x x x x x 试给出乘法次数尽可能少的计算形式。2.设有递推公式01361,1,2,n n y y y n 如果取*003 1.732y y 作实际计算,问计算到10y 时误差为初始误差*00y y 的多少倍?这个计算过程数值稳定吗?二. (14分)解答下列问题_____________________姓名学 号学 院专 业任 课教师(密封线内不答题)…… … … … …………密………………………………………………封………………………………………线…………………………………………………

1. 若2() 63f x x +,则[1,2,3]f 和[1,2,34]f ,的值分别是多少?2. 已知100101211114412===,,,试利用二次插值方法求115的近似值,并估计误差。 三. (10分) 设f 在互易节点 i x 上的值0,1,....i i f f x i n 。试证明:f 在节点i x 上

的n 次最小二乘拟合多项式n p x 与f 在节点i x 上的n 次Lagrange 插值多项式n L x 一致,即=n n p x L x 。 四. (12分) 按代数精度的定义,构造下列形式的求积公式(即确定参数,A B ,):

数值分析第四版习题及答案

第四版 数值分析习题 第一章绪论 1.设x>0,x得相对误差为δ,求得误差、 2.设x得相对误差为2%,求得相对误差、 3.下列各数都就是经过四舍五入得到得近似数,即误差限不超过最后一位得半个单位,试指 出它们就是几位有效数字: 4.利用公式(3、3)求下列各近似值得误差限: 其中均为第3题所给得数、 5.计算球体积要使相对误差限为1%,问度量半径R时允许得相对误差限就是多少? 6.设按递推公式 ( n=1,2,…) 计算到、若取≈27、982(五位有效数字),试问计算将有多大误差? 7.求方程得两个根,使它至少具有四位有效数字(≈27、982)、 8.当N充分大时,怎样求? 9.正方形得边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝? 10.设假定g就是准确得,而对t得测量有±0、1秒得误差,证明当t增加时S得绝对误差增 加,而相对误差却减小、 11.序列满足递推关系(n=1,2,…),若(三位有效数字),计算到时误差有多大?这个计算过程 稳定吗? 12.计算,取,利用下列等式计算,哪一个得到得结果最好? 13.,求f(30)得值、若开平方用六位函数表,问求对数时误差有多大?若改用另一等价公式 计算,求对数时误差有多大? 14.试用消元法解方程组假定只用三位数计算,问结果就是否可靠? 15.已知三角形面积其中c为弧度,,且测量a ,b ,c得误差分别为证明面积得误差满足 第二章插值法 1.根据(2、2)定义得范德蒙行列式,令 证明就是n次多项式,它得根就是,且 、 2.当x= 1 , -1 , 2 时, f(x)= 0 , -3 , 4 ,求f(x)得二次插值多项式、 3. 4., 研究用线性插值求cos x 近似值时得总误差界、

数值分析整理版试题及答案

数值分析整理版试题及答案

例1、 已知函数表 x -1 1 2 ()f x -3 0 4 求()f x 的Lagrange 二次插值多项式和Newton 二次插值多项式。 解: (1)k x -1 1 2 k y -3 0 4 插值基函数分别为 ()()()()()()()()()() 1200102121()1211126 x x x x x x l x x x x x x x ----= ==-------- ()()()()()()()() ()()021******* ()1211122x x x x x x l x x x x x x x --+-= ==-+---+- ()()()()()()()()()()0122021111 ()1121213 x x x x x x l x x x x x x x --+-= ==-+--+- 故所求二次拉格朗日插值多项式为 () ()()()()()()()()()()2 20 2()11131201241162314 121123537623k k k L x y l x x x x x x x x x x x x x ==?? =-? --+?-+-+?+-????=---++-=+-∑ (2)一阶均差、二阶均差分别为

[]()()[]()()[][][]010********* 011201202303 ,11204 ,412 3 4,,5 2,,126 f x f x f x x x x f x f x f x x x x f x x f x x f x x x x x ---===-----= = =----=== --- k x ()k f x 一阶 二阶 -1 -3 1 0 3/ 2 2 4 4 5/6 故所求Newton 二次插值多项式为 ()()[]()[]()() ()()()20010012012,,,35 311126537623P x f x f x x x x f x x x x x x x x x x x x =+-+--=-+ +++-=+- 例2、 设2 ()32f x x x =++,[0,1]x ∈,试求()f x 在[0, 1]上关于()1x ρ=,{} span 1,x Φ=的最佳平方逼近多项式。 解: 若{}span 1,x Φ=,则0()1x ?=,1()x x ?=,且()1x ρ=,这样,有

数值分析试题

文档来源为:从网络收集整理.word 版本可编辑.欢迎下载支持. 1文档来源为:从网络收集整理.word 版本可编辑. 华南理工大学研究生课程考试 《数值分析》试卷A (2015年1月9日) 1. 考前请将密封线内各项信息填写清楚; 所有答案请按要求填写在本试卷上; 课程代码:S0003004; 4. 考试形式:闭卷; 5. 考生类别:硕士研究生; 本试卷共八大题,满分100分,考试时间为150分钟。 一.(12分)解答下列问题 1.欲计算下式: ()13(1)2(1)(2)7(1)(2)(3)6(1)(2)(3)(4),P x x x x x x x x x x x =+-+------+---- 2.设有递推公式 0161,1,2,n n y y y n -?=??=-=?? *001.732y y ≈= 作实际计算,问计算到10y 时误差为初始误差*00y y -的多少 这个计算过程数值稳定吗 ? . (14分)解答下列问题 1. 若2()63f x x =+,则[1,2,3]f 和[1,2,34]f ,的值分别是多少? 2. 1012 . (10分) 设f 在互易节点i x 上的值()()0,1,....i i f f x i n ==。试证明:f 在节点i x 上n 次最小二乘拟合多项式()n p x 与f 在节点i x 上的n 次Lagrange 插值多项式()n L x 一致,()()=n n p x L x 。 . (12分) 按代数精度的定义,构造下列形式的求积公式(即确定参数,A B ,α): Gauss 型求积公式。 . (14分) 已知线性代数方程组Ax=b 为: (1) 用顺序高斯消去法求解方程组Ax=b ; ????????????????=????????????????????????????????-----n n n n n n n n b b b b x x x x d u u u v d v d v d 121121121112211000000 .0)/(,0,11,,,≠-≠∑-=n i i i i n i i i i i d v u d d b v u d 已知其中

数值分析实验题目及解答

内容包括: 实验题目1:算法的数值稳定性实验 实验题目2:LU分解实验 实验题目3:三次样条插值外推样条实验 实验题目4:第二类Fredholm 积分方程实验实验题目5:M级显式R_K法

实验题目:算法的数值稳定性实验 实验内容:计算积分()1 0()d 1515n x I n x a x ==+? (n=1,2,…,20) 易得到下面递推公式 ()()1 1I n aI n n =--+ 并有估计式 ()() ()() 1 1 111I n a n a n << +++ 计算方法: 算法一:采用下面递推公式计算: ()()1 1I n aI n n =--+ ()1,2,,20 n = 取初值()116 0ln ln 15a I a +== 算法二: 采用下面递推公式计算: ()()111I n I n a n ??-= -+???? ()20,19,,1 n =

结果分析:(分析哪个好哪个不好,原因是什么) 我觉得算法二比较好, 原因一:根据式 ()() ()() 1 1 111I n a n a n << +++得知,I(n)不可能小于 零,而算法一的计算结果有部分结果小于零。原因二:对算法一记初始误差 ε0=/I 0-I(0)/>0; 则εn =/I n -I(n)/=a/I n-1-I(n-1)/=a n *ε0 由此可知,当n=20时, ε20把ε0放大了a 20倍,其结果造成严重的。 而对于算法二^ ^ 11n n a εε-= ,…, ^ ^ 01 n n a εε=,尽管有初始误差^ 20ε,但随着计算的进程,这个误差的影响不断减小。 附:源程序:(把源程序附上) 算法一程序: >> format long >> a=15;I=log(16/15); for n=1:20 n I=-a*I+1/n end 算法二程序: >> format long >> a=15;I=31/10080; >> for n=20:-1:1 n I I=1/a*(-I+1/n); End

数值分析试卷及其答案1

1. 已知325413.0,325413*2*1==X X 都有6位有效数字,求绝对误差限。(4分) 解: 由已知可知6 5.0102 1 ,0,6,10325413.0016*1=?= =-=?=ε绝对误差限n k k X 2分 620*2102 1 ,6,0,10325413.0-?= -=-=?=ε绝对误差限n k k X 2分 2. 已知?? ???=0 01 A 220- ?????440求21,,A A A ∞ (6分) 解: {}, 88,4,1max 1==A 1分 {}, 66,6,1max ==∞A 1分 () A A A T max 2λ= 1分 ?????=0 1 A A T 4 2 ???? ? -420?????0 01 2 20 - ???? ?440= ?????0 01 80 ???? ?3200 2分 {}32 32,8,1max )(max ==A A T λ

1分 24322==A 3. 设32)()(a x x f -= (6分) ① 写出f(x)=0解的迭代格式 ② 当a 为何值时,)(1k k x x ?=+ (0,1……)产生的序列{}k x 收敛于 2 解: ①迭代格式为: x a x x x a x a x x a x x x f x f x x k k k k k k k k k k 665)(665)(6)()(')(2 2 32 1 += +=---=-=+? 3 分 ②时迭代收敛即当222,112 10)2(',665)('2<<-<-=-= a a x a x ?? 3分 4. 给定线性方程组,其中:?? ?=13A ?? ?2 2,?? ? ???-=13b 用迭代公式 )()()()1(k k k Ax b x x -+=+α(0,1……)求解,问取什么实数α ,可使 迭代收敛 (8分) 解: 所给迭代公式的迭代矩阵为?? ? --???--=-=ααααα21231A I B 2分

数值分析整理版试题及答案

例1、 已知函数表 求()f x 的Lagrange 二次插值多项式和Newton 二次插值多项式。 解: (1) 故所求二次拉格朗日插值多项式为 (2)一阶均差、二阶均差分别为 例2、 设2 ()32f x x x =++,[0,1]x ∈,试求()f x 在[0,1]上关于()1x ρ=,{}span 1,x Φ=的最佳平 方逼近多项式。 解: 若{}span 1,x Φ=,则0()1x ?=,1()x x ?=,且()1x ρ=,这样,有 所以,法方程为

011231261192 34a a ??????????=?????????? ?????????? ,经过消元得012311 62110123a a ??? ???????=???????????????????? 再回代解该方程,得到14a =,011 6 a = 故,所求最佳平方逼近多项式为* 111()46S x x =+ 例3、 设()x f x e =,[0,1]x ∈,试求()f x 在[0,1]上关于()1x ρ=,{}span 1,x Φ=的最佳平方逼近 多项式。 解: 若{}span 1,x Φ=,则0()1x ?=,1()x x ?=,这样,有 所以,法方程为 解法方程,得到00.8732a =,1 1.6902a =, 故,所求最佳平方逼近多项式为 例4、 用4n = 的复合梯形和复合辛普森公式计算积分1 ? 。 解: (1)用4n =的复合梯形公式 由于2h =,( )f x =()121,2,3k x k k =+=,所以,有 (2)用4n =的复合辛普森公式 由于2h =,( )f x =()121,2,3k x k k =+=,()12 220,1,2,3k x k k + =+=,所以,有 例5、 用列主元消去法求解下列线性方程组的解。 解:先消元 再回代,得到33x =,22x =,11x = 所以,线性方程组的解为11x =,22x =,33x = 例6、 用直接三角分解法求下列线性方程组的解。 解: 设 则由A LU =的对应元素相等,有 1114u = ,1215u =,1316u =, 2111211433l u l =?=,3111311 22 l u l =?=, 2112222211460l u u u +=?=-,2113232311 545l u u u +=?=-,

数值计算方法期末考精彩试题

1. 已知函数 21 1y x = +的一组数据: 求分段 线性插值函数,并计算 () 1.5f 的近似值. 计算题1.答案 1. 解 [] 0,1x ∈, ()1010.510.50110x x L x x --= ?+?=---% []1,2x ∈,()210.50.20.30.81221x x L x x --=?+?=-+--% 所以分段线性插值函数为 ()[][]10.50,10.80.31,2x x L x x x ?-∈?=? -∈??% ()1.50.80.3 1.50.35 L =-?=% 4. 写出梯形公式和辛卜生公式,并用来分别计算积分1 01 1dx x +?. 计算题4.答案 4 解 梯形公式 ()()()2b a b a f x dx f a f b -≈ ?+???? 应用梯形公式得 1 01111 []0.75121011dx x ≈+=+++? 辛卜生公式为

确定下列求积公式中的待定系数,并证明确定后的求积公式具有3次代数精确度 ()()()() 1010h h f x dx A f h A f A f h --=-++? 证明题答案

故 ( )()()()40333h h h h f x dx f h f f h -= -++? 具有三次代数精确度。 1.设 3 2 01219 (), , 1, 44f x x x x x ==== (1)试求()f x 在 19,44???? ??上的三次Hermite 插值多项式()x H 使满足''11()(), 0,1,2,... ()()j j H x f x j H x f x === () x H 以升幂形式给出。 (2)写出余项()()()R x f x H x =-的表达式 计算题1.答案 1、(1) ()32142632331 22545045025x x x x H =- ++- (2) ()522191919()(1)(),()(,) 4!164444R x x x x x ξξξ-=---=∈ 3.试确定常数A ,B ,C 和 a ,使得数值积分公式 有尽可能高的代数精度。试问所得的数值积分公式代数精度是多少?它是否为Gauss 型的? 计算题3.答案

华南理工大学数值分析试题

华南理工大学研究生课程考试 《数值分析》试卷B (2015 年 1 月 9 日) 师教课任 注意事项: 1. 考前请将密封线内各项信息填写清楚; 2. 所有答案请按要求填写在本试卷上; 3. 课程代码:S0003004; 4. 考试形式:闭卷; 5. 考生类别:硕士研究生; 6. 本试卷共八大题,满分 100分,考试时间为150分钟。 线 一?单项选择题(每小题2分,共10分) 1 ?设有某数x,则x的具有四位有效数字且绝对误差限是0. 5 10 5的近似值 应是( )° (A) 0.693 (B) 0.6930 (C) 0.0693 (D) 0.06930 业专 院学 ) 题 答 不 内 线 封 密 { 2 ?选择数值稳定的算法是为了() (A)简化计算步骤 (C)节省存储空间 (B)控制舍入误差的积累 (D)减小截断误差 3.如果对不超过m次的多项式,求积公式 式具有( )次代数精度。 (A)至少 m (B) m b f (x)dx a (C) 不足m A k f (x k)精确成立,则该求积公 k 0 (D)多于m 号学名姓 4.为使两点数值求积公式 1 1 f(x)dx f(X。) f (x1)具有最高次代数精度, 则求积节点应为( )° (A)x°,X1 任意(B) X。1,X1 1 (C) X。- ,x1 3 _3 3 (D) x o 1 1 ,X1 2 2 密 5.在下列求解常微分方程初值问题的数值方法中, (A) Euler 公式(B) (C) 3 阶 Runge— Kutta 公式(D) 4 () 的局部截断误差为 梯形公式 阶 Runge— Kutta 公式 O(h 3)。

数值分析第五版答案

第一章 绪论 3.下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指 出它们是几位有效数字:*1 1.1021x =,*20.031x =, *3385.6x =, * 456.430x =,*57 1.0.x =? 解:*1 1.1021x =是五位有效数字; *20.031x =是二位有效数字; *3385.6x =是四位有效数字; *456.430x =是五位有效数字; *57 1.0.x =?是二位有效数字。 4.利用公式(2.3)求下列各近似值的误差限:(1) ***124x x x ++,(2) ***123x x x ,(3) **24/x x . 其中****1234 ,,,x x x x 均为第3题所给的数。 解: *4 1* 3 2* 13* 3 4* 1 51 ()1021()1021()1021()1021()102 x x x x x εεεεε-----=?=?=?=?=? *** 124***1244333 (1)()()()() 1111010102221.0510x x x x x x εεεε----++=++=?+?+?=? *** 123*********123231132143 (2)() ()()() 111 1.10210.031100.031385.610 1.1021385.610222 0.215 x x x x x x x x x x x x εεεε---=++=???+???+???≈

** 24**** 24422 * 4 33 5 (3)(/) ()() 11 0.0311056.430102256.43056.430 10x x x x x x x εεε---+≈ ??+??= ?= 5计算球体积要使相对误差限为1,问度量半径R 时允许的相对误差限是多少? 解:球体体积为34 3 V R π= 则何种函数的条件数为 2 3'4343 p R V R R C V R ππ=== (*)(*)3(*)r p r r V C R R εεε∴≈= 又(*)1r V ε= 故度量半径R 时允许的相对误差限为1 (*)10.333 r R ε=?≈ 6.设028Y = ,按递推公式1n n Y Y -= (n=1,2,…) 计算到100Y 27.982≈(5位有效数字),试问计算100Y 将有多大误差? 解:1n n Y Y -=- 10099Y Y ∴=- 9998Y Y = 9897Y Y =-…… 10Y Y =- 依次代入后,有1000100Y Y =- 即1000Y Y = 27.982, 100027.982Y Y ∴=-

数值分析期末试题

数值分析期末试题 一、填空题(20102=?分) (1)设??? ? ? ??? ??---=28 3 012 251A ,则=∞ A ______13_______。 (2)对于方程组?? ?=-=-3 4101522121x x x x ,Jacobi 迭代法的迭代矩阵是=J B ?? ? ? ??05.25.20。 (3)3*x 的相对误差约是*x 的相对误差的 3 1倍。 (4)求方程)(x f x =根的牛顿迭代公式是) ('1)(1n n n n n x f x f x x x +-- =+。 (5)设1)(3 -+=x x x f ,则差商=]3,2,1,0[f 1 。 (6)设n n ?矩阵G 的特征值是n λλλ,,,21 ,则矩阵G 的谱半径=)(G ρi n i λ≤≤1max 。 (7)已知?? ? ? ??=1021 A ,则条件数=∞ )(A Cond 9 (8)为了提高数值计算精度,当正数x 充分大时,应将)1ln(2 -- x x 改写为 )1ln(2 ++ -x x 。 (9)n 个求积节点的插值型求积公式的代数精确度至少为1-n 次。 (10)拟合三点))(,(11x f x ,))(,(22x f x ,))(,(33x f x 的水平直线是)(3 1 3 1 ∑== i i x f y 。 二、(10分)证明:方程组? ?? ??=-+=++=+-1 211 2321321321x x x x x x x x x 使用Jacobi 迭代法求解不收敛性。 证明:Jacobi 迭代法的迭代矩阵为 ???? ? ?????---=05 .05 .01015.05.00J B J B 的特征多项式为

华南理工大学数值分析试题-14年下-C

华南理工大学研究生课程考试 《数值分析》试卷C (2015年1月9日) 1. 考前请将密封线内各项信息填写清楚; 所有答案请按要求填写在本试卷上; 课程代码:S0003004; 4. 考试形式:闭卷; 5. 考生类别:硕士研究生; 本试卷共八大题,满分100分,考试时间为150分钟。 一、(12分)解答下列问题: 1)设近似值0x >,x 的相对误差为δ,试证明ln x 的绝对误差近似为δ。 2)利用秦九韶算法求多项式 542()681p x x x x x =-+-+ 在3x =时的值(须写出计算形式),并统计乘法次数。 (12分)解答下列问题: 1)设()235f x x =+,求[]0,1,2f 和[]0,1,2,3f 。 2)利用插值方法推导出恒等式: 33220,0[]j j i i x j i x i j =≠=-=-∑∏ 。

(1)设{}∞ =0)(k k x q 是区间[]1,0上带权1=ρ而最高次项系数为1的正交多项式族,其中1)(0=x q ,求1()q x 和2()q x 。 (2)求形如2y a bx =+的经验公式,使它与下列数据拟合: 四、(14分)对积分()10I f x dx = ?,试 (1)构造一个以012113,,424 x x x ===为节点的插值型求积公式; (2)指出所构造公式的代数精度; (3)用所得数值求积公式计算积分1 203x dx ?的精确值; (4)指出所得公式与一般的Newton-Cotes 型公式在形式上的重要区别。

(1)设?? ????=4321A ,计算1A 、()Cond A ∞和()A ρ。 (2)用列主元Gauss 消去法解方程组: 12312315410030.112x x x ????????????=????????????-?????? 六、(13分)对2阶线性方程组 11112212112222 a x a x b a x a x b +=??+=? (11220a a ≠ ) (1)证明求解此方程组的Jacobi 迭代与Gauss-Seidel 迭代同时收敛或同时发散; (2)当同时收敛时,试比较它们的收敛速度。

相关文档
最新文档