实验三 编程实现列主元高斯消去法.docVIP

实验三 编程实现列主元高斯消去法 实验目的：实现高斯主消元法，对计算过程加深理解。 实验内容： 编写c++程序，实现对角占优方程组的编程求解。 3. 实验步骤 1、设计方程组的存储为二位数组，最大方程组数为100，第i行j列的元素值代表第i个方程的第j个系数，输入时没有的系数项填0。 2 对于每个方程按主消元法从0~n-1依次消元。 3迭代求解，对于第i个未知数的值，依次迭代i+1~n-1已求出的结果 4．实验结果分析： 对于除数很小的情况程序不能很好的解决，对于不是对角占优的也不能搜索出主元素，以后在进一步解决上述问题。 #includeiostream using namespace std; const int MAX=100; double str[MAX][MAX]; int main() { while(1) { cout输入方程组个数 0n100endl; int n,m,i,j,k,flag=0; cinn; if(n=0||n=100)break; cout输入方程组的系数endl; for(i=0;in;i++) { for(j=0;j=n;j++) { cinstr[i][j]; if(i==jstr[i][j]==0)flag=1; } } if(flag==1){cout方程组不是对角占优的，此程序不能解决endl;system(pause);return 0;} for(i=0;in;i++) { for(j=n;j=i;j--)//对角变一 str[i][j]/=str[i][i]; for(k=i+1;kn;k++)//方程消元 { double tem=str[k][i]/str[i][i]; for(j=i;j=n;j++) { str[k][j]-=str[i][j]*tem; } } } double ans[MAX]={0}; ans[n-1]=str[n-1][n]; for(i=n-2;i=0;i--) { ans[i]=str[i][n]; for(j=n-1;ji;j--) ans[i]-=str[i][j]*ans[j]; } for(i=0;in;i++) coutXi+1 = ans[i]endl; system(pause); } }