clcclearclose?all;
img?=?imread（‘sjtu.jpg‘）;%讀取原圖像
gray?=?rgb2gray（img）;
[heightwidth]?=?size（gray）;???%輸入圖像是height?x?width的且hight>nwidth>n
x1?=?double（gray）;??
x2?=?x1;
name={‘3x3?Laplace‘‘5x5?Laplace‘‘7x7?Laplace‘‘9x9?Laplace‘}
sort=1;
n1=[3?5?7?9]
for?n?=?n1;
a（1:n1:n）?=?1;???%n×n模板元素全是1
for?i?=?1:height-n+1??
????for?j?=?1:width-n+1??
????????c?=?x1（i:i+（n-1）j:j+（n-1））.*a;?%取出x1中從（ij）開始的n行n列元素與模板相乘??
????????s?=?sum（sum（c））;?????????????????%求c矩陣中各元素之和??
????????x2（i+（n-1）/2j+（n-1）/2）?=?s/（n*n）;?%將與模板運算后的各元素的均值賦給模板中心位置的元素??
????end?
end??
final?=?uint8（x2）;
subplot（22sort）;
%?Laplace
H?=?fspecial（‘unsharp‘）;
sharpened?=?imfilter（finalH‘replicate‘）;
imshow（sharpened）;
title（name{sort}）;
sort=sort+1;
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????827??2018-04-10?12:00??exp2_Ave.m
?????文件?????????749??2018-04-10?12:00??exp2_Gauss.m
?????文件?????????691??2018-04-10?12:00??exp2_Midfilter.m