clear;
[fnamepname]?=?uigetfile（‘*.jpg‘‘Sample?Dialog?Box（原始圖像）‘）;
fullname=strcat（pnamefname）;
yuanshi=imread（fullname）;?
%yuanshi=rgb2gray（yuanshi）;%轉化為灰度圖像
figure（1）;
%subplot（121）;
imshow（uint8（yuanshi））;
title（‘原始圖像‘）;
[heightwidth]=size（yuanshi）;%原始圖像大小
I=im2double（yuanshi）;%圖像存儲類型轉換
T=dctmtx（8）;%離散余弦變換矩陣
V=blkproc（I[88]‘P1*x*P2‘TT‘）;%對原圖像進行DCT變換
V1=V*255;%圖像存儲類型轉換
V0=round（V1/5）;%DCT系數量化
for?m=1:height/8;
????for?n=1:width/8;
????????x=（m-1）*8+1;
????????y=（n-1）*8+1;
????????M11（mn）=V0（xy）;%每個DCT塊第（1，1）個點組成新矩陣
????end
end
%一次量化直方圖
[mn]=size（M11）;
MAX=max（max（M11））;
SS=zeros（1MAX+50）;
for?k=0:MAX+49????
???p=length（find（M11==k））;
???SS（k+1）=p/（m*n）;
end
%figure;
%subplot（122）;
%imshow（uint8（M11））;
figure;
%subplot（121）;
bar（0:MAX+

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2012-11-09?23:15??新建文件夾\
?????文件????????3269??2012-05-28?10:01??新建文件夾\xinxinxinxin1.asv
?????文件????????3267??2012-05-28?10:49??新建文件夾\xinxinxinxin1.m