拉普拉斯差分滤波器模板实现+MATLAB

clear all;
close all;
clc;

k = 5/1; 
sigma = 1;
G1=fspecial('gauss',[round(k*sigma), round(k*sigma)], sigma);
% 使用diff函数做离散差分，会导致降低唯独
Gdx= diff(G1, 1, 1);
Gdy= diff(G1, 2, 1);
Gdxx= diff(G1, 1, 2);
Gdyy= diff(G1, 2, 2);

% 使用gradient保证差分后维度不改变，将原矩阵先插值，再进行差分
[Gx,Gy] = gradient(G1);   
[Gxx,Gxy] = gradient(Gx);
[Gyx,Gyy] = gradient(Gy);