(Matlab源码)Matlab实现算术编码(Arithmetic coding)超级详解(每一段代码都可以看懂)

1、代码功能

  • 输入:一个字符串
  • 输出:
    •   codeword(码值)
    •  codeword所占的位数

2、代码框图

3、代码超详解

  •  统计字符串中的字符种类,调用函数,放入数组b中。
b=unique(str1);
  • 统计每种字符的个数,放入数组a中。
for i=1:length(b)
    a(i)=length(strfind(str1,b(i)));
end
  •  排序,形成对应关系的数组。排序后,若数组b中的第一个元素为字符‘a’, 则数组a中的第一个元素为字符‘a’,出现的次数
for i=1:length(a)-1
    for j=i+1:length(a)
        if(a(i)<a(j))
            t1=a(i);
            a(i)=a(j);
            a(j)=t1;
            t2=b(i);
            b(i)=b(j);
            b(j)=t2;
        end
    end
end
  • 统计每种字符出现的频率,放入数组frequency中。
for i=1:length(a)
    frequency(i) = a(i)/sum(a);
end

 根据频率分配概率范围,范围高值放入数组high中,范围低值放入数组low中。

low(1)=0;
for i = 1:length(frequency)
    if i==1
        high(i)=frequency(i);
    else
        low(i)=frequency(i-1)+low(i-1);
        high(i)=low(i)+frequency(i);
    end
end

编码

%编码
low_value=0; %初始化
high_value=1; %初始化
range=high_value-low_value; %初始化
for i = 1:length(str1) %遍历字符串中的每一个元素,使用for循环语句。
    
     d=strfind(b,str1(i)); %找到字符串中的字符,在数组d中的角标
     
     %算术编码
     high_value = low_value+ range*high(d);
     low_value = low_value+range*low(d);
     range=high_value-low_value;
     
     %编码后数据处理
     if i==length(str1)
         ave=(high_value+low_value)/2;
         ave_str=num2str(ave,length(str1)); %把数字转成字符串
         low_value_str=num2str(low_value,length(str1)); %把数字转成字符串
         high_value_str=num2str(high_value,length(str1)); %把数字转成字符串
         
         %当ave_str与low_value_str对应位的元素不一样时
         %截去ave_str不一样元素后面的字符
         for ii= 3: length(low_value_str)
             if low_value_str(ii) ~= high_value_str(ii)
                  ave_str((ii+1):end)=[];
                  break;
             end
             ii=ii+1;           
         end
     end
end

%把ave_str值赋给codeword
codeword=ave_str
%计算codeword占用的位数
bits=-log2(range)

 

  • 整个代码地址

       Matlab实现算术编码

4、代码缺点

  • 这个代码可以实现短字符串的编码,对于长字符串的编码,会出现underflow。如何解决underflow, 这篇博客有专题研究

       算术编码(arithmetic coding)的underflow问题

  • 6
    点赞
  • 85
    收藏
    觉得还不错? 一键收藏
  • 5
    评论
好的,以下是一个简单的算术编码和译码的 MATLAB 代码实现,供参考: 算术编码: ```matlab function encoded = arithmetic_encode(msg, prob) % 计算每个符号的累积概率 cumprob = cumsum(prob); % 初始化区间 low = 0; high = 1; % 对于每个符号,更新区间 for i=1:length(msg) symbol = msg(i); range = high - low; high = low + range * cumprob(symbol+1); low = low + range * cumprob(symbol); end % 返回编码后的值 encoded = (low + high) / 2; end ``` 译码: ```matlab function decoded = arithmetic_decode(encoded, prob, msg_len) % 计算每个符号的累积概率 cumprob = cumsum(prob); % 初始化区间 low = 0; high = 1; % 初始化译码结果 decoded = zeros(1, msg_len); % 对于每个符号,更新区间并译码 for i=1:msg_len range = high - low; cumulative = (encoded - low) / range; % 找到对应的符号 for j=1:length(cumprob) if cumulative < cumprob(j) decoded(i) = j - 1; high = low + range * cumprob(j); low = low + range * cumprob(j-1); break; end end end end ``` 其中,`msg` 为待编码的消息,`prob` 为每个符号的概率,`encoded` 为编码后的值,`msg_len` 为消息长度。具体使用方式可以参考以下示例: ```matlab % 测试数据 msg = [2 1 3 0 2 3 1 2]; prob = [0.25 0.25 0.25 0.25]; % 编码 encoded = arithmetic_encode(msg, prob); % 输出编码结果 disp(['Encoded value: ' num2str(encoded)]); % 译码 decoded = arithmetic_decode(encoded, prob, length(msg)); % 输出译码结果 disp(['Decoded message: ' num2str(decoded)]); ```
评论 5
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值