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

%1.產生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.產生測量矩陣，兩種方法：toeplitz矩陣與隨機高斯矩陣；
????%toeplitz矩陣
?%???X1=toeplitz_1（MN）;
?%???X2=toeplitz_1（MN）;
????
????%隨機高斯矩陣
????X1=randn（MN）;
????X2=randn（MN）;
????
%3.獲得信道沖擊響應觀測值，暫時未加入噪聲；
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.重構算法，預計使用傳統的信道估計方法（LSMMSE），及壓縮感知估計（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（‘發射天線1到接收天線1的信道沖激響應‘）;
subplot（422）;plot（h11）;title（‘omp重構的響應‘）;
subplot（423）;plot（h（2:））;title（‘發射天線1到接收天線2的信道沖激響應‘）;
subplot（424）;plot（h12）;title（‘omp重構的響應‘）;
subplot（425）;plot（h（3:））;title（‘發射天線2到接收天線1的信道沖激響應‘）;
subplot（426）;plot（h21）;title（‘omp重構的響應‘）;
subplot（427）;plot（h（4:））;title（‘發射天線2到接收天線2的信道沖激響應‘）;
subplot（428）;plot（h22）;title（‘omp重構的響應‘）;

figure;??%1發1收的原始與重構信號
hold?on;
plot（h11‘k.-‘）?????????????????????????????????%??重建信號
plot（h（1:）‘r*‘）???????????????????????????????%??原始信號
legend（‘Recovery‘‘Original‘）

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