资源简介
本例包含OMP、SP、IRLS、CoSaMP等重构算法,比较性能优劣。
代码片段和文件信息
function [Sestd]=cosamp2(PhiuutrueK)
% Cosamp algorithm
% Input
% K : sparsity of Sest
% Phi : measurement matrix
% u: measured vector
% tol1 : tolerance for approximation between successive solutions.
% Output
% Sest: Solution found by the algorithm
% d : success index (d=1 is success d = 0 is no convergence)
%
% Algorithm as described in “CoSaMP: Iterative signal recovery from
% incomplete and inaccurate samples“ by Deanna Needell and Joel Tropp.
%
% This implementation was written by David Mary
%
% This script/program is released under the Commons Creative Licence
% with Attribution Non-commercial Share Alike (by-nc-sa)
% http://creativecommons.org/licenses/by-nc-sa/3.0/
% Short Disclaimer: this script is for educational purpose only.
% Longer Disclaimer see http://igorcarron.googlepages.com/disclaimer
% Initialization
Sest=zeros(size(utrue));%不明白不知道utrue是什么
v=u;
t=1; T2=[];
while t < 101
[kz]=sort(abs(Phi‘*v));k=flipud(k);z=flipud(z);%
Omega=z(1:2*K);%取了前2K个。
T=sort(union(OmegaT2));phit=Phi(:T);%T为支撑集合并后的集合,PHIT为观察矩阵phi只取下标为T的列后所组成的矩阵
%keyboard
b=abs(pinv(phit)*u);%通过最小二乘法进行信号估计,pinv表示伪逆
[k3z3]=sort((b));k3=flipud(k3);z3=flipud(z3);%把得到的信号的新的估计值进行从大到小的排序
Sest=zeros(size(utrue));
Sest(T(z3(1:K)))=abs(b(z3(1:K))); %信号估计后进行删减,只取了前K个,并取绝对值后赋给sestsest此时为信号估计
[k2z2]=sort(abs(Sest));k2=flipud(k2);z2=flipud(z2);
T2=z2(1:K);%T2为信号估计值的最大K个下标集合。
v=u-Phi*Sest;%残差进行更新
d=0;n2=norm(abs(Sest-utrue)‘inf‘);
if n2 <1e-3
d=1;t=1e10;
disp(‘CoSaMP: success‘);
end
t=t+1;
end
if d==0
disp(‘CoSaMP: failed‘)
end
属性 大小 日期 时间 名称
----------- --------- ---------- ----- ----
文件 1780 2009-12-17 21:01 OMP_SP\cosamp2.m
文件 1369 2008-04-07 06:59 OMP_SP\IRLSregcomp.m
文件 767 2008-04-07 06:59 OMP_SP\krored.m
文件 1105 2013-07-04 10:06 OMP_SP\OMP.m
文件 2142 2008-04-07 06:58 OMP_SP\OMPmod3.m
文件 1256 2009-12-22 11:27 OMP_SP\SPA.m
文件 4932 2013-07-04 09:30 OMP_SP\test.m
目录 0 2013-07-04 10:23 OMP_SP
----------- --------- ---------- ----- ----
13351 8
相关资源
-
SPWM变频调速系统仿真Simuli
nk模型 - 数字信号处理-基于计算机的方法课本
- spiht改进算法matlab
- hopfield神经网络解决TSP问题
- NSP algorithm 彭思龙教授提出的零空间追
- Compressed Sensing 压缩感知书籍
- pso with wind 40uint
- 斯坦福2011年CVPR发表文章
- CS Primary tutorial CS压缩传感的初级教学
- 51个城市的TSP问题
- 交流电机spwm调速
-
simuli
nk仿真spwm逆变器 - 邮政运输网络中的邮路规划和邮车调
- 三电平逆变器的spwm仿真
- SPAMS的原始代码
- Crack dsp builder 11.1
-
Matlab ba
sed DIC code Version 1 130814 - SPIN_Matlab
- MRAS方法无速度传感器控制异步电机模
- 四种压缩传感的重构算法
- 语音信号时频分析 (the basic analysis
- KKSVD
- 1.5维谱
- 机器学习、人工智能、数据挖掘中经
- 石墨烯的电导率计算 (Complex conduc
- 提取地震子波的完整代码(Wave Extra
- SOMP算法代码
- LTE energy and spectrum efficiency
- 信号的功率谱香农熵和功率谱指数熵
- 基于字典学习的语音增强中字典更新
评论
共有 条评论