Capacidad de canal

Cálculo de Capacidad de Canal
Introducción
La capacidad de canal discreto, determinado por , se determina conociendo la máxima información mutua entre dos variables aleatorias, , es decir,donde , se escribe como sigue,

Lo que a continuación se expone, es precisamente la capacidad de canal, para lo cual se considera la siguiente matriz de transición:

Sin embargo, las probabilidadesiniciales, o a priori, se desconocen, éstas serán propuestas a fin de que la información mutua, , sea la máxima posible, y en consecuencia conocer la capacidad de canal. Aquí conviene precisar que, sila matriz de transición P fuese simétrica y las probabilidades iniciales fuesen equivalentes, se tendría de manera adyacente un máximo en la entropía, por tanto en la información mutua, y sin más, sesabría de buena tinta la capacidad de canal discreto; éste no es el caso, sin embargo.
Se ha hecho uso de una herramienta computacional para realizar el cálculo de la capacidad de canal, correspondeel caso a MATLAB, el código es el siguiente:
% -----------------------------

clc;
clear all;
close all;

% -----------------------------

P = [2/10 3/10 1/2; 1/10 1/2 2/5; 4/10 1/40.35];

% -----------------------------



dprob = 0.1:0.1:1;
v = zeros(3,power(length(dprob),3));
infm = zeros(1,power(length(dprob),3));

m = 1;
for n1 = 1:length(dprob)for n2 = 1 : length(dprob)
for n3 = 1 : length(dprob)
if ( ( dprob(n1) + dprob(n2) + dprob(n3) ) == 1 )
v(:,m) = [dprob(n1),dprob(n2),dprob(n3)];m = m + 1;
end
end
end
end

% -----------------------------

for h = 1:length(v);
for i = 1:3
for j = 1:3
s = 0;for k = 1:3
s = v(k,h)*P(k,j) + s;
end
infm(h) = v(i,h)*P(i,j)*log2( P(i,j) / s ) + infm(h);
end
end
end

M = max(infm);...