北京科技大学应用计算方法课后复习与标准答案.docx

、第一次作业 (一) 2-6 计算下列向量的 1-范数、∞ -范数、 2-范数 (1) x=(12,-4,-6,2) T A=[12,-4,-6,2] A = 12 -4 -6 2 norm(A,1) ans = 24 norm(A,inf) ans = 12 norm(A,2) ans = 14.1421 (2) x=(1,3,-4) T A=[1,3,-4] A = 1 3 -4 norm(A,1) ans = 8 norm(A,inf) ans = 4 norm(A,2) ans = F 范数 F 范数 (二) 2-9 计算下列矩阵的行范数、列范数、谱范数、 3 1 1 (1) A 1 1 1 2 1 1 A=[3,-1,1;1,1,1;2,1,-1] A = TOC \o 1-5 \h \z 3 -1 1 1 1 1 2 1 -1 norm(A,1) ans = 6 norm(A,inf) ans = 5 norm(A,2) ans = 3.7888 norm(A,fro) ans = 4.4721 0a (2) A ,a R a0 A=[0,1;-1,0] A = 0 1 -1 0 norm(A,1) ans = 1 norm(A,inf) ans = 1 norm(A,2) ans = 1 norm(A,fro) ans = 1.4142 第二次作业 用牛顿迭代法求方程 x3 3x 1 0在 x0 2附近的根。要求：给成程序和运行结 果. 1、牛顿法的基本原理 在求解非线性方程 在求解非线性方程 f (x) 0时，它的困难在于 f (x) 是非线性函数，为克服这一困难， 考虑它的线性展开。设当前点为 xk ，在 xk 处的 Taylor 展开式为 f(x) f(xk) f (xk )(x xk) (2.1) 令 f (x) 0 ，解其方程得到 xk 1 xk f(xk )f (xk ) ,(k 0,1, ) (2.2) (2.2)式为牛顿迭代公式， 用牛顿迭代公式求方程 f (x) 0根的方法称为牛顿迭代法。 此即 牛顿迭代法的设计原理。 2、 matlab 程序代码 function root=NewtonRoot(f,a,b,eps) if (nargin==3) eps=1.0e-4; end f1=subs(sym(f),findsym(sym(f)),a); f2=subs(sym(f),findsym(sym(f)),b); if (f1==0) root=a; end if (f2==0) roor=b; end if (f1*f20) disp( 两端点函数值乘积大于 0！ ) return ; else tol=1; fun=diff(sym(f)); fa=subs(sym(f),findsym(sym(f)),a); fb=subs(sym(f),findsym(sym(f)),b); dfa=subs(sym(fun),findsym(sym(fun)),a); dfb=subs(sym(fun),findsym(sym(fun)),b); if (dfadfb) root=a-fa/dfa; else root=b-fb/dfb; end while (toleps) r1=root; fx=subs(sym(f),findsym(sym(f)),r1); dfx=subs(sym(fun),findsym(sym(fun)),r1); root=r1-fx/dfx; tol=abs(root-r1); end end 3、运行结果截图 结论：通过计算可以看出 x3 3x 1 0在 x0 2附近的根为 1.8794. 三、第三次作业 编写高斯顺序消元法求解下面方程组的程序并计算结果 10x1x 10x1 x2 x3 6.2 x1 10x2 2x3 8.5 2x1 x2 5x3 3.2 1、高斯顺序消元法的设计原理 高斯顺序消元法的基本思想是将线性方程组 a11x1 a12 x2 a1n xn b1 a21x1 a22x2 a2n xn b2 an1x1 an2x2 ann xn bn (3.1) (3.2) 通过消元，逐步转化为等价的上(或下)三角形方程组，然后用回代法求解。 2、 matlab 程序代码 function [x,XA]=GaussXQByOrder(A,b) %高斯顺序消元法 N = size(A); n = N(1); for i=1:(n-1) for j=(i+1):n if (A(i,i)==0) disp( 对角元素为 0！ ); %防止对角元素为 0 return ; end l = A(j,i); %消元方