计算数值方法实验报告

实验地点 指导教师 于亚男 太原理工大学学生实验报告

学院名称 学生姓名 课程名称 计算机数值方法 专业班级 实验日期 实验题目 4月26日 学号 成绩 求解线性方程组 一、实验目的和要求 (1)了解直接法解线性方程组 (2)对Gauss消元法、LU分解法、追赶法有更深刻的理解 二、实验内容和原理 合理利用Gauss消元法、LU分解法、追赶法求解下列方程组： ?123??x1??14???x???8? 012① ????2?????241????x3????13???0.3?10?15?5.291②??11.2?1???4?8③??4??1227121??x1??59.17???????6.130?12??x2??46.78?? ?????x19523?????x2211???4???59.1435??x1???2??x???7?10???2???? 836??x3???7??????61120??x4???3??21??x1???7??121??x???5????2???④????????????（n=5,10,100,…） ??????x121???n?1???5??12??????5???xn???三．主要仪器设备 HP笔记本，VC++6.0 四．操作方法与实验步骤 1.高斯： #include void main() { int i,j,f,n; float a[5][5],b[5],x[5],l,k; printf(\ scanf(\ printf(\ for(i=0;i=0;i--) { x[i]=b[i]; k=0; for(f=i+1;f #include #define L 30 double a[L][L],b[L],l[L][L],u[L][L],x[L],y[L]; int main() { int n,i,j,k,r; printf(\ scanf(\ printf(\ for(i=1;i<=n;++i) { for(j=1;j<=n;++j) { scanf(\ } } printf(\ for(i=1;i<=n;++i) { scanf(\ } for(i=1;i<=n;++i) { for(j=1;j<=n;++j) { l[i][j]=0; u[i][j]=0.0; } } for(k=1;k<=n;++k) { for(j=k;j<=n;++j) { u[k][j]=a[k][j]; for(r=1;r0;--i) { x[i] = y[i]; for(j=i+1;j<=n;++j) { x[i]-=u[i][j]*x[j]; } x[i]/= u[i][i]; } for(i=1;i<=n;++i) { printf(\ } return 0; } 五．实验数据记录和处理 1.高斯