Pn=ceil(sqrt(K)); %K表示字典的大小,K取1024;
%初始化一个扁矩阵,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字典;
希望对你有帮组!!!
Pn=ceil(sqrt(K)); %K表示字典的大小,K取1024;
%初始化一个扁矩阵,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字典;
希望对你有帮组!!!