计算方法第三章新.pdfVIP

第3章 线性代数方程组的数值解法 《九章算术》方程术： 今有上禾三秉，中禾二秉，下禾一秉，实三十九斗；上禾二 秉，中禾三秉，下禾一秉，实三十四斗；上禾一秉，中禾二 秉，下禾三秉，实二十六斗。问上、中、下禾实一秉各几何？ 答曰： 上禾一秉，九斗、四分斗之一， 中禾一秉，四斗、四分斗之一， 下禾一秉，二斗、四分斗之三。 3x+2y+z=39 2x+3y+z=34 X+2y+3z=26 问题的提出: … a11 x1 + a12 x2 + + a1n xn = b1 a x + a x + … + a x = b 21 1 22 2 2n n 2 Ax=b …… … an1 x1 + an2 x2 + + ann xn = bn 克莱姆法则 ( Cramer ): D的第i 列换成b D D= det(A), X = i (i = 1,2,…….,n ) i D 解n阶方程组要算（n+1）个n阶行列式， 计算一个n阶行列式需要 n！=n*(n-1)*…*2*1 次乘法。 例如: n=20，工作量:（20+1）*20！=21！=9.7*1020次乘法， 假定1秒钟做3千亿次乘法，需要： 9.7× 1020 ≈100 年。 310×××××11 606024365 1949年夏末，哈佛大学教授列昂惕 夫把美国经济分解为500个部门，例如煤炭 工业、汽车工业、交通系统等。对每个部 门，他写出了一个描述该部门的产出如何 分配给其他经济部门的线性方程。 由于当时最大的计算机Mark Ⅱ还不 能处理包含500个未知数的500个方程的方 程组，列昂惕夫只好把问题简化为包含42 个未知数的42个方程的方程组，计算机运 算了56个小时，才得到最后的答案。 列昂惕夫获得了1973年诺贝尔经济 学奖，他打开了研究经济数学模型的新时 代的大门。 §3.1 消去法 一、三角形方程组的求解方法 设给定非奇异线性方程组： ?ax++???+= ax ax b ? 11 1 12 2 1nn 1 ? ax22 2+ ???+ ax 2nn = b 2 ? ? axnn n= b n xnnnn= ba/ n xbii=?()/1,,1∑ axa ikkii in =???? ki=+1 递推公式称回代法。 回代法 function X=backsub(A,B) %Input -A is an nxn upper-triangular nonsingular matrix ％ - B is an n x 1 matrix %Output - X is the solution to the linear system AX=B %Find the dimension of B and initialize X n=length(B); X=zeros(n,1); X(n)=B(n)/A(n,n); for k=n-1:-1:1 X(k)=(______________________)/A(k,k); end 二．Gauss消去法 由两个过程组成: a）消元－－通过变换，把原方程变为等价的三角形方程组。 b）回代－－解三角形方程组。 例 解方程组： ?2431xxx123?+= ? ?421127xx12+ += x 3 ? ??+64xxx123 ?= 13 ?2 - 4 3 1 ? 2 ? - 4 3 1 ? 2 ? - 4 3