K=8;??????%?信道非零抽頭稀疏，即信道稀疏度；
N=128;????%?信道長度
M=64;?????%?輸出序列長度，即測量數(shù)（M>=K*log（N/K）至少40但有出錯(cuò)的概率）

%1.產(chǎn)生2*2稀疏多徑信道，共4條路徑，每條路徑稀疏度為K=8；
h=zeros（4N）;
p=zeros（1K）;
for?i=1:4
????????for?n=1:K
??????????p（1n）=randint（11[1?N]）;
??????????h（ip（1n））=2*rand（11）-1;
????????end
end

%2.產(chǎn)生測量矩陣，兩種方法：toeplitz矩陣與隨機(jī)高斯矩陣；
????%toeplitz矩陣
?%???X1=toeplitz_1（MN）;
?%???X2=toeplitz_1（MN）;
????
????%隨機(jī)高斯矩陣
????X1=randn（MN）;
????X2=randn（MN）;
????
%3.獲得信道沖擊響應(yīng)觀測值，暫時(shí)未加入噪聲；
X=zeros（64256）;
X（:1:128）=X1;
X（:129:256）=X2;
y1=X1*h（1:）‘+X2*h（3:）‘;
y2=X1*h（2:）‘+X2*h（4:）‘;

%4.重構(gòu)算法，預(yù)計(jì)使用傳統(tǒng)的信道估計(jì)方法（LSMMSE），及壓縮感知估計(jì)（BPOMPCoSaMP）作圖并對比
h1=omp（y1X2*N）;
h11=h1（11:128）;
h21=h1（1129:256）;
h2=omp（y2X2*N）;
h12=h2（11:128）;
h22=h2（1129:256）;

subplot（421）;plot（h（1:））;title（‘發(fā)射天線1到接收天線1的信道沖激響應(yīng)‘）;
subplot（422）;plot（h11）;title（‘omp重構(gòu)的響應(yīng)‘）;
subplot（423）;plot（h（2:））;title（‘發(fā)射天線1到接收天線2的信道沖激響應(yīng)‘）;
subplot（424）;plot（h12）;title（‘omp重構(gòu)的響應(yīng)‘）;
subplot（425）;plot（h（3:））;title（‘發(fā)射天線2到接收天線1的信道沖激響應(yīng)‘）;
subplot（426）;plot（h21）;title（‘omp重構(gòu)的響應(yīng)‘）;
subplot（427）;plot（h（4:））;title（‘發(fā)射天線2到接收天線2的信道沖激響應(yīng)‘）;
subplot（428）;plot（h22）;title（‘omp重構(gòu)的響應(yīng)‘）;

figure;??%1發(fā)1收的原始與重構(gòu)信號(hào)
hold?on;
plot（h11‘k.-‘）?????????????????????????????????%??重建信號(hào)
plot（h（1:）‘r*‘）???????????????????????????????%??原始信號(hào)
legend（‘Recovery‘‘Original‘）

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件????????1621??2015-05-07?14:18??回心轉(zhuǎn)意的程序\mimo.m
?????文件????????1376??2015-04-11?20:02??回心轉(zhuǎn)意的程序\omp.m
?????文件?????????173??2015-05-06?20:07??回心轉(zhuǎn)意的程序\toeplitz_1.m
?????目錄???????????0??2015-05-07?13:35??回心轉(zhuǎn)意的程序\