%初始化壹個扁矩陣,bb為字典原子維數開平方,如妳要產生256為的,則設置bb=16;
DCT=zeros(bb,Pn);
for k=0:1:Pn-1,
V=cos([0:1:bb-1]'*k*pi/Pn);
if k>0, V=V-mean(V); end;
DCT(:,k+1)=V/norm(V); end;
G=kron(DCT,DCT);%G就是最終的DCT字典;
希望對妳有幫組!!!