直接上程式碼：

#include<iostream>
#include<vector>
using namespace std;

int main()
{
    //定義被卷積的矩陣（其實是一個數組，陣列元素的個數8*8）
    int const map = 8;
    float  A[map*map] =
    {
        1, 2, 3, 4, 5, 6, 7, 8,
        1, 2, 3, 4, 5, 6, 7, 8,
        1, 2, 3, 4, 5, 6, 7, 8,
        1, 2, 3, 4, 5, 6, 7
 
, 8,
        1, 2, 3, 4, 5, 6, 7, 8,
        1, 2, 3, 4, 5, 6, 7, 8,
        1, 2, 3, 4, 5, 6, 7, 8,
        1, 2, 3, 4, 5, 6, 7, 8
    };

    //定義卷積核矩陣（其實也是一個數組，陣列元素的個數3*3）
    int const kernel = 3;
    float B[kernel*kernel] =
    {
      1, 1, 1,
      1, 1, 1,
      1, 1, 1
    };

    //計算卷積輸出矩陣的維數（其實是輸出陣列元素個數的開根號）
 

    int const outm = map - kernel + 1;   //被卷積矩陣的維數-卷積核的維數+1  即8-3+1=6

    //計算卷積過程中的被卷積矩陣的寬和高（就是把寬拉成和卷積核的高一樣，這樣才好對應相乘）
    int const convAw = kernel*kernel;//3*3=9
    int const convAh = map*map;//8*8=64

    float A_convert[convAh*convAw] = { 0 };//定義一個卷積過程中的矩陣（也就是被拉長過後的矩陣）

    for (int i = 0; i < outm; i++)
    {
        for
 
 (int j = 0; j < outm; j++)
        {
            int wh = i * outm * convAw + j * convAw; 

            int col1 = i * map + j;
            A_convert[wh] = A[col1];        //第一次迴圈時把A[0] 的值賦給 A_convert[0]
            A_convert[wh + 1] = A[col1 + 1];//第一次迴圈時把A[1] 的值賦給 A_convert[1]
            A_convert[wh + 2] = A[col1 + 2];//第一次迴圈時把A[2] 的值賦給 A_convert[2]

            int col2 = (i + 1) * map + j;
            A_convert[wh + 3] = A[col2];    //第一次迴圈時把A[8] 的值賦給 A_convert[3]
            A_convert[wh + 4] = A[col2 + 1];//第一次迴圈時把A[9] 的值賦給 A_convert[4]
            A_convert[wh + 5] = A[col2 + 2];//第一次迴圈時把A[10] 的值賦給 A_convert[5]

            int col3 = (i + 2) * map + j;
            A_convert[wh + 6] = A[col3];     //第一次迴圈時把A[16] 的值賦給 A_convert[6]
            A_convert[wh + 7] = A[col3 + 1]; //第一次迴圈時把A[17] 的值賦給 A_convert[7]
            A_convert[wh + 8] = A[col3 + 2]; //第一次迴圈時把A[18] 的值賦給 A_convert[8]
        }

    }
    vector<int> C;
    for (int i = 0; i < outm; i++)
    {

        for (int j = 0; j < outm; j++)
        {
            int a = 0;
            int wh = i * outm * convAw + j * convAw;
            for (int m =0; m < convAw; m++)
            {
                a += A_convert[wh + m] * B[m] ;

            }
            C.push_back(a); //在C中新增資料a

        }

    }


    //輸出被卷積矩陣
    cout << "被卷積矩陣 :" << endl;
    for (int i = 0; i < map; i++)
    {
        for (int j = 0; j < map; j++)
        {
            cout << A[i*map + j] << " ";
        }
        cout << endl;
    }
    cout << endl;


    //輸出卷積核矩陣
    cout << "卷積核矩陣:" << endl;
    for (int i = 0; i < kernel; i++)
    {
        for (int j = 0; j < kernel; j++)
        {
            cout << B[i*kernel + j] << " ";
        }
        cout << endl;
    }
    cout << endl;


    //輸出卷積後輸出矩陣
    cout << "卷積後輸出矩陣：" << endl;
    for (int i = 0; i < outm; i++)
    {
        for (int j = 0; j < outm; j++)
        {
            cout << C[i*outm + j] << " ";
        }
        cout << endl;
    }

    system("pause");
    return 0;
}