线性方程组的迭代解法
线性方程组的迭代解法 C语言实现
至少选择雅可比迭代法、高斯-赛德尔迭代法上述方法中的一种方法(鼓励多选)对下列方程组进行求解。
精度要求10-3
//雅可比迭代法:
#include<stdio.h>
#include<math.h>
double fx1(double x2,double x3){
double x1=0.1*x2+0.2*x3+0.72;
return x1;
}
double fx2(double x1,double x3){
double x2=0.1*x1+0.2*x3+0.83;
return x2;
}
double fx3(double x1,double x2){
double x3=0.2*x1+0.2*x2+0.84;
return x3;
}
double max(double x1,double x2,double x3,double a,double b,double c){
double max=0,m;
if(fabs(a-fx1(x2,x3))>fabs(b-fx2(x1,x3))){
m=fabs(a-fx1(x2,x3));
} else m=fabs(b-fx2(x1,x3));
if(m>fabs(c-fx3(x1,x2))){
max=m;
}else max=fabs(c-fx3(x1,x2));
return max;
}
int main(){
double x1,x2,x3,a,b,c;
int n=0;
printf("输入初值x1,x2,x3\n");
scanf("%f %f %f",&x1,&x2,&x3);
a=x1;
b=x2;
c=x3;
while(max(x1,x2,x3,a,b,c)>0.001){
a=fx1(x2,x3);
b=fx2(x1,x3);
c=fx3(x1,x2);
x1=a;
x2=b;
x3=c;
n++;
printf("%f\t%f\t%f\t\n",x1,x2,x3);
}
printf("迭代次数:%d\n",n);
printf("结果为 x1=%f\t x2=%f\t x3=%f\t",x1,x2,x3);
}
//高斯-赛德尔迭代法:
#include<stdio.h>
#include<math.h>
double fx1(double x2,double x3){
double x1=0.1*x2+0.2*x3+0.72;
return x1;
}
double fx2(double x1,double x3){
double x2=0.1*x1+0.2*x3+0.83;
return x2;
}
double fx3(double x1,double x2){
double x3=0.2*x1+0.2*x2+0.84;
return x3;
}
double max(double x1,double x2,double x3){
double max=0,m;
if(fabs(x1-fx1(x2,x3))>fabs(x2-fx2(x1,x3))){
m=fabs(x1-fx1(x2,x3));
} else m=fabs(x2-fx2(x1,x3));
if(m>fabs(x3-fx3(x1,x2))){
max=m;
}else max=fabs(x3-fx3(x1,x2));
return max;
}
int main(){
double x1,x2,x3,a,b,c;
int n=0;
printf("输入初值x1,x2,x3\n");
scanf("%f %f %f",&x1,&x2,&x3);
while(max(x1,x2,x3)>0.001){
x1=fx1(x2,x3);
x2=fx2(x1,x3);
x3=fx3(x1,x2);
n++;
printf("%f\t%f\t%f\t\n",x1,x2,x3);
}
printf("迭代次数:%d\n",n);
printf("结果为 x1=%f\t x2=%f\t x3=%f\t",x1,x2,x3);
}
第一次写博客,代码是自己编写的,仅供参考。