%%%%%%%%%%%%%初始化%%%%%%%%%%%%%%
N=40;
t=[1317192531];
f=[2515203010];
H=zeros（NN）;
for?k=1:5
????H（f（k）t（k））=1;
end
figure（1）;
pcolor（H）;
%%%%%%%%%%%%%%y&A%%%%%%%%%%%%%%%
fA=zeros（N1）;
for?k=1:N
????fA（k）=（1/N^0.5）*exp（2*pi*k^3*j/N）;
end
y=H*fA;
W=zeros（NN）;
for?p=0:N-1
????for?q=0:N-1
????????W（p+1q+1）=exp（2*pi*p*q*j/N）;
????end
end
A=zeros（NN*N）;
fa=fA;
for?k=0:N-1
????D=diag（fa）;
????A（1:N（k*N+1）:（k*N+N））=D*W;
????for?l=1:N
????????if?l==N
????????????fa（l）=fa（1）;
????????else
????????????fa（l）=fa（l+1）;
????????end
????end
end
%%%%%%%%%%%%l1范數求解%%%%%%%%%%%%%%
z=zeros（N*N1）;
s=zeros（N*N1）;
K=5;
u=zeros（K1）;
v=zeros（K1）;
for?k=1:10000
????a=s+A‘*（y-A*s）;
%%??a=abs（a）;
????[mn]=sort（a1）;
????u=m（N-4:N）‘;
????v=n（N-4:N）‘;
????for?i=1:5
????????v1=v（i）;
????????z（v1）=u（i）;
????end
????s=z;
????z=zeros（N*N1）;
end
%%%%%%%%%%%%仿真圖%%%%%%%%%%%%%%%%%
H1=zeros（NN）;
for?k=1:N
????H1（k:）=s（（k-1）*N+1:k*N）‘;
end
H1=abs（H1）;
figure（2）;
pcolor（H1）;

????
????????
????

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1069??2013-06-08?19:53??CSRsim.m