(合集)西南交通大学2005-2013研究生数值分析试题

西南交大 数值分析题库

考试目标及考试大纲 本题库的编纂目的旨在给出多套试题，每套试题的考查范围及难度配置均基于“水平测试”原则，按照教学大纲和教学内容的要求，通过对每套试题的解答，可以客观公正的评定出学生对本课程理论体系和应用方法等主要内容的掌握水平。通过它可以有效鉴别和分离不同层次的学习水平，从而可以对学生的学习成绩给出客观的综合评定结果。 本题库力求作到能够较为全面的覆盖教学内容，同时突显对重点概念、重点内容和重要方法的考查。考试内容包括以下部分： 绪论与误差：绝对误差与相对误差、有效数字、误差传播分析的全微分法、相对误差估计的条件数方法、数值运算的若干原则、数值稳定的算法、常用数值稳定技术。 非线性方程求解：方程的近似解之二分法、迭代法全局收敛性和局部收敛定理、迭代法误差的事前估计法和事后估计法、迭代过程的收敛速度、r 阶收敛定理、Aitken加速法、Ne w to n法与弦截法、牛顿局部收敛性、Ne w to n收敛的充分条件、单双点割线法（弦截法）、重根加速收敛法。 解线性方程组的直接法：高斯消元法极其充分条件、全主元消去法、列主元消去法、高斯-若当消元法、求逆阵、各种消元运算的数量级估计与比较、矩阵三角分解法、Doolittle 和Crout三角分解的充分条件、分解法的手工操作、平方根法、Cholesky分解、改进的平方根法(免去开方)、可追赶的充分条件及适用范围、计算复杂性比较、严格对角占优阵。 解线性方程组迭代法：向量和矩阵的范数、常用向量范数的计算、范数的等价性、矩阵的相容范数、诱导范数、常用范数的计算；方程组的性态和条件数、基于条件数误差估计与迭代精度改善方法；雅可比（Jacobi）迭代法、Gauss-Seidel迭代法、迭代收敛与谱半径的关系、谱判别法、基于范数的迭代判敛法和误差估计、迭代法误差的事前估计法和事后估计法；严格对角占优阵迭代收敛的有关结论；松弛法及其迭代判敛法。 插值法：插值问题和插值法概念、插值多项式的存在性和唯一性、插值余项定理；Lagrange插值多项式；差商的概念和性质、差商与导数之间的关系、差商表的计算、牛顿（Newton）插值多项式；差分、差分表、等距节点插值公式；Hermite插值及其插值基函数、误差估计、插值龙格（Runge）现象；分段线性插值、分段抛物插值、分段插值的余项及收敛性和稳定性；样条曲线与样条函数、三次样条插值函数的三转角法和三弯矩法。 曲线拟合和函数逼近：最小二乘法原理和多项式拟合、函数线性无关概念、法方程有唯一解的条件、一般最小二乘法问题、最小二乘拟合函数定理、可化为线性拟合问题的常见函数类；正交多项式曲线拟合、离散正交多项式的三项递推法。最佳一致逼近问题、最佳一致逼近多项式、切比雪夫多项式、切比雪夫最小偏差定理、切比雪夫多项式的应用(插值余项近似极小化、多项式降幂)。本段加黑斜体内容理论推导可以淡化，但概念需要理解。 数值积分与微分：求积公式代数精度、代数精度的简单判法、插值型求积公式、插值型求积公式的代数精度；牛顿一柯特斯(Newton-Cotes)公式、辛卜生（Simpson）公式、几种低价牛顿一柯特斯求积公式的余项；牛顿一柯特斯公式的和收敛性、复化梯形公式及其截断误差、复化Simpson公式及其截断误差、龙贝格（Romberg）求积法、外推加速法、高斯型求积公式、插值型求积公式的最高代数精度、高斯点的充分必要条件。正交多项式的构造方法、高斯公式权系数的建立、Gauss-Legendre公式的节点和系数。本段加黑斜体内容理论推导可以淡化，但概念需要理解。 常微分方程数值解：常微分方程初值问题数值解法之欧拉及其改进法、龙格—库塔法、阿当姆斯方法。

最新第六章习题答案-数值分析

第六章习题解答 2、利用梯形公式和Simpson 公式求积分2 1 ln xdx ? 的近似值，并估计两种方法计算值的最大 误差限。 解：①由梯形公式： 21ln 2 ()[()()][ln1ln 2]0.3466222 b a T f f a f b --= +=+=≈ 最大误差限 3''2 ()111 ()()0.0833******** T b a R f f ηη-=-=≤=≈ 其中，(1,2)η∈ ②由梯形公式： 13()[()4()()][ln14ln()ln 2]0.38586262 b a b a S f f a f f b -+= ++=++≈ 最大误差限 5(4)4()66 ()()0.0021288028802880 S b a R f f ηη-=-=≤≈， 其中，(1,2)η∈。 4、推导中点求积公式 3''()()()()() ()224 b a a b b a f x dx b a f f a b ξξ+-=-+<

2009哈工大级研究生《数值分析》试卷

2009级研究生《数值分析》试卷 一.(6分) 已知描述某实际问题的数学模型为x y y x y x u 223),(+=，其中，y x ,由 统计方法得到，分别为4,2==y x ,统计方法的误差限为0.01,试求出u 的误差限 )(u ε和相对误差限)(u r ε. 二.(6分) 已知函数13)(3+=x x f 计算函数)(x f 的2阶均差]2,1,0[f ,和4阶均差]4,3,2,1,0[f . 三.(6分)试确定求积公式: )]1(')0('[12 1 )]1()0([21)(10f f f f dx x f -++≈?的代数精 度. 四.(12分) 已知函数122)(2 3 -++=x x x x f 定义在区间[-1,1]上,在空间 },,1{)(2x x Span x =Φ上求函数)(x f 的最佳平方逼近多项式. 其中,权函数1)(=x ρ,15 4 ))(),((,1532))(),((,34))(),((210-==-=x x f x x f x x f ???. 五.(16分) 设函数)(x f 满足表中条件: (1) 填写均差计算表(标有*号处不填): (2) 分别求出满足条件)2,1,0(),()(),()(22===k x f x N x f x L k k k k 的 2次 Lagrange 和 Newton 差值多项式.

(3) 求出一个四次插值多项式)(4x H ,使其满足表中所有条件.并用多项式降幂形式表示. 六.(16分) (1). 用Romberg 方法计算?3 1 dx x ,将计算结果填入下表(*号处不填). (2). 试确定三点 Gauss-Legender 求积公式?∑-=≈1 1 2 )()(k k k x f A dx x f 的Gauss 点k x 与系数 k A ,并用三点 Gauss-Legender 求积公式计算积分: ?3 1dx x . 七.(14分) (1) 证明方程02ln =--x x 在区间(1,∞)有一个单根.并大致估计单根的取值范围. (2) 写出Newton 迭代公式,并计算此单根的近似值.(要求精度满足: 5 110||-+<-k k x x ). 八. (12分) 用追赶法求解方程组: ???? ?? ? ??=??????? ????????? ??022112111131124321x x x x 的解. 九. (12分) 设求解初值问题???==0 0)() ,('y x y y x f y 的计算格式为: )],(),([111--+++=n n n n n n y x bf y x af h y y ,假设11)(,)(--==n n n n y x y y x y ,试确定参数b a ,的值,使该计算格式的局部截断误差为二阶,即截断部分为: )(3h o .

研究生《数值分析》练习题

硕士研究生 《数值分析》练习题 一、判断题 1、用Newton 切线法求解非线性线性方程可以任选初值。 （ ） 2、求解非线性线性方程，Newton 切线法比弦截法迭代次数多。 （ ） 3、若n n A R ?∈非奇异，用Jacobi 迭代法求解线性方程组Ax b =必收敛。（ ） 4、Lagrange 插值法与Newton 插值法得到同一个插值多项式。 （ ） 二、填空题 1、近似数 3.14108937a =关于π具 位有效数字。 2、双点弦截法具有 阶收敛速度。 3、求方程x x e =根的单点弦截法迭代公式是 。 4、设2112A ?? = ? ?? ? ，则()A ρ= 。 5、设,0,1,2,3i x i =是插值基点，,0,1,2,3i l i =是对应的三次Lagrange 插值基函数，则()()3 3012i i i x l =-=∑ 。 6、由下数据表确定的代数插值多项式的不超过 次。 7、若()8754321f x x x x =+-+，则差商[]0,1,2,,8f = 。 8、拟合三点()()()0,1,1,3,2,2A B C 的直线是y = 。 三、分析与计算题 1、设()14,2,3515T A x -??==-?? -?? ，求∞=,2,1,,p x A p p 和()1A cond 。

2、1001012,20253A x -???? ? ? == ? ? ? ?-???? ，试计算p p x A ,，p=1,2,∞，和1)(A cond 。 3、线性方程组,0Ax b b =≠，用Jacobi 迭代法是否收敛，为什么？其中 122111221A -?? ?=-- ? ?--?? 。 4、线性方程组,0Ax b b =≠，用Jacobi 迭代法是否收敛，为什么？其中 2-11=11111-2A ?? ???? ???? 。 5、已知函数表如下: ⑴ ()111.75ln11.75L ≈、估计截断误差并说明结果有几位有效数字； ⑵ ()211.75ln11.75N ≈、估计截断误差并说明结果有几位有效数字。 6、已知函数表 如下: ⑴用Lagrange 插值法求ln0.55的近似值()10.55N 、估计截断误差并说明结果的有效数字； ⑵用 Newton 插值法求ln0.55的近似值()20.55N 、估计截断误差并说明结果的有效数字。 7、已知数据如下，求满足条件的Hermite 插值多项式。

西南交大数值分析题库填空

一. 填空 2．Gauss型求积公式不是插值型求积公式。（限填“是”或“不是”） 3. 设l k(x)是关于互异节点x0, x1,…, x n, 的Lagrange 插值基函数，则 0 m=1,2,…,n 5．用个不同节点作不超过次的多项式插值，分别采用Lagrange插值方法与Newton插值方法所得多项式相等(相等, 不相等)。 。 7. n个不同节点的插值型求积公式的代数精度一定会超过n-1次 8.f(x)=ax7+x4+3x+1,f[20, 21,…,27]= a，f [20, 21,…,28]= 0 10设 (i=0,1,…,n)，则＝ _x_ , 这里（x i x j,ij, n2）11.设称为柯特斯系数 则=______1____ 12采用正交多项式拟合可避免最小二乘或最佳平方逼近中常见的_法方程组病态___问题。 13辛卜生（Simpson）公式具有___3____次代数精度。 14 牛顿插商与导数之间的关系式为： 15试确定[0,1]区间上2x3的不超过二次的最佳一致逼近多项式p(x), 该多项式唯一否？答：p(x)=(3/2)x, ; 唯一。 17.给定方程组记此方程组的Jacobi迭代矩阵为B J=(a ij)33，则a23= -1； ,且相应的Jacobi迭代序列是__发散_____的。 18.欧拉预报--校正公式求解初值问题的迭代格式(步长为h) ，此方法是阶方法。 ，此方法是 2阶方法。 19. 2n阶Newton-Cotes公式至少具有2n+1次代数精度。 20.设，则关于的 ||f|| =1 21矩阵的LU分解中L是一个 _为单位下三角阵，而U是一个上三角阵____。 22.设y=f (x1,x2) 若x1,x2,的近似值分别为x1*, x2*,令y*=f(x1*,x2*)作为y的近似值,其绝对误差限的估计式为: ||f(x1*,x2*)|x1-x*1|+ |f(x1*,x2*)|x2- x*2| 23设迭代函数(x)在x*邻近有r（1）阶连续导数，且x* = (x*)，并且有(k) (x*)=0 (k=1,…,r-1)，但(r) (x*)0，则x n+1=(x n)产生的序列{ x n }的收敛阶数为___r___ 24设公式为插值型求积公式，则, 且=b-a 25称微分方程的某种数值解法为p阶方法指的是其局部截断误差 为O(h p+1)。 26.设x0, x1,x2是区间[a, b]上的互异节点，f(x)在[a, b]上具有各阶导数，过

电子科技大学数值分析研究生期末考试习题一

习 题 请尽可能提供程序 1.用二分法求方程012=--x x 的正根，要求误差05.0<。 2. 为求方程0123=--x x 在5.10=x 附近的一个根，设将方程改写成下列等价形式，并建立相应的迭代公式： 1）2/11x x +=，迭代公式21/11k k x x +=+；2）231x x +=，迭代公式3211k k x x +=+； 3）1 12-=x x ，迭代公式1/11-=+k k x x ；4）132-=x x ，迭代公式131-=+k k x x 。 试分析每种迭代公式的收敛性。 3. 给定函数)(x f ，设对一切x ，)(x f '存在且M x f m ≤'≤<)(0，证明对于范围M /20<<λ内的任意定数λ，迭代过程)(1k k k x f x x λ-=+均收敛于)(x f 的根*x 。 4.设a 为正整数，试建立一个求 a 1的牛顿迭代公式，要求在迭代公式中不含有除法运算，并考虑公式的收敛性。请提供程序。 5.用Gauss 消去法求解方程组： ???? ? ??-=????? ??????? ??----50312131 2111321x x x （请提供程序） 用列主元Gauss 消去法求解下列方程组： （1）???? ? ??=????? ??????? ??13814142210321321x x x （请提供程序） 6.用追赶法解三对角方程组b Ax =，其中 ????????????????--------=210001 2100012100012100012A ，??????? ?????????=00001b 。 7.设n n R P ?∈且非奇异，又设x 为n R 上一向量范数，定义Px x p =。试证明p x 是n R 上向量的一种范数。 8.用平方根法（Cholesky 分解）求解方程组：

西南交大数值分析题库积分微分方程

用复化梯形公式计算积分 1 ()f x dx ?,要把区间[0,1]一般要等分 41 份才能保 证满足误差小于0.00005的要求（这里(2) () 1f x ∞ ≤） ；如果知道(2) ()0f x >，则 用复化梯形公式计算积分1 ()f x dx ? 此实际值 大 (大，小)。 在以1 0((),())()(),(),()[0,1]g x f x xf x g x dx f x g x C = ∈?为内积的空间C[0,1] 中，与非零常数正交的最高项系数为1的一次多项式是 2 3 x - 3. (15分)导出用Euler 法求解 (0)1y y y λ'=??=? 的公式, 并证明它收敛于初值问题的精确解 解 Euler 公式 11,1,,,k k k x y y h y k n h n λ--=+== L -----------(5分) ()()1011k k k y h y h y λλ-=+==+L ------------------- (10分) 若用复化梯形求积公式计算积分1 x I e dx = ? 区间[0,1]应分 2129 等分，即要 计算个 2130 点的函数值才能使截断误差不超过 71 102 -?；若改用复化Simpson 公式，要达到同样精度区间[0,1]应分12 等分，即要计算个 25 点的函数值 1．用Romberg 法计算积分 2 3 2 x e dx -? 解 []02()()2b a T f a f b -= += 9.6410430E-003 10221()222 b a a b T T f -+=+= 5.1319070E-003 10 022243 T T S -= = 4.6288616E-003 22T = 4.4998E-003 21 122243 T T S -= = 4.E-003 10 02221615 S S C -= = 4.6588636E-003 32T = 4.7817699E-003 32 222243 T T S -= = 4.1067038E-003

第六章习题答案数值分析.docx

第六章习题解答 2 2、利用梯形公式和 Simpson 公式求积分 ln xdx 的近似值， 并估计两种方法计算值的最大 1 误差限。 解：①由梯形公式： T ( f ) b a [ f (a) f (b)] 2 1 [ln1 ln 2] ln 2 0.3466 2 2 2 最大误差限 R ( f ) (b a)3 f '' ( ) 1 1 1 0.0833 T 12 12 2 12 12 其中， (1,2) ②由梯形公式： b a 4 f ( b a f (b)] 1 4ln( 3 ln 2] 0.3858 S( f ) [ f (a) ) [ln1 ) 6 2 6 2 最大误差限 R S ( f ) (b a)5 f (4) ( ) 6 6 0.0021， 2880 2880 4 2880 其中， (1,2) 。 4、推导中点求积公式 f ( x)dx (b a) f ( a b ) (b a) 3 (a b) b a 2 24 证明： 构造一次函数 P （ x ），使 P a 2 b f a b , P ' ( a b ) f ' ( a b ), P '' ( x) 0 2 2 2 则，易求得 P( x) f ' ( a b )( x a b ) f ( a b ) 2 2 2 且 P(x)dx f ' ( a b )( x a b ) f ( a b ) dx b b a a 2 2 2 f ( a b )dx (b a) f ( a b ) ，令 P(x)dx I ( f ) b b a 2 2 a 现分析截断误差：令 r ( x) f ( x) P(x) f ( x) f ' ( a b )( x a b ) f ( a b ) 2 2 2 由 r ' ( x) f ' (x) f ' ( a b ) 易知 x a 2 b 为 r (x) 的二重零点， 2 a b )2 ， 所以可令 r (x) ( x)( x 2

北航2010-2011年研究生数值分析期末模拟试卷1-3

数值分析模拟试卷1 一、填空（共30分，每空3分） 1 设??? ? ??-=1511A ，则A 的谱半径=)(a ρ______，A 的条件数=________. 2 设 ,2,1,0,,53)(2==+=k kh x x x f k ，则],,[21++n n n x x x f ＝________, ],,[321+++n n n n x x x x f ，＝________. 3 设?????≤≤-++≤≤+=2 1,121 0,)(2 323x cx bx x x x x x S ，是以0，1，2为节点的三次样条函数，则b=________,c=________. 4 设∞=0)]([k k x q 是区间［0，1］上权函数为x x =)(ρ的最高项系数为1的正交多项式族，其中1)(0=x q ，则 ?=1 )(dx x xq k ________，=)(2 x q ________. 5 设???? ??????=11001a a a a A ，当∈a ________时，必有分解式，其中L 为下三角阵，当 其对角线元素)3,2,1(=i L ii 满足条件________时，这种分解是唯一的. 二、（14分）设4 9,1,41,)(2102 3 === =x x x x x f , （1）试求)(x f 在]4 9,41[上的三次Hermite 插值多项式)(x H 使满足 2,1,0),()(==i x f x H i i ，)()(11x f x H '='. （2）写出余项)()()(x H x f x R -=的表达式. 三、（14分）设有解方程0cos 2312=+-x x 的迭代公式为n n x x cos 3 2 41+ =+， （1） 证明R x ∈?0均有? ∞ →=x x n x lim （? x 为方程的根）； （2） 取40=x ，用此迭代法求方程根的近似值，误差不超过，列出各次迭代值； （3）此迭代的收敛阶是多少？证明你的结论. 四、(16分) 试确定常数A ，B ，C 和，使得数值积分公式 有尽可能高的代数精度. 试问所得的数值积分公式代数精度是多少？它是否为Gauss 型的？

2012研究生数值分析课期末考试复习题及答案

一、填空 1. 设 2.3149541...x * =，取5位有效数字，则所得的近似值x= 2.3150 . 2.设一阶差商 ()()()21122114 ,321f x f x f x x x x --= = =---， ()()()322332 615 ,422f x f x f x x x x --= = =-- 则二阶差商 ()123,,______ f x x x =11/6 3. 设(2,3,1)T X =--, 则2||||X = 14 ，=∞||||X 3 。p49 4. 4．求方程 2 1.250x x --= 的近似根，用迭代公式 1.25x x =+，取初始值 01 x =， 那么 1______x =。 1.5 5．解初始值问题 00 '(,)()y f x y y x y =?? =?近似解的梯形公式是 1______k y +≈。 ()()[]11,,2 ++++k k k k k y x f y x f h y 6、 1151A ??= ? -??，则A 的谱半径 ＝ 6 。 7、设 2()35, , 0,1,2,... , k f x x x kh k =+== ，则 []12,,n n n f x x x ++= —————— ————3 和 []123,,,n n n n f x x x x +++= _______________0_____ 。 8、 若线性代数方程组AX=b 的系数矩阵A 为严格对角占优阵，则雅可比迭代和高斯-塞德尔迭代都 收敛 。 9、解常微分方程初值问题的欧拉（Euler ）方法的局部截断误差为_______O(h ） ___。

西南交通大学2018-2019数值分析Matlab上机实习题

数值分析2018-2019第1学期上机实习题 f x，隔根第1题.给出牛顿法求函数零点的程序。调用条件：输入函数表达式() a b，输出结果：零点的值x和精度e，试取函数 区间[,] ，用牛顿法计算附近的根，判断相应的收敛速度，并给出数学解释。 1.1程序代码： f=input('输入函数表达式：y=','s'); a=input('输入迭代初始值：a='); delta=input('输入截止误差：delta='); f=sym(f); f_=diff(f); %求导 f=inline(f); f_=inline(f_); c0=a; c=c0-f(c0)/f_(c0); n=1; while abs(c-c0)>delta c0=c; c=c0-f(c0)/f_(c0); n=n+1; end err=abs(c-c0); yc=f(c); disp(strcat('用牛顿法求得零点为',num2str(c))); disp(strcat('迭代次数为',num2str(n))); disp(strcat('精度为',num2str(err))); 1.2运行结果： run('H:\Adocument\matlab\1牛顿迭代法求零点\newtondiedai.m') 输入函数表达式：y=x^4-1.4*x^3-0.48*x^2+1.408*x-0.512 输入迭代初始值：a=1 输入截止误差：delta=0.0005 用牛顿法求得零点为0.80072 迭代次数为14 精度为0.00036062 牛顿迭代法通过一系列的迭代操作使得到的结果不断逼近方程的实根，给定一个初值，每经过一次牛顿迭代，曲线上一点的切线与x轴交点就会在区间[a,b]上逐步逼近于根。上述例子中，通过给定初值x=1，经过14次迭代后，得到根为0.80072，精度为0.00036062。

研究生数值分析试卷

2005~2006学年第一学期硕士研究生期末考试试题（A 卷） 科目名称：数值分析 学生所在院： 学号： 姓名： 注意：所有的答题内容必须答在答题纸上，凡答在试题或草稿纸上的一律无效。 一、（15分）设求方程 0cos 2312=+-x x 根的迭代法 k k x x cos 3 2 41+=+ (1) 证明对R x ∈?0,均有*lim x x k k =∞ →,其中*x 为方程的根. (2) 此迭代法收敛阶是多少? 证明你的结论. 二、（12分）讨论分别用Jacobi 迭代法和Gauss-Seidel 迭代法求解下列方程组的收敛性。 ??? ??=++-=++=-+. 022,1, 122321 321321x x x x x x x x x 三、（8分）若矩阵??? ? ? ??=a a a a A 000002，说明对任意实数0≠a ，方程组b AX =都是非病态的。（范数用∞?） 四、（ 求)(x f 的Hermite 插值多项式)(3x H ，并给出截断误差)()()(3x H x f x R -=。 五、（10分）在某个低温过程中，函数 y 依赖于温度x （℃）的试验数据

为 已知经验公式的形式为 2bx ax y += ，试用最小二乘法求出 a ，b 。 六、（12分）确定常数 a ，b 的值，使积分 [ ] dx x b ax b a I 2 1 1 2 ),(?--+= 取得最小值。 七、（14分）已知Legendre(勒让德)正交多项式)(x L n 有递推关系式： ?? ? ? ???=+-++===-+),2,1()(1)(112)()(, 1)(1110 n x L n n x xL n n x L x x L x L n n n 试确定两点的高斯—勒让德（G —L ）求积公式 ? -+≈1 1 2211)()()(x f A x f A dx x f 的求积系数和节点，并用此公式近似计算积分 ?=2 11 dx e I x 八、（14分）对于下面求解常微分方程初值问题 ?????==0 0)() ,(y x y y x f dx dy 的单步法： ??? ? ??? ++==++=+) ,() ,()2 121(1 21211 hk y h x f k y x f k k k h y y n n n n n n

数值分析最佳习题(含答案)

第一章 绪论 姓名 学号 班级 习题主要考察点：有效数字的计算、计算方法的比较选择、误差和误差限的计算。 1 若误差限为5105.0-?,那么近似数有几位有效数字（有效数字的计算） 解：2*103400.0-?=x ，325*102 1102 1---?=?≤-x x 故具有3位有效数字。 2 14159.3=π具有4位有效数字的近似值是多少（有效数字的计算） 解：10314159.0?= π，欲使其近似值*π具有4位有效数字，必需 41*102 1 -?≤-ππ，3*3102 1102 1--?+≤≤?-πππ，即14209.314109.3*≤≤π 3 已知2031.1=a ，978.0=b 是经过四舍五入后得到的近似值，问b a +， b a ?有几位有效数字（有效数字的计算） 解：3*1021 -?≤-a a ，2*102 1-?≤-b b ，而1811.2=+b a ，1766.1=?b a 2123****102 1 10211021)()(---?≤?+?≤ -+-≤+-+b b a a b a b a 故b a +至少具有2位有效数字。 2123*****102 1 0065.01022031.1102978.0)()(---?≤=?+?≤ -+-≤-b b a a a b b a ab

故b a ?至少具有2位有效数字。 4 设0>x ，x 的相对误差为δ，求x ln 的误差和相对误差（误差的计算） 解：已知δ=-* *x x x ，则误差为 δ=-= -* **ln ln x x x x x 则相对误差为 * * ** * * ln ln 1ln ln ln x x x x x x x x δ = -= - 5测得某圆柱体高度h 的值为cm h 20*=，底面半径r 的值为cm r 5*=， 已知cm h h 2.0||*≤-，cm r r 1.0||*≤-，求圆柱体体积h r v 2π=的绝对误差 限与相对误差限。（误差限的计算） 解：*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 相对误差限为 %420 1 20525) 5,20() 5,20(),(2 ==??≤ -ππv v r h v 6 设x 的相对误差为%a ,求n x y =的相对误差。（函数误差的计算） 解：%* *a x x x =-， )%(* **** *na x x x n x x x y y y n n n =-≤-= - 7计算球的体积，为了使体积的相对误差限为%1，问度量半径r 时允许的相对误差限为多大（函数误差的计算）

数值分析 第六章 习题

第六章 习 题 1. 计算下列矩阵的1A ，2A ，A ∞三种范数。 （1）1101A ???=????，（2）312020116A ????=??????? . 2. 用Jacobi 方法和Gauss-Seidel 迭代求解方程组 1231231 238322041133631236x x x x x x x x x ?+=??+?=??++=? 要求取(0)(0,0,0)T x =计算到(5)x ，并分别与精确解(3,2,1)T x =比较。 3. 用Gauss-Seidel 迭代求解 12312312 35163621122x x x x x x x x x ??=??++=???+=?? 以(0)(1,1,1)T x =?为初值，当(1)() 310k k x x +?∞?<时，迭代终止。 4. 已知方程组121122,2,x x b tx x b +=?? +=? （1）写出解方程组的Jacobi 迭代矩阵，并讨论迭代收敛条件。 （2）写出解方程组的Gauss-Seidel 迭代矩阵，并讨论迭代收敛条件. 5. 设有系数矩阵 122111221A ?????=?????? ， 211111112B ?????=??????? ， 证明：（1）对于系数矩阵A ，Jacobi 迭代收敛，而Gauss-Seidel 迭代不收敛. （2）对于矩阵B ，. 6. 讨论方程组 112233302021212x b x b x b ?????????????=??????????????????? 用Jacobi 迭代和Gauss-Seidel 迭代的收敛性；如果都收敛，比较哪种方法收敛更快.

2008级研究生数值分析试题

太原科技大学 2008级硕士研究生08/09学年第一学期 《数值分析》考试试卷 说明：1、Legendre 正交多项式)(x L n 有三项递推关系式： ?? ?? ???=+-++===-+ ,2,1)(1)(112)()(,1)(1110n x L n n x xL n n x L x x L x L n n n 2、Chebyshev 多项式)(x T n 有三项递推关系式： ?? ? ??=-===-+ ,2,1)()(2)()(,1)(1110n x T x xT x T x x T x T n n n 一、填空题：（每题4分，共20分） 1、设??? ? ??-=1511A ，则=∞)(A Cond 2、为提高数值计算精度，当x 充分小时，应将 x x sin cos 1-改写为 3、设)5()(2 -+=x a x x ?，要使)(1k k x x ?=+局部收敛到5* = x ，则a 的取值范围为 4、近似数235.0* =x 关于真值229.0=x 有 位有效数字。 5、设,1)(3 -+=x x x f 则差商=]3,2,1,0[f 二、（本题满分10分）用数值积分的方法建立求解初值问题b x a y a y y x f y a ≤≤==',)(),,(的Simpson 公式： )4(3 1111-+-++++=n n n n n f f f h y y 其中1,,1),,(+-==n n n i y x f f i i i ，11-+-=-=n n n n x x x x h . 三、（本题满分15分）设要用Gauss-Seidel 迭代法求解下列线性方程组

数值分析2016上机实验报告

序言 数值分析是计算数学的范畴，有时也称它为计算数学、计算方法、数值方法等，其研究对象是各种数学问题的数值方法的设计、分析及其有关的数学理论和具体实现的一门学科，它是一个数学分支。是科学与工程计算（科学计算）的理论支持。许多科学与工程实际问题（核武器的研制、导弹的发射、气象预报）的解决都离不开科学计算。目前，试验、理论、计算已成为人类进行科学活动的三大方法。 数值分析是计算数学的一个主要部分,计算数学是数学科学的一个分支,它研究用计算机求解各种数学问题的数值计算方法及其理论与软件实现。现在面向数值分析问题的计算机软件有：C,C++,MATLAB,Python,Fortran等。 MATLAB是matrix laboratory的英文缩写，它是由美国Mathwork公司于1967年推出的适合用于不同规格计算机和各种操纵系统的数学软件包，现已发展成为一种功能强大的计算机语言，特别适合用于科学和工程计算。目前，MATLAB应用非常广泛，主要用于算法开发、数据可视化、数值计算和数据分析等，除具备卓越的数值计算能力外，它还提供了专业水平的符号计算，文字处理，可视化建模仿真和实时控制等功能。 本实验报告使用了MATLAB软件。对不动点迭代，函数逼近（lagrange插值，三次样条插值，最小二乘拟合），追赶法求解矩阵的解，4RungeKutta方法求解，欧拉法及改进欧拉法等算法做了简单的计算模拟实践。并比较了各种算法的优劣性，得到了对数值分析这们学科良好的理解，对以后的科研数值分析能力有了极大的提高。

目录 序言 (1) 问题一非线性方程数值解法 (3) 1.1 计算题目 (3) 1.2 迭代法分析 (3) 1.3计算结果分析及结论 (4) 问题二追赶法解三对角矩阵 (5) 2.1 问题 (5) 2.2 问题分析（追赶法） (6) 2.3 计算结果 (7) 问题三函数拟合 (7) 3.1 计算题目 (7) 3.2 题目分析 (7) 3.3 结果比较 (12) 问题四欧拉法解微分方程 (14) 4.1 计算题目 (14) 4.2.1 方程的准确解 (14) 4.2.2 Euler方法求解 (14) 4.2.3改进欧拉方法 (16) 问题五四阶龙格-库塔计算常微分方程初值问题 (17) 5.1 计算题目 (17) 5.2 四阶龙格-库塔方法分析 (18) 5.3 程序流程图 (18) 5.4 标准四阶Runge-Kutta法Matlab实现 (19) 5.5 计算结果及比较 (20) 问题六舍入误差观察 (22) 6.1 计算题目 (22) 6.2 计算结果 (22) 6.3 结论 (23) 7 总结 (24) 附录

研究生数值分析习题

1. 五个节点的Newton-Cotes 求积公式的代数精度为______，五个节点的求积公式最高代数精度为___________。（即Gauss 型求积公式） 2. 已知数值求积公式为3 11 ()[(1)4(2)(3)]3 f x dx f f f ≈++? ， 则其代数精度为______。 3. 数值积分公式1 '12 ()[(1)8(0)(1)]9 f x dx f f f -≈-++?的代数 精度为_________。 4. 要使求积公式1 110 1 ()(0)()4 f x dx f A f x ≈ +?具有2次代数精度，则1x =___,1A =___。 5. 在Newton-Cotes 求积公式：() ()()()n b n i i a i f x dx b a C f x =≈-∑? 中，当系数()n i C 是负值时，公式的稳定性不能保证，所以实际应用中，当___________时的Newton-Cotes 求积公式不能使用。 ()8()7()10()6A n B n C n D n ≥≥≥≥ 6. 若用复化梯形公式计算1 0x e dx ?，要求误差不超过6 10-，利 用余项公式估计，至少用______个求积节点。 7. 对于Gauss 型求积公式3 1 ()()()b k k a k f x x dx A f x ρ=≈∑?，其中 ()x ρ为权函数，下列说法错误的是_________。

（A ）该求积公式一定是稳定的； （B ）3 1()k k k A f x b a ==-∑； （C ）该求积公式的代数精度为5； （D ）2 (35)()()0b a x x x x dx ωρ-=? ，其中3 1 ()()k k x x x ω==∏-。 8. 0{()}k k x ?∞ =是区间[0,1]上权函数 ()x x ρ=的最高系数为1的正交多项式族，其中0()1x ?=，则1 40()_______x x dx ?=?。 9. 构造代数精度最高的如下形式的求积公式，并求出其代数精度： 1 010 1 ()()(1)2 xf x dx A f A f ≈+? 10. 数值积分公式形如 1 ()()(0)(1)(0)(1)xf x dx S x Af Bf Cf Df ''≈=+++? （1）试确定参数A 、B 、C 、D ，使公式的代数精度尽量高； （2）设4 ()[0,1]f x C ∈，推导余项公式1 0()()()R x xf x dx S x =-?， 并估计误差。 11. 用8n =的复化梯形公式和复化Simpson 公式计算 1 x e d x -? 时， （1）试用余项估计其误差； （2）计算积分的近似值。

研究生数值分析试卷

1 I(a,b) 2 ax 2 b x dx 2005~2006学年第一学期硕士研究生期末考试试题（A 卷） 科目名称：数值分析 学生所在院： ________ 学号： ________ 姓名： ______ 注意：所有的答题内容必须答在答题纸上，凡答在试题或草稿纸上的一律无效。 一、 （15分）设求方程12 3x 2cosx 0根的迭代法 / 2 X ki 4 cosx k 3 （1） 证明对X o R ,均有lim X k x *,其中X *为方程的根. k （2） 此迭代法收敛阶是多少？证明你的结论. 二、 （12分）讨论分别用Jacobi 迭代法和Gauss-Seidel 迭代法求解下列方程组的 收敛性。 x 1 2x 2 2x 3 1, X 1 X 2 X 3 1, 2x 1 2x 2 x 3 0. 0 0a 非病态的。（范数用HI ） 求f （X ）的Hermite 插值多项式H 3（x ），并给出截断误差R （x ） f （x ） H 3（x ） 五、（10分）在某个低温过程中，函数 y 依赖于温度x （T ）的试验数据为 已知经验公式的形式为 y ax bx 2，试用最小二乘法求出 a ， b 、（8分）若矩阵A 2a a 0 0 a 0，说明对任意实数a 0，方程组AX b 都是 四、（15六、（12分）确定常数 a ，b 的值，使积分

、（15分）设求方程 12 3x 2cosx 0根的迭代法 取得最小值。 七、（14分）已知Legendre 勒让德）正交多项式L n （x ）有递推关系式: L o (x) 1, L i (x) x (n 1, 2,) 试确定两点的咼斯一勒让德（G — L ）求积公式 1 1 f （x ）dx 入仁花）A 2f （x 2） 的求积系数和节点，并用此公式近似计算积分 1 2 一 e x dx 1 八、（14分）对于下面求解常微分方程初值冋题 dx f （x,y ）的单步法: y （x 。） y 。 1 1 y n 1 y n h(?k 1 - k 2) k 1 f(X n ,y n ) k 2 f(X n h, y n hkj （1） 验证它是二阶方法； （2） 确定此单步法的绝对稳定域。 2005~2006学年第一学期硕士研究生期末考试试题（B 卷） 科目名称：数值分析 学生所在院： _______ 学号： _________ 姓名： ______ 注意：所有的答题内容必须答在答题纸上，凡答在试题或草稿纸上的一律无效。 一、（12分）讨论分别用Jacobi 迭代法和Gauss-Seidel 迭代法求解下列方程组的 收敛性。 X 1 2x 2 2x 3 1, X 1 X 2 X 3 1, 2x 1 2x 2 x 3 0. L n 1(X ) 2n 1 n 1 xL n (x) L n 1（X ）

数值分析上机实验

目录 1 绪论 (1) 2 实验题目(一) (2) 2.1 题目要求 (2) 2.2 NEWTON插值多项式 (3) 2.3 数据分析 (4) 2.3.1 NEWTON插值多项式数据分析 (4) 2.3.2 NEWTON插值多项式数据分析 (6) 2.4 问答题 (6) 2.5 总结 (7) 3 实验题目(二) (8) 3.1 题目要求 (8) 3.2 高斯-塞德尔迭代法 (8) 3.3 高斯-塞德尔改进法—松弛法 (9) 3.4 松弛法的程序设计与分析 (9) 3.4.1 算法实现 (9) 3.4.2 运算结果 (9) 3.4.3 数据分析 (11) 4 实验题目(三) (13) 4.1 题目要求 (13) 4.2 RUNGE-KUTTA 4阶算法 (13) 4.3 RUNGE-KUTTA 4阶算法运算结果及数值分析 (14) 总结 (16) 附录A (17)

1绪论 数值分析是计算数学的一个主要部分，它主要研究各类数学问题的数值解法，以及分析所用数值解法在理论上的合理性。实际工程中的数学问题非常复杂，所以往往需要借助计算机进行计算。运用数值分析解决问题的过程：分析实际问题，构建数学模型，运用数值计算方法，进行程序设计，最后上机计算求出结果。 数值分析这门学科具有面向计算机、可靠的理论分析、好的计算复杂性、数值实验、对算法进行误差分析等特点。 本学期开设了数值分析课程，该课程讲授了数值分析绪论、非线性方程的求解、线性方程组的直接接法、线性方程组的迭代法、插值法、函数逼近与曲线拟合、数值积分和数值微分、常微分方程初值问题的数值解法等内容。其为我们解决实际数学问题提供了理论基础，同时我们也发现课程中很多问题的求解必须借助计算机运算，人工计算量太大甚至无法操作。所以学好数值分析的关键是要加强上机操作，即利用计算机程序语言实现数值分析的算法。本报告就是基于此目的完成的。 本上机实验是通过用计算机来解答数值分析问题的过程，所用的计算工具是比较成熟的数学软件MATLAB。MATLAB是Matrix Laboratory的缩写，是以矩阵为基础的交互式程序计算语言。MATLAB是一款具有强大的矩阵运算、数据处理和图形显示功能的软件，其输出结果可视化，编程效率极高，用极少的代码即可实现复杂的运行，因此它使工程技术人员摆脱了繁琐的程序代码，以便快速地验证自己的模型和算法。其主要特点包括：强大的数值运算功能；先进的资料视觉化功能高阶但简单的程序环境；开方及可延展的构架；丰富的程式工具箱。 在科学研究和工程计算领域经常会遇到一些非常复杂的计算问题，利用计算器或手工计算是相当困难或无法实现的，只能借助计算机编程来实现。MATLAB将高性能的数值计算和可视化的图形工具集成在一起，提供了大量的内置函数，使其在科学计算领域具有独特的优势。 最后感谢数值分析课程任课教师赵海良老师的悉心指导！