function dct8x8(I,m,s) %定義DCT量化子程序

T=inline('dctmtx(8)');

y=blkproc(I,[8 8],T);%進行DCT變換，求DCT變換系數
y=round(y./m);% DCT系數除以量化表中的值，并且進行四舍五入取整數

p=1;te=1;  
while(p<=64)   
    for q=1:te
    y1(s+p)=y(te-q+1,q);p=p+1;
    end

    for q=te:-1:1
    y1(s+p)=y(te-q+1,q);p=p+1;
     end
end

f=haffman(y1);
c(s:s+64,1)=f(:,1);c(s:s+64,2)=f(:,2);c(s:s+64,3)=f(:,3)