数值计算方法实验报告.doc

《数值计算方法实验报告.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、数值计算方法实验报告班级：姓名：序号：1.实验内容（1）高斯消去法（2）全主元消去法（3）列主元消去法（4）Doolittle分解法（5）追赶法（6）Jacobi迭代法（7）Gauss-Seidel迭代法2.实验过程（源程序代码）（1）高斯消去法/*求解线性方程组--高斯消去法*/#includeusingnamespacestd;classCGauss{private:float**a;float*x;intn;public:CGauss(introw);~CGauss();voidInput();voidOutput();vo

2、idCalculate();};CGauss::CGauss(introw){n=row;a=newfloat*[n];x=newfloat[n];for(inti=0;i

3、>>a[i][j];}/*输出过程*/voidCGauss::Output(){cout<<"解向量为："<

4、p;a[n-1][n]=a[n-1][n]/a[n-1][n-1];for(i=0,k=n-2;k>=0;i++,k--){temp=0;for(j=k+1;j>n;CGaussobj(n);obj.Input();obj.Calculate();obj.Ou

5、tput();return0;}（1）全主元消去法/*求解线性方程组--全主元消去法*/#includeusingnamespacestd;classCGauss{private:float**a;float*x;int*z;/*z为记录列交换过程的数组下标，用在最后理顺解向量下标*/intn;public:CGauss(introw);~CGauss();voidInput();voidOutput();voidCalculate();};CGauss::CGauss(introw){n=row;a=newfloat*[n];x

6、=newfloat[n];z=newint[n];for(inti=0;i>a[i][j];}/*输出过程*

7、/voidCGauss::Output(){cout<<"解向量为："<temp){temp=a[i][j];row=

8、i;col=j;}}}if(0==temp){cout<<"系数矩阵为奇异阵！/n";return;}/*行