欢迎来到天天文库
浏览记录
ID:14039084
大小:187.05 KB
页数:24页
时间:2018-07-25
《数值计算方法实验报告》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
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();vo
2、idInput();voidOutput();voidCalculate();};CGauss::CGauss(introw){n=row;a=newfloat*[n];x=newfloat[n];for(inti=0;i3、<>a[i][j];}/*输出过程*/voidCGauss::Output(){cout<<"解向量为:"<4、(j=k+1;j=0;i++,k--){temp=0;for(j=k+1;j5、<>n;CGaussobj(n);obj.Input();obj.Calculate();obj.Output();return0;}(1)全主元消去法/*求解线性方程组--全主元消去法*/#includeusingnamespacestd;classCGauss{private:float**a;float*x;int*z;/*z为记录列交换过程的数组下标,用在最后理顺解向量下标*/6、intn;public:CGauss(introw);~CGauss();voidInput();voidOutput();voidCalculate();};CGauss::CGauss(introw){n=row;a=newfloat*[n];x=newfloat[n];z=newint[n];for(inti=0;i7、nti=0;i>a[i][j];}/*输出过程*/voidCGauss::Output(){cout<<"解向量为:"<8、auss::Calculate(){/*循环变量*/inti,j,k;//消去过程floattemp;introw,col;for(k=0;ktemp){temp=a[i][j];row=i;col=j;}}}if(0==temp){cout<<"系数矩阵为奇异阵!/n";return;}/*行
3、<>a[i][j];}/*输出过程*/voidCGauss::Output(){cout<<"解向量为:"<4、(j=k+1;j=0;i++,k--){temp=0;for(j=k+1;j5、<>n;CGaussobj(n);obj.Input();obj.Calculate();obj.Output();return0;}(1)全主元消去法/*求解线性方程组--全主元消去法*/#includeusingnamespacestd;classCGauss{private:float**a;float*x;int*z;/*z为记录列交换过程的数组下标,用在最后理顺解向量下标*/6、intn;public:CGauss(introw);~CGauss();voidInput();voidOutput();voidCalculate();};CGauss::CGauss(introw){n=row;a=newfloat*[n];x=newfloat[n];z=newint[n];for(inti=0;i7、nti=0;i>a[i][j];}/*输出过程*/voidCGauss::Output(){cout<<"解向量为:"<8、auss::Calculate(){/*循环变量*/inti,j,k;//消去过程floattemp;introw,col;for(k=0;ktemp){temp=a[i][j];row=i;col=j;}}}if(0==temp){cout<<"系数矩阵为奇异阵!/n";return;}/*行
4、(j=k+1;j=0;i++,k--){temp=0;for(j=k+1;j5、<>n;CGaussobj(n);obj.Input();obj.Calculate();obj.Output();return0;}(1)全主元消去法/*求解线性方程组--全主元消去法*/#includeusingnamespacestd;classCGauss{private:float**a;float*x;int*z;/*z为记录列交换过程的数组下标,用在最后理顺解向量下标*/6、intn;public:CGauss(introw);~CGauss();voidInput();voidOutput();voidCalculate();};CGauss::CGauss(introw){n=row;a=newfloat*[n];x=newfloat[n];z=newint[n];for(inti=0;i7、nti=0;i>a[i][j];}/*输出过程*/voidCGauss::Output(){cout<<"解向量为:"<8、auss::Calculate(){/*循环变量*/inti,j,k;//消去过程floattemp;introw,col;for(k=0;ktemp){temp=a[i][j];row=i;col=j;}}}if(0==temp){cout<<"系数矩阵为奇异阵!/n";return;}/*行
5、<>n;CGaussobj(n);obj.Input();obj.Calculate();obj.Output();return0;}(1)全主元消去法/*求解线性方程组--全主元消去法*/#includeusingnamespacestd;classCGauss{private:float**a;float*x;int*z;/*z为记录列交换过程的数组下标,用在最后理顺解向量下标*/
6、intn;public:CGauss(introw);~CGauss();voidInput();voidOutput();voidCalculate();};CGauss::CGauss(introw){n=row;a=newfloat*[n];x=newfloat[n];z=newint[n];for(inti=0;i7、nti=0;i>a[i][j];}/*输出过程*/voidCGauss::Output(){cout<<"解向量为:"<8、auss::Calculate(){/*循环变量*/inti,j,k;//消去过程floattemp;introw,col;for(k=0;ktemp){temp=a[i][j];row=i;col=j;}}}if(0==temp){cout<<"系数矩阵为奇异阵!/n";return;}/*行
7、nti=0;i>a[i][j];}/*输出过程*/voidCGauss::Output(){cout<<"解向量为:"<8、auss::Calculate(){/*循环变量*/inti,j,k;//消去过程floattemp;introw,col;for(k=0;ktemp){temp=a[i][j];row=i;col=j;}}}if(0==temp){cout<<"系数矩阵为奇异阵!/n";return;}/*行
8、auss::Calculate(){/*循环变量*/inti,j,k;//消去过程floattemp;introw,col;for(k=0;ktemp){temp=a[i][j];row=i;col=j;}}}if(0==temp){cout<<"系数矩阵为奇异阵!/n";return;}/*行
此文档下载收益归作者所有