高斯约当消去法.docx

《高斯约当消去法.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、姓名吴娴学号2011210927班级2011211123高斯约当消去法#include#include#include/*形成R矩阵函数*/voidxz(intm,floatres[10],intm1[],intm2[],floatr[][10]){inti,j,k;for(i=0;i0){r[k][k]=r[k][k]+res[i];//自电阻形成R值r[j][k]=r[j][k]-res[i];//互

2、电阻形成R值r[k][j]=r[j][k];//矩阵为对称矩阵}r[j][j]=r[j][j]+res[i];//网孔独有电阻仅算自电阻}}/*选列主元函数*/floatzu(floatalfa[],floateps,intny,intnp,inti,floatra[][11]){intir/*计数器*/,m,ic;floatp;alfa[0]=fabs(ra[i][i]);m=i;for(ir=i;ir<=ny;ir++){if(fabs(ra[ir][i])<=alfa[0])continue;alfa[0]=fabs(ra[ir][i])

3、;//如果该元素大于原主元，则更新主元m=ir;//确定主元的所在行数}if(alfa[0]<=eps)return(0.0);//若该列所有主元均近似为零，则该方程无唯一解alfa[0]=ra[m][i];//更新归一因子if(m==1)return(1.0);for(ic=i;ic<=np;ic++){p=ra[m][ic];ra[m][ic]=ra[i][ic];ra[i][ic]=p;}//交换原主元与现主元所在行的所有电阻系数return(1.0);}/*求解方程的Gauss函数*/floatgauss(floatra[][11],f

4、loatvn[],intny){inti,j,k,np;floatdeta,eps,beta/*消元因子*/,alfa[2];//归一因子eps=1.e-10;np=ny+1;//包括增广部分for(i=1;i<=ny;i++){deta=1.0;deta=zu(alfa,eps,ny,np,i,ra);if(deta==0.0)return(deta);//电路无唯一解for(j=i;j<=np;j++)ra[i][j]=ra[i][j]/alfa[0];//主元归一for(k=1;k<=ny;k++){if(k==i)continue;be

5、ta=ra[k][i];//对于非对角线电阻系数来确定消元因子for(j=i;j<=np;j++)ra[k][j]=ra[k][j]-beta*ra[i][j];//消元}}for(i=1;i<=ny;i++)vn[i]=ra[i][np];//电流解即为增广部分的值return(deta);}intmain(){floatres[20]/*存储第i个电阻值*/,r[10][10]/*存储电阻矩阵rij*/,ra[10][11];//增广矩阵，最大限度10列floatamps[10]/*电流矩阵*/,v[10]/*增广（电压值）部分*/,det

6、a;//消元因子inti,j,m,n,m1[20]/*与第i个网孔关联的电阻*/,m2[20];//与第i个电阻相关联的相邻网孔的序号FILE*InputFile,*OutputFile;//文件流v[0]=0;/*输入数据*/InputFile=fopen("Input.txt","r");//打开文件if(InputFile==NULL){printf("CannotopentheInputFile!");return0;}fscanf(InputFile,"%d%d",&n,&m);//读取网孔数，电阻数for(i=1;i

7、++)fscanf(InputFile,"%f",&v[i]);//读取增广部分，即等式右边部分for(i=0;i

8、++)ra[i][j]=r[i][j];ra[i][n+1]=v[i];//补上增广部分}//形成增广矩阵/*打印RA矩阵*/OutputFile=f