雅可比迭代法.docVIP

雅可比迭代法

班级：计算062姓名：王保翔3060811028

目的意义：设方程组Ax=b的系数矩阵A非奇异而且，将A分裂为A=D+L+U,可以使计算简便。其中

，

思想：A=D+L+U，其中

，

将方程组乘以，得到等价的方程组，i=1，2，…n

算法：

Step1：给定一组x，即初值。

Step2：用for循环计算：

x[k+1]=(b[i]-)/a[k][k].

Step3:当abs(x[k+1]-x[k])eps时停止。

程序：

#includestdio.h

#definep0.1

#defineN5

voidmain()

{

inti,j,k,s,n;

doublex[N],a[N][N],b[N],q[N],c[N],d[N];

doublet,u;

printf(Inputn:);

scanf(%d,n);

for(i=0;i=n;i++)

{

printf(Inputx[%d],b[%d]:,i,i);

scanf(%lf,%lf,x[i],b[i]);

}

for(i=0;i=n;i++)

for(j=0;j=n;j++)

{

printf(Inputa[%d][%d]:,i,j);

scanf(%lf,a[i][j]);

}

for(i=0;i=n;i++)

for(j=0;j=n;j++)

printf(a=%lf\n,a[i][j]);

for(k=0;;k++)

{

t=0;

for(s=0;s=n;s++)

d[s]=x[s];

for(i=0;i=n;i++)

{

u=0;

d[i]=x[i];

for(j=0;j=n;j++)

{

if(j!=i)

u=u+a[i][j]*d[j];

elsecontinue;

}

printf(u=%lf\n,u);

x[i]=(b[i]-u)/a[i][i];

q[i]=x[i]-d[i];

printf(x[%d]=%lf\n,i,x[i]);

}

for(i=0;i=n;i++)

{

if(q[i]=0)

c[i]=q[i];

elsec[i]=-q[i];

t=t+c[i];

}

if(tp)

break;

}

printf(k=%d\n,k);

for(i=0;i=n;i++)

printf(x[%d]=%lf\n,i,x[i]);

}算例及运算结果：

分析：输入初始向量和矩阵为数经过k次迭代后输出x的值，结果正确可靠。

参考文献

[1]刑志栋，矩阵数值分析，陕西：陕西科学技术出版社，2005。

[2]谭浩强，C语言程序设计，北京：清华大学出版社，2005。

[3]翁惠玉，c语言程序设计思想与方法，北京：人民邮电出版社，2008