用SOR迭代法解法方程① SOR.m文件为 function [x,n]=SOR(a,b,E,p) temp1=size(a);temp2=size(b); if temp1(1)~=temp1(2)|temp1(1)~=temp2(1)|temp2(2)~=1|det(a)==0 error('矩阵或向量的大小不对应或矩阵为奇异矩阵') return; end %%化为三角矩阵过程 leng=temp1(1);N=-tril(a,-1)-triu(a,1);M=a+N;N=-a; t=1;
高斯迭代法解矩阵方程 function [a,b] = gaussj(a,n,np,b,m,mp) for j=1:n ipiv(j)=0; end for i=1:n big=0; for j=1:n if(ipiv(j)~=1) for k=1:n if (ipiv(k)==0) if (abs(a(j,k))>=big) big=abs(a(j,k)); irow=j; icol=k; end elseif (ipiv(k)>1) pause, dips('gaussj')
