经典polar编码译码程序

经典polar编码译码程序

clc
clear

n = 12;
N = 2^n;

Channel = ‘BEC‘;

p = 0.5;

if strcmp（Channel‘BSC‘）
    I_w = 1 + p*log2（p） + （1-p）*log2（1-p）;
elseif strcmp（Channel‘BEC‘）
    I_w = 1 - p;
end
Z_w = 1- I_w ;

temp_Z = Z_w;
for i = 1 : n
    temp_Z = [2*temp_Z - temp_Z.^2; temp_Z.^2];
    temp_Z = reshape（temp_Z[12^i]）;
end
[Val ind] = sort（temp_Z‘descend‘）;
Mutual_information = 1 - Val; % for BEC
plot（indMutual_information ‘sk‘‘MarkerFaceColor‘‘k‘）
ylim（[0 1]）;
xlabel（‘Bit channel index‘）; ylabel（‘capacity‘）;
title（[‘Capacity of bit channels N = ‘num2str（N）]）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2013-06-21 23:53  Polar_HW（1）\
     文件        6527  2013-06-21 23:53  Polar_HW（1）\BEC（R_is_0.7）.fig
     文件        2267  2013-06-21 23:53  Polar_HW（1）\BEC.mat
     文件        7471  2013-06-21 23:53  Polar_HW（1）\BSC（R_is_0.2）.fig
     文件        2169  2013-06-21 23:53  Polar_HW（1）\BSC.mat
     文件         586  2013-06-21 23:53  Polar_HW（1）\Compare_z_parameter.m
     文件         368  2013-06-21 23:53  Polar_HW（1）\Gen_G.m
     文件        3053  2013-06-21 23:53  Polar_HW（1）\Main_Homework.asv
     文件        2940  2013-06-21 23:53  Polar_HW（1）\Main_Homework.m
     文件        9368  2013-06-21 23:53  Polar_HW（1）\Polarization（1024）.fig
     文件        2982  2013-06-21 23:53  Polar_HW（1）\Polarization（16）.fig
     文件       17582  2013-06-21 23:53  Polar_HW（1）\Polarization（2048）.fig
     文件        4990  2013-06-21 23:53  Polar_HW（1）\Polarization（256）.fig
     文件       23100  2013-06-21 23:53  Polar_HW（1）\Polarization（4096）.fig
     文件        3377  2013-06-21 23:53  Polar_HW（1）\Polarization（64）.fig
     文件        1113  2013-06-21 23:53  Polar_HW（1）\prob.m
     文件        1235  2013-06-21 23:53  Polar_HW（1）\W.m