ppm_soft_decode.m

这是关于PPM信号解调方法的研究

function [rxbits llr]=PPM_soft_decode(r_PPMseq,numbits)
rxbits=zeros(1,numbits);
for nb=1:numbits/4
    %mxk=r_PPMseq(1+(nb-1)*160:160+(nb-1)*160);
    zp=r_PPMseq(1+(nb-1)*16:16+(nb-1)*16);
    %for np=1:16
        %mxkp=mxk(1+(np-1)*10:10+(np-1)*10);
        %zp(np)=sum(mxkp);
    %end
    zp_sum=sum(zp);
    p=[zp(1)/zp_sum zp(2)/zp_sum zp(3)/zp_sum zp(4)/zp_sum zp(5)/zp_sum zp(6)/zp_sum zp(7)/zp_sum zp(8)/zp_sum zp(9)/zp_sum zp(10)/zp_sum zp(11)/zp_sum zp(12)/zp_sum zp(13)/zp_sum zp(14)/zp_sum zp(15)/zp_sum zp(16)/zp_sum];
    %L1=log((p(9)+p(10)+p(11)+p(12)+p(13)+p(14)+p(15)+p(16))/(1-(p(9)+p(10)+p(11)+p(12)+p(13)+p(14)+p(15)+p(16))));
    %L1=(p(9)+p(10)+p(11)+p(12)+p(13)+p(14)+p(15)+p(16));%/(1-(p(9)+p(10)+p(11)+p(12)+p(13)+p(14)+p(15)+p(16)));
    L1 = log((1-(p(9)+p(10)+p(11)+p(12)+p(13)+p(14)+p(15)+p(16)))/(p(9)+p(10)+p(11)+p(12)+p(13)+p(14)+p(15)+p(16)));
    %L2=log((p(5)+p(6)+p(7)+p(8)+p(13)+p(14)+p(15)+p(16))/(1-(p(5)+p(6)+p(7)+p(8)+p(13)+p(14)+p(15)+p(16))));
    %L2=(p(5)+p(6)+p(7)+p(8)+p(13)+p(14)+p(15)+p(16));%/(1-(p(5)+p(6)+p(7)+p(8)+p(13)+p(14)+p(15)+p(16)));
    L2 = log((1-(p(5)+p(6)+p(7)+p(8)+p(13)+p(14)+p(15)+p(16)))/(p(5)+p(6)+p(7)+p(8)+p(13)+p(14)+p(15)+p(16)));
    %L3=log((p(3)+p(4)+p(7)+p(8)+p(11)+p(12)+p(15)+p(16))/(1-(p(3)+p(4)+p(7)+p(8)+p(11)+p(12)+p(15)+p(16))));
    %L3=(p(3)+p(4)+p(7)+p(8)+p(11)+p(12)+p(15)+p(16));%/(1-(p(3)+p(4)+p(7)+p(8)+p(11)+p(12)+p(15)+p(16)));
    L3 = ((1-(p(3)+p(4)+p(7)+p(8)+p(11)+p(12)+p(15)+p(16)))/(p(3)+p(4)+p(7)+p(8)+p(11)+p(12)+p(15)+p(16)));
    %L4=log((p(2)+p(4)+p(6)+p(8)+p(10)+p(12)+p(14)+p(16))/(1-(p(2)+p(4)+p(6)+p(8)+p(10)+p(12)+p(14)+p(16))));
    %L4=(p(2)+p(4)+p(6)+p(8)+p(10)+p(12)+p(14)+p(16));%/(1-(p(2)+p(4)+p(6)+p(8)+p(10)+p(12)+p(14)+p(16)));
    L4 = ((1-(p(2)+p(4)+p(6)+p(8)+p(10)+p(12)+p(14)+p(16)))/(p(2)+p(4)+p(6)+p(8)+p(10)+p(12)+p(14)+p(16)));
   % if L1>0.5
   %     a1=1; 
   % else 
   %     a1=0;
   % end
   % if L2>0.5
   %     a2=1; 
   % else 
   %     a2=0;
   % end
   % if L3>0.5
   %     a3=1; 
   % else 
   %     a3=0;
    %end
    %if L4>0.5 
    %    a4=1; 
   % else 
    %    a4=0;
   % end
    %rxbits(1+(nb-1)*4:4+(nb-1)*4) = [a1 a2 a3 a4]
    llr(1+(nb-1)*4:4+(nb-1)*4) = [L1 L2 L3 L4]
end