gusucode.com > Diabetic retinopathy > compute_cosine_modified_basis.m
function B = compute_cosine_modified_basis(M,N) B=zeros(M*N,M*N); for k=1:M*N for t=1:M*N p=floor((k-1)/N); q=floor((t-1)/N); i=mod((k-1),N); x=mod((t-1),N); if(p==q) B(k,t)=cos(i*pi*((2*x+1)/2)/N); else B(k,t)=0; end end end return