function X_dct=mydct_good(x,N)
%計(jì)算dct變換
%N必須是2的整數(shù)次冪,表示零填充后的矩陣大小
[row,col]=size(x);
%如果設(shè)定的填充大小比原始圖像的尺寸小,那么設(shè)定的填充尺寸不合理
if(N<row)|(N<col)
    disp('setted range isnot reasonable');
end
%補(bǔ)零
x(row+1:N,:)=0;
x(1:row,col+1:N)=0;

%求T矩陣
Tpq=zeros(N,N);
Tpq(1,:)=1/sqrt(N);
p=1:N-1;
q=0:N-1;
Tpq(p+1,q+1)=sqrt(2/N)*cos(p'*(2*q+1)*pi/(2*N));

%由T*x*T'得到DCT變換
X_dct=Tpq*x*Tpq';