CRC算法原理及C语言实现
摘 要 本文从理论上推导出CRC算法实现原理,给出三种分别适应不同计算机或微控制器硬件环境的C语言程序。读者更能根据本算法原理,用不同的语言编写出独特风格更加实用的CRC计算程序。
关键词 CRC 算法 C语言
1 引言
循环冗余码CRC检验技术广泛应用于测控及通信领域。CRC计算可以靠专用的硬件来实现,但是对于低成本的微控制器系统,在没有硬件支持下实现CRC检验,关键的问题就是如何通过软件来完成CRC计算,也就是CRC算法的问题。
这里将提供三种算法,它们稍有不同,一种适用于程序空间十分苛刻但CRC计算速度要求不高的微控制器系统,另一种适用于程序空间较大且CRC计算速度要求较高的计算机或微控制器系统,最后一种是适用于程序空间不太大,且CRC计算速度又不可以太慢的微控制器系统。
2 CRC简介
CRC校验的基本思想是利用线性编码理论,在发送端根据要传送的k位二进制码序列,以一定的规则产生一个校验用的监督码(既CRC码)r位,并附在信息后边,构成一个新的二进制码序列数共(k+r)位,最后发送出去。在接收端,则根据信息码和CRC码之间所遵循的规则进行检验,以确定传送中是否出错。
16位的CRC码产生的规则是先将要发送的二进制序列数左移16位(既乘以 )后,再除以一个多项式,最后所得到的余数既是CRC码,如式(2-1)式所示,其中B(X)表示n位的二进制序列数,G(X)为多项式,Q(X)为整数,R(X)是余数(既CRC码)。
(2-1)
求CRC码所采用模2加减运算法则,既是不带进位和借位的按位加减,这种加减运算实际上就是逻辑上的异或运算,加法和减法等价,乘法和除法运算与普通代数式的乘除法运算是一样,符合同样的规律。生成CRC码的多项式如下,其中CRC-16和CRC-CCITT产生16位的CRC码,而CRC-32则产生的是32位的CRC码。本文不讨论32位的CRC算法,有兴趣的朋友可以根据本文的思路自己去推导计算方法。
CRC-16:(美国二进制同步系统中采用)
CRC-CCITT:(由欧洲CCITT推荐)
CRC-32:
接收方将接收到的二进制序列数(包括信息码和CRC码)除以多项式,如果余数为0,则说明传输中无错误发生,否则说明传输有误,关于其原理这里不再多述。用软件计算CRC码时,接收方可以将接收到的信息码求CRC码,比较结果和接收到的CRC码是否相同。
3 按位计算CRC
对于一个二进制序列数可以表示为式(3-1):
(3-1)
求此二进制序列数的CRC码时,先乘以 后(既左移16位),再除以多项式G(X),所得的余数既是所要求的CRC码。如式(3-2)所示:
(3-2)
可以设: (3-3)
其中 为整数, 为16位二进制余数。将式(3-3)代入式(3-2)得:
(3-4)
再设: (3-5)
其中 为整数, 为16位二进制余数,将式(3-5)代入式(3-4),如上类推,最后得到:
(3-6)
根据CRC的定义,很显然,十六位二进制数 既是我们要求的CRC码。
式(3-5)是编程计算CRC的关键,它说明计算本位后的CRC码等于上一位CRC码乘以2后除以多项式,所得的余数再加上本位值除以多项式所得的余数。由此不难理解下面求CRC码的C语言程序。*ptr指向发送缓冲区的首字节,len是要发送的总字节数,0x1021与多项式有关。
unsigned int cal_crc(unsigned char *ptr, unsigned char len) {
unsigned char i;
unsigned int crc=0;
while(len--!=0) {
for(i=0x80; i!=0; i/=2) {
if((crc&;0x8000)!=0) {crc*=2; crc^=0x1021;} /* 余式CRC乘以2再求CRC */
else crc*=2;
if((*ptr&;i)!=0) crc^=0x1021; /* 再加上本位的CRC */
}
ptr++;
}
return(crc);
}
按位计算CRC虽然代码简单,所占用的内存比较少,但其最大的缺点就是一位一位地计算会占用很多的处理器处理时间,尤其在高速通讯的场合,这个缺点更是不可容忍。因此下面再介绍一种按字节查表快速计算CRC的方法。
4 按字节计算CRC
不难理解,对于一个二进制序列数可以按字节表示为式(4-1),其中 为一个字节(共8位)。
(4-1)
求此二进制序列数的CRC码时,先乘以 后(既左移16位),再除以多项式G(X),所得的余数既是所要求的CRC码。如式(4-2)所示:
(4-2)
可以设: (4-3)
其中 为整数, 为16位二进制余数。将式(4-3)代入式(4-2)得:
(4-4)
因为:
(4-5)
其中 是 的高八位, 是 的低八位。将式(4-5)代入式(4-4),经整理后得:
(4-6)
再设: (4-7)
其中 为整数, 为16位二进制余数。将式(4-7)代入式(4-6),如上类推,最后得:
(4-8)
很显然,十六位二进制数 既是我们要求的CRC码。
式(4-7)是编写按字节计算CRC程序的关键,它说明计算本字节后的CRC码等于上一字节余式CRC码的低8位左移8位后,再加上上一字节CRC右移8位(也既取高8位)和本字节之和后所求得的CRC码,如果我们把8位二进制序列数的CRC全部计算出来,放如一个表里,采用查表法,可以大大提高计算速度。由此不难理解下面按字节求CRC码的C语言程序。*ptr指向发送缓冲区的首字节,len是要发送的总字节数,CRC余式表是按0x11021多项式求出的。
unsigned int cal_crc(unsigned char *ptr, unsigned char len) {
unsigned int crc;
unsigned char da;
unsigned int crc_ta[256]={ /* CRC余式表 */
0x0000, 0x1021, 0x2042, 0x3063, 0x4084, 0x50a5, 0x60c6, 0x70e7,
0x8108, 0x9129, 0xa14a, 0xb16b, 0xc18c, 0xd1ad, 0xe1ce, 0xf1ef,
0x 1231, 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=0;
while(len--!=0) {
da=(uchar) (crc/256); /* 以8位二进制数的形式暂存CRC的高8位 */
crc<<=8; /* 左移8位,相当于CRC的低8位乘以 */
crc^=crc_ta[da^*ptr]; /* 高8位和当前字节相加后再查表求CRC ,再加上以前的CRC */
ptr++;
}
return(crc);
}
很显然,按字节求CRC时,由于采用了查表法,大大提高了计算速度。但对于广泛运用的8位微处理器,代码空间有限,对于要求256个CRC余式表(共512字节的内存)已经显得捉襟见肘了,但CRC的计算速度又不可以太慢,因此再介绍下面一种按半字节求CRC的算法。
5 按半字节计算CRC
同样道理,对于一个二进制序列数可以按字节表示为式(5-1),其中 为半个字节(共4位)。
(5-1)
求此二进制序列数的CRC码时,先乘以 后(既左移16位),再除以多项式G(X),所得的余数既是所要求的CRC码。如式(4-2)所示:
(5-2)
可以设: (5-3)
其中 为整数, 为16位二进制余数。将式(5-3)代入式(5-2)得:
(5-4)
因为:
(5-5)
其中 是 的高4位, 是 的低12位。将式(5-5)代入式(5-4),经整理后得:
(5-6)
再设: (5-7)
其中 为整数, 为16位二进制余数。将式(5-7)代入式(5-6),如上类推,最后得:
(5-8)
很显然,十六位二进制数 既是我们要求的CRC码。
式(5-7)是编写按字节计算CRC程序的关键,它说明计算本字节后的CRC码等于上一字节CRC码的低12位左移4位后,再加上上一字节余式CRC右移4位(也既取高4位)和本字节之和后所求得的CRC码,如果我们把4位二进制序列数的CRC全部计算出来,放在一个表里,采用查表法,每个字节算两次(半字节算一次),可以在速度和内存空间取得均衡。由此不难理解下面按半字节求CRC码的C语言程序。*ptr指向发送缓冲区的首字节,len是要发送的总字节数,CRC余式表是按0x11021多项式求出的。
unsigned cal_crc(unsigned char *ptr, unsigned char len) {
unsigned int crc;
unsigned char da;
unsigned int crc_ta[16]={ /* CRC余式表 */
0x0000,0x1021,0x2042,0x3063,0x4084,0x50a5,0x60c6,0x70e7,
0x8108,0x9129,0xa14a,0xb16b,0xc18c,0xd1ad,0xe1ce,0xf1ef,
}
crc=0;
while(len--!=0) {
da=((uchar)(crc/256))/16; /* 暂存CRC的高四位 */
crc<<=4; /* CRC右移4位,相当于取CRC的低12位)*/
crc^=crc_ta[da^(*ptr/16)]; /* CRC的高4位和本字节的前半字节相加后查表计算CRC,
然后加上上一次CRC的余数 */
da=((uchar)(crc/256))/16; /* 暂存CRC的高4位 */
crc<<=4; /* CRC右移4位, 相当于CRC的低12位) */
crc^=crc_ta[da^(*ptr&;0x0f)]; /* CRC的高4位和本字节的后半字节相加后查表计算CRC,
然后再加上上一次CRC的余数 */
ptr++;
}
return(crc);
}
2. CRC验证可以用在单片机里吗
可以是可以,不过这段程序在8位的51单片机里直接编译的话效率太低。
一般都是用等价8位查表法来实现CRC16的。
3. 超声波测距离
一种用于汽车倒车避撞的超声波无线距离测量系统
Research of Ultrasonic Distance Measurement System
Abstract: A kind of ultrasonic distance measurement system used in the car is designed in this paper. The system includes the lower microcomputer system and the upper microcomputer system. The lower microcomputer system is mainly composed of ultrasonic transmitting circuit, receiving circuit ,wireless communicating mole and microcomputer. The data from the lower microcomputer system is transmitted to the upper microcomputer system by the wireless way. The design principle of ultrasonic distance measurement circuit is analyzed. The design method that the data is transmitted is also introced. The system is of the characteristics of measurement convenience, fast response and stability.
Key words : wireless communicating;microcomputer; ultrasonic;distance measurement;temperature compensation
摘 要:本文介绍一种用于汽车倒车避撞的超声波无线距离测量系统。系统由下位机与上位机两部分组成,下位机主要由超声波发射电路、超声波接收电路、无线收发模块及单片机组成,上位机由单片机、无线收发模块、显示电路等组成,下位机与上位机之间通过无线收发模块传输信息。文中分析了超声波测距电路的设计方法,叙述了采用无线通信技术实现数据远程传输的设计思路。该系统测量距离方便、灵活、稳定。
关键词:无线通信;单片机;超声波;距离测量;温度补偿
1. 引言
随着经济的发展,人们的生活水平越来越高。当今,对许多人来说,汽车进入家庭已不再是奢望,但随之而来的事情就是如何保证汽车使用过程中的安全问题,特别是如何防止汽车与其他物体碰撞的事情发生。据初步调查统计,l5%的汽车事故是由汽车倒车“后视”不良造成的。因此,增强汽车的后视能力,对于提高行车安全,减轻司机的劳动强度和心理压力,是十分重要的。如果车辆能适时检测与周围障碍物的距离并给出警告信息,使司机及早采取行动,可避免车辆相撞事故的发生。
随着科学技术的发展,用超声波进行无接触测量得到了广泛的应用。超声波是由机械振动产生的,可在不同介质中以不同的速度传播,它具有定向性好、能量集中、在传输过程中衰减较小,反射能力较强,在恶劣工作环境下具有一定的适应能力等优点。因此可用于液位测量、车辆自动导航[2]等领域。本文介绍一种基于无线数据传输方式的超声波车辆倒车避撞预警系统。
2. 超声波测距原理
发射的超声波遇到障碍物时就会发生反射,反射波可由接收器接收,这样只要测出超声波从发送点到反射回来的时间间隔Δt,然后根据公式(1)即可求出超生波从发射处到障碍物之间的距离。
S=CΔt/2 (1)
式中:S—超生波发射处与障碍物间的距离
C—超声波在介质中的传播速度
由于超声波是一种声波,其声速C受环境温度的影响,关系如式(2),因此使用超生波测量距离时应该采用温度补偿的方法对式(1)中的声速值加以校正。
C=331.4+ 0.61×T (2)
式中:T—环境温度
3. 硬件电路设计
如图1,硬件电路主要由单片机、超声波传感器、温度测量电路、无线收发模块等组成。
系统中单片机均采用ATMEL公司的AT89S51作为核心控制芯片,它与MCS-51的指令和引脚兼容[1],并且具有ISP在线编程功能,便于系统的设计和调试。
超声波传感器是超声波测距电路中的重要元件,其性能优劣直接影响到测距准确度和可靠性。通常超声波传感器有两类:一类是发射电路和接收电路互相独立的分体式超声波传感器,此类传感器测距有效范围比较大,但不具备防尘、防水性能。另一类是同时具有发射与接收功能的收发一体式超声波传感器,此类超声波测距有效范围比较小,但防尘、防水性能好。该系统选择分体式超声波传感器。
考虑到超声波具有指向性,本系统在汽车尾部左、右两个部位各安装一个超声波传
感器,适当调整安装位置,可准确测量汽
车后部障碍物。
如图1所示,下位机的P1.1、P1.2引脚分别用于控制两路超声波发射,INT0,INT1分别用于两路超声波信号检测,P1.3用于温度检测,串口RXD、TXD分别连接无线收发模块A的输入、输出端。同样,上位机串口RXD、TXD分别连接无线收发模块B的输入、输出端,当接收到下位机发送的测量数据时,下位机进行处理,然后显示测量结果,当车辆离障碍物的距离超过安全警戒线时发出报警信号。
实际安装时,该系统的下位机部分安装在汽车的尾部,上位机部分安装于驾驶室内。
3.1 超声波发射电路
超声波发射电路由超声波换能器(或称超声波振头)和超声波发生器两部分组成,电路如图2所示。系统中,超声波换能器的型号为CSB40T,它将超声波发生器提供的电信号转换为机械振动并发射出去。40KHz的超声波信号是利用NE555时基电路振荡产生的,振荡频率f ≈1.44/((R22+2×R23)×C21),通过R23调节信号频率,使之与换能器的40KHz固有频率一致。工作时,下位机通过P1.1口定时向超声波发生电路发出控制信号,超声波发生电路产生40KHz的调制脉冲,经换能器转换为超声波信号向前方空间发射。
3.2 超声波接收电路
超声波接收电路采用了集成电路CX20106A,CX20106A是日本索尼公司生产的红外遥控信号接收集成电路,它由前置放大、自动偏压控制、振幅放大、峰值检波和整形电路组成。该集成电路红外发射的频率38KHZ,超声波换能器的固有频率是40KHz,适当设计CX20106A外围电路参数,就可以将其用于超声波的接收放大电路,如图3所示,引脚1为CX20106A信号输入端,引脚2为CX20106A的RC网络连接端,引脚3为CX20106A检波电容连接端,
引脚4为CX20106A的接地端,引脚5为CX20106A带通滤波器中心设置端,引脚6为CX20106A积分电容连接端,引脚7为CX20106A信号输出端,引脚8为CX20106A供电电源端。
工作时,换能器CSB40T将所接收到的微弱声波振动信号转化成为电信号,送给CX20106A的输入端1,当CX20106A接收到信号时,7脚就会输出一个低电平,可用于下位机的中断信号源。当下位机接收到中断信号时,说明检测到了反射回来的超声波,下位机就进入中断状态,开始距离计算,并将计算结果发送给上位机。
3.3温度检测电路
温度检测电路采用DALLS公司的1-WIRE式总线器件DS18B20数字温度传感器,电路连接非常简单,但是必须保证时序与单片机严格同步。DS18B20具有9,10,11,和12位转换精度,未编程时默认精度为12位,测量精度一般为0.5℃,软件处理后可达0.1℃。温度输出以16位符号扩展的二进制数形式提供,低位在先,以0.0625℃/ LSB形式表达,高五位为扩展符号位。转换周期与转换精度设定有关,9位精度时,最大转换时间为93.75ms;12位精度时,最大转化时间为750ms。在本系统中采用默认的12位精度。关于DS18B20的使用方法可参考有关书籍。
3.4 数据无线收发模块
为避免在车内铺设电缆,系统的上位机部分与下位机部分采用无线的方式进行通信。
无线通信模块采用PTR2000,它是收发一体的工作在国际通用数传频段433MHz的无线通信模块,最高传输速率可以达到20Kbit/s,功耗低,待机状态下仅为8μA,可以直接与单片机的串口连接使用。PTR2000的引脚定义如下:TXE是发送控制端;PWR是节能控制端;DI是数据输入端;DO是数据输出端;CS是频道选择端。
硬件连接时,由单片机3个通用I/O口分别控制TXE、PWR、CS,单片机的串口与DI,DO连接。TXE为1时,为发送状态,TXE为0时,为接收状态。状态转换需要5毫秒。PWR为0时,为节电待机状态,此时模块无法进行接收或者发送。
无线通信具有无需布线、便于安装、灵活性强等诸多优点,但是数据在传输过程中难以避免的会产生误码,而且产生误码的几率要远远大于有线网络,并且误码的产生与多方面的因素有关,因此有很大的不确定性。所以必须采用一种差错控制机制,该系统采用停止等待协议来实现差错控制。此外,还采用校验机制以确定何时需要重传,CRC校验码的检错能力很强,它除了能检查出离散传输错误外,还能检查出突发传输错误。考虑到硬件和传输的开销问题,使用CRC16校验码。
PTR2000灵敏性很高,在无载波的情况下在接收端会产生随机的数据,因此需制定传输协议,格式如表1所示。通信协议中,必须在有效数据前加上两个或多个固定的前导字符作为同步信号,使得接收端能够辨别出有效数据的开始。
前导字符采用0xAA、0xAA、0xFF、0x00共4字节,其中前两个字节为同步信号,后两个字节为帧开始标志,接收端只要能够接收到0xAA、0xAA、0xFF与0x00,就可以认为新的一帧开始了。帧类型分为数据帧、有序数据帧、控制命令帧、确认帧等多种帧类型。帧编号为可选项,与帧类型相关,只有帧类型是有序数据帧时才有效。校验为2字节CRC16校验码。帧结束标志:为0x00。
4.软件设计
4.1下位机程序设计
下位机程序主要由数据通信程序、距离计算程序、温度补偿程序等组成。
距离计算程序流程图如图4所示。
温度补偿通常有两种方法:一种方法是每次按照公式C=331.4+ 0.61×T计算当前声速C,进行温度补偿。其特点是:根据当时的温度得到精确声速,从而计算得到的距离值也比较精确,但程序中牵涉到浮点数运算,由微处理器系统实现,难度较大。另一种方法是将温度与声速的对应关系列成温度---声速二维表,固化到系统中。温度补偿时,根据温度---声速表,查取最接近当前温度的那个温度所对应的声速值,此声速值即作为当前声速。其特点是:避开了复杂的浮点运算和浮点运算后各字节的提取操作,这样既保证了一定的精度要求,又可以避免浮点运算。因此本系统采用方法二进行温度补偿。程序流程图略。
4.2 上位机主程序设计
上位机与下位机通信时,上位机按照通信协议格式将开始测量命令发送给下位机,下位机接收到命令后就开始测量汽车离障碍物的距离,然后将测量结果发送给上位机,上位机先判断前导字符来确定是否为有效数据,若是有效数据,则解开封包进行相应操作,否则丢弃该数据包,上位机再按照同样的方式继续发命令、接收数据,直到接收到正确的数据为止。程序流程图如图5所示。
5 结束语
通过对系统硬件电路和软件的合理设计,本系统能在-20℃到50℃之间正常工作,三位数码管以厘米为单位显示距离,能准确判断距离汽车1.5米内的物体并及时报警,提高了汽车倒车的安全性。本文的创新点是在汽车防撞系统中采用了数据无线通信策略,减少了车内布线。
参考文献:
[1]MCS-51系列单片机应用系统设计,何立民,北京航空航天大学出版社,1990年.
[2]高准确度超声波测距仪的研制,赵珂,向瑛等,传感器技术,2003年第2期.
[3]无线通信在嵌入式系统中的应用,曹玲芝,石军等,微计算机信息,2005年第11期
作者简介:
曹玲芝,女,1965年生,硕士,副教授,主要从事远程测控技术研究。
联系方式:郑州轻工业学院电气信息工程学院办公室 邮编 450002
Email: [email protected]
任亚萍,女,1973年生,硕士生,主要从事计算机技术研究.
4. 用单片机汇编语言或C51语言或实现CRC校验码
Uint16 Crc16(unsigned char *puchMsg, int usDataLen)
{
unsigned int uchCRCHi = 0xFF ; /* 高CRC字节初始化 */
unsigned int uchCRCLo = 0xFF ; /* 低CRC 字节初始化 */
unsigned int temp16;
Uint32 uIndex ; /* CRC循环中的索引 */
while (usDataLen--) /* 传输消息缓冲区 */
{
temp16=*puchMsg++;
uIndex = uchCRCHi ^ temp16 ; /* 计算CRC */
uchCRCHi = uchCRCLo ^ auchCRCHi[uIndex] ;
uchCRCLo = auchCRCLo[uIndex] ;
}
return (uchCRCHi << 8 | uchCRCLo) ;
}
5. 几种CRC16算法
一. CRC16算法首先在源文件头文件加入表值:[c] view plain ////////////////////////////////////////////////////////////////////////// // CRC16码表 static WORD const wCRC16Table[256] = { 0x0000, 0xC0C1, 0xC181, 0x0140, 0xC301, 0x03C0, 0x0280, 0xC241, 0xC601, 0x06C0, 0x0780, 0xC741, 0x0500, 0xC5C1, 0xC481, 0x0440, 0xCC01, 0x0CC0, 0x0D80, 0xCD41, 0x0F00, 0xCFC1, 0xCE81, 0x0E40, 0x0A00, 0xCAC1, 0xCB81, 0x0B40, 0xC901, 0x09C0, 0x0880, 0xC841, 0xD801, 0x18C0, 0x1980, 0xD941, 0x1B00, 0xDBC1, 0xDA81, 0x1A40, 0x1E00, 0xDEC1, 0xDF81, 0x1F40, 0xDD01, 0x1DC0, 0x1C80, 0xDC41, 0x1400, 0xD4C1, 0xD581, 0x1540, 0xD701, 0x17C0, 0x1680, 0xD641, 0xD201, 0x12C0, 0x1380, 0xD341, 0x1100, 0xD1C1, 0xD081, 0x1040, 0xF001, 0x30C0, 0x3180, 0xF141, 0x3300, 0xF3C1, 0xF281, 0x3240, 0x3600, 0xF6C1, 0xF781, 0x3740, 0xF501, 0x35C0, 0x3480, 0xF441, 0x3C00, 0xFCC1, 0xFD81, 0x3D40, 0xFF01, 0x3FC0, 0x3E80, 0xFE41, 0xFA01, 0x3AC0, 0x3B80, 0xFB41, 0x3900, 0xF9C1, 0xF881, 0x3840, 0x2800, 0xE8C1, 0xE981, 0x2940, 0xEB01, 0x2BC0, 0x2A80, 0xEA41, 0xEE01, 0x2EC0, 0x2F80, 0xEF41, 0x2D00, 0xEDC1, 0xEC81, 0x2C40, 0xE401, 0x24C0, 0x2580, 0xE541, 0x2700, 0xE7C1, 0xE681, 0x2640, 0x2200, 0xE2C1, 0xE381, 0x2340, 0xE101, 0x21C0, 0x2080, 0xE041, 0xA001, 0x60C0, 0x6180, 0xA141, 0x6300, 0xA3C1, 0xA281, 0x6240, 0x6600, 0xA6C1, 0xA781, 0x6740, 0xA501, 0x65C0, 0x6480, 0xA441, 0x6C00, 0xACC1, 0xAD81, 0x6D40, 0xAF01, 0x6FC0, 0x6E80, 0xAE41, 0xAA01, 0x6AC0, 0x6B80, 0xAB41, 0x6900, 0xA9C1, 0xA881, 0x6840, 0x7800, 0xB8C1, 0xB981, 0x7940, 0xBB01, 0x7BC0, 0x7A80, 0xBA41, 0xBE01, 0x7EC0, 0x7F80, 0xBF41, 0x7D00, 0xBDC1, 0xBC81, 0x7C40, 0xB401, 0x74C0, 0x7580, 0xB541, 0x7700, 0xB7C1, 0xB681, 0x7640, 0x7200, 0xB2C1, 0xB381, 0x7340, 0xB101, 0x71C0, 0x7080, 0xB041, 0x5000, 0x90C1, 0x9181, 0x5140, 0x9301, 0x53C0, 0x5280, 0x9241, 0x9601, 0x56C0, 0x5780, 0x9741, 0x5500, 0x95C1, 0x9481, 0x5440, 0x9C01, 0x5CC0, 0x5D80, 0x9D41, 0x5F00, 0x9FC1, 0x9E81, 0x5E40, 0x5A00, 0x9AC1, 0x9B81, 0x5B40, 0x9901, 0x59C0, 0x5880, 0x9841, 0x8801, 0x48C0, 0x4980, 0x8941, 0x4B00, 0x8BC1, 0x8A81, 0x4A40, 0x4E00, 0x8EC1, 0x8F81, 0x4F40, 0x8D01, 0x4DC0, 0x4C80, 0x8C41, 0x4400, 0x84C1, 0x8581, 0x4540, 0x8701, 0x47C0, 0x4680, 0x8641, 0x8201, 0x42C0, 0x4380, 0x8341, 0x4100, 0x81C1, 0x8081, 0x4040}; 然后在文件中加入下列函数:[c] view plain ////////////////////////////////////////////////////////////////////////// // 函数功能: CRC16效验 // 输入参数: pDataIn: 数据地址 // iLenIn: 数据长度 // 输出参数: pCRCOut: 2字节校验值 void CCRCDlg::CRC16(const CHAR* pDataIn, int iLenIn, WORD* pCRCOut) { WORD wResult = 0; WORD wTableNo = 0; for(int i = 0; i < iLenIn; i++) { wTableNo = ((wResult & 0xff) ^ (pDataIn[i] & 0xff)); wResult = ((wResult >> 8) & 0xff) ^ wCRC16Table[wTableNo]; } *pCRCOut = wResult; } 二.CRC16(MODBUS)[c] view plain ////////////////////////////////////////////////////////////////////////// // CRC MODBUS 效验 // 输入参数: pDataIn: 数据地址 // iLenIn: 数据长度 // 输出参数: pCRCOut: 2字节校验值 void CCRCDlg::CheckCRCModBus(const CHAR* pDataIn, int iLenIn, WORD* pCRCOut) { WORD wHi = 0; WORD wLo = 0; WORD wCRC; wCRC = 0xFFFF; for (int i = 0; i < iLenIn; i++) { wCRC = CalcCRCModBus(*pDataIn, wCRC); pDataIn++; } wHi = wCRC / 256; wLo = wCRC % 256; wCRC = (wHi > 1; wCRCIn = wCRCIn & 0x7fff; if(wCheck == 1) { wCRCIn = wCRCIn ^ 0xa001; } wCRCIn = wCRCIn & 0xffff; } return wCRCIn; } 三.CRC16(CCITT的0XFFFF)[c] view plain ////////////////////////////////////////////////////////////////////////// // 函数功能: CRC16效验(CCITT的0XFFFF效验) // 输入参数: pDataIn: 数据地址 // iLenIn: 数据长度 // 输出参数: pCRCOut: 2字节校验值 void CCRCDlg::CRCCCITT(const CHAR* pDataIn, int iLenIn, WORD* pCRCOut) { WORD wTemp = 0; WORD wCRC = 0xffff; for(int i = 0; i < iLenIn; i++) { for(int j = 0; j < 8; j++) { wTemp = ((pDataIn[i] > 8); wCRC
6. crc16校验的c语言程序
下面我们以CRC-16为例来说明任意长度数据流的CRC校验码生成过程。我们采用将数据流分成若干个8bit字符,并由低字节到高字节传送的并行方法来求CRC校验码。具体计算过程为:用一个16bit的寄存器来存放CRC校验值,且设定其初值为0x0000;将数据流的第一个8bit与16bit的CRC寄存器的高字节相异或,并将结果存入CRC寄存器高字节;CRC寄存器左移一位,最低1bit补零,同时检查移出的最高1bit,若移出的最高1bit为0,则继续按上述过程左移,若最高1bit为1,则将CRC寄存器中的值与生成多项式码相异或,结果存入CRC寄存器值;继续左移并重复上述处理方法,直到将8bit数据处理完为止,则此时CRC寄存器中的值就是第一个8bit数据对应的CRC校验码;然后将此时CRC寄存器的值作为初值,用同样的处理方法重复上述步骤来处理下一个8bit数据流,直到将所有的8bit字符都处理完后,此刻CRC寄存器中的值即为整个数据流对应的CRC校验码。
下面示出了其计算过程的流程图:
在用C语言编写CRC校验码的实现程序时我们应该注意,生成多项式 对应的十六进制数为0x18005,由于CRC寄存器左移过程中,移出的最高位为1时与 相异或,所以与16bit的CRC寄存器对应的生成多项式的十六进制数可用0x8005表示。下面给出并行处理8bit数据流的C源程序:
unsigned short crc_dsp(unsigned short reg, unsigned char data_crc)
//reg为crc寄存器, data_crc为将要处理的8bit数据流
{
unsigned short msb; //crc寄存器将移出的最高1bit
unsigned short data;
unsigned short gx = 0x8005, i = 0; //i为左移次数, gx为生成多项式
data = (unsigned short)data_crc;
data = data << 8;
reg = reg ^ data;
do
{
msb = reg & 0x8000;
reg = reg << 1;
if(msb == 0x8000)
{
reg = reg ^ gx;
}
i++;
}
while(i < 8);
return (reg);
}
以上为处理每一个8bit数据流的子程序,在计算整个数据流的CRC校验码时,我们只需将CRC_reg的初值置为0x0000,求第一个8bit的CRC值,之后,即可将上次求得的CRC值和本次将要处理的8bit数据作为函数实参传递给上述子程序的形参进行处理即可,最终返回的reg值便是我们所想得到的整个数据流的CRC校验值。
7. 求教 C#语言编写 的CRC16的校验程序 (多项式为:CRC-16/X25 x16+x12+x5+1)
using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;
namespace ConvertToCRC16
{
public static class CRC16Util
{
// CRC高位字节表
private static readonly byte[] m_CRCHighOrderByteTable = new byte[]
{
0x00, 0xC1, 0x81, 0x40, 0x01, 0xC0, 0x80, 0x41, 0x01, 0xC0,
0x80, 0x41, 0x00, 0xC1, 0x81, 0x40, 0x01, 0xC0, 0x80, 0x41,
0x00, 0xC1, 0x81, 0x40, 0x00, 0xC1, 0x81, 0x40, 0x01, 0xC0,
0x80, 0x41, 0x01, 0xC0, 0x80, 0x41, 0x00, 0xC1, 0x81, 0x40,
0x00, 0xC1, 0x81, 0x40, 0x01, 0xC0, 0x80, 0x41, 0x00, 0xC1,
0x81, 0x40, 0x01, 0xC0, 0x80, 0x41, 0x01, 0xC0, 0x80, 0x41,
0x00, 0xC1, 0x81, 0x40, 0x01, 0xC0, 0x80, 0x41, 0x00, 0xC1,
0x81, 0x40, 0x00, 0xC1, 0x81, 0x40, 0x01, 0xC0, 0x80, 0x41,
0x00, 0xC1, 0x81, 0x40, 0x01, 0xC0, 0x80, 0x41, 0x01, 0xC0,
0x80, 0x41, 0x00, 0xC1, 0x81, 0x40, 0x00, 0xC1, 0x81, 0x40,
0x01, 0xC0, 0x80, 0x41, 0x01, 0xC0, 0x80, 0x41, 0x00, 0xC1,
0x81, 0x40, 0x01, 0xC0, 0x80, 0x41, 0x00, 0xC1, 0x81, 0x40,
0x00, 0xC1, 0x81, 0x40, 0x01, 0xC0, 0x80, 0x41, 0x01, 0xC0,
0x80, 0x41, 0x00, 0xC1, 0x81, 0x40, 0x00, 0xC1, 0x81, 0x40,
0x01, 0xC0, 0x80, 0x41, 0x00, 0xC1, 0x81, 0x40, 0x01, 0xC0,
0x80, 0x41, 0x01, 0xC0, 0x80, 0x41, 0x00, 0xC1, 0x81, 0x40,
0x00, 0xC1, 0x81, 0x40, 0x01, 0xC0, 0x80, 0x41, 0x01, 0xC0,
0x80, 0x41, 0x00, 0xC1, 0x81, 0x40, 0x01, 0xC0, 0x80, 0x41,
0x00, 0xC1, 0x81, 0x40, 0x00, 0xC1, 0x81, 0x40, 0x01, 0xC0,
0x80, 0x41, 0x00, 0xC1, 0x81, 0x40, 0x01, 0xC0, 0x80, 0x41,
0x01, 0xC0, 0x80, 0x41, 0x00, 0xC1, 0x81, 0x40, 0x01, 0xC0,
0x80, 0x41, 0x00, 0xC1, 0x81, 0x40, 0x00, 0xC1, 0x81, 0x40,
0x01, 0xC0, 0x80, 0x41, 0x01, 0xC0, 0x80, 0x41, 0x00, 0xC1,
0x81, 0x40, 0x00, 0xC1, 0x81, 0x40, 0x01, 0xC0, 0x80, 0x41,
0x00, 0xC1, 0x81, 0x40, 0x01, 0xC0, 0x80, 0x41, 0x01, 0xC0,
0x80, 0x41, 0x00, 0xC1, 0x81, 0x40
};
// CRC低位字节表
private static readonly byte[] m_CRCLowOrderByteTable = new byte[]
{
0x00, 0xC0, 0xC1, 0x01, 0xC3, 0x03, 0x02, 0xC2, 0xC6, 0x06,
0x07, 0xC7, 0x05, 0xC5, 0xC4, 0x04, 0xCC, 0x0C, 0x0D, 0xCD,
0x0F, 0xCF, 0xCE, 0x0E, 0x0A, 0xCA, 0xCB, 0x0B, 0xC9, 0x09,
0x08, 0xC8, 0xD8, 0x18, 0x19, 0xD9, 0x1B, 0xDB, 0xDA, 0x1A,
0x1E, 0xDE, 0xDF, 0x1F, 0xDD, 0x1D, 0x1C, 0xDC, 0x14, 0xD4,
0xD5, 0x15, 0xD7, 0x17, 0x16, 0xD6, 0xD2, 0x12, 0x13, 0xD3,
0x11, 0xD1, 0xD0, 0x10, 0xF0, 0x30, 0x31, 0xF1, 0x33, 0xF3,
0xF2, 0x32, 0x36, 0xF6, 0xF7, 0x37, 0xF5, 0x35, 0x34, 0xF4,
0x3C, 0xFC, 0xFD, 0x3D, 0xFF, 0x3F, 0x3E, 0xFE, 0xFA, 0x3A,
0x3B, 0xFB, 0x39, 0xF9, 0xF8, 0x38, 0x28, 0xE8, 0xE9, 0x29,
0xEB, 0x2B, 0x2A, 0xEA, 0xEE, 0x2E, 0x2F, 0xEF, 0x2D, 0xED,
0xEC, 0x2C, 0xE4, 0x24, 0x25, 0xE5, 0x27, 0xE7, 0xE6, 0x26,
0x22, 0xE2, 0xE3, 0x23, 0xE1, 0x21, 0x20, 0xE0, 0xA0, 0x60,
0x61, 0xA1, 0x63, 0xA3, 0xA2, 0x62, 0x66, 0xA6, 0xA7, 0x67,
0xA5, 0x65, 0x64, 0xA4, 0x6C, 0xAC, 0xAD, 0x6D, 0xAF, 0x6F,
0x6E, 0xAE, 0xAA, 0x6A, 0x6B, 0xAB, 0x69, 0xA9, 0xA8, 0x68,
0x78, 0xB8, 0xB9, 0x79, 0xBB, 0x7B, 0x7A, 0xBA, 0xBE, 0x7E,
0x7F, 0xBF, 0x7D, 0xBD, 0xBC, 0x7C, 0xB4, 0x74, 0x75, 0xB5,
0x77, 0xB7, 0xB6, 0x76, 0x72, 0xB2, 0xB3, 0x73, 0xB1, 0x71,
0x70, 0xB0, 0x50, 0x90, 0x91, 0x51, 0x93, 0x53, 0x52, 0x92,
0x96, 0x56, 0x57, 0x97, 0x55, 0x95, 0x94, 0x54, 0x9C, 0x5C,
0x5D, 0x9D, 0x5F, 0x9F, 0x9E, 0x5E, 0x5A, 0x9A, 0x9B, 0x5B,
0x99, 0x59, 0x58, 0x98, 0x88, 0x48, 0x49, 0x89, 0x4B, 0x8B,
0x8A, 0x4A, 0x4E, 0x8E, 0x8F, 0x4F, 0x8D, 0x4D, 0x4C, 0x8C,
0x44, 0x84, 0x85, 0x45, 0x87, 0x47, 0x46, 0x86, 0x82, 0x42,
0x43, 0x83, 0x41, 0x81, 0x80, 0x40
};
/// <summary>
/// 获得CRC16效验码
/// </summary>
/// <param name="buffer"></param>
/// <returns></returns>
public static void CalculateCrc16(byte[] buffer, out byte crcHighOrderByte, out byte crcLowOrderByte)
{
// CRC高位字节/低位字节
crcHighOrderByte = crcLowOrderByte = 0xff;
for (int i = 0; i < buffer.Length; i++)
{
// 计算CRC查找索引
int crcIndex = crcHighOrderByte ^ buffer[i];
crcHighOrderByte = (byte)(crcLowOrderByte ^ m_CRCHighOrderByteTable[crcIndex]);
crcLowOrderByte = (byte)m_CRCLowOrderByteTable[crcIndex];
}
}
}
}
// 工具函数如上