信息论与编码实验报告

《信息论与编码实验报告》由会员分享，可在线阅读，更多相关《信息论与编码实验报告（26页珍藏版）》请在装配图网上搜索。

1、实用文档 实验报告 课程名称： 信息论与编码 姓 名： 系： 专 业： 年 级： 学 号： 指导教师： 职 称： 年 月 日 . 实用文档 目 录 实验一 信源熵值的计算 1 实验二 Huf

2、fman 信源编码 5 实验三 Shannon编码 9 实验四 信道容量的迭代算法 12 实验五 率失真函数 14 实验六 差错控制方法 19 实验七 汉明编码 22 . 实用文档 实验一 信源熵值的计算 一、实验

3、目的 1 进一步熟悉信源熵值的计算 2 熟悉 Matlab 编程 二、实验原理 熵（平均自信息）的计算公式 q 1 q H (x) pi log 2 pi pi log 2 i 1 pi i 1 MATLAB实现： HX sum( x.* log 2 ( x)) ；或者 h h x(i ) * log 2 ( x(i )) 流程：第一步：打开一个名为“ nan311”的 TXT文档，读入一篇英文文章存入一个数组 temp，为了程序准确性将所读内容转存到另一个数组 S，计算该数组中每个字母与空格的

4、出现次数 ( 遇到小写字母都将其转化为大写字母进行计数 ) ，每出现一次该字符的计数器 +1； 第二步：计算信源总大小计算出每个字母和空格出现的概率； 最后，通过统计数据和信息熵公式计算出所求信源熵值（本程序中单位为奈特 nat ）。 程序流程图： . 实用文档 三、实验内容 1、写出计算自信息量的 Matlab 程序 2、已知：

5、信源符号为英文字母（不区分大小写）和空格。 输入：一篇英文的信源文档。 输出：给出该信源文档的中各个字母与空格的概率分布，以及该信源的熵。 四、实验环境 Microsoft Windows 7 Matlab 6.5 五、编码程序 #include"stdio.h" #include #include #define N 1000 int main(void) { char s[N]; int i,n=0; float num[27]={0}; double r

6、esult=0,p[27]={0}; FILE *f; char *temp=new char[485]; f=fopen("nan311.txt","r"); while (!feof(f)) { fread(temp,1, 486, f);} fclose(f); s[0]=*temp; for(i=0;i

7、; else if(s[i]>='a'&&s[i]<='z') num[s[i]-97]++; else if(s[i]>='A'&&s[i]<='Z') . 实用文档 num[s[i]-65]++; } printf(" 文档中各个字母出现的频率 :\n"); for(i=0;i<26;i++) { p[i]=num[i]/strlen(s); printf("%3c:%f\t",i+65,p[i]); n++; if(n==3) { printf("\n"); n=0;

8、 } } p[26]=num[26]/strlen(s); printf(" 空格 :%f\t",p[26]); printf("\n"); for(i=0;i<27;i++) { if (p[i]!=0) result=result+p[i]*log(p[i]); } result=-result; printf(" 信息熵为 :%f",result); printf("\n"); return 0; } 六、求解结果

9、 其中 nan311.txt 中的文档如下： There is no hate without fear. Hate is crystallized fear, fear ’dividend,s fear objectivized. We hate what we fear and so where hate is, fear is lurking. Thus we hate . 实用文档 what threatens our person, our vanity and our dreams and plans for

10、 ourselves. If we can isolate this element in what we hate we may be able to cease from hating. 七、实验总结 通过这次实验， 我们懂得了不必运行程序时重新输入文档就可以对文档进行 统计，既节省了时间而且也规避了一些输入错误。 在实验中，我们进一步了解到 信源熵的计算，理论和实践的结合让我们对这个知识点了解的更加深刻了。

11、 . 实用文档 实验二 Huffman 信源编码 一、实验目的 1．理解信源的最优变长编码的基本思想。 2．熟练掌握 Huffman 信源编码方法。 二、设计原理 设信源 S={s1,s2, ..,sq}，其对应的概率分布为 P(si)={p1,p2,p3, .，,pq}则其编 码步骤如下： (1)将

12、 q 个信源符号按递减方式排列。 (2)用 0、1 码符分别表示概率最小的两个信源符号，并将这两个符号合并成 一个新的符号，从而得到 q-1 个符号的新信源成为 S 信源的缩减信源 S1。 (3)将缩减信源 S1 中的符号仍按递减顺序排列，再将最小两个概率相加，合 并成一个符号，并分别用 0、1 码表示，这样有形成了 q-2 个缩减信源 S2。 (4)依次继续下去，直到缩减信源只剩下两个符号为止， 将最后两个符号用 0、 1 分别表示。 (5)从最后一次缩减信源开始，向前返回，沿信源缩减过程的反方向取出所 编的马元。 三、实验内容 计

13、算定信源和输入信号字母表的 Huffman 编码，并计算 Huffman 编码的平均 码长。实验具体要求如下： 信源字母表的概率分布为： P={ 0.15,0.12,0.2,0.08,0.04,0.18,0.02,0.09,0.04,0.02,0.06} 输入信号字母表： U={0,1,2} ； 1. 独立设计信源和输入信号字母表进行Huffman 编码 , 其中信源字母表元 素个数要求是 8 以上，信号字母表元素个数是 2 以上； 2. 输出 Huffman 编码的平均码长。 四、实验环境 Microsoft Windows 7

14、 Matlab 6.5 五、编码程序 MATLAB编码： . 实用文档 function[h,L]=huffman(p,r) %变量 p 为符号出现概率所组成的概率向量 %返回值 h 为利用 Huffman 编码算法编码后最后得到编码结果 %返回值 L 为进行 Huffman 编码后所得编码的码字长度 if length(find(p<0))~=0 error('Not a prob.vector,negative component(s)'); end %判断概率向量中是否有 0 元素，有 0 元素程

15、序显示出错，终止运行 if (sum(p,2)>1) error('Not a prob.vector,components do not add up to 1'); end %判断所有符号出现概率之和是否大于 1，如果大于 1 程序显示出错，终止 运行 a=length(p); %测定概率向量长度，将长度值赋给变量 n k=fix((a-1)/(r-1)); l1=a-k*r+k; q=zeros(1,a); m=zeros(k+1,a); mp=m; q=p; [m(1,:),mp(1,:)]=sort(q);

16、 if (l1>1) s=sum(m(1,1:l1),2); q=[s,m(1,(l1+1):a),ones(1,l1-1)]; [m(2,:),mp(2,:)]=sort(q); else m(2,:)=m(1,:); mp(2,:)=1:1:a; end for i=3:k+1 s=sum(m(i-1,1:r),2); q=[s,m(i-1,r+1:a),ones(1,r-1)]; [m(i,:),mp(i,:)]=sort(q); end n1=m; n2=mp; for i=1:k+1 n

17、1(i,:)=m(k+2-i,:); n2(i,:)=mp(k+2-i,:); end m=n1; mp=n2; c=cell(k+1,a); for j=1:r . 实用文档 c{1,j}=num2str(j-1); end for i=2:k p1=find(mp(i-1,:)==1); for j=1:r c{i,j}=strcat(c{i-1,p1},int2str(j-1)); end for j=(r+1):(p1+r-1) c{i,j}=c{i-1,j-r

18、}; end for j=(p1+r):a c{i,j}=c{i-1,j-r+1}; end end if l1==1 for j=1:a c{k+1,j}=c{k,j}; end else p1=find(mp(k,:)==1); for j=1:l1 c{k+1,j}=strcat(c(k,p1),int2str(j-1)); end for j=(l1+1):(p1+l1) c{k+1,j}=c{k,mp(1,j-l1)}; end for j=(p1(1)+l1+1):a

19、c{k+1,j}=c{k,mp(1,j-l1+1)}; end end for j=1:a l(j)=length(c{k+1,j}); end h=cell(1,a); for j=1:a h{1,j}=c{k+1,j}; end L=sum(l.*m(k+1,:)); % 求平均码长 2、在 MATLAB 命令窗口中输入： p=[0.15,0.12,0.2,0.08,0.04,0.18,0.02,0.09,0.04,0.02,0.06]; r=3; [h,L]=huffman(p,r).

20、 . 实用文档 六、运行结果 得出的结论为： 概率 编码 概率 编码 0.15 2120 0.02 11 0.12 2121 0.09 12 0.2 2122 0.04 20 0.08 210 0.02 22 0.04 211 0.06 0 0.18 10 L=2.0600 七、实验总结

21、 在 huffman 编码的过程中，我们运用了平时熟悉的数学软件 MATLAB的运行来实现，把书本上 huffman 的算法运用编程来实现。通过这次实验，使我更加清晰地理解 huffman 编码的原理及实现过程， 并且能够在 MATLAB中熟练地进行编码运行。 . 实用文档 实验三 Shannon编码 一、实验目的 1、熟悉离散信源的特点； 2、学习仿真离散信源的方法 3、学习离散信源平均信息量的计算方法 4、熟悉 Matlab 编程 二、实验原理 给定某个信源符号的概率分布，通过以下的步骤进行香农编

22、码 1、信源符号按概率从大到小排列； p1 p2 ....... pn 2、确定满足下列不等式的整数码长 K i 为 lb ( pi ) Ki lb ( pi ) 1 3、为了编成唯一可译码，计算第 i 个消息的累加概率： i 1 Pi p(ak ) k 1 4、将累加概率 Pi 变换成二进制数； 5、取 Pi 二进制数的小数点后 Ki 位即为该消息符号的二进制码字。 三、实验内容 1、写出计算自信息量的 Matlab 程序 2、写出计算离散信源平均信息量的 Matlab 程序。 3、将程序在计算机上仿真

23、实现，验证程序的正确性并完成习题。 四、实验环境 Microsoft Windows 7 Matlab 6.5 五、编码程序 计算如下信源进行香农编码，并计算编码效率： X a0 a1 a2 a3 a4 a5 a6 P 0.2 0.19 0.18 0.17 0.15 0.1 0.01 MATLAB程序： (1) a=[0.2 0.18 0.19 0.15 0.17 0.1 0.01]; k=length(a);y=0; for i=1:k-1 . 实用文档 for n=

24、i+1:k if (a(i)

25、'),disp(a(m)) disp('? oíó?á1?'),disp(s(m)) disp(' à±????êy'),disp(b(m)) disp(' ???×±à??'),disp(z') end (2) function y=b2d10(x) for i=1:8 temp=x.*2; if(temp<1) y(i)=0; x=temp; else x=temp-1; y(i)=1; end end (3) p=[0.2 0.19 0.18 0.17 0.15 0.1 0.01]; sum=0;sum1=0;

26、 for i=1:7 a(i)=-log2(p(i)); K(i)=ceil(a(i)); R(i)=p(i)*K(i); sum=sum+R(i); c(i)=a(i)*p(i); . 实用文档 sum1=sum1+c(i); end K1=sum; H=sum1; Y=H/K1; disp('??? ùD??￠á?'),disp(H) disp('??? ù??3¤'),disp(K1) disp(' à±??D§?ê'),disp(Y) 六、实验结果 输出结果为： 出事概率 0.2000

27、，求和结果 0，编码位数 3，最终编码 000 出事概率 0.1900 ，求和结果 0.2000 ，编码位数 3，最终编码 001 出事概率 0.1800 ，求和结果 0.3900 ，编码位数 3，最终编码 011 出事概率 0.1700 ，求和结果 0.5700 ，编码位数 3，最终编码 100 出事概率 0.1500 ，求和结果 0.7400 ，编码位数 3，最终编码 101 出事概率 0.1000 ，求和结果 0.8900 ，编码位数 4，最终编码 1110 出事概率 0.0100 ，求和结果 0.9900 ，编码位数 7，最

28、终编码 1111110 编码效率： 平均信息量 2.6087 平均码长 3.1400 编码效率 0.8308 七、实验总结 通过本次的实验，掌握了 Shannon 编码的实验原理以及编码过程。 Shannon 编码中，对概率的排序是最基本的， 如果没有将其按照从大到小的顺序排序， 则经过 MATLAB的程序运行后，将出现错误。在运用 MATLAB编程的过程中，调用 了各种函数，实现了编程。通过与队友的讨论，不但让我们更快的完成 MATLAB 编码，也深深体会到只有将大家的智慧融合起来，才能更快更好的解决难题。

29、 . 实用文档 实验四 信道容量的迭代算法 一、实验目的 1、进一步熟悉信道容量的迭代算法； 2、学习如何将复杂的公式转化为程序； 3、熟悉程序设计语言的数值计算程序和调试技术。 二、实验原理 (1)初始化信源分布 P(0 ) P1,P2 , , Pr (一般初始化为均匀分布 ),置迭代计数 器 k=0，设信道容量相对误差门限为 ， >0，可设 ; (2) ( k ) Pij Pi (k ) i 1,2, , s

30、 ji Pij ( k) i Pi exp Pij ln ( k) ji (3) Pi( k 1) j i 1,2, r exp Pij ln ( k ) ji

31、 i j (4) C ( k 1) ln exp Pij ln ( k ) ji i j (5)如果 C ( k 1) C ( k) ，转向 (7)； C ( k 1) (6)置迭代序号 k 1 k ，转向 (2)； (7)输出 Pi (k 1) 和 C (k 1) 的结果； (8)停止。 三

32、、实验内容 1、已知：信源符号个数 r 、新宿符号个数 s、信道转移概率矩阵 P； 2、输入：任意的一个信道转移概率矩阵，信源符号个数、信宿符号个数和 每一个具体的转移概率在运行时从键盘输入； 3、 输出：最佳信源分布 P*,信道容量 C。 四、实验环境 Microsoft Windows 7 、 Matlab 6.5 五、编码程序 . 实用文档 aa.m文件： clear; r=input('输入信源个数： '); s=input('输入信宿个数： '); deta=inp

33、ut('输入信道容量的精度： '); Q=rand(r,s); %创建 m*n 随机分布矩阵 A=sum(Q,2); B=repmat(A,1,s); disp('信源转移概率矩阵 :'),p=Q./B %信源转移概率矩阵 i=1:1:r; q(i)=1/r; disp('原始信源分布： '),q c=-10e-8; C=repmat(q',1,s); for k=1:1:100000 m=p.*C; %后验概率的分子部分 a=sum(m); %后验概率的分母部分 su1=repmat(a,r,1); t=m.

34、/su1; %后验概率矩阵 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<=0.000001) break; end end disp('最大信道容量时的信源分布： q='),disp(q') disp('最大信道容量： c=')

35、,disp(c(k+1)) 六、实验结果结果 1）检验：运行 aa.m 输入信源的个数： 2 输入信宿的个数： 3 输入信道容量的精度： 0.000001 信宿转移概率矩阵 :p =0.5000 0.30000.2000 0.3000 0.5000 0.2000 原始信源分布： q = 0.5000 0.5000 最佳信源分布： q= 0.5000 0.5000 . 实用文档 最大信道容量： c= 0.0365 2）计算信源个数为 3，

36、信宿个数为 5 的信道容量： 运行 aa.m 输入信源的个数： 3 输入信宿的个数： 5 输入信道容量的精度： 0.000001 信宿转移概率矩阵 :p =0.0484 0.1385 0.3058 0.2845 0.2227 0.2104 0.2471 0.1077 0.3762 0.0585 0.3430 0.0800 0.1808 0.3428 0.0534 原始信源分布： q = 0.3333 0.3333 0.3333 最佳信源分布： q =0.4691 0.17

37、94 0.3515 最大信道容量： c =0.1559 七、实验总结 通过实验，我们对信道容量的理解更加深刻了。 信道容量是指信道能无错误传送的最大信息率。 信道的输入、 输出都取值于离散符号集， 且都用一个随机变量来表示的信道就是离散单符号信道。 由于信道中存在干扰， 因此输入符号在传输中将会产生错误， 这种信道干扰对传输的影响可用传递概率来描述。 为了评价实际信道的利用率， 应具体计算已给信道的容量。 这是一个求最大值的问题。 由于互信息对输入符号概率而言是凸函数， 其极值将为最大值， 因此这也就是求极值的问题。对于离散信道， P(x)是一组数，满足

38、非负性和归一性等条件，可用拉格朗日乘子法求得条件极值。对于连续信道， P(x)是一函数，须用变分法求条件极值。 实验过程中，我们虽然也遇到了很多困难，但也正是因为如此，我们才能发现自己基础的薄弱点，学的更有方向。对于编程方面，我们也有了很大的提升。 实验五 率失真函数 . 实用文档 一、实验目的 验证率失真函数的极值特性，理解相关参数的变化对率失真函数的影 响。 二、实验原理 （1）输入 S， d 的初始值、条件概率、输出

39、分布等值； M （2）计算输出分布 q j pii pji j ； i 1 （3）进入迭代，标志为 0 或者误差大于指定 eps 则迭代，否则退出迭代； M M log( pji ij ) ； （4）计算一个互信息 I (qj ; pji ) pii pji ij i 1 j 1 q j （5）计算一个条件概率分布 P( j; i) q j eSdij ； M qk eSdik k 1 M （6）重算一次（ 4），并计算

40、D pii pji ij dij ； i 1 （ 7）重算（ 3） -（ 7）步骤，直到退出迭代； 三、实验环境 Microsoft Windows 7 、 Visual Studio 2005 profession 四、编码程序 #include #include #include using namespace std; //Define some global var const int M = 10; //M 元信源 const doub

41、le S = -50; //迭代算法中的中间量， S 越小，允许最大 失真度 D 越小，当 S 很小时（例如 -100），R(D)=H(X) static int d[M][M]; //失真函数 static double q[M], Pji[M][M]; //输出分布和条件概率分布 static double Pi[M] = {0.4, 0.1, 0.25, 0.1, 0.05, 0.05, 0.01, 0.02, 0.005, 0.015}; //初始化信源的概率分布 const int systemDefine = 2; //定义进制（默认为 2 进制

42、，结果为 bit， 为 e 时，结果为 nat） const double eps = 1e-8; //允许误差 . 实用文档 //计算输出分布 (qj) void calcOutDistribution() { int i, j; for(j=0; j

43、robabilityDistribution() { int i, j, k; double temp = 0; for(i=0; i

44、r)=I(qj;Pji) 【实际上就是根据互信息量公式求互信息】 double getSelfInformation() { int i, j; double I=0; for(i=0; i

45、c, char *argv[]) { double probabilityCount = 0.0; //概率和 for(int k=0; k eps) { cout<<"概率和不为 1，程序异常退出 !"<

46、D 代表限定失真 double mutualInformation1, mutualInformation2, D; int i, j, flag, nCount; //初始值 mutualInformation1 = 0; mutualInformation2 = 0; D=0; flag = 0; nCount = 0; //迭代次数指示器 //init mothod //输出分布的初始化 for(i=0; i

47、始化 . 实用文档 for(i=0; i

48、OutDistribution(); //迭代算法 cout<<"误差精度： "< eps) { cout<

49、率分布 calcProbabilityDistribution(); //在上面的原来的输出分布 q 和新生成的条件概率分布 Pji 的基础上获得新的互信息 R(r, r+1) mutualInformation2 = getSelfInformation(); . 实用文档 //再计算条件概率分布 calcOutDistribution(); cout<<" 互 信息 1 ： "<

50、<

51、 } return 0; } 五、实验结果 运行实验结果如下： 六、实验总结 通过这次实验，让我们更好的掌握了率失真的求解方法，而且通过计算机解 决问题效率提高了很多， 节省了很多繁琐的步骤， 更加直观和方便的让我们了解 到相关参数变化对率失真的影响。 实验六 差错控制方法 . 实用文档 一、实验目的 1、 了解纠错编码的基本原理 2、了解几种常用编码：奇偶校验码、

52、正反码等，线性分组码、循环码、卷 积码的编解码原理 3、 重点掌握线性分组码、循环码、卷积码的编解码原理。 二、实验原理 N 个重复码是一种将输入比特重复 n 遍的编码，假设信道的错误率为 p， 接收端收到 n 个比特后进行译码，如果 n 个接收比特的“ 1”的个数多于” 0 “的个数，则译码为“ 1”反之为“ 0”，假设编码输入时等概的。 （1）计算 n=5 的信道错误率与译码的错误率的关系； （2）用 matlab 仿真得到上述的曲线。 三、实验内容 n 重复码是一种将输入比特重复 n 遍的编码，假设信道的错误率为

53、 p ，接收端收到 n 个比特后进行译码，如果 n 个接收比特的“ 1”的个数多于“ 0”的个数，则译码为“ 1”，反之为“ 0”。假设编码输入时等概的。 （1）计算 n =5 时信道错误率与译码错误率的关系； （2）用 Matlab 仿真得到上述的曲线； 实验步骤： （ 1）令 n1，n2 分别表示接收到的 n 个比特中“ 0”和“ 1”的个数，则误码率可以写成 Pb=P（n1n0|”0”)P(0) 当 n=5 时，编码时“ 1”被映射成“ 11111”；“ 0”映射成“ 00000”，信道错误率为 p，则 P(

54、 n1 n2 "1") C50 pe5 C51 (1 pe ) pe4 C52 (1 pe ) 2 pe3 P( n n 2 "0" ) C 0 p 5 C1 (1 p ) p 4 C 2 (1 p ) 2 p 3 1 5 e 5 e e 5 e e 因此 Pb pe5 5 pe4 (1 pe ) 10 pe3 (1 pe ) 2 四、实验环境 Microsoft Windows 7 Matlab 6.5 五、编码程序

55、 . 实用文档 MATLAB编码： n=5; m=0:-0.5:-3; pe=10.^m; Data d=(sign(randn(1,100000))+1)/2; s=[d;d;d;d;d]; s=reshape(s,1,5*length(d)); for k=1:length(pe) err=rand(1,length(d)*5); err=err2;

56、 error(k)=sum(abs(dd-d))/length(d); end loglog(pe,error) 六、实验结果 七、实验总结 通过本次实验，掌握了差错控制编码的实验原理与编码过程。同时通过实验解决了书 本上的例题，学会了计算信道的错误率与译码错误率的关系，能更好的理解编码、解码原 理。 . 实用文档 实验七 汉明编

57、码 一、实验目的 1、掌握线性分组码的编码原理 2、掌握汉明码编码方法 3、了解编码对误码性能的改善 二、实验原理 （n，k）线性分组码的 H 矩阵是一个 (n k ) n r n 阶矩阵，这里 r n k 是校验元的数目。显然， r 个校验元能组成 2 r 列互不相同的 r 重矢量，其中非全 零矢量有 2 r 1个。如果用这 2r 1个非全零矢量作为 H 矩阵的全部列，即令 H 矩 阵的列数 n 2 r 1，则此 H 矩阵的各列均不相同，且无全零列，由此可构造一 个纠正单个错误的（ n，k）线性分

58、组码。 同时， 2r 1是 n 所能取的最大值，因为如果 n 2r 1，那么 H 矩阵的 n 列中 必会出现相同的两列， 这样就不能满足对 H 矩阵的要求。而由于 n 2r 1是 n 所 能取的最大值，也就意味着码率 R 取得了最大值，即 k n r r r R 1 n 1 1 n n 2r 这样设计出来的码是符合我们的要求的，这样的码就是汉明码。 定义 若 H 矩阵的列是由非全零且互不相同的所有二进制 r 重矢量组

59、成，则 由此得到的线性分组码，称为 GF(2)上的（ 2r 1， 2r 1 r ）汉明码。 三、实验内容 1、写出产生汉明码的的 Matlab 程序 2、将程序在计算机上仿真实现，验证程序的正确性并完成习题。 四、实验环境 Microsoft Windows 7 Matlab 6.5 五、编码程序 本实验要求写出产生（ 3， 1）汉明码的生成矩阵，由上述可知，我们的 n=2 r 1=3, 而 k= 2r 1 r =1, 由此可得出 r=2. 当 r=2 时，有 3 个非全零的二重矢量： （01），（10），（1

60、1） 构成矩阵： . 实用文档 0 1 1 H 1 0 1 由此得到一个能纠正单个错的（ 3，1）汉明码。若码字传输中左边第一位出 错，则相应的伴随式 s (01) 就是 H 矩阵的第一列，也正好是“ 1”的二进制表示。 同理可知，无论哪一位出错， 它对应的伴随式就是该位的二进制表示， 故译码十 分方便，特别适用于计算机内部运算和记忆系统中的纠错。 如果要得到系统码形式的 H 矩阵，只需对上述矩阵进行初等变换 交换列即 可。 H 1 1 0 1 0 1 相

61、应地，生成矩阵 G 为： G 1 1 1 根据 C G 生成码字。 由此构成的（ 3， 1）汉明码如表 3-1 所示。 表 3-1 （3， 1）系统码 信息组 码字 0 000 1 111 MATLAB程序： %function f=humm_enc(a); G=[1 1 1];%(3,1)汉明码的生成矩阵 t=input( 输‘入 0或1:');%t=0则产生 (3,1)汉明码， t=1则对输入序列进行编码 if

62、 t==1 a=input('输入信息元序列 :'); c=mod(a*G,2);%编码的码字 c disp('编码后序列 :');%显示“编码后序列 :”字样 disp(c); %显示编码的码字 else disp('(3,1)汉明系统码为： ');%显示“ (3,1)汉明系统码为：”字样 for i=0:1; %循环变量的取值为 0,1 a=dec2bin(i,1); %系统自动给出两信息元从 0到 1 c=mod(a*G,2) ; %输出编码的码字 disp(a); %显示系统自动给出两信息元从 0到 1 d

63、isp('对应码字为： '); %显示“对应码字为 :”字样 . 实用文档 disp(c); %显示编出的码字 end %结束 for循环 end %结束整个程序 六、实验结果 ① Untitled3 输入 0或1:1 输入信息元序列 :0 编码后序列 : 0 0 0 ② Untitled3 输入 0或1:0 (3,1) 汉明系统码为： 0 对应码字为： 0 0 0 汉明系统码为： 1 对应码字为： 1 1 1 七、实验总结 通过这次实验我对汉明编码有了进一步的认识，对生成矩阵的求法更熟悉 了，同时对仿真模块的应用也有了深刻的认识， 认真检验了各参数对编码的影响， 此实验将为以后相类似实验的基础，以后相类似实验将更容易更快捷的被做出， 可达到了学以致用的效果。 .