拉普拉斯差分滤波器模板实现+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);