联合稀疏SOMP代码
来源:互联网 发布:淘宝账号和密码全忘记 编辑:程序博客网 时间:2024/05/28 23:10
原理网上有好多,在此就不写啦,这个是我自己在OMP基础上稍微改进的,作为一个函数可以调用,大家转载时要注明出处哟~~
function A=SOMP1(D,X,L)
% 输入参数:
% D - 过完备字典,注意:必须字典的各列必须经过了规范化
% X - 信号
% L - 系数中非零元个数的最大值(可选,默认为D的列数,速度可能慢)
% 输出参数:
% A - 稀疏系数
%如果函数的输入参数只有两个时即D和X,则L默认为D的列数
P=size(X,2);%返回X的列数,如果是size(X,1)返回的是X的行数
K=size(D,2);
temp=zeros(K,P);
a=zeros(K,P);
x=X;
residual=x; %残差
indx=zeros(L,1);%索引集
for j=1:1:L
proj=sum(residual'*D);
[maxVal,pos]=max(abs(proj));
pos=pos(1);
indx(j)=pos;
a=inv(D(:,indx(1:j))'*D(:,indx(1:j)))*D(:,indx(1:j))'*X;
residual=x-D(:,indx(1:j))*a;
temp(indx(1:j),:)=a;
A=sparse(temp);
if sum(sum(residual.^2)) < 1e-6
break;
end
end
return;
- 联合稀疏SOMP代码
- 稀疏表示之OMP,SOMP算法及openCV实现
- 联合稀疏-阅读笔记
- 联合查询代码示例
- 稀疏
- 两表联合查询插入代码:
- vcs联合编译v/sv/c++代码
- UNION联合多表查询代码举例
- ThinkPHP多表联合查询实例代码
- 深度学习笔记:稀疏自编码器(4)——稀疏自编码器代码练习
- 联合
- 联合
- 联合
- 联合
- 联合
- 联合
- 联合
- 稀疏矩阵【压缩算法】,具体代码没有,只是思想
- MySQL递归查询
- 面向对象——类设计(三)
- 遍历聚合对象中的元素——迭代器模式
- 杭电1180——诡异的楼梯(BFS+优先队列)
- PAT (Top Level) Practise 1005Programming Pattern (35)
- 联合稀疏SOMP代码
- php使用explode分割字符串新手容易忽略的问题
- 通过汇编一个简单的C程序,分析汇编代码理解计算机工作原理
- Log4J日志配置详解
- 系统框架 springmvc mybatis Bootstrap html5 shiro maven
- WebStorm快捷键
- PAT (Top Level) Practise 1006 Tree Traversals - Hard Version (35)
- 单例模式的三种常用的形式
- Android: 解决动画完成后位置恢复到初始位置的问题