N=10000; E=1; snr=10^(snr_in_dB/10); sgma=sqrt(E/(3*2*snr)); s000=[1 0]; s001=[cos(pi/4) sin(pi/4)]; s011=[0 1]; s010=[cos(3*pi/4) sin(3*pi/4)]; s110=[-1 0]; s111=[cos(5*pi/4) sin(5*pi/4)]; s101=[0 -1]; s100=[cos(7*pi/4) sin(7*pi/4)]; for i=1:N, temp=rand; if(temp