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

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

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

%--------------------------------------------------------------------------
%3.計(jì)算四個(gè)共生矩陣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;
?????????