急问一个问题

odbo1984

我编了一个高斯塞得尔的程序:
function x=gs(A,b)
%高斯塞德尔迭代法
%A*x=b
numer=length(A);
D=zeros(numer);
L=zeros(numer);
U=zeros(numer);
x0=zeros(numer,1);
for i=1:numer
D(i,i)=A(i,i);%对角矩阵
x0(i,1)=1;%初始迭代向量
end
for i=2:3
for j=1:i-1
U(j,i)=-A(j,i);%U上半部矩阵
end
for j=1:i-1
L(i,j)=-A(i,j); %L下半部矩阵
end
end
errs=1;
while errs>0.001
x1=inv(D-L)*(U*x0+b);
errs=norm(x1-x0);
x0=x1;
end
x=x0;
但是结果却是:>> a=[2 -1 1;-1 -2 3;1 3 1];
>> b=[4;5;6];
>> x=gs(a,b)

x =

NaN
Inf
-Inf
为什么呢?急!

AaronSpark

没太看懂你的程序，下面是C语言写的，可以参考一下

#include<iostream.h>
#include<math.h>
const int m=10;

void main()
{
int choice=1;

while(choice==1)
{

double a[m][m],b[m],e,x[m],y[m],w,se,max;
int n,i,j,N,k;
cout<<"\n请输入方程的个数：";
cin>>n;cout<<endl;
for(i=1;i<=n;i++)
{
cout<<"请输入第"<<i<<"个方程的各项系数：";
for(j=1;j<=n;j++)
cin>>a[j];
}
cout<<"\n请输入各个方程等号右边的常数项。\n";
for(i=1;i<=n;i++)
{
cin>>b;
}
cout<<"\n请输入最大迭代次数：";
cin>>N;
cout<<"\n请输入最大偏差：";
cin>>e;
for(i=1;i<=n;i++)
{
x=0;
y=x;
}
k=0;
while(k!=N)
{
k++;
for(i=1;i<=n;i++)
{
w=0;
for(j=1;j<=n;j++)
{
if(j!=i)
w=w+a[j]*y[j];
}
y=(b-w)/double(a);
}

max=fabs(x[1]-y[1]);
for(i=1;i<=n;i++)
{
se=fabs(x-y);
if(se>max)
max=se;
}
if(max<e)
{
cout<<endl;
for(i=1;i<=n;i++)
cout<<"x"<<i<<"="<<y<<endl;
break;
}
for(i=1;i<=n;i++)
{
x=y;
}

}
if(k==N)
cout<<"迭代失败！！"<<endl;
cout<<endl;
cout<<"--------------------------------------------------------------"<<endl;
cout<<endl;
cout<<" 1、继续运算。"<<endl;

cout<<" 2、结束运算，退出本程序。"<<endl;
cout<<endl;
cout<<"请输入你的选择：";
cin>>choice;

}
}

suffer

形成LU矩阵的时候错了

不想飞的鸟

不是太明白