作业帮 > 综合 > 作业

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

来源:学生作业帮 编辑:灵鹊做题网作业帮 分类:综合作业 时间:2024/05/08 14:39:12
急求信道容量matlab编程代码!已知一个信道的信道转移矩阵为 p用Matlab编写函数求信道容量.
已知一个信道的信道转移矩阵为 p,p为一个两行三列矩阵
急求信道容量matlab编程代码!已知一个信道的信道转移矩阵为 p用Matlab编写函数求信道容量.
%亲自验证,绝对可用!
% 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)