channelestimation_lmmse.m

关于MATLAB中实现802.16e的物理层仿真,希望对大家有帮助!

function Hlmmse=channelestimation_LMMSE(P1,Pilot,trms,t_max,snr,Hp_Ls)

%LMMSE信道估计算法准则
N=length(P1);
beta=17/9; % 仅仅对于16QAM
Rhh=zeros(N,N);

%  计算Rhh --频域信道矢量的自协方差矩阵
for k=1:N
    for l=1:N
        Rhh(k,l)=(1-exp((-1)*t_max*((1/trms)+j*2*pi*(k-l)/N)))./(trms*(1-exp((-1)*t_max/trms))*((1/trms)+j*2*pi*(k-l)/N));
    end
end

Hlmmse=(Rhh*inv(Rhh+(beta/snr)*eye(N))*Hp_Ls')';