%??1-D信號壓縮傳感的實現（正交匹配追蹤法Orthogonal?Matching?Pursuit）
%??測量數M>=K*log（N/K）K是稀疏度N信號長度可以近乎完全重構
clc;clear

%%??1.?時域測試信號生成
K=7;??????%??稀疏度（做FFT可以看出來）
N=256;????%??信號長度
M=64;?????%??測量數（M>=K*log（N/K）至少40但有出錯的概率）
f1=50;????%??信號頻率1
f2=100;???%??信號頻率2
f3=200;???%??信號頻率3
f4=400;???%??信號頻率4
fs=800;???%??采樣頻率
ts=1/fs;??%??采樣間隔
Ts=1:N;???%??采樣序列
x=0.3*cos（2*pi*f1*Ts*ts）+0.6*cos（2*pi*f2*Ts*ts）+0.1*cos（2*pi*f3*Ts*ts）+0.9*cos（2*pi*f4*Ts*ts）;??%??完整信號

%%??2.??時域信號壓縮傳感
Phi=randn（MN）;???????????????????????????????????%??測量矩陣（高斯分布白噪聲）
s=Phi*x‘;????????????????????????????????????????%??獲得線性測量?
%%??3.??正交匹配追蹤法重構信號（本質上是L_1范數最優化問題）
%%??%??算法迭代次數（m>=K）
hat_x=csomp（sPhiNK）;


%%??4.??恢復信號和原始信號對比
figure（1）;
hold?on;
plot（hat_x‘k.-‘）?????????????????????????????????%??重建信號
p

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2013-10-17?12:48??Wavelet_OMP\
?????文件????????1135??2013-05-18?11:34??Wavelet_OMP\CS.asv
?????文件????????1122??2013-05-18?11:48??Wavelet_OMP\CS.m
?????文件????????1374??2013-05-18?11:39??Wavelet_OMP\csomp.asv
?????文件????????1304??2013-05-18?11:48??Wavelet_OMP\csomp.m
?????文件????????1065??2009-10-14?18:27??Wavelet_OMP\DWT.m
?????文件???????66614??2009-06-08?18:00??Wavelet_OMP\lena256.bmp
?????文件????????2287??2013-05-18?12:53??Wavelet_OMP\main1.m
?????文件????????1329??2013-04-10?09:02??Wavelet_OMP\omp.m
?????文件????????5220??2010-04-04?21:50??Wavelet_OMP\SL0.m
?????文件????????1228??2013-05-18?12:29??Wavelet_OMP\test.asv
?????文件????????1227??2013-05-18?12:51??Wavelet_OMP\test.m
?????文件????????2497??2013-05-18?10:48??Wavelet_OMP\Wavelet_OMP.m
?????文件?????????912??2010-04-01?16:44??Wavelet_OMP\wavelet_OMP_colour.m