hamming.c 5.0 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300
  1. #include "hamming.h"
  2. #include "math.h"
  3. #include <stdbool.h>
  4. #include <string.h>
  5. #include "includes.h"
  6. //编码
  7. unsigned char HMBM_74(unsigned char a)
  8. {
  9. unsigned char b =0;
  10. b = a & 0x0F;
  11. b = b <<0x03;
  12. if(b&0x08){b=b^0x03;}
  13. if(b&0x10){b=b^0x05;}
  14. if(b&0x20){b=b^0x06;}
  15. if(b&0x40){b=b^0x07;}
  16. return(b & 0x7F);
  17. }
  18. //解码
  19. unsigned char HMYM_74(unsigned char a, uint8_t *err)
  20. {
  21. unsigned char b=0,c=0;
  22. b = a;
  23. c = 0;
  24. if(b&0x40) {c = c^0x07 ;}
  25. if(b&0x20) {c = c^0x06 ;}
  26. if(b&0x10) {c = c^0x05 ;}
  27. if(b&0x08) {c = c^0x03 ;}
  28. if(b&0x04) {c = c^0x04 ;}
  29. if(b&0x02) {c = c^0x02 ;}
  30. if(b&0x01) {c = c^0x01 ;}
  31. // printf("c = %d\n",c);
  32. *err = c;
  33. switch(c){
  34. case 0 :b = b >> 3;break ;
  35. case 1 :b = b >> 3;break ;
  36. case 2 :b = b >> 3;break ;
  37. case 3 :b = (b^0x08) >> 3;break ;
  38. case 4 :b = b >> 3;break ;
  39. case 5 :b = (b^0x10) >> 3;break ;
  40. case 6 :b = (b^0x20) >> 3;break ;
  41. case 7 :b = (b^0x40) >> 3;break ;
  42. }
  43. return (b & 0x0F);
  44. }
  45. //对u32数据编码
  46. //data:输入的待编码数据
  47. //code:编码后生成的数据(1字节的原始数据生成2字节的编码数据)
  48. void HM_encode_u32(uint32_t data, uint32_t *code)
  49. {
  50. int i,j;
  51. uint32_t temp = 0;
  52. uint8_t data_char[8] = {0};
  53. uint8_t data_encode[8] = {0};
  54. memset(code,0,2);
  55. // 原始数据拆分:
  56. for(i = 0; i < 8; i++)
  57. {
  58. data_char[i] = (data>>(4*(7-i))) & 0x0f;
  59. }
  60. //编码
  61. for(i = 0; i < 8; i++)
  62. {
  63. data_encode[i] = HMBM_74(data_char[i]);
  64. }
  65. //编码后数据合成
  66. for(i = 0; i < 2; i++)
  67. {
  68. for(j = 0; j < 4; j++)
  69. {
  70. temp = data_encode[i*4+j];
  71. code[i] |= temp<<(3-j)*8;
  72. }
  73. }
  74. }
  75. //对u32数据解码
  76. uint32_t HM_decode_u32(uint32_t *data)
  77. {
  78. int i,j;
  79. uint32_t temp = 0;
  80. uint32_t data_de=0;
  81. uint8_t data_temp[8] = {0};
  82. uint8_t data_decode[8] = {0};
  83. uint8_t err;
  84. // 原始数据拆分:
  85. for(i = 0; i < 2; i++)
  86. {
  87. for(j = 0; j < 4; j++)
  88. {
  89. data_temp[i*4+j] = (data[i]>>(3-j)*8) & 0xff;
  90. }
  91. }
  92. //解码
  93. for(i = 0; i < 8; i++)
  94. {
  95. data_decode[i] = HMYM_74(data_temp[i],&err);
  96. }
  97. //解码后数据合成
  98. for(i = 0; i < 8; i++)
  99. {
  100. temp = data_decode[i];
  101. data_de |= temp<<4*(7-i);
  102. }
  103. return data_de;
  104. }
  105. //汉明码编码(一个字节)
  106. void HM_encode_char(uint8_t data, uint8_t *code)
  107. {
  108. int i;
  109. uint8_t data_char[2] = {0};
  110. memset(code,0,2);
  111. // 原始数据拆分:
  112. data_char[0] = data & 0x0f;
  113. data_char[1] = (data>>4) & 0x0f;
  114. //编码
  115. for(i = 0; i < 2; i++)
  116. {
  117. code[i] = HMBM_74(data_char[i]);
  118. }
  119. }
  120. //汉明码解码(一个字节)
  121. uint8_t HM_decode_char(uint8_t *code, uint8_t *c_err)
  122. {
  123. int i;
  124. uint8_t data = 0;
  125. uint8_t err;
  126. uint8_t data_decode[2] = {0};
  127. *c_err = 0;
  128. //解码
  129. for(i = 0; i < 2; i++)
  130. {
  131. data_decode[i] = HMYM_74(code[i],&err);
  132. if(err)
  133. {
  134. *c_err = 1;
  135. }
  136. }
  137. //解码后数据合成
  138. data |= data_decode[0];
  139. data |= data_decode[1]<<4;
  140. return data;
  141. }
  142. //汉明码编码(多字节)
  143. //返回编码后的字节数
  144. void HM_encode(uint8_t *data, uint8_t *code, uint16_t size)
  145. {
  146. int i;
  147. memset(code,0,size*2);
  148. for(i = 0; i < size; i++)
  149. {
  150. HM_encode_char(data[i],code+2*i);
  151. }
  152. }
  153. //汉明码解码(多个字节)
  154. uint8_t HM_decode(uint8_t *code, uint8_t *data, uint16_t size)
  155. {
  156. uint8_t err,ret=0;
  157. int i;
  158. memset(data,0,size/2);
  159. for(i = 0; i < size; i=i+2)
  160. {
  161. data[i/2] = HM_decode_char(code+i,&err);
  162. if(err)
  163. ret = 1;
  164. }
  165. return ret;
  166. }
  167. #define NUM 8
  168. void hamming_test(uint8_t num)
  169. {
  170. #if 1
  171. int i,j,temp,temp_bit;
  172. unsigned char coder_T[NUM] = {0};
  173. unsigned char encoder_T[NUM*2] = {0};
  174. unsigned char decoder_T[NUM] = {0};
  175. uint8_t err;
  176. srand(TickCounter);
  177. printf("汉明码测试:\n");
  178. printf("原始数据:");
  179. //生成原数据
  180. for(i = 0; i < NUM; i++)
  181. {
  182. coder_T[i] = rand()%0x100;
  183. printf("%02x ",coder_T[i]);
  184. }
  185. printf("\n");
  186. //编码
  187. // for(i = 0; i < NUM; i++)
  188. // {
  189. // encoder_T[i] = HMBM_74(coder_T[i]);
  190. // }
  191. HM_encode(coder_T,encoder_T,NUM);
  192. printf("编码数据:");
  193. for(i = 0; i < NUM*2; i++)
  194. {
  195. printf("%02x ",encoder_T[i]);
  196. }
  197. printf("\n");
  198. printf("共修改%d位\n",num);
  199. for(i = 0; i < num; i++)
  200. {
  201. temp = rand()%(NUM*8*2);
  202. j = temp/8;
  203. temp_bit = temp%8;
  204. encoder_T[j] ^= 1<<temp_bit;
  205. printf("修改第%02d位\n",temp);
  206. }
  207. // for(i = 0; i < NUM; i++)
  208. // {
  209. // for(j = 0; j < num; j++)
  210. // {
  211. // temp = rand()%7;
  212. // encoder_T[i] ^= 1<<temp;
  213. // printf("修改第%02d位\n",temp);
  214. // }
  215. // }
  216. printf("修改数据:");
  217. for(i = 0; i < NUM*2; i++)
  218. {
  219. printf("%02x ",encoder_T[i]);
  220. }
  221. printf("\n");
  222. //解码
  223. // for(i = 0; i < NUM; i++)
  224. // {
  225. // decoder_T[i] = HMYM_74(encoder_T[i],&err);
  226. // printf("err=%d\n",i);
  227. // if(err)
  228. // {
  229. // printf("%d:err\n",i);
  230. // }
  231. // }
  232. err=HM_decode(encoder_T,decoder_T,NUM*2);
  233. printf("err=%d\n",err);
  234. printf("解码数据:");
  235. for(i = 0; i < NUM; i++)
  236. {
  237. printf("%02x ",decoder_T[i]);
  238. }
  239. printf("\n");
  240. #endif
  241. #if 0
  242. uint32_t data = 0, data_de = 0;
  243. uint32_t data_en[2] = {0};
  244. srand(TickCounter);
  245. printf("汉明码测试:\n");
  246. data = rand()%0xffffffff;
  247. printf("原始数据:%08x\n",data);
  248. HM_encode_u32(data,data_en);
  249. printf("编码数据:%08x, %08x\n",data_en[0],data_en[1]);
  250. data_de = HM_decode_u32(data_en);
  251. printf("解码数据:%08x\n",data_de);
  252. #endif
  253. }