数值分析上机作业.docVIP

《数值分析》上机作业 姓名： 学号： 学院： 专业： 第一章 4题 设计程序如下： clear; %定义x，y的区间 x=-10:0.2:10; %步长为100 y=-10:0.2:10; %在XOY平面生成二维网格数据 [X,Y]=meshgrid(x,y); %对二元函数进行表达 a=-abs(X);b=X+Y;c=X.^2+Y.^2+1; Z=exp(a)+cos(b)+1./c; %绘制三维图形 mesh(X,Y,Z); 图1. 区间等分100份 图2.区间等分200份 图3.区间等分400份 第二章 2题 设计程序如下： Clear n=8; V=220; R=27; a=[0 -2 -2 -2 -2 -2 -2 -2]; b=[2 5 5 5 5 5 5 5]; c=[-2 -2 -2 -2 -2 -2 -2]; I=[V/R 0 0 0 0 0 0 0]; for i=2:n a(i)=a(i)/b(i-1); b(i)=b(i)-c(i-1)*a(i); I(i)=I(i)-a(i)*I(i-1); end; I(n)=I(n)/b(n); for i=n-1:-1:1 I(i)=(I(i)-c(i)*I(i+1))/b(i); end; I; fprintf(“各电路的电流量I为”); 运行结果如下： 各电路的电流量I为 I=8.1478 4.0737 2.0365 1.0175 0.5073 0.2506 0.1194 0.0477 第三章 2题 高斯赛德尔迭代 设计程序如下： ：A=[10,1,2,3,4;1,9,-1,2,-3;2,-1,7,3,-5;3,2,3,12,-1;4,-3,-5,-1,15;] x=[0;0;0;0;0]; b=[12;-27;14;-17;12] c=0.000001 L=-tril(A,-1) U=-triu(A,1) D=(diag(diag(A))) X=inv(D-L)*U*x+inv(D-L)*b; k=1; while norm(X-x,inf)= c x=X; X=inv(D-L)*U*x+inv(D-L)*b; k=k+1; end X k 计算结果：X = 1.0000 -2.0000 3.0000 -2.0000 1.0000 k =37 Jacobi迭代法： 设计程序如下： Jacobi：b=[12;-27;14;-17;12] x = [0;0;0;0;0;] k = 0; r = 1; e=0.000001 A=[10,1,2,3,4;1,9,-1,2,-3;2,-1,7,3,-5;3,2,3,12,-1;4,-3,-5,-1,15;] D = diag(diag(A)); B = inv(D)*(D-A); f = inv(D)*b; p = max(abs(eig(B))); if p = 1 迭代法不收敛 return end while r e x0 = x; x = B*x0 + f; k = k + 1; r = norm (x-x0,inf); end x k 计算结果：x = 1.0000 -2.0000 3.0000 -2.0000 1.0000 k =65 SOR：A=[10,1,2,3,4;1,9,-1,2,-3;2,-1,7,3,-5;3,2,3,12,-1;4,-3,-5,-1,15] x=[0;0;0;0;0]; b=[12;-27;14;-17;12] e=0.000001 w=1.44; L=-tril(A,-1) U=-triu(A,1) D=(diag(diag(A))) X=inv(D-w*L)*((1-w)*D+w*U)*x+w*inv(D-w*L)*b n=1; while norm(X-x,inf)=e x=X; X=inv(D-w*L)*((1-w)*D+w*U)*x+w*inv(D-w*L)*b; n=n+1; end X n 计算结果：X = 1.0000 -2.0000 3.0000 -2.0000 1.0000 n =22 结果分析：由迭代次数的比较可知该情况下Jacobi迭代法比Gauss-Seide