c语言实现矩阵的加减乘除求逆运算

1、#include#include#definecol3#definerow3classmatrix/类的定义(private:doublemcolrow;/矩阵设置为私有的，public:matrix()/无参数的构造函数matrix(doubleacolrow);/有参数的构造函数matrixAdd(matrix&b);/加法运算声明matrixSub(matrix&b);/减法运算声明matrixMul(matrix&b);/乘法运算声明matrixDiv(matrix&b);/除法运算声明matrixInverse。;/求逆运算声明matrix();/析构

2、函数声明 voiddisplay。;/显示函数声明;matrix:matrix(doubleacolrow)/构造函数的定义inti,j;for(i=0;icol;i+)for(j=0;jrow;j+)mij=aij;matrixmatrix:Add(matrix&b)/加法运算inti,j;matrix*c=(matrix*)malloc(sizeof(matrix);for(i=0;icol;i+)for(j=0;jmij=mij+b.mij;return(*c);matrixmatrix:Sub(matrix&b)/减法运算inti,j;matrix*c=(matrix*

3、)malloc(sizeof(matrix);for(i=0;icol;i+)for(j=0;jmij=mij-b.mij;return*c;matrixmatrix:Mul(matrix&b)/乘法运算(inti,j,k;doublesum=0;matrix*c=(matrix*)malloc(sizeof(matrix);for(i=0;icol;i+)(for(j=0;jrow;j+)(for(k=0;kmij=sum;sum=0;)return(*c);)matrixmatrix:Div(matrix&b)/除法运算(除法直接求解，参见主函数matrixc;return

4、(c);)matrixmatrix:Inverse()求逆运算(http:/ bkk,貌似是不正确的，因为这对比如说是0,0,1,1,0,1,0,1,1的矩阵就会判断为不可逆，而实际上该矩阵是可逆的，这里应该是作者笔误，待进一步求证for(intj=0;j=i;j-)bij/=bii;for(k=0;kM;k+)if(k!=i)doubletemp=bki;for(j=0;jN;j+)bkj-=temp*bij;/*导出结果*/for(i=0;iM;i+)for(j=3;jmij-3=bij;return(*c);matrix:matrix()()voidmatrix:display()(in

5、ti,j;for(i=0;icol;i+)(for(j=0;jrow;j+)printf(%f,mij);printf(n);)voidmain()(doublea33=1,0,1,0,1,1,0,3,1;doubleb33=0,0,1,1,0,1,0,1,0;matrixma(a),mb(b),mc;intflag;printf(n 请选择要进行的操作：n1、打印t2、加法);printf(t3、减法n4、乘法t5、除法t6、求逆n7、退出n);printf(n);scanf(%d,&flag);while(flag=1)|(flag=2)|(flag=3)|(flag=4)|(fl

6、ag=5)|(flag=6)|(flag=7)if(flag=1)printf(矩阵 a 为：n);ma.display();printf(矩阵 b 为：n);mb.display();if(flag=2)/矩阵加法运算printf(矩阵加法运算结果：n);mc=ma.Add(mb);mc.display();elseif(flag=3)/矩阵减法运算printf(矩阵减法运算结果：n);mc=ma.Sub(mb);mc.display();)elseif(flag=4)/矩阵乘法运算(printf(矩阵乘法运算结果：n);mc=ma.Mul(mb);mc.display();)elseif(

7、flag=5)/矩阵除法运算(printf(矩阵除法运算结果：n);printf(矩阵的除法分成两类：n1、AB=inverse(A)*Bn2、B/A=B*inverse(A)n);printf(采用第 1 类，则 ab 的结果为：n);mc=ma.Inverse();mc=mc.Mul(mb);mc.display();printf(采用第 2 类，则 a/b 的结果为：n);mc=mb.Inverse();mc=ma.Mul(mc);mc.display();)elseif(flag=6)/矩阵求逆运算(printf(矩阵 a 求逆运算结果为：n);mc=ma.Inverse();mc.display();printf(矩阵 b 求逆运算结果为：n);mc=mb.Inverse();mc.displa