数值计算方法
?随着科学技术的飞速发展,科学计算愈来 愈显示出其重要性。科学计算的应用之广 已遍及各行各业,例如:气象资料的分析 图像,飞机、汽车及轮船的外形设计,高 科技研究等都离不开科学计算。因此,作 为科学计算的数学工具数值计算方法已成 为各咼輕院校数学、物理和井算*几应用专 ?工科本科生的专业基础课,也是工科 矗究生的学位必修课。
?数值分析或数值计算方法主要是研究如何 运用计算机去获得数学问题的数值解的理 论和方法?对那些在经典数学中,用解析方法 在理论上已作出解的存在,但要求出他的解 析解又十分困难,甚至是不可能的这类数学 问题,数值解法就显得不可缺少,同时有十分
?计算机解决科学计算问题时经历的几个过 程
-实际问题——> 数学模型——> 数值计算方法——> 程序设计——> 上机运行求出解
-实际问题——〉数学模型:由实际问题应用科
学知识和数学理论建立数学模型的过程,是应
-数值计算方法——> 程序设计——> 计算结果:
根据数学模型提出求解的数值计算方法,直到 编出程序上机算出解,是计算数学的任务。
?数值计算方法重点研究:求解的数值方法
及与此有关的理论
-包括:方法的收敛性,稳定性,误差分析,计 q寸间的最小(也就是计算费用),占用内存空
数学问题的数值解法例示
数学问题的数值解法例示
数学问题的数值解法例示
数学问题的数值解法例示
?例仁1.1试求函数方程x=cosx在区间(0,彳)内的 —个根。
解 令/■(%) r - COSX,易知/'(兀)在[0,Q上是连续函数且
7171
f(0)f(-) = (-i)*-<0
本题用解析法求解较为困难?若用图解 法,可大致判定此零点位置作图像
J “兀
[y = cosx
取两曲线交点/的横坐标/为所求方程的解,从图中可以
例1.12计算定积分
rlJodx ⑵
rl
Jo
dx ⑵ I2 =
e~x dx
a)i1= —4
1 + x2
角军:(1)由牛顿一莱布尼兹公式
法有多种,如== 被积函数厶=4arctanx !{)= 4arctanl -4arctan0 =
法有多种,如== 被积函数
y 4—、
帀庄的复 Rimpson公式有 捋t 牌尿《
2 賁/(0) + 4/(|) + 2f(|) + 4f(|) + /(I)]
6424
= 3.141568627
二无原函数,因此,2 ▲
二无原函数,因此,
(2) /2 = £ edx,由 刁(兀)
2q 0.746855379。- Leibniz公式无法求解,仅可fl?值方 仍选扌韌=2, 二丄,的复jhsimpson^式进 幽
2
q 0.746855379。
例1.L3求解初值问题
dy
< dx
2xy(0) = 1
2x
该方程是5emoullij^程,令ii = y$解得
= V2x + lo本题数值方法很多,如I 择经典的四阶R-K方法:
H 一斗(rx)
々2 h
VF+F y+e)
2 y
p———-A 涉炜木。
A
现取h=0.05,其结果见下表:
XnYnyxnyny
01-00000 1.00000 1.21.84931 1.84931
0-21.18322 1.18322 1-41.94396 1.94396
0-4134164 1-34164 1.62.04939 2.04939
32051245483241.48324 1.81.61245 2.01.73205 ...
3205
1245
48324
1.48324 1.8
1.61245 2.0
1.73205 ...
2-23607 2.23607
?在任何科学计算中其解的精确性总是相对 的,而误差则是绝对的?我们从下面这个例子 就可以了解误差产生的原因.
例1.2.1试求摆长为L的单摆运动周期.
在物理学中我们知道单罢周期T = X
在物理学中我们知道单罢周期T = X
其中:/为摆长;g为自由落体加速度;加是质点
的质量。如图所示:由牛顿定律
in & 二 ma - -ml
d20
所以
ml
d20
dt2
=-mg sin 3
即等子in心
当倂艮小时,sin。a &令
CD1
咚+皿0
dt2
解微分方程得1],2 =故有
0 = cx cos劲+ C2 sin cot = Jef +c;?sin伽+ ?)
1°模型误差忽略空气阻力忽略0
1°模型误差
忽略空气阻力
忽略0点处的摩擦力
a 声
Iot 展式:sinP = & + [1
3!5!
差:g=9.8米/秒2长度
差:开方
详差的分类
?模型误差从实际问题建立的数学模型往往 都忽略了许多次要的因素,因此产生的误差 称为模型误差.
?观测误差一般数学问