BPSK.zip_BPSK generation_Eb/No_NOISE_bpsk noise_ebno

function decoded_signal = demodulation(mode, rx_signal) lengthmode = length(rx_signal)/mode; length_rx_signal = length(rx_signal); decoded_signal = zeros(length_rx_signal,1); temp = zeros(mode,length_rx_signal); switch mode; %BPSK case 1; decoded_signal = (sign(real(rx_signal))+1)./2; %QPSK case 2; temp(1,:) = (sign(real(rx_signal))+1)./2; temp(2,:) = (sign(imag(rx_signal))+1)./2; decoded_signal = temp(:); % %16QAM case 4; % 0010 0110 1110 1010 % 0011 0111 1111 1011 % 0001 0101 1101 1001 % 0000 0100 1100 1000 xr = sqrt(10).*real(rx_signal)'; xi = sqrt(10).*imag(rx_signal)'; temp(1,:) = (0<=xr); temp(2,:) = (-2<=xr).*(xr<=2); temp(3,:) = (0<=xi); temp(4,:) = (-2<=xi).*(xi<=2); decoded_signal = temp(:); %64QAM case 6; for i=1:lengthmode temp_re = real(rx_signal(i,1)); temp_im = imag(rx_signal(i,1)); n = (i-1)*mode+1; if temp_re>0 temp(n,1) = 0; if temp_re<(4/sqrt(42)) temp(n+2,1) = 0; if temp_re>(2/sqrt(42)) temp(n+4,1) = 0; else temp(n+4,1) = 1; end else temp(n+2,1) = 1; if temp_re<(6/sqrt(42)) temp(n+4,1) = 0; else temp(n+4,1) = 1; end end else temp(n,1) = 1; if temp_re>(-4/sqrt(42)) temp(n+2,1) = 0; if temp_re<(-2/sqrt(42)) temp(n+4,1) = 0; else temp(n+4,1) = 1; end else temp(n+2,1) = 1; if temp_re>(-6/sqrt(42)) temp(n+4,1) = 0; else temp(n+4,1) = 1; end end end if (temp_im>0) temp(n+1,1) = 0; if temp_im<(4/sqrt(42)) temp(n+3,1) = 0; if temp_im>(2/sqrt(42)) temp(n+5,1) = 0; else temp(n+5,1) = 1; end else temp(n+3,1) = 1; if temp_im<(6/sqrt(42)) temp(n+5,1) = 0; else temp(n+5,1) = 1; end end else temp(n+1,1) = 1; if temp_im>(-4/sqrt(42)) temp(n+3,1) = 0; if temp_im<(-2/sqrt(42)) temp(n+5,1) = 0; else temp(n+5,1) = 1; end else temp(n+3,1) = 1; if temp_im>(-6/sqrt(42)) temp(n+5,1) = 0; else temp(n+5,1) = 1; end end end end decoded_signal = temp; end