symbol_decision.m

B3g_phase2_C语言_Matlab程序及说明

function  [x_symbol,x_variance]=Symbol_Decision(x_rake,rou,var)

N=length(x_rake);

% 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);

for n=1:N
    
    P =exp(-abs(x_rake(n)-rou*MapMatrix16QAM).^2/2/var(n))/sum(exp(-abs(x_rake(n)-rou*MapMatrix16QAM).^2/2/var(n)));
    x_symbol(n)=P*MapMatrix16QAM.';
    x_variance(n)=P*abs(MapMatrix16QAM.^2)'-abs(x_symbol(n))^2;
end