smmse_mse_calc.m

ls算法smmse算法仿真和对比

%Function Declaration:
function ms_error=SMMSE_MSE_calc(X,H,Y,Rgg,SNR_send);
%This function generates mean squared error for the the SMMSE estimator..
%EVALUATION OF Hsmmse
u=rand(64,64);
F=fft(u)*inv(u);%The 64 X 64 twiddle factor matrix..
I=eye(64,64);
Rhh=F*Rgg*F';
Hls =(inv(X)) * Y;
Hsmmse=Rhh*inv(Rhh+(1/SNR_send)*I)*Hls;
ms_error_mat=mean(((abs(H)-abs(Hsmmse))/abs(H)).^2);
for i=1:64
    if(ms_error_mat(i)~=0)
        ms_error=ms_error_mat(i);
    end
end