前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >香农编码c++实现_香农费诺编码例子

香农编码c++实现_香农费诺编码例子

作者头像
全栈程序员站长
发布2022-09-30 20:31:06
4460
发布2022-09-30 20:31:06
举报
文章被收录于专栏:全栈程序员必看

大家好,又见面了,我是你们的朋友全栈君。

实验三香农编码

信息论与编码实验报告

院系: 哈尔滨理工大学荣成校区

专业: 电子信息工程

学号:

姓名:

日期: 2015年6月16日

香农编码

信息论与编码第三次实验报告

一、实验目的和任务?

1、?理解信源编码的意义;?

2、?熟悉?MATLAB程序设计;??

3、?掌握香农编码的方法及计算机实现;??

4、?对给定信源进行香农编码,并计算编码效率;?

二、实验原理介绍?

给定某个信源符号的概率分布,通过以下的步骤进行香农编码

?1、信源符号按概率从大到小排列;?

???

2、确定满足下列不等式的整数码长为?

????

3、为了编成唯一可译码,计算第i个消息的累加概率:?

?

?4、将累加概率变换成二进制数;?

取二进制数的小数点后位即为该消息符号的二进制码字。

三、实验设备

计算机

Matlab2014

四、编码程序

N=input(‘N=’); %输入信源符号的个数

s=0; %对行的概率进行初始化

l=0; %对列的概率初始化

H=0; %对信息熵初始化

for i=1:N

p(i)=input(‘p=’);%输入信源符号概率分布矢量

p(i)<1; %输入概率小于1

s=s+p(i); %行概率求和

H=H+(-p(i)*log2(p(i)));%计算信源信息熵

end

if abs(s-1)>0 %若行概率和>1则报错

error(‘不符合概率分布’)

end

for i=1:N-1 %for循环按概率分布大小对信源排序

for j=i+1:N

if p(i)

m=p(j);

p(j)=p(i);

p(i)=m;

end

end

end

for i=1:N

a=-log2(p(i)); %确定满足不等式的整数a,并领a为第i个码字的长度

if mod(a,1)==0 %a对1取余后赋值给a

w=a;

else

w=fix(a+1); %计算各信源符号的码长

end

l=l+p(i)*w; %计算平均码长

end

r=l; %香农编码的信息率(这里L=1,m=2)

n=H/r; %计算编码效率

P(1)=0; %初始1列为0;

for i=2:N

P(i)=0; %初始1列为0;累加

for j=1:i-1

P(i)=P(i)+p(j); %从2列开始累加

end

end

%计算累加概率

for i=1:N

for j=1:w

W(i,j)=fix(P(i)*2); %对每次的P(i)取码字,长度与a相同

P(i)=P(i)*2-fix(P(i)*2); %将累加概率转化为L(i)位二进制码字

end

end

W %显示码字

r %显示平均码长

n %显示编码效率

H %显示信息熵H

实验结果

>> xiangnong

N=6

p=0.25

p=0.25

p=0.20

p=0.15

p=0.10

p=0.05

W =

0 0 0 0 0

0 1 0 0 0

1 0 0 0 0

1 0 1 1 0

1 1 0 1 1

1 1 1 1 0

r =

2.7000

n =

0.8975

H =

2.4232

六、实验心得

此次实验让我明白了做程序需要先进行需求分析,然后画出流程图,根据流程图把每部分模块完成,然后再进行调试和编辑,提高了我做

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员栈长,转载请注明出处:https://javaforall.cn/195272.html原文链接:https://javaforall.cn

本文参与 腾讯云自媒体同步曝光计划,分享自作者个人站点/博客。
原始发表:2022年9月10日 ,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 作者个人站点/博客 前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体同步曝光计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档