turbo码的matlab仿真

turbo码在matlab中的编码和仿真，学习的好资料，可供参考

% funcion principal donde se setean los par醡etros de la simulaci髇

function bermaxrand

clear all

% parametros de poly2trellis
% K: Longitud restringida
% G: Codegenerator en octal
% F: Feedback connection en octal

% para el ejemplo de Casti K=3 G=7 F=5
%K=3 ; G=7; F=5;

% Paper de Berrou
K=5; G=21; F=37; % este codificador tiene muy buena performance

% asigno el vector de realimentacion para terminar el trellis del primer
% codificador
assignin（‘base‘‘F‘F）;

trellis = poly2trellis（KGF）;
% asigno la estructura al workspace ‘base‘
assignin（‘base‘‘trellis‘trellis）;

% Genero la matriz de conexiones y salidas del trellis

trellismat;


EbNovector = 0:0.5:2;

iteracionesvector = [8 10];

% Trama de longitud N
N=53*53;
assignin（‘base‘‘N‘N）;

CantTramas=100; % cantidad de tramas que se transmiten por cada una se 
% realiza el proceso de decodificacion

% distintos trazos para diferentes iteraciones en un mismo grafico
trazo = [‘-‘ ‘:‘ ‘-.‘ ‘--‘];

contador = 0;

% vector con el valor de BER para cada EbNo
bervector = zeros（length（EbNovector）1）;

for j = 1:length（iteracionesvector）

    iteraciones = iteracionesvector（j）;
    assignin（‘base‘‘iteraciones‘iteraciones）;
    
    for i = 1:length（EbNovector）

        EbNo =   EbNovector（i）;
        assignin（‘base‘‘EbNo‘EbNo）;
        
        totalerrors = 0;
        
        for k=1:CantTramas
            
            % funcion que simula el esquema Tx-canal-Rx
            turbo;
                        
            % errors es una funcion
            totalerrors = totalerrors + errors; 
            
            contador = contador + 1;
            progreso = 100*contador/（length（iteracionesvector）*...
                length（EbNovector）*CantTramas）
       
        end

        bervector（i） = totalerrors/（N*CantTramas）;

    end

    semilogy（EbNovectorbervector trazo（j））;
    xlabel（‘Eb/No [dB]‘）;
    ylabel（‘BER‘）;
    hold on
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        2011  2006-04-20 01:40  bermaxrand.m
     文件         650  2006-04-11 16:48  beruncoded.m
     文件          37  2005-09-12 17:06  bin2polar.m
     文件        4801  2006-04-11 16:45  decmax.m
     文件         162  2006-03-31 06:28  errors.m
     文件        1332  2009-05-22 21:53  license.txt
     文件          43  2005-09-16 16:19  polar2bin.m
     文件        1200  2006-04-04 06:09  trellismat.m
     文件        4516  2006-04-11 16:46  turbo.m