polar代码MATLAB (1).rar_Matlab代码_Polar码_polar_polar coding_sc译码

% debugging function PolarCode_Simulation %% 编码参数 clc;clear;clear all; L =16; R=1/2; K=L*R; n=log2(L); TM = TransformMatrixGeneration(n);%有比特翻转 SNR=[0:0.5:4]; % generator = [1 1 0 0 0 0 1 1 0 0 1 0 0 1 1 0 0 1 1 1 1 1 0 1 1]; % c = [1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]; % x^k generator = [1 1 0 0 1]; c = [1 0 0 0 0 ]; % x^k CRC_len=length(generator)-1; r_len=K-CRC_len; ChannelType='BAWGNC'; %% 仿真结果记录 bit_count_SC = zeros(size(SNR)); frame_count_SC= zeros(size(SNR)); bler_SC= zeros(size(SNR)); %% 最大蒙特卡洛次数 MaxCircle = 1000; for i=1:length(SNR) sigma = sqrt(0.5*10^(-0.1*SNR(i))/R); [ChannelSet, LLR]=guassian_approximation(L,K,sigma); G = TM(:,ChannelSet); % 生成矩阵 for cir = 1:MaxCircle x_m = randi([0, 1], r_len,1); multip = conv(c,x_m); [divid, remainder]=deconv(multip,generator); remainder=(mod(remainder(end-3:end),2))'; % 对除法的余数做模二运算 code=[x_m;remainder]; % 只是在原序列后面加上CRC码 x_s = BinaryMatrixProduct(G,code); x_a = 1-2*x_s; %% 信道 noise = sigma*randn(L,1); y= x_a+ noise; %% 解调 x_m_SC = SC_Decode(y, ChannelSet, sigma, ChannelType); % Path Metric x_m_sc=x_m_SC(1:r_len); [divid, remainder]=deconv(x_m_SC,generator); % 对多项式的除，就等于做反卷积 remainder=mod(remainder(end-3:end),2); % 对除法的余数做模二运算 % if isequal(remainder,[0;0;0;0]); % msgbox(['CRC校验正确！']); % else % msgbox(['CRC校验错误']); % end frame_count_SC(i) = frame_count_SC(i) + any(x_m_sc ~=x_m); if frame_count_SC(i) > 20 break; end end bler_SC(i) = frame_count_SC(i)/cir % ber_SC(i) = bit_count_SC(i)/cir /K; end semilogy(SNR, bler_SC, '-+'); hold on end