首页> R语言教程

用R语言实现霍夫曼编码的示例代码

withpy 2021-05-25

简介这篇文章主要介绍了用R语言实现霍夫曼编码的示例代码,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧

可读性极低,而且其实也没必要用R语言写,图个乐罢了 

p=c(0.4,0.2,0.2,0.1,0.1)###输入形如c(0.4,0.2,0.2,0.1,0.1)的概率向量,即每个待编码消息的发生概率
p1=p###将概率向量另存,最后计算编码效率要用
mazijuzhen=matrix(,nrow=length(p),ncol=length(p)-1)###码字矩阵:第i行对应向量p的第i个分量所对应的那个待编码消息的编码后的码字
group=matrix(c(1:length(p),rep(NA,length(p)*(length(p)-1))),nrow=length(p),ncol=length(p))###初始分组:每一行代表一组,每个行向量的所有分量代表此组的所有元素,初始时,有多少个待编码消息就分多少个组,每组只有一个待编码消息,以整数i代表向量p的第i个分量所对应的那个待编码消息
i=1###开始编码
for(i in 1:(length(p)-1))
{
 orderp=order(p,decreasing = FALSE)###orderp的分量依次是:p的最小分量的下标,p的第二小分量的下标。。。
 mazijuzhen[group[orderp[1],],i]=0###给概率最小的两个消息组编上0和1
 mazijuzhen[group[orderp[2],],i]=1
 group[min(c(orderp[1],orderp[2])),]=c(na.omit(group[min(c(orderp[1],orderp[2])),]),na.omit(group[max(c(orderp[1],orderp[2])),]),rep(NA,length(p)-length(c(na.omit(group[min(c(orderp[1],orderp[2])),]),na.omit(group[max(c(orderp[1],orderp[2])),])))))###把此次迭代的两个消息组中组编号较大的分到组编号较小的组里去。
 group[max(c(orderp[1],orderp[2])),]=NA###删除组编号较大的组
 p[min(c(orderp[1],orderp[2]))]=p[orderp[1]]+p[orderp[2]]###计算本次迭代得到的新的消息组的发生概率
 p[max(c(orderp[1],orderp[2]))]=NA###由于组编号较大的组被删除,所以相应删除它所对应的概率
 print("当前迭代次数")###本次迭代的结果总结
 print(i)
 print("概率向量")
 print(p)
 print("分组矩阵")
 print(group)
 print("码字矩阵")
 print(mazijuzhen)
}
i=1###由霍夫曼编码的特性,将所有编码倒转得到最终编码
for (i in 1:length(p)) 
{
 mazijuzhen[i,]=rev(mazijuzhen[i,]) 
}
i=1###构建码长向量
machang=c()
for (i in 1:length(p))
{
 machang=c(machang,length(na.omit(mazijuzhen[i,])))
}
xiaolv=-p1%*%log(p1,2)/mean(machang)###计算编码效率
print("最终的码字矩阵和编码效率")
mazijuzhen
xiaolv

到此这篇关于用R语言实现霍夫曼编码的示例代码的文章就介绍到这了,更多相关R语言 霍夫曼编码内容请搜索本站以前的文章或继续浏览下面的相关文章希望大家以后多多支持本站!

相关文章

  • Python实现Kerberos用户的增删改查操作

    这篇文章主要介绍了Python实现Kerberos用户的增删改查操作,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧

  • 用C语言实现五子棋游戏

    这篇文章主要为大家详细介绍了用C语言实现五子棋游戏,文中示例代码介绍的非常详细,具有一定的参考价值,感兴趣的小伙伴们可以参考一下

  • php使用pecl方式安装扩展操作示例

    这篇文章主要介绍了php使用pecl方式安装扩展操作,结合实例形式分析了pecl的安装、以及基于pecl的扩展查找、安装,并给出了swoole与xdebug扩展安装的操作步骤,需要的朋友可以参考下

  • 详解java为什么不允许类多重继承却允许接口多重继承

    这篇文章主要介绍了java为什么不允许类多重继承却允许接口多重继承,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧

  • python argparse命令行参数解析(推荐)

    Python argparse模块是解析命令行参数的首选方法。解析命令行参数是一个非常常见的任务,Python脚本根据传递的值来执行和操作

  • AngularJS基于http请求实现下载php生成的excel文件功能示例

    这篇文章主要介绍了AngularJS基于http请求实现下载php生成的excel文件功能,结合实例形式分析了AngularJS http请求及文件下载等相关操作技巧,需要的朋友可以参考下