2022年2022年實(shí)驗三香農(nóng)編碼

1、精選學(xué)習(xí)資料 - - - 歡迎下載信息論與編碼試驗報告院系：哈爾濱理工高校榮成校區(qū)專業(yè)：電子信息工程學(xué)號：姓名：日期： 2021年 6 月 16 日精品學(xué)習(xí)資料精選學(xué)習(xí)資料 - - - 歡迎下載香農(nóng)編碼信息論與編碼第三次試驗報告一.試驗?zāi)康暮腿蝿?wù)1. 懂得信源編碼的意義；2. 熟識 matlab程序設(shè)計；3. 把握香農(nóng)編碼的方法及運(yùn)算機(jī)實(shí)現(xiàn)；4. 對給定信源進(jìn)行香農(nóng)編碼,并運(yùn)算編碼效率；二.試驗原理介紹給定某個信源符號的概率分布,通過以下的步驟進(jìn)行香農(nóng)編碼1.信源符號按概率從大到小排列；p1p2.pn精品學(xué)習(xí)資料精選學(xué)習(xí)資料 - - - 歡迎下載2.確定滿意以下不等式的整數(shù)碼長k i 為精品學(xué)

2、習(xí)資料精選學(xué)習(xí)資料 - - - 歡迎下載- lbp i k ilb pi 13.為了編成唯獨(dú)可譯碼,運(yùn)算第i 個消息的累加概率：l 1pip ak k 1精品學(xué)習(xí)資料精選學(xué)習(xí)資料 - - - 歡迎下載4.將累加概率pi 變換成二進(jìn)制數(shù)；精品學(xué)習(xí)資料精選學(xué)習(xí)資料 - - - 歡迎下載精品學(xué)習(xí)資料精選學(xué)習(xí)資料 - - - 歡迎下載5.取pi 二進(jìn)制數(shù)的小數(shù)點(diǎn)后k i 位即為該消息符號的二進(jìn)制碼字；精品學(xué)習(xí)資料精選學(xué)習(xí)資料 - - - 歡迎下載三.試驗設(shè)備1.運(yùn)算機(jī)2.matlab2021四.編碼程序n=input'n='%輸入信源符號的個數(shù)s=0;l=0;%對行的概率進(jìn)行初始化%對

3、列的概率初始化精品學(xué)習(xí)資料精選學(xué)習(xí)資料 - - - 歡迎下載h=0;%對信息熵初始化 fori=1:npi=input'p='%輸入信源符號概率分布矢量 pi<1;%輸入概率小于1s=s+pi;%行概率求和h=h+-pi*log2pi;%運(yùn)算信源信息熵 endifabss-1>0%如行概率和>1 就報錯error' 不符合概率分布' endfori=1:n-1%for循環(huán)按概率分布大小對信源排序forj=i+1:nifpi<pj%冒泡排序 m=pj;pj=pi; pi=m;endend endfori=1:na=-log2pi;%確定滿意

4、不等式的整數(shù)a,并領(lǐng) a 為第 i個碼字的長度ifmoda、1=0%a 對 1 取余后賦值給a w=a;elsew=fixa+1;%運(yùn)算各信源符號的碼長 endl=l+pi*w;%運(yùn)算平均碼長 endr=l;%香農(nóng)編碼的信息率（這里l=1、m=2）n=h/r;%運(yùn)算編碼效率p1=0;%初始 1 列為 0；fori=2:npi=0;%初始 1 列為 0 ；累加forj=1:i-1pi=pi+pj;%從 2 列開頭累加end end%運(yùn)算累加概率fori=1:nforj=1:wwi、j=fixpi*2;%對每次的p（ i ）取碼字 、 長度與 a 相同 pi=pi*2-fixpi*2;%將累加概率轉(zhuǎn)化為li位二進(jìn)制碼字 endend精品學(xué)習(xí)資料精選學(xué)習(xí)資料 - - - 歡迎下載w%顯示碼字r%顯示平均碼長n%顯示編碼效率h%顯示信息熵 h五.試驗結(jié)果>> xiangnong n=6p=0.25 p=0.25 p=0.20 p=0.15 p=0.10 p=0.05w =000000100010000101101101111110r =2.7000n =0.8975h =2.4232精品學(xué)習(xí)資料精選學(xué)習(xí)資料 - - - 歡迎下載六.試