get_awgn_Pb.m

function [Pb, M]=get_awgn_Pb(N,var)
%函数用来计算BEC信道的Bhattacharyya参数
%输入为信道个数N,BEC擦除概率e
%输出为Bhattacharyya参数矩阵
load('phi_inv_in_out_chart.mat')
n=log2(N);
temp=2/var;
for i=1:n
    for ii=1:2^(i-1)
        temp1(2*ii)=2*temp(ii);
        temp1(2*ii-1)=f2(temp(ii), phi_inv_in, phi_inv_out, length(phi_inv_in), phi_inv_out(1));
    end
    temp=temp1;
end
M=temp;
Pb=1-qfunc(sqrt(M./2)*(-1));

end