mean_var.m

B3g_phase2_C语言_Matlab程序及说明

function [x_rec, x_variance]=Mean_Var(LLR_P);

[M,N]=size(LLR_P);

% Gray coded 16-QAM mapping vector
MapMatrix16QAM = [-1-j -1-3*j -1+j -1+3*j -3-j -3-3*j -3+j -3+3*j 1-j 1-3*j 1+j 1+3*j 3-j 3-3*j 3+j 3+3*j].';
MapMatrix16QAM = MapMatrix16QAM/sqrt(10);
matTxBits=de2bi(0:15);

vecOrder = (2.^(3:-1:0)).';
vecTxBits = matTxBits * vecOrder + 1;
Symbols = MapMatrix16QAM(vecTxBits);

All_d=2*matTxBits'-1;
%All_d=2*All_d-1;
for n=1:N
    
    P=prod((1+diag(tanh(LLR_P(:,n)/2))*All_d)/2);
    
    x_rec(n)=P*Symbols;
    x_var(n)=P*(abs(Symbols).^2)-abs(x_rec(n))^2;
    
end
x_variance=x_var;