二元HammingCode的Matlab实现（转贴）

suffer

1. function [index,H,s]=jghostyhammingcode(f,r)
   
2. %---------HammingCode 译码---------------------------
   
3. %---------HammingCode By jghosty 2006-11-29----------
   
4. %f为输入二元码，码长n=2^r-1,index为标志符，H是校验矩阵，s为最终输出码
   
5. lf=length(f);
   
6. hamlength=2^r-1;
   
7. %判断码长是否符合
   
8. if lf<hamlength
   
9.     index='输入错误！'
   
10.     return;
    
11. end
    
12. for i=1:15
    
13.     x(i)=i;
    
14. end
    
15. h=ones(1,hamlength);
    
16. for i=1:hamlength
    
17.     h(i)=i;
    
18. end
    
19. %生成校验矩阵
    
20. H1=dec2bin(h);
    
21. H=H1';
    
22. S=mod(f*(H1),2);
    
23. %判断码字是否正确
    
24. if mod(S,2)==0
    
25.     index='无错误';
    
26.     s=f;
    
27.     return;
    
28. else
    
29.     index='译码出错,纠错处理';
    
30.     K=ones(r,1);
    
31.     for j=1:hamlength
    
32. for l=1:r
    
33.         K(l)=H1(j,l);
    
34. end
    
35.      G=mod(S*K,2);
    
36.      if G~=0;
    
37.          break;
    
38.         end
    
39.     end
    
40. s=mod((f+[zeros(1,j-1),1,zeros(1,hamlength-j)]),2);
    
41. end

复制代码