hamming.c 4.4 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269
  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)
  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. switch(c){
  32. case 0 :b = b >> 3 ;break ;
  33. case 1 :b = b >> 3 ;break ;
  34. case 2 :b = b >> 3 ;break ;
  35. case 3 :b = (b^0x08) >> 3 ;break ;
  36. case 4 :b = b >> 3 ;break ;
  37. case 5 :b = (b^0x10) >> 3 ;break ;
  38. case 6 :b = (b^0x20) >> 3 ;break ;
  39. case 7 :b = (b^0x40) >> 3 ;break ;
  40. }
  41. return (b & 0x0F);
  42. }
  43. //对u32数据编码
  44. //data:输入的待编码数据
  45. //code:编码后生成的数据(1字节的原始数据生成2字节的编码数据)
  46. void HM_encode_u32(uint32_t data, uint32_t *code)
  47. {
  48. int i,j;
  49. uint32_t temp = 0;
  50. uint8_t data_char[8] = {0};
  51. uint8_t data_encode[8] = {0};
  52. memset(code,0,2);
  53. // 原始数据拆分:
  54. for(i = 0; i < 8; i++)
  55. {
  56. data_char[i] = (data>>(4*(7-i))) & 0x0f;
  57. }
  58. //编码
  59. for(i = 0; i < 8; i++)
  60. {
  61. data_encode[i] = HMBM_74(data_char[i]);
  62. }
  63. //编码后数据合成
  64. for(i = 0; i < 2; i++)
  65. {
  66. for(j = 0; j < 4; j++)
  67. {
  68. temp = data_encode[i*4+j];
  69. code[i] |= temp<<(3-j)*8;
  70. }
  71. }
  72. }
  73. //对u32数据解码
  74. uint32_t HM_decode_u32(uint32_t *data)
  75. {
  76. int i,j;
  77. uint32_t temp = 0;
  78. uint32_t data_de=0;
  79. uint8_t data_temp[8] = {0};
  80. uint8_t data_decode[8] = {0};
  81. // 原始数据拆分:
  82. for(i = 0; i < 2; i++)
  83. {
  84. for(j = 0; j < 4; j++)
  85. {
  86. data_temp[i*4+j] = (data[i]>>(3-j)*8) & 0xff;
  87. }
  88. }
  89. //解码
  90. for(i = 0; i < 8; i++)
  91. {
  92. data_decode[i] = HMYM_74(data_temp[i]);
  93. }
  94. //解码后数据合成
  95. for(i = 0; i < 8; i++)
  96. {
  97. temp = data_decode[i];
  98. data_de |= temp<<4*(7-i);
  99. }
  100. return data_de;
  101. }
  102. //汉明码编码(一个字节)
  103. void HM_encode_char(uint8_t data, uint8_t *code)
  104. {
  105. int i;
  106. uint8_t data_char[2] = {0};
  107. memset(code,0,2);
  108. // 原始数据拆分:
  109. data_char[0] = data & 0x0f;
  110. data_char[1] = (data>>4) & 0x0f;
  111. //编码
  112. for(i = 0; i < 2; i++)
  113. {
  114. code[i] = HMBM_74(data_char[i]);
  115. }
  116. }
  117. //汉明码解码(一个字节)
  118. uint8_t HM_decode_char(uint8_t *code)
  119. {
  120. int i;
  121. uint8_t data = 0;
  122. uint8_t data_decode[2] = {0};
  123. //解码
  124. for(i = 0; i < 2; i++)
  125. {
  126. data_decode[i] = HMYM_74(code[i]);
  127. }
  128. //解码后数据合成
  129. data |= data_decode[0];
  130. data |= data_decode[1]<<4;
  131. return data;
  132. }
  133. //汉明码编码(多字节)
  134. //返回编码后的字节数
  135. void HM_encode(uint8_t *src, uint8_t *des, uint16_t size)
  136. {
  137. int i;
  138. memset(des,0,size*2);
  139. for(i = 0; i < size; i++)
  140. {
  141. HM_encode_char(src[i],des+2*i);
  142. }
  143. }
  144. //汉明码解码(多个字节)
  145. void HM_decode(uint8_t *src, uint8_t *des, uint16_t size)
  146. {
  147. int i;
  148. memset(des,0,size/2);
  149. for(i = 0; i < size; i=i+2)
  150. {
  151. des[i/2] = HM_decode_char(src+i);
  152. }
  153. }
  154. #define NUM 8
  155. void hamming_test(uint8_t num)
  156. {
  157. #if 0
  158. int i,j,temp,temp_bit;
  159. unsigned char coder_T[NUM] = {0};
  160. unsigned char encoder_T[NUM] = {0};
  161. unsigned char decoder_T[NUM] = {0};
  162. srand(TickCounter);
  163. printf("汉明码测试:\n");
  164. printf("原始数据:");
  165. //生成原数据
  166. for(i = 0; i < NUM; i++)
  167. {
  168. coder_T[i] = rand()%0x10;
  169. printf("%02x ",coder_T[i]);
  170. }
  171. printf("\n");
  172. //编码
  173. for(i = 0; i < NUM; i++)
  174. {
  175. encoder_T[i] = HMBM_74(coder_T[i]);
  176. }
  177. printf("编码数据:");
  178. for(i = 0; i < NUM; i++)
  179. {
  180. printf("%02x ",encoder_T[i]);
  181. }
  182. printf("\n");
  183. printf("共修改%d位\n",num);
  184. for(i = 0; i < num; i++)
  185. {
  186. temp = rand()%(NUM*8);
  187. j = temp/8;
  188. temp_bit = temp%8;
  189. encoder_T[j] ^= 1<<temp_bit;
  190. printf("修改第%02d位\n",temp);
  191. }
  192. printf("修改数据:");
  193. for(i = 0; i < NUM; i++)
  194. {
  195. printf("%02x ",encoder_T[i]);
  196. }
  197. printf("\n");
  198. //解码
  199. for(i = 0; i < NUM; i++)
  200. {
  201. decoder_T[i] = HMYM_74(encoder_T[i]);
  202. }
  203. printf("解码数据:");
  204. for(i = 0; i < NUM; i++)
  205. {
  206. printf("%02x ",decoder_T[i]);
  207. }
  208. printf("\n");
  209. #endif
  210. uint32_t data = 0, data_de = 0;
  211. uint32_t data_en[2] = {0};
  212. srand(TickCounter);
  213. printf("汉明码测试:\n");
  214. data = rand()%0xffffffff;
  215. printf("原始数据:%08x\n",data);
  216. HM_encode_u32(data,data_en);
  217. printf("编码数据:%08x, %08x\n",data_en[0],data_en[1]);
  218. data_de = HM_decode_u32(data_en);
  219. printf("解码数据:%08x\n",data_de);
  220. }