matlab 哈夫曼编码

最强废铁h

1. clc;
   
2. clear;
   
3. n=10;
   
4. % c=[];
   
5. c=randsample(100,10);    %随机生成10个数字
   
6. % c=rand(1,n);%********随机生成数据********************
   
7. c=c/sum(c);
   
8. % c=[0.1546 0.0948 0.1361 0.0186 0.1258 0.0474 0.1897 0.0351 0.1423 0.0557];
   
9. p=c;
   
10. disp(p);
    
11. w=[;];%***************位置矩阵，记录数据的变化过程*****
    
12. for i=1:n-1
    
13.     [p,l]=sort(p);
    
14.     w(i,:)=[l(1:n-i+1),zeros(1,i-1)]; %存放每一次排序的顺序
    
15.     p=[p(1)+p(2),p(3:n),1];
    
16. end
    
17. a={;};%**********字符数组，生成霍夫曼码****************
    
18. for i=1:2*(n-1)
    
19.     for j=1:n
    
20.         
    
21.         a{i,j}=' ';
    
22.     end
    
23. end%*****************初始化*********************
    
24. a{2*n-2,1}='0';
    
25. a{2*n-2,2}='1';
    
26. for i=1:n-1 % 1-9
    
27.     for j=1:i+1  %
    
28.         a{2*(n-i)-1,j}=a{2*(n-i),find(w(n-i,:)==j)};%恢复未排序前  奇数行
    
29.     end
    
30.     if((n-i-1)<=0)
    
31.         break;
    
32.     end
    
33.     a{2*(n-i-1),1}=strcat(a{2*(n-i)-1,1},'0');   %每次都是最小两个概率合并，这是恢复到合并前
    
34.     a{2*(n-i-1),2}=strcat(a{2*(n-i)-1,1},'1');   %前两个每次加一个码字，分别为0和1
    
35.     for j=3:i+2
    
36.         a{2*(n-i-1),j}=a{2*(n-i)-1,j-1};   %除取前两个的，后面和上一行结果相同
    
37.     end
    
38. end  
    
39. h=-c.*log2(c);
    
40. hx=sum(h);
    
41. for i=1:n
    
42.     len(i)=length(a{1,i});
    
43. end
    
44. plen=len.*c';
    
45. averagelen=sum(plen);
    
46. fprintf('霍夫曼编码为:\n');
    
47. % disp(a);
    
48. for i=1:n
    
49.     disp(a(1,i));
    
50. end
    
51. fprintf('平均码长为:\n');
    
52. disp(averagelen);
    
53. fprintf('信源熵为:\n');
    
54. disp(hx);
    
55. fprintf('编码效率为:\n');
    
56. disp(hx/averagelen);

复制代码

1. 错误使用 horzcat
   
2. 串联的矩阵的维度不一致。
   
3. 
   
4. 出错 huffman_1 (line 15)
   
5.     p=[p(1)+p(2),p(3:n),1];

复制代码

烟雨任平生

请问解决了吗？

最强废铁h

烟雨任平生 发表于 2022-5-13 16:47
请问解决了吗？

我换代码了，