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）]）;

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????目錄???????????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