计算方法作业_最终版.docxVIP

计算方法作业 姓名： 蒋姚亮 学号： 专业： 控制工程 学院： 自动化学院 成绩：__________________ 任课教师： 卫宏儒 2012年11月 作业一： 1、计算下列向量的1-范数、-范数、2-范数。 (1) (2) 解：(1) 程序：x=[12,-4,-6,2]; norm(x,1) norm(x,inf) norm(x,2) 运行结果： 得到的1-范数、-范数、2-范数分别为：24、12、14.1421。 (2) 程序：x=[1,3,4]; norm(x,1) norm(x,inf) norm(x,2) 运行结果： 得到的1-范数、-范数、2-范数分别为：8、4、5.0990。 2、计算下列矩阵的行范数、列范数、谱范数、F范数。 (1) (2) , 解：(1) 程序：clc x=[3 -1 1;1 1 1;2 1 -1]; norm(x,inf) norm(x,1) norm(x,2) norm(x,fro) 运行结果： (2) , 取a=1; 程序：clc x=[0 1;-1 0]; norm(x,inf) norm(x,1) norm(x,2) norm(x,fro) 得到：矩阵的行范数、列范数、谱范数、F范数分别为：5、6、3.7888、4.4721； ，当a=1时，其值分别为：1、1、1、1.4142。 作业二： 1、用牛顿迭代法求方程在附近的根。要求:给出程序和运行结果，计算结果保留4位有效数字。 解：设定初值p0，最大迭代步数N，误差限Tol，近似根p1，迭代步数k，则 程序：syms x; f=x^3-3*x-1; %定义f以及df函数 df=3*x^2-3; p0=2; N=1000; Tol=1e-5; for k=1:N p1=p0-subs(f/df,x,p0); % subs为置换函数，把f/df表达式中的x用p0替代 if abs(p1-p0)Tol break end p0=p1; end disp(p1);disp(k); %迭代完毕后输出p1和k的值 运行结果： 有运行结果可知：当迭代4次时，方程的解为1.879，已经满足精度要求。 作业三： 1、编写高斯顺序消元法求解下面方程组的程序及并计算结果。 解： 原理： 高斯消去法的求解过程：分为两个阶段:首先，把原方程组化为上三角形方程组，称之为“消去”过程；然后，用逆次序逐一求出三角方程组(原方程组的等价方程组)的解，并称之为“回代”过程。 程序： % A为系数矩阵，B为增广矩阵，n为方程中未知数的个数 A=[10 -1 -1;-1 10 -2;-2 -1 5]; b=[6.2 8.5 3.2]; n=3; B=[A,b]; for k=1:n-1 if B(k,k)==0 disp(Gauss fails to solve this); return end for i=k+1:n m=B(i,k)/B(k,k); B(i,k+1:n+1)=B(i,k+1:n+1)-m*B(k,k+1:n+1); % 顺序消元过程 end end X(n)=B(n,n+1)/B(n,n); % 回代过程 for i=n-1:-1:1 X(i)=(B(i,n+1)-B(i,i+1:n)*X(i+1:n))/B(i,i); end disp(X) 运行结果： 结果显示：经高斯顺序消元法得到线性方程组的解为：，，。 作业四： 1、编写Jacobi迭代法和Seidel迭代法求解上述方程组的程序，并计算出结果。精度要求： 解：(1) Jacobi迭代程序： A=[10 -1 -1;-1 10 -2;-2 -1 5]; b=[6.2 8.5 3.2]; n=3; Tol=1e