现在的位置: 首页 > 综合 > 正文

【线性代数】矩阵相乘

2012年11月20日 ⁄ 综合 ⁄ 共 820字 ⁄ 字号 评论关闭
【线性代数】矩阵相乘
 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 }


抱歉!评论已关闭.