《科学与工程计算方法及其应用》课程上机实验报告.docVIP

《科学与工程计算方法及其应用》课程 上机实验报告 姓名： 班级：自动化二班 学号：20110703217 日期：2013-5-15 指导老师： 刘小燕 本次实验题号：实验二 一 程序流程图 程序代码及注释 迭代法： function [K,X]=Jacobi(A,B,X0,delta,Nmax) N=length(B); for K=1:Nmax for i=1:N X(i)=(B(i)-A(i,[1:i-1,i+1:N])*X0([1:i-1,i+1:N]))/A(i,i); end err=abs(norm(X-X0)); relerr=err/(norm(X)+eps); X0=X; if(errdelta)|(relerrdelta) break end end X=X; 高斯Seidel法： function[K,X]=GSeidel(A,B,X0,delta,Nmax) N=length(B); for K=1:Nmax for i=1:N if i==1 X(1)=(B(1)-A(1,2:N)*X0(2:N))/A(1,1); elseif i==N X(N)=(B(N)-A(N,1:N-1)*(X(1:N-1)))/A(N,N); else X(i)=(B(i)-A(i,1:i-1)*X(1:i-1)-A(i,i+1:N)*X0(i+1:N))/A(i,i); end end err=abs(norm(X-X0)); relerr=err/(norm(X)+eps); X0=X; if(errdelta)|(relerrdelta) break end end X=X; 运算结果 （1）矩阵的运算 1）矩阵乘法 A=[1 2;3 4] A = 1 2 3 4 B=[5 5;7,8] B = 5 5 7 8 A^2*B ans = 105 115 229 251 2）矩阵除法 A=[1 2 3;4 5 6;7 8 9] A = 1 2 3 4 5 6 7 8 9 B=[1 0 0;0 2 0;0 0 3] B = 1 0 0 0 2 0 0 0 3 A\B ans = 1.0e+016 * -0.4504 1.8014 -1.3511 0.9007 -3.6029 2.7022 -0.4504 1.8014 -1.3511 A/B ans = 1.0000 1.0000 1.0000 4.0000 2.5000 2.0000 7.0000 4.0000 3.0000 3）使用冒号选出指定元素 A=[1 2 3;4 5 6;7 8 9] A = 1 2 3 4 5 6 7 8 9 A(1:2,3) ans = 3 6 A(2:3,:) ans = 4 5 6 7 8 9 迭代法 [K,X]=Jacobi([4 2 1 ;2 4 1;1 2 4],[39;34;26],[0;0;0],0.00001,1000) K = 43 X = 6.9048 4.4048 2.5714 高斯Seidel [K,X]=GSeidel([4 2 1 ;2 4 1;1 2 4],[39;34;26],[0;0;0],0.00001,1000) K = 9 X = 6.9048 4.4048 2.5714 四．分析及总结 此次上机实验内容是用Jacobi迭代法进行迭代运算，但是Jacobi迭代法相较于Gua