仿真极化码在不同码长下的极化效果图，包括散点图和柱状图 matlab代码.rar

clc clear n = 8; 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,[1,2^i]); end [Val, ind] = sort(temp_Z,'descend'); Mutual_information = 1 - Val; % for BEC plot(ind,Mutual_information ,'sk','MarkerFaceColor','b') ylim([0, 1]); xlabel('Bit channel index'); ylabel('capacity'); title(['Capacity of bit channels, N = ',num2str(N)]); figure; hist(Mutual_information,10); xlabel('Bit channel index'); ylabel('capacity'); title(['Capacity of bit channels, N = ',num2str(N)]);