计算物理课件 Chapter0.pptVIP

计算物理方法是研究并解决物理问题的数值近似解的方法，简称计算物理，也叫算法物理等。算法：从给定的已知量出发，经过有限次四则运算及规定的运算顺序，最后求出未知量的数值解，这样构成的完整计算步骤称为算法。 运算量： 一个算法所需的乘除运算总次数. 计算量是衡量一个算法好坏的重要指标。 算法优劣的重要性 例1 已知 a0, a1, a2 ,…, an, x, 计算多项式： 直接计算：运算量（乘） 秦九韶算法（1247年）：运算量：例 2 解线形方程组 其中， ? 克兰姆(Cramer)法则：运算量(乘除)： ?高斯消元法(Gauss)：运算量(乘除)Gauss: 3060次； Cramer: 本课程的任务： ● 建立各种物理问题的计算物理算法的方法和理论。 通俗地讲，就是为各种实际物理问题提供有效的数值近似解方法。 ● 提供在计算机上实际可行的、理论可靠的、计算复杂 性好的各种常用物理算法。 先修课程和后续课程：先修课程：基础物理,高等数学，线性代数，常微分方程，计算机语言等。后继课程：数值代数，数值逼近，最优化方法等。 0.2 误差与有效数字 误差的来源 1、模型误差指数学模型与实际问题之间出现的误差。 2、观测误差 由观测产生的误差。3、截断误差当数学模型不能得到精确解时，通常要用数值方法求它的近似解，其近似解与精确解之间的误差称为截断误差。如函数 f (x)用Taylor展式的有限项来近似代替。 4、舍入误差 由于受计算机字长的限制，计算时只能取有限位数进行运算，由此产生的误差称舍入误差。 例如 近似值 x1= 123.456 和 x2 = 0.0123456 。它们的绝对误差限分别为： 相对误差 先看两种产品的不合格率：8/1000=0.8% ； 12/2000=0.6% （相对误差的值） 定义：设 x* 为精确值，x 为近似值，相对误差为当绝对误差 e 较小时，相对误差可写为 相对误差限 有效数字： 设 x* 的近似值为p为整数，若 x 的绝对误差限不超过末位的半个单位，即 则称用 x 近似 x*是具有“ n 位有效数字。” 例0.1 x=12.34；y=0.004067 （均为4位）；3.00（3位），3.0000（5位）。| π－3.14|≤10－(n－p)/2=10－2 /2 ( n = 3, p = 1)。 0.3 约束误差 计算物理的选择在计算物理中，要讨论计算公式的收敛性、稳定性和截断误差。约束误差就是要选择收敛性要求低、稳定性好的计算物理。 提高计算物理精度数值在计算机中存放的位数称为字长。有限位的字长是带来舍入误差和拟制计算物理精度的根源。因此，选用字长长的计算机和增加数的有效位数可以提高计算物理精度。 0.4 范数 在实轴上，任意两点a,b之间的距离可以用绝对值| a－b |来度量。三维空间R３中向量的长度为 向量范数是用来衡量向量的大小，它是三维空间中向量长度概念的推广。 向量范数 矩阵范数 向量范数的定义 定义 任意向量x=(x1, x2,…, xn)T∈Rn，按照一定规则对应一个非负实数，记为 ‖x‖,如果满足： （1）非负性 ‖x‖≥0，而且‖x‖=0 当且仅当x=0； （2）齐次性 x∈Rn ，α∈R，有‖αx‖= |a|‖x‖； （3）三角不等式 则称‖x‖为向量 x 的范数。 Rn中常用范数形式：设x=(x1,x2,…,xn)T， 1-范数： ‖x‖1=|x1|+|x2|+…+|xn|= 2-范数： ‖x‖2= ∞-范数:‖x‖∞=max{|x1|,|x2|,…,|xn|}= p-范数： ‖x‖p= 注：定义中的规则十分重要 若定义‖x‖=x1 + x2+…+xn,就构不成范数，因为由‖x‖=0不能推出x≠0。例如，x=（1，-1，0）， 显然，‖x‖=0，但是，x≠0。 例0.2 已知 x=(1, 3,-5)，求1，2，∞-范数。‖x‖1=1+3+5=9;‖x‖2=(12+32+52)1/2;‖x‖∞ =max{1,3,|-5|}=5. 范数的等价性和等价性定理 定义 设‖?‖α 和‖?‖β 为Rn 上定义的两种范数，若存在非负常数C1, C2,使?x∈ Rn 成立C1‖?‖α≤‖?‖β≤ C2‖?‖α，则称‖?‖α 和‖?‖β 为Rn 上的等价范数。 定理 0.1 有限维空间Rn上的所有范数都等价。 容易验证：（1）‖x‖2≤‖x‖1≤ n1/2‖x‖2；（2）‖x‖∞≤‖x‖2≤ n1/2‖x‖∞；（3）‖x‖∞≤‖x‖1≤ n‖x‖∞。 证明提示：（1）由 x12+ x22+ … + xn2≤(|x1|+ |x2|+… +| xn|) 2， 两边开方