用C语言实现高斯赛德尔迭代方法Word格式.docx
- 文档编号:17244783
- 上传时间:2022-11-29
- 格式:DOCX
- 页数:6
- 大小:14.98KB
用C语言实现高斯赛德尔迭代方法Word格式.docx
《用C语言实现高斯赛德尔迭代方法Word格式.docx》由会员分享,可在线阅读,更多相关《用C语言实现高斯赛德尔迭代方法Word格式.docx(6页珍藏版)》请在冰豆网上搜索。
printf("
PleaseinputthematrixAbyrow!
\n"
);
label:
for(i=0;
i<
n;
i++)
Row%d:
"
i+1);
for(j=0;
j<
j++)
scanf("
%f"
&
a[i][j]);
}
Pleaseinputthevectorb:
b[i]);
if(i==j)
c[i][j]=0;
continue;
c[i][j]=-a[i][j]/a[i][i];
\nThematrixAandvectorb:
%10.5f"
a[i][j]);
b[i]);
\nTheGauss-Seideliterativescheme:
c[i][j]);
b[i]/a[i][i]);
floatinit_vec(intn,floatx[N])
inti;
Pleaseinputtheinitialiterationvectorx:
x[i]);
\nTheinitialiterationvectorx:
x[i]);
floatgs(intn,floata[N][N],floatb[N],floatx[N])
inti,j,k;
floattmp1,tmp2,x2[N];
for(k=0;
k<
10001;
k++)
x2[i]=x[i];
tmp1=0.0;
tmp2=0.0;
i;
tmp1+=a[i][j]*x[j];
for(j=i+1;
tmp2+=a[i][j]*x2[j];
x[i]=(b[i]-tmp1-tmp2)/a[i][i];
for(i=0,j=0;
if(fabs(x2[i]-x[i])<
0.00001)j++;
if(j==n)
\nThisGauss-Seideliterativeschemeisconvergent!
\nNumberofiterations:
%d"
k+1);
break;
if(k==499)
\nThisJacobiiterativeschememaybenotconvergent!
\nTheresults:
%12.7f"
}
intmain()
intn;
floatx[N],a[N][N],b[N];
Inputn:
%d"
n);
Table(n,a,b);
init_vec(n,x);
gs(n,a,b,x);
getch();
另一种
/*seideliterativemathod*/#include<
stdio.h>
#include<
math.h>
voidload();
#definem10/*定义矩阵*/
floata[m][m];
floatc[m];
main()
floatx[m],x1[m],eps[m];
floats=0;
floatt=0;
intp=1;
intq=1;
intk=0;
floatepsi;
load();
=m-1;
s=s+fabs(a[i][j]);
t=t+fabs(a[j][i]);
q=q&
&
(s<
2*fabs(a[i][i]));
p=p&
(t<
s=0;
t=0;
if((p+q)==0)/*判断是否收敛*/
seideliterativeisnotretrain"
else
{x[i]=0;
x1[i]=0;
do
epsi=x[0]-x1[0];
j++)s=s+a[i][j]*x[j];
x[i]=(c[i]+a[i][i]*x[i]-s)/a[i][i];
eps[i]=fabs(x[i]-x1[i]);
x1[i]=x[i];
epsi=(epsi>
eps[i])?
epsi:
eps[i];
x%d=%f"
i,x[i]);
k=k+1;
while(epsi>
=1e-3);
iteratived%dtimes"
k);
voidload()
floatb[m*m];
\npleaseinputmatrixa\n"
b[j+i*m]);
a[i][j]=b[j+i*m];
pleaseinputconstvectorc\n"
for(i=0;
i++)scanf("
c[i]);
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 语言 实现 高斯赛德尔迭代 方法