jointeqdec.m

safsdrghfydjg kukugilihg b nmvhjtdyewrtwert

function [LDM]=JointEqDec(r,MPI,MPO,SP,SM,PS,NS,OP,PO,PI,SIG0,SIGR,Lc)

w=3; % Number of code bits per data bit.
J=2;
s=size(r);
M=s(1);
La1=zeros(1,M);
LEI=zeros(M/w,J);

for j=1:PO
    
    [LE]=BCJRTurbo(r,La1,1,PS,NS,OP,SIG0,SIGR,Lc);
    DilrLE=Deinter(MPO,LE-Lc*r(:,1)'-La1);
    La=Row2mtrx(DilrLE,w);
        
    LM=La;
    for x=1:PI
        [LEO1,LD1]=MoMAPDecoder([LM(:,1) LM(:,2)],LEI(:,2),SP,SM,PS,NS,OP,MPI);
        [LEO2,LD2]=MoMAPDecoder([LM(:,1) LM(:,3)],LEI(:,1),SP,SM,PS,NS,OP,MPI);
        LEI=[LEO1(:,1) LEO2(:,1)];
    end
    LETC=[LEO1(:,1)+LEO2(:,1) LEO1(:,2) LEO2(:,2)];
    LD=(LD1+LD2)*0.5;
    La=LETC;
    LDM=LD;
        
    DilrLE=Mtrx2row(La);
    La1=Interleaver(MPO,DilrLE);
            
end