太原理工大学数值计算方法实验报告
- 格式:doc
- 大小:161.00 KB
- 文档页数:10
printf("%lf ",x[i]);
printf("\n");
}
int _tmain(int argc, _TCHAR* argv[])
{
double a[3][4],x[3]={0,0,0},d[3];
for(int i=0;i<3;i++)
for(int j=0;j<4;j++)
scanf("%lf",&a[i][j]);
shuchu(x);
do
{
x[0]=(a[0][3]-a[0][1]*x[1]-a[0][2]*x[2])/a[0][0];
x[1]=(a[1][3]-a[1][0]*x[0]-a[1][2]*x[2])/a[1][1];
x[2]=(a[2][3]-a[2][0]*x[0]-a[2][1]*x[1])/a[2][2];
d[0]=abs((a[0][3]-a[0][1]*x[1]-a[0][2]*x[2])/a[0][0]-x[0]);
d[1]=abs((a[1][3]-a[1][0]*(a[0][3]-a[0][1]*x[1]-a[0][2]*x[2])/a[0][0]-a[1][2]*x[2]) /a[1][1]-x[1]);
d[2]=abs((a[2][3]-a[2][0]*(a[0][3]-a[0][1]*x[1]-a[0][2]*x[2])/a[0][0]-a[2][1]*(a[1] [3]-a[1][0]*(a[0][3]-a[0][1]*x[1]-a[0][2]*x[2])/a[0][0]-a[1][2]*x[2])/a[1][1])/a[2][2]-x[2]);
shuchu(x);
}while(d[0]>0.5e-5&&d[1]>0.5e-5&&d[2]>0.5e-5);
system("pause");
return 0;
}
实验结果与分析
1.列主元素消元法
2.完全组元素消元法
3.LU分解法
4.高斯-赛德尔迭代法
讨论、心得(可选):
了解Gauss消元法、LU分解法、追赶法等线性方程组直接求解的基本方法、基本原理;能够按照工程实际要求,选择适当的算法;通过编写程序,进行算法设计和数值求解,了解雅可比迭代法、高斯-赛德尔迭代法等线性方程组迭代求解的基本方法、基本原理,能够按照工程实际要求,选择适当的算法,通过编写程序,进行算法设计和数值求解。