当前位置：文档之家 › 数值分析2-1-0

数值分析2-1-0

  • 格式：ppt
  • 大小：937.00 KB
  • 文档页数：21

下载文档原格式

  / 21

合集下载

数值分析实验报告2

数值分析实验报告2

实验报告实验项目名称函数逼近与快速傅里叶变换实验室数学实验室所属课程名称数值逼近实验类型算法设计实验日期班级学号姓名成绩512*x^10 - 1280*x^8 + 1120*x^6 - 400*x^4 + 50*x^2 - 1并得到Figure，图像如下：实验二：编写程序实现[-1,1]上n阶勒让德多项式，并作画（n=0,1，…，10 在一个figure中）。

要求：输入Legendre(-1,1,n)，输出如a n x n+a n-1x n-1+…多项式。

在MATLAB的Editor中建立一个M-文件，输入程序代码，实现勒让德多项式的程序代码如下：function Pn=Legendre(n,x)syms x;if n==0Pn=1;else if n==1Pn=x;else Pn=expand((2*n-1)*x*Legendre(n-1)-(n-1)*Legendre(n-2))/(n);endx=[-1:0.1:1];A=sym2poly(Pn);yn=polyval(A,x);plot (x,yn,'-o');hold onend在command Windows中输入命令：Legendre（10），得出的结果为：Legendre(10)ans =(46189*x^10)/256 - (109395*x^8)/256 + (45045*x^6)/128 - (15015*x^4)/128 + (3465*x^2)/256 - 63/256并得到Figure，图像如下：实验三：利用切比雪夫零点做拉格朗日插值，并与以前拉格朗日插值结果比较。

在MATLAB的Editor中建立一个M-文件，输入程序代码，实现拉格朗日插值多项式的程序代码如下：function [C,D]=lagr1(X,Y)n=length(X);D=zeros(n,n);D(:,1)=Y';for j=2:nfor k=j:nD(k,j)=(D(k,j-1)- D(k-1,j-1))/(X(k)-X(k-j+1));endendC=D(n,n);for k=(n-1):-1:1C=conv(C,poly(X(k)));m=length(C);C(m)= C(m)+D(k,k);end在command Windows 中输入如下命令：clear,clf,hold on;k=0:10;X=cos(((21-2*k)*pi)./22); %这是切比雪夫的零点Y=1./(1+25*X.^2);[C,D]=lagr1(X,Y);x=-1:0.01:1;y=polyval(C,x);plot(x,y,X,Y,'.');grid on;xp=-1:0.01:1;z=1./(1+25*xp.^2);plot(xp,z,'r')得到Figure ，图像如下所示：比较后发现，使用切比雪夫零点做拉格朗日插值不会发生龙格现象。

数值分析课件1

数值分析课件1

提出数值问题
数值问题是指有限个输入数据（问题 的自变量、原始数据）与有限个输出数据 （待求解数据）之间函数关系的一个明确 无歧义的描述。这正是数值分析所研究的 对象。
数值问题举例
dy = x +y2 dx y ( 0) = y 0 x ∈ [0, 1]
是用一阶常微分方程初值问题表示的 数学模型，要求无穷多个输出，因而它不是 数值问题 。但当我们要求出有限个点处函 数值的近似值时，便成为一数值问题。
设计高效可靠的算法
计算方法的任务之一就是提供求得数值问 题近似解的方法—算法。
算法：指把对数学问题的解法归结为只有 加、减、乘、除等基本运算，并确定运算次序 的完整而准确的描述。
算法的可靠性：算法的可靠性包括算法的收 敛性、稳定性、误差估计等几个方面。这些是
数值分析研究的第二个任务。
一个算法在保证可靠的大前提下再评价其 优劣才是有价值的。 算法的优劣评价：可靠算法的优劣，应该考 虑其时间复杂度（计算机运行时间）、空间 复杂度（占据计算机存储空间的多少）以及 逻辑复杂度（影响程序开发的周期以及维护 ）。这是数值分析研究的第三个任务。
e − e = (e − en ) + (en − e)
* *
二、误差的度量
1) 2) 3) 4)
绝对误差 相对误差 有效数字 各种度量之间的关系
1. 绝对误差

绝对误差定义：近似值减准确值
* ∆
x − x= e( x * ) * * e ( x ) 在不引起混淆时， 简记 为 e 。
• 绝对误差限：
位有效数字。如 A = sin 29 20′ = 0.4900 设其近似值a=0.484，其相对误差为：
0.4900 − 0.484 1 = 0.012397 < 0.0125 = × 101− 2 0.484 2× 4

《数值分析》所有参考答案

《数值分析》所有参考答案
解：
等价三角方程组
， ，
11.设计算机具有4位字长。分别用Gauss消去法和列主元Gauss消去法解下列方程组，并比较所得的结果。
解：Gauss消去法
回代
列主元Gauss消去
15.用列主元三角分解法求解方程组。其中
A= ,
解：
等价三角方程组
回代得
， ， ，
16.已知 ，求 ， ， 。
解：
, ,
17.设 。证明
，(II)
，
当 时
当 时
迭代格式(II)对任意 均收敛
3) ，
构造迭代格式 (III)
，
当 时
当 时
迭代格式(III)对任意 均收敛
4)
取格式(III)
， ， ，
4.用简单迭代格式求方程 的所有实根，精确至有3位有效数。
解：
当 时， ，
1 2
当 时
，
，
， ，
1)
迭代格式 ，
,
当 时， ，
任取 迭代格式收敛于
是中的一种向量范数。
解：
当 时存在 使得
，
,
所给 为 上的一个范数
18.设 。证明
(1) ;
(2) ;
(3) 。
解：(1)
(2)
(3)
19.设
A=
求 ， ， 及 ， 。
解： ,
Newton迭代格式
,
20.设 为 上任意两种矩阵(算子)范数，证明存在常数
， 使得
对一切 均成立。
解：由向量范数的等价性知道存在正常数 使得
,
=0.187622
[23.015625 , 23.015625+0.187622]

数值分析_第2章

数值分析_第2章

证：由1。 f '( x) C[a, b],由2。 f '( x)不变号，故f ( x) 知 知 单调，再由3。 唯一的 [a, b]，使f ( ) 0. 知
由1 3 知f ( x)在[a, b]上必属于下列四种情形之一：
。 。
f ''( x) 0 f (a) 0, f (b) 0, f '( x) 0(增) f ''( x) 0
二．收敛性：
mn . n .
◆判定二分次数：
1 lim n 1 b0 a0 0 n 2
1 对 0，若要求 mn n 1 b0 a0 2
b0 a0 则2 n log 2 1与取整的 1抵消 .
定理1.(单点法收敛的充分条件) 设f ( x)在[a, b]上二阶 可导，且满足：
。 1. f ''( x)在[a, b]上不变号（凹凸不变性）；
2。 f '( x)在[a, b]上不为0（单调性）； . 3。 f (a) f (b) 0; . 4。取x0 [a, b], 使f ( x0 ) f ''( x0 ) 0.x1 [a, b], f ( x1 ) f ( x0 ) 0. . 则由(6)所得 xn 单调收敛于f ( x) 0在[a, b]上的唯一根。
列表计算：
n
0 1 2 3 4 5
xn
2 1 1.33333 1.40000 1.41176 1.40378
2
f ( xn )
2 -1 -0.22223 -0.04000 -0.00692
hn

数值分析 第1章 绪论 张铁版

数值分析 第1章 绪论 张铁版

3.绝对值太小的数不宜作除数 例7 仿计算机，采用3位十进制，用消元法求解方程组
1.00105 x 1.00y 1.00 1.00x 1.00y 2.00
105 x 1105
1.00001 0.9999899
(2) (1) 10
1.00105 x 1.00 y 1.00 解： x得， 消 5 5 5 (1.00 1.0010 ) y (2.00 1.0010 )
算法1:直接计算 n(n 1) 乘法次数：1+2+ +n= 2 加法次数：n
算法2：秦九韶算法(Hernor算法)：
S n an , S k xS k 1 ak , (k n - 1, ,0) P ( x) S . 0 n
乘法次数：n，加法次数：n
( ) n1 Rn ( x) x (n 1)! f
( n1)
截断误差：
舍入误差 R 3.14159 0.0000026. 数制转换、机器数.
§1.3 绝对误差、相对误差与有效数字
定义1 绝对误差，简称误差：
e x * x, 其中x为准确值x *的近似值.
误差限： | e | 的一个上界,即 x * x .
5
y
2 105 1105
1.00 10 x 1.00 y 1.00 x* 0.00, y* 1.00 y 1.00
错.为什么,怎么办？
４.简化计算程序，减少运算次数 减少运算次数可以不但节省时间，而且减少舍入误差. 例8 计算多项式的值 Pn ( x) an x n an1x n1 a1x a0 .
1 1 e1 * I 9 0.0684, ( I 9 ( ) 0.0684) 2 10 10 ( B) * * I n1 1 (1 I n ), n 9,8,,1. n

数值分析(清华大学出版社)第1,2章

数值分析(清华大学出版社)第1,2章
20
2.
x 的相对误差是
x x d x er ( x ) d ln x x x

它是对数函数的微分。
设 u = xy , 则 lnu=lnx+lny , 因而 dlnu = dlnx + dlny
e r ( u ) e r ( x ) e r ( y ) r ( u ) r ( x ) r ( y )

即 m- n = - 2, m=1, n = 3, 所以 x = 3.14 作为 近似值 时, 就有3 位有效数字。

16
四、 相对误差限与有效数字的关系
定理1
设近似值
x 0.a1a2 an 10m
有n 位有效数字， a1 0 。则其相对误差限为 1 n 1 r (x ) 10 2a1 x 0.a1a2 an 10m 故 证明
a1 10
m 1
| x | (a1 1) 10

m 1
r (x )
x x x


0 .5 10m n 1 10 n1 a1 10m 1 2a1
17
定理2 设近似值 x 0.a1a2 an 10 的相对误差限
m

1 10 n1 ,则它至少有n 位有效数字。 不大于 2( a 1) 1

25
对多元函数 y f ( x1 , x2 , , xn ), 自变量的近似值为 x1 , x2 ,, xn , y 的近似值为 y f ( x1 , x2 , , xn ),
y 的运算误差为 函数值
e ( y ) e[ f ( x1 , x2 , , xn )] df ( x1 , x2 , , xn )

数值分析课后习题答案

数值分析课后习题答案

x2 6.6667x2 8.205
再解
1
15 56
x31.785,7得 x35.769
1 25069x4 0.47847x4 1.4872
1 x5 5.3718 x5 5.3718
2-10.证明下列不等式:
(1)x-yx-z+z-y; (2)|x-y|x-y;
证明 (1)x-y=(x-z)+(z-y)x-z+z-y
b.用Gauss消元法
102 x y 1 x y 2
回代得解: y=1, x=0.
102 x Байду номын сангаасy 1
100y 100
再用列主元Gauss消元法
102 x y 1 x y 2
回代得解: y=1, x=1.
x y
y 1
2
2-8.用追赶法求解方程组:
4 1
x1 100
1 4 1
x2 0
3-8.判定求解下列方程组的SOR方法的收敛性.
2 1 0 0 x1 1
1
0 0
2 1 0
1 2 1
0 12
x2 x3 x4
0 00
解 直接可验证系数矩阵A是负定矩阵,所以-A是对称
1-3.为了使101/2的相对误差小于0.01%,试问应取几位 有效数字?
解 因为101/2=3.162…=0.3162…10,若具有n位有效 数字,则其绝对误差限为0.5 101-n ,于是有
r=0.5101-n/3.162…<0.5101-n/3<0.01% 因此只需n=5.即取101/2=3.1623
1 2
0
12 1,
1 2
1 2
0
12

数值分析课程第五版课后习题答案(李庆扬等)(OCR)

数值分析课程第五版课后习题答案(李庆扬等)(OCR)

根是ｘ，，２…，ｘ－，且Ｖ。ｘ，ｘ…·，ｘ）＝Ｖ，Ｃｘ６，ｘ…·）（ｘ－ｘ）…（ｘ－ｘ）。
Ｖ，（ｘｏ，ｘ，…ｘ－ｘ）＝１１】 －ｘ，）用ａ－ｘ，）
［证明］由
可得求证。
＝Ｖ，（Ｃｘ８，ｘ，…，ｘＸ））１１（ｘ－ｘ）
２、当ｘ＝１－１，２时，ｆ（ｘ）＝０，－３．４，求ｆ（ｘ）的二次插值多项式。
Ｌ，（ｘ）＝ｙ％（（ｘｘ６－－ｘｘ，）（（ｘｘ－２ｘ－ｘ２２））
ｙ＝ｆ（ｘ）＝ｆ０．５）＝－０．６９３１４７，ｙ２＝ｆ（ｘ）＝ｆ（０．６）＝－０．５１０８２６，则
Ｌ２（ｘ）＝ｙ。 （ｘ－ｘ）（ｘ－ｘ２）
（ｘ６－ｘ）ｘ－ｘ）
（ｘ－ｘ）（ｘ－ｘ）
（ｘ－ｘ）（ｘ－ｘ２）
（ｘ－ｘｏ）（ｘ－ｘ） （ｘ２－ｘｏ）（ｘ２－ｘ）
＝－０．９１６２９１×．（０（．ｘ４－０－．０５．）５（）ｘ（－００．．６４）－０．６－．
３０—＋２—９．ｘ９５８３ｘ３１￣０２＇＝０．８３３６×１０４
１４、试用消元法解方程ｘ组１＋１０＂ｘ＝１００
ｘ＋ｘ２＝２
，假定只有三位数计算，问结果是否
可靠？
［解］精确解为ｘ１＝０１００－＊１ １０＂－２ ，当使用三位数运 算时，得到
ｘ ＝１，ｘ２＝１，结果可靠。
１５、已知三角形面积ｓ＝ｓ去＝ ａｂｓｉｎｃ，其中ｃ为弧度，０＜ｃ＜ 且测量ａ，ｂ，ｃ
位有效数字；ｘ＝５６．４３０有５位有效数字；ｘ＝７×１０有２位有效数字。 ４、利用公式（３．３）求下列各近似值的误差限，其中ｘ，ｘ；，ｘ，ｘ；均为第３题所给
的数。
（１）ｘ＋ｘ２＋ｘ：
ｅ（ｘ＋ｘ写＋ｘ）＝＞
［解］
Ｅ（ｘ）＝Ｅ（ｘ）＋Ｅ（ｘ）＋Ｅ（ｘ；）
３＋ｔｘ１０＝１．０５×１０３
（２）ｘｘｘ；

《数值分析简明教程》第二版（王能超编著）课后习题答案高等教育出版社

《数值分析简明教程》第二版（王能超编著）课后习题答案高等教育出版社

《数值分析简明教程》第⼆版（王能超编著）课后习题答案⾼等教育出版社0.1算法1、（p.11，题1）⽤⼆分法求⽅程013=--x x 在[1,2]内的近似根，要求误差不超过10-3.【解】由⼆分法的误差估计式311*10212||-++=≤=-≤-εk k k a b x x ，得到100021≥+k .两端取⾃然对数得96.812ln 10ln 3≈-≥k ，因此取9=k ，即⾄少需2、（p.11，题2）证明⽅程210)(-+=x e x f x在区间[0,1]内有唯⼀个实根；使⽤⼆分法求这⼀实根，要求误差不超过21021-?。

【解】由于210)(-+=x e x f x ，则)(x f 在区间[0,1]上连续，且012010)0(0<-=-?+=e f ，082110)1(1>+=-?+=e e f ，即0)1()0(⼜010)('>+=x e x f ，即)(x f 在区间[0,1]上是单调的，故)(x f 在区间[0,1]内有唯⼀实根.由⼆分法的误差估计式211*1021212||-++?=≤=-≤-εk k k a b x x ，得到1002≥k .两端取⾃然对数得6438.63219.322ln 10ln 2=?≈≥k ，因此取7=k ，即⾄少需⼆分0.2误差1．（p.12，题8）已知e=2.71828…，试问其近似值7.21=x ，71.22=x ，x 2=2.71，718.23=x 各有⼏位有效数字？并给出它们的相对误差限。

【解】有效数字：因为11102105.001828.0||-?=<=- x e ，所以7.21=x 有两位有效数字；因为12102105.000828.0||-?=<=- x e ，所以71.22=x 亦有两位有效数字；因为3310210005.000028.0||-?=<=- x e ，所以718.23=x 有四位有效数字；%85.17.205.0||111=<-=x x e r ε； %85.171.205.0||222=<-=x x e r ε； %0184.0718.20005.0||333=<-=x x e r ε。

数值分析实验报告2

数值分析实验报告2

实验报告一、实验名称复合梯形求积公式、复合辛普森求积公式、龙贝格求积公式及自适应辛普森积分。

二、实验目的及要求1. 掌握复合梯形求积计算积分、复合辛普森求积计算积分、龙贝格求积计算积分和自适应辛普森积分的基本思路和步骤.2. 培养Matlab 编程与上机调试能力. 三、实验环境计算机，MATLAB 软件 四、实验内容1.用不同数值方法计算积分94ln 10-=⎰xdx x 。

（1）取不同的步长h 。

分别用复合梯形及复合辛普森求积计算积分，给出误差中关于h 的函数，并与积分精确指比较两个公式的精度，是否存在一个最小的h ，使得精度不能再被改善。

（2）用龙贝格求积计算完成问题（1）。

（3）用自适应辛普森积分，使其精度达到10-4。

五、算法描述及实验步骤1.复合梯形公式将区间[a,b]划分为n 等份，分点x k =a+ah,h=(b-a)/h,k=0,1,...,n ，在每个子区间[x k ,x k +1](k=0,1,...,n-1)上采用梯形公式（1.1），得)]()([2)(b f a f ab dx x f b a+-≈⎰ （1.1） )]()(2)([2)]()([211110b f x f b f hx f x f h T n k k k n k k n ++=+=∑∑-=+-= （1.2）),(),(12)(''2b a f h a b f R n ∈--=ηη（1.3） 其中Tn 称为复合梯形公式，Rn 为复合梯形公式的余项。

2.复合辛普森求积公式将区间[a,b]划分为n 等份，在每个子区间[x k ,x k +1](k=0,1,...,n-1)上采用辛普森公式（1.4），得)]()2(4)([6b f ba f a f ab S +++-=（1.4） )]()(2)(4)([6)]()()([611102/112/11b f x f x f b f hx f x f x f h S n k k n k k k k n k k n +++=++=∑∑∑-=-=+++-= （1.5） ),(),()2(180)()4(4b a f h a b f R n ∈-=ηη （1.6）其中Sn 称为复合辛普森求积公式，Rn 为复合辛普森求积公式的余项。

数值分析参考答案第二章

数值分析参考答案第二章

第二章插值法解:X 0 1,X ! 1,X 2 2,f(X °)0,f (X 1)3,f(X 2) 4 7 l 0(x) (X X 1)(x X 2)1(X 1)(X 2)(X 0 X 1）（X X 2)2h(x) (X X 0）（X X 2)-(x 1)(x 2)(X 1 X0XX1 X 2) 6l 2（X(X x °)(x X 1) 1 “1)(x 1) (x (X 2 X °)(X 2 X i ) 3则二次拉格朗日插值多项式为2L 2（X ） yh （x ）k 02.给出f （x ） In x 的数值表X 0.4 0.5 0.6 0.7 0.8 Inx -0.916291-0.693147-0.510826-0.356675-0.223144用线性插值及二次插值计算 In0.54的近似值。

解：由表格知，x 0 0.4, x 1 0.5, x 2 0.6, x 3 0.7, x 4 f （x 。

0.916291, f (x 1) 0.693147 f(X 2) 0.510826, f (x 3) 0.356675 f(X 4)0.223144若采用线性插值法计算In0.54即f （0.54）， 则 0.5 0.54 0.61 当 x 1, 1,2 时，f (x)0, 3,4，求f （x ）的二次插值多项式。

(X4 一2)(XX /VX2l 1(x) X X I X 2 X 2 10(x 0.6) l 2（XX X 2 X 1X 110(x 0.5)L i (x) f(xOh(x) f(X 2)〔2(x)6.93147(x 0.6) 5.10826( x 0.5) ^(0.54)0.6202186 0.620219L 2(X ) f(X 0)I °(X ) f(xj 1(x) f(X 2)J(x)L 2(0.54) 0.61531984 0.6153203.给全cosxO x 90'的函数表，步长h 1 究用线性插值求cosx 近似值时的总误差界。

《数值分析》2019-2020学年第一学期期末试卷

《数值分析》2019-2020学年第一学期期末试卷

吉林大学《数值分析》2019-2020学年第一学期期末试卷一、单项选择题（每小题3分，共15分） 1. 已知近似值1x ，2x ，则()12,x x ()=A. ()()2112x x x x +B. ()()12x x +C. ()()1122x x x x +D. ()()12x x2. 已知求积公式()()211211()(6362)f x dx f Af f ≈++∫，则A ＝（ ） A ． 16 B. 13 C. 12 D. 233. 已知，则化为2112A ⎡⎤=⎢⎣⎦⎥A 为对角阵的平面旋转变换角θ＝（ ） A ．6πB.4πC.3πD.2π4. 设求方程()0f x =的根的切线法收敛，则它具有（ ）敛速。

A ． 线性 B. 超越性 C. 平方 D. 三次5. 改进欧拉法的局部截断误差为（ ）A ． B. ()5O h ()4O h C. ()3O h D. ()2O h二、填空题（每小题3分，共15分）1. π的近似值3.1428是准确到 近似值。

2. 满足()a a f x x =，()b b x x =，()c f c f x x =的拉格朗日插值余项为 。

3. 用列主元法解方程组时，已知第2列主元为()142a 则()142a ＝ 。

4．乘幂法师求实方阵 的一种迭代方法。

5. 欧拉法的绝对稳定实区间为 。

三、计算题（每小题12分，共60分） 1. 用已知函数表x 0 1 2y 1 2 5求抛物插值多项式，并求1()2f 的近似值。

2. 用紧凑格式解方程组 123410114130141x x x −⎡⎤⎡⎤⎢⎥⎢⎥−−=⎢⎥⎢⎥⎢⎥⎢⎥−⎣⎦⎣⎦⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦⎤⎥3. 已知方程组123210113110121x x x ⎡⎤⎡⎤⎡⎢⎥⎢⎥⎢=−⎢⎥⎢⎥⎢⎢⎥⎢⎥⎢⎣⎦⎣⎦⎣⎥⎥⎦) (1) 证明高斯－塞德尔法收敛；(2)写出高斯－塞德尔法迭代公式； (3) 取初始值，求出()(00,0,0TX=()1X4. 用复化辛卜公式计算积分4n =1011dx x +∫，并估计误差。

第一章数值分析(误差分析)

第一章数值分析(误差分析)
*
* e x x * e r * * x x x x* er 则称η 为 x* 的相对误差限。 x
如果
这时 x=10,
x*=10±1;
2019/3/13
第一章 绪论与误差分析
2
本章内容安排
1. 目的意义：了解计算数学的背景知识；掌握误 差的基本知识 2.重 点：误差来源、误差表示、误差传播 及算法设计原则 3.难 点：有效数字 4.内容分配： 第 1 次：§1 计算数学研究的对象和内容 第
§2 误差的来源和分类 2 次：§3 误差的表示 §4 误差的传播 §5 算法设计的若干原则
由于计算机的字长有限，参加运算的数据以及计算结 果在计算机上存放时，计算机会按舍入原则舍去每个数据 字长之外的数字，从而产生误差，这种误差称为舍入误差 或计算误差。 例如，在十进制十位的限制下，会出现 (1.000002)2-1.000004=0
这个结果是不准确的，准确的结果应是 (1.000002)2-1.000004 =1.000004000004-1.000004=4×10-12 这里所产生的误差就是计算舍入误差。 在数值分析中，一般总假定数学模型是准确的，因而 不考虑模型误差和观测误差，主要研究截断误差和舍入误 差对计算结果的影响。
则有误差限 |x-x*|≤1= εx ,
虽然εy是εx 的3倍，但在1000内差3显然比10内差1更精确 些。这说明一个近似值的精确程度除了与绝对误差有关 外，还与精确值的大小有关，所以这时可以用相对误差 来比较这两个近似数的准确度。
2019/3/13 19
第一章 绪论与误差分析
e x x 定义1 .2 记 er x x 则称其为近似值 x *的相对误差。 由于 x 未知， 实际使用时总是将 x * 的相对误差取为

《数值分析》第二章答案

《数值分析》第二章答案

习题21. 分析下列方程各存在几个根，并找出每个根的含根区间：(1) 0cos =+x x ； (2) 0cos 3=-x x ； (3) 0sin =--x e x ； (4) 02=--x e x 。

解：(1) 0cos =+x x (A) x x x f cos )(+= ，0sin1)(≥-='x x f ,),(∞-∞∈x10cos 0)0(=+=f ,01cos 1)1cos(1)1(<+-=-+-=-f ∴ 方程(A) 有唯一根 ]0,1[*-∈x (2) 0cos 3=-x x (B) x x x f c o s 3)(-=,0sin 3)(>+='x x f , ),(∞-∞∈x 时010c o s03)0(<-=-⨯=f ，01cos 31cos 13)1(>-=-⨯=f ∴ 方程(B) 有唯一根 ]1,0[*∈x (3)sin =--xex (C)xex -=sinx x f sin )(1=, xex f -=)(2方程(C)有无穷个正根，无负根 在[22,2πππ+k k ] 内有一根 )(1k x ，且0]2[lim )(1=-∞→πk x k k在[ππππ++k k 2,22]内有一根)(2k x ，且0])12([lim )(2=+-∞→πk x k k (示图如下) 3,2,1,0=k)(2x f x(4)02=--xex(D) xex-=2,)(21x x f = xex f -=)(2方程(D) 有唯一根 ]1,0[*∈x 当 0<x 时 (D)与方程2x ex -=- (E) 同解 当 0<x 时 (E)无根 2. 给定方程 012=--x x ； (1)(2)若在[0 , 2]上用二分法求根，要使精确度达到6位有效数，需二分几次？ 解：012=--x x1) 01)(2=--=x x x f 1)1(-=f , 025.0)5.1(<-=f ,1)2(=f]2,5.1[*∈x, 618034.1251*=+=x)(5.1- 1.75(+) 2(+) )(5.1- 1.625(+) 1.75(+) )(5.1-1.5625(+) 1.625(+))(5625.1- )(59375.1-1.625(+)1102103125.02)5625.1625.1(-⨯<=-6.159375.1*≈≈x2位有效近似值为 1.6 2)00==a a ， 20==b b)(21k k k b a c +=kk k a b c x 2121*=-≤-+5102121-⨯≤k，51102≥-k60.162ln 10ln 51=≥-k∴ 只要2等分18次3. 为求0353=--x x 的正根，试构造3种简单迭代格式，判断它们是否收敛，且选择一种较快的迭代格式求出具有3位有效数的近似根。

(仅供参考)数值分析课后第二章习题解答

(仅供参考)数值分析课后第二章习题解答

6第二章 习题解答一、习题解答1 用二分法求解下列方程，要求误差不超过10-5 （1）x – ln x=2 在区间[2，4]内的根；（2）x e x– 1 = 0 在区间[0，1]内的根；（3）x 3 + 4x 2– 10 = 0 在区间[1，2]内的根。

解：（1）x 18 = 3.14618。

二分法程序2-1a=2;b=4;k=0;f=inline('x-log(x)-2'); ya=f(a);while (b-a)>0.00001x0=.5*(a+b);y0=f(x0); if ya*y0<0 b=x0; elsea=x0;ya=y0; end k=k+1; endformat long disp([x0,k])（2） x 17 = 0.567146。

（利用上面程序修改前两行） （3） x 17 = 1.365226。

2 证明方程 1 – x – sin x = 0 在区间[0，1]上有一根。

使用二分法求误差不大于41021−×的根需二分多少次？证明 令f (x ) = 1 – x – sin x ，则f (0) = 1，f (1)= – sin 1。

于是 f (0) f (1)< 0，故所给方程在区间[0，1]上必有根。

又因为x x f cos 1)(−−=′ 所以，函数 f (x ) 在区间[0，1]内单减。

故，方程在区间[0，1]内只有一个根。

利用二分法收敛定理，由411021201−+×≤−n 得 2n ≥ 104，所以二分法求根至少需14次二分计算能满足误差要求。

3 比较以下两种方法求 e x+ 10 x – 2 =0 的根到三位小数所需要的计算量。

（1） 在区间[0，1]内用二分法；（2） 用迭代法)2(1011n x n e x −=+，取初值 x 0 = 0。

解：（1）二分法迭代11次，x 11 = 0.090。

（利用第1题程序修改前两行） （2）不动点迭代5次，x 5 = 0.0905。

数值分析2(二分法)

数值分析2(二分法)

)
函数在[0,1]内有唯一零点, 故[0,1]是隔根区间。
n 1
函数值符号 f(0.5)= - 0.1006<0
区间 [0,0.5]
2
3
f(0.25)=0.3961>0
f(0.375)=0.1317>0
[0.25,0.5]
[0.375,0.5]
4
f(0.4375)=0.0113>0 [0.4375,0.5]
其中( x, y, z, t )表示点的当前位置。
参考: 《数值分析计算引论》P. 183
非线性科学是当今科学发展的一个重要研究 方向,而非线性方程的求根也成了一个不可缺 的内容。但是非线性方程的求根非常复杂。
高于４次的代数方程,不存在通用的求根公式
f ( x ) a0 a1 x an x (an 0)
性质:1. f(an)f(bn)<0; 2. bn – an = (b0 – a0)/ 2n
定理2.2 设x*是 f(x)=0在[a 0, b 0]内的唯一根,且 f(a 0)f(b 0)<0,则二分计算过程中, 各区间的中点数 列 1 xn (an bn )( n 0,1,2,) 2 满足 | xn – x*|≤ (b 0 – a 0)/ 2n+1 1 log2 10 3.3219 * 3 注记: 若要 | x n x | 10 2 b0 a0 1 b0 a0 3 只需 10 n log 2 n 1 3 2 2 10
《数值分析》2
非线性方程
二分法迭代
GPS (Global Positioning System) LBS (Location Based Service)
  1. 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
设 (t ) f (t ) 1 (t ) K (t x0 )(t x1 ) (1) ( x0 ) ( x1 ) 0, R( x) 取 K (2) ( x) 0 ( x x0 )( x x1 ) (1 ) ( 2 ) 0, ( ) 0, R( x) d2 (1), (2) (t ) f (t ) 1(t ) (t x0 )(t x1 ) 2 ( x x0 )( x x1 ) dt 2! R( x) (t ) f (t ) (3) ( x x0 )( x x1 ) 2! R( x) f ( ) ( ) 0, a b 12 ( x x0 )( x x1 )
（3）Newton插值公式
2 ( x) f ( x0 ) f ( x0 , x1 )( x x0 )
f ( x0 , x1 , x2 )( x x0 )( x x1 )
18
三、二次插值的三种插值形式
3. 逐次线性插值
f ( x0 ) 1 1 ( x ) x 0 x1 f ( x1 )
( x x0 )( x x2 ) l1 ( x) ( x1 x0 )( x1 x2 )
(2) 二次Lagrange插值公式
2 ( x) y0l0 ( x) y1l1 ( x) y2l2 ( x)
17
三、二次插值的三种插值形式
1 ( x) f ( x0 ) f ( x0 , x1 )( x x0 )
称 li(x) (i 0,1 ) 为线性Lagrange插值基函数 (3)Lagrange插值公式 1 ( x) y0l0 ( x) y1l1 ( x)
8
三、线性插值的三种插值形式
2. Newton插值（点斜式）
(1) Newton插值公式
f ( x1 ) f ( x0 ) 1 ( x) f ( x0 ) ( x x0 ) x1 x0
3、行列式形式：
f ( x0 ) 1 1 ( x) x0 x1 f ( x1 )
x x0 x x1
四、线性插值的误差
设f(x) c1[a,b],f 在区间[a,b]上存在，则对于x [a,b]有 (x)
f ( ) 1、 ( x) f ( x) 1 ( x) R ( x x0 )( x x1 ) , a b 2!
A B( x x0 )
f ( x0 ) f ( x1 ) f ( x0 , x1 ) x0 x1
x x0 x x1
3、行列式形式：
线性插值的误差
f ( x0 ) 1 1 ( x ) x 0 x1 f ( x1 )

1、 ( x) f ( x) 1 ( x) f ( ) ( x x0 )( x x1 ) , a b R 2! ( x1 x0 ) 2 max f ( x) , a x b 2、R( x) 8 15
3
问题的引入：
并且用P(x)近似代替f(x)
这就是插值问题, (1)式为插值条件,
称函数P(x)为函数f(x)的插值函数 如果P(x)为多项式函数, 则称之为插值多项式
点 xi , i 0,1,2, ,n,称为插值节点
区间[a,b]称为插值区间
例如：函数y sin x,若给定[0,π]上5个等分点 其插值函数的图象如图
x : x0 , x1, y : y0, y1 , yi f ( xi )
--------( )
构造y f(x)的插值函数(x)
使 (x)是不高与于1次的多项式
(xi ) yi f(xi ), i 0,1
二、方法：
令 且满足
1 ( x) ax b
--------(1) --------(2)
( x1 x0 ) max f ( x) , a x b 2、R( x) 8
2
10
四、线性插值的误差
f ( ) R ( x x0 )( x x1 ) , a b 分析：1、 ( x) f ( x) 1 ( x) 2! f ( ) R( x) f ( x) 1 ( x) ( x x0 )( x x1 ) , a b 2! 2! R( x) f ( ) ( ) 0, a b ( x x0 )( x x1 ) 2! R( x) (t ), 使 ( ) 0, ( ) f ( ) ( x x0 )( x x1 )
f ( x0 ) f ( x1 ) f ( x0 , x1 ) x0 x1
（2）1阶均差商
（3）Newton插值公式又可表示为
1 ( x) f ( x0 ) f ( x0 , x1 )( x x0 )
A B( x x0 )
9
三、线性插值的三种插值形式
x x0 x x1 1 ( x) y0 y1 x0 x1 x1 x0
四、线性插值的误差
f ( ) R( x) f ( x) 1 ( x) ( x x0 )( x x1 ) 2!
( x1 x0 ) 2 证明 ：2、 R( x) max f ( x) , a x b 8 f ( )
R( x) 2! max f ( x) ( x x0 )( x x1 ) ( x x0 )( x1 x )
解：
R1 (0.3367) 0.92 105
14
线性插值的三种插值形式
（小结）
1.拉格朗尔插值（两点式） 1 x xi i 0,1 1 ( x) y0l0 ( x) y1l1 ( x) li (x )
0 x xi
2. Newton插值（点斜式）
1 ( x) f ( x0 ) f ( x0 , x1 )( x x0 )
1 x xi i 0,1 li (x ) 0 x xi x xi i 0,1,2 x xi
( x x1 )( x x2 ) l0 ( x ) ( x0 x1 )( x0 x2 ) ( x x0 )( x x1 ) l2 ( x ) ( x2 x0 )( x2 x1 )
用线性插值计算 sin 0.3367的值，并估计误差。
x0 0.32, x1 0.34 x x0 x x1 sin x y0 y1 x0 x1 x1 x0 0.3367 0.34 0.3367 0.32 sin 0.3367 0.314567 0.333487 0.32 0.34 0.34 0.32 0.330356 1 R1 ( x) M 2 ( x x0 )( x x1 ) M 2 max sin x 0.3335 2
--------(1)
2 ( x) a0 a1 x a2 x 2
2 ( xi ) yi
i 0,1,2
--------(2)
16
1 ( x) y0l0 ( x) y1l1 ( x)
三、二次插值的三种插值形式
1.拉格朗尔插值 （1）构造插值基
1 li (x ) 0
解决方法：
假如可以通过实验或测量, 可以获得f(x)在区间[a,b] 上的一组n 1个不同的点
a x0 x1 x2 xn b
上的函数值 yi f ( xi ),
i 0,1,2,, n
能否存在一个性能优良、便于计算的函数
比如多项式函数P(x),满足
P( xi ) yi i 0,1,2 ,, n ------(1)
x x0 x x1
f ( x0 ) 1 (1) 01 ( x) x0 x1 f ( x1 ) f ( x0 ) 1 (2) 02 ( x) x0 x2 f ( x2 )
x x0 x x1 x x0 x x2
01 ( x) 1 (3) 012 ( x) x1 x 2 02 ( x)
2 max f ( x) ( x x0 ) ( x1 x) 2 [ ] 2 2
ab ab 2
( x1 x0 ) 2 max f ( x) 8
13
五、例题：
已知 sin 0.32 0.314567, sin 0.34 0.333487, sin 0.36 0.35227
4
sinxµ å IJ Öµ
1
yy
0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0
0
0.5
1
1.5 1.5

x x x
对于被插函数f(x)和插值函数P(x)
在节点xi 处的函数值必然相等
但在节点外P(x)的值可能就会偏离f(x) 因此P(x)近似代替f(x)必然存在着误差
计算物理
插值法
湖北大学
物理学与电子技术学院
1
问题的引入：




1 函数表达式过于复杂不便于计算, 而又需 要计算许多点处的函数值 2 仅有几个采样点处的函数值, 而又需要知 道非采样点处的函数值 …… 上述问题的一种解决思路：建立复杂函数或 者未知函数的一个便于计算的近似表达式.
2
问题的引入：
5
本章要点
用简单的函数(如多项式函数)作为一个 复杂函数的近似,最简单实用的方法就是插 值。 本章主要介绍有关插值法的一些基本 概念，及多项式插值的基础理论和几个常 用的插值方法:Lagrange插值、分段线性插 值、Newton插值、Hermite插值和三次样条 插。
6
§1. 线性插值
一、问题： 已知
附：Rolle定理
设f(x)在[a, b] 上连续，f ( x)在(a, b)上存在，且f (a) f (b), 则在[a, b]内至少有一点c, 使f (c) 0