vec/2.cpp

#include <stdio.h>
#include <math.h>
#include <time.h>

void applySeparableGaussianBlur(float src[5][5], float dst[5][5], int h, int w, float kx[3], float ky[3]) {
    float buf[5][5] = {0};
    
	//开始进行可分离卷积 
    for (int i = 1; i < h - 1; i++) {
    	// 计算当前行的行内卷积 
        for (int j = 1; j < w - 1; j++) {
            buf[i][j] = src[i][j-1]*kx[0] + src[i][j]*kx[1] + src[i][j+1]*kx[2];
        }
    }
	
	// 进行行间的卷积得到最终像素值 
    for (int i = 1; i < h - 1; i++) {
        for (int j = 1; j < w - 1; j++) {
            dst[i][j] = buf[i-1][j]*ky[0] + buf[i][j]*ky[1] + buf[i+1][j]*ky[2];
        }
    }
}

int main() {
    float src[5][5] = {
        {0, 0, 0, 0, 0},
        {0, 1, 2, 3, 0},
        {0, 4, 5, 6, 0},
        {0, 7, 8, 9, 0},
        {0, 0, 0, 0, 0}
    };
    float dst[5][5] = {0};
    float kx[3] = {0.25, 0.5, 0.25};
    float ky[3] = {0.25, 0.5, 0.25};
	

    clock_t start = clock();
    applySeparableGaussianBlur(src, dst, 5, 5, kx, ky);
    clock_t end = clock();

    printf("矩阵结果：\n");
    for (int i = 0; i < 5; i++) {
        for (int j = 0; j < 5; j++) {
            printf("%.2f ", dst[i][j]);
        }
        printf("\n");
    }
    printf("步骤2运行时间：%f秒\n", (double)(end - start) / CLOCKS_PER_SEC);

    return 0;
}