数值分析考试复习题(附matlab程序实现).docVIP

模式识别_周晓勇_Z201102021 PAGE PAGE 11 期末作业 给定一个方程组试建立一个收敛的迭代格式，并说明收敛的理由 解： A=[2 1 4;1 4 1;3 1 1]; b=[6 3 2]; %对A|b进行行变换,使A为严格对角占优阵 %如果A为严格对角占优阵,则解Ax=b的Jacobi迭代法和Gauss-Seidel迭代法均收敛 A=[3 1 1;1 4 1;2 1 4]; b=[2 3 6]; A\b%正确答案 w=[1:20]*0.1; for i=1:20 [x,k,flag]=SOR(A,b,1e-5,w(i),100);%函数见SOR.m result(:,i)=[w(i);flag;k;x]; end [C I]=min(result(3,:)); plot(w,result(2,:).*result(3,:),w(I),C,r*);xlabel(w);ylabel(k); legend(k-w曲线,[k=,num2str(C), ,\omega=,num2str(w(I))]); title([方程组解为：,num2str(result(4:end,I),%10.4f)]); 结果： 证明方程在区间内有根，并建立收敛的迭代格式，使用任意初值都收敛，并说明理由 解： [x,err,k,y]=NewtonDownhill(x^3+