matlab 非锐化掩蔽、高提升滤波

clc;
clear;

pic = imread（‘../../pic/Fig0340（a）（dipxe_text）.tif‘）;
figure（1）;
imshow（pic[0255]）;
pic = double（pic）;
[mn]=size（pic）;

h = zeros（55）;
sum = 0;
%ρ = 3 %高斯滤波
for i = 1: 5
    for j = 1 : 5
        h（ij） = floor（exp（-（（i-3）^2 + （j-3）^2）/（2*3*3））*1024）;
        sum = sum + h（ij）;
    end
end


pic_exp = zeros（m+4n+4）;
pic_exp（3:end-23:end-2） = pic;
pic_exp（1:） = pic_exp（5:）;
pic_exp（2:） = pic_exp（4:）;
pic_exp（:1） = pic