?m=30:50

L=100;%信號的長度

signal=zeros（L1）;
p=L+m-1;
T=5;
q=randperm（L）;
signal（q（1:T））?=?sign（randn（T1））;%產生稀疏信號
for?u=1:21
?X=zeros（p（u）+L-1L）;
b=rand（1p（u））;
for?i=1:p（u）
????if?b（i）>=0.5?b（i）=1;end
????if?b（i）<0.5?b（i）=-1;end
end
for?i=1:p（u）
X（i1）=b（1i）;
end
for?i=1:p（u）+L-2
????for?j=1:L-1
????????X（i+1j+1）=X（ij）;
????end
end




Xp=zeros（m（u）L）;
for?j=1:m（u）
???Xp（j:）=X（L-1+j:）;
end?
%提取其中m行
z=randn（1m（u））;
v=Xp*signal;
d=1;


yp=v+z‘;


signal_ls=Xp‘*inv（Xp*Xp‘）*yp;
delt1（u）=10*log10（norm（signal_ls-signal）/norm（signal））;
mm=2*T;????????????????????????????????????????????%??算法迭代次數（m>=K）
hat_y=zeros（1L）;?????????????????????????????????%??待重構的譜域（變換域）向量?????????????????????
Aug_t=[];?????????????????????????????????????????%??增量

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1879??2010-12-04?13:53??daopin.m

-----------?---------??----------?-----??----

?????????????????1879????????????????????1