C语言纯软件计算任意多项式CRC7、CRC8、CRC16和CRC32的代码

第一节 完全按照列竖式笔算模二除法的方式计算CRC7、CRC8、CRC16和CRC32

条件:初始值INIT=0,结果异或值XOROUT=0,输入数据不反转,输出数据不反转。

cpp 复制代码
/* CRC(循环冗余校验)在线计算: http://www.ip33.com/crc.html */
// 初始值INIT=0, 结果异或值XOROUT=0, 输入数据不反转, 输出数据不反转
#include <stdint.h>
#include <stdio.h>
#include <string.h>
 
#define POLYNOMIAL_CRC7 0x89ul
#define POLYNOMIAL_CRC8 0x107ul
#define POLYNOMIAL_CRC16 0x11021ul
#define POLYNOMIAL_CRC32 0x104c11db7ull
 
/* 计算CRC7校验码 */
uint8_t calc_crc7(const void *data, int len)
{
	const uint8_t *p = data;
	int i, j;
	uint16_t temp = 0;
 
	if (len != 0)
		temp |= p[0] << 8;
 
	for (i = 1; i <= len; i++)
	{
		if (i != len)
			temp |= p[i];
		for (j = 0; j < 8; j++)
		{
			if (temp & 0x8000)
				temp ^= POLYNOMIAL_CRC7 << 8;
			temp <<= 1;
		}
	}
	return temp >> 9;
}
 
/* 计算CRC8校验码 */
uint8_t calc_crc8(const void *data, int len)
{
	const uint8_t *p = data;
	int i, j;
	uint16_t temp = 0;
 
	if (len != 0)
		temp |= p[0] << 8;
 
	for (i = 1; i <= len; i++)
	{
		if (i != len)
			temp |= p[i];
		for (j = 0; j < 8; j++)
		{
			if (temp & 0x8000)
				temp ^= POLYNOMIAL_CRC8 << 7;
			temp <<= 1;
		}
	}
	return temp >> 8;
}
 
/* 计算CRC16校验码 */
uint16_t calc_crc16(const void *data, int len)
{
	const uint8_t *p = data;
	int i, j;
	uint32_t temp = 0;
 
	if (len != 0)
		temp |= p[0] << 24;
	if (len > 1)
		temp |= p[1] << 16; // 填充前二分之一
	if (len > 2)
		temp |= p[2] << 8; // 填充到四分之三
 
	for (i = 3; i <= len + 2; i++)
	{
		if (i < len)
			temp |= p[i]; // 每次都填充最后四分之一的空间
		
		// 从左数第0~7位计算到左数第16~23位
		// 当j=0时是从左数第0位计算到第16位, 共计17位
		// 当j=1时是从左数第1位计算到第17位, 共计17位
		// ...
		// 当j=7时是从左数第7位计算到第23位, 共计17位
		for (j = 0; j < 8; j++)
		{
			if (temp & 0x80000000)
				temp ^= POLYNOMIAL_CRC16 << 15;
			temp <<= 1;
		}
	}
	return temp >> 16;
}
 
/* 计算CRC32校验码 */
uint32_t calc_crc32(const void *data, int len)
{
	const uint8_t *p = data;
	int i, j;
	uint64_t temp = 0;
 
	if (len != 0)
		temp |= (uint64_t)p[0] << 56;
	if (len > 1)
		temp |= (uint64_t)p[1] << 48;
	if (len > 2)
		temp |= (uint64_t)p[2] << 40;
	if (len > 3)
		temp |= (uint64_t)p[3] << 32; // 填充前二分之一
	if (len > 4)
		temp |= (uint64_t)p[4] << 24;
	if (len > 5)
		temp |= (uint64_t)p[5] << 16;
	if (len > 6)
		temp |= (uint64_t)p[6] << 8; // 填充到八分之七
 
	for (i = 7; i <= len + 6; i++)
	{
		if (i < len)
			temp |= p[i]; // 每次都填充最后八分之一的空间
		
		// 从左数第0~7位计算到左数第32~39位
		// 当j=0时是从左数第0位计算到第32位, 共计33位
		// 当j=1时是从左数第1位计算到第33位, 共计33位
		// ...
		// 当j=7时是从左数第7位计算到第39位, 共计33位
		for (j = 0; j < 8; j++)
		{
			if (temp & 0x8000000000000000)
				temp ^= POLYNOMIAL_CRC32 << 31;
			temp <<= 1;
		}
	}
	return temp >> 32;
}
 
void main()
{
	char *str = "https://zh.purasbar.com/";
	int i, n;
 
	n = strlen(str);
	for (i = 0; i <= n; i++)
		printf("len=%d, crc7=0x%02x, crc8=0x%02x, crc16=0x%04x, crc32=0x%08x\n", i, calc_crc7(str, i), calc_crc8(str, i), calc_crc16(str, i), calc_crc32(str, i));
}

程序运行结果:


第二节 计算反转式CRC32

条件:初始值0xffffffff、结果异或值0xffffffff、输入数据反转、输出数据反转。

cpp 复制代码
#include <stdint.h>
#include <stdio.h>
#include <string.h>

#define INIT_CRC32 0xffffffffull
#define POLYNOMIAL_CRC32 0x104c11db7ull

uint8_t reverse(uint8_t data)
{
	int i;
	uint8_t result = 0;

	for (i = 0; i < 8; i++)
	{
		result <<= 1;
		if (data & 1)
			result |= 1;
		data >>= 1;
	}
	return result;
}

uint32_t reverse32(uint32_t data)
{
	uint32_t result;

	result = reverse((data >> 24) & 0xff);
	result |= reverse((data >> 16) & 0xff) << 8;
	result |= reverse((data >> 8) & 0xff) << 16;
	result |= reverse(data & 0xff) << 24;
	return result;
}

/* 计算CRC32校验码(初始值0xffffffff、结果异或值0xffffffff、输入数据反转、输出数据反转) */
uint32_t calc_crc32_rev(const void *data, int len)
{
	const uint8_t *p = data;
	int i, j;
	uint64_t temp = INIT_CRC32 << 32;
 
	if (len != 0)
		temp ^= (uint64_t)reverse(p[0]) << 56;
	if (len > 1)
		temp ^= (uint64_t)reverse(p[1]) << 48;
	if (len > 2)
		temp ^= (uint64_t)reverse(p[2]) << 40;
	if (len > 3)
		temp ^= (uint64_t)reverse(p[3]) << 32; // 填充前二分之一
	if (len > 4)
		temp |= (uint64_t)reverse(p[4]) << 24;
	if (len > 5)
		temp |= (uint64_t)reverse(p[5]) << 16;
	if (len > 6)
		temp |= (uint64_t)reverse(p[6]) << 8; // 填充到八分之七
 
	for (i = 7; i <= len + 6; i++)
	{
		if (i < len)
			temp |= reverse(p[i]); // 每次都填充最后八分之一的空间
		
		// 从左数第0~7位计算到左数第32~39位
		// 当j=0时是从左数第0位计算到第32位, 共计33位
		// 当j=1时是从左数第1位计算到第33位, 共计33位
		// ...
		// 当j=7时是从左数第7位计算到第39位, 共计33位
		for (j = 0; j < 8; j++)
		{
			if (temp & 0x8000000000000000)
				temp ^= POLYNOMIAL_CRC32 << 31;
			temp <<= 1;
		}
	}
	return reverse32(temp >> 32) ^ 0xffffffff;
}

void main()
{
	char *str = "https://zh.purasbar.com/";
	int i, n;

	n = strlen(str);
	for (i = 0; i <= n; i++)
		printf("len=%d, crc32_rev=0x%08x\n", i, calc_crc32_rev(str, i));
}

程序运行结果:


第三节 不使用uint64_t变量计算CRC32

有的单片机不支持uint64_t类型,无法运行calc_crc32函数和calc_crc32_rev函数。

注意到

if (temp & 0x8000000000000000)

temp ^= POLYNOMIAL_CRC32 << 31;

temp <<= 1;

也就是

if (temp & 0x8000000000000000)

{

temp ^= POLYNOMIAL_CRC32 << 31;

temp <<= 1;

}

else

temp <<= 1;

相当于

if (temp & 0x8000000000000000)

{

temp <<= 1;

temp ^= POLYNOMIAL_CRC32 << 32;

}

else

temp <<= 1;

这样的话,我们就可以把一个uint64_t temp变量拆成两个uint32_t变量,使calc_crc32函数和calc_crc32_rev函数可以在不支持uint64_t类型的平台上运行。

因为异或运算由33位变成了32位,所以模二除法的除数POLYNOMIAL_CRC32=0x104c11db7最高位的1就不再需要了,可以删掉,变成0x4c11db7。

cpp 复制代码
#include <stdint.h>
#include <stdio.h>
#include <string.h>
 
#define POLYNOMIAL_CRC32 0x4c11db7
 
/* 计算CRC32校验码 */
uint32_t calc_crc32(const void *data, int len)
{
	const uint8_t *p = data;
	int i, j;
	uint32_t temp = 0;
	uint32_t temp2 = 0;
	uint32_t msb;
 
	// 填充前二分之一
	if (len != 0)
		temp |= p[0] << 24;
	if (len > 1)
		temp |= p[1] << 16;
	if (len > 2)
		temp |= p[2] << 8;
	if (len > 3)
		temp |= p[3];
 
	// 填充到八分之七
	if (len > 4)
		temp2 |= p[4] << 24;
	if (len > 5)
		temp2 |= p[5] << 16;
	if (len > 6)
		temp2 |= p[6] << 8;
 
	for (i = 7; i <= len + 6; i++)
	{
		if (i < len)
			temp2 |= p[i]; // 每次都填充最后八分之一的空间
		
		// 从左数第0~7位计算到左数第32~39位
		// 当j=0时是从左数第0位计算到第32位, 共计33位
		// 当j=1时是从左数第1位计算到第33位, 共计33位
		// ...
		// 当j=7时是从左数第7位计算到第39位, 共计33位
		for (j = 0; j < 8; j++)
		{
			msb = (temp & 0x80000000);
			
			temp <<= 1;
			if (temp2 & 0x80000000)
				temp |= 1;
			temp2 <<= 1;
 
			if (msb)
				temp ^= POLYNOMIAL_CRC32;
		}
	}
	return temp;
}
 
uint8_t reverse(uint8_t data)
{
	int i;
	uint8_t result = 0;
 
	for (i = 0; i < 8; i++)
	{
		result <<= 1;
		if (data & 1)
			result |= 1;
		data >>= 1;
	}
	return result;
}
 
uint32_t reverse32(uint32_t data)
{
	uint32_t result;
 
	result = reverse((data >> 24) & 0xff);
	result |= reverse((data >> 16) & 0xff) << 8;
	result |= reverse((data >> 8) & 0xff) << 16;
	result |= reverse(data & 0xff) << 24;
	return result;
}
 
/* 计算CRC32校验码(初始值0xffffffff、结果异或值0xffffffff、输入数据反转、输出数据反转) */
uint32_t calc_crc32_rev(const void *data, int len)
{
	const uint8_t *p = data;
	int i, j;
	uint32_t temp = 0xffffffff;
	uint32_t temp2 = 0;
	uint32_t msb;
 
	// 填充前二分之一
	if (len != 0)
		temp ^= reverse(p[0]) << 24;
	if (len > 1)
		temp ^= reverse(p[1]) << 16;
	if (len > 2)
		temp ^= reverse(p[2]) << 8;
	if (len > 3)
		temp ^= reverse(p[3]);
 
	// 填充到八分之七
	if (len > 4)
		temp2 |= reverse(p[4]) << 24;
	if (len > 5)
		temp2 |= reverse(p[5]) << 16;
	if (len > 6)
		temp2 |= reverse(p[6]) << 8;
 
	for (i = 7; i <= len + 6; i++)
	{
		if (i < len)
			temp2 |= reverse(p[i]); // 每次都填充最后八分之一的空间
		
		// 从左数第0~7位计算到左数第32~39位
		// 当j=0时是从左数第0位计算到第32位, 共计33位
		// 当j=1时是从左数第1位计算到第33位, 共计33位
		// ...
		// 当j=7时是从左数第7位计算到第39位, 共计33位
		for (j = 0; j < 8; j++)
		{
			msb = (temp & 0x80000000);
			
			temp <<= 1;
			if (temp2 & 0x80000000)
				temp |= 1;
			temp2 <<= 1;
 
			if (msb)
				temp ^= POLYNOMIAL_CRC32;
		}
	}
	return reverse32(temp) ^ 0xffffffff;
}
 
void main()
{
	char *str = "https://zh.purasbar.com/";
	int i, n;
 
	n = strlen(str);
	for (i = 0; i <= n; i++)
		printf("len=%d, crc32=0x%08x, crc32_rev=0x%08x\n", i, calc_crc32(str, i), calc_crc32_rev(str, i));
}

程序运行结果:

相关推荐
恣艺1 分钟前
Python 游戏开发与文件处理:PyGame + Turtle + openpyxl + python-docx + PyPDF2
开发语言·python·pygame
yu85939581 分钟前
STM32 控制 W5500 以太网传输程序
stm32·单片机·嵌入式硬件
LCG元8 分钟前
STM32实战:基于STM32F103的车内防窒息系统(红外检测+GSM报警)
stm32·单片机·嵌入式硬件
高林雨露9 分钟前
kotlin 相关code
开发语言·kotlin
我还记得那天12 分钟前
函数的递归调用
c语言·开发语言·visualstudio
zhangfeng113312 分钟前
ThinkPHP5 事件系统的标准最佳实践 事件系统的完整设计逻辑tags.php tags.php(事件地图)
android·开发语言·php
xyq202416 分钟前
HTML 标签简写及全称
开发语言
tongluowan00717 分钟前
数据结构 Bitmap(位图)示例 - 用户签到系统
开发语言·数据结构·bitmap·用户签到系统
就叫_这个吧17 分钟前
Java线程池应用的四种方式+线程池底层实现原理
java·开发语言
Rust研习社19 分钟前
Rust 官方拟定 LLM 政策,防止 LLM 污染开源社区?
开发语言·后端·ai·rust·开源