c語(yǔ)言實(shí)現(xiàn)矩陣的加減乘除求逆運(yùn)算

1、#include#include#definecol3#definerow3classmatrix/類(lèi)的定義(private:doublemcolrow;/矩陣設(shè)置為私有的，public:matrix()/無(wú)參數(shù)的構(gòu)造函數(shù)matrix(doubleacolrow);/有參數(shù)的構(gòu)造函數(shù)matrixAdd(matrix&b);/加法運(yùn)算聲明matrixSub(matrix&b);/減法運(yùn)算聲明matrixMul(matrix&b);/乘法運(yùn)算聲明matrixDiv(matrix&b);/除法運(yùn)算聲明matrixInverse。;/求逆運(yùn)算聲明matrix();/析構(gòu)

2、函數(shù)聲明 voiddisplay。;/顯示函數(shù)聲明;matrix:matrix(doubleacolrow)/構(gòu)造函數(shù)的定義inti,j;for(i=0;icol;i+)for(j=0;jrow;j+)mij=aij;matrixmatrix:Add(matrix&b)/加法運(yùn)算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)/減法運(yùn)算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)/乘法運(yùn)算(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)/除法運(yùn)算(除法直接求解，參見(jiàn)主函數(shù)matrixc;return

4、(c);)matrixmatrix:Inverse()求逆運(yùn)算(http:/ bkk,貌似是不正確的，因?yàn)檫@對(duì)比如說(shuō)是0,0,1,1,0,1,0,1,1的矩陣就會(huì)判斷為不可逆，而實(shí)際上該矩陣是可逆的，這里應(yīng)該是作者筆誤，待進(jìn)一步求證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;/*導(dǎo)出結(jié)果*/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 請(qǐng)選擇要進(jì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)/矩陣加法運(yùn)算printf(矩陣加法運(yùn)算結(jié)果：n);mc=ma.Add(mb);mc.display();elseif(flag=3)/矩陣減法運(yùn)算printf(矩陣減法運(yùn)算結(jié)果：n);mc=ma.Sub(mb);mc.display();)elseif(flag=4)/矩陣乘法運(yùn)算(printf(矩陣乘法運(yùn)算結(jié)果：n);mc=ma.Mul(mb);mc.display();)elseif(

7、flag=5)/矩陣除法運(yùn)算(printf(矩陣除法運(yùn)算結(jié)果：n);printf(矩陣的除法分成兩類(lèi)：n1、AB=inverse(A)*Bn2、B/A=B*inverse(A)n);printf(采用第 1 類(lèi)，則 ab 的結(jié)果為：n);mc=ma.Inverse();mc=mc.Mul(mb);mc.display();printf(采用第 2 類(lèi)，則 a/b 的結(jié)果為：n);mc=mb.Inverse();mc=ma.Mul(mc);mc.display();)elseif(flag=6)/矩陣求逆運(yùn)算(printf(矩陣 a 求逆運(yùn)算結(jié)果為：n);mc=ma.Inverse();mc.display();printf(矩陣 b 求逆運(yùn)算結(jié)果為：n);mc=mb.Inverse();mc.displa