计算机专业---数值计算方法.docxVIP

数值计算方法 ?随着科学技术的飞速发展，科学计算愈来 愈显示出其重要性。科学计算的应用之广 已遍及各行各业，例如：气象资料的分析 图像，飞机、汽车及轮船的外形设计，高 科技研究等都离不开科学计算。因此，作 为科学计算的数学工具数值计算方法已成 为各咼輕院校数学、物理和井算*几应用专 ?工科本科生的专业基础课，也是工科 矗究生的学位必修课。 ?数值分析或数值计算方法主要是研究如何 运用计算机去获得数学问题的数值解的理 论和方法?对那些在经典数学中，用解析方法 在理论上已作出解的存在，但要求出他的解 析解又十分困难，甚至是不可能的这类数学 问题，数值解法就显得不可缺少，同时有十分 ?计算机解决科学计算问题时经历的几个过 程 -实际问题—— 数学模型—— 数值计算方法—— 程序设计—— 上机运行求出解 -实际问题——〉数学模型：由实际问题应用科 学知识和数学理论建立数学模型的过程，是应 -数值计算方法—— 程序设计—— 计算结果: 根据数学模型提出求解的数值计算方法，直到 编出程序上机算出解，是计算数学的任务。 ?数值计算方法重点研究：求解的数值方法 及与此有关的理论 -包括：方法的收敛性，稳定性，误差分析，计 q寸间的最小（也就是计算费用），占用内存空 数学问题的数值解法例示 数学问题的数值解法例示 数学问题的数值解法例示 数学问题的数值解法例示 ?例仁1.1试求函数方程x=cosx在区间(0,彳)内的 —个根。 解 令/■(%) r - COSX,易知/(兀)在［0,Q上是连续函数且 71 71 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)] 6 4 2 4 = 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 hf+ay + ) ^n^KF+営++) VF+F y+e) 2 y p———-A 涉炜木。 A 现取h=0.05,其结果见下表: Xn Yn y xn yn y 0 1-00000 1.00000 1.2 1.84931 1.84931 0-2 1.18322 1.18322 1-4 1.94396 1.94396 0-4 134164 1-34164 1.6 2.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长度 差:开方 详差的分类 ?模型误差从实际问题建立的数学模型往往 都忽略了许多次要的因素，因此产生的误差 称为模型误差. ?观测误差一般数学问题包含若干参数，他们 是通过