矩陣的運算程序C語言版

1、#include<stdio.h>#defineM20#defineN20floatAMN;floatBMN;floatCMN;inti,j,m,n,p,q;floaty=1.0;voidmain()(printf("#'n”);printf("#歡迎您使用矩陣函數(shù)包系統(tǒng)。print*”#系統(tǒng)功能：printf("#計算兩個矩陣相加，相減；printf("#數(shù)乘矩陣；矩陣轉(zhuǎn)置；求逆#榭#n");矩printf("#陣制用的M”)；printf("#、n”);printf("請選擇您需要的運算，按

2、回車鍵結(jié)束n”)；printf*n").printf("1,矩陣相加n”);,printf("2矩陣相減n")；,printf("3數(shù)乘矩陣n")；,printf("4矩陣轉(zhuǎn)置n");.printf("5逆矩陣n”)；printf("6,退出n”)；print%'intx;seanf("%d",&x);switch(x)(case 1: printf("請輸入矩陣A的行數(shù)和列數(shù)(用逗號隔開):")；seanf("%d,%d&quo

3、t;,&i,&j);printf("請輸入矩陣B的行數(shù)和列數(shù)(用逗號隔開):")；seanf(”d,%d”,&m,&n);if(i!=m|j!=n)printf("*對不起，您輸入的兩個矩陣不能相加，請重試.*S”)；elseprintf("請輸入矩陣A:fi");for(p=0;p<i;p+)for(q=0;qvj;q+)seanf("%f',&Apq);printf("輸出矩陣A：n");for(p=0;p<i;p+)for(q=0;q<j;q+

4、)(printf("%10.2r,Apq);if(q+1)%j=0)printf("n");printf("請輸入矩陣B：n");for(p=0;p<i;p+)for(q=0;q<j;q+)scanf(“r,&Bpq);printf("輸出矩陣B：n");for(p=0;p<i;p+)for(q=0;q<j;q+)pnntf(H%10.2r,Bpq);if(q+1)%j=0)printf(nnH);>printf(H矩陣A+矩陣B為:n“);計算兩個矩陣相加for(p=0;p<i;

5、p+)for(q=0;q<j;q+)Cpq=Apq+Bpq;for(p=0;p<i;p+)for(q=0;q<j;q+)printf("%10.2f',Cpq);if(q+1)%j=0)printf("n");>break;case 2:(printf("請輸入矩陣A的行數(shù)和列數(shù)(用逗號隔開):")；scanf("%d,%d",&i,&j);printf("請輸入矩陣B的行數(shù)和列數(shù)(用逗號隔開):")；scanf("%d,%d",&

6、;m,&n);if(i!=m|j!=n)printf("*對不起，您輸入的兩個矩陣不能相減，請重試elseprintf("請輸入矩陣A:n");for(p=0;p<i;p+)for(q=0;q<j;q+)scanf("%f',&Apq);printf("輸出矩陣A:n");for(p=0;p<i;p+)for(q=0;q<j;q+)(printf("%10.2f,Apq);if(q+1)%j=0)printf("n");printf("請輸入矩陣B

7、：n");for(p=0;p<i;p+)for(q=0;q<j;q+)scanf("%f',&Bpq);printf("輸出第矩陣B:n");for(p=0;p<i;p+)for(q=0;q<j;q+)(printf("%10.2f,Bpq);if(q+1)%j=0)printf("n");>printf("矩陣A-矩陣B為：n");計算兩個矩陣相減for(p=0;p<i;p+)for(q=0;q<j;q+)Cpq=Apq-Bpq;for(p=0;

8、p<i;p+)for(q=0;q<j;q+)<printf("%10.2f',Cpq);if(q+1)%j=0)printf("n");break;case 3:floatk;printf("請輸入矩陣A的行數(shù)和列數(shù)(用逗號隔開):")；scanf("%d,%d",&i,&j);printf("請輸入矩陣An");for(p=0;p<i;p+)for(q=0;q<j;q+)scanf("%f',&Apq);printf(&qu

9、ot;輸出矩陣An");for(p=0;p<i;p+)for(q=0;qvj;q+)(printf("%10.2f',Apq);if(q+1)%j=O)printf("n");)printf("請輸入一個實數(shù):n”)；seanf("%f',&k);for(p=0;p<i;p+)/數(shù)乘矩陣for(q=0;q<j;q+)Bpq=k*Apq;printf("輸出k乘矩陣A的結(jié)果山”)；for(p=0;p<i;p+)for(q=0;q<j;q+)(printf("%10

10、.2f',Bpq);if(q+1)%j=O)printf("n"););break;case 4: printfd請輸入矩陣A的行數(shù)和列數(shù)(用逗號隔開)seanf(“d,%d”，&printf("請輸入矩陣A:n");for(p=0;p<i;p+)for(q=0;q<j;q+)seanf("%f',&Apq);printf(輸出矩陣AnM);for(p=0;p<i;p+)for(q=0;q<j;q+)printf(H%10.2f',Apq);if(q+1)%j=O)printf(H

11、nM);)for(p=0;p<i;p+)/轉(zhuǎn)置for(q=0;q<j;q+)Bqp=Apq;printf("輸出矩陣A的轉(zhuǎn)置矩陣:n");for(p=0;p<j;p+)for(q=0;q<i;q+)(printfr%10.2r,Bpq);if(q+1)%i=0)printfCAn'1);/;break;case 5:floataM2*M;floatbN2wN;floatt,x;intk,T;printf("輸入方陣的維數(shù)：n");請輸入方陣，即行和列相等的矩陣。scanf(H%dH,&T);printf("

12、;請輸入矩陣:n");for(i=0;i<T;i+)for(j=O;j<T;j+)scanf("%f',&biO)；printf(H原矩陣為：nH);for(i=0;i<T;i+)(for(j=O;j<T;j+)printf(M%1O.3fbi0)；printfCAn'1);/for(i=0;i<T;i+)for(j=0;j<(2*T);j+)ifG<T)ai0=bi0；elseif(j=T+i)ai0=1.O;elseai0=O.O;for(i=0;i<T;i+)(for(k=0;k<T;k+)(if(k!=i)t=aki/aii;for(j=0;jv(2T);j+)x=aij*t;akO=akj-x;for(i=0;i<T;i+)4=aii;ford=0;j<(2*T);j+)aiO=aiO/t;for(i=0;i<T;i+)y=y)皿i；if(y=O)printf("對不起，您輸入的矩陣沒有逆矩陣，請重新輸入。n")；else(for(i=0;i<T;i+)forQ=0;j<T;j+)bij=aiU+T;printf("逆矩陣為:n");f