算法仿真代码极化码.zip

tic clc clear all close all %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%全局定义%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% global N global Channel_Parameter global Generated_Matrix global glo_Det global saturate_LLR global saturate_PM global saturate_Mode %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 参数设置 %%%%%%%%%%%%%%%%%%%%%%%%%%%%% N = 1024; %设置码长 L = 16; %设置搜索路径 Block = 1; %设置仿真块数 EbN0 =[1 1.5 2 2.25]; %设置位信噪比范围 load('Channel_Parameter_1024','Channel_Parameter'); %设置对应的信道参数 load('Generated_Matrix_N=1024','Generated_Matrix'); %编码矩阵 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%定点参数设置%%%%%%%%%%%%%%%%%%%%%%%%%%%%% saturate_LLR = 8; %设置LLR字长 saturate_PM = 13; %设置PM字长 Fra_Shift =4 ; %设置小数点的移位运算系数 saturate_Mode = 'symmetric'; %设置定点模式 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% k = N/2; Info_Length = N-k-24; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% CRC-24校验 %%%%%%%%%%%%%%%%%%%%%%%%%% gen = crc.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]); glo_Det=crc.detector([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]); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Error_Rate = zeros(1,length(EbN0)); Block_Error = zeros(1,length(EbN0)); BLER = zeros(1,length(EbN0)); Info_index = find(Channel_Parameter==1); %%%%%%%%%%%%%%%%%%%%%%% 在编码矩阵中挑选出与信息位下标对应的行 %%%%%%%%%%%%%%% Info_Generated_Matrix = Generated_Matrix(Info_index,:); processing = 0; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% for jj = 1:length(EbN0) All_Num = 0; ii = 0; %%%%%%%%%%%%%%%%%%%%%%%% 设置需要仿真的帧数 %%%%%%%%%%%%%%%%%%%%%%%%%%%% while ((ii <=Block)) ii = ii + 1; Msg = round(rand(1,Info_Length)); %%%%%%%%%%%%%%%%%% 进行 polar 编码 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Info = (generate(gen,Msg'))'; Code = mod(Info*Info_Generated_Matrix,2); %%%%%%%%%%%%%%%%%%%%%基带处理加噪声%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Code = 1 - 2*Code; Delte_Square = (N/(2*(Info_Length)))/(10^((EbN0(jj)/10))); y = Code +sqrt(Delte_Square)*randn(1,N); LLR = (2*y)/Delte_Square ; %%%%%%%%%%%%%%%% LLR 定点%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% LLR=round(Fra_Shift*LLR'); %小数点移位 for i=1:length(LLR) LLR(i)=saturate(LLR(i),saturate_LLR,saturate_Mode); %浮点饱和 end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%译码开始%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Data_Output = SCL_CRC_Decoder_FP(LLR,L); %%%%%%%%%%%%%%%%%%%%误码统计%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Num = sum(bitxor(Msg,Data_Output)); if Num ~=0 Block_Error(jj) = Block_Error(jj) +1; end All_Num = Num + All_Num; processing = processing + 1; if mod(processing,100)==0 disp('目前正在进行中....'); processing/(100) disp('还要继续....') end end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% BLER(jj) = Block_Error(jj)/ii; Error_Rate(jj) = All_Num/(Info_Length*ii); end save('1024BLER_L=32(D0802)','BLER'); %每次仿真前记住更改保存目录，不然以前仿真的结果被覆盖了！！切记！！ save('1024Error_Rate_L=32(D0802)','Error_Rate'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%绘图%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% figure(1) semilogy(EbN0,Error_Rate,'-.rd'); title('码率0.5，码长256，SCL-CRC译码'); xlabel('EbN0'),ylabel('BER'); grid on; figure(2) semilogy(EbN0,BLER,'-.rd'); title('码率0.5，码长256，SCL-CRC译码'); xlabel('EbN0'),ylabel('BLER'); grid on; toc