优秀的编程知识分享平台

网站首页 > 技术文章 正文

史上解释CRC最清楚的文章

nanyue 2024-11-27 18:16:06 技术文章 1 ℃

CRC的全称是循环冗余校验(Cyclic Redundancy Check),具体的描述可以参考:百度百科:CRC (循环冗余校验),地址为:https://baike.baidu.com/item/CRC/1453359

网上各种各样的讲CRC的内容够多了,本篇文章的目的在于用最简单的方式讲清楚CRC,重在科普、实战和应用,并给出代码实现,数学好的同学不要找我抬杠,谢谢!

CRC的本质是什么?

CRC的目的是保证数据的完整性,其方法是在发送数据的后面再增加多余的若干位数据,接收方使用同样的CRC计算方法,检查接收到的数据CRC是否为0:

  • 如果为0,则表示数据是完整的,接收方可以开开心心的去处理这个数据。
  • 如果不为0,则表示数据不完整/出错,接收方就需要处理下这个数据(一般是丢弃/要求重发)。

那么CRC的核心算法就是如何通过一堆“发送数据”,来计算“多余的若干数据”。

CRC采用的策略是用除法求余数,这就是CRC算法的本质。

不同于十进制的除法运算,CRC用的是二进制的除法运算,二进制的除法运算,属于模二运算,模二运算的特点是:没有进位

模二加法:0+0=0 0+1=1 1+0=1 1+1=0

模二减法:0-0=0 0-1=1 1-0=1 1-1=0

模二乘法:0x0=0 0x1=0 1x0=0 1x1=1

模二除法:是模二乘法的逆运算,参见下图示例

好了,这里我们先给一个例子,我们来计算0x1C的CRC8的校验结果:

如上图,在这个示例中:

  • CRC8的多项式是x8+x2+x+1,对应的除数就是二进制100000111
  • 被除数是0x1C,转化成二进制就是00011100
  • CRC8为8位,被除数后面补8个0
  • 最后的计算结果是0x54

CRC的细节

在上一小节中,我们给出了一个简单的CRC8的例子,细心的同学可能看到了,这里面还有初始值结果异或值输入数据反转输出数据反转,都是什么呢?

简单而言:

  • 初始值是给CRC计算一个初始值,可以是0,也可以是其他值;结果异或值是把计算结果再异或某一个值;这么做的目的是防止全0数据的CRC一直为0
  • 输入数据反转是指输入数据以字节为单位按位逆序处理;输出数据反转是指CRC计算结果整体按位逆序处理;这么做的目的我看到的一个合理解释是右移比左移更容易计算,效率高,它跟大小端无关。

我们再把上面CRC8的例子分解开来看,把CRC校验码如何计算的细节讲清楚:

第一步,我们列出除数,被除数。

第二步,CRC8的输入数据反转为False,所以0x1C仍然是:00011100。

第三步,CRC8补8个0。

第四步,CRC8先计算初始值(0x00),此时被除数保持不变。

第五步,进行模二除法运算,得到结果(0101 0100)。

第六步,CRC8的输出数据反转为False,所以计算结果仍然是:0101 0100。

第七步,CRC的输出异或值(0x00),结果保持不变,最终结果为0x54。

因为模二加减运算的特性,一般计算数包括全0时,步骤可以省略:

  • 初始值为全0时,被除数不变。
  • 结果异或值为0时,结果不变。

如果大家学习过异或操作,可以看到模二加减法实际上就是异或操作,这也是CRC很容易在实际中应用的数学基础。

下面给出了其他的几个示例,包括了数据反转,初始化值和结果异或值,弄清楚了这些例子,CRC就完全搞清楚了。

CRC4-itu,多项式:x4+x+1,输入数据反转,输出数据反转,初始值0x00,输出异或值0x00。

0x1C的CRC4-itu校验码为0x2。

CRC16-modbus,多项式:x16+x15+x2+1,输入数据反转,输出数据反转,初始值0xFFFF,输出异或值0x0000。

0x1C的CRC16-modbus校验码为0x89EB。

CRC5-usb,多项式:x5+x2+1,输入数据反转,输出数据反转,初始值0x1F,输出异或值0x1F。

0x1C的CRC5-usb校验码为0x0D。

总结一下计算CRC的二进制除法过程:

  1. 如果输入数据反转为True,对输入数据按字节进行反转作为被除数;否则直接将输入数据作为被除数。
  2. 被除数补0,CRC码为几位就补几个0,CRC4补4个0,CRC5补5个0,...,CRC32补32个0。
  3. 如果初始值不为0,第一步计算需要先把被除数同初始值做模二加法运算(异或)。
  4. 进行二进制除法运算,求得余数,得到结果。
  5. 如果结果异或值不为0,那么最后的结果需要同结果异或值做异或运算。
  6. 如果输出数据反转为True,那么最终的结果还需要按位进行反转。

最后给一个0x1C2D(2字节)的CRC4-itu示例,供参考。

CRC优化-crc表

由于CRC的计算过程中需要不停的循环做异或运算,占用CPU较多,算法上有一种空间换时间的做法:提前把0x00-0xFF共256个数据的CRC码提前算好保存,那么计算时可以节省CPU,这个提前算好的表叫CRC表(CRC表仅与多项式有关)。

这里特别需要说明下,为什么CRC表仅与多项式有关

如果你理解了CRC的原理,会发现:输入初始值,输入数据反转,输出数据反转,结果异或值都是跟输入输出数据相关的,那么CRC表(指通用的CRC表)仅仅只与多项式(Poly)有关。

但我们在应用中受资源限制,一般会做一些处理,即:把输入数据反转和输出数据反转体现到CRC表中。

我们先来说说CRC表的生成示例,网上给的CRC16表一般是2种,一种最后两个分别是0x0ED1和0x1EF0,另外一种最后两个分别是0x8081和0x4040,实际上这两个分别对应的是CRC16/CCITT_FALSE和CRC16/MODBUS。

看懂了如上的计算原理,那么CRC表的原理就很清楚了;修改算法参数,就可以生成任意的CRC表。

下面给出常用的CRC16的示例和生成CRC的代码,供参考:

private static ushort[] crc16_table_0x1021 =
{
    0x0000, 0x1021, 0x2042, 0x3063, 0x4084, 0x50A5, 0x60C6, 0x70E7,
    0x8108, 0x9129, 0xA14A, 0xB16B, 0xC18C, 0xD1AD, 0xE1CE, 0xF1EF,
    0x1231, 0x0210, 0x3273, 0x2252, 0x52B5, 0x4294, 0x72F7, 0x62D6,
    0x9339, 0x8318, 0xB37B, 0xA35A, 0xD3BD, 0xC39C, 0xF3FF, 0xE3DE,
    0x2462, 0x3443, 0x0420, 0x1401, 0x64E6, 0x74C7, 0x44A4, 0x5485,
    0xA56A, 0xB54B, 0x8528, 0x9509, 0xE5EE, 0xF5CF, 0xC5AC, 0xD58D,
    0x3653, 0x2672, 0x1611, 0x0630, 0x76D7, 0x66F6, 0x5695, 0x46B4,
    0xB75B, 0xA77A, 0x9719, 0x8738, 0xF7DF, 0xE7FE, 0xD79D, 0xC7BC,
    0x48C4, 0x58E5, 0x6886, 0x78A7, 0x0840, 0x1861, 0x2802, 0x3823,
    0xC9CC, 0xD9ED, 0xE98E, 0xF9AF, 0x8948, 0x9969, 0xA90A, 0xB92B,
    0x5AF5, 0x4AD4, 0x7AB7, 0x6A96, 0x1A71, 0x0A50, 0x3A33, 0x2A12,
    0xDBFD, 0xCBDC, 0xFBBF, 0xEB9E, 0x9B79, 0x8B58, 0xBB3B, 0xAB1A,
    0x6CA6, 0x7C87, 0x4CE4, 0x5CC5, 0x2C22, 0x3C03, 0x0C60, 0x1C41,
    0xEDAE, 0xFD8F, 0xCDEC, 0xDDCD, 0xAD2A, 0xBD0B, 0x8D68, 0x9D49,
    0x7E97, 0x6EB6, 0x5ED5, 0x4EF4, 0x3E13, 0x2E32, 0x1E51, 0x0E70,
    0xFF9F, 0xEFBE, 0xDFDD, 0xCFFC, 0xBF1B, 0xAF3A, 0x9F59, 0x8F78,
    0x9188, 0x81A9, 0xB1CA, 0xA1EB, 0xD10C, 0xC12D, 0xF14E, 0xE16F,
    0x1080, 0x00A1, 0x30C2, 0x20E3, 0x5004, 0x4025, 0x7046, 0x6067,
    0x83B9, 0x9398, 0xA3FB, 0xB3DA, 0xC33D, 0xD31C, 0xE37F, 0xF35E,
    0x02B1, 0x1290, 0x22F3, 0x32D2, 0x4235, 0x5214, 0x6277, 0x7256,
    0xB5EA, 0xA5CB, 0x95A8, 0x8589, 0xF56E, 0xE54F, 0xD52C, 0xC50D,
    0x34E2, 0x24C3, 0x14A0, 0x0481, 0x7466, 0x6447, 0x5424, 0x4405,
    0xA7DB, 0xB7FA, 0x8799, 0x97B8, 0xE75F, 0xF77E, 0xC71D, 0xD73C,
    0x26D3, 0x36F2, 0x0691, 0x16B0, 0x6657, 0x7676, 0x4615, 0x5634,
    0xD94C, 0xC96D, 0xF90E, 0xE92F, 0x99C8, 0x89E9, 0xB98A, 0xA9AB,
    0x5844, 0x4865, 0x7806, 0x6827, 0x18C0, 0x08E1, 0x3882, 0x28A3,
    0xCB7D, 0xDB5C, 0xEB3F, 0xFB1E, 0x8BF9, 0x9BD8, 0xABBB, 0xBB9A,
    0x4A75, 0x5A54, 0x6A37, 0x7A16, 0x0AF1, 0x1AD0, 0x2AB3, 0x3A92,
    0xFD2E, 0xED0F, 0xDD6C, 0xCD4D, 0xBDAA, 0xAD8B, 0x9DE8, 0x8DC9,
    0x7C26, 0x6C07, 0x5C64, 0x4C45, 0x3CA2, 0x2C83, 0x1CE0, 0x0CC1,
    0xEF1F, 0xFF3E, 0xCF5D, 0xDF7C, 0xAF9B, 0xBFBA, 0x8FD9, 0x9FF8,
    0x6E17, 0x7E36, 0x4E55, 0x5E74, 0x2E93, 0x3EB2, 0x0ED1, 0x1EF0,
};

生成CRC表的代码如下:

private static ushort[] generate_crc16_table(ushort poly)
{
    ushort crc;
    ushort[] table = new ushort[256];
    for (int index = 0; index < 256; index++)
    {
        crc = (ushort)(index << 8);
        for (int i = 0; i < 8; i++)
        {
            if ((crc & 0x8000) != 0)
            {
                crc = (ushort)((crc << 1) ^ poly);
            }
            else
            {
                crc = (ushort)(crc << 1);
            }
        }
        table[index] = crc;
    }
    return table;
}

代码地址

gitee地址:https://gitee.com/anyangchina/crc_all

参考文献

  • http://www.sunshine2k.de/articles/coding/crc/understanding_crc.html
  • http://users.ece.cmu.edu/~koopman/pubs/KoopmanCRCWebinar9May2012.pdf
  • http://srecord.sourceforge.net/crc16-ccitt.html
  • https://reveng.sourceforge.io/crc-catalogue/all.htm

Tags:

最近发表
标签列表