jgraving

function [C r] = BlahutArimoto(p)

disp('BlahutArimoto')

% Capacity of discrete memoryless channel
% Blahut-Arimoto algorithm

% Input
% p: m x n matrix
% p is the transition matrix for a channel with m inputs and n outputs

jgraving

/*
2D Angle Interpolation (shortest distance)

Parameters:
a0 = start angle
a1 = end angle
t = interpolation factor (0.0=start, 1.0=end)

Benefits:
1. Angles do NOT need to be normalized.