水情遥测系统中快速差错校验的软件方法

House

    水情遥测系统中快速差错校验的软件方法

0 B6 Q. F6 p) s" C    摘要：讨论了在自报式水情无线遥测系统数据通信中进行快速差错校验的必要性,给出了经过实验检验、可行的软件快速校验方法,并比较了它们的优劣与适应的场合。
7 t/ r$ a+ H- ^& j. W" |$ w    关键词：CRC 汉明码 水情遥测 数据通信 差错校验
将测站的实时水情数据（水位、闸位、雨量等）准确无误地发送到中心站,提供给水文洪水预报、洪水调度、防洪排涝决策等高一级系统,是水情遥测系统最基本、最重要的功能之一。水情遥测系统是一个软硬件综合系统。其基本工作流程是：在测量端（测站）完成水情参数的采集与处理（信源编码、存储记录、信道编码等）,然后将处理过的数据通过无线或有线信道直接或经中继发送至远端的中心站,由中心站进行接收解码并作进一步处理。图1为水情无线遥测系统结构示意图。测站和中继站的主控设备一般采用单片机,用汇编语言编程;而中心接收端主机一般采用微型机,用C语言（或其它高级语言）编程。
$ o/ ]$ e: W; r: {/ u2 Q水情遥测系统的报汛方式一般有三种：定时自报式、查询-应答多和混合式（自报和查询-应答兼容）。三种报汛方式兼有长短。
+ a/ q3 @; z5 H. k$ y

综合考虑系统功耗、可靠性、复杂性等要素,定时自报方式在水情遥测系统中仍占主流。其优点是：功耗极低（值守状态<50μA@+12V,大多以蓄电池供电）,系统结构简单,可靠性较高;缺点是无法实现反馈重发、反传校验等差错控制。显然,需要选择一种合适的有差错校验方法。
# H+ t0 M1 _2 i3 o水情遥测系统的数据通信方式可分为超短波通信、微波通信、卫星通信、移动通信、有线通信等。因其遥测站点常建于交通不便、供电及有线通信条件不足的地区,基于建设成本及运行费用等考虑,数据通信仍以无线超短波通信为主要方式。
9 ]* `! \& s& D2 c/ e5 X/ x在无线数据通信过程中,即使信道质量良好,但由于信号衰减、失真,特别是某些突发性的干扰（如雷电,电磁辐射）不可避免地会发生数据误传,即误码。根据水情遥测系统的相关规范,超短波数据传输的误码率应小于10 -4,以及在每个数据收集周期平均应有90%以上测站（重点控制站必须包括在内）能准确传送数据至中心站。因此采取适当的差错控制方法,提高数据传输的可靠性很有必要。常用的差错控制方法分软件和硬件方式。最简单的是由器件直接实现奇偶校验方式,它占用10%的时间,只检出奇数个位出错。据检测,在电话网中以1200波特率传输数据时,若采用奇偶校验方式,仍会有40%的错误不能检出,这对水情遥测显示是不够的。欲对包括中继在内的每一个站实行码校验,还要求所选校验方式具有高检出率、速度快、编码简单等特点。常见的方式有汉明码、循环冗余校验（CRC）等,虽然这些校验方式也可由硬件实现,但人们角倾向于采用简单经济又具灵活性的软件校验。以下结合工作实际给出经验证可行的快速校验方式,并比较了它们的优劣。文中所涉及到的程序算法均以C语言的形式给出,而将其转变成单片机的算法也不难。
, [9 L) v! {4 X! \$ V) O1 CRC校验
/ c7 N, E; _, D2 E5 U/ uCRC（Cyclical Redundancy Check）校验,又称循环冗余校验,具有极强的检错能力（不能纠错）,算法简单。早期用硬件电路直接搭成,但软件方法成本更低,实现更简单,运算速度也很快。16位的CRC检错率如表1所示[1]。
0 e) a9 o& T* n
: r. o& B6 D+ d- b表1 16位的CRC检错率
6 ~5 f9 c5 q3 O  g0 {4 r

单位个位错误	双位错误	奇数个位错误	比16位短的突出性错误	恰好17位的突发性错误	其他所有突发性错误
全部	全部	全部	全部	99.9969%	99.9984%

常用的16位CRC多项式有两种：一种是CRC-CCITT标准,在微机通信的XMODEM协议中得到了应用;另一种是CRC-16标准,它实际捕获错误的能力不如CRC-CCITT,在IBM的二进制同步协议（BYSYNC）的数据传送中应用已久。两者采用的多项式如表2所示,本文采用前者。

表2 常用的16位CRC多项式

生成多项式的值（genpoly）	本原多项式表示	标  准
1021H F005H	X 16+X 12+X 5+1 X 16+X 15+X 2+1	CRC-CCITT 0 b% K7 u! s7 E& z* ~1 E& vCRC-16

注：genpoly为generator polynomial的合成调,在程序中用作“生成多项式”寄存器。
1.1 直接模2除法CRC实现方式
对16位的CRC而言,用信息段作被除数,生成多项式（本文1021H,CCITT标准）作除数,进行模2除法所产生的余数（2字节）即为CRC校验值,且CRC校验只间余数而不管商是多少。发送时将校验值连在信息段的后面一起发送。在接收端,接收方只需把接收到的CRC校验值连同信息一,作为新的信息段并对其进行相同的CRC运算（只比发送时多2字节）。若得到的新余数（校验值）为0,则表明接收到的信息段和CRC都无差错;反之,说明信息段或CRC有错,应做相应处理。所以CRC的编码和译码并没有本质的区别。程序如下：USHORT crc(USHORT data,USHORT genpoly,USHORT accum)
4 T' i: U3 p, k2 l6 c{// data：数据,所用信息字的一个字节;genpoly:CRC多项式,如1021H;accum：累加器的值,一次赋0,以后放每次校验结果。
  i( M1 X1 q' ?, B  Rdata<<=8; //信息字节左移到高字节
9 j, S3 z% Y+ c; E5 Pfor(int i=8;i>0;i--){
& K2 C0 e' ^4 ?2 a3 B8 Aif((data^accum)&0x800) //如果（data异或accum）的较高位是1
2 R. r. @  N9 X' c% W" s8 h5 Z0 Baccum=(accum<<1)^genpoly; //移位与genpoly异域
else accum<<=1; //否则仅移位
, a6 u& ]/ o% Z1 F! zdata<<=1; //将信息字的下一位升格
$ y2 `2 r8 A$ L7 _9 y8 `}
/ E% z+ q3 k4 A& Zreturn accum; //返回用作下一个信息字校验的累加器值
}
1.2 快速CRC实现方式
直接模2除法CRC方式虽编程简单,但效率不高。采用方式方式,要使用16位的多项式及两字节的累加器,对每一信息位（bit）累加器都要移位一次,再根据移位结果判断是否作异式;每一字节重复8位,运算速度相对较慢,不符合计算机按比特进行计算的规律。但如果采用微机通信中XMODEM协议所使用的CRC查询方式,则比直接CRC模2除法方式快4～10倍。查询方法实施过程：首先用信息字节与累加器的高字节进行异或,并将其结果作初始累加器为0的CRC;然后与原累加器的低字节再作一次异或。一步只有256个样式,可以构造一个256个双字节的查询表,一步实现。这样对每一字节只要作两次操作就可完成。以下是具体步骤。
; D6 @, w) ?* B1 A) ]. ]' f, F（1）构造查询表,运行直接模2除法CRC函数CRC（i,1021,0）,用I从0～255代入,将结果按序排列可得到一个256个样式的双字节查询表。该表只作一次,可以先用C语言微型机上作好,然后再移到单片机上,留作以后查询使用。
（2）取一个双字节累加器accum,赋初值0,将信息流的一个字节赋给另一双字节变量data(accum和data都是双字节变量,以下步骤也是作双字节运算)。
. R- y1 [; p; j8 _1 v* ?) D（3）将accum>>8（也即取原累加器accum的高字节）的值与信息字data相异或,所得结果（是一个<256的值）查上述构造好的查询表,得到一个16位的暂存值。
# `" i2 _$ l9 \% y2 ~) {4 C（4）将accum<<8（即原累加器accum的低字节左移成高字节,低位补0）,与上一步得到的暂存值（16位的值）相异或,结果作为新的累加器值,赋值给accum。
（5）取信息流的下一字节赋给data,重复进行第（3）步和第（4）步,直至所有的信息字节寝用完为止,较后累加器的值就是余数。
, Y: i: _; k* c+ K( ?2 a8 m2 扩展汉明码
2.1 编码方法
CRC校验只能检错但不能纠错。而1949年提出的汉明码是一种能纠正单个错误的线性分组码。其中,既是线性分组码同时也是循环码的（7,4）码有两种。其生成矩阵和校验矩阵分列如下：

5 H4 [% H- ]0 U# R
: `+ o4 i+ H, o
) m7 F$ ~8 N( f4 q
两者使用效果等价。
汉明码是纠正单个错误的完备码,所有的接收码都可对应到一个信息（多一对应）,要么是正确信息,要么是发生单个错误的情形。当有两个错误时,会把它当成另一个码的单个错误加以纠正,导致误码。
* G  R+ l+ ]2 L. k9 J1 k( u, `/ V扩展汉明码在此基础上引入一个校验和,即在编码在时候增加第8位偶校验位,构成（8,4）线性分组码,因而可以纠正一位错误同时检出两位错误。事件上,在发生错误时就是这个偶校验位确定了是错一位还是错两位。若错一位则可以纠正,错两位就只能检出但不能纠正。
$ i( d4 G: H/ J0 \" ]编、译码均以扩展汉明码（8,4）线性分组码为例。为了方便单片机的运算,实现快速编码,可采用查询法。因为信息是一个4位的矢量,记作C,共有16个可能值。为了构成8位发射码矢量,可以建立16个一字节的查询作为8位的发送码。以生成矩阵G1为例,用信息矢量C乘以成矩阵G1再加上一位偶校验就得到了生成码（发送码）。查询表为：
: R! R4 r+ v1 m" _

信息	生成码	信息	生成码	信息	生成三	信息	生成码
0000   {) B$ m( d; z: v6 T0001 # F+ p( |0 d/ V2 y6 S( J) X0010 0011	00000000 ; n) t) [0 [* `00010111 9 L; @. o# _* P7 M- I# @00101101 ) e$ L* v) l, `9 O5 Q& r00111010	0100 ' z( K* h/ V; m* |0 @$ [! q: {3 F) A0101 & S1 v, Z# k' |/ ]0110 - P+ P- Z- Y9 D7 b2 F; s! F; D1 L0111	01001110 ! X5 _/ O, W* }& @01011001 01100011 01110100	1000 1001 . }* o' r9 \, u' ~0 o  L0 Q* ]1010 + m6 I; u9 B! H1011	10001011 10011100 10100110 2 w! Q, G# _! G3 r+ w! Y7 }10110001	1100 # z9 R( `! N  f) F3 N1001 1110 8 |& r+ t8 O4 w2 I( o1111	11000101 * w* f3 @) z9 c( @# |. P10011100 2 r) k7 _  h+ a- |11101000 ( C  _* |0 n2 d11111111

2.2 译码方法
用查询法对（8,4）码进行译码,需要建造有256个值的查询表。按照译码编写查询表。先定出扩展汉明码的校验矩阵,实际上就是将原校验矩阵H1扩展,记为H1,

+ S- x9 \4 P. K, j

: `5 G. m( Z* W1 k' ]  {9 C% {对于作一8位的接收码矢量R,进行RH1 T运算,得到一个4位的伴随矢量,再按如下步骤比较确定原信息。
5 ]4 G! K+ ^; x* `6 E9 r（1）如果伴随式矢量是全0矢量,接收码是正确的,码的前（低）4位就是信息。
$ s, p% @- f; x! b: V, ]7 t! H（2）如果伴随矢量的较后一位是1,则有一位错,可纠正。将伴随矢量与矩 阵H1的每一列相比较,找出相同的那一列,记下列号,再将接收码与该列号相对应的那一位变号（1变0,0变1）,得到的码就是纠正后的原码,信息取码的前（低）4位。
（3）否则,是一位以上的错码且不能纠正。
将一个字节可能出现的所有0～255个可能值值都按上面的译码步骤做一遍得到查询表,留作译码用。另外译码和编码还可以对整个信息字节作一字节的垂直校验以增强校验能力。
上述检验方式已在江苏、宁夏、福建等地的实际工作中得到了验证。CRC校验虽不具备纠错功能但有很高的检错率,应用面也很广。其中,直接模2除法CRC方式因编程简单、占用程序空间少（不用查询表）,适合于数据通信量不大且程序及内存空间有限的场合,反之可选用快速CRC方式。在对数据完整性要求高的场合,可根据具体情况考虑使用汉明码呀扩展汉明码。某些要求更高的特殊情况下,则可选用更复杂一些的校验码,同时通信条件的好坏也是影响校验方式选用的因素之一。
) a* W7 u9 G) G: U2 A, \& x
; m$ T* U5 U# l
8 n" @0 g& M$ g* _+ c+ M. E

sunygd

谢谢分享