tcturbomud.m

safsdrghfydjg kukugilihg b nmvhjtdyewrtwert

function [LDM]=TCTurboMUD(SP,SM,PS,NS,OP,PO,PI,r,S,A,psd,MPI,MPO)

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

for j=1:PO
    [LE]=SISOMUD(La1,A,S,r,psd);
    for k=1:K
        DilrLE(k,:)=Deinter(MPO,LE(k,:));
    end
    for k=1:K
        LaUsers(k,:,:)=Row2mtrx(DilrLE(k,:),w);
    end
    for k=1:K
        LM=TD2mtrx(LaUsers,k,(M/w),w);
        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;
        LaUsers(k,:,:)=LETC;
        LDM(k,:)=LD;
    end
    for k=1:K
        DilrLE(k,:)=Mtrx2row(TD2mtrx(LaUsers,k,(M/w),w));
    end
    for k=1:K
        La1(k,:)=Interleaver(MPO,DilrLE(k,:));
    end
end