研究生作业封面.docVIP

2012级研究生 《计算方法》作业 姓名： 高强 学号： 专业： 材料科学与工程 学院： 新金属材料国家重点实验室 成绩：__________________ 任课教师：数理学院卫宏儒 2012年11月9日 作业一 算法：||||1= ||||∞= ||||2= 程序：x=[12;-4;-6;2] norm(x,1) norm(x,2) norm(x,inf) 结果截图： 程序：x=[1;3;-4] norm(x,1) norm(x,2) norm(x,inf) 结果截图： 算法：||A||1= ||A||∞= ||A||2=,是ATA最大特征值 程序：A=[3,-1,-1;1,1,1;2,1,1]; norm(A,1) norm(A,2) norm(A,inf) norm(A, fro) 结果截图： 程序：A=[0,1;-1,0]; norm(A,1) norm(A,2) norm(A,inf) norm(A, fro) 结果截图： 作业二 算法：设pn是方程f(x)=0的一个近似根，将f(x)在pn处作泰勒展开，=+( -)+()2 取展开式的线性部分近似(称为的线性化)，得近似的线性方程≈+( -)=0 设≠0，令上式的解为，得=-，n=0,1, … 程序：p0=2;f=@(x)x^3-3*x-1;tol=1e-3; g=@(x)3*x^2-3;N=2000; for k=1:N p1=p0-f(p0)/g(p0) if abs(p1-p0)tol break end p0=p1; end disp(p1);disp(k) 结果截图： 作业三 算法：对于线性方程组x=b,如果a11≠0，则保留第一个方程，利用a11≠0分别消去其下面的n-1方程中的变量x1,即分别用乘数mi1=-ai1/a11 (i=2,3,…,n)去乘第一个方程，然后依次加到第个方程上就完成了第一步消元。第k步消元时，假设≠0,保留前面第1到第k个方程，利用≠0,消去下面的n-k个方程中的变量xk。依此类推，经过n-1次消元后，方程组就化为同届的上三角方程组x=b(n-1)。利用回代法求解上三角方程组，即得到方程组的解。解为 xn=/ xk=,k=n-1,n-2, …,1 程序：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(高斯消元法失败) return end for i=k+1:n m=B(i,k)/B(k,k) B(i,k:n+1)=B(i,k:n+1)-m*B(k,k: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 结果截图： 作业四 算法: 设A=L+D+U,其中D=diag[a11,a22, …ann], 0 0 a12 a13 … a1n a21 0 0 a23 … a2n L＝ a31 a22 0 U= 0 ┇ ┇ ┇ ┇ an-11 an1 an2 an3 ann-1 0 0 x(k+1)=D-1(b－Lx(k+1)－Ux(k)),k=0,1,2… 整理后可得：x(k+1)=－(L+D)-1Ux(k)+(L+D)-1b 上式可简化为x(k+1)=BGS x(k)＋dGS 程序：function [y,n]=gauseidel(A,b,x0,eps) if nargin==3 eps=1.0e-3; elseif nargin3