数值计算基础实验指导书及代码2013.docVIP

数值计算基础 实验指导书 2013年 目录 实验一 直接法解线性方程组的 1 实验二 插值方法 4 实验三 数值积分 6 实验四 常微分方程的数值解 8 实验五 迭代法解线性方程组与非线性方程 10 实验一 直接法解线性方程组 一、实验目的 掌握列选主元消去法与追赶法解线性方程组。 二、实验内容 分别写出Guass列选主元消去法与追赶法的算法，编写程序上机调试出结果，要求所编程序适用于任何一解线性方程组问题，即能解决这一类问题，而不是某一个问题。实验中以下列数据验证程序的正确性。 1、用Guass列选主元消去法求解方程组 2、用追赶法求解方程组 三、实验仪器设备与材料 主流微型计算机 四、实验原理 1、Guass列选主元消去法 对于AX =B 1）、消元过程：将（A|B）进行变换为，其中是上三角矩阵。即： k从1到n-1 列选主元 选取第k列中绝对值最大元素作为主元。 换行 归一化 消元 2）、回代过程：由解出。 #includeiostream.h #include stdio.h #include stdlib.h #include math.h #define N 3 #define EPS 1.0e-4 double a[N][N],b[N],c[N][N],d[N][N]; void swap(double a,double b)//换行 { double temp; temp=a; a=b; a=temp; } void gauss(double a[][N],double b[N])//gauss求解 { int i,j,k,l; double temp,m; for(k=0;k=N-2;k++) { temp=fabs(a[k][k]); //选主元 l=k; for(i=k+1;iN;i++) { if(fabs(a[i][k]temp))//fabs为取绝对值函数。包含在头文件math.h中 { temp=fabs(a[i][k]); l=i; } } if(fabs(a[l][k])EPS) { printf(方程有无数组解!); exit(0); } if(l!=k) //换行 { for(j=k;j=N-1;j++) swap(a[l][j],a[k][j]); swap(b[l],b[k]); } for(i=k+1;iN;i++) //按行进行消元 { m=a[i][k]/a[k][k];//行乘数 for(j=k;jN;j++) a[i][j]-=m*a[k][j]; b[i]-=m*b[k]; } } if(a[N-1][N-1]==0) //矩阵行列式值为0 { printf(矩阵的行列式为0!\n); exit(-1); } b[N-1]/=a[N-1][N-1]; for(i=N-2;i=0;i--) { for(j=N-1;ji;j--) b[i]-=a[i][j]*b[j]; b[i]/=a[i][i]; } } void main() { int i; double a[][N]={{2.5,2.3,-5.1},{5.3,9.6,1.5},{8.1,1.7,-4.3}}; double b[]={3.7,3.8,5.5}; gauss(a,b); //高斯求解法 for(i=0;iN;i++) //输出结果 printf(X%d= %.4f\n,i+1,b[i]); } 2、追赶法 线性方程组为： 做LU分解为： 分解公式： 则 回代公式： #includestdio.h #define N 5 double a[N][N],f[N],u[N][N],l[N][N],y[N]; void zhuigan(double a[][N],double f[N]){ int i;double x[N]={0,0,0,0}; l[0][0]=a[0][0]; for(i=0;iN;i++){ if(l[i][i]!=0) break; else printf(无解！); } for(i=1;i=N-1;i++){ l[i][i-1]=a[i][i-1]; l[i][i]=a[i][i]-l[i][i-1]*u[i-1][i];