欢迎来到天天文库
浏览记录
ID:44819752
大小:62.82 KB
页数:5页
时间:2019-10-30
《实验三 高斯消去法》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、实验三高斯消去法一、实验目的(1)熟悉求解线性方程组的有关理论和方法;(2)能编程实现列主元高斯消去法;(4)根据不同类型的方程组,选择合适的数值方法。二、实验内容用选主元高斯消去求方程组三、算法基本原理A.;B.消元结果直接存储在系数矩阵中;C.当消元过程发生两行对调的情况为偶数次时,行列式值为对角线乘积,否则为对角线乘积的相反数。四、算法设计与实现列主元高斯消去法:列主元消元回代图4.2列主元的约当消去约当消去关键点://选主元素for(k=0;k2、;for(i=k+1;ifabs(max)){max=a[i][k];max_i=i;}if(max==0)break;if(max_i!=k)//交换两行for(j=k;j3、][j];}//消元}for(k=n-1;k>=0;k--){sum=0;for(j=k+1;j#includeusingnamespacestd;#defineMAX50voidinput(doublea[MAX][MAX+1],intn){cout<<"输入原方程组的增广矩阵"<4、or(inti=0;i>a[i][j];}voidoutput(doublex[],intn){cout<<"Gauss消去法得到的原方程组的解为"<5、的阶"<>n;input(a,n);for(k=0;kfabs(max)){max=a[i][k];max_i=i;}if(max==0)break;if(max_i!=k)//交换两行for(j=k;j6、a[i][k]=a[i][k]/-a[k][k];for(j=k+1;j=0;k--){sum=0;for(j=k+1;j
2、;for(i=k+1;ifabs(max)){max=a[i][k];max_i=i;}if(max==0)break;if(max_i!=k)//交换两行for(j=k;j3、][j];}//消元}for(k=n-1;k>=0;k--){sum=0;for(j=k+1;j#includeusingnamespacestd;#defineMAX50voidinput(doublea[MAX][MAX+1],intn){cout<<"输入原方程组的增广矩阵"<4、or(inti=0;i>a[i][j];}voidoutput(doublex[],intn){cout<<"Gauss消去法得到的原方程组的解为"<5、的阶"<>n;input(a,n);for(k=0;kfabs(max)){max=a[i][k];max_i=i;}if(max==0)break;if(max_i!=k)//交换两行for(j=k;j6、a[i][k]=a[i][k]/-a[k][k];for(j=k+1;j=0;k--){sum=0;for(j=k+1;j
3、][j];}//消元}for(k=n-1;k>=0;k--){sum=0;for(j=k+1;j#includeusingnamespacestd;#defineMAX50voidinput(doublea[MAX][MAX+1],intn){cout<<"输入原方程组的增广矩阵"<4、or(inti=0;i>a[i][j];}voidoutput(doublex[],intn){cout<<"Gauss消去法得到的原方程组的解为"<5、的阶"<>n;input(a,n);for(k=0;kfabs(max)){max=a[i][k];max_i=i;}if(max==0)break;if(max_i!=k)//交换两行for(j=k;j6、a[i][k]=a[i][k]/-a[k][k];for(j=k+1;j=0;k--){sum=0;for(j=k+1;j
4、or(inti=0;i>a[i][j];}voidoutput(doublex[],intn){cout<<"Gauss消去法得到的原方程组的解为"<5、的阶"<>n;input(a,n);for(k=0;kfabs(max)){max=a[i][k];max_i=i;}if(max==0)break;if(max_i!=k)//交换两行for(j=k;j6、a[i][k]=a[i][k]/-a[k][k];for(j=k+1;j=0;k--){sum=0;for(j=k+1;j
5、的阶"<>n;input(a,n);for(k=0;kfabs(max)){max=a[i][k];max_i=i;}if(max==0)break;if(max_i!=k)//交换两行for(j=k;j6、a[i][k]=a[i][k]/-a[k][k];for(j=k+1;j=0;k--){sum=0;for(j=k+1;j
6、a[i][k]=a[i][k]/-a[k][k];for(j=k+1;j=0;k--){sum=0;for(j=k+1;j
此文档下载收益归作者所有