I=imread（‘lena.bmp‘）;?%讀取lena圖像
J=imnoise（I‘salt?&?pepper‘0.02）;?%在原圖中添加椒鹽噪聲
imwrite（J‘lena1.bmp‘）;?%將添加椒鹽噪聲的圖像寫入文件夾
[Amap]=imread（‘lena1.bmp‘）;?%讀取加入噪聲后的索引圖像
image=double（A）;?%使索引圖片像素的灰度值為雙精度浮點數(shù)

prompt=（‘請輸入窗口大小對應(yīng)n值?例3*3窗口，n=3‘）;%設(shè)置對話框中顯示的文字
name=‘輸入窗口大小‘;%設(shè)置對話框的名字
numlines=1;%設(shè)置對話框顯示文字的行數(shù)
defaultanswer={‘3‘};%設(shè)定編輯欄的默認值
answer=inputdlg（promptnamenumlinesdefaultanswer）;%將編輯欄的輸入值賦給cell型變量answer

n=eval（answer{1}）;%將cell型變量answer轉(zhuǎn)換為整型變量賦給n，作為窗口邊長
z=n*n;%取出窗口中像素個數(shù)賦給z
u=zeros（1z）;%定義一個內(nèi)容全零的1xz數(shù)組，用于放需要排列的數(shù)
m=（n-1）/2;%定義窗口半徑

for?x=m+1:255-m;
????for?y=m+1:255-m;%進入循環(huán)，