基本矩阵运算.doc

/* 这个类提供的是一些基本的矩阵运算 可以满足+,-,*,det,R以及取反矩阵的运算 但是+,-,*,det,都没有对运算的可行性进行检查 同时，注意不要让秩超过整数 */ //CMatrix:: class CMatrix { public: CMatrix(); CMatrix(int xM,int xN); virtual ~CMatrix(); int getm()const; int getn()const; int getij(int i,int j)const;//下标从0,0开始！ void setij(int xData,int i,int j); void reset(int xM,int xN); void operator=(const CMatrix xa); CMatrix operator-()const; CMatrix operator+(const CMatrix xb)const;//不进行合法性检查 CMatrix operator-(const CMatrix xb)const;//不进行合法性检查 CMatrix operator*(const CMatrix xb)const;//不进行合法性检查 int R()const; int det()const;//不定义m!=n时的行列式，返回未知值 CMatrix Anti()const; private: int mM,mN; int **mData; }; //The math function that has been used //****************************************************************** int min(int a,int b); int gcd(int a,int b); int abs(int a); int swap(int a,int b); void FMswap(double *a,double *b); void FMbeishu(double *a,double k,int n); void FMxiao(double *a,double *b,double k,int n); int Fdet(int **a,int n); int FR(int **a,int m,int n); int min(int a,int b) { if(ab) return b; return a; } int gcd(int a,int b) { if(b==0) return a; else if(a==1) return b; else return (b,a%b); } int abs(int a) { if(a0)return -a; return a; } int swap(int a,int b) { int t; t=a; a=b; b=t; return 0; } void swap(double a,double b) { double t; t=a; a=b; b=t; } void FMswap(double *a,double *b) { double *t; t=a; a=b; b=t; } void FMbeishu(double *a,double k,int n) { for(int i=0; in; i++) a[i]=a[i]*k; } void FMxiao(double *a,double *b,double k,int n) { for(int i=0; in; i++) b[i]+=a[i]*k; } int Fdet(int **a,int n) { bool ansget=false; int i,j,k; double ans; double **x; x=new double*[n]; for(i=0; in; i++) x[i]=new double[n]; for(i=0; in; i++) for(j=0; jn; j++) x[i][j]=double(a[i][j]); for(j=0; jn; j++) { int notzerox=j; for(i=j; in; i++) if(x[i][j]!=0){ notzerox=i; break; } for(k=j; kn; k++) swap(x[notzerox][k],x[j][k]); if(x[j][j]==0){ ansget=true; ans