信道容量迭代計(jì)算實(shí)驗(yàn)報(bào)告

信道容量迭代計(jì)算實(shí)驗(yàn)報(bào)告（總3頁(yè)）--本頁(yè)僅作為文檔封面，使用時(shí)請(qǐng)直接刪除即可----內(nèi)頁(yè)可以根據(jù)需求調(diào)整合適字體及大小-

信道容量迭代計(jì)算實(shí)驗(yàn)報(bào)告王升信科1002

信道容量迭代計(jì)算實(shí)驗(yàn)報(bào)告一、實(shí)驗(yàn)?zāi)康模毫私庑诺廊萘康亩x和計(jì)算方法，能編寫出正確的程序進(jìn)行迭代計(jì)算得出信道容量。二、實(shí)驗(yàn)要求：1）輸入：輸入信源個(gè)數(shù)、信宿個(gè)數(shù)和信道容量的精度，程序能任意生成隨機(jī)的信道轉(zhuǎn)移概率矩陣。2）輸出：輸出最佳信源分布和信道容量。三、實(shí)驗(yàn)環(huán)境：Matlab、實(shí)驗(yàn)原理：p（硝"3輸入廠 p（bt\a.）悶=^rLp?t 乞訊叩鞏眄丨碼）Jd……『）Ciit十L/i)二】n工p(at疣C\n+],ra}-ln(ttiska)C(n11.in匚J(inC=C(]i>l五、源程序代碼：clear;r=input（'輸入信源個(gè)數(shù)：’）；s=input（'輸入信宿個(gè)數(shù)：'）;deta=input（'輸入信道容量的精度：’）；Q=rand（r,s）;%創(chuàng)建m*n隨機(jī)分布矩陣A=sum（Q,2）;B=repmat（A,1,s）;

disp('信源轉(zhuǎn)移概率矩陣：’)，p=Q./B %信源轉(zhuǎn)移概率矩陣i=1:1:r;q(i)=1/r;disp('原始信源分布:'),qc=-10e-8;C=repmat(q',1,s);fork=1:1:100000m=p.*C;%后驗(yàn)概率的分子部分a=sum(m);%后驗(yàn)概率的分母部分su1=repmat(a,r,1);t=m./su1; %后驗(yàn)概率矩陣D=exp(sum(p.*log(t),2));%信源分布的分子部分su2=sum(D); %信源分布的分母部分q=D/su2; %信源分布C=repmat(q,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<=break;endenddisp('最大信道容量時(shí)的信源分布：q='),disp(q')disp('最大信道容量：c='),disp(c(k+1))六、實(shí)驗(yàn)結(jié)果：輸入ft源個(gè)埶：2輸入信宿亍埶：3輸入信道容里的精度：0.000001信源轉(zhuǎn)移拯率矩蹲：P=0.4392 0.3572 0.19350.3146 □.0320 □.5335原始信源分布：0.=0.50