%?已知?A?和?y（=Ax）?求?x
%?聯系方式：sy1133@163.com
%===========================================


clear;?clc

n?=?500;?%?信號長度
m?=?100;?%?測量次數?
k?=?20;??%?稀疏度

%Generate?sparse?signal?
z?=?randperm（n）;
x?=?zeros（n?1）;
x（z（1:k））?=?sign（randn（k1））;

A?=?randn（mn）;?%?用于重構的高斯隨機矩陣
y?=?A*x;????????%?投影信息：y

alpha?=?1;?%?如果alpha取2就是CoSaMP算法

r?=?y;?L?=?[];?a_index=[];a_index2=[];
a?=?zeros（size（x））;
iter?=?1;
err?=?1e-5;?%?信號誤差

tic
while?（iter?err）
????
????%?Identify
????h?=?A‘*?r;
????[h_newh_index]?=?sort（abs（h）‘descend‘）;
????
????%?Merge
????L?=?union（a_index2h_index（1:alpha*k））;
????
????%?Estimate
????a（L）?=?A（:L）\y;??
????
????%Prune?
????[a_newa_index]?=?sort（abs（a）‘descend‘）;
????a（a_index（k+1:end））=0;?
????a_index2?=?a_index（1

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件??????299441??2010-11-08?13:37??CoSaMP.pdf
?????文件??????455354??2010-12-18?09:25??Somp.pdf
?????文件????????1097??2010-12-20?08:41??Somp.m