%下面涉及到的主要還是公式編程，具體的算法和你的資料一樣，我就不詳細介紹算法了，如有對函數不理解的再聯系我
Gray?=?imread（‘main.jpg‘）;
[MNO]?=?size（Gray）;
%下面2行用于對待大圖片，強制截取其一部分，減少運算量，由于是基于統計，所以對結果影響不大
%M?=?256;????????????%此2行與19、20行相關，詳見19、20行
%N?=?256;

%--------------------------------------------------------------------------
%1.將各顏色分量轉化為灰度，如果圖像是灰度圖像，就注釋掉此段
%??如果未注釋此行，運行會出現報錯“Index?exceeds?matrix?dimensions.”
%--------------------------------------------------------------------------
Gray?=?double（0.3*Gray（::1）+0.59*Gray（::2）+0.11*Gray（::3））;

%--------------------------------------------------------------------------
%2.為了減少計算量，對原始圖像灰度級壓縮，將Gray量化成16級
%--------------------------------------------------------------------------
for?i?=?1:M
????for?j?=?1:N
????????%for?n=1:floor（M+N/2）/16??%如果5、6行被注釋，則使用此行。floor函數為向負無窮大方向近似，即近似為等于或小于自己的整數
????????for?n?=?1:256/16????????%如果5、6行未注釋，則使用此行（其實5、6行注釋與否使用上行結果都一樣，因為M+N/2==256）
????????????if?（n-1）*16<=Gray（ij）&Gray（ij）<=（n-1）*16+15
????????????????Gray（ij）?=?n-1;
????????????end
????????end
????end
end

%--------------------------------------------------------------------------
%3.計算四個共生矩陣P取距離為1，角度分別為04590135
%--------------------------------------------------------------------------
P?=?zeros（16164）;
for?m?=?1:16
????for?n?=?1:16
????????for?i?=?1:M
????????????for?j?=?1:N
????????????????if?j????????????????????P（mn1）?=?P（mn1）+1;
????????????????????P（nm1）?=?P（mn1）;
????????????????end
????????????????if?i>1&j????????????????????P（mn2）?=?P（mn2）+1;
????????????????????P（nm2）?=?P（mn2）;
????????????????end
????????????????if?i????????????????????P（mn3）?=?P（mn3）+1;
????????????????????P（nm3）?=?P（mn3）;
????????????????end
????????????????if?i????????????????????P（mn4）?=?P（mn4）+1;
?????????