Huffman编码的Matlab实现_信息论基础huffman编码matlab上机作业-程序员宅基地

Huffman编码的Matlab实现

实现代码

信息论与编码大作业

function [ Huff_chart, Ave_length, yita] = huff_code( c )

%统计字符串中字符出现频率
    t=unique(c);
    for i=1:length(t)
        Num(i)=length(strfind(c,t(i)));
    end
    Num=Num/sum(Num);

%构建Huffman编码表
    Huff_chart=cell(3,length(t));
    for i=1:length(t)
        Huff_chart(1,i)={t(i)};
        Huff_chart(2,i)={Num(i)};
        Huff_chart(3,i)={0};
    end

%生成待编码矩阵，每列最后一个元素为新求和得到元素的位置
    n=length(t);   
    A=sort(Num,'descend');%按降序排列
    T=A;
    A=A/sum(A);
    B=zeros(n,n);
    %生成第一列
    for i=1:n
        B(i,1)=T(i);
    end
    r=B(i,1)+B(i-1,1);
    T(n-1)=r;T(n)=0;
    T=sort(T,'descend');
    t=n-1;
    %生成编码表的其他各列
    for j=2:n
        for i=1:t
            B(i,j)=T(i);
        end
        if t>1
            K=find(T==r);
            B(n,j)=K(end);
            r=(B(t-1,j)+B(t,j));
            T(t-1)=r;
            T(t)=0;
            T=sort(T,'descend');
            t=t-1;
        else
            B(n,j)=1;
        end
    end
    
%对编码矩阵进行Huffman编码编码格式为：
    %2--->0
    %1--->1
    
    Huff_code=zeros(1,n);
    Huff_code(1,1)=2;
    Huff_code(1,2)=1;
     for i=n-2:-1:1
         t1=Huff_code(B(n,i+1))*10+1;
         t0=Huff_code(B(n,i+1))*10+2;
         Huff_code(B(n,i+1))=[];
         Huff_code(n-i+1)=t1;
         Huff_code(n-i)=t0;
     end
     
%将Huffman编码对应输出到Huffman编码表格中  
    A1=A;
    for i=1:n
        t=find(abs(A-Num(i))<10.^-3);
        t=t(1);
        A(t)=0;
        Huff_chart(3,i)={strrep(num2str(Huff_code(1,t)),'2','0')};
    end
    
%相关计算   
    Ave_length=zeros(1,n);
    for i=1:n
        Ave_length(i)=length(num2str(Huff_code(i)));
    end
    Ave_length=sum(Ave_length.*A1);
    H1=log2(Num);
    H=-Num*(H1');
    yita=H/Ave_length;
    
end

本文链接：https://blog.csdn.net/qq_43633720/article/details/104075508

原作者删帖不实内容删帖广告或垃圾文章投诉

智能推荐

linux devkmem 源码,linux dev/mem dev/kmem实现访问物理/虚拟内存-程序员宅基地

文章浏览阅读451次。dev/mem: 物理内存的全镜像。可以用来访问物理内存。/dev/kmem: kernel看到的虚拟内存的全镜像。可以用来访问kernel的内容。调试嵌入式Linux内核时，可能需要查看某个内核变量的值。/dev/kmem正好提供了访问内核虚拟内存的途径。现在的内核大都默认禁用了/dev/kmem，打开的方法是在 make menuconfig中选中 device drivers --> ..._dev/mem 源码实现

vxe-table 小众但功能齐全的vue表格组件-程序员宅基地

文章浏览阅读7.1k次，点赞2次，收藏19次。vxe-table，一个小众但功能齐全并支持excel操作的vue表格组件_vxe-table

（开发）bable - es6转码-程序员宅基地

文章浏览阅读62次。参考：http://www.ruanyifeng.com/blog/2016/01/babel.htmlBabelBabel是一个广泛使用的转码器，可以将ES6代码转为ES5代码，从而在现有环境执行// 转码前input.map(item => item + 1);// 转码后input.map(function (item) { return item..._让开发环境支持bable

FPGA 视频处理 FIFO 的典型应用_fpga 频分复用视频-程序员宅基地

文章浏览阅读2.8k次，点赞6次，收藏29次。摘要：FPGA视频处理FIFO的典型应用，视频输入FIFO的作用，视频输出FIFO的作用，视频数据跨时钟域FIFO，视频缩放FIFO的作用_fpga 频分复用视频

R语言：设置工作路径为当前文件存储路径_r语言设置工作目录到目标文件夹-程序员宅基地

文章浏览阅读575次。【代码】R语言：设置工作路径为当前文件存储路径。_r语言设置工作目录到目标文件夹

background 线性渐变-程序员宅基地

文章浏览阅读452次。格式：background: linear-gradient(direction, color-stop1, color-stop2, ...);<linear-gradient> = linear-gradient([ [ <angle> | to <side-or-corner>] ,]? &l..._background线性渐变

随便推点

【蓝桥杯省赛真题39】python输出最大的数中小学青少年组蓝桥杯比赛算法思维python编程省赛真题解析-程序员宅基地

文章浏览阅读1k次，点赞26次，收藏8次。第十三届蓝桥杯青少年组python编程省赛真题一、题目要求（注：input（）输入函数的括号中不允许添加任何信息）1、编程实现给定一个正整数N，输出正整数N中各数位最大的那个数字。例如:N=132，则输出3。2、输入输出输入描述：只有一行，输入一个正整数N输出描述：只有一行，输出正整数N中各数位最大的那个数字输入样例：

网络协议的三要素-程序员宅基地

文章浏览阅读2.2k次。一个网络协议主要由以下三个要素组成：1.语法数据与控制信息的结构或格式，包括数据的组织方式、编码方式、信号电平的表示方式等。2.语义即需要发出何种控制信息，完成何种动作，以及做出何种应答，以实现数据交换的协调和差错处理。3.时序即事件实现顺序的详细说明，以实现速率匹配和排序。不完整理解：语法表示长什么样，语义表示能干什么，时序表示排序。转载于:https://blog.51cto.com/98..._网络协议三要素csdn

The Log: What every software engineer should know about real-time data's unifying abstraction-程序员宅基地

文章浏览阅读153次。主要的思想，将所有的系统都可以看作两部分，真正的数据log系统和各种各样的query engine所有的一致性由log系统来保证，其他各种query engine不需要考虑一致性，安全性，只需要不停的从log系统来同步数据，如果数据丢失或crash可以从log系统replay来恢复可以看出kafka系统在linkedin中的重要地位，不光是d..._the log: what every software engineer should know about real-time data's uni