计算方法绪论与误差等精要.ppt

计算方法研究的对象 研究数值方法的设计、分析和有关理论基础与软件实现。 计算方法的内容 连续系统的离散化 离散性方程的数值求解 数值方法和数值软件过去50年的主要进展 计算数学未来50年的展望 1.1 误差的来源 1.1.1 舍入误差(Round-off Errors) 1.2 误差 误差限 有效数字 1.3 相对误差和绝对误差 [例1.5]测得某桌面的长a的近似值a*=120cm,宽b的近似值b*=60cm。若已知|e(a*)|≤0.2cm, |e(b*)|≤0.1cm。 试求近似面积s*=a*b* 的绝对误差限与相对误差限。 1.5 在近似计算中需要注意的问题 1.尽量简化计算步骤，减少乘除运算的次数 例如，计算多项式 通常运算的乘法次数为 若采用递推算法, 则乘法次数仅为n. 又如 2.防止大数“吃掉”小数 当|a||b|时，尽量避免a+b 。例如，假设计算机 只能存放10位尾数的十进制数，则 108+0.04=108 3.尽量避免相近数相减 例如，当x很大时，应 4.避免绝对值很小的数做分母 当|b||a|时，应尽量避免a/b 病态方程组： 如果方程组AX=b由于A或b的小扰动而导致解严重失真，则此方程组称为病态方程组，否则称为良态方程组。 判定一个病态方程组的简单方法； 病态方程组一般不能用解方程组的常有方法求解，而采用“迭代求精法”来计算。（列主元消元法的应用） 追赶法 追赶法仍然保持LU分解特性，它是一种特殊的LU分解。充分利用了系数矩阵的特点，而且使之分解更简单，得到对三对角线性方程组的快速解法。 因三对角矩阵的非零元素呈“带状”，我们也因此将它叫做带状矩阵。 [解]先LU分解系数矩阵，由2.17式得 再由2.18式得 对第一个方程组，由2.15式得 再由2.16式得 对第二个方程组，由2.15式得 再由2.16式得 2．平方根法 定理 设A对称正定，则有非奇异下三角阵L，使 ---- 理论基础 (证略） 分解方法：设 ( choleskey分解) 将例2.1用增广矩阵的变化写成 高斯消去法的求解过程分为两个阶段: 首先，把原方程组化为上三角形方程组，称之为“消元”过程； 然后，用逆次序逐一求出三角方程组(原方程组的等价方程组)的解，并称之为“回代”过程。 2.2 高斯消去法(Gaussian elimination) 消元 ：将(2.1)式写成 矩阵形式 (2.3) (2.4) 第1步:若a11 (1) ≠0，用第二个方程减去第一个方程乘以a21(1)/ a11(1),用第三个方程减去第一个方程乘以a31(1)/ a11(1)… 则有 矩阵形式 (2.5) (2.6) 第2步:若a22 (2) ≠0，用第三个方程减去第二个方程乘以a32(2)/ a22(2),用第四个方程减去第二个方程乘以a42(2)/ a22(2)… 则有 矩阵形式 (2.7) (2.8) 第k步:若akk (k) ≠0，用第k+1个方程减去第k个方程乘以ak+1k(k)/ akk(k) … 则有 矩阵形式 (2.9) (2.10) 重复n－1次，得到等价的上三角形方程组 矩阵形式 (2.11) (2.12) 以上过程把系数矩阵A(1)变成上三角矩阵A(n)，称之为消元，计算公式可归纳为 (2.13) 回代 (2.12) 2.3 主元素消去法 因此，x1＝0 x2＝1 因此，x1＝1 x2＝1 精确解，x1＝10000/9999 x2＝9998/9999 在做除法运算时，选取绝对值大的作分母。 ——主元素消去法的基本思路。 列主元素消去法 列主元素消去法基本思想 1.用高斯消去法求解线性方程组时,应避免小的主元.在实际计算中,进行第k步消去前,应该在第k列元素aik (i=k,…,n)中找出绝对值最大者,例如 ∣a ∣= max ∣a ∣ 2.再把第p个方程与第k个方程组进行交换,使apk(k-1)成为主元.我们称这个过程为选主元素.由于只在第k列元素中选主元素,通常也称为按列选主元素(或称部分选主元). 3.如果在第k步消去前,在第k个方程到第n个方程所有的xk到xn的系数a (i=k,…,n;j=k,…,n)中,找出绝对值最大者,例如 ∣a ∣=max∣a ∣ 再交换第k,p两个方程和第k,q两个未知量的次序,使a 成为主元素. 称这个过程为完全选主元素。 4.不论是哪种方式选出主元素,而后再按上面介绍的计算步骤进行消去的计算,一般都称为主元素高斯消去法。在实际计算中,常用按列选主元素的