%%%%%%%%%%%%%?若干頻率正弦疊加信號(hào)?隨機(jī)欠采樣?壓縮感知恢復(fù)（OMP算法）?%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%?國防科學(xué)技術(shù)大學(xué)?楊俊剛?2011.03?%%%%%%%%%%%%%%%%%%%%

clear;

Tend=6;?%?信號(hào)持續(xù)時(shí)間?0-Tend
fs=200;?%?原始波形采樣頻率
tt1=0:1/fs:Tend;?%?原始波形采樣時(shí)間點(diǎn)
N=size（tt12）;?%?采樣點(diǎn)數(shù)?
f1=13;?%?正弦信號(hào)頻率
f2=5;
f3=2;
y1=cos（2*pi*f1*tt1）+cos（2*pi*f2*tt1）+cos（2*pi*f3*tt1）;?%?波形

M=50;??%?隨機(jī)欠采樣信號(hào)的點(diǎn)數(shù)

%%%%%%%%?生成隨機(jī)欠采樣位置?%%%%%%%%%%%%%%%%
indexM=fix（rand（1M）*N）;?
indexM=sort（indexM）;
if?（indexM（1）==0）
????indexM（1）=1;
end
?for?kk=1:M-1
?????while（indexM（kk+1）<=indexM（kk））
?????????indexM（kk+1）=indexM（kk+1）+1;
?????end
?end
?%%%%%%%%?生成隨機(jī)欠采樣位置?end?%%%%%%%%%%%%%%%%
?
?y2=y1（indexM）;?%?生成欠采樣序列

DCT_Matrix=（dct（eye（N）））.‘;?%?生成稀疏變換矩陣?采用DCT矩陣

Sense=DCT_Matrix（indexM:）;?%?生成對(duì)應(yīng)的測量矩陣

%%%%%%%%?OMP?算法恢復(fù)
%%%%%%%%?OMP?start?%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Aug=[];
cor

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件????????1557??2020-09-29?05:26??CS_Examples\CS_OMP.m
?????文件????????1220??2020-09-29?05:26??CS_Examples\CS_SPGL1.m
?????文件????????3861??2020-09-29?05:26??CS_Examples\spgl1_1_7\ChangeLog
?????文件?????????697??2020-09-29?05:26??CS_Examples\spgl1_1_7\Contents.m
?????文件???????26436??2020-09-29?05:26??CS_Examples\spgl1_1_7\COPYING
?????文件?????????184??2020-09-29?05:26??CS_Examples\spgl1_1_7\NormGroupL2_dual.m
?????文件?????????176??2020-09-29?05:26??CS_Examples\spgl1_1_7\NormGroupL2_primal.m
?????文件?????????375??2020-09-29?05:26??CS_Examples\spgl1_1_7\NormGroupL2_project.m
?????文件?????????221??2020-09-29?05:26??CS_Examples\spgl1_1_7\NormL12_dual.m
?????文件?????????209??2020-09-29?05:26??CS_Examples\spgl1_1_7\NormL12_primal.m
?????文件?????????463??2020-09-29?05:26??CS_Examples\spgl1_1_7\NormL12_project.m
?????文件??????????63??2020-09-29?05:26??CS_Examples\spgl1_1_7\NormL1_dual.m
?????文件??????????63??2020-09-29?05:26??CS_Examples\spgl1_1_7\NormL1_primal.m
?????文件?????????227??2020-09-29?05:26??CS_Examples\spgl1_1_7\NormL1_project.m
?????文件????????2783??2020-09-29?05:26??CS_Examples\spgl1_1_7\README
?????文件???????16195??2020-09-29?05:26??CS_Examples\spgl1_1_7\spgdemo.m
?????文件???????30273??2020-09-29?05:26??CS_Examples\spgl1_1_7\spgl1.m
?????文件????????4725??2020-09-29?05:26??CS_Examples\spgl1_1_7\spgSetParms.m
?????文件????????1608??2020-09-29?05:26??CS_Examples\spgl1_1_7\spgsetup.m
?????文件????????1594??2020-09-29?05:26??CS_Examples\spgl1_1_7\spg_bp.m
?????文件????????1814??2020-09-29?05:26??CS_Examples\spgl1_1_7\spg_bpdn.m
?????文件????????2513??2020-09-29?05:26??CS_Examples\spgl1_1_7\spg_group.m
?????文件????????1626??2020-09-29?05:26??CS_Examples\spgl1_1_7\spg_lasso.m
?????文件????????2853??2020-09-29?05:26??CS_Examples\spgl1_1_7\spg_mmv.m
?????文件????????1561??2020-09-29?05:26??CS_Examples\spgl1_1_7\private\ensure.m
?????文件????????6398??2020-09-29?05:26??CS_Examples\spgl1_1_7\private\heap.c
?????文件????????3660??2020-09-29?05:26??CS_Examples\spgl1_1_7\private\heap.h
?????文件???????11849??2020-09-29?05:26??CS_Examples\spgl1_1_7\private\lsqr.m
?????文件????????2838??2020-09-29?05:26??CS_Examples\spgl1_1_7\private\oneProjector.m
?????文件????????5704??2020-09-29?05:26??CS_Examples\spgl1_1_7\private\oneProjectorCore.c
?????文件????????1485??2020-09-29?05:26??CS_Examples\spgl1_1_7\private\oneProjectorCore.h
............此處省略7個(gè)文件信息