这次的代码解决了先前代码中程序无法自行判断误差并跳出迭代的问题,代码实用性更强。
雅可比迭代源代码
#include<stdio.h>
int main()
{
int n,i,j,k,b;
double num1;
double num2;
printf("输入方程组的维数:\n");
scanf("%d",&n);//方程组维数
double arrD[n][n]={0};//系数矩阵D
double arrB[n][n]={0};//矩阵D-1*(l+U)
printf("输入方程组各个系数:\n");
for(i=0;i<n;i++)
{ for(j=0;j<n;j++)
{ scanf("%lf",&num1);
arrD[i][j]=num1;
}
}
for(i=0;i<n;i++)
{ for(j=0;j<n;j++)
{ if(i==j)
{
arrB[i][j]=0;
}
else{
arrB[i][j]=-arrD[i][j]/arrD[i][i];
}
}
}
printf("输入矩阵b:\n");
double arrf[n]={0};//矩阵f
for(b=0;b<n;b++)
{
scanf("%lf",&num2);
arrf[b]=num2/arrD[b][b];
}
double sum=0;
double arrx[n]={0};//初始X向量
double arrxd[n]={0};
double wc;//误差限
int m;
int g;
double sum0=0;
double sum1=0;
double wc_=0;
int maxdd=0;
int count=0;
printf("输入精度:\n");
scanf("%lf",&wc);
printf("输入最大迭代次数:\n");
scanf("%d",&maxdd);
while(1)
{
count++;
for(k=0;k<n;k++)//迭代过程
{
for(j=0;j<n;j++)
{
sum=sum+arrB[k][j]*arrx[j];
}
arrxd[k]=sum+arrf[k];
sum=0;
}//迭代过程
for(g=0;g<n;g++)
{
sum0+=arrx[g];
sum1+=arrxd[g];
}
wc_=(sum1-sum0);
if((wc_<wc&&wc_>(-wc))||count>maxdd)
{
break;
}
sum0=sum1=0;
for(m=0;m<n;m++)
{
arrx[m]=arrxd[m];
}
}
if(count>maxdd)
{
printf("超出最大迭代次数\n");
}
else
{
printf("迭代次数为:%d\n",count);
printf("Jacobi迭代的结果为:");
for(k=0;k<n;k++)
{
printf("%lf ",arrxd[k]);
if(k==(n-1))
{
printf("\n");
}
}
}
return 0;
}
结果展示
高斯赛德尔迭代源代码
#include<stdio.h>
int main()
{
int n,i,j,k,b;
double num1;
double num2;
printf("输入方程组的维数:\n");
scanf("%d",&n);//方程组维数
double arrD[n][n]={0};//系数矩阵D
double arrB[n][n]={0};//矩阵D-1*(l+U)
printf("输入方程组各个系数:\n");
for(i=0;i<n;i++)
{ for(j=0;j<n;j++)
{ scanf("%lf",&num1);
arrD[i][j]=num1;
}
}
for(i=0;i<n;i++)
{ for(j=0;j<n;j++)
{ if(i==j)
{
arrB[i][j]=0;
}
else{
arrB[i][j]=-arrD[i][j]/arrD[i][i];
}
}
}
printf("输入矩阵b:\n");
double arrf[n]={0};//矩阵f
for(b=0;b<n;b++)
{
scanf("%lf",&num2);
arrf[b]=num2/arrD[b][b];
}
double sum=0;
double arrx[n]={0};//初始X向量
double arrxd[n]={0};
double wc;//误差限
int m;
int g;
double sum0=0;
double sum1=0;
double wc_=0;
int maxdd=0;
int count=0;
printf("输入精度:\n");
scanf("%lf",&wc);
printf("输入最大迭代次数:\n");
scanf("%d",&maxdd);
while(1)
{
count++;
for(k=0;k<n;k++)//迭代过程
{
for(j=0;j<n;j++)
{
sum=sum+arrB[k][j]*arrxd[j];
}
arrxd[k]=sum+arrf[k];
sum=0;
}//迭代过程
for(g=0;g<n;g++)
{
sum0+=arrx[g];
sum1+=arrxd[g];
}
wc_=(sum1-sum0);
if((wc_<wc&&wc_>(-wc))||count>maxdd)
{
break;
}
sum0=sum1=0;
for(m=0;m<n;m++)
{
arrx[m]=arrxd[m];
}
}
if(count>maxdd)
{
printf("超出最大迭代次数\n");
}
else
{
printf("迭代次数为:%d\n",count);
printf("高斯赛德尔迭代的结果为:");
for(k=0;k<n;k++)
{
printf("%lf ",arrxd[k]);
if(k==(n-1))
{
printf("\n");
}
}
}
return 0;
}
结果展示
代码改进的过程中我充分认识到了自身代码水平的不足,还得练!
————如有错误,欢迎指正(#^.^#)