基于压缩感知的MIMO信道估计算法

基于压缩感知的MIMO信道估计算法，使用了CoSaMP算法，GBP算法，OMP算法，CE算法，还有ls信道估计算法，MMSE算法，MIMO信道模型，等很全面

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  回心转意的程序\