一般离散无记忆信道容量的迭代计算_第1页
一般离散无记忆信道容量的迭代计算_第2页
一般离散无记忆信道容量的迭代计算_第3页
一般离散无记忆信道容量的迭代计算_第4页
一般离散无记忆信道容量的迭代计算_第5页
已阅读5页，还剩1页未读， 继续免费阅读

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

一般离散无记忆信道容量的迭代计算信道容量的迭代算法1信道容量的迭代算法的步骤一、用了matlab实现DMC容量迭代的算法如下：第一步:首先要初始化信源分布：即选取一个精度，本次中我选deta=0.000001。第二步：。第三步：第四步：第五步：若，则执行k=k+1，然后转第二步。直至转移条件不成立，接着执行下面的程序。第六步：输出迭代次数k和和，程序终止。i=1:1:r;%设置循环首项为1，公差为1，末项为r（Q的行数）的循环p(i)=1/r;%原始信源分布r个信源,等概率分布disp('原始信源分布：')p(i)E=repmat(p',1,s);%把r个等概率元素组成一列，复制为s列fork=1:1:1/detam=E.*P;%m=p.*E;%后验概率的分子部分a=sum(m);%把得到的矩阵m每列相加之和构成一行su1=repmat(a,r,1);%把得到的行矩阵a复制r行,成一新矩阵sul,后验概率的分母部分t=m./su1;%后验概率矩阵n=exp(sum(P.*log(t),2));%信源分布的分子部分su2=sum(n);%信源分布的分母部分p=n/su2;%信源分布E=repmat(p,1,s);C(k+1)=log(sum(exp(sum(P.*log(t),2))))/log(2);kk=abs(C(k+1)-C(k))/C(k+1);if(kk<=deta)break;enddisp('迭代次数：k='),disp(k)enddisp('最大信道容量时的信源分布：p='),disp(p')disp('最大信道容量：C='),disp(C(k+1))运行结果及分析结果分析：这两组数据都是我随机选的，都是选的信源个数为2，信宿的个数为3，选用的精度为0.000001。然后输入信道转移矩阵P，执行，将得到的结果代入，得到后，再进行的判断，这个条件满足时继续求，再依次往下计算，直至这个条件不满足，然后输出k和和。总的来说这不过是将矩阵不断的代入公式，当最后的精度不大于0.000001.输出k和和。第二组数据：4.心得体会

在此次matlab编程实现一般离散无记忆信道容量的迭代算法的过程中，开始我觉得这十分的困难，经过网上查阅资料和请教同学，最后才对本次设计有了一定的理解，详细理解了信道容量

人人文库> 全部分类> 行业资料 > 信息产业

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论