高斯消元法c++代码

function.h
#ifndef FUNCTION_H
#define FUNCTION_H
#include<iostream>
#include<iomanip>
#include<cmath>
using namespace std;
class juZhen{
private:
double **m_p;
cin>>p[i];
return p;
}
void main(){
int pp(1);
while(pp){
cout<<"**********************************************************"<<endl;
cout<<"* 高斯消元法解方程组 *"<<endl;
}
}
void juZhen::display(){ //显示矩阵
cout.setf(ios::right);
for(int x=0;x<m_sizex;x++){
for(int y=0;y<m_sizey;y++){
cout<<setw(10)<<m_p[x][y]<<setw(10);
cout<<"* ⑤退出当前矩阵运算 *"<<endl;
cout<<"* *"<<endl;
cout<<"* ⑥退出程序 *"<<endl;
p={
cout<<"方程组无确定解！"<<endl;
break;
}
else{
double t;
for(int y=n;y<m_sizey;y++){
t=m_p[n][y];
m_p[n][y]=m_p[p][y];
void columnPrincipalComponentEliminationMethod(); //列主元消去法函数
void principalComponentEliminationMethod(); //全主元消去法
void result(); //输出结果函数
for(int i=0;i<m_sizex;i++)
delete[] m_p[i];
}
void display(); //显示函数
void eliminationGaussian(); //高斯消去法函数
}
cout<<endl;
}
}
void juZhen::eliminationGaussian(){ //高斯直接消去法
for(int n=0;n<m_sizey-2;n++){
for(int x=n;x<m_sizex-1;x++){
double d1=m_p[x+1][n]/m_p[n][n];
}
m=(m_p[k][m_sizey-1]-sum2)/m_p[k][k];
sum1[k]=m;
cout<<"**********************************************************"<<endl;
cout<<"X("<<k+1<<")="<<m<<endl;
}
}
}
}
display();
result();
}
void juZhen::principalComponentEliminationMethod(){ //全主元消去法
for(int n=0;n<m_sizey-2;n++){
}
for(int x=0;x<m_sizex;x++){ //列变换
t=m_p[x][n];
m_p[x][n]=m_p[x][p2];
m_p[x][p2]=t;
}
for(int x=n;x<m_sizex-1;x++){
double d1=m_p[x+1][n]/m_p[n][n];
}
delete[] sum1;
}
}
main.cpp
#include"function.h"
using namespace std;
double* input(){
double *p;
int size;
cout<<"请输入矩阵元素的总和："<<endl;
for(int y=n;y<m_sizey;y++){
m_p[x+1][y]=m_p[x+1][y]-d1*m_p[n][y];
}
}
}
display();
result();
}
void juZhen::result(){
int m_sizex; //矩阵的行数
int m_sizey; //矩阵的列数
public:
juZhen(double *p=0,int sizex=0,int sizey=0); //构造函数
~juZhen(){ //析构函数
for(int y=n;y<m_sizey;y++){
m_p[x+1][y]=m_p[x+1][y]-d1*m_p[n][y];
}
}
}
display();
result();
}
void juZhen::columnPrincipalComponentEliminationMethod(){ //列主元消去法
cin>>size;
p=new double[size];
cout<<"请输入矩阵元素："<<endl;
for(int i=0;i<size;i++)
cin>>p[i];
return p;
}
#include"function.h"
using namespace std;
};
#endif
function.cpp
#include"function.h"
juZhen::juZhen(double *p,int sizex,int sizey){ //构造函数
m_sizex=sizex;
m_sizey=sizey;
m_p=new double*[m_sizex];
cout<<"X("<<m_sizey-1<<")="<<m<<endl;
for(int k=m_sizex-2;k>=0;k--){
double sum2(0);
for(int y=k+1;y<m_sizey-1;y++){
sum2+=m_p[k][y]*sum1[y];
if(m_p[m_sizex-1][m_sizey-2]==0)
cout<<"方程组无确定解！"<<endl;
else{
cout<<"**********************************************************"<<endl;
cout<<"该方程组解如下："<<endl;
cout<<"* *"<<endl;
cout<<"* ①直接消元法 *"<<endl;
cout<<"* *"<<endl;
m_p[p][y]=t;
}
for(int x=n;x<m_sizex-1;x++){
double d1=m_p[x+1][n]/m_p[n][n];
for(int y=n;y<m_sizey;y++){
m_p[x+1][y]=m_p[x+1][y]-d1*m_p[n][y];
for(int i=0;i<sizex;i++)
m_p[i]=new double[m_sizey];
for(int x=0;x<m_sizex;x++){
for(int y=0;y<m_sizey;y++){
m_p[x][y]=p[x*m_sizey+y];
}
p1=x;
p2=y;
}
}
}
double t;
for(int y=n;y<m_sizey;y++){ //行变换