map feed back

WAclear all disp('This program only consider the puncturing patterns according'); disp(' to CDMA-2000 standard, rc=1/4,1/3,1/2 cases.'); L_info = input(' Please enter the infomation size (frame= info + tail(2*m bits), default: 378)\n '); if isempty(L_info) L_info = 378; % L_frame=infomation bits plus tail bits end % decode channel condition disp(' Please choose channel condition. in AWGN channel, press 0'); disp(' uncorelatived Rayleigh fading channel press 1'); disp(' corelatived Rayleigh fading or Jack fading channel, press 2'); channel = input('default=0,AWGN channel '); if isempty(channel) channel=0; %AWGN case end % encoder generator %input form is "[gd; g1; g2;...gm]". gd is feed back generator, g1...gm are forward generator. disp(' Please enter code generator: ( default:'); g = input('g =[1 0 1 1; 1 1 0 1; 1 1 1 1]\n' ); if isempty(g) g = [1 0 1 1; 1 1 0 1; 1 1 1 1]; end %other generator patterns %[1 1 1;1 0 1] %[1 1 0 1;1 1 1 1] %[1 1 1 1 1;1 0 0 0 1] [n,K] = size(g); m = K - 1; %shift register memory max_states = 2^m; %trellis states L_frame=L_info+2*m; % z RSC encoder need z*m bits tail bit to flash the memory. %input the p=1/code rate p = input(' Please input an integer p (p=1/rate, only consider p=2 or 3 or 4 cases): default p=3\n'); if isempty(p) p=3; end % input the number of iterations niter = input(' Please enter number of iterations for each frame: default 4\n'); if isempty(niter) niter = 4; end % Number of frame errors to count as a stop criterior. ferrlim = input(' Please enter number of frame errors to terminate: default 5\n '); if isempty(ferrlim) ferrlim = 5; end EbN0db = input(' Please enter Eb/N0 in dB : default [0.5:0.5:2.0]\n'); if isempty(EbN0db) EbN0db = [0.5:0.5:2.0]; end fprintf('\n\n----------------------------------------------------\n'); fprintf(' === Modified Addition MAP decoder === \n'); fprintf(' Frame size = %6d\n',L_frame); fprintf(' Channel condition = %6d\n',channel); fprintf(' co