PolarCodesMatLab/lr_new.m at master · zking77/PolarCodesMatLab · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
function L=lr_new(y,u_hat,N,i,var)
%函数用来计算最大似然比,函数使用的递归的思路
%参数y,为译码器输入,u_hat为前面的译码结果,N为码长,i为当前要计算的译码输出编号

if N~=1                         %没有满足递归终止条件,递归继续
    if mod(i,2)==1              %i为奇数时
        if i==1
            L=(lr_new(y(1:N/2),u_hat,N/2,(i+1)/2,var)*lr_new(y(N/2+1:N),u_hat,N/2,(i+1)/2,var)+1)/(lr_new(y(1:N/2),u_hat,N/2,(i+1)/2,var)+lr_new(y(N/2+1:N),u_hat,N/2,(i+1)/2,var));    %i=1时
            if L>100
                L=100;
            else if L<1/100
                    L=1/100;
                end
            end
        else
            uoe=mod(u_hat(1:2:i-1)+u_hat(2:2:i-1),2);                                                                                                               %i~=1时
            ue=u_hat(2:2:i-1);
            L=(lr_new(y(1:N/2),uoe,N/2,(i+1)/2,var)*lr_new(y(N/2+1:N),ue,N/2,(i+1)/2,var)+1)/(lr_new(y(1:N/2),uoe,N/2,(i+1)/2,var)+lr_new(y(N/2+1:N),ue,N/2,(i+1)/2,var));
            if L>100
                L=100;
            else
                if L<1/100
                    L=1/100;
                end
            end
        end
    else                        %i为偶数时
        uoe=mod(u_hat(1:2:i-2)+u_hat(2:2:i-2),2);
        ue=u_hat(2:2:i-2);
        L=lr_new(y(1:N/2),uoe,N/2,i/2,var)^(1-2*u_hat(i-1))*lr_new(y(N/2+1:N),ue,N/2,i/2,var);
        if L>100
            L=100;
        else
            if L<1/100
                L=1/100;
            end
        end
    end
else
    L=exp(-2*y/var);

end


end