clear?all;
I1?=?imread（‘try1.jpg‘）;
S?=?size（I1）;
M?=?S（1）;
N?=?S（2）;
%考慮M和N不能被8整除，需要先做判斷用0補足
complementM?=?M+mod（8-mod（M8）8）;
complementN?=?N+mod（8-mod（N8）8）;
I?=?zeros（complementMcomplementN）;
I（1:M1:N）=I1;
%分8*8小塊完成DCT變換
%直接截取前20（實際上是21）個數據（對角線掃描），其余變0
compressI?=?zeros（complementMcomplementN）;
tmp?=?zeros（88）;
for?i?=?1?:?（complementM/8）
????for?j?=?1?:?（complementN/8）
????????tmp?=?imDCT（?I（8*i-7:8*i??8*j-7:8*j）?）;
????????for?m?=?1:8
????????????for?n?=?1:8
????????????????if?（m+n）>7?
????????????????????tmp（mn）=0;
????????????????end
????????????end
????????end
????????compressI（8*i-7:8*i??8*j-7:8*j）?=?iimDCT（tmp）;
????end
end

%計算PSNR
tmp2?=?（double（I1）-compressI（1:M1:N））.^2;
PSNR?=?255^2/mean（tmp2（:））;
PSNRdB?=?10*log10（PSNR）;

subplot（211）;imshow（I1）;title（‘原圖‘）;
subplot（212）;imshow（uint8（compressI））;title（‘壓縮圖‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2015-04-23?15:08??實驗1（張鑫）\
?????文件?????????946??2015-04-01?13:04??實驗1（張鑫）\Homework1_DCT.m
?????文件?????????410??2015-03-30?00:01??實驗1（張鑫）\Homework1_DFT.m
?????文件?????????489??2015-04-15?18:09??實驗1（張鑫）\iimDCT.m
?????文件?????????555??2015-03-30?00:02??實驗1（張鑫）\iimFFT.m
?????文件?????????510??2015-04-15?18:08??實驗1（張鑫）\imDCT.m
?????文件?????????531??2015-03-29?23:27??實驗1（張鑫）\imFFT.m
?????文件????????8632??2015-03-29?08:32??實驗1（張鑫）\try1.jpg
?????文件????????7983??2015-03-29?08:32??實驗1（張鑫）\try2.jpg
?????文件?????????202??2015-04-22?19:07??實驗1（張鑫）\說明.txt