dctm_two.m

生成DCT矩阵,只要简单输入矩阵的阶数,就可以返回一个用语DCT变换的矩阵

%this program is for the call of the dct transform matrix of type -11
function dctm=dctm(n)
for p=0:(n-1)
    for q=0:(n-1)
    t(p+1,q+1)=1/n^(1/2);
    end
end
for p=1:(n-1)
    for q=0:(n-1)
        t(p+1,q+1)=(2/n)^(1/2)*cos(pi*(2*q+1)*p/(2*n));
    end
end
dctm=t