高斯消元解方程代码,C代码库,德仔网

代码语言

代码分类

【C】高斯消元解方程

作者: / 发布于2011/1/12/ 698

<div> 矩阵法求多元一次方程的解</div>
<div>
#include <iostream>
using namespace std;
int main()
{
 void input(float**p,int m,int n);
 int fun(float**p,int m,int n);
 float **q,b[30],w[30][30],x=1,d[30];
 int i,j,k,m,n;
 cout<<"请输入方程式的元数"<<endl;
 cin>>m;
 n=m;
 q=new float*[m];
 for(i=0;i<m;i++)
 q[i]=new float[m];
 cout<<"请输入方程的矩阵"<<endl;//输入各个元的系数
 input(q,m,n);
 cout<<"请输入各方程的值"<<endl;
 for(i=0;i<m;i++)
 cin>>d[i];
 for(i=0;i<m;i++)
 for(j=0;j<n;j++)
 w[i][j]=q[i][j];//设置一个替换数组，备用
 x=fun(q,m,n);//调用求矩阵的值
 cout<<"矩阵的值为"<<x<<endl;
 for(k=0;k<m;k++)
 {
 for(i=0;i<m;i++)
 for(j=0;j<n;j++)
 q[i][j]=w[i][j];//返回原数组
 for(i=0;i<m;i++)
 q[i][k]=d[i];//把第i列替换为方程的右得数
 b[k]=fun(q,m,n);
 }
 if(x==0) cout<<"此方程无解"<<endl;
 else 
 {
 for(i=0;i<m;i++)
 cout<<b[i]/x<<" ";//输出解
 cout<<endl;
 } 
 return 0;
}
<div>void input(float**p,int m,int n)
{
 for(int i=0;i<m;i++)
 for (int j=0;j<n;j++)
 cin>>p[i][j];
}
int fun(float**p,int m,int n)
{
 int i,j,k=1,t,c,e=0,x;
 int sort(float**p,int x,int m);
 for(i=0;i<m;i++)
 {
 if(p[i][i]==0)//调用查看对角线是否有为0的，如果有则与非0行替换
 {
 e++;//每两行交换则矩阵的值乘-1 见第72行
 x=i;
 sort(p,x,m);
 }
 for(j=i+1;j<n;j++)
 {
 c=-p[j][i]/p[i][i];
 for(t=0;t<m;t++)
 p[j][t]=p[j][t]+p[i][t]*c;
 }
 }
 for(i=0;i<m;i++)
 k=k*p[i][i];
 if(e%2==0) return k;
 else return -k;
}
int sort(float**p,int x,int m)
{
 int i,j,temp;
 for(i=0;i<m;i++)
 if(p[i][x]!=0)
 {
 for(j=0;j<m;j++)
 temp=p[x][j],p[x][j]=p[i][j],p[i][j]=temp;//交换
 break;
 }
 return **p;
}
<div>此代码适合与就任何方程，主要思想就是把矩阵先化为上三角，或下三角
</div>
<div></div>

试试其它关键字

同语言下

可能有用的

贡献的其它代码Label