迭代算法计算信道容量m=8;n=5;%n行m列信道矩阵 e=0.0001;%精度 c1=1;% c2=0;% B=zeros(n,m);%计算 D=zeros(n,m);%计算 P=abs(rand(n,m));%信道矩阵元素全为正 [n,m]=size(P); A=ones(1,n)./n;%初始输入概率分布 %------生成一个 随机的信道矩阵,满足 , for i=1:1:n if norm(P(i,:),1)>0 P(i,:)=P(i,:)/norm(P(i,:),1); else P(i,:)=1/n; end end %...
m=8;n=5;%n行m列信道矩阵 e=0.0001;%精度 c1=1;% c2=0;% B=zeros(n,m);%计算 D=zeros(n,m);%计算 P=abs(rand(n,m));%信道矩阵元素全为正 [n,m]=size(P); A=ones(1,n)./n;%初始输入概率分布 %------生成一个 随机的信道矩阵,满足 , for i=1:1:n if norm(P(i,:),1)>0 P(i,:)=P(i,:)/norm(P(i,:),1); else P(i,:)=1/n; end end %------迭代算法 while abs(c1-c2)>e %------计算 for i=1:1:m E=P(:,i)./(A*P(:,i));%计算 D(:,i)=log(E);%计算 End %-----对D矩阵中的-Inf元素全取为0，E中可能有零元素存在，0的对数是负无穷，%-----没有意义 for i=1:1:n for j=1:1:m if D(i,j)==-Inf D(i,j)=0 ; end end end %----------------- B=P*D';% 在B的对角线上 a=exp(diag(B));%提取对角线元素,计算% ， % 矩阵 c1=log(A*a);%计算 c2=log(max(a));%计算 A=A.*a'./(A*a);%计算 End c1/log(2)%转换单位为比特/符号
本文档为【迭代算法计算信道容量】，请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑，
图片更改请在作品中右键图片并更换，文字修改请直接点击文字进行修改，也可以新增和删除文档中的内容。
该文档来自用户分享，如有侵权行为请发邮件ishare@vip.sina.com联系网站客服，我们会及时删除。
[版权声明] 本站所有资料为用户分享产生，若发现您的权利被侵害，请联系客服邮件isharekefu@iask.cn，我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权，请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传，仅限个人学习分享使用，禁止用于任何广告和商用目的。
浙公网安备 33021202002483