實(shí)驗(yàn)一信道容量的迭代算法

1、陜西科技大學(xué)實(shí)驗(yàn)報(bào)告班級(jí)： 學(xué)號(hào)： 姓名： 實(shí)驗(yàn)組別： 實(shí)驗(yàn)日期： 2012-12-18 報(bào)告日期： 2012-12-18 成績： 報(bào)告內(nèi)容：（目的和要求、原理、步驟、數(shù)據(jù)、計(jì)算、小結(jié)等）實(shí)驗(yàn)名稱 ：信道容量的迭代算法一、實(shí)驗(yàn)?zāi)康?） 進(jìn)一步熟悉信道容量的迭代算法。2） 學(xué)習(xí)如何將復(fù)雜的公式轉(zhuǎn)化為程序。3） 掌握高級(jí)語言數(shù)值計(jì)算程序的設(shè)計(jì)和調(diào)試技術(shù)。二、實(shí)驗(yàn)要求（1）已知：信源符號(hào)個(gè)數(shù) r 、信宿符號(hào)個(gè)數(shù) s 、信道轉(zhuǎn)移概率矩陣 P。 （2）輸入：任意的一個(gè)信道轉(zhuǎn)移概率矩陣。信源符號(hào)個(gè)數(shù)、信宿符號(hào)個(gè)數(shù)和每 個(gè)具體的轉(zhuǎn)移概率在運(yùn)行時(shí)從鍵盤輸入。（ 3）輸出：信道容量 C。三、實(shí)驗(yàn)原理（流程圖）

2、：1: procedure ChannelCapacity( r, s, ( pji )12: initialize: 信源分布 pi 1, 相對誤差門限 , C r3: repeat4:pi pjirpi pjii15:exp pji log ijj1pirrj1s6:exppji log iji 1 j1rC log expi1j1pjilog ij7: until8: output C9: end procedure流程圖 :四、實(shí)驗(yàn)數(shù)據(jù)(源代碼) ：#include #include #define MAX 100 double Calculate_a(int k,double pa)

3、; double Calculate_C1(double pa,double a); double Calculate_C2(double a);int r,s; double pba33; void main() int i,j;double C1,C2,E; double aMAX,pa3; E=0.000001;r=3;word 文檔 可自由復(fù)制編輯 s=3;pa0=0.2;pa1=0.3; pa2=0.5; pba00=0.5; pba01=0.33333333; pba02=0.16666666; pba10=0.16666666; pba11=0.5; pba12=0.333333

4、333; pba20=0.33333333; pba21=0.16666666; pba22=0.5;dofor(i=0;i=E)double sum=0; for(i=0;ir;i+) sum+=pai*ai;for(i=0;ir;i+) pai=pai*ai/sum;elseprintf( 最佳信源概率 :n); for(i=0;i=E);printf( 信道容量為 :%lfn,C1/log(2); double Calculate_a(int k,double pa)int i,j;double temp,sum2=0; for(j=0;js;j+) double sum1=0; for

5、(i=0;ir;i+)word 文檔 可自由復(fù)制編輯 sum1+=pai*pbaij; temp=pbakj/sum1; temp=log(temp);sum2+=pbakj*temp; return exp(sum2);double Calculate_C1(double pa,double a) int i;double sum=0; for(i=0;ir;i+) sum+=pai*ai; return log(sum);double Calculate_C2(double a)int i;double max=a0; for(i=0;ir;i+) if(maxai) max=ai; return log(max);五、程序運(yùn)行結(jié)