赞
踩
function [x1] = OMP(K,A,y) %UNTITLED 此处显示有关此函数的摘要 % 此处显示详细说明 %Initialize x=0; r=y;%残差向量 I=[];%索引集合 B=[];%储存更新之后的A列向量 k=0;%循环索引 while(k<K) in_product=A'*r; [val pos]=max(abs(in_product));%返回内积绝对值中,最大元素的值及其位置 I=[I;pos];%更新索引 B=[B A(:,pos)];%采样矩阵中的列 %size(B) x=pinv(B)*y;%最小二乘法求目标信号 r=y-B*x;%更新残差向量 k=k+1; end [M N]=size(A); x1=zeros(N,1); x1(I)=x; end
Copyright © 2003-2013 www.wpsshop.cn 版权所有,并保留所有权利。