function d=highf(a,n); % 高提升濾波函數，a是圖像，n是乘法系數
w1=[1 1 1;1 -8 1;1 1 1];
g1=fil(a,w1); %由于老師不許使用Matlab中現有的函數，是以自己編寫fil代替濾波函數
g2=(n*a-g1); %由于-8<0,是以是減号
immax=max(max(g2)); %找到g2矩陣的最大值
immin=min(min(g2)); %找到g2矩陣的最小值
d=(g2-immin)./(immax-immin); %原圖檔與拉普拉斯濾波後的圖像相加減，
%結果可能超出[0，1]的範圍，是以進行歸一化将immin，immax經過線性變換y=b*x+c成0，1
end
結果顯示
