metric_acc_cmp_sel.m

本人编写的TCM解调源程序共享

function [sur_acc_metric, sur_path_index] = metric_acc_cmp_sel(acc_metric, branch_metric, state_trans, stat_trans_input)

% sqwang's source code, used for academic research, not commerce.

[num_row, num_col] = size(state_trans);

stat_trans_input = stat_trans_input + 1;  % Add 1 for matlab index convenience

tentative_acc_metric = zeros(num_col, num_row);
sur_acc_metric = zeros(num_row, 1);
sur_path_index = zeros(num_row, 1);

for k = 1 : num_row
    % State metric accumulation
    for l = 1 : num_col  % num_col branch to each state
        tentative_acc_metric(l, k) = branch_metric(state_trans(k, l), stat_trans_input(k, l)) + acc_metric(state_trans(k, l));  
    end
end

[max_acc_metric, sur_path_index] = max(tentative_acc_metric);
sur_acc_metric = max_acc_metric';