%??1-D信號壓縮傳感的實現(xiàn)（正交匹配追蹤法Orthogonal?Matching?Pursuit）
%??測量數(shù)M>=K*log（N/K）K是稀疏度N信號長度可以近乎完全重構(gòu)
%??編程人--香港大學(xué)電子工程系?沙威??Email:?wsha@eee.hku.hk
%??編程時間：2008年11月18日
%??文檔下載:?http://www.eee.hku.hk/~wsha/Freecode/freecode.htm?

clc;clear

%%??1.?時域測試信號生成
K=8;??????%??稀疏度??
N=256;????%??信號長度
M=64;?????%??測量數(shù)（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）;??%??完整信號
%S?=?zeros（500003）;
%[S（:1）?Fs?bits]?=?wavread（‘source4.wav‘）;
%[S（:2）?Fs?bits]?=?wavread（‘source2.wav‘）;
%[S（:3）?Fs?bits]?=?wavread（‘source3.wav‘）;
%S?=?S‘;
%for?i=10001:10256
%s（i-10000）=S（1i）;
%en

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

?????文件???????2787??2012-04-14?13:59??CS_OMP.m

?????文件???????1083??2012-04-10?22:41??OMPerr.m

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

?????????????????3870????????????????????2