Question

急求信道容量matlab编程代码！！已知一个信道的信道转移矩阵为 p用Matlab编写函数求信道容量.

已知一个信道的信道转移矩阵为 p，p为一个两行三列矩阵

Answer 1

%亲自验证，绝对可用！
% Matlab实现离散信道容量的迭代算法% 功能：利用迭代算法计算离散信道的容量
%参数解释
%C：信道容量
%P：转移概率矩阵
%B：中间变量矩阵
%e: 信道容限，一般选0.00001
%X：输入概率分布
%n：迭代次数

function channel_cap(P, e)
n=0;
C=0;
C_0=0;
C_1=0;
[r,s]=size(P);
for i=1:r
if(sum(P(i,:))~=1)%检测概率转移矩阵是否行和为1.
error('概率转移矩阵输入有误！！')
return;
end
for j=1:s
if(P(i,j)<0||P(i,j)>1)%检测概率转移矩阵是否负值或大于1
error('概率转移矩阵输入有误！！')
return;
end
end
end
X=ones(1,r)/r;
A=zeros(1,r);
B=zeros(r,s);
while(1)
n=n+1;
for i=1:r
for j=1:s
B(i,j)=log(P(i,j)/(X*P(:,j))+eps);
end
A(1,i)=exp(P(i,:)*B(i,:)');
end
C_0=log2(X*A');
C_1=log2(max(A));
if (abs(C_0-C_1)<e) %满足迭代终止条件停止迭代
C=C_0;
fprintf('迭代次数: n=%d\n',n)
fprintf('信道容量: C=%f比特/符号\n',C)
break; %满足后输出结果并退出
else
X=(X.*A)/(X*A');
continue;
end
end

Answer 2

这个东西自己写