【文件属性】:
文件名称:压缩传感OMP程序,实现压缩数据重构
文件大小:2KB
文件格式:NONE
更新时间:2012-10-29 07:43:23
OMP
压缩传感中OMP算法
for times=1:m; % 迭代次数
for col=1:N; % 恢复矩阵的所有列向量
product(col)=abs(T(:,col)'*r_n); % 恢复矩阵的列向量和残差的投影系数(内积值)
end
[val,pos]=max(product); % 最大投影系数对应的位置
Aug_t=[Aug_t,T(:,pos)]; % 矩阵扩充
T(:,pos)=zeros(M,1); % 选中的列置零(实质上应该去掉,为了简单我把它置零)
aug_y=(Aug_t'*Aug_t)^(-1)*Aug_t'*s; % 最小二乘,使残差最小
r_n=s-Aug_t*aug_y; % 残差
pos_array(times)=pos; % 纪录最大投影系数的位置
网友评论
- 很基础易懂。
- 简单有效,用的很顺手,感谢楼主
- 指令、语句的注解很详细,楼主很用心帮我解决了一个问题,完成了一个任务。感谢楼主
- 能够简单的对OMP算法有基础的了解,例子还不错。
- 很基础易懂,学习了
- 基础的OMP原理性代码范例,值得学习