数值计算方法课程的设计.docVIP

数值计算方法课程设计 姓名 学号 成绩 课程实际报告 实验一：秦九韶算法 题目 用选列主元高斯消去法解线性方程组 算法语言： 利用c语言的知识编写该算法程序 算法步骤叙述： 秦九昭算法的基思路是v[0]=a[0]*x+a[1] v[i]=v[i-1]*x+a[i+1];利用秦九昭算法计算多项式函数。 程序清单： #includestdio.h void main() { float a[5],x,sum; int i; printf(presase input the value of x=); scanf(%f,x); for(i=5;i=0;i--) { printf(please input the value of a%d=,i); scanf(%f,a[i]); } sum=a[5]; for(i=5;i=1;i--) {sum=sum*x+a[i-1]; } printf(f(x)=%f/n,sum); } 输出结果计算： 实验总结： 通过运用C语言，解决了秦九韶算法手写的复杂。为以后的雪地打下基础。 实验二：用选列主元高斯消去法解线性方程组 题目 用选列主元高斯消去法解线性方程组 算法步骤叙述 第一步消元——在增广矩阵（A,b）第一列中找到绝对值最大的元素，将其所在行与第一行交换，再对（A,b）做初等行变换使原方程组的第一列元素除了第一行的全变为0； 第二步消元——在增广矩阵（A,b）中第二列中（从第二行开始）找到绝对值最大的元素，将其所在行与第二行交换，再对（A,b）做初等行变换使原方程组的第二列元素除了第一和第二行的全变为0； 第三步消元——在增广矩阵（A,b）中第三列中（从第三行开始）找到绝对值最大的元素，将其所在行与第三行交换，再对（A,b）做初等行变换使原方程组的第三列第四行元素为0； 第四，按x4-x3-x2-x1的顺序回代求解出方程组的解,x[n]=b[n]/a[n][n],x[i]=(b[i]-Σa[i][j]x[j])/a[i][i],i=n-1,…,2,1 程序清单： #includeiostream #includemath #define N 4 static double A[N][N] = {-3，-1,0,0，-1,2，-1,0，0，-1,2，-1，0,0，-1,2}; static double B[N]={1,0,0,0}; static double X[N]; int i,j,k; void main() { for(k = 0; k N-1 ;k++) { int index = k; for(i = k; i N ;i++) { if(fabs(A[index][k]) fabs(A[i][k])) { index = i; } } double temp; for( i = k ; i N ;i++ ) { temp = A[index][i]; A[index][i] = A[k][i]; A[k][i] = temp; } temp = B[index]; B[index] = B[k]; B[k] = temp; for(i = k+1; iN; i++) { double T = A[i][k]/A[k][k]; B[i] = B[i] - T * B[k]; for ( j = k+1 ; j N ; j++ ) { A[i][j] = A[i][j] - T * A[k][j]; } } } X[N-1] = B[N-1]/A[N-1][N-1]; for (i = N-2; i =0 ; i--) { double Temp = 0; for (int j = i+1; jN ;j++) Temp = Temp + A[i][j] * X[j]; X[i] = (B[i] - Temp) /A[i][i]; } cout 线性方程组的解（X1，X2，X3......Xn）为：endl; for( i = 0; i N ;i++) { cout X[i] ; } } 实验总结： 通过c++语言的编写过程掌握高斯消元法及选列主元元素的技术，掌握了简单的c++程序编辑语言编写算法程序 实验五：二分法与牛顿法 题目 用二分法和Newton迭代法求下列方程的正根：要求结果的误差限为 1.二分法 算法语言：C语言 算法思路