数值分析第一章思考题
数值计算方法思考题
数值计算方法思考题第一章 预篇1．什么是数值分析？它与数学科学和计算机的关系如何？2．何谓算法？如何判断数值算法的优劣？3．列出科学计算中误差的三个来源，并说出截断误差与舍入误差的区别。
4．什么是绝对误差与相对误差？什么是近似数的有效数字？它与绝对误差和相对误差有何关系？5．什么是算法的稳定性？如何判断算法稳定？为什么不稳定算法不能使用？6．判断如下命题是否正确：（1）一个问题的病态性如何，与求解它的算法有关系。
（2）无论问题是否病态，好的算法都会得到好的近似解。
（3）解对数据的微小变化高度敏感是病态的。
（4）高精度运算可以改善问题的病态性。
（5）用一个稳定的算法计算良态问题一定会得到好的近似值。
（6）用一个收敛的迭代法计算良态问题一定会得到好的近似值。
（7）两个相近数相减必然会使有效数字损失。
（8）计算机上将1000个数量级不同的数相加，不管次序如何结果都是一样的。
7．考虑二次代数方程的求解问题ax 2 + bx + c = 0.下面的公式是熟知的aac b b x 242-±-=. 与之等价地有ac b b c x 422--=.对于 a = 1, b = -100 000 000 , c = 1应当如何选择算法？8．指数函数有著名的级数展开++++=!3!2132x x x e x如果对x < 0用上述的级数近似计算指数函数的值，这样的算法结果是否会好？为什么？9．考虑数列x i , i = 1,…, n , 它的统计平均值定义为∑==n i i x x x 11 它的标准差2112)(11⎥⎦⎤⎢⎣⎡--=∑-n i i x x n σ 数学上它等价于2112211⎥⎥⎦⎤⎢⎢⎣⎡⎪⎪⎭⎫ ⎝⎛--=∑=n i i x n x n σ 作为标准差的两种算法，你如何评价它们的得与失？第二章 非线性方程求根1．判断如下命题是否正确：(a) 非线性方程的解通常不是唯一的；(b) Newton 法的收敛阶高于割线法；(c) 任何方法的收敛阶都不可能高于Newton 法；(d) Newton 法总是比割线法更节省计算时间；(e) 如果函数的导数难于计算，则应当考虑选择割线法；(f) Newton 法是有可能不收敛；(g) 考虑简单迭代法x k +1 = g (x k )，其中x * = g (x *)。
数值分析思考题1
数值分析思考题11、 讨论绝对误差（限）、相对误差（限）与有效数字之间的关系。
答：（1）绝对误差（限）与有效数字：若*120....10m n x ααα=⨯（a 1≠0，m 为整数） 绝对误差：*1*102m n e x x -=-≤⨯，那么*x 就有 n 个有效数字。
因此，从有效数字可以算出近似数的绝对误差限；有效数字位数越多，其绝对误差限也越小。
（2）相对误差限与有效数字：*120....10m n x ααα=⨯（a1≠0，m 为整数）相对误差限:*1111110*1210*102m n n r m x x e x αα--+-⨯-=≤=⨯⨯，*1*102m n e x x -=-≤⨯，11*10m x α-≥⨯可见*x 至少有n 位有效数字。
2、相对误差在什么情况下可以用下式代替？答：实际情况下真实值 x 是无法得到的，当测量值与真实值之间的误差可以忽略不计时，可用下式代替。
3、 查阅何谓问题的“病态性”，并区分与“数值稳定性”的不同点。
r e x x e x x *****-==答：病态性：数学问题本身性质所决定的，与算法无关，却能引起问题真解很大变化。
同:都是输入数据的微小误差导致输出数据误差的增大。
异：数值稳定性是相对于算法而言的，算法的不同直接影响结果的不同；而病态性是数学模型本身的问题，与算法无关。
4、 取，计算，下列方法中哪种最好？为什么？（1）(33-，（2）(27-，（3）()313+，（4）()611，（5）99-答：)631 5.05110-≈⨯ （1）(()333332 1.41 5.83210--≈-⨯≈⨯（2）223(7(75 1.41) 2.510--≈-⨯=⨯（3331 5.07310(32 1.41)-≈≈⨯+⨯（4361 5.10410(1.411)-≈≈⨯+（5）9999700.3-≈-=方法3最好，误差最小141.≈)61。
《数值分析》第一章答案
《数值分析》第⼀章答案习题11．以下各表⽰的近似数，问具有⼏位有效数字？并将它舍⼊成有效数。
（1）*1x ＝451.023， 1x ＝451.01；（2）*2x ＝－0.045 113， 2x ＝－0.045 18；（3）*3x ＝23.421 3， 3x ＝23.460 4；（4）*4x ＝31， 4x ＝0.333 3；（5）*5x ＝23.496， 5x ＝23.494；（6）*6x ＝96×510， 6x ＝96.1×510；（7）*7x ＝0.000 96， 7x ＝0.96×310-；（8）*8x ＝－8 700， 8x ＝－8 700.3。
解：(1) =*1x 451.023 =1x 451.01=-1*1x x 0.01311021-?≤，1x 具有4位有效数字。
→1x 451.0(2) -=*2x 0.045 113 -=2x 0.045 18=-241021x x 0.045 18045113.0-=0.000 06731021-?<2x 具有2位有效数字，045.02-→x(3)=*3x x =-4604.234213.23=-4213.234604.231 10210391.0-?≤3x 具有3位有效数字，4.233→x (不能写为23.5) (4) =*4x 31 ，=4x 0.3333=-4*4x x 41021000033.0-?<,4x 具有4位有效数字，=4x 0.3333(5) =*5x 23.496，=5x 23.494=-5*5x x =-494.23496.2321021002.0-?<5x具有4位有效数字，→5x 23.50 (不能写为23.49)(6) =*6x 51096?710961.0?==-6*6x x 710001.0-?72101021--??≤6x 具有2位有效数字，57610961096.0?=?=x(7) =*7x 0.00096 371096.0-?=x3*71096.0-?=x =-7*7x x 0 7x 精确(8) 8700*8-=x 8x 3.8700-=8*8x x -010213.0?≤=8x 具有4位有效数字，8x 8700-=精确2.以下各数均为有效数字： (1) 0.1062 + 0.947; (3)2.747?6.83; (2)23.46―12.753; (4)1.473 / 0.064 。
数值分析第一章思考题
数值分析第一章思考题第一章思考题(2012级本科学生作品)1、什么样的算法被称为不稳定算法？试列举一个例子进行说明。
在算法执行过程中，舍入算法对计算结果影响大的一类算法被称为数值不稳定的一种算法。
例如，假设初始数据有一点微小误差，就会对一个算法的数据结构产生很大的影响，造成误差扩散。
用计算公式ln 1ln n n =-，构造出的递推算法是一个数值不稳定的算法；而另一公式ln 1(1ln)/n -=-则可以构造出一个数值稳定的算法。
2、我们都知道秦九韶算法能够减少运算次数，高中也学过他的具体过程，请举出一个例子并用秦九韶算法计算。
答；一般的，一元n 次多项式的求值需要经过(1)/2n n +次乘法和n 次加法，而秦九韶算法只需要n 次乘法和n 次加法。
具体的不太会了。
3、为什么要设立相对误差的概念？答：相对误差是近似值误差与精确值的比值，用来衡量近似值的近似程度。
x=10±1，y=1000±5。
虽然x 的误差比y 的误差小，但y 的近似程度比x 更好。
这单用误差无法表现出来，而相对误差可以解决这个问题。
4、误差在生活中有什么作用？答：误差的作用不仅仅体现在数学课题研究中，在生活中误差的作用也非常大，比如在建筑行业中，设计图纸时必须要达到一定的精确度才行。
5、有效数字以及计算规则答：有效数字是指实际上能测量到的数值，在该数值中只有最后一位是可疑数字，其余的均为可靠数字。
它的实际意义在于有效数字能反映出测量时的准确程度。
例如，用最小刻度为0.1cm 的直尺量出某物体的长度为11.23cm ，显然这个数值的前3位数是准确的，而最后一位数字就不是那么可靠，医|学教育网搜集整理因为它是测试者估计出来的，这个物体的长度可能是11.24cm ，亦可能是11.22cm ，测量的结果有±0.01cm 的误差。
我们把这个数值的前面3位可靠数字和最后一位可疑数字称为有效数字。
这个数值就是四位有效数字。
数值分析思考题
数值分析重点考察内容第一章：基本概念第二章：Gauss消去法，Lu分解法第三章：题型：具体题＋证明，误差分析三个主要迭代法，条件误差估计，范数的小证明第四章：掌握三种插值方法：拉格朗日，牛顿，厄尔米特，误差简单证明，构造复合函数第五章：最小二乘法计算第六章：梯形公式，辛普森（抛物线）公式，高斯公式三个重要公式，误差分析。
高斯求积公式的构造第七章：几种常用的迭代格式构造，收敛性证明。
第九章：基本概念（收敛阶，收敛条件，收敛区域等）, 简单欧拉法。
第一章 误差1. 科学计算中的误差来源有4个，分别是________，________，________，________。
2. 用Taylor 展开近似计算函数000()()'()()f x f x f x x x ≈+-，这里产生是什么误差？3. 0.7499作34的近似值，是______位有效数字，65.380是舍入得到的近似值，有____几位有效数字，相对误差限为_______. 0.0032581是四舍五入得到的近似值，有_______位有效数字.4. 改变下列表达式，使计算结果比较精确：（1）11,||1121x x x x --++ （2）||1x （3） 1cos ,0,|| 1.x x x x -≠ （4） sin sin ,αβαβ-≈5.采用下列各式计算61)时，哪个计算效果最好？并说明理由。
（1）（2）99-（3）6(3- （46. 已知近似数*x 有4位有效数字，求其相对误差限。
上机实验题：1、利用Taylor 展开公式计算 0!kx k x e k ∞==∑，编一段小程序，上机用单精度计算x e 的函数值. 分别取 x =1，5，10，20，-1，-5，-10，-15，-20，观察所得结果是否合理，如不合理请分析原因并给出解决方法.2、已知定积分10,0,1,2,,206n n 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. 在下列各对数中，X 是精确值ａ的近似值（1） ａ=π，x=3.1 （2） ａ=1/7，x=0.143 （3） ａ=π/1000，x=0.0031 （4） ａ=100/7，x=14.3 试估计x 的绝对误差和相对误差。
解：（1） e=∣3.1-π∣≈0.0416， δr = e/∣x ∣≈0.0143 （2） e=∣0.143-1/7∣≈0.0143 δr = e/∣x ∣≈0.1 （3） e=∣0.0031-π/1000∣≈0.0279 δr = e/∣x ∣≈0.9 (4) e=∣14.3-100/7∣≈0.0143 δr = e/∣x ∣≈0.0012. 已知四个数：x 1=26.3，x 2=0.0250， x 3= 134.25，x 4=0.001。
试估计各近似数的有效位数和误差限，并估计运算μ1= x 1 x 2 x 3和μ1= x 3 x 4 /x 1的相对误差限。
解：x 1=26.3 ｎ=3 δx 1=0.05 δr x 1=δx 1/∣x 1∣=0.19011×10-2x 2=0.0250 ｎ=3 δx 2=0.00005 δr x 2=δx 2/∣x 2∣=0.2×10-2x 3= 134.25 ｎ=5 δx 3=0.005 δr x 3=δx 3/∣x 3∣=0.372×10-4x 4=0.001 ｎ=1 δx 4=0.0005 δr x 4=δx 4/∣x 4∣=0.5由公式：e r （μ）= e （μ）/∣μ∣≦1/∣μ∣Σni=1∣∂f/∂x i ∣δx ie r （μ1）≦1/∣μ1∣［x 2 x 3δx 1+ x 1 x 3δx 2 +x 1x 2δx 3］ =0.34468/88.269275 =0.0039049e r （μ2）≦1/∣μ2∣［-x 3 x 4/ x 21δx 1+ x 4/ x 1δx 3 + x 3/ x 1δx 4］ =0.497073. 设精确数ａ>0，x 是ａ的近似值，x 的相对误差限是0.2，求㏑x 的相对误差限。
数值分析课后习题与解答
课后习题解答第一章绪论习题一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),使它满足解：这种题目可以有很多方法去做，但应以简单为宜。
数值分析作业思考题
数值分析思考题11、讨论绝对误差（限）、相对误差（限）与有效数字之间的关系。
2、相对误差在什么情况下可以用下式代替？3、查阅何谓问题的“病态性”，并区分与“数值稳定性”的不同点。
4、取，计算，下列方法中哪种最好？为什么？（1）(33-，（2）(27-，（3）()313+，（4）()611，（5）99-数值实验数值实验综述：线性代数方程组的解法是一切科学计算的基础与核心问题。
求解方法大致可分为直接法和迭代法两大类。
直接法——指在没有舍入误差的情况下经过有限次运算可求得方程组的精确解的方法，因此也称为精确法。
当系数矩阵是方的、稠密的、无任何特殊结构的中小规模线性方程组时，Gauss消去法是目前最基本和常用的方法。
如若系数矩阵具有某种特殊形式，则为了尽可能地减少计算量与存储量，需采用其他专门的方法来求解。
Gauss消去等同于矩阵的三角分解，但它存在潜在的不稳定性，故需要选主元素。
对正定对称矩阵，采用平方根方法无需选主元。
方程组的性态与方程组的条件数有关，对于病态的方程组必须采用特殊的方法进行求解。
数值计算方法上机题目11、实验1. 病态问题实验目的：算法有“优”与“劣”之分，问题也有“好”和“坏”之别。
所谓坏问题就是问题本身的解对数据变化的比较敏感，反之属于好问题。
希望读者通过本实验对此有一个初步的体会。
数值分析的大部分研究课题中，如线性代数方程组、矩阵特征值问题、非线性方程及方程组等都存在病态的问题。
病态问题要通过研究和构造特殊的算法来解决，当然一般要付出一些代价（如耗用更多的机器时间、占用更多的存储空间等）。
问题提出：考虑一个高次的代数多项式re x xex x*****-==141.≈)61∏=-=---=201)()20)...(2)(1()(k k x x x x x p （E1-1）显然该多项式的全部根为l ，2，…，20，共计20个，且每个根都是单重的（也称为简单的）。
现考虑该多项式方程的一个扰动0)(19=+xx p ε (E1-2)其中ε是一个非常小的数。
(完整版)数值分析第四版习题和答案解析
h 应取多少 ?
9. 若 yn 2 n , 求 4 yn 及 4 yn .
10. 如 果 f ( x) 是 m 次 多 项 式 , 记 f (x) f (x h) f ( x) , 证 明 f (x) 的 k 阶 差 分
k f (x)(0 k m) 是 m k 次多项式 , 并且 m l f ( x) 0(l 为正整数 ).
.专业资料 . 整理分享 .
.WORD. 格式 .
11. 证明 ( f k g k ) fk g k gk 1 f k .
n1
fk gk
12. 证明 k 0
fngn
f0 g0
n1
gk 1 f k .
k0
n1
2 yj
13. 证明 j 0
14. 若 f (x) a0
yn y0. a1 x L an 1 xn 1
.
.专业资料 . 整理分享 .
.WORD. 格式 .
18. f ( x) 、 g( x) C1 a,b , 定义
b
b
( a)( f , g) f (x) g (x)dx;( b)( f , g ) f ( x) g ( x) dx f (a) g (a);
a
a
问它们是否构成内积 ?
6
1 x dx
19. 用许瓦兹不等式 (4.5) 估计 0 1 x 的上界 , 并用积分中值定理估计同一积分的上下界
5. 计算球体积要使相对误差限为 1％ , 问度量半径 R时允许的相对误差限是多少 ?
6. 设 Y0 28, 按递推公式
1
Yn Yn 1
783
100
( n=1,2, … )
计算到 Y100 . 若取 783 ≈ 27.982( 五位有效数字 ), 试问计算 Y100 将有多大误差 ?
第一章复习与思考题
第一章复习与思考题1. 什么是数值分析？它与数学科学和计算机的关系如何？答：数值分析也称计算数学，是数学科学的一个分支，主要研究的是用计算机求解各种数学问题的数值计算方法及其理论与软件实现.数值分析以数学问题为研究对象，但它并不像纯数学那样只研究数学本身的理论，而是把理论与计算紧密结合，着重研究数学问题的数值方法及其理论.2. 何谓算法？如何判断数值算法的优劣？答：一个数值问题的算法是指按规定顺序执行一个或多个完整的进程，通过算法将输入元变换成输出元.一个面向计算机，有可靠理论分析且计算复杂性好的算法就是一个好算法. 因此判断一个算法的优劣应从算法的可靠性、准确性、时间复杂性和空间复杂性几个方面考虑.3. 列出科学计算中误差的三个来源，并说出截断误差与舍入误差的区别.答：用计算机解决实际问题首先要建立数学模型，它是对被描述的实际问题进行抽象、简化而得到的，因而是近似的，数学模型与实际问题之间出现的误差叫做模型误差.在数学模型中往往还有一些根据观测得到的物理量，如温度、长度等，这些参量显然也包含误差，这种由观测产生的误差称为观测误差.当数学模型不能得到精确解时，通常要用数值方法求它的近似解，其近似解和精确解之间的误差称为截断误差或方法误差.有了求解数学问题的计算公式以后，用计算机做数值计算时，由于计算机字长有限，原始数据在计算机上表示时会产生误差，计算过程又可能产生新的误差，这种误差称为舍入误差.截断误差和舍入误差是两个不同的概念，截断误差是由所采用的数值方法而产生的，因而也称方法误差，舍入误差是由数值计算而产生的.4. 什么是绝对误差与相对误差？什么是近似数的有效数字？它与绝对误差和相对误差有何关系？答：设x 为准确值，*x 为x 的一个近似值，称x x e -=**为近似值*x 的绝对误差，简称误差. 近似值的误差*e 与准确值x 的比值xx x x e -=**称为近似值*x 的相对误差，记作*r e .通常我们无法知道误差的准确值，只能根据测量工具或计算情况估计出误差绝对值的一个上界*ε，*ε叫做近似值的误差限.若近似值*x 的误差限是某一位的半个单位，该位到*x 的第一位非零数字共有n 位，就说*x 有n 位有效数字.有效数位越多，绝对误差限越小，相对误差限也越小.5. 什么是算法的稳定性？如何判断算法稳定？为什么不稳定算法不能使用？答：一个算法如果输入数据有误差，而在计算中舍入误差不增长，则称此算法是数值稳定的；否则称为不稳定的.判断一个算法是否稳定主要是看初始数据误差在计算中的传播速度，如果传播速度很快就是数值不稳定的.对于不稳定的算法来说，由于其误差传播是逐步扩大的，因而计算结果不可靠，所以不稳定的算法是不能使用的.6. 什么是问题的病态性？它是否受所用算法的影响?答：对一个数值问题本身来说，如果输入数据有微小扰动（即误差），引起输出数据（即问题解）相对误差很大，这就是病态问题.病态问题是数值问题本身固有的，不是由计算方法引起的，病态性并不受所用算法的影响，对病态问题必须采用特殊的方法以减少误差危害.7. 什么是迭代法？试利用03=-a x 构造计算3a 的迭代公式.答：迭代法是一种按同一公式重复计算逐次逼近真值的算法，是数值计算普遍使用的重要方法. 在计算3a 时，可从等价的方程求根问题03=-a x 出发，利用方程的等价形式)2(312xa x x +=即可得到计算3a 的迭代公式)2(3121k k k x a x x +=+. 8. 直接利用以直代曲的原则构造求方程02=-a x 的根a x =*的迭代法.答：求方程0)(=x f 的根在几何上就是求曲线)(x f y =与x 轴交点*x 的横坐标. 假如已给出一个近似值k x ，用该点))(,(k k x f x 处的切线逼近曲线，令1+k x 为该切线与x 轴交点的横坐标，一般情况下，1+k x 近似方程根*x 的程度比k x 近似*x 的程度要好，这就是以直代曲的思想. 曲线a x y -=2在点))(,(k k x f x 处的切线方程为a x x x y k k --=22，切线方程的根)(21kk x a x x +=，以此作为新的近似值，就得到了求方程02=-a x 的根a x =*的迭代法)(211k k k x a x x +=+. 9. 举例说明什么是松弛技术. 答：在积分近似计算的梯形公式∑=-+=ni i i n x f x f h T 11)]()([2中，取2,1=n 可分别得 .2)],()(2)([4)],()([221b a c b f c f a f a b T b f a f a b T +=++-=+-=令121221)1()(T T T T T S ωωω-+=-+=,若取3/1=ω，则得)],()(4)([63134121b f c f a f a b T T S ++-=-= 这就是松弛技术，ω称为松弛因子.10. 考虑无穷级数∑∞=11n n,它是发散的，在计算机上计算它的部分和，会得到什么结果？为什么？答：虽然在理论上无穷级数∑∞=11n n是发散的，但在计算机上计算时，由于计算机只能进行有限数的计算，所以无论n取多大的值，级数的和都是有限数.11. 判断下列命题的正确性：(1) 解对数据的微小变化高度敏感是病态的.(2) 高精度运算可以改善问题的病态性.(3) 无论问题是否病态，只要算法稳定都能得到好的近似值.(4) 用一个稳定的算法计算良态问题一定会得到好的近似值.(5) 用一个收敛的迭代法计算良态问题一定会得到好的近似值.(6) 两个相近数相减必然会使有效数字损失.(7) 计算机上将1000个数量级不同的数相加，不管次序如何结果都是一样的.答：(1)对.(2)错.(3)错.(4)错.(5)对.(6)错.(7)错.。
数值分析思考题
数值分析思考题61、数值计算中迭代法与直接法的区别是什么(D直接法是指在没有舍入误差的情况下经过有限次运算求得方程组的精确解的方法。
直接法又称为精确法。
(2)迭代法是采取逐次逼近的方法，即从一个初始向量出发，按照一定的计算格式，构造一个向量的无穷序列，其极限才是方程组的精确解，只经过有限次运算得不得精确解。
迭代法是一种逐次逼近的方法，与直接法比较，具有程序简单，存储量小的优点。
2、详述你所知道的线性方程组的迭代法的收敛性定理。
迭代公式X(_I)二Bx(k)+ g(k二0,1,2, ?)收敛的充分必要条件是M k->0.假设矩阵M的谱半径p(B),可知MkTO的充分必要条件是p (B) < 1 o 迭代公式x(k I)二Bx(k)+ g(k 二0,1,2, ?) 和x(k + 1)二Bix(k + 1) + B2x(k) + g(k 二0,1,2,?),收敛。
严格对角占优线性方程组Ax二b(其中A e R m x n,b e L)的Jacobi 迭代公式x(k + 1) = Bx(k)+ g(k = 0,1,2,?),收敛。
Gauss-Seidel 迭代公式x(k + 1)二Bix(k + 1) + B2x(k) + g(k = 0,1,2,?),收敛。
3、详述你所知道的非线性方程(组)的迭代法以及收敛性结果。
(1)不动点迭代法：不一定收敛，若存在常数L<1 ,使得I 4> (x) - 0 (y) I W L|x 一y|,?x, y G [a, b],则收敛于x*。
(2)斯蒂芬森迭代法：若不动点迭代公式的迭代函数e(x)在不动点X*的某邻域内具有二阶连续导数，e'(x*)二A工1且A工0,则二阶收敛，极限是X*。
(3)牛顿迭代法：收敛4、举例说明解线性方程组的S0R方法的最佳松弛因子与何种因素有解线性方程组的S0R方法的最佳松弛因子与迭代矩阵的谱半径有关，是单峰关系。
数值分析思考题
数值分析复习思考题（2006-12-28）这几天的答疑时间中，解答了部分同学的问题，更多是作为教师的深入思考。
而共同探讨问题是非常重要的。
由于时间有限，这个文档中提出问题的深度可能不够，有些问题还没给出解答，希望研究生同学一起来思考，提出更多的问题。
我会在以后的时间中形成新的文档。
第一章 思考题1．在科学计算中，一般认为误差的来源有几种？列举在数值分析课中主要讨论误差。
数值计算中一个基本的手段是近似，所以就有了各种误差。
误差来源有四种：模型误差，观测误差，截断误差，舍入误差。
一般分为两类，第一类是固有误差（包括模型误差和观测误差），第二类是计算误差（包括截断误差和舍入误差）。
计算方法课中主要讨论计算误差。
这是因为在用计算机解决数学问题时，常常用“有限代替无穷，用近似代替准确”。
例如，解决连续性问题时通常要将其转化为离散问题求解，这将引起截断（方法）误差；由于机器数的位数有限，计算机表示数据时一般带有舍入误差。
下面不全面列举出本课程内容涉及的误差线性方程组直接求解方法——舍入误差多项式插值方法——插值误差数据拟合方法——残差数值积分方法——求积误差微分方程数值解方法——局部截断误差………………………………………………2．有效数字的概念是如何抽象而来的，请简单给予叙述。
有效数字位数与计算近似值x的误差这两个概念是通过末位数半个单位相联系的。
由于计算机的机器数只能表示有限位浮点数，对于很多数据只能近似表示，近似采用“四舍五入”的原则进行。
有效数字概念正是根据日常生活中的“四舍五入”原则抽象而来的。
若近似值x的绝对误差限是某一位上半个单位，该位到x的第一位非零数字一共有n位，则称这一近似数具有n位有效数字。
而相对误差则与有效数位数基本一致。
3．什么样的算法被称为是不稳定的算法？试举一个例子说明在算法执行过程中，舍入误差对计算结果影响大的一类算法被称为数值不稳定的算法。
例如初始数有一点微小的误差，就会对一个算法的数据结果产生较大的影响，造成误差扩散，用计算公式I n = 1 – n I n-1构造出的递推算法是一个数值不稳定的算法；而另一个公式I n-1= ( 1 – I n )/n则可以构造出一个数值稳定的算法。
数值分析思考题答案
数值分析思考题答案数值分析课程思考题1．叙述拉格朗⽇插值法的设计思想。
Lagrange插值是把函数y=f(x)⽤代数多项式pn(x)代替，构造出⼀组n次差值基函数；将待求得n次多项式插值函数pn(x)改写成另⼀种表⽰⽅式，再利⽤插值条件确定其中的待定函数，从⽽求出插值多项式。
2．函数插值问题的提出以及插值法发展的脉络。
问题的提出：实际问题中常遇到这样的函数y=f(x),其在某个区间[a,b]上是存在的。
但是，通过观察或测量或试验只能得到在[a,b]区间上有限个离散点x0,x1,…,xn上的函数值y=f(xi),(i=0,…,n)或者f(x)函数表达式是已知的，但却很复杂⽽不便于计算希望⽤⼀个简单的函数描述它。
发展脉络：在⼯程中⽤的多的是多项式插值和分段多项式插值。
在多项式插值中，⾸先谈到的是Lagrange插值，其成功地⽤构造插值基函数的⽅法解决了求n次多项式插值函数的问题，但是其⾼次插值基函数计算复杂，且次数增加后，插值多项式需要重新计算，所以在此基础上提出Newton插值，它是另⼀种构造插值多项式的⽅法，与Lagrange插值相⽐，具有承袭性和易于变动节点的特点。
如果对插值函数，不仅要求他在节点处与函数同值，还要求它与函数有相同的⼀阶，⼆阶甚⾄更⾼阶的导数值，这就提出了Hermite插值，它是利⽤未知函数f(x)在插值节点上的函数值及导数值来构造插值多项式的。
为了提⾼精度，加密节点时把节点分成若⼲段，分段⽤低次多项式近似函数，由此提出了分段多项式插值。
最后，由于许多⼯程中对插值函数的光滑性有较⾼的要求，就产⽣了样条插值。
3．描述数值积分算法发展和完善的脉络。
数值积分主要采⽤插值多项式来代替函数构造插值型求积公式。
通常采⽤Lagrange插值。
如果取等距节点，则得到Newton-Cotes公式，其中，当n=1时，得到梯形公式；当n=2时，得到Simpson公式；当n=4时，得到Cotes公式。
常州大学数值分析第一章习题解答
1.1解：m=3;f=@(x)digit(digit(x^4,m)- digit(x^3,m)+ digit(3*x^2,m)+ digit(x-2,m),m);g=@(x)digit(digit(digit( digit(digit(digit( (x-1)*x,m)+3,m)*x,m)+1,m)*x,m)-2,m);f(3.33)g(3.33)有ans = 121ans =121实际上，当m=2时，就可以看出这两种算法在计算的精确度上的区别：m=2;f=@(x)digit(digit(x^4,m)- digit(x^3,m)+ digit(3*x^2,m)+ digit(x-2,m),m);g=@(x)digit(digit(digit( digit(digit(digit( (x-1)*x,m)+3,m)*x,m)+1,m)*x,m)-2,m);f(3.33)g(3.33)有ans = 120ans =130，可以看出，两者在计算精度上的不同区别，数学上恒等，在数值上不一定恒等。
1.2解：（1）精确到小数点后第三位，故有4位有效数字（2）精确到小数点后第三位，故有2位有效数字（3）精确到小数点后第三位，故有0位有效数字1.3 解;记圆的面积为S，由题意有|e(S)|≤1%。
由S=πr2知：dS=2πrdr所以dS/S=(2πrdr)/(πr2)=2(dr/r)∴|e(r)|≈1／2|e(S)|≤0.5×1%=0.5%1.4 解：由题有：|e(x)|≤1/2×10^-2 ; |e(y)|≤1/2×10^-2; |e(z)||≤1/2×10^-2∴|e（S）|≈|xe(x)+ye(y)|+ |ze(z)|^2≈x|e(x)|+y|e(y)|+z^2|z(z)|^2≤4.21×0.005+1.79×1.005+2.11×2.11×0.005^2=0.03≤1/2×10^-1又S=4.21*1.79+2.11^2=11.988∴S至少具有3位有效数字。
常州大学数值分析第一章习题解答
1.1解：m=3;f=@(x)digit(digit(x^4,m)- digit(x^3,m)+ digit(3*x^2,m)+ digit(x-2,m),m);g=@(x)digit(digit(digit( digit(digit(digit( (x-1)*x,m)+3,m)*x,m)+1,m)*x,m)-2,m);f(3.33)g(3.33)有ans = 121ans =121实际上，当m=2时，就可以看出这两种算法在计算的精确度上的区别：m=2;f=@(x)digit(digit(x^4,m)- digit(x^3,m)+ digit(3*x^2,m)+ digit(x-2,m),m);g=@(x)digit(digit(digit( digit(digit(digit( (x-1)*x,m)+3,m)*x,m)+1,m)*x,m)-2,m);f(3.33)g(3.33)有ans = 120ans =130，可以看出，两者在计算精度上的不同区别，数学上恒等，在数值上不一定恒等。
1.2解：（1）精确到小数点后第三位，故有4位有效数字（2）精确到小数点后第三位，故有2位有效数字（3）精确到小数点后第三位，故有0位有效数字1.3 解;记圆的面积为S，由题意有|e(S)|≤1%。
由S=πr2知：dS=2πrdr所以dS/S=(2πrdr)/(πr2)=2(dr/r)∴|e(r)|≈1／2|e(S)|≤0.5×1%=0.5%1.4 解：由题有：|e(x)|≤1/2×10^-2 ; |e(y)|≤1/2×10^-2; |e(z)||≤1/2×10^-2∴|e（S）|≈|xe(x)+ye(y)|+ |ze(z)|^2≈x|e(x)|+y|e(y)|+z^2|z(z)|^2≤4.21×0.005+1.79×1.005+2.11×2.11×0.005^2=0.03≤1/2×10^-1又S=4.21*1.79+2.11^2=11.988∴S至少具有3位有效数字。
数值分析第五版第1章习题答案-李庆扬
第1章复习与思考题习题0，依据定义|X*X|/X*|X*-X|=X***X X：按定义求解E(lnX)=|lnX-LnX*|=ln(X/X*)=ln(1-)X|lnX ln*|/lnX*ln(1)/lnX*ln(1)/lnX*按泰勒展开求解，2(x*)=(x*)(x x*)"()(x x*)f f f1(lnX*)|lnX lnX*|lnX*(X X*)**r(lnX*)|lnX lnX*|/ln */ln *X X X X问题是解法1错了吗？ 很小时，ln(1-)=，求n X 的相对误差*0.02*X X ： 按照定义：{(10.02)X*}*(1.02)n n X (X )(X X*)/X* 1.021(10.02)1n n n n n多项式展开，有100.02n n i i：按照泰勒展开11(X X*)*(X X*)nX*(0.02X*)0.02X*)(XX*)/X*0.02*/*0.02n n n n n nn nn nnX nX X n问题是解法1错了吗？10.02n n i i 收敛于（0.02/(1-0.02)= 0.002004008016032064128256释？应该没有错，按照泰勒展开，相当于将误差限放大了。
3、下列各数都是经过四舍五入得到的近似数，即误差限不超过最后一位的半个单位，试指出它们是几位有效数字* 1.1021x ，0.031x ，385.6，56.430，7 1.0。
解：有4位有效数字；X2有1位有效数字；X3有3位有效数字；X4有1位有效数字**24x x ，*1x x 其中***123x x x ，，，公式（2.3）：1(A*)|()*|(X )n k k kfX 解：******124124()()()()0.510x x x x x x ************1232311321234313()()()()0.031385.60.510 1.1021385.60.510 1.10210.0310.510(0.59768212.48488 1.708255)100.214790815x x x x x x x x x x x x******2*2442244323335(/)1/()/()()1/56.4300.5100.031/(56.430)0.510(10.031/56.430)/56.4300.510=(0.99945064681906787169945064681907)0.5/56.430100.885610x x x x x x x5、计算求体积要使相对误差限为1%，问度量半径R 所允许的相对误差是多少？解：球的体积公式3V R23(V)3R (R)/R 3(R)/R 0.01r有(R)0.01/3R (R)(R)/1/300r R6、设Y0=28，按递推公式11783100n nY Y ，1,2,3.....n 计算到Y100，若取27.982（解：n n n (Y )(Y Y )，所以(Y )0n 。
数值分析思考题1
（1） ，（2） ，（3） ，（4） ，（；
（2） =0.52=0.25；
（3） =0.0050726；
（4） =0.00510385；
（5） =99-98.70=0.3；
2、相对误差在什么情况下可以用下式代替？
答：在实际计算时，由于真值常常是未知的，当 较小时，通常用 代替。
3、查阅何谓问题的“病态性”，并区分与“数值稳定性”的不同点。
答：（1）病态问题：对于数学问题本身，如果输入数据有微小变化，就会引起输出数据（即问题真解）的很大变化，这就是病态问题。
（2）不同点：数值稳定性是相对于算法而言的，算法的不同直接影响结果的不同；而病态性是数学问题本身性质所决定的，与算法无关，也就是说对病态问题，用任何算法（或方法）直接计算都将产生不稳定性。
数值分析思考题1
1、讨论绝对误差（限）、相对误差（限）与有效数字之间的关系。
答：（1）绝对误差（限）与有效数字：将x的近似值x*表示成
x*=±10m×(a1×10﹣1+a2×10﹣2+…an×10﹣n+…+ak×10﹣k+…)，其中m是整数，a1≠0，a1，a2，…，ak是0到9中的一个数字。若绝对误差 ，那么x*至少有n个有效数字，即a1，a2，…，an为有效数字，而an+1，…，ak，…不一定是有效数字。因此，从有效数字可以算出近似数的绝对误差限；有效数字位数越多，其绝对误差限也越小。
（2）相对误差（限）与有效数字：将x的近似值x*表示成
x*=±10m×(a1×10﹣1+a2×10﹣2+…an×10﹣n+…+ak×10﹣k+…)，其中m是整数，a1≠0，a1，a2，…，ak是0到9中的一个数字。若ak是有效数字，那么相对误差不超过 ;反之，如果已知相对误差r，且有 ，那么ak必为有效数字。
第一章习题解答数值分析
第一章习题解答3、设精准数ａ>0，x 是ａ的近似值，x 的相对误差限是，求㏑x 的相对误差限。
解：设=()u f x ，()()()()()()||||||||||()||()||||()||()||||r r r x e u df x e x df x e x e u u dx u dx u x df x x df x x e x x dx u dx u δ=≈==≤()||10.2(())||()||ln ln ln r r r r df x x x x f x x x dx u x x x xδδδδ==⋅⋅==4、长方体的长宽高别离为50cm ，20cm 和10cm ，试求测量误差知足什么条件时其表面积的误差不超过1cm 2。
解：设2()S xy yz zx =++{}[]{}(,,)(,,)(,,)()||()||()||()(,,)(,,)(,,)||||||max (),(),()2()2()2()max (),(),()1S x y z S x y z S x y z e S e x e y e z x y zS x y z S x y z S x y z e x e y e z x y z y z z x x y e x e y e z ∂∂∂≤++∂∂∂⎛⎫∂∂∂≤++ ⎪∂∂∂⎝⎭=+++++< {}[]11max (),(),()2()2()2()4()110.0031254(502010)320e x e y e z y z z x x y x y z <=+++++++===++测量误差小于时其表面积的误差不超过1cm 2。
7、计算61)1.414≈。
利用以下四种计算格式，试问哪一种算法误差最小。
（1（2）3(3- （3（4）99- 解：计算各项的条件数'()(())||()xf x cond f x f x = 11 1.41461(),(())| 3.5147(1)x f x cond f x x ===+ 322 1.414()(32),(())|49.3256x f x x cond f x ==-=33 1.41431(),(())| 1.4557(32)x f x cond f x x ===+ 44 1.414()9970,(())|4949x f x x cond f x ==-= 由计算知，第三种算法误差最小。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第一章思考题
(2012级本科学生作品)
1、什么样的算法被称为不稳定算法？试列举一个例子进行说明。
在算法执行过程中，舍入算法对计算结果影响大的一类算法被称为数值不稳定的一种算法。
例如，假设初始数据有一点微小误差，就会对一个算法的数据结构产生很大的影响，造成误差扩散。
用计算公式ln 1ln n n =-，构造出的递推算法是一个数值不稳定的算法；而另一公式ln 1(1ln)/n -=-则可以构造出一个数值稳定的算法。
2、我们都知道秦九韶算法能够减少运算次数，高中也学过他的具体过程，请举出一个例子并用秦九韶算法计算。
答；一般的，一元n 次多项式的求值需要经过(1)/2n n +次乘法和n 次加法，而秦九韶算法只需要n 次乘法和n 次加法。
具体的不太会了。
3、为什么要设立相对误差的概念？
答：相对误差是近似值误差与精确值的比值，用来衡量近似值的近似程度。
x=10±1，y=1000±5。
虽然x 的误差比y 的误差小，但y 的近似程度比x 更好。
这单用误差无法表现出来，而相对误差可以解决这个问题。
4、误差在生活中有什么作用？
答：误差的作用不仅仅体现在数学课题研究中，在生活中误差的作用也非常大，比如在建筑行业中，设计图纸时必须要达到一定的精确度才行。
5、有效数字以及计算规则
答：有效数字是指实际上能测量到的数值，在该数值中只有最后一位是可疑数字，其余的均为可靠数字。
它的实际意义在于有效数字能反映出测量时的准确程度。
例如，用最小刻度为0.1cm 的直尺量出某物体的长度为11.23cm ，显然这个数值的前3位数是准确的，而最后一位数字就不是那么可靠，医|学教育网搜集整理因为它是测试者估计出来的，这个物体的长度可能是11.24cm ，亦可能是11.22cm ，测量的结果有±0.01cm 的误差。
我们把这个数值的前面3位可靠数字和最后一位可疑数字称为有效数字。
这个数值就是四位有效数字。
在确定有效数字位数时，特别需要指出的是数字“0”来表示实际测量结果时，它便是
有效数字。
例如，分析天平称得的物体质量为7.1560g，滴定时滴定管读数为20.05mL，这两个数值中的“0”都是有效数字。
在0.006g中的“0”只起到定位作用，不是有效数字。
在计算中常会遇到下列两种情况：一是化学计量关系中的分数和倍数，这些数不是测量所得，它们的有效数字位数可视为无限多位；另一种情况是关于pH、pK和lgK等对数值，其有效数字的位数仅取决于小数部分的位数，因为整数部分只与该真数中的10的方次有关。
例如，pH=13.15为两位有效数字，整数部分13不是有效数字。
若将其表示成[H+]=7.1×10-14，就可以看出13的作用仅是确定了[H+]在10-14数量级上，其数学意义与确定小数点位置的“0”相同。
在滴定分析中，实验数据的记录只应保留一位可疑数字，结果的计算和数据处理均应按有效数字的计算规则进行。
在进行加减运算时，有效数字取舍以小数点后位数最少的数值为准。
例如，0.0231、24.57和1.16832三个数相加，24.57的数值小数点后位数最少，故其他数值也应取小数点后两位，医学教育网搜集整理其结果是：0.02+24.57+1.17=25.76
在乘除运算中，应以有效数字最少的为准。
例如，0.0231、24.57和1.16832三个数相乘，0.0231的有效数字最少，只有3位，故其他数字也只取3位。
运算的结果也保留3位有效数字：
0.0231×24.6×1.17=0.665
在对数运算中，所取对数的位数应与真数的有效数字位数相同。
例如，lg9.6的真数有两位有效数字，则对数应为0.98，不应该是0.982或0.9823.又如[H+]为3.0×10-2mol·L-1时，PH应为1.52.
正确运用有效数字规则进行运算，不但能够反映出计算结果的可信程度，而且能大大简化计算过程。
在滴定分析中一般常采用四位有效数字。
6、在科学计算中，一般认为误差的来源有几种？列举在数值分析课中主要讨论误差。
数值计算中一个基本的手段是近似，所以就有了各种误差。
误差来源有四种：模型误差，观测误差，截断误差，舍入误差。
一般分为两类，第一类是固有误差（包括模型误差和观测误差），第二类是计算误差（包括截断误差和舍入误差）。
计算方法课中主要讨论计算误差。
这是因为在用计算机解决数学问题时，常常用“有限代替无穷，用近似代替准确”。
例如，解决连续性问题时通常要将其转
化为离散问题求解，这将引起截断（方法）误差；由于机器数的位数有限，计算机表示数据时一般带有舍入误差。
下面不全面列举出本课程内容涉及的误差
线性方程组直接求解方法——舍入误差
多项式插值方法——插值误差
数据拟合方法——残差
数值积分方法——求积误差
微分方程数值解方法——局部截断误差
7、有效数字的意义是什么？都有哪些作用？
我的思考：有效数字是一个数据的在保证最小误差的情况下，取的一个能够在计算中发挥其有效作用的近似值。
有效数字的作用在于，最大精度地去发挥这个数值在计算中的作用，而又不会对计算结果造成太大影响，是计算过程较为简单。
8、什么是算法？怎样判断一个算法是否稳定？
答：（1）由基本运算及运算顺序的规定所构成的完整解题步骤，称为算法；（2）研究一个计算公式是否稳定，只要假设初始值有误差ɛ，中间不再产生新误差，考察由ɛ引起的误差积累是否增长，如不增长就认为是稳定的，否则是不稳定的。
9、有效数字的概念是如何抽象而来的？
答：有效数字位数与计算近似值x的误差这两个概念是通过末位数半个单位相联系的。
有效数字概念正是根据日常生活中的“四舍五入”原则抽象而来的。
若近似值x的绝对误差限是某一位上半个单位，该位到x的第一位非零数字一共有n位，则称这一近似值具有n 位有效数字。
10、什么是误差和误差限?他们有什么区别?
答:误差是近似值与精确值的差,可以是正的,也可以是负的,它是一个不确定的数.误差限有误差下线和误差下限之分,误差上限是误差的最大值,误差下限是误差的最小值,他有一个范围,不是一个数,也是不确定的.
11、在一项浩大的工程中，要进行成千上万的数值分析，也许一处误差就会造成整个庞大工程的失误，除了课本中介绍的几种预防误差的方法，还有哪些方法？有没有一种可以预防数值误差的合理计算系统？
12、什么是截断误差和舍入误差？它们分别对应算法的那种性质？
截断误差：是指计算某个算式时没有精确的计算结果，如积分计算，无穷级数计算等，使用极限的形式表达的，显然我们只能截取有限项进行计算，此时必定会有误差存在。
舍入误差：是指由于计算机表示位数的有限，很难表示位数很长的数字，这时计算机就会将其舍成一定的位数，引起舍入误差，每一步的舍入误差是微不足道的，但经过传播和积累，舍入误差可能会淹没所要的真解。
13、如何评价一个数值算法的好坏？
答：设计数值算法主要是处理计算精度与计算速度两个问题，计算精度就是计算结果的可靠性问题，一般应具有收敛性与数值稳定性，计算速度就是算法的效率问退，它包含计算量和存储量两个方面。
14、利用级数公式π/4=1－1/3+1/5－1/7...........可计算出无理数π的近似值，由于交错级数的部分和数列Sn在其极限上下值摆动，故截断误差将小于第一个被舍去的项的绝对值|an+1|。
试分析，为得到级数三位有效数字近似值。
应取多少项求和？
15、在实际生活中有些测量值不是经过四舍五入得来的，如何判断它们的相对误差？比如有一种四舍六入五保留的取值方法应该怎样计算它的相对误差？
16、在计算机软件编程中应在怎样设计一种软件将误差降到最小？
答：首先，要注意运算次序，防止大数“吃掉”小数；其次，要尽量减少运算次数，节省计算时间，减小所占空间，提高精确率。