【线性代数】矩阵相乘

【线性代数】矩阵相乘

 1 //矩阵相乘
 2 #include<iostream>
 3 using namespace std;
 4  
 5 
 6 int main()
 7 {
 8     int am,an;
 9     int bm,bn;
10     int a[100][100];//矩阵A
11     int b[100][100];//矩阵B
12     int c[100][100]={0};//矩阵C 结果矩阵
13     cout<<"请输入矩阵A的行数和列数";
14     cin>>am>>an;
15     cout<<"请输入矩阵B的行数和列数";
16     cin>>bm>>bn;
17     if(an!=bm)//若A矩阵的列数须与B矩阵的行数不相等则无法进行乘法运算
18         cout<<"该两矩阵无法进行乘法运算！";
19     else
20     {//A矩阵的列数须与B矩阵的行数相等才能进行乘法运算
21         cout<<"请输入矩阵A"<<endl;
22         for(int i=0;i<am;i++)
23             for(int j=0;j<an;j++)
24             {
25                 cin>>a[i][j];
26             }
27         cout<<"请输入矩阵B"<<endl;
28         for(int i=0;i<bm;i++)
29             for(int j=0;j<bn;j++)
30             {
31                 cin>>b[i][j];
32             }
33         //下面进行矩阵的乘法计算
34 //结果矩阵的行列数由矩阵AB决定 结果矩阵为C am*bn
35         for(int i=0;i<am;i++)
36             for(int j=0;j<bn;j++)
37             {
38                 for(int ii=0;ii<an;ii++)
39                 {
40                     c[i][j]+=a[i][ii]*b[ii][j];
41                 }
42             }
43         cout<<"矩阵A*矩阵B的结果为："<<endl;
44         for(int i=0;i<am;i++)
45         {
46             for(int j=0;j<bn;j++)
47             {
48                 cout<<c[i][j]<<"";
49             }
50             cout<<endl;
51         }
52     }
53     system("pause");
54     return 0;
55 }

转载于:https://www.cnblogs.com/cplushua/archive/2011/11/15/2249630.html