Mathematics of CRC

From Wikipedia, the free encyclopedia

Cyclic Redundancy Check (CRC) is based on division in the ring of polynomials over the finite field GF(2) (the integers modulo 2), that is, the set of polynomials where each coefficient is either zero or one, and arithmetic operations wrap around (due to the nature of binary arithmetic).

Any string of bits can be interpreted as the coefficients of a message polynomial of this sort, and to find the CRC, we multiply the message polynomial by xn and then find the remainder when dividing by the degree-n generator polynomial. The coefficients of the remainder polynomial are the bits of the CRC.

In general form:

M(x) \cdot x^{n} = Q(x) \cdot G(x) + R (x)

Here M(x) is the original message polynomial and G(x) is the degree-n generator polynomial. The bits of M(x) \cdot x^{n} are the original message with n zeroes added at the end. R(x) is the remainder polynomial, which is the CRC 'checksum'. The quotient polynomial Q(x) is of no interest.

In communication, the sender attaches the n bits of R after the original message bits of M, sending out M(x) \cdot x^{n} - R(x) (the codeword.) The receiver, knowing G(x) and therefore n, separates M from R and repeats the calculation, verifying that the received and computed R are equal. If they are, then the receiver assumes the received message bits are correct.

In practice CRC calculations resemble long division in binary, except that the subtractions involved do not borrow from more significant digits, and thus become exclusive or operations.

A CRC is a checksum in a strict mathematical sense, as it can be expressed as the weighted modulo-2 sum of per-bit syndromes, but that word is generally reserved more specifically for sums computed using larger moduli, such as 10, 256, or 65535.

CRCs can also be used as part of error-correcting codes, which allow not only the detection of transmission errors, but the reconstruction of the correct message. These codes are based on closely related mathematical principles.

Contents

[edit] Polynomial arithmetic modulo 2

Since the coefficients are constrained to a single bit, any math operation on CRC polynomials must map the coefficients of the result to either zero or one. For example in addition:

(x^3 + x) + (x + 1) = x^3 + 2x + 1 \equiv x^3 + 1 \pmod 2

Note that 2x becomes zero in the above equation because addition of coefficients is performed modulo 2:

2x = x + x = x\times(1 + 1) = x\times0 \equiv 0 \pmod 2

Multiplication is similar:

(x^2 + x)(x + 1) = x^3 + 2x^2 + x \equiv x^3 + x \pmod 2

We can also divide polynomials mod 2 and find the quotient and remainder. For example, suppose we're dividing x³ + x² + x by x + 1. We would find that

\frac{(x^3 + x^2 + x)}{(x+1)} = (x^2 + 1) - \frac{1}{(x+1)}

In other words,

(x^3 + x^2 + x) = (x^2 + 1)(x + 1) - 1 \equiv (x^2 + 1)(x + 1) + 1 \pmod 2

The division yields a quotient of x² + 1 with a remainder of -1, which, since it is odd, has a last bit of 1.

In the above equations, x2 + x + 1 represents the original message bits 111, x + 1 is the generator polynomial, and the remainder 1 (equivalently, x0) is the CRC. The degree of the generator polynomial is 1, so we first multiplied the message by x1 to get x3 + x2 + x.

[edit] Variations

There are several standard variations on CRCs, any or all of which may be used with any CRC polynomial. Implementation variations such as endianness and CRC presentation only affect the mapping of bit strings to the coefficients of M(x) and R(x), and do not impact the properties of the algorithm.

  • To check the CRC, instead of calculating the CRC on the message and comparing it to the CRC, a CRC calculation may be run on the entire codeword. If the result is zero, the check passes. This works because the codeword is M(x) \cdot x^{n} - R(x) = Q(x) \cdot G(x), which is always divisible by G(x).
This simplifies many implementations by avoiding the need to treat the last few bytes of the message specially when checking CRCs.
  • The shift register may be initialized with ones instead of zeroes. This is equivalent to inverting the first n bits of the message feeding them into the algorithm. The CRC equation becomes M(x) \cdot x^{n} + \sum_{i=m}^{m+n-1} x^{i} = Q(x) \cdot G(x) + R (x), where m > deg(M(x)) is the length of the message in bits. The change this imposes on R(x) is a function of the generating polynomial and the message length, \left ( \sum_{i=m}^{m+n-1} x^{i} \right ) \bmod G(x).
The reason this method is used is because an unmodified CRC does not distinguish between two messages which differ only in the number of leading zeroes, because leading zeroes do not affect the value of M(x). When this inversion is done, the CRC does distinguish between such messages.
  • The CRC may be inverted before being appended to the message stream. While an unmodified CRC distinguishes between messages M(x) with different numbers of trailing zeroes, it does not detect trailing zeroes appended after ther CRC remainder itself. This is because all valid codewords are multiples of G(x), so x times that codeword is also a multiple. (In fact, this is precisely why the first variant described above works.)

In practice, the last two variations are invariably used together. They change the transmitted CRC, so must be implemented at both the transmitter and the receiver. While presetting the shift register to ones is straightforward to do at both ends, inverting affects receivers implementing the first varaition, because the CRC of a full codeword that already includes a CRC is no longer zero. Instead, it is a fixed non-zero pattern, the CRC of the inversion pattern of n ones.

The CRC thus may be checked either by the obvious method of computing the CRC on the message, inverting it, and comparing with the CRC in the message stream, or by calculating the CRC on the entire codeword and comparing it with an expected fixed value called the check polynomial C(x).

This result is sometimes called the corresponding magic number, and may be computed as C(x) = \left ( \sum_{i=n}^{2n-1} x^{i} \right ) \bmod G(x), or equivalently by computing the unmodified CRC of a message consisting of n ones, M(x) = \sum_{i=0}^{n-1} x^{i}

These inversions are extremely common, so much so that they may not be mentioned at all in a brief description of a communication protocol and simply assumed. In particular, they are by convention always performed when using the CRC-32 or CRC-16-CCITT polynomials.

[edit] Reversed representations and reciprocal polynomials

[edit] Polynomial representations

All the well-known CRC generator polynomials of degree n have two common hexadecimal representations. In both cases, the coefficient of xn is omitted and understood to be 1.

  • The msbit-first representation is a hexadecimal number with n bits, the least significant bit of which is always 1. The most significant bit represents the coefficient of xn − 1 and the least significant bit represents the coefficient of x0.
  • The lsbit-first representation is a hexadecimal number with n bits, the most significant bit of which is always 1. The most significant bit represents the coefficient of x0 and the least significant bit represents the coefficient of xn − 1.

The msbit-first form is often referred to in the literature as the normal representation, while the lsbit-first is called the reversed representation. It is essential to use the correct form when implementing a CRC. If the coefficient of xn − 1 happens to be zero, the forms can be distinguished at a glance by seeing which end has the bit set.

To further confuse the matter, the excellent paper by P. Koopman and T. Chakravarty [1][2] converts CRC generator polynomials to hexadecimal numbers in yet another way: msbit-first, but including the xn coefficient and omitting the x0 coefficient. This "Koopman" representation has the advantage that the degree can be determined from the hexadecimal form and the coefficients are easy to read off in left-to-right order. However, it is not used anywhere else and is not recommended due to the risk of confusion.

[edit] Reciprocal polynomials

A reciprocal polynomial is created by assigning the xn through x0 coefficients of one polynomial to the x0 through xn coefficients of a new polynomial. That is, the reciprocal of the n bit polynomial G(x) is xnG(x − 1). Example: The reverse of x16 + x15 + x2 + 1 is x16 + x14 + x1 + 1. The most interesting property of reciprocal polynomials, when used in CRCs, is that they have the exact same error-detecting strength as the polynomials they are reciprocals of. The reciprocal of a polynomial generates the same codewords, only bit reversed — that is, if all but the first n bits of a codeword under the original polynomial are taken, reversed and used as a new message, the CRC of that message under the reciprocal polynomial equals the reverse of the first n bits of the original codeword. But the reciprocal polynomial is not the same as the original polynomial, and the CRCs generated using it are not the same (even modulo bit reversal) as those generated by the original polynomial.

Another interesting property of reciprocal polynomials is related to their representation. Consider again the CRC-16-IBM polynomial x16 + x15 + x2 + 1 and its reciprocal. These are its representations

  • Normal original: 0x8005
  • Reversed original: 0xA001
  • Koopman original: 0xC002
  • Normal reciprocal: 0x4003
  • Reversed reciprocal: 0xC002
  • Koopman reciprocal: 0xA001

Note that the Koopman representation of the reciprocal polynomial is the same as the reversed representation of the original polynomial. This means that if you mistakenly use the Koopman representation of a polynomial in a right-shift algorithm, you get a CRC that is just as strong as (but different from) the one you intended. This holds only for polynomials with a degree which is a multiple of 4.[citation needed]

[edit] Error detection strength

The error-detection ability of a CRC depends on the degree of its key polynomial and on the specific key polynomial used. The "error polynomial" E(x) is the symmetric difference of the received message codeword and the correct message codeword. An error will go undetected by a CRC algorithm if and only if the error polynomial is divisible by the CRC polynomial.

  • Because a CRC is based on division, no polynomial can detect errors consisting of a string of zeroes prepended to the data, or of missing leading zeroes. However, see Variations.
  • All single bit errors will be detected by any polynomial with at least two terms with non-zero coefficients. The error polynomial is xk, and xk is divisible only by polynomials xi where i \le k.
  • All two bit errors separated by a distance less than the order of the primitive polynomial which is a factor of the generator polynomial will be detected. The error polynomial in the two bit case is E(x) = x^i + x^k = x^k \cdot (x^{i-k} + 1), \; i > k. As noted above, the xk term will not be divisible by the CRC polynomial, which leaves the xi − k + 1 term. By definition, the smallest value of i − k such that a polynomial divides xi − k + 1 is the polynomial's order or exponent. The polynomials with the largest order are called primitive polynomials, and for polynomials of degree n with binary coefficients, have order 2n − 1.
  • All errors in an odd number of bits will be detected by a polynomial which is a multiple of x + 1. This is equivalent to the polynomial having an even number of terms with non-zero coefficients. This capacity assumes that the generator polynomial is the product of x + 1 and a primitive polynomial of degree n − i since all primitive polynomials except x + 1 have an odd number of non-zero coefficients.
  • All burst errors of length n will be detected by any polynomial of degree n or greater which has a non-zero x0 term.

(as an aside, there is never reason to use a polynomial with a zero x0 term. Recall that a CRC is the remainder of the message polynomial times x^n divided by the CRC polynomial. A polynomial with a zero x0 term always has x as a factor. So if K(x) is the original CRC polynomial and K(x) = x \cdot K'(x), then

M(x) \cdot x^{n-1} = Q(x) \cdot K'(x) + R(x)

M(x) \cdot x^n = Q(x) \cdot x \cdot K'(x^n) + x \cdot R(x)

M(x) \cdot x^n = Q(x) \cdot K(x^n) + x \cdot R(x)

That is, the CRC of any message with the K(x) polynomial is the same as that of the same message with the K'(x) polynomial with a zero appended. It is just a waste of a bit.)

The combination of these factors means that good CRC polynomials are often primitive polynomials (which have the best 2-bit error detection) or polynomials whose factors are a primitive polynomial of degree n − 1 and x + 1 (which detects all odd numbers of bit errors, and has half the two-bit error detection ability of a primitive polynomial of degree n)[1].

[edit] Bitfilters

Analysis Technique using bitfilters[1] allows to determine quite more efficient the properties of a given generator polynomial. The results are the following:

  1. All burst errors (but one) with length no longer than the generator polynomial can be detected by any generator polynomial 1 + ... + Xn. This includes 1-bit errors (burst of length 1). The maximum length is n + 1, when n is the degree of the generator polynomial (which itself has a length of n + 1). The exception to this result is a bit pattern the same as that of the generator polynomial.
  2. All uneven bit errors are detected by generator polynomials with even number of terms.
  3. 2-Bit errors in a (multiple) distance of the longest bitfilter of even parity to a generator polynomial are not detected; all others are detected. For degrees up to 32 there is an optimal generator polynomial with that degree and even number of terms; in this case the period mentioned above is 2n − 1 − 1. For n = 16 this means that blocks of 32767 bits length do not contain undiscovered 2-Bit errors. For uneven number of terms in the generator polynomial there can be a period of 2n − 1; however, these generator polynomials (with odd number of terms) do not discover all odd number of errors, so they should be avoided. A list of the corresponding generators with even number of terms can be found in the link mentioned at the beginning of this section.
  4. All single bit errors within the bitfilter period mentioned above (for even terms in the generator polynomial) can be identified uniquely by their residual. So CRC method can be used for error correction as well (within those limits, e.g. 32767 bits with optimal generator polynomials of degree 16). Since all odd errors leave an odd residual, all even an even residual, 1-bit errors and 2-bit errors can be distinguished, however not necessarily 1-bit errors and 3-bit errors. Thus bit error correction might be erroneous itself and produce more errors.

[edit] See also

General category

Specific Technological References

[edit] References

  1. ^ a b c Koopman, P. (June 2002). "32-Bit Cyclic Redundancy Codes for Internet Applications". The International Conference on Dependable Systems and Networks: 459. doi:10.1109/DSN.2002.1028931.  - verification of Castagnoli's results by exhaustive search and some new good polynomials
  2. ^ Koopman, P. and Chakravarty, T. (2004). "Cyclic Redundancy Code (CRC) Polynomial Selection For Embedded Networks". The International Conference on Dependable Systems and Networks: 131. doi:10.1109/DSN.2004.1311885.  - analysis of short CRC polynomials for embedded applications