作业帮急求信道容量matlab编程代码!已知一个信道的信道转移矩阵为 p用Matlab编写函数求信道容量.已知一个信道的信道转移矩阵为 p,p为一个两行三列矩阵
来源：学生作业帮助网 编辑：作业帮 时间：2024/05/06 20:30:49
急求信道容量matlab编程代码!已知一个信道的信道转移矩阵为 p用Matlab编写函数求信道容量.已知一个信道的信道转移矩阵为 p,p为一个两行三列矩阵
急求信道容量matlab编程代码!已知一个信道的信道转移矩阵为 p用Matlab编写函数求信道容量.
已知一个信道的信道转移矩阵为 p,p为一个两行三列矩阵
急求信道容量matlab编程代码!已知一个信道的信道转移矩阵为 p用Matlab编写函数求信道容量.已知一个信道的信道转移矩阵为 p,p为一个两行三列矩阵
%亲自验证,绝对可用!
% 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)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)