极化码仿真代码

极化码的编译码， matlab开发工具进行仿真验证，能快速上手了解

function initPC（NKEcN0designSNRdBsilentflag） %designSNRdB optional
%
% This prepares the collection of all implicit parameters related to
% polar coding & SC decoding; to be used by all subsequent routines later.
%   （Including the memory resources to be used by the polar SC decoding）
%
%        Usage: initPC（NKEcN0designSNRdB）
% 
%            N  -  Blocklength‘）;
%            K  -  Message length （Rate = K/N）‘）;
%            Ec -  the BPSK symbol energy （linear scale）;
%            N0 （optional） -  Noise power spectral density （default N0/2 = 1） 
%                               （used for Monte-Carlo simulate etc）‘）;
%            designSNRdB （optional） - the SNR at which the polar code construction （PCC） should be performed‘）;
%                                      （As far as PCC is concerned SNR=Ec/N0 by definition & defaults to “0dB“）‘）;
% 
%            silentflag （optional） - whether to print the last result or not
%                                       （defaults to 0）
% 
%        Note: This routine must be called before we use any other utility around here.\n\n‘）;
% 

if nargin==3 %Normalize noise-power-spectral-density when not supplied: --- N0/2=1 （default）
    N0=2;
    designSNRdB=0;
    silentflag=0;
elseif nargin==4 
    designSNRdB=0;
    silentflag=0;
elseif nargin==5
    silentflag=0;
elseif （nargin<3 || nargin>6）
    fprintf（‘\n   Usage: initPC（NKEcN0designSNRdB）\n‘）;
    fprintf（‘\n       N  -  Blocklength‘）;
    fprintf（‘\n       K  -  Message length （Rate = K/N）‘）;
    fprintf（‘\n       Ec -  the BPSK symbol energy‘）;
    fprintf（‘\n       N0 （optional） -  Noise power spectral density （default N0/2 = 1）‘）;
    fprintf（‘\n                           （used for Monte-Carlo simulate etc）‘）;
    fprintf（‘\n       designSNRdB （optional） - the SNR at which the polar code should be constructed‘）;
    fprintf（‘\n                                 （Here SNR=Ec/N0 - by definition for a PCC defaults to “0dB“）‘）;
    fprintf（‘\n\n   Note: This routine must be called before we use any other utility around here.\n\n‘）;
    return;
end

addpath（genpath（‘./functions‘））; %All helping routines.

global PCparams;

PCparams = struct（‘N‘ N ...
                  ‘K‘ K ...
                  ‘n‘ 0 ...
                  ‘FZlookup‘ zeros（N1） ...
                  ‘Ec‘ Ec ...
                  ‘N0‘ N0 ...        %N0/2 =1 normalization
                  ‘LLR‘ zeros（12*N-1） ...
                  ‘BITS‘ zeros（2N-1） ...
                  ‘designSNRdB‘ designSNRdB）;

PCparams.n = log2（N）;
PCparams.FZlookup = pcc（NKdesignSNRdB）; %polar code construction at the given design-SNR （default 0dB）

if（~silentflag）
fprintf（‘\n All polar coding parameters/resources initialized. （in a structure - “PCparams“） \n‘）;
disp（PCparams）;
end

end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2015-10-28 01:04  polarcoding_matlab\
     文件        2086  2015-10-28 01:04  polarcoding_matlab\Sample session.txt
     文件       25634  2015-10-28 01:02  polarcoding_matlab\diary
     文件        1500  2015-10-27 13:01  polarcoding_matlab\systematic_pencode.m
     文件        3776  2015-10-28 01:02  polarcoding_matlab\plotPC.m
     文件        3776  2015-10-28 01:02  polarcoding_matlab\plotPC_systematic.m
     文件        2249  2015-10-22 04:05  polarcoding_matlab\systematic_pdecode.m
     文件        3601  2015-10-13 08:19  polarcoding_matlab\pdecode.m
     文件        1025  2015-01-28 01:21  polarcoding_matlab\pencode.m
     文件        2801  2015-10-28 00:15  polarcoding_matlab\initPC.m
     目录           0  2015-10-27 12:56  polarcoding_matlab\functions\
     文件         926  2015-10-27 12:58  polarcoding_matlab\functions\pcc.m
     文件        1281  2015-10-27 12:52  polarcoding_matlab\functions\test_systematic.m
     文件         393  2015-10-22 03:09  polarcoding_matlab\functions\spencode_BPSK_AWGN.m
     文件         381  2015-10-22 03:04  polarcoding_matlab\functions\pencode_BPSK_AWGN.m
     文件         413  2015-10-24 00:53  polarcoding_matlab\functions\bitreversed.m
     文件         363  2015-01-28 01:21  polarcoding_matlab\functions\setdesignSNRdB.m
     文件         102  2015-01-28 01:21  polarcoding_matlab\functions\logdomain_diff.m
     文件         106  2015-01-28 01:20  polarcoding_matlab\functions\logdomain_sum.m
     文件         576  2015-10-22 01:19  polarcoding_matlab\functions\FN_transform.m
     文件         177  2015-01-28 01:20  polarcoding_matlab\functions\upperconv.m
     文件        1198  2015-01-29 07:08  polarcoding_matlab\functions\updateBITS.m
     文件         304  2015-01-28 01:20  polarcoding_matlab\functions\lowerconv.m
     文件         331  2015-10-22 03:02  polarcoding_matlab\functions\BPSK_AWGN.m
     文件        1298  2015-01-29 06:17  polarcoding_matlab\functions\updateLLR.m
     文件         518  2015-10-13 11:30  polarcoding_matlab\functions\EncoderC.m
     文件         646  2015-10-13 11:14  polarcoding_matlab\functions\EncoderB.m
     文件         719  2015-10-13 11:54  polarcoding_matlab\functions\EncoderA.m
     文件         248  2015-10-22 02:57  polarcoding_matlab\functions\BPSK.m