lmmse.m

来自「ＭＩＭＯ－ＯＦＤＭ系统仿真」· M 代码 · 共 17 行

function X_Esti=Lmmse(Y_FFT,N_FFT,H_FFT,GAUSS_NOISE,SNR)

for j=1:2
   for i=1:N_FFT
           if GAUSS_NOISE==1
              X_Esti(:,(j-1)*N_FFT+i) =  inv( H_FFT(:,:,i)'*H_FFT(:,:,i)+2/10^(SNR/10)*eye(2) )*H_FFT(:,:,i)' * Y_FFT(:,(j-1)*N_FFT+i);
%               X_Esti(:,(j-1)*N_FFT+i) = H_FFT(:,:,i)' / ( H_FFT(:,:,i)'*H_FFT(:,:,i)+1/10^(SNR/10) ) * Y_FFT(:,(j-1)*N_FFT+i);             
%               X_Esti(:,(j-1)*N_FFT+i) = inv( H_FFT(:,:,i)) * Y_FFT(:,(j-1)*N_FFT+i);      %ZF
           else %无噪声
              % X_Esti(:,(j-1)*N_CBPS+i) = H_FFT(:,:,i)' * inv( H_FFT(:,:,i)'*H_FFT(:,:,i) ) * Data_FFT(:,(j-1)*N_CBPS+i);
              % X_Esti(:,(j-1)*N_CBPS+i) = H_FFT(:,:,i)' / (H_FFT(:,:,i)'*H_FFT(:,:,i)) * Data_FFT(:,(j-1)*N_CBPS+i);
              X_Esti(:,(j-1)*N_FFT+i) = inv( H_FFT(:,:,i)) * Y_FFT(:,(j-1)*N_FFT+i);      %ZF
           end
   end
end