计算方法实验指导书(有例程).docVIP

计算 理学院实验中心 数学专业实验室编写 实验一 【实验类型】 性 【实验学时】 学时 【实验内容】 能够利用列主元的高斯消去法解任意阶数的线性方程组 【实验前的预备知识】 1．计算机基础知识； 2．熟悉编程基本思想； 3．熟悉常见数学函数【实验方法或步骤】 列主元消去法基本思路 设有线性方程组，设是可逆矩阵。列主元消去法的基本思想就是通过列主元的选取将初等行变换作用于方程组的增广矩阵，将其中的变换成一个上三角矩阵，然后求解这个三角形方程组。 列主元高斯消去法算法描述 将方程组用增广矩阵表示。 步骤1：消元过程，对 选主元，找使得 如果，则矩阵奇异，程序结束；否则执行（3）； 如果，则交换第行与第行对应元素位置，， ； 消元，对，计算对，计算 步骤 2：回代过程： 若则矩阵奇异，程序结束；否则执行（2）； 对，计算 [实验内容] 在课后习题中选择一个求解线性方程组的题编程计算。交回实验报告与计算结果 [实验程序] (由学生独立完成) 参考程序： 列主元消去法 #include math.h void zhuyuan (k,n,a) int k,n; float a[5][6]; {int t,i,j; float p,q; p=fabs(a[k][k]); t=k; for (i=k+1;i=n;i++) if (fabs(a[i][k])p) {p=fabs(a[i][k]); t=i; } for (j=k;j=n+1;j++) {q=a[k][j]; a[k][j]=a[t][j]; a[t][j]=q; } } void shuru(n,a) int n; float a[5][6]; {int i,j; printf(input a[i][j]:\n); for(i=0;i=n;i++) {for(j=0;j=n+1;j++) {scanf(%f,%a[i][j]); printf(%8.5f,a[i][j]); } printf(\n); } } void xiaoyuan(n,a) int n; float a[5][6]; {int(k=0;kn;k++) {zhuyuan(k,n,a); for (i=k+1;i=n;i++) for (j=k+1;j=n+1;j++) a[i][j]=a[i][j]-a[i][k]*a[k][j]/a[k][k]; } } void huidai(n,a,x) int n; float a[5][6],x[5]; {int k,j; x[n]=a[n][n+1]/a[n][n]; for (k=n-1;k=0;k--) {x[k]=a[k][n+1]; for (j=k+1;j=n;j++) x[k]=x[k]-a[k][j]*x[j]; x[k]=x[k]/a[k][k]; } } main() {int n,i; float a[5][6],x[5]; printf(input n:\nn=); scanf(%d,n); printf(%d\n,n); n=n-1; shuru(n,a); huidai(n,a,x); for (i=0;i=n;i++) printf(x[%d]=%f\n,i+1,x[i]); printf(\n); } /*说明：本程序只能求解方程个数与未知量个数相同，且系数行列式不为零的线形方程组。*/ 例 求解线形方程组： (书118页 2(1) ) 此程序运行结果如下： input n: n=4 input a[i][j]: 1.00000 0.83240 0.76750 0.98310 8.89970 0.83240 0.69300 0.64000 0.81900 7.41440 0.76750 0.64000 0.59110 0.75800 6.85280 0.98310 0.81900 0.75800 0.00550 4.91710 x[1]=1.000094 x[2]=1.999953 x[3]=2.999929 x[4]=4.000000 实验二 最小二乘法 【实验类型】性 【实验学时】 学时 【实验内容】 曲线拟合的最小二乘法的能根据给定的函数值表构造出次数不相同的拟合多项式。 【实验