โปรแกรม C++ เพื่อทำการคูณเมทริกซ์

เมทริกซ์คืออาร์เรย์ของตัวเลขสี่เหลี่ยมที่จัดเรียงในรูปแบบของแถวและคอลัมน์

ตัวอย่างของเมทริกซ์มีดังนี้

เมทริกซ์ 3*2 มี 3 แถวและ 2 คอลัมน์ดังแสดงด้านล่าง -

8 1
4 9
5 6

โปรแกรมทำการคูณเมทริกซ์มีดังนี้

ตัวอย่าง

#include<iostream>
using namespace std;
int main() {
   int product[10][10], r1=3, c1=3, r2=3, c2=3, i, j, k;
   int a[3][3] = { {2, 4, 1} , {2, 3, 9} , {3, 1, 8} };
   int b[3][3] = { {1, 2, 3} , {3, 6, 1} , {2, 4, 7} };
   if (c1 != r2) {
      cout<<"Column of first matrix should be equal to row of second matrix";
   } else {
      cout<<"The first matrix is:"<<endl;
      for(i=0; i<r1; ++i) {
         for(j=0; j<c1; ++j)
         cout<<a[i][j]<<" ";
         cout<<endl;
      }
      cout<<endl;
      cout<<"The second matrix is:"<<endl;
      for(i=0; i<r2; ++i) {
         for(j=0; j<c2; ++j)
         cout<<b[i][j]<<" ";
         cout<<endl;
      }
      cout<<endl;
      for(i=0; i<r1; ++i)
      for(j=0; j<c2; ++j) {
         product[i][j] = 0;
      }
      for(i=0; i<r1; ++i)
      for(j=0; j<c2; ++j)
      for(k=0; k<c1; ++k) {
         product[i][j]+=a[i][k]*b[k][j];
      }
      cout<<"Product of the two matrices is:"<<endl;
      for(i=0; i<r1; ++i) {
         for(j=0; j<c2; ++j)
         cout<<product[i][j]<<" ";
         cout<<endl;
      }
   }
   return 0;
}

ผลลัพธ์

The first matrix is:
2 4 1
2 3 9
3 1 8
The second matrix is:
1 2 3
3 6 1
2 4 7
Product of the two matrices is:
16 32 17
29 58 72
22 44 66

ในโปรแกรมข้างต้น เมทริกซ์ a และ b ทั้งสองถูกกำหนดค่าเริ่มต้นดังนี้ −

int a[3][3] = { {2, 4, 1} , {2, 3, 9} , {3, 1, 8} };
int b[3][3] = { {1, 2, 3} , {3, 6, 1} , {2, 4, 7} };

หากจำนวนคอลัมน์ในเมทริกซ์แรกไม่เท่ากับจำนวนแถวในเมทริกซ์ที่สอง จะไม่สามารถทำการคูณได้ ในกรณีนี้ ข้อความแสดงข้อผิดพลาดจะพิมพ์ออกมา ได้ดังนี้

if (c1 != r2) {
   cout<<"Column of first matrix should be equal to row of second matrix";
}

ทั้งเมทริกซ์ a และ b จะแสดงโดยใช้ nested for loop สิ่งนี้แสดงให้เห็นโดยข้อมูลโค้ดต่อไปนี้

cout<<"The first matrix is:"<<endl;
for(i=0; i<r1; ++i) {
   for(j=0; j<c1; ++j)
   cout<<a[i][j]<<" ";
   cout<<endl;
}
cout<<endl;
cout<<"The second matrix is:"<<endl;
for(i=0; i<r2; ++i) {
   for(j=0; j<c2; ++j)
   cout<<b[i][j]<<" ";
   cout<<endl;
}
cout<<endl;

หลังจากนี้ เมทริกซ์ผลคูณ [][] จะถูกเตรียมข้อมูลเบื้องต้นเป็น 0 จากนั้นจึงใช้ลูปที่ซ้อนกันเพื่อค้นหาผลคูณของเมทริกซ์ 2 เมทริกซ์ a และ b ซึ่งแสดงให้เห็นในข้อมูลโค้ดด้านล่าง

for(i=0; i<r1; ++i)
for(j=0; j<c2; ++j) {
   product[i][j] = 0;
}
for(i=0; i<r1; ++i)
for(j=0; j<c2; ++j)
for(k=0; k<c1; ++k) {
   product[i][j]+=a[i][k]*b[k][j];
}

หลังจากได้รับผลิตภัณฑ์แล้วจะมีการพิมพ์ ดังแสดงไว้ดังนี้

cout<<Product of the two matrices is:"<<endl;
for(i=0; i<r1; ++i) {
   for(j=0; j<c2; ++j)
   cout<<product[i][j]<<" ";
   cout<<endl;
}