二维卷积c代码

二维信号的卷积原理请参考另外一篇文章：http://blog.csdn.net/carson2005/article/details/43702241

这里直接给出参考代码：

void Conv2(int** filter, int** arr, int** res, int filterW, int filterH, int arrW, int arrH){int temp;for (int i=0; i<filterH+arrH-1; i++){for (int j=0; j<filterW+arrW-1; j++){temp = 0;for (int m=0; m<filterH; m++){for (int n=0; n<filterW; n++){if ((i-m)>=0 && (i-m)<arrH && (j-n)>=0 && (j-n)<arrW){temp += filter[m][n]*arr[i-m][j-n];}}}res[i][j] = temp;}}}void Conv2Test(){int filterW = 3;int filterH = 3;int arrW = 5;int arrH = 5;int resW = filterW + arrW - 1;int resH = filterH + arrH - 1;int** pFilter = new int*[filterH];for (int i=0; i<filterH; i++){pFilter[i] = new int[filterW];}int** arr = new int*[arrH];for (int i=0; i<arrH; i++){arr[i] = new int[arrW];}int** res = new int*[resH];for (int i=0; i<resH; i++){res[i] = new int[resW];}pFilter[0][0] = 1;pFilter[0][1] = 3;pFilter[0][2] = 1;pFilter[1][0] = 0;pFilter[1][1] = 5;pFilter[1][2] = 0;pFilter[2][0] = 2;pFilter[2][1] = 1;pFilter[2][2] = 2;arr[0][0] = 17;arr[0][1] = 24;arr[0][2] = 1;arr[0][3] = 8;arr[0][4] = 15;arr[1][0] = 23;arr[1][1] = 5;arr[1][2] = 7;arr[1][3] = 14;arr[1][4] = 16;arr[2][0] = 4;arr[2][1] = 6;arr[2][2] = 13;arr[2][3] = 20;arr[2][4] = 22;arr[3][0] = 10;arr[3][1] = 12;arr[3][2] = 19;arr[3][3] = 21;arr[3][4] = 3;arr[4][0] = 11;arr[4][1] = 18;arr[4][2] = 25;arr[4][3] = 2;arr[4][4] = 9;printf("pFilter: \n");for (int i=0; i<filterH; i++){for (int j=0; j<filterW; j++){printf("%d ", pFilter[i][j]);}printf("\n");}printf("*********************************************** \n");printf("arr: \n");for (int i=0; i<arrH; i++){for (int j=0; j<arrW; j++){printf("%d ", arr[i][j]);}printf("\n");}printf("*********************************************** \n");Conv2(pFilter, arr, res, filterW, filterH, arrW, arrH);printf("res: \n");for (int i=0; i<resH; i++){for (int j=0; j<resW; j++){printf("%d ", res[i][j]);}printf("\n");}printf("*********************************************** \n");}