| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 2614 人关注过本帖
标题:matlab实现huffman编码
只看楼主 加入收藏
fangvsyu
Rank: 1
等 级:新手上路
帖 子:2
专家分:0
注 册:2010-5-7
收藏
 问题点数:0 回复次数:3 
matlab实现huffman编码
如何用matlab实现图象的huffman编解码,最好有详细的解释.谢谢!
搜索更多相关主题的帖子: matlab huffman 编码 
2010-05-07 18:15
ahyshong
Rank: 2
等 级:论坛游民
帖 子:28
专家分:86
注 册:2009-2-14
收藏
得分:0 
function A=f(a,b1)
    A=zeros(1,length(b1))
for i=1:length(b1)
    A(i)=length(strfind(a,b1(i)));
end

function END=huffman(A)
A=A/sum(A);
A=fliplr(sort(A));
 T=A;
[m,n]=size(A);
B=zeros(n,n-1);
 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=fliplr(sort(T));
t=n-1;
for j=2:n-1
    for i=1:t
        B(i,j)=T(i);
    end
        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=fliplr(sort(T));
        t=t-1;
end
B;
END1=sym('[0,1]');
END=END1;
t=3;
d=1;
for j=n-2:-1:1
    for i=1:t-2
        if i>1&&B(i,j)==B(i-1,j)
            d=d+1;
        else
            d=1;
        end
        B(B(n,j+1),j+1)=-1;
        temp=B(:,j+1);
 
        x=find(temp==B(i,j));
        END(i)=END1(x(d));
    end
    y=B(n,j+1);
    END(t-1)=[char(END1(y)),'0'];
    END(t)=[char(END1(y)),'1'];
    t=t+1;
    END1=END;
end
这个是我以前用过的程序,你自己看着修改修改吧。
不明白的命令在help里面找。
2010-05-24 09:35
ahyshong
Rank: 2
等 级:论坛游民
帖 子:28
专家分:86
注 册:2009-2-14
收藏
得分:0 
A为需要编码的那个序列
2010-05-24 09:36
AntsiLynn
Rank: 1
等 级:新手上路
帖 子:1
专家分:0
注 册:2017-5-29
收藏
得分:0 
回复 2楼 ahyshong
在Matlab上跑了一遍 好像不对?在n>10后会出现编码效率大于1 而且能看出来编码出来的不是唯一码了
程序如下
N=10;%N位的信源
SignalSourcewithoutNormilization=unifrnd(0,1,1,N);%产生随机数
sum=0;
for n=1:N
    sum=sum+SignalSourcewithoutNormilization(n);
end
SignalSource=SignalSourcewithoutNormilization./sum;%归一化
SignalSource=sort(SignalSource,'descend')
%排序
SortingMatrix=zeros(N,N);
for i=1:N
    P=sort(SignalSource,'descend');%降序排列各概率
    SortingMatrix(i,1)=P(i);%将信源各概率输入排列矩阵第一列
end
rear=SortingMatrix(i,1)+SortingMatrix(i-1,1);%合并概率最小两项
P(N-1)=rear;P(N)=0;
P=sort(P,'descend');
t=N-1;

for j=2:n%生成编码表的其他各列
        for i=1:t
            SortingMatrix(i,j)=P(i);
        end
        if t>1
            K=find(P==rear);
            SortingMatrix(N,j)=K(end);%每列最后一个元素用于记录合并发生位置
            rear=(SortingMatrix(t-1,j)+SortingMatrix(t,j));%合并最后两项(即概率最小两项)
            P(t-1)=rear;
            P(t)=0;
            P=sort(P,'descend');
            t=t-1;
        else
            SortingMatrix(n,j)=1;
        end
    end

    %生成Huffman码字矩阵和排序后元素的码表
    m=3;s1=sym('[2,1]');s2=s1;%码表,用1表示0,2表示1,因为0会因为被省略导致无法记录
    %最后一列和倒数第二列的情况可以直接推出故从倒数第三列开始建表
    for i=N-2:-1:1
        p=SortingMatrix(N,i+1);%p是每列保存上一列最后两个数的合并项的位置
        if p==1
            s2(1:m-2)=s1(2:m-1);
        elseif p==m
            s2(1:m-2)=s1(1:m-2);
        elseif p==2
            s2(1)=s1(1);
            s2(2:m-2)=s1(3:m-1);
        else
            s2(1:p-1)=s1(1:p-1);
            s2(p+1:m-2)=s1(p+1:m-2);
        end
        s2(m-1)=[char(s1(p)),'2'];
        s2(m)=[char(s1(p)),'1'];
        m=m+1;
        s1=s2;
    end
    L=zeros(1,N);
    for i=1:N
        [~,rear]=size(char(s2(i)));
        L(i)=rear;
    end

    %将码表转换成0和1输出
    array=zeros(1,N);
    array(1:n)=s2(1:N);
    String=[];
    for i=1:n
        s=num2str(array(i));
        s=strrep(s,'1','0');
        s=strrep(s,'2','1');
        if i==1
            String=s;
        else
            String=[String,' ',s];
        end
    end
   
   String=regexp(String,' ','split')%编码后的输出码表
   
   AverageLength=0;
   for n=1:N
       AverageLength=AverageLength+L(n)*SignalSource(n);
   end
    AverageLength%平均码长
   
    H=0;
    for n=1:N
        H=H-SignalSource(n)*log2(SignalSource(n));
    end
    H
   
    Efficiency=H/AverageLength%编码效率

出来的结果:
>> Huffman

SignalSource =

    0.2088    0.1922    0.1591    0.1555    0.1081    0.0495    0.0409    0.0382    0.0247    0.0230


String =

    '01'    '00'    '111'    '110'    '101'    '1111'    '11111'    '11111'    '111111'    '111110'


AverageLength =

    2.9497


H =

    2.9558


Efficiency =

    1.0021
2017-05-29 15:16
快速回复:matlab实现huffman编码
数据加载中...
 
   



关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.037714 second(s), 7 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved