CRC16并行计算的Matlab推导

1、CRC16并行计算的Matlab推导本文使用的CRC16的生成多项式为：x16 + x15 + x2 + 1其对应的串行编码图如以下图所示。假设输入数据的位宽为8比特，即17,16,15,14,13,12,11,10, I为In put的首字母。I0表示最低比特位，I7表示最高比特位。在串行模式下，I0先输入CRC16计算模块，于是I0输入后各个存放 器的状态变化如下：D15 = D14 + D15 + I0D14= D13fD13= D12D12 = D11D11 = D10fD10 = D9D9= DgD8 = D7D7= D6D6 = D5D5 = D4D4= D3D3 = D2D2 =

2、 D15 + D 1 + I0D1 = D0D0 = D15 + I0可以将以上表达式组成矩阵乘法的形式，那么有：D T?D S?I01其中，D为DoDi5构成的列向量，用转置矩阵的形式表示为：T D D0 D1 D2 D3 D4 D5 D6 D7 D8 D9 D10 D11 D12 D13 D14 D15同理，D是D0Di5构成的列向量，用转置矩阵的形式表示为：T DD0D1D2D3D4D5D6D7D8D9D10D11D12D13D14D15表达式 1中的矩阵T，表示为00000000000000011000000000000000010000000000000100100000000000

3、000001000000000000000010000000000000000100000000000000001000000000T00000001000000000000000010000000000000000100000000000000001000000000000000010000000000000000100000000000000001000000000000000011表达式 1中的矩阵 S 也是一个列向量，表示为：TS 1 0 1 0 0 0 0 0 0 0 0 0 0 0 0 1T当 I1 输入后各个存放器的状态为：D 2 T ?D S?I1 T T?D S?I0 S?I

4、0 T2?D T ?S?I0 S?I1以此类推，当 I7 输入后各存放器的状态为：D(8) T8 ?D T7?S?I0 T6?S?I1 T5 ?S?I2 T4?S?I3 T3?S?I4 T2?S?I5 T1?S?I6 T0?S?I7上式中T7?S?I0 T6?S?I1 T5?S?I2 T4?S?I3 T3?S?I 4 T2?S?I5 T1?S?I6 T0?S?I7可以看成两个矩阵A和B的相乘A*B,其中A = T7S T6S T5S T4S T3S T2S T1S T0SI000000000I100000000I200000000I30000B0000I400000000I500000000I

5、600000000I7于是 D(8) 可以简化为D(8) T8 ?D A?B而此时各个存放器中的值就是最后要求的CRC16的值。对应的 matlab 程序如下：T8二mod(D8),2); %计算结果对2取模，是因为我们执行的是二进制%加减运算T7 二mod(D7),2);T6 =mod(TA6),2);T5 =mod(TA5),2);T4 =mod(TA4),2);T3 =mod(TA3),2);T2 =mod(TA2),2);T1 =mod(TA1),2);T0 =mod(TA0),2);S=1,0,0,0,0,0,0,0,0,0,0,0,0,1,0,1;T7S =mod(T7*S),2)

6、;T6S =mod(T6*S),2);T5S =mod(T5*S),2);T4S =mod(T4*S),2);T3S =mod(T3*S),2);T2S =mod(T2*S),2);T1S =mod(T1*S),2);T0S =mod(T0*S),2);A= T7S, T6S, T5S, T4S, T3S, T2S, T1S, T0S;B=eye(8);C=A*B;运算结果:TS豪IiT11111110e90A90nA9AeQflQa1ARnn90000&0e0Q00100090nfliaARa0annt1nnnfl09aaDQ00ea90i00uftnaRD0QQaaRA1R01D09Q9B

7、0QDog9a1011anG0Q00a0aahn1I91iQ000G00aaa9Qa101GGG0aa&aeaElGg0a1100e0aa0aCa0a11QGaa9a0GHsuUU11Uu00HaUU0Qi0BEl611G003a001111111tiunu0HUUi1111II1100e000c 1111111100000Q(I0U000u0cI00Q000Q)c000000CI0U000u0cI01Q0009(I611008S(I0011Q00(I00Q1100EIe000110(I0U00a1TfI00Q000T16000Q0S111111111011111111所以将T8和C带入D(8

8、) T8?D A?B，就能得到各个存放器中最后的值，这也就是CRC16的并行计算实现，最后用 Verilog写出CRC16并行计算的代码:auEign npi(r_crr_e(yL ii assign next_crc_reg1 jGEign npxt_crc_reg2j hssIijii _i：n：_rpij3 jssiqn next crc reg it atfLgn next_crc_rpg5 design next_crt_reg6 dHmiqn next crc rpg|；7 Rzlgn npm crr_rpgR next crc ireq91 (rrr lnLinj (nrc rp

9、gilbj; -(*crc_ini: * (*crc_rea15;?)!=crc_inZj Cri:_lnft -crc in乍 crc_iniij =crc_in3 -crc imCS =crc_lnl -crc inGcrc_in6亠 crcZlREi * crc_inf *i crc_ln(aj crcZin23 * crc in1 crc21nOj * crc req门刃crc_rtg9Jcrc_rc 叭创 1 crc_regl2 crc_rcg13crc_rsg0J; i:fR_rpy9； crc rrq1; crc_rQgH ;crc_reg1 E;crc F?q门片D crc reql31;crc r* crc_rpg 14J 亠 crc_req|11;crr_rpgO;jfsign ne?ttcrc_rpgLlii