北京科技大学应用计算方法作业与答案.doc

一、第一次作业 （一）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 = 5.0990 （二）2-9 计算下列矩阵的行范数、列范数、谱范数、范数。 （1） A=[3,-1,1;1,1,1;2,1,-1] A = 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 （2） 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 二、第二次作业 用牛顿迭代法求方程在附近的根。要求：给成程序和运行结果. 牛顿法的基本原理 在求解非线性方程时，它的困难在于是非线性函数，为克服这一困难，考虑它的线性展开。设当前点为，在处的展开式为 令，解其方程得到 式为牛顿迭代公式，用牛顿迭代公式求方程根的方法称为牛顿迭代法。此即牛顿迭代法的设计原理。 2、程序代码 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、运行结果截图 结论：通过计算可以看出在附近的根为1.8794. 三、第三次作业 编写高斯顺序消元法求解下面方程组的程序并计算结果。 高斯顺序消元法的设计原理 高斯顺序消元法的基本思想是将线性方程组 通过消元，逐步转化为等价的上（或下）三角形方程组，然后用回代法求解。 程序代码 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