工程硕士 数值分析 第一章 绪论.pptVIP

§1.1 计算方法概论 应用问题举例 2、天体力学中的Kepler方程 全球定位系统：在地球的任何一个位置，至少可以同时收到4颗以上卫星发射的信号 3、全球定位系统（Global Positioning System, GPS) 非线性方程组的数值方法 4、已经测得在某处海洋不同深度处的水温如下： 深度（M） 466 741 950 1422 1634 水温（oC）7.04 4.28 3.40 2.54 2.13 根据这些数据，希望合理地估计出其它深度（如500米，600米，1000米…）处的水温 5、人口预测 下面给出的是中国1900年到2000年的人口数，我们的目标是预测未来的人口数（数据量较大时） 本课程第四章的内容：曲线拟合 6、铝制波纹瓦的长度问题 这个问题就是要求由函数f(x)=sin x给定的曲线从x=0到x=48英寸间的弧长L. 由微积分学我们知道,所求的弧长可表示为: 7、生物化学反应的例子 A，B，C是三种蛋白质，其反应如下： 我们通过建模可以得到如下方程组 本课程第六章的内容： 常微分方程的数值方法 数值分析的特点 1、方法是近似的； 2、与计算机不能分离：上机实习 （掌握一门语言：C语言或Fortran语言， 会用一种数学软件：Matlab或Mathematica，Maple） 数值算法的要求: 1. 仿真性：模型尽可能仿真实际问题; §1.3 误差基础知识 舍入误差很小，本课程将研究它在 运算过程中是否能被有效控制。 二．误差基本概念 1．绝对误差。设x为准确值， 为近似值。则称 　 e = x - x*为x*的绝对误差。如果 则称 为x*的绝对误差限。 例1.3.3 测得某桌面的长a的近似值a*=120cm,宽b的近似值b* = 60cm。若已知|e(a*)|≤0.2cm, |e(b*)|≤0.1cm。 试求近似面积s*= a*b*的绝对误差限与相对误差限。 §1.4 选用算法应遵循的原则(补充内容) 1.尽量简化计算步骤，减少乘除运算的次数。 当|a||b|时，尽量避免a+b 。例如，假设计算机只能存放10位尾数的十进制数，则 本章基本要求： 1.熟悉计算方法在解决实际问题中所处的地位； 2.熟悉绝对误差（限），相对误差（限）及有效数字概念； 3.熟悉定理1.3.1； 4.了解公式(1.3)--(1.9)； 5.熟悉选用算法应遵循的原则。 注：? Frobenius 范数不是算子范数。 ? 我们只关心有相容性的范数，算子范数总是相容的。 ? 即使 A中元素全为实数，其特征值和相应特征向量 仍可能是复数。将上述定义中绝对值换成复数模均成立。 若不然，则必存在某个向量范数 || · ||v 使得 对任意A 成立。 反例 ? ? 谱半径 (定义1.2.5) 矩阵A的谱半径记为 ? (A) = ，其中?i 为 A 的特征值。 Re Im ? ? ? ? ? ? ? ? ? (A) 证明： 由算子范数的相容性，得到 将任意一个特征值 ? 所对应的特征向量 代入 证明： A对称 若? 是 A 的一个特征值，则?2 必是 A2 的特征值。 又：对称矩阵的特征值为实数，即 ?2(A) 为非负实数， 故得证。 对某个 A 的特征值? 成立 所以2-范数亦称为谱范数。 对任意算子范数 || · || 有 定理1.2.2 若A对称，则有 定理1.2.3 定理1.2.4 对任意 ? 0, 存在算子范数 || · || 使得 || A || ? ? (A) + ? 。 证明：省略 设： A A k k =  ? ? lim 是指 ij k ij k a a =  ? ? ) ( lim 对所有 1? i, j ? n 成立。 等价于对 任何算子范数有 定义1.2.5 证明： “?” 若 ? 是 A 的特征值, 则?k 是 Ak 的特征值。 则 [? (A)]k = [ max | ? | ]k = | ?mk | ? ? ( Ak ) ? ? (A) 1 “?” 由定理1.2.4和 ? (A) 1 可知 存在算子范数|| · || 使得 || A || 1。 || Ak || ? || A ||k ? 0 （ k ? ?） Ak ? 0 ? ? ( A ) 1 定理1.2.5 ? || Ak || ? 0 （由定理1.2.2得） 证明： ① 若不然，则 有非零解，即存在非零向量 使得 ? ② 若矩阵