信息論上機實驗報告

PAGEI1-華中科技大學文華學院信息論上機實驗報告專業(yè)年級：10級通信工程**學號：1001********姓名：**指導教師：***2012年目錄TOC\o"1-3"\h\u85831等長編碼 227361.1求信源的熵 2861.2求等長碼碼長 285831.3求等長長等編碼 385831.4等長編碼全部源代碼 3151381.5MATLAB仿真及其截圖 5177822霍夫曼編碼 6180442.1霍夫曼編碼源代碼 6130182.2MATLAB仿真及其截圖 93游程編碼23073 9120523.1游程編碼源代碼 9246013.2MATLAB仿真及其截圖 11196834學習心得 12272865參考書籍 141、用MATLAB實現等長編碼。1.1求信源的熵symsh;s=[0.25,0.25,0.25,0.25];len=size(s,2);ifsum(s)~=1error('sisnotacompleteset');endh=0;fori=1:1:lenh=h-s(i)*log(s(i))/log(2);end1.2求等長碼長r=2;l1=log(len)/log(r);lc=0;ifl1-round(l1)>0lc=round(11)+1;elselc=round(l1);end1.3等長編碼A=zeros(4,lc);fori=1:1:4;s=i;t=lc;forj=1:1:lc;ifs>power(2,t)/2;A(i,lc-t+1)=1;s=s-power(2,t)/2;elseA(i,lc-t+1)=0;endt=t-1;endend1.4全部源代碼s=[0.25,0.25,0.25,0.25];len=size(s,2);ifsum(s)~=1error('sisnotacompleteset');endh=0;fori=1:1:lenh=h-s(i)*log(s(i))/log(2);endr=2;l1=log(len)/log(r);lc=0;ifl1-round(l1)>0lc=round(11)+1;elselc=round(l1);endA=zeros(4,lc);fori=1:1:4;s=i;t=lc;forj=1:1:lc;ifs>power(2,t)/2;A(i,lc-t+1)=1;s=s-power(2,t)/2;elseA(i,lc-t+1)=0;endt=t-1;endend結果：h=2、lc=2、A=000110111.5MATLAB仿真及其截圖2、霍夫曼編碼2.1霍夫曼編碼源代碼A=[0.3,0.2,0.1,0.2,0.2];%信源消息的概率序列A=fliplr(sort(A));%按降序排列T=A;[m,n]=size(A);B=zeros(n,n-1);%空的編碼表（矩陣）fori=1:nB(i,1)=T(i);%生成編碼表的第一列endr=B(i,1)+B(i-1,1);%最后兩個元素相加T(n-1)=r;T(n)=0;T=fliplr(sort(T));t=n-1;forj=2:n-1%生成編碼表的其他各列fori=1:tB(i,j)=T(i);endK=find(T==r);B(n,j)=K(end);%從第二列開始，每列的最后一個元素記錄特征元素在%該列的位置r=(B(t-1,j)+B(t,j));%最后兩個元素相加T(t-1)=r;T(t)=0;T=fliplr(sort(T));t=t-1;endB;%輸出編碼表END1=sym('[0,1]');%給最后一列的元素編碼END=END1;t=3;d=1;forj=n-2:-1:1%從倒數第二列開始依次對各列元素編碼fori=1:t-2ifi>1&B(i,j)==B(i-1,j)d=d+1;elsed=1;endB(B(n,j+1),j+1)=-1;temp=B(:,j+1);x=find(temp==B(i,j));END(i)=END1(x(d));endy=B(n,j+1);END(t-1)=[char(END1(y)),'0'];END(t)=[char(END1(y)),'1'];t=t+1;END1=END;endA%排序后的原概率序列END%編碼結果fori=1:n[a,b]=size(char(END(i)));L(i)=b;endavlen=sum(L.*A)%平均碼長H1=log2(A);H=-A*(H1')%熵P=H/avlen%編碼效率結果：A=0.30000.20000.20000.20000.1000END=[01,10,11,000,001]avlen=2.3000H=2.2464P=0.97672.2MATLAB仿真及其截圖3、游程編碼3.1、游程編碼源代碼I=checkerboard(10,2);[m,n]=size(I);J=[];fori=1:mvalue=I(i,1);num=1;forj=2:nifI(i,j)==valuenum=num+1;elseJ=[Jnumvalue];num=1;value=I(i,j);endendJ=[Jnumvalue00];enddisp('原圖像大?。?)whos('I');disp('壓縮圖像大?。?)whos('J')disp('圖像壓縮比：')disp(m*n/length(J))結果：原圖像大小：NameSizeBytesClassI40x4012800doublearray壓縮圖像大?。篘ameSizeBytesClassJ1x4003200doublearray圖像壓縮比：43.2MATLAB仿真及其截圖信息論學習感想對于信息論上機學習來說，我覺得的一些技巧對自己的幫助還是很大的，特別是Matlab中的一些基礎知識。我自我感覺Matlab學的不好，在開始的時候我只做了一個題目，而且還是百度上拷下來的。不過在后面的時間我發(fā)現自己越做越快，首先可能克服了難的心理，然后跟自己好好看參考書籍也有一定關系。所以凡事主要靠自己，要給自己充裕的時間處理事情。不要等到快要交作業(yè)了，快要答辯了，才發(fā)現自己沒做，那時就會手忙腳亂，叫天不應，叫地不靈。只能祈求哪位老師或大哥大姐給個例子或寫個程序吧，百度上有這個題目。不能把太多的希望都寄托在別人身上，也不是每次都那么幸運的能百度出來。其次自己手里還有幾本參考資料吧，不說非常好的，學校圖書館借的就行了，我們的課本講的東西太少了，舉得例子也太少了，做我們的參考遠遠不夠。在做題之前一定要先把基礎問題先搞懂，不能一開始就想著怎么把任務完成，我開始做的時候，怎么定義函數都不知道，而且好多函數都不了解怎么用。其實如果要做的快一點可以先把一些常用的函數記住，然后做哪個題就重點看哪一個，可以避免無目的的看書浪費時間。我在做題時遇到了很多問題，首先現在是網絡時代了，我第一個想到的就是百度，那個函數不會用，那個函數不會定義，以至于出錯了為什么會出現錯誤都百度過。另外一個就是找同學幫忙，我學的只能算一個中等偏下的水平，比我學的好的同學很多，出現了一些問題就直接向他們求教，雖然有些時候他們也不一定能解答出來，但最少也能給我一些提示。最后有時候也靠自己，慢慢的想，慢慢的翻資料一遍一遍的輸入，校正，最后就做出來了。雖然有些題目我也不知道自己的對不對，但是通過自己動手把它做出結果了我就覺得有收獲了。Matlab的課比較少，課堂上學到的東西也不是特別的多，但是通過實驗，我自學了許多東西，也完成了