文档介绍：.1压缩感知部分
压缩感知算法主要可分为三类：贪婪迭代算法、凸凸优化（或最优化逼近方法）和基于贝叶 斯框架提出的重构算法。由于第三类方法注重信号的时间相关性，不适合图像处理问题，故 目前的研究成果主要集中在前两类中。目前已实现6中算法，分匹配追踪(Stagewise OMP)也是由OMP改进而来的一种贪心算法，与CoSaMP、 SP算法类似，不同之处在于CoSaMP、SP算法在迭代过程中选择的是与信号内积最大的2K 或K个原子，而StOMP是通过门限阈值来确定原子。此算法的输入参数中没有信号稀疏度 K，因此相比于ROMP及CoSaMP有独到的优势。
StOMP的算法流程：
ere^ sere t sea
输入；
CDMxy的传感矩阵A = W
⑵MM观测向筮,
⑴送代族数4默认为ID
〔■on限参数默认为m
输出；
⑴信号■稀疏表示系数估计&
Q）Nxl推残差令=7-4息
以下流程中，弓表示残差，成示这代被数，。表示空集，劣表示每枕迭代找到曲索 引例序号），4表示嵌迭代的索引例序与集宜注意；设元素个裁为々，一般有££3, 因为每欢迭代找到的索引为一般并非只含一个列序号），/表示矩阵4的第/列，4表示 搜索引4选出的矩阵刀的列集舍，执为队Q的列向堂，符与U表示集合并运算，6・）表 示求向堂内积，皿知［・］表示事模值使对值:、
■阅始化与二入4=。4 = 0』=»
［亨计箕i<=ciag日「弓_顼即计蔑（弓\1冬JW月）』选择中中大于门跟75?的值/将 这些值对应思的列序与/构成集合4（列序号集合）；
⑶令4 = 4-2禹，河=42"徵& 4）：若4 = 4』（即无籍列被谕中） 则停止遥代进入第（乃步,
［4｝求y=-M的最小二乘解；0 = arcmin|,-耳切 11 ='4X「在？ F
⑴更新殃差号=尹-40 =7-_4,4二4「1#> S
〔珈。+1,如果L<s^回第⑵=0则停止迭代进 入葡7｝堤
〔7.•重构所得白在4处有非零项，其值分别为最后一校送代所爵a・
注1；得到&后，利用稀疏矩阵可得重构信号x=^f
注土从输入参数中可以知道琴;算法并不需要已知信号稀疏度（文献［2］ 二段中提到藉Q'lP算法对稀疏度K依籁性很大，只有正瞒估计了信号的稀疏度，才能精 确重构原始信导％本人对此双点保留意见）
注％在浏量矩阵为随机高斯矩阵的情况下，第⑵步中的门服Nep,其中与取值 范困一股为2冬£3,冬=时|/灰，弓表示上次遥代计尊出的残差，||也表示？范熟 注4；随机高所测童矩阵如文献［3］（2-7｝所示，这里注意方差一定要是
（标准差为1/皿7 ）,即Matab生成代码为randn（M-^钩政明，这是因为第心步的原子选 择是将传感矩阵列向童与残差的内萩向童元素与残差向童2范歌的某倍薮相比较，对于前面 的几篇介藉的其它重构尊法此处无所谓，为保持一致此后尽重全臻以r观血GL©钩以来 生成测筮矩阵.
经StOMP算法重构后的结果如下所示：
该算法的用时情况如下:
探查摘要
基于M「立「伊时间于1 9桐a-20TE?2;28;2西生成
蹈官称 调用透总时间目用时同*忍甘亘图
〔深色箓育=白用时同」
STOMP 1 ；
imsho't
2
s
s
nmvvpl nt
4
s
£
new口1 ol:》。日:；更rveiFiq unJNextPlot
4
§
s
dr
2
s
5
mDrph 逐
4
s
s
■
1
£
s
■
imdi玳己
3
s
5
.
in'tSize
2
SO s
s
■
imope 门
1
s
s
■
(SwOMP)
分段弱正交匹配追踪(Stagewise Weak OMP)可以说是StOMP的一种修改算法，它们的 唯一不同是选择原子时的门限设置，这可以降低对测量矩阵的要求。我们称这里的原子选择 方式为〃弱选择〃(WeakSelection)，StOMP的门限设置由残差决定，这对测量矩阵(原子 选择)提出了要求，而SWOMP的门限设置则对测量矩阵要求较低(原子选择相对简单、粗 糙)。
SWOMP的算法流程：
输入
Mx V的传感坦阵4 = @妙
^x1维湿测向童卜
迭代决敬5,默认为抻
输出，
〔1