《计算机数学基础(下)》数值分析部分辅导(1).

合集下载

计算机数学基础 ppt课件

计算机数学基础 ppt课件

王雷@版权所有
1. 2., E},{a},P,A),其 中P = {Aa,AaE,EaA}。 在许多的文法中,有多条产生式的左部相 同,可以将左部相同的产生式写成合并的 产生式形式。在此例文法G中,P中的前两 个产生式的左部相同,都是A,可以合并 为A a | aE,这样一来,P = { A a | aE, EaA}。 在许多情况下,只需要将文法的产生式写 出就可以表明该文法了。
文法所生成的语言,根据四种类型文法,也分为四种,即: 0型语言、1型语言、2型语言和3型语言。
Chomsky建立的形式语言理论对计算机科学的发展规律有 着深刻的影响,特别是对计算机程序设计语言的设计、编 译方法和计算复杂性等方面具有更大的作用。
王雷@版权所有
1.2.5 文法和语言的类型
定义1.11 设文法G = (V, T, P, S),如果,对 于P,满足(V∪T)+且中至少含 有一个非终结符,(V∪T)*,则G称为0型 文法(或短语结构文法,简记为PSG)或者 无约束文法(Unrestricted Grammar)。
• “不怕太阳晒,也不怕那风雨狂,只怕先生骂我 笨,没有学问无颜见爹娘 ……”
• “太阳当空照,花儿对我笑,小鸟说早早早……”
王雷@版权所有
1.1 符号、符号串及其运算
符号串的联结:联结是符号串的基本运算。两个符号串X和Y的 联结,记为XY,就是把Y跟随在X的后面形成的符号串。 例1.1:设 = {1, 2}是一个字母表。设X = 11、Y = 22分别是 上的两个符号串。则: XY = 1122是X、Y两个符号串的联结,XY是上的一符号串。 YX = 2211是Y、X两个符号串的联结,YX也是上的一符号串。
符号串的方幂:设X是符号 串,把X自身联结n次后,得 到的符号串Z,即Z = XX…XX = Xn,称为X的方幂。 我们约定X0 = 。这个定义 可以递归地表示为:

数值分析第一章PPT

数值分析第一章PPT

1.1.2 计算数学与科学计算 现代科学的三个组成部分: 科学理论, 科学实验, 科学计算 科学计算 的核心内容是以现代化的计算机及数学软件 (Matlab, Mathematica, Maple, MathCAD etc. )为工具,以数学 模型为基础进行模拟研究。
一些边缘学科的相继出现:
计算数学,计算物理学,计算力学,计算化学,计算生物学, 计算地质学,计算经济学,等等

取 0 e
1
x2
dx S4 ,
S4
R4
/* Remainder */
1 1 1 1 由留下部分 称为截断误差 /* Truncation Error */ 4! 9 5! 11 /* included terms */ 1 1 这里 R4 引起.005 0 由截去部分 4! 9 /* excluded terms */ 1 1 1 S4 1 1 0 .333 0 .1 0 .024 0 .743 引起 3 10 42 | 舍入误差 /* Roundoff Error */ | 0.0005 2 0.001
数值分析
第1章
数值分析与科学计算引论
§1.1 数值分析的对象、作用与特点
1.1.1 什么是数值分析 数值分析是计算数学的主要部分,计算数学是数学 科学的一个分支,它研究用计算机求解各种数学问题的 数值计算方法及其理论与软件实现.这门课程又称为(数 值)计算方法、科学与工程计算等。

在电子计算机成为数值计算的主要工具的今天, 需要研究适合计算机使用的数值计算方法。使用计 算机解决科学计算问题时大致要经历如下几个过程:
造成这种情况的是不稳定的算法 /* unstable algorithm */ 我们有责任改变。

计算机数学基础(二)教学辅导_01

计算机数学基础(二)教学辅导_01

计算机数学基础(二)教学辅导_01辅导日期:2002年8月28日(星期三)辅导提纲:1.《数值分析》课程主要内容介绍,重点掌握的内容概要。

2.教材内容:第9章数值分析中的误差(全部讲完)。

辅导内容一.《数值分析》与《离散数学》教材体系的不同特点1.《离散数学》是计算机数学基础上册的内容,教材体系的特点是概念多,定义多,知识点多。

但都比较好理解,各章之间的联系不是十分紧密。

学习的困难之处,不在于知识的深度,而在于知识的广度和对基本概念掌握的准确程度。

2.《数值分析》是计算机数学基础下册的主要内容,教材体系的特点是知识点集中,不是太多,但每个知识点的难度都较大,掌握起来都不是很轻松。

对同一个数值计算的问题,从不同的角度来解决,就产生了不同的计算方法。

这些方法总是和程序设计的流程相对应的。

3.学习策略的应变:加强对同一数值分析问题的各种不同的算法设计原理和计算方法的理解,注意它们相互之间的区别与联系。

有条件的话,最好能编程来实现之。

二.《数值分析》课程的主要知识点1.数值分析中的误差(第9章)(1)基本概念:绝对误差、绝对误差限;相对误差、相对误差限;近似数的有效数字。

(2)分析原理:数值计算中误差的传播及在数值计算中的若干准则。

(3)主要运算①给定近似值和有效数字的位数,求其相对误差限;②结定近似值和相对误差限,求它应当有几位有效数字。

2.线性方程组的数值解法(第10章)(1)高斯顺序消去法(2)列主元消去法(3)雅可比迭代法(4)高斯——赛德尔迭代法(5)迭代法的收敛性。

3.函数插值与最小二乘拟合(第11章)(1)从测试样点值决定函数曲线的基本原理(2)通过n+1个样点的拉格朗日n次插值多项式(3)通过n+1个样点的牛顿n次插值多项式(4)分段线性插值(5)三次样条插值(6)用最小二乘曲线拟合法从测试样点值求函数曲线的基本原理(7)直线拟合(8)多项式拟合(9)指数拟合。

28208272.doc 第1 页共2 页打印:203/28/201328208272.doc 第 2 页 共 2 页 打印:203/28/2013 4. 数值积分与微分(第12章)(1) 代数精度(2) 牛顿——科茨系数与牛顿——科茨公式(3) 高斯求积公式(4) 数值微分的二点求导公式与三点求导公式5. 方程求根(第13章)(1) 二分法求根(2) 迭代法求根:一般迭代法与牛顿迭代法(3) 弦截法求根6. 常微分方程的数值解法(第14章)(1) 欧拉法与改进的欧拉法(2) 二阶、三阶、四阶龙格——库塔法三. 第9章 数值分析中的误差(一) 知识点1. 基本概念: x ——近似值; *x ——精确值(1) 绝对误差:e = x - *x (注意:近似值在先,精确值在后。

数值分析ppt课件

数值分析ppt课件

数值积分与微分
数值积分
通过数值方法近似计算定积 分,如梯形法则、辛普森法 则等。
数值微分
通过数值方法近似计算函数 的导数,如差分法、中心差 分法等。
常微分方程的数值解法
通过数值方法求解常微分方 程,如欧拉方法、龙格-库塔 方法等。
03
数值分析的稳定性与误差分析
误差的来源与分类
模型误差
由于数学模型本身的近 似性和简化,与真实系
非线性代数方法
非线性方程组的求解
通过迭代法、直接法等求解非线性方程组,如牛顿法、拟牛顿法 等。
非线性最小二乘问题
通过迭代法、直接法等求解非线性最小二乘问题,如GaussNewton方法、Levenberg-Marquardt方法等。
多项式插值与逼近
通过多项式插值与逼近方法对函数进行近似,如拉格朗日插值、 样条插值等。
机器学习与数值分析的交叉研究
机器学习算法
利用数值分析方法优化和改进机器学 习模型的训练和预测过程,提高模型 的准确性和效率。
数据驱动的模型
通过数值分析方法处理大规模数据集 ,提取有用的特征和模式,为机器学 习模型提供更好的输入和输出。
大数据与数值分析的结合
大数据处理
利用数值分析方法处理和分析大规模数 据集,挖掘其中的规律、趋势和关联信 息。
数值分析PPT课件
contents
目录
• 引言 • 数值分析的基本方法 • 数值分析的稳定性与误差分析 • 数值分析的优化方法 • 数值分析的未来发展与挑战
01
引言
数值分析的定义
数值分析
数值分析是一门研究数值计算方法及 其应用的学科,旨在解决各种数学问 题,如微积分、线性代数、微分方程 等。

数值分析试题(卷)与答案解析

数值分析试题(卷)与答案解析

数值分析试题(卷)与答案解析数值分析试题一、填空题(2 0×2′)1.322A1, X23设 x=0.231 是精确值 x*=0.229 的近似值,则 x 有 2 位有效数字。

2. 若 f(x)=x7- x3+ 1 ,则 f[20 ,21,2 2,23 ,24,25,26,2 7]= 1 ,f[2 0,2 1,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(xi)的系数 ai(x)的特点是:a i ( x)1;所以i当系数 ai(x)满足ai(x)>1 ,计算时不会放大f(xi)的误差。

8.要使 20 的近似值的相对误差小于0.1% ,至少要取 4 位有效数字。

9.对任意初始向量(0) 及任意向量g ,线性方程组的迭代公式x(k+1)=Bx(k)+ (=0,1, ?)X g k收敛于方程组的精确解x*的充分必要条件是(B)<1 。

10.由下列数据所确定的插值多项式的次数最高是 5 。

x 0 0.5 1 1.5 2 2.5y=f(x) -2 -1.75 -1 0.25 2 4.2511.牛顿下山法的下山条件为|f(xn+1)|<|f(xn)| 。

12. 线性方程组的松弛迭代法是通过逐渐减少残差ri ( i=0,1, ? ,n)来实现的,其中的残差r = (b -a x-a x-? -a x )/aii, (i=0,1,,n)。

计算机数学基础(2)--误差分析(02-09)

计算机数学基础(2)--误差分析(02-09)

《计算机数学基础(2)》辅导第9章 数值分析中的误差 (2002级(秋季)用) 中央电大 冯 泰 《计算机数学基础》是中央广播电视大学开放本科教育计算机科学与技术专业教学中重要的核心基础课程,它是学习专业理论不可少的数学工具. 通过本课程的学习,要使学生具有现代数学的观点和方法,初步掌握处理离散结构所必须的描述工具和方法以及计算机上常用数值分析的构造思想和计算方法. 同时,也要培养学生抽象思维和慎密概括的能力,使学生具有良好的开拓专业理论的素质和使用所学知识,分析和解决实际问题的能力.本学期讲授数值分析部分,包括数值分析中的误差、线性方程组的数值解法、函数插值和最小二乘拟合、数值积分与微分、方程求根和常微分方程的数值解法. 通过本课程的学习,使学生熟悉数值计算方法的基本原理,掌握常见数值计算的方法. 依据教学大纲,我们对本学期的教学内容,逐章进行辅导,供师生学习参考.第9章 数值分析中的误差一、重点内容绝对误差-设精确值x *的近似值x , 差e =x -x *称为近似值x 的绝对误差(误差). 绝对误差限―绝对误差限ε是绝对误差e 绝对值的一个上界,即ε≤-=*x x e . 相对误差e r ―绝对误差e 与精确值x *的比值,***-==xx x xe e r .常用xe e r =计算.相对误差限r ε―相对误差e r 绝对值的一个上界,r r e ≥ε,常用xε计算.绝对误差限的估计式:)()()(2121x x x x εεε+=±)()()(122121x x x x x x εεε+≈22122121+=x x x x x x x )()()(εεε相对误差限的估计式:⎪⎪⎭⎫⎝⎛≠-±±≤±21212212121121)()()(x x x x x x x x x x x x x x r r r 时εεε112221)()()(x x x x x x r r r εεε+≤,221121)()()(x x x x x x r r r εεε+≤有效数字―如果近似值x 的绝对误差限ε是它某一个数位的半个单位,我们就说x 准确到该位. 从这一位起到前面第一个非0数字为止的所有数字称为x 的有效数字.关于有效数字的结论有: (1)设精确值x *的近似值x ,若mn a a a x 10.021⨯±=a 1,a 2,…,a n 是0~9之中的自然数,且a 1≠0,n l x x l m ≤≤110⨯50=≤--,.*ε 则x 有l 位有效数字.(2)设近似值mn a a a x 10.021⨯±= 有l 位有效数字,则其相对误差限111021+-⨯≤l r a ε(3) 设近似值m n a a a x 10.021⨯±= 的相对误差限不大于1110)1(21+-⨯+l a则它至少有l 位有效数字.(4) 要求精确到10-k(k 为正整数),则该数的近似值应保留k 位小数. 二、实例例1 设x *= π=3.1415926…,求x *的近似值及有效数字.解 若取x *的近似值x =3.14=0.314×101, 即m =1,它的绝对误差是-0.001 592 6…,有31105.06592001.0-*⨯≤=- x x ,即l =3,故近似值x =3.14有3位有效数字.或x =3.14的绝对误差限0.005,它是x *的小数后第2位的半个单位,故近似值x =3.14准确到小数点后第2位,有3位有效数字. 若取近似值x =3.1416,绝对误差是0.0000074…,有5-1*10⨯50≤00000740=-.. x x ,即m =1,l =5,故近似值x =3.1416有5位有效数字.或x =3.1416的绝对误差限0.00005,它是x *的小数后第4位的半个单位,故近似值x =3.1416准确到小数点后第4位,亦即有4位有效数字.若取近似值x =3.1415,绝对误差是0.0000926…,有 0000926.0=-*x x 41105.0-⨯≤,即m =1,l =4,故近似值x =3.1415只有4位有效数字.或x =3.1415的绝对误差限0.0005,它是x *的小数后第3位的半个单位,故近似值x =3.1415准确到小数点后第3位.注意:这就是说某数有s 位数,若末位数字是四舍五入得到的,那么该数有s 位有效数字.若末位数不是四舍五入得到的,那末它就不一定有s 位有效数字,必须用其绝对误差限来确定.绝对误差限是哪一位的半个单位,也就是精确到该位,从而确定有效数字. 例2 指出下列各数具有几位有效数字,及其绝对误差限和相对误差限: 2.000 4 -0.002 00 9 000 9 000.00解 因为x 1=2.000 4=0.200 04×101, 它的绝对误差限0.000 05=0.5×101―5,即m =1,l =5,故x =2.000 4有5位有效数字. a 1=2,相对误差限025000.01021511=⨯⨯=-a r εx 2=-0.002 00,绝对误差限0.000 005,因为m =-2,l =3,x 2=-0.002 00有3位有效数字. a 1=2,相对误差限εr =3110221-⨯⨯=0.002 5x 3=9 000,绝对误差限为0.5×100,因为m =4, l =4, x 3=9 000有4位有效数字,a =9,相对误差限εr =4110921-⨯=0.000 056x 4=9 000.00,绝对误差限0.005,因为m =4,l =6,x 4=9 000.00有6位有效数字,相对误差限为εr =6110921-⨯=0.000 000 56由x 3与x 4可以看到小数点之后的0,不是可有可无的,它是有实际意义的.例3 ln2=0.69314718…,精确到10-3的近似值是多少?解 精确到10-3=0.001,即绝对误差限是ε=0.0005, 故至少要保留小数点后三位才可以.ln2≈0.693例4 数值x *=2.197224577…的六位有效数字的近似值x =2.19722,而不是2.19723.注意:取一个数的近似数,若取5位有效数字,则只看该数第6位数,采取四舍五入的方法处理.与第7位,第8位的数值大小无关.本例取6位有效数字,左起第6个数是2,而第7个数是4,故应舍去,得到x=2.19722.本例第8个数,第9个数都是大于或等于5的数,再入上去,就得到x=2.19723,是不对的. 我们计算一下它们的误差. 取x=2.19722,e=x -x*=-0.000 004 577…,∣e ∣=∣x -x*∣=0.000 004 577…<0.000 005=0.5×101-6取x=2.19723,e=x -x*=0.000 005 423…,∣e ∣=∣x -x*∣=0.000 005 423…<0.000 05=0.5×101-5 即x=2.19723只有五位有效数字. 例5 设近似值x 1,x 2满足ε(x 1)=0.05,ε(x 2)=0.005,那么ε(x 1x 2)=?解 已知x 1,x 2的绝对误差限,求x 1x 2的绝对误差限.由绝对误差限的传播公式)()()(211221x x x x x x εεε+==1221005.005.0)(x x x x +=ε注:该传播公式也可以用于多个数的积, 213312321321)()()()(x x x x x x x x x x x x εεεε++=)(3)(),(2)(232x x x x x x εεεε==等.三、练习题1.下列各数中,绝对误差限为0.000 05的有效近似数是( B ) (A)-2.180 (B) 2.1200 (C) -123.000 (D) 2.120 2. 数8.000033的5位有效数字的近似值是多少? 答案:8.000 03. 若误差限为0.5×10-5,那么近似数0.003400有( B )位有效数字. (A) 2 (B) 3 (C) 4 (D) 64. 若近似值x 的绝对误差限为ε=0.5×10-2,那么以下有4位有效数字的x 值是( B ).(A) 0.934 4 (B) 9.344 (C) 93.44 (D)934.4 5. 已知准确值x *与其有t 位有效数字的近似值x =0.0a 1a 2…a n ×10s (a 1≠0)的绝对误差∣x *-x ∣≤( A ). (A) 0.5×10 s -1-t (B) 0.5×10 s -t (C) 0.5×10s +1-t (D) 0.5×10 s +t6. 已知x *1=x 1±0.5×10-3,x *2=x 2±0.5×10-2,那么近似值x 1,x 2之差的误差限是多少?答案:0.55×10-2. 7. 设近似值x =-9.73421的相对误差限是0.0005,则x 至少有几位有效数字. 答案:38. 用四舍五入的方法得到近似值x =0.0514,那么x 的绝对误差限和相对误差限各是几? 答案:0.000 05,0.0019. 设近似值x 1,x 2满足ε(x 1)=0.05,ε(x 2)=0.005,那么ε(x 1+x 2)=? 答案:0.05510. 设近似值x =±0.a 1a 2…a n ×10m ,具有l 位有效数字,则其相对误差限为( B ).(A) 1110121+-⨯+l a (B)1110)1(21+-⨯+l a(C)111021+-⨯l a (D) la -⨯1021111. 测量长度为x =10m 的正方形,若ε(x )=0.05m ,则该正方形的面积S 的绝对误差限是多少?答案:1(m)12.数值x*=2.197224577…的六位有效数字的近似值x=( B ).(A) 2.19723 (B) 2.19722 (C) 2.19720 (D) 2.19722513. 将下列各数舍入成三位有效数字,并确定近似值的绝对误差和相对误差. (1) 2.1514 (2) -392.85 (3) 0.00392214. 已知各近似值的相对误差,试确定其绝对误差:(1) 13267 e r=0.1% (2) 0.896 e r=10%四、练习题答案1. B2. 8.000 03. B4. B .5. A6. 0.55×10-2.7. 38. 0.000 05,0.0019. 0.05510. B11. 1(m)12. B13. (1)2.15, e=-0.001 4, e r=-0.000 65;(2) -393 , e=-0.15, e r=-0.00038;(3)0.00392, e=-0.000 002, e r=0.0005114. (1) e=0.13×102 (2) 0.9×10-1。

数值分析 (1)

数值分析 (1)

e * − e = (e * − en ) + (en − e )
2009-09-26 zhwang@ 17
2. 误差的度量
1) 2) 3) 4)
绝对误差 相对误差 有效数字 度量间的关系
2009-09-26
zhwang@
18
1)绝对误差
绝对误差定义:
zhwang@
22
相对误差(续2)
* * e ε ( x r 相对误差限: 的上界,记为 r ) 。 相对误差限:数值
相对误差限也可以通过
ε r* = ε * / x*
来计算。
Remark: 当要求计算相对误差,是指估计一个尽 可能小的相对误差限。 相对误差与相对误差限没有量纲。
分类方法1:若算法包含 有一个进程则称其为串行算法, 否则为并行算法。 分类方法2:从算法执行所 花费的时间角度来讲,若算术运 算占绝大多数时间则称其为数值 型算法,否则为非数值型算法。 本课程介绍数值型串行算 法。(其它类型算法参阅数据结 构、并行算法等课程)
2009-09-26
zhwang@
19
绝对误差(续)
•绝对误差限:
* * 如果存在正数 ε = ε(x ) ,使得有绝对误差
e * = x * − x ≤ ε* ,
则称 ε* 为 x*近似 x 的一个绝对误差限 绝对误差限。 绝对误差限
x ∈ [x * − ε * , x * + ε * ] , x = x * ± ε * 。
Remark: 通常计算中所要求的误差,是指 估计一个尽可能小的绝对误差限。 绝对误差与绝对误差限有量纲。
2009-09-26 zhwang@ 10
算法应用状态
数值分析研究对象以及解决问题方法的 广泛适用性,著名流行软件如Maple、Matlab、 Mathematica等已将其绝大多数内容设计成函 数,简单调用之后便可以得到运行结果。 但由于实际问题的具体特征、复杂性, 以及算法自身的适用范围决定了应用中必须选 择、设计适合于自己特定问题的算法,因而掌 握数值方法的思想和内容是至关重要的。

《计算机数学基础》数值部分第二单元辅导

《计算机数学基础》数值部分第二单元辅导
其中基函数 , 。
当n=2时,得到二次多项式,就是二次插值。
拉格朗日插值多项式的余项为
其中,
注意:过n+1个互异节点,所得插值多项式应该是次数不超过n的多项式。
3. 均差与牛顿插值多项式
函数值之差与自变量之差的商就是均差,
一阶均差
二阶均差
... ... ... ...
k xk yk xkyk 1 1 4 1 4 2 2 4.5 4 9 3 3 6 9 18 4 4 8 16 32 5 5 8.5 25 42.5 ? 15 31 55 105.5
解得a0=2.45, a1=1.25。所求拟合直线方程为 y=2.45+1.25x
(*)
其中: ,
(k=1,2,...,n-1)
附加条件:(1) 当已知S?(x0)=y?0 ,S?(xn)=y?n时,(*)式中
最小,求?(x)的方法,称为最小二乘法。
(1) 直线拟合 若,a0,a1满足法方程组
即a0, a1是法方程组的解。
(2) 二次多项式拟合 若满足法方程组
即a0, a1,a2是法方程组的解。
二、实例
例1 已知函数y=f(x)的观察数据为
上式正是Qn(x)的拉格朗日插值多项式。可见,Qn(x)的拉格朗日插值多项式就是它自身,即次数不超过n的多项式在n+1个互异节点处的拉格朗日插值多项式就是它自身。
例4 已知函数e-x的下
列数据用分段线性插值法 x 0.10 0.15 0.25 0.30 求x=0.2的近似值。 e-x 0.904 837 0.860 708 0.778 805 0.740 818 解 用分段线性插值,先求基函数。

(完整版)数值分析每节课的教学重点、难点

(完整版)数值分析每节课的教学重点、难点

计算方法教案新疆医科大学数学教研室张利萍一、课程基本信息1、课程英文名称:Numerical Analysis2、课程类别:专业基础课程3、课程学时:总学时544、学分:45、先修课程:《高等数学》、《线性代数》、《Matlab 语言》二、课程的目的与任务:计算方法是信息管理与信息系统专业的重要理论基础课程,是现代数学的一个重要分支。

其主要任务是介绍进行科学计算的理论方法,即在计算机上对来自科学研究和工程实际中的数学问题进行数值计算和分析的理论和方法。

通过本课程的学习,不仅使学生初步掌握数值分析的基本理论知识,而且使学生具备一定的科学计算的能力、分析问题和解决问题的能力,为学习后继课程以及将来从事科学计算、计算机应用和科学研究等工作奠定必要的数学基础。

三、课程的基本要求:1.掌握计算方法的常用的基本的数值计算方法2.掌握计算方法的基本理论、分析方法和原理3.能利用计算机解决科学和工程中的某些数值计算应用问题,增强学生综合运用知识的能力4.了解科学计算的发展方向和应用前景四、教学内容、要求及学时分配:(一) 理论教学:引论(2学时)第一讲(1-2节)1.教学内容:计算方法(数值分析)这门课程的形成背景及主要研究内容、研究方法、主要特点;算法的有关概念及要求;误差的来源、意义、及其有关概念。

数值计算中应注意的一些问题。

2.重点难点:算法设计及其表达法;误差的基本概念。

数值计算中应注意的一些问题。

3.教学目标:了解数值分析的基本概念;掌握误差的基本概念:误差、相对误差、误差限、相对误差限、有效数字;理解有效数字与误差的关系。

学会选用相对较好的数值计算方法。

A 算法B 误差第二讲典型例题第二章线性方程组的直接法(4学时)第三讲1.教学内容:线性方程组的消去法、Gauss消去法及其Gauss列主元素消去法的计算过程;三种消去法的程序设计。

2.重点难点:约当消去法,Gauss消去法,Gauss列主元素消去法3.教学目标:了解线性方程组的解法;掌握约当消去法、Gauss消去法、Gauss列主元素消去的基本思想;能利用这三种消去法对线性方程组进行求解,并编制相应的应用程序。

计算机数学基础-数值分析与组合数学第二版下册教学设计

计算机数学基础-数值分析与组合数学第二版下册教学设计

计算机数学基础-数值分析与组合数学第二版下册教学设计一、教学目标本教材是计算机数学基础的一部分,主要介绍数值分析与组合数学的相关知识。

通过本课程的学习,学生应该能够:1.掌握数值分析的基本原理和方法;2.熟练应用常用数值分析算法,如牛顿迭代法、二分法等;3.熟悉组合数学的基本概念和方法;4.掌握排列、组合、多重集合等概念,并能够应用到实际问题中。

二、教学内容与进度安排1. 数值分析第1章引言第2章非线性方程的数值解法第3章线性方程组的数值解法第4章差分法第5章插值法第6章数值微积分第7章数值积分第8章常微分方程的初值问题2. 组合数学第1章引言第2章排列、组合与离散概率第3章二项式系数第4章多重集合第5章递推关系与生成函数第6章偏序与格3. 实践环节在课程中,我们将通过各种实例和案例让学生更好地理解和掌握所学知识。

同时,实践环节也是促进学生思考和创新的一个重要机会。

三、教学方法与评价方法1. 教学方法本课程主要采用“讲授+实践”的教学模式。

具体包括:1.讲授:由教师进行讲述和解释;2.实践:通过项目、应用案例等实践环节,让学生积累实际经验,加深自己的理解;3.评价:通过每次课程作业和期末考试等多个环节对学生进行评价和反馈。

2. 评价方法对于学生的学习成果,我们将主要采用以下方法进行评价:1.课程作业:每周布置一定量的作业,检验学生掌握情况;2.课堂测试:随时进行小测验,帮助学生巩固所学知识;3.期末考试:检验学生对整个课程的综合掌握情况。

四、参考书目1.《数值分析》(第2版),北京:高等教育出版社,2014年;2.《组合数学》(第2版),北京:高等教育出版社,2011年;3.《计算机数学基础》(第3版),北京:清华大学出版社,2015年。

五、教学团队本课程由清华大学计算机科学与技术系教师组成的团队授课,教师有多年数值分析和组合数学的教学和研究经验,充分保证了教学质量和水平。

计算机数学基础下第5编数值分析PPT课件

计算机数学基础下第5编数值分析PPT课件

例1
考察求积公式 1 f (x)dx 1[ f (1) 2 f (0) f (1)]
1
2
具有几次代数精度。
解:考察 f (x) 为零次多项式的情况。例如 f (x) 1
1
左边 1dx 2 1
公式成立,
右边 1 [1 2 1] 2 2
再考察 f (x)为一次多项式时的情况,取 f (x) x
C (n) k
(1)nk n k!(n k)!
n 0
n
(t
j0
j)dt
jk
Ck(n) 称为科茨系数
牛顿 科茨公式
则Ak (b a)Ck(n)
b
f (x)dx (b a)
a
n
C (n) k
f
(xk
)
k 0
该公式称为n阶牛顿—科茨公式
是数值积分的基 本公式之一, 应当记住
科茨系数
Ck(
Ak
n 0
t(t k(k
11))((tknn))hdt
(分子中没有 t k ,分母中没有 k k )
考察分母
k(k 1)(k 2)1 (1)(2)[(n k)] k!(1)nk (n k)!
Ak
(1)nk b a k!(n k)! n
n 0
n
(t
j0
j)dt
jk
若记
n)具有下列性质:
n
性质(1)
C (n) k
1
(归一性)
k 0
性质(2)
C (n) k
C (n) nk
(对称性)
常见的科茨系数表见 P.119
需记住的有:
当n=1时,C0(1)
1 2
,

《计算机数学基础》(第二版)习题参考答案

《计算机数学基础》(第二版)习题参考答案

《计算机数学基础》(第二版)习题参考答案习题1.11.42,23,42---x x ,1722++x x ,4682-+x x ,h x 234++。

2. (1)]14,6[,]3,2[-=-=R D 。

(2)]1,0[,]1,1[=-=R D 。

(3)),0[,),(∞+=∞+-∞=R D 。

(4)),0[,),(∞+=∞+-∞=R D 。

(5)]1,1[,),(-=∞+-∞=R D 。

3.(1)不同,因为定义域不同。

(2)不同,因为对应规则不同。

(3)相同,因为定义域和对应规则均相同。

4.(1)]2,2[-=D 。

(2)}1|{≠=x x D 。

(3)),(D ∞+-∞=。

(4)),(D ∞+-∞=。

图略5.(1)2010h T +-=。

(2)10k =。

(3)C 5︒-。

6.(1)有界;(2)有界;(3)无界;(4)有界。

7.(1)非奇非偶函数;(2)奇函数;(3)偶函数;(4)偶函数。

8.(1)周期函数,周期是π2;(2)非周期函数;(3)周期函数,周期是π。

习题1.21.(1)),(,)13(2))((223∞+-∞=-±+=±±g f D x x x x g f ; ),(,263))((2345∞+-∞=--+=∙fg D x x x x x g f ;}33|{,132))(/(/223±≠=-+=x x D x x x x g f g f 。

(2)]1,1[,11))((-=-±+=±±g f D x x x g f ;]1,1[,1))((2-=-=∙fg D x x g f ;)1,1[,11))(/(/-=-+=g f D xx x g f 。

2.(1)),(,62118))((2∞+-∞=++=g f D x x x g f , ),(,236))((2∞+-∞=+-=f g D x x x f g , ),(,88))((34∞+-∞=+-=f f D x x x x f f ,),(,89))((∞+-∞=+=g g D x x g g 。

数值分析全套课件

数值分析全套课件

Ln n si n

ˆ L2n (4L2n Ln ) / 3
n L error 192 3.1414524 1.4e-004 384 3.1415576 3.5e-005 3.1415926 4.6e-010
3/16
通信卫星覆盖地球面积
将地球考虑成一 个球体, 设R为地 球半径,h为卫星 高度,D为覆盖面 在切痕平面上的 投影(积分区域)
( x1 x2 ) | x1 | ( x2 ) | x2 | ( x1 )
15/16
例3.二次方程 x2 – 16 x + 1 = 0, 取
求 x1 8 63 使具有4位有效数
63 7.937
解:直接计算 x1≈8 – 7.937 = 0.063
( x1 ) (8) (7.937) 0.0005
5/16
误差的有关概念
假设某一数据的准确值为 x*,其近似值 为 x,则称
e(x)= x - x*
为 x 的绝对误差 而称
e( x) x x er ( x ) , x x
*
( x 0)

为 x 的相对误差
6/16
如果存在一个适当小的正数ε

,使得
e( x) x x
计算出的x1 具有两位有效数
1 0.062747 修改算法 x1 8 63 15.937 4位有效数 (15.937) 0.0005 ( x1 ) 0.000005 2 2 (15.937) (15.937)
16/16
1
参考文献
[1]李庆扬 关治 白峰杉, 数值计算原理(清华) [2]蔡大用 白峰杉, 现代科学计算 [3]蔡大用, 数值分析与实验学习指导 [4]孙志忠,计算方法典型例题分析 [5]车刚明等, 数值分析典型题解析(西北工大) [6]David Kincaid,数值分析(第三版) [7] John H. Mathews,数值方法(MATLAB版)

《计算机数学基础(下)》期未复习辅导.doc

《计算机数学基础(下)》期未复习辅导.doc

《计算机数学基础》数值分析期末复习提纲中央电大数理教研室《计算机数学基础》数值分析部分是中央广播电视大学本科开放教育计算机科学与技术专业学生必修的一门专业基础课程,使用教材是任现淼主编、吴裕树副主编的《计算机数学基础(下册7值分析与组合数学》(中央电大出版社出版)。

期末考试全国统一命题。

一、期末考试试题期末考试的试卷有单项选择题、填空题和解答题。

单项选择题和填空题各5个题,分数约占30%。

解答题共5个题,包括计算题、化简题和证明题等,分数约占70%。

各章分数的分布为第9章约6分,第10〜14各章有选择题、填空题和解答题,分数分配大致与所用课时成比例。

期末考试的内容和要求以中央电大编发的《计算机数学基础(下)数值分析部分考核说明》为准。

主要考核基本概念、基本原理和基本运算。

可以带简易计算器。

二、考核知识点、要求、例题与参考练习题以下分章给出期末考试的考核知识点、复习要求、例题与参考练习题,供期末复习和考试参考。

第9章数值分析中的误差(一)考核知识点误差的来源类型;绝对误差和绝对误差限,相对误差和相对误差限,有效数字;绝对误差的传播。

(二)复习要求1.知道产生误差的主要来源。

2.了解绝对误差和绝对误差限、相对误差和相对误差限和有效数字等概念以及它们之间的关系。

3.知道四则运算中的误差传播公式。

(三)例题例1 指出下列各数具有几位有效数字,及其绝对误差限和相对误差限:2.000 4 一0.002 00 9 000.00解因为由=2.000 4= 0.200 04X 101,绝对误差限0.000 05=0.5X 10 1-5,即m=l,/=5,故x=2.000 4有5位有效数字.相对误差限& = —x 10* = 0.000 0252x2x2= - 0.002 00 ,绝对误差限0.000 005 , 3位有效数字。

相对误差限沪一*—x lO-3+i = 0.002 52x2X3=9 000.00,绝对误差限0.005, 6位有效数字,相对误差限为£户—^― x 10一6袒=0.000 0005 62x9例2所2=0.69314718...,精确到1。

《计算机数学基础》辅导(3)

《计算机数学基础》辅导(3)

《计算机数学基础》辅导(3)−−集合及其运算本章重点:集合概念,集合的运算,集合恒等式的证明,笛卡儿积.一、重点内容1. 集合的概念集合与元素,具有确定的,可以区分的若干事物的全体称为集合,其中的事物叫元素.集合A中所含元素的个数记作A.集合中的元素不能重复出现,集合中的元素无序之分. 集合与其元素之间有属于“∈”或不属于“∉”之分.集合的表示方法:列举法和描述法.2. 特殊集合:全集、空集和幂集全集合E,在一个具体问题中,所涉及的集合都是某个集合的子集,该集合为全集;空集∅,不含任何元素的集合为空集. 空集是惟一的,它是任何集合的子集.集合A的幂集P(A),有集合A的所有子集构成的集合. 若∣A∣=n, 则∣P(A)∣=2n.3. 集合的关系:包含,子集,集合相等.包含(子集),若B∈⇒∀,,则B包含A(或A包含于B),A是B的aa∈A子集,记BA⊆又A≠B,则A是B的真子集,记A⊂B.集合相等,若A⊆B,B⊆A,则A=B.注意:在集合概念部分要特别注意:元素与集合,集合与子集,子集与幂集,∈与⊂(⊆),空集∅与所有集合等的关系.4. 集合的运算集合A和B的并,由集合A和B的所有元素组成的集合,A⋃B集合A和B的交,由集合A和B的公共元素组成的集合,A⋂B集合A的补集~A,由不属于集合A的元素组成的集合,~A. 补集总相对于一个全集.集合A与B的差集,由属于A,而不属于B的所有元素组成的集合,A-B.集合A与B的对称差,A⊕B=(A-B)⋃(B-A),也有A⊕B=)A⋃B〕-(A⋂B)应该很好地掌握10条运算律(运算的性质),即交换律、结合律、分配律、幂等律、同一律、零律、补余律、吸收律、摩根律和双补律等.5. 恒等式证明集合的运算部分有三个方面的问题:其一是进行集合的运算;其二是集合运算式的化简;其三是集合恒等式的推理证明.集合恒等式的证明方法通常有二:其一,要证明A=B,只需要证明A⊆B,又A⊇B;其二,通过运算律进行等式推导.6. 有序对与笛卡儿积有序对,就是有顺序的数组,如<x,y>,x,y的位置是确定的,不能随意放置.注意:有序对<a,b>≠<b,a>,以a,b为元素的集合{a,b}={b,a};有序对(a,a)有意义,而集合{a,a}不成立,因为它只是单元素集合,应记作{a}.笛卡儿积,是一种集合合成的方法,把集合A,B合成集合A×B,规定A×B={<x,y>∣x∈A,y∈B}由于有序对<x,y>中x,y的位置是确定的,因此A×B的记法也是确定的,不能写成B×A.笛卡儿积也可以多个集合合成,A1×A2×…×A n.笛卡儿积的运算性质.二、实例例3.1已知S={2,a,{3},4},R={{a},3,4,1},指出下列命题的真值.(1) {a}∈S; (2) {a}∈R;(3) {a,4,{3}}⊆S; (4) {{a},1,3,4}⊆R;(5) R=S; (6) {a}⊆S(7) {a}⊆R (8) ∅⊂R(9) ∅⊆{{a}}⊆R (10) {∅}⊆S(11) ∅∈R (12) ∅⊆{{3},4}解集合S有四个元素组成:2,a,{3},4,而元素{3}又是集合. 集合R 类似.(1) {a},这是单元素的集合,{a}不是集合S的元素. 故命题A:{a}∈S的真值为0;(2) {a}是R的元素,故命题B:{a}∈R的真值为1.(3) a,4,{3}都是集合S的元素,它们可以构成S的子集. 故命题C:{a,4,{3}}⊆S的真值为1(4) {a},1,3,4都是R的元素,它们可以构成R的子集,故命题D:{{a},1,3,4}⊆R的真值为1.(6)和(8),(9)和(12)相应题号的命题,其真值为1;而(5),(7),(10)相应题号的命题,其真值为0.例3.2设A={=,∈,∉,⊂, ⊃}选择适当的符号填在各小题的横线上.(1) (1,2,3,4) N ; (2) Z Q Q ,2 (3) },056{}5,1{2R x x x x ∈=+-∅ (4) },3{},2{22R y y R x x ∈<∈< (5) }},{{}{a a a(6) {正方形} {菱形} {四边形}(7) {(1,2,3)} {1,2,3,{(1,2,3)}}解 (1) ⊂ (2) ∉, ⊃ (3) ⊂ , = (4) ⊂ (5) ∈或⊂(6) ⊂ ⊂ (7) ∈例3.3 写出下列集合的子集:(1) A ={a ,{b },c }(2) B ={∅}(3) C =∅解 (1)因为∅是任何集合的子集,所以∅是集合A 的子集;由A 的任何一个元素构成的集合,都是A 的子集,所以{a },{{b }},{c }是A 的子集;由A 的任何两个元素构成的集合,都是A 的子集,所以{a ,{b }},{{b },{c }},{a , c }是A 的子集;由A 的任何三个元素构成的集合,也是A 的子集,所以{a ,{b },c }=A 是A 的子集;于是集合A 的所有子集为∅,{a },{{b }},{c },{a ,{b }},{{b },{c }},{a , c },{a ,{b },c }=A(2) 同(1),B 的子集有:∅,{∅}.(3) 因为∅是任何集合的子集,故∅也是C 的子集. 因为C 中没有元素,因此C 就没有其它子集,所以C 的子集只有:∅.说明:(1) 以集合A 的8个子集为元素的集合,就是集合A 的幂集,即P (A )={ ∅,{a },{{b }},{c },{a ,{b }},{{b },{c }},{a , c },{a ,{b },c }}那么集合B 的幂集为;P (B )={∅,{∅}};集合C 的幂集:P (C )={∅}.一般地,如果集合A ,有,n A =那么P (A )有2n 个元素.(2) 根据真子集的定义,对于任何集合A ,除了集合A 本身不是A 的真子集外,其它子集均是A 的真子集. 于是本例集合A 有7个真子集:∅,{a },{{b }},{c },{a ,{b }},{{b },{c }},{a , c }集合B 只有1个真子集:∅集合C 没有真子集.例3.4设集合A ={1,2,3,4},B ={2,3,5},求B A A B B A B A B A ⊕--⋂⋃,,,,. 解 };5,4,3,2,1{=⋃B A}5,4,1{}5{}4,1{)()(}5{};4,1{};3,2{=⋃=-⋃-=⊕=-=-=⋂A B B A B A A B B A B A例3.5 试证A -(B -C )=(A -B )⋃(A ⋂C ) 证明 [方法1] 对任意x ,)()()()()()()()()~()(C A B A C B A C A B A x C x A x B x A x C x B x A x C B x A x C B A x ⋂⋃-⊆--∴⋂⋃-∈⇒∈∧∈∨∉∧∈⇒∈∨∉∧∈⇒⋂∉∧∈⇒--∈ 同理,有)()()(C B A x C A B A x --∈⇒⋂⋃-∈所以,A -(B -C )=(A -B )⋃(A ⋂C )说明:事实上,方法1的证明,完全是等值过程,可以写作)()()()()()~()(C A B A x C x A x B x A x C x B x A x C B x A x C B A x ⋂⋃-∈⇔∈∧∈∨∉∧∈⇔∈∨∉∧∈⇔⋂∉∧∈⇔--∈[方法2] 进行恒等推导.A -(B -C )=)~(~C B A ⋂⋂(分配律)摩根律)()~()()(~C A B A C B A ⋂⋃⋂=⋃⋂==(A -B )⋃(A ⋂C )例3.6 化简))(()))(((A B B A C B A --⋃⋂-⋃解 ))(()))(((A B B A C B A --⋃⋂-⋃=))~())(((A B B C B E A ⋃⋃⋃-⋃⋂)=A A E A =⋃⋂(例3.7 设集合 A ={a ,b },B ={1,2,3},C ={d },求A ×B ×C ,B ×A.解 先计算A ×B ={<a ,1>,<a ,2>,<a ,3>,<b ,1>,<b ,2>,<b ,3>}A ×B ×C ={<a ,1>,<a ,2>,<a ,3>,<b ,1>,<b ,2>,<b ,3>}×{d }={<<a ,1>,d >,<<a ,2>,d >,<<a ,3>,d >,<<b ,1>,d >,<<b ,2>,d >,<<b ,3>,d >}B ×A ={<1,a >,<2,a >,<3,a >,<1,b >,<2,b >,<3,b >}例3.8 设集合A ={1,2},求A ×P (A ).解P(A)={∅,{1},{2},{1,2}}A×P(A)={1,2}×{∅,{1},}{2},{1,2}={<1,∅>,<2,∅>,<1,{1}>,<2,{1}>,<1,{2}>,<2,{2}>,<1,{1,2}>,<2,{1,2}>} 例3.9 单项选择题1. 若集合A={a,b,c},∅为空集合,则下列表示正确的是( )(A) {a}∈A(B){a}⊂A(C) a⊂A(D) ∅∈A答案:(B)解答:由集合A的元素构成的集合是A的子集,{a}是A的子集,故选择(B)正确.2.对任意集合S,S⋃∅=S,满足( )(A) 幂等律 (B) 零一律 (C) 同一律 (D) 互补律答案:{C}解答:见集合的运算性质,A⋃∅=A和E⋂A=A称为同一律.例3.10 填空题1设全集合E={1,2,3,4,5},A={1,2,3},B={2,5},A⋂B= ,~B= .~A⋃~B=答案:{2},{1,3,4},{1,3,4,5}解答:A⋂B是由集合A,B的公共元素构成的新集合. 此处A,B公共元素只有2,故A⋂B={2},~B是全集合中除去B的元素所剩余元素构成的新集合,全集合E有1,2,3,4,5,除去B的元素2,5,余下有1,3,4. 故~B={1,3,4}. ~A={4,5},于是~A⋃~B={1,3,4,5}2.设集合A={a,b,c},B={a,b},那么P(A)-P(B)= P(B)-P(A)=答案:{c},{a,c},{b,c},{a,b,c};∅解答:P(A)={∅,{a},{b},{c},{a,b},{b,c},{a,c},{a,b,c}}P(B)={∅,{a},{b},{a,b}}所以P(A)-P(B)={ {c},{a,c},{b,c},{a,b,c}}.∵P(A) ⊂P(B),∴P(B)-P(A)=∅三、练习题1.设S,T,M为任意集合,判定下列命题的真假:(1)∅是∅的子集;(2)如果S⋃T=S⋃M,则T=M;(3)如果S-T=∅,则S=T;(4)如果~S⋃T=E,则S⊆T;(5) S⊕S=S2. 用列举法表示以下集合:(1) }7xxAN{2≤∧∈=x(2) }3Nx=xxA-3{<∧∈(3) }0)1({2≤+∧∈=x R x x A3. 求使得下列集合等式成立时,a , b , c , d 应该满足的条件:(1) {a , b }={a , b , c }(2) {a , b , a }={a ,b }(3) {{a , ∅}, b , {c }}={{∅}}4. 求幂集P (A ),设集合A 为(1) A ={{1, 1 }, {2, 1 },{1, 2, 1} };(2) A =P (A )5. 设A ,B 为任意集合,试证明B A A B B A =⇔-=-6 设集合A ={1,2,{1,2},∅}, 试求:(1) A -{1,2};(2) A -∅;(3)A -{∅};(4){{1,2}}-A ;(5)∅-A ;(6) {∅}-A7. 试证对任意集合A ,B ,C ,等式(A -B )⋃(A -C )=A 成立的充分必要条件是A ⋂B ⋂C =∅四、练习题答案1. (1),(4)为真,其余为假.2. (1) A ={0,1,2}(2) A ={1,2,3,4,5}(3) A ={-1}3. (1) a =c 或c =b(2) 任意a , b(3) a =c =∅,且b ={∅}4. (1) P (A )={∅, {{1}}, {{1,2}}, {{1}, {1, 2}}}先将集合A 化简为{{1},{1,2}},再求幂集.(2) P (A )={∅, {∅}, {{1}}, {{2}}, {{1, 2}},{∅,{1}}, {∅, {2}}, {∅, {1,2}}, {{1}, {2}}, {{1}, {1, 2}}, {{2}, {1,2}}, {∅, {1}, {2} }, {∅, {1},{1,2}}, {∅, {2}, {1,2 }}, {{1}, {2}, {1, 2 }}, {∅, {1}, {2}, {1,2}}}先求P (A ),再求幂集.5. 当A =B 时,必有A -B =B -A ;反之,由A -B =B -A ,得到B A B B B A ⋂-=⋂-)()(化简后得到∅=-A B ,即A B ⊆;同理,由A -B =B -A ,得到A AB A B A ⋂-=⋂-)()(化简后得到∅=-B A ,即B A ⊆.A =B6.(1) {{1,2},∅}. (2) A ; (3) {1,2,{1,2}};(4) ∅; (5) ∅; (6)∅提示:(1)此处{1,2}是以1,2为元素的A 的子集. 属于A ,而不属于{1,2}的元素有{1,2}和∅,故A -{1,2}={{1,2},∅}.此处把{1,2}理解为A 的元素,所求集合A 减去一个元素是无意义的. 也就是说,集合之间可以进行并、交、补、差等运算,一个集合与一个元素之间不能进行运算.(2) 此处的∅是空集合,不能理解为集合A 的元素. 从集合A 减去一个没有元素的集合,结果还是A.注意:A 中有元素∅,如果理解为元素∅,也就出现了集合减元素的错误.(3) 此处{∅}是A 的子集,结果为从A 中除去元素∅,为{1,2,{1,2}}(4) 集合{1,2}是集合A 的以1,2为元素的子集,属于{1,2}而不属于A 是不可能的,故其结果为∅.(5) 属于空集合∅而不属于A 这是不可能的,故结果为∅.(6)以A 的元素∅为元素的A 的子集{∅}减去A ,结果为∅.7. 必要性设(A -B )⋃(A -C )=A ,因为(A -B )⋃(A -C )=)~()~(C A B A ⋂⋂)()(~)~(~(C B A C B A C B A -==⋂= 所以 A C B A =-)(于是对于任意,A x ∈必有)(C B A x -∈,而必有C B x ∉,故有ΦC B A =)(充分性设ΦC B A =)( ,则对于任意A x ∈,必有C B x ∉,即)(~C B x ∈,因此)~C B A ⊆ 于是,AC B A C B A C A B A C A B A =⋂===--)(~)~(~)~()~()()(。

《计算机数学基础 下》数值分析部分辅... - 中华残疾

《计算机数学基础 下》数值分析部分辅... - 中华残疾

《计算机数学基础(下)》数值部分辅导(5)中央电大 冯泰第13章 方程求根一、重点内容1. 二分法: 设方程f (x )=0在区间[a ,b ]内有根,用二分有根区间的方法,得到有根区间序列:[a ,b ]⊇ [a 1,b 1] ⊇ [a 2,b 2] ⊇---⊇ [a n ,b n ] ⊇…x *≈x n =)(21n n b a +(a 0=a ,b 0=b ),n =0,1,2,… 有误差估计式: ∣x *-x n ∣≤1+2-n a b ,n =0,1,2,… 二分法区间次数: 2--≥1+ln ln )ln(εa b n 2. 简单迭代法: 若方程f (x )=0表成x =ϕ(x ),于是有迭代格式:x n =ϕ(x n -1) (n=1,2,…))()lim (lim *1*x x x x n n n n ϕϕ===-∞→∞→ x *≈x n若存在0<<1,(x ),在区间[a ,b ]内任一点为初始值进行迭代,迭代数列收敛.快速迭代法: ⎪⎪⎩⎪⎪⎨⎧ 3. 牛顿法:用切线与x 轴的交点,逼近曲线f (x )与x 轴的交点.迭代公式为 )()(111---'-=n n n n x f x f x x (n=1,2,…) 选初始值x 0满足f (x 0)f (x 0)>0,迭代解数列一定收敛.4. 弦截法: 用两点连线与x 轴交点逼近曲线f (x )与x 轴的交点.迭代公式为)()()()(1-1-1+---=n n n n n n n x x x f x f x f x x (n =1,2,…) 二、实例例1 证明方程1-x -sin x =0在区间[0,1]内有一个根,使用二分法求误差不超过0.5×10-4的根要迭代多少次?证明 令f (x )=1-x -sin x ,∵ f (0)=1>0,f (1)=-sin1<0∴ f (x )=1-x -sin x =0在[0,1]有根.又f (x )=1-c os x >0(x [0.1]),故f (x )=0在区间[0,1]内有唯一实根.给定误差限=0.5×10-4,有728713=1-2104+50-=1-2--≥.ln ln .ln ln ln )ln(εa b n只要取n =14.例2 用迭代法求方程x 5-4x -2=0的最小正根.计算过程保留4位小数. [分析] 容易判断[1,2]是方程的有根区间.若建立迭代格式校正值 )(~1n n x x ϕ=+ 再校正值 )~(11++=n n x x ϕ 改进值 n n n n n n n xx x x x x x +---=++++++1121111~2)~()),(()(,)(,21∈1>45='42-=42-=454x x x x x x x ϕϕ即,此时迭代发散. 建立迭代格式)),(()(,)(,21∈54<2+454='2+4=2+4=455x x x x x x x ϕϕ,此时迭代收敛.解 建立迭代格式552+4=2+4=x x x x )(,ϕ1=21∈54<2+454='04x x x x 取初始值)),,(()(ϕ ≈6=2+4=5501x x 1.431 0 ≈7247=2+4=5512.x x 1.505 1 ≈02048=2+4=5523.x x 1.516 5 ≈0668=2+4=5534.x x 1.518 2≈07288=2+4=5545.x x 1.5185取≈*x 1.5185例3 试建立计算3a 的牛顿迭代格式,并求3791411.的近似值,要求迭代误差不超过10-6.[分析]首先建立迭代格式.确定取几位小数,求到两个近似解之差的绝对值不超过10-6.解 令0=-==33a x x f a x )(,,求x 的值.牛顿迭代格式为),...,,()()(10=3+32=3--='-=2231+k x a x x a x x x f x f x x k k k k k k k k k迭代误差不超过10-6,计算结果应保留小数点后6位.当x =7或8时,x 3=343或512,0>8''80<7''7)()(,)()(f f f f 而,取x 0=8,有 ≈8⨯3791411+8⨯32=3+32=22001.x a x x 7.478 078 ≈4780787⨯3791411+4780787⨯32=3+32=22212...x a x x 7.439 9560381220=-21.x x≈4399567⨯3791411+4399567⨯32=3+32=22223...x a x x 7.439760 0001960=-32.x x≈4397607⨯3791411+4397607⨯32=3+32=22334...x a x x 7.439760 于是,取≈*x 7.439760例4 用弦截法求方程x 3-x 2-1=0,在x[分析] 先确定有根区间.再代公式.解 f (x )= x 3-x 2-1,f (1)=-1,f (2)=3,有根区间取[1,2]. 取x 1=1, 迭代公式为 )()()()(1-1-1+---=n n n n n n n x x x f x f x f x x (n =1,2,…)251≈1⨯43-2=-+--1---=0120302131213112.)(x x x x x x x x x x ≈2-251⨯2+2-251-2511-251-251-251=2323233).(.....x 1.37662 ≈251-376621⨯251+251-376621-3766211-376621-376621-376621=2323234)..(.......x 148881 ≈376621-488811⨯376621+376621-488811-4888111-488811-488811-488811=2323235)..(.......x 146348 ≈488811-463481⨯488811+488811-463481-4634811-463481-463481-463481=2323236)..(.......x 1.46553 取≈*x 1.46553,f (1.46553)-0.000145例4 选择填空题1. 设函数f (x )在区间[a ,b ]上连续,若满足 ,则方程f (x )=0在区间[a ,b ]一定有实根.答案:f (a )f (b )<0解答:因为f (x )在区间[a ,b ]上连续,在两端点函数值异号,由连续函数的介值定理,必存在c ,使得f (c )=0,故f (x )=0一定有根.2. 用简单迭代法求方程f (x )=0的实根,把方程(x )=0表成x =(x ),则f (x )=0的根是( )(A)y =x 与y =(x )的交点 (B) y =x 与y =(x )交点的横坐标(C) y =x 与x 轴的交点的横坐标 (D) y =(x )与x 轴交点的横坐标答案:(B)解答:把f (x )=0表成x =(x ), 满足x =(x )的x 是方程的解,它正是y =x 与y =(x )的交点的横坐标.3.为求方程x 3―x 2―1=0在区间[1.3,1.6]内的一个根,把方程改写成下列形式,并建立相应的迭代公式,迭代公式不收敛的是( )(A) 1-1=1-1=1+2k k x x x x :,迭代公式 (B) 21+21+1=1+1=k k x x x x :,迭代公式(C) 3121+23+1=+1=/)(:,k k x x x x 迭代公式(D) 1+++1==1-221+23k k k k x x x x x x :,迭代公式答案:(A) 解答:在(A)中,.).()()(,)(,//0758281=1-6121>1-21-='1-1=1-1=23232x x x x x x ϕϕ 故迭代发散.在(B)中1<9010=311<2-='1+1=1+1=3322..),)(,x x x x x x ϕϕ,故迭代收敛. 在(C)中,1<55150≈31+1361⨯2<+132='+1=32232232.).(.)(),,)(//x x x x x ϕϕ,故迭代收敛.在(D)中,类似证明,迭代收敛.4.牛顿切线法是用曲线f (x )上的 与x 轴的交点的横坐标逐步逼近f (x )=0的解;而弦截法是用曲线f (x )上的 与x 轴的交点的横坐标逐步逼近f (x )=0的解. 答案:点的切线;两点的连线解答:见它们的公式推导.三、练习题1. 用二分法求方程f (x )=0在区间[a ,b ]内的根x n ,已知误差限,确定二分的次数n 是使( )(A) b -a (B) f (x ) (C)x *-x n (D)x *-x n b -a2. 设方程f (x )=x -4+2x =0,在区间[1,2]上满足 ,所以f (x ) =0在区间[1,2]内有根.建立迭代公式x x 2-4==ϕ(x ),因为 ,此迭代公式发散.3. 牛顿切线法求解方程f (x )=0的近似根,若初始值x 0满足( ),则解的迭代数列一定收敛.(A))()(00''x f x f <0 (B) )()(00''x f x f >0(C))()(00''x f x f 0 (D))()(00''x f x f 04. 设函数f (x )区间[a ,b ]内有二阶连续导数,且f (a )f (b )<0, 当 时,则用弦截法产生的解数列收敛到方程f (x )=0的根.5. 用二分法求方程x 3-x -1=0在区间[1.0,1.5]内的实根,要求准确到小数点后第2位.6. 试用牛顿切线法导出下列各式的迭代格式:(1) c 1不使用除法运算; (2) c1不使用开方和除法运算. 四、练习题答案 1.(C) 2.0>20<1)(,)(f f ; 3861≈22>22-='.ln ln )(x x ϕ>13.(B )4. f (x )05. 1.326. (1) 31+21+50-51=2-2=nn n n n n cx x x cx x x ..)(,附录:教材中练习与习题答案练习13.1 (A )1. 0.9212.3.1455(B )1.B2.C3. ln()ln ln b a --ε24.[1.5,2] 练习13.2 (A )1.(1) 0.7391 (2) 0.20391 (3)2.09455(B )1.A2.x x =ϕ()的形式3.B4.Cϕϕ()(~)(~)~x x x x x x x x n n n n n n n+++++---+1112112 练习13.3 (A ) 1.x 02= x 1=1.5829 x 2=1.5009 x 3=x 4=1.49732. x 1=1.5970149 x 2=1.5945637 x 3=x 4=1.59456213. 3.317(B)1.线性化2.x f x f x n n n n ----'=11112()()(,,) ,0≠'1-)(n x f 3.C 4.B 5.,...),,()(210=5+231=21-1-n x x x n n n 6.D 练习13.4 (A )1. 2.09455 2. 3.14619(B ) 1. ,...),()()()()(32=---=2-1-2-1-1-1-n x x x f x f x f x x n n n n n n n 2.A 3.D 4.C练习13.5 (A)1. x ***..=⎡⎣⎢⎤⎦⎥=-⎡⎣⎢⎤⎦⎥x x 12025409836049726776 2.x x 1121125225()()..⎡⎣⎢⎤⎦⎥=⎡⎣⎢⎤⎦⎥,x x 1222120277778()().⎡⎣⎢⎤⎦⎥=⎡⎣⎢⎤⎦⎥ (B)1. 线性化2.B3.C习题131.(1) -0.56714 (2)1.267172.(1)1.9405 (2) 0.4502 (3)-1.154173.(1) 3.63197 (2) 0.011754.(1)0.33767 1.30749 (2)1.386165. (1) (0.8260,0.5636)T (2) (0.9991,-0.0888)T。

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

《计算机数学基础(下)》数值分析部分辅导(1)
中央电大 冯 泰
第9章 数值分析中的误差
一、重点内容
误差 设精确值x *的近似值x ,差e =x -x *称为近似值x 的误差(绝对误差)。

误差限 近似值x 的误差限ε是误差e 的一个上界,即ε≤-=*
x x e 。

相对误差e r 是误差e 与精确值x *
的比值,*
*
*-==x x x x e e r 。

常用x e e r
=计算。

相对误差限r ε 是相对误差的最大限度,r r e ≥ε,常用x
ε
计算相对误差限。

绝对误差的运算:
)()()(2121x x x x εεε+=± )()()(122121x x x x x x εεε+≈
22
122121
+=x x x x x x x )()()(
εεε 有效数字 如果近似值x 的误差限ε 是它某一个数位的半个单位,我们就说x 准确到该
位. 从这一位起到前面第一个非0数字为止的所有数字称为x 的有效数字.
关于有效数字:
(1) 设精确值x *的近似值x ,
m n a a a x 10.021⨯±=
a 1,a 2,…,a n 是0~9之中的自然数,且a 1≠0,
n l x x l m ≤≤110⨯50=≤--,.*ε
则x 有l 位有效数字.
(2) 设近似值m n a a a x 10.021⨯±= 有n 位有效数字,则其相对误差限
1+-1
10⨯21

n r a ε
(3) 设近似值m n a a a x 10.021⨯±= 的相对误差限不大于
1110)
1(21
+-⨯+n a
则它至少有n 位有效数字.
(4) 要求精确到10-
3,取该数的近似值应保留4位小数。

一个近似值的相对误差是与准确数字有关系的,准确数字是从一个数的第一位有效数字一直数到它的绝对误差的第一位有效数字的前一位,例如具有绝对误差e=0.0926的数x =20.7426只有三位准确数字2,0,7。

一般粗略地说,具有一位准确数字,相对于其相对误差为10%的量级;有二位准确数字,相对于其相对误差为1%的量级;有三位准确数字,相对于其相对误差为0.1%的量级。

二、实例
例1 设x *= π=3.1415926…
近似值x =3.14=0.314×101, 即m=1,它的误差是 0.0015926…,有
3-1*10⨯50≤0015260=-.. x x
即l =3,故x =3.14有3为有效数字。

x =3.14准确到小数点后第2位。

又近似值x =3.1416,它的误差是0.0000074…,有
5-1*10⨯50≤00000740=-.. x x
即m =1,l =5,x =3.1416有5位有效数字。

而近似值x =3.1415,它的误差是0.0000926…,有
4-1*10⨯50≤00009260=-.. x x
即m =1,l =4,x =3.1415有4位有效数字。

这就是说某数有s 位数,若末位数字是四舍五入得到的,那么该数有s 位有效数字;若末位数字不是四舍五入得到的,那么该数有s 位或s -1位有效数字。

例2 指出下列各数具有几位有效数字,及其绝对误差限和相对误差限: 2.000 4 -0.002 00 9 000 9 000.00
解 因为x 1=2.000 4=0.200 04×101, 它的误差限0.000 05=0.5×10 1―
5,即m=1,l=5,故
x =2.000 4有5位有效数字. 相对误差限%...50020=4
000205
0000=
r ε
x 2=-0.002 00,误差限0.000 005,因为m=-2,l =3,x 2=-0.002 00有3位有效数字。

相对误差限εr =0.000 05/0.002 00=0.25%。

x 3=9 000,绝对误差限为0.5,因为m=4, l=4, x 3=9 000有4位有效数字,相对误差限
εr =0.5/9 000=0.005 6%
x 4=9 000.00,绝对误差限0.005,因为m=4,l =6,x 4=9 000.00有6位有效数字,相对误差限为εr =0.005/9 000.00=0.000 056% 由x 3与x 4可以看到小数点之后的0,不是可有可无的,它是有实际意义的。

例3 ln2=0.69314718…,精确到10-
3的近似值是多少?
解 精确到10-
3=0.001,即绝对误差限是ε=0.05,故至少要保留小数点后三位才可以。

Ln2≈0.693。

三、练习题 1. 设某数x *,它的保留三位有效数字的近似值的绝对误差是 。

2. 设某数x *,它的精确到10-
4的近似值应取小数点后 位。

3. ( )的3位有效数字是0.236×102。

(A) 235.54×10-1 (B) 235.418 (C) 2354.82×10-
2 (D) 0.0023549×10
3 4. 设a *=2.718181828…,取a=2.718,则有( ),称a 有四位有效数字。

(A) 4
-*10⨯50≤-.a a (B)4-1*10⨯50≤-.a a (C) 4
-*10≤-a a (D) 00030≤-*.a a
5. 设某数x *,对其进行四舍五入的近似值是( ),则它有3位有效数字,绝对误差限是
4-10⨯2
1。

(A) 0.315 (B) 0.031 50 (C) 0.0315 (D) 0.00315
6. 以下近似值中,保留四位有效数字,相对误差限为3
-10⨯250.。

(A) 0.01234 (B) –12.34 (C) –2.20 (D) 0.2200
7. 将下列各数舍入成三位有效数字,并确定近似值的绝对误差和相对误差。

(1) 2.1514 (2) -392.85 (3) 0.003922
8. 已知各近似值的相对误差,试确定其绝对误差:
(1) 13267 e r =0.1% (2) 0.896 e r =10%
9. 已知各近似值及其绝对误差,试确定各数的有效位数。

(1) 0.3941 e=0.25×10-2 (2)293.481 e=0.1 (3) 0.00381 e=0.1×10-
4 10. 已知各近似值及其相对误差,试确定各数的有效位数。

(1) 1.8921 e r =0.1×10-
2 (2) 22.351 e r =0.15 (3) 48361 e r =1% 四、练习题答案
1.该数有效数字第四位的一半。

2 .四 3.(A) 4.(B) 5.(C) 6.(D)
7. (1)2.15, e = -0.14*10-2, e r =0.65*10-3;(2) -393 , e=-0.15, e r =0.38*10-
3;
(3)0.00392, e = -0.2*10-5, e r =0.51*10-
3
8.(1) e=0.13×102 (2) 0.9×10-
1 9. (1)
2 (2)
3 (3)2 10.(1) 3 (2)1 (3)2
附录:教材中练习与习题答案
练习9.1 (B) 1.B 2. A 3.量纲 4. 半个 5.D 6. 3位 7.C 8. B 9. ≤ 10. 0.5mm 练习9.2 (B) 1. )()(21x x εε+ )()(21x x εε+
2. )()(1221x x x x εε+,
2
2
1221)
()(x x x x x εε-
3. 舍入误差不增加
4. )()(x e x f '
5. 使用数值稳定的算法;防止两个相近数相减;简化计算步骤,减少运算次数;避免除数的绝对值远小于被除数的绝对值;防止大数“吃掉”小数 习题9 1. 0.00005 0.017% 四位有效数字 0.005 0.017% 四位有效数字 0.0005 0.0017% 五位有效数字 0.5 0.0017% 五位有效数字 0.00005×105 0.017% 四位有效数字 2. 5.5(Ω) 0.2375(Ω) 4.32% 3. (1) B (2)A (3)A
4.
)()()(2121+=±x x x x εεε )()()(22
12
121121±±±=
±x x x x x x x x x x r r r εεε
)()()(122121+≈x x x x x x εεε )()()(1221+≈x x x x r r r εεε
22
122121+=x x x x x x x )()()(
εεε )()()(2121
+=x x x x r r r εεε 5. 0.005
6. 0.00333…
7. 取四位. 利用定理2.。

相关文档
最新文档