stomp算法matlab实现

function [omp_result support residual_error] =StOMP（Sc dictionary_matrix S Pfa）
% Sc表观测
% dictionary_matrix表超完备矩阵
% iteration表迭代次数
% S稀疏度
residual_error=Sc;%残差
support = [];
omp_result=zeros（size（dictionary_matrix2）1）;
Tfa = sqrt（2*log10（1/Pfa））;
for k=1:10  %迭代10次                                
    support_0 = [];
    inner_product=dictionary_matrix‘*residual_error;
    nrm = norm（inner_product）;
%    nrm = norm（residual_error）;
    Threshold = Tfa*nrm/sqr