正在加载图片...
Jacobi迭代法的 MATLAB函数文件 Jacobi m如下: function y, n=jacobi(A, b, XO, eps) if nargin==3 eps=10e-6; elseifnargin<3 error return end D=diag(diag(A);%求A的对角矩阵 L=tri(A,-1);%求A的下三角阵 U=triu(A,1);%求A的上三角阵 B=DL+U) =D\b y=B X0+f; n %迭代次数 while norm(y-x0)>=eps XO y=BX0+f: n=n+1; endJacobi迭代法的MATLAB函数文件Jacobi.m如下: function [y,n]=jacobi(A,b,x0,eps) if nargin==3 eps=1.0e-6; elseif nargin<3 error return end D=diag(diag(A)); %求A的对角矩阵 L=-tril(A,-1); %求A的下三角阵 U=-triu(A,1); %求A的上三角阵 B=D\(L+U); f=D\b; y=B*x0+f; n=1; %迭代次数 while norm(y-x0)>=eps x0=y; y=B*x0+f; n=n+1; end
<<向上翻页向下翻页>>
©2008-现在 cucdc.com 高等教育资讯网 版权所有