%??1-D信號壓縮傳感的實現（正交匹配追蹤法Orthogonal?Matching?Pursuit）
%??測量數M>=K*log（N/K）K是稀疏度N信號長度可以近乎完全重構
%??編程時間：2008年11月18日


clc;clear

%%??1.?時域測試信號生成
K=8;??????%??稀疏度??
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*sin（2*pi*f1*Ts*ts）+0.6*sin（2*pi*f2*Ts*ts）+0.1*sin（2*pi*f3*Ts*ts）+0.9*sin（2*pi*f4*Ts*ts）;??%??完整信號

%%??2.??時域信號壓縮傳感
Phi=randn（MN）;???????????????????????????????????%??測量矩陣（高斯分布白噪聲）
s=Phi*x.‘;????????????????????????????????????????%??獲得線性測量?

%%??3.??正交匹配追蹤法重構信號（本質上是L_1范數最優化問題）
m=2*K;????????????????????????????????????????????%??算法迭代次數（m>=K）
Psi=fft（eye（NN））/sqrt（N）;????????????????????????%??傅里葉正變換矩陣
T=Ph

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????2422??2010-04-21?21:08??OMP\CS_OMP.m

?????文件???????1065??2009-10-14?18:27??OMP\Wavelet_OMP\DWT.m

?????文件??????66614??2009-06-08?18:00??OMP\Wavelet_OMP\lena256.bmp

?????文件???????2499??2009-10-14?18:13??OMP\Wavelet_OMP\Wavelet_OMP.m

?????目錄??????????0??2010-06-21?20:16??OMP\Wavelet_OMP

?????目錄??????????0??2010-06-21?20:16??OMP

-----------?---------??----------?-----??----

????????????????72600????????????????????6