www.gusucode.com > Diabetic retinopathy > cosine_dct_basis.m
function B = cosine_dct_basis(N) B=zeros(N,N); for k=1:N for t=1:N B(k,t)=cos(((2*(t-1)+1)*(k-1)*pi)/(2*N)); end end return