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