信息论-----哈夫曼编码实验

clear
p = [0.2 0.4 0.1 0.15 0.15];
disp(p);
len = size(p,2);
rank = zeros(len,len+1);%定义一个矩阵用来存储每次的概率排序及合并的结果
rank(:,1:2) = [p' (1:len)'];%定义一个矩阵用来存储每次的概率排序及合并的结果，第1列为概率，从第2列开始为此概率对应的符号标号  
cw = ones(len,len-1);%定义一个矩阵用于记录码字中0，1分配过程
cw = 2*cw;
for i=1:len-1  %编码次数比符号个数少1
    %%%%%%%%（1）概率排序%%%%%%%%%%%%%
    rank =
sortrows(rank,-1);%概率降序排序
    %%%%%%%%%%（2）给信源符号分配0，1码元%%%%%%%%%%%%
    R_L_0 =
rank(len-i,2:len);%确定分配0码元的符号标号所在行
    R_L_1 =
rank(len+1-i,2:len);%确定分配1码元的符号标号所在行
    num_0 = R_L_0(R_L_0~=0);%num表示本次分配0的概率对应哪些信源符号。rank中的0元素不是符号标号，去除0。num_0存储的数据为此次分配0的符号标号
    num_1 =
R_L_1(R_L_1~=0);
cw(num_0,len-i) = 0 ;%在cw矩阵中给相应的符号分配码元0
cw(num_1,len-i) = 1;
    %%%%%%%%%%%（3）将概率最小的两个概率相加%%%%%%%%%%%%%%%%%%
    rank(len-i,1)
= rank(len+1-i,1)+rank(len-i,1);%概率相加
rank(len+1-i,:) = 0;%两个概率合并后，其中下面的概率所在行要置零，防止下一轮概率排序时产生干扰。
rank(len-i,2:length([num_1 num_0])+1) = [num_1 num_0];%将合并的概率（即前面分配0和1码元的符号）所对应的符号标号记录在rank矩阵中
end
cw2 = cell(len,1);%定义一个len*1的元胞数组，它的每个单元可以储存不同的数据类型。
for i=1:len
    tmp =
cw(i,:);
    cw2{i} =
num2str(tmp(tmp~=2));%将cw中的0和1取出来（2不是分配的码元），并转换为字符串类型
disp(cw2{i});
end