View source for Burst error-correcting code

{{Short description|Codes intended to correct short, contiguous errors in a communications channel}}
In [[coding theory]], '''burst error-correcting codes''' employ methods of correcting [[burst error]]s, which are errors that occur in many consecutive bits rather than occurring in bits independently of each other.

Many codes have been designed to correct [[random error]]s. Sometimes, however, channels may introduce errors which are localized in a short interval. Such errors occur in a burst (called [[burst error]]s) because they occur in many consecutive bits. Examples of burst errors can be found extensively in storage mediums. These errors may be due to physical damage such as scratch on a disc or a stroke of lightning in case of wireless channels. They are not independent; they tend to be spatially concentrated. If one bit has an error, it is likely that the adjacent bits could also be corrupted. The methods used to correct random errors are inefficient to correct burst errors.

==Definitions==
[[File:Burst Error,Length 5,Description.jpg|thumb|400px|A burst of length 5]]
'''A burst of length {{math|ℓ}}<ref name = "Coding Bounds">{{cite arXiv |last=Fong |first=W.H. |date=2011|title=Coding Bounds for Multiple Phased-Burst Correction and Single Burst Correction Codes |class=cs.IT |eprint=1104.1408}}</ref>'''

Say a codeword <math>C</math> is transmitted, and it is received as <math>Y = C + E.</math> Then, the error vector <math>E</math> is called a burst of length <math>\ell</math> if the nonzero components of <math>E</math> are confined to <math>\ell</math> consecutive components. For example, <math>E = (0\textbf{1000011}0)</math> is a burst of length <math>\ell = 7.</math>

Although this definition is sufficient to describe what a burst error is, the majority of the tools developed for burst error correction rely on cyclic codes. This motivates our next definition.

'''A cyclic burst of length {{math|ℓ}}<ref name = "Coding Bounds" />'''

An error vector <math>E</math> is called a cyclic burst error of length <math>\ell</math> if its nonzero components are confined to <math>\ell</math> cyclically consecutive components. For example, the previously considered error vector <math>E = (010000110)</math>, is a cyclic burst of length <math>\ell = 5</math>, since we consider the error starting at position <math>6</math> and ending at position <math>1</math>. Notice the indices are <math>0</math>-based, that is, the first element is at position <math>0</math>.

For the remainder of this article, we will use the term burst to refer to a cyclic burst, unless noted otherwise.

=== Burst description ===
It is often useful to have a compact definition of a burst error, that encompasses not only its length, but also the pattern, and location of such error. We define a burst description to be a tuple <math>(P,L)</math> where <math>P</math> is the pattern of the error (that is the string of symbols beginning with the first nonzero entry in the error pattern, and ending with the last nonzero symbol), and <math>L</math> is the location, on the codeword, where the burst can be found.<ref name = "Coding Bounds" />

For example, the burst description of the error pattern <math>E = (010000110)</math> is <math>D = (1000011,1)</math>. Notice that such description is not unique, because <math>D' = (11001,6)</math> describes the same burst error. In general, if the number of nonzero components in <math>E</math> is <math>w</math>, then <math>E</math> will have <math>w</math> different burst descriptions each starting at a different nonzero entry of <math>E</math>. To remedy the issues that arise by the ambiguity of burst descriptions with the theorem below, however before doing so we need a definition first.

'''Definition.''' The number of symbols in a given error pattern <math>y,</math> is denoted by <math>\mathrm{length}(y).</math>

{{math theorem | name = Theorem (Uniqueness of burst descriptions) | math_statement = Suppose <math>E</math> is an error vector of length <math>n</math> with two burst descriptions <math>(P_1,L_1)</math> and <math>(P_2,L_2)</math>. If <math>\mathrm{length}(P_1) + \mathrm{length}(P_2) \leqslant n + 1,</math> then the two descriptions are identical that is, their components are equivalent.<ref>{{cite book |first=R.J. |last=McEliece |title=The Theory of Information and Coding |publisher=Cambridge University Press |edition=Student |date=2004 |isbn=978-0-521-83185-7 |pages= |url={{GBurl|bQ9fWH7fh3IC|pg=PR5}}}}</ref>

{{math proof | proof = Let <math>w</math> be the [[hamming weight]] (or the number of nonzero entries) of <math>E</math>. Then <math>E</math> has exactly <math>w</math> error descriptions. For <math>w = 0, 1,</math> there is nothing to prove. So we assume that <math>w \geqslant 2</math> and that the descriptions are not identical. We notice that each nonzero entry of <math>E</math> will appear in the pattern, and so, the components of <math>E</math> not included in the pattern will form a cyclic run of zeros, beginning after the last nonzero entry, and continuing just before the first nonzero entry of the pattern. We call the set of indices corresponding to this run as the zero run. We immediately observe that each burst description has a zero run associated with it and that each zero run is disjoint. Since we have <math>w</math> zero runs, and each is disjoint, we have a total of <math>n-w</math> distinct elements in all the zero runs. On the other hand we have: <math display="block">\begin{align}
n-w = \text{number of zeros in } E &\geqslant (n - \mathrm{length}(P_1)) + (n - \mathrm{length}(P_2))\\ 
&= 2n - \left ( \mathrm{length}(P_1) +  \mathrm{length}(P_2) \right ) \\
&\geqslant 2n - (n+1) && \mathrm{length}(P_1) + \mathrm{length}(P_2) \leqslant n+1 \\
&= n-1
\end{align}</math>
This contradicts <math>w \geqslant 2.</math> Thus, the burst error descriptions are identical.}}}}

A [[corollary]] of the above theorem is that we cannot have two distinct burst descriptions for bursts of length <math>\tfrac{1}{2}(n+1).</math>

==Cyclic codes for burst error correction==

[[Cyclic code]]s are defined as follows: think of the <math>q</math> symbols as elements in <math>\mathbb{F}_q</math>. Now, we can think of words as polynomials over <math>\mathbb{F}_q,</math> where the individual symbols of a word correspond to the different coefficients of the polynomial. To define a cyclic code, we pick a fixed polynomial, called [[generator polynomial]]. The codewords of this cyclic code are all the polynomials that are divisible by this generator polynomial.

Codewords are polynomials of degree <math>\leqslant n-1</math>. Suppose that the generator polynomial <math>g(x)</math> has degree <math>r</math>. Polynomials of degree <math>\leqslant n-1</math> that are divisible by <math>g(x)</math> result from multiplying <math>g(x)</math> by polynomials of degree <math>\leqslant n-1-r</math>. We have <math>q^{n-r}</math> such polynomials. Each one of them corresponds to a codeword. Therefore, <math>k=n-r</math> for cyclic codes.

Cyclic codes can detect all bursts of length up to <math>\ell = n-k = r</math>. We will see later that the burst error detection ability of any <math>(n, k)</math> code is bounded from above by <math>\ell \leqslant n-k</math>. Cyclic codes are considered optimal for burst error detection since they meet this upper bound:

{{math theorem | name = Theorem (Cyclic burst correction capability) | math_statement = Every cyclic code with generator polynomial of degree <math>r</math> can detect all bursts of length <math>\leqslant r.</math>

{{math proof | proof = We need to prove that if you add a burst of length <math>\leqslant r</math> to a codeword (i.e. to a polynomial that is divisible by <math>g(x)</math>), then the result is not going to be a codeword (i.e. the corresponding polynomial is not divisible by <math>g(x)</math>). It suffices to show that no burst of length <math>\leqslant r</math> is divisible by <math>g(x)</math>. Such a burst has the form <math>x^i b(x)</math>, where <math>\deg(b(x)) < r.</math> Therefore, <math>b(x)</math> is not divisible by <math>g(x)</math> (because the latter has degree <math>r</math>). <math>g(x)</math> is not divisible by <math>x</math> (Otherwise, all codewords would start with <math>0</math>). Therefore, <math>x^i</math> is not divisible by <math>g(x)</math> as well.}}
}}
The above proof suggests a simple algorithm for burst error detection/correction in cyclic codes: given a transmitted word (i.e. a polynomial of degree <math>\leqslant n-1</math>), compute the remainder of this word when divided by <math>g(x)</math>. If the remainder is zero (i.e. if the word is divisible by <math>g(x)</math>), then it is a valid codeword. Otherwise, report an error. To correct this error, subtract this remainder from the transmitted word. The subtraction result is going to be divisible by <math>g(x)</math> (i.e. it is going to be a valid codeword).

By the upper bound on burst error detection (<math>\ell \leqslant n-k = r</math>), we know that a cyclic code can not detect ''all'' bursts of length <math>\ell > r</math>. However cyclic codes can indeed detect ''most'' bursts of length <math>> r</math>. The reason is that detection fails only when the burst is divisible by <math>g(x)</math>. Over binary alphabets, there exist <math>2^{\ell-2}</math> bursts of length <math>\ell</math>. Out of those, only <math>2^{\ell-2-r}</math> are divisible by <math>g(x)</math>. Therefore, the detection failure probability is very small (<math>2^{-r}</math>) assuming a uniform distribution over all bursts of length <math>\ell</math>.

We now consider a fundamental theorem about cyclic codes that will aid in designing efficient burst-error correcting codes, by categorizing bursts into different cosets.

{{math theorem | name = Theorem (Distinct [[Coset]]s) | math_statement = A linear code <math>C</math> is an <math>\ell</math>-burst-error-correcting code if all the burst errors of length <math>\leqslant \ell</math> lie in distinct [[coset]]s of <math>C</math>.

{{math proof | proof = Let <math>\mathbf{e}_1, \mathbf{e}_2</math> be distinct burst errors of length <math> \leqslant \ell</math> which lie in same coset of code <math>C</math>. Then <math>\mathbf{c} = \mathbf{e}_1 - \mathbf{e}_2</math> is a codeword. Hence, if we receive <math>\mathbf{e}_1,</math> we can decode it either to <math>\mathbf{0}</math> or <math>\mathbf{c}</math>. In contrast, if all the burst errors <math>\mathbf{e}_1</math> and <math>\mathbf{e}_2</math> do not lie in same coset, then each burst error is determined by its syndrome. The error can then be corrected through its syndrome. Thus, a linear code <math>C</math> is an <math>\ell</math>-burst-error-correcting code if and only if all the burst errors of length <math>\leqslant \ell</math> lie in distinct cosets of <math>C</math>.}}}}

{{math theorem | name = Theorem (Burst error codeword classification) | math_statement = Let <math>C</math> be a linear <math>\ell</math>-burst-error-correcting code. Then no nonzero burst of length <math>\leqslant 2\ell</math> can be a codeword.

{{math proof | proof = Let <math>c</math> be a codeword with a burst of length <math>\leqslant 2\ell</math>. Thus it has the pattern <math>(0,1,u,v,1,0)</math>, where <math>u</math> and <math>v</math> are words of length <math>\leqslant \ell -1.</math> Hence, the words <math>w = (0,1,u,0,0,0)</math> and <math>c-w = (0,0,0,v,1,0)</math> are two bursts of length <math> \leqslant \ell</math>. For binary linear codes, they belong to the same coset. This contradicts the Distinct Cosets Theorem, therefore no nonzero burst of length <math>\leqslant 2\ell</math> can be a codeword.}}
}}
==Burst error correction bounds==

===Upper bounds on burst error detection and correction===

By upper bound, we mean a limit on our error detection ability that we can never go beyond. Suppose that we want to design an <math>(n, k)</math> code that can detect all burst errors of length <math>\leqslant \ell.</math> A natural question to ask is: given <math>n</math> and <math>k</math>, what is the maximum <math>\ell</math> that we can never achieve beyond? In other words, what is the upper bound on the length <math>\ell</math> of bursts that we can detect using any <math>(n, k)</math> code? The following theorem provides an answer to this question.

{{math theorem | name = Theorem (Burst error detection ability) | math_statement = The burst error detection ability of any <math>(n, k)</math> code is <math>\ell \leqslant  n-k.</math>

{{math proof | proof = First we observe that a code can detect all bursts of length <math>\leqslant  \ell</math> if and only if no two codewords differ by a burst of length <math>\leqslant  \ell</math>. Suppose that we have two code words <math>\mathbf{c}_1</math> and <math>\mathbf{c}_2</math> that differ by a burst <math>\mathbf{b}</math> of length <math>\leqslant  \ell</math>. Upon receiving <math>\mathbf{c}_1</math>, we can not tell whether the transmitted word is indeed <math>\mathbf{c}_1</math> with no transmission errors, or whether it is <math>\mathbf{c}_2</math> with a burst error <math>\mathbf{b}</math> that occurred during transmission. Now, suppose that every two codewords differ by more than a burst of length <math>\ell.</math> Even if the transmitted codeword <math>\mathbf{c}_1</math> is hit by a burst <math>\mathbf{b}</math> of length <math>\ell</math>, it is not going to change into another valid codeword. Upon receiving it, we can tell that this is <math>\mathbf{c}_1</math> with a burst <math>\mathbf{b}.</math> By the above observation, we know that no two codewords can share the first <math>n-\ell</math> symbols. The reason is that even if they differ in all the other <math>\ell</math> symbols, they are still going to be different by a burst of length <math>\ell.</math> Therefore, the number of codewords <math>q^k</math> satisfies <math>q^k \leqslant  q^{n-\ell}.</math> Applying <math>\log_q</math> to both sides and rearranging, we can see that <math>\ell \leqslant  n-k</math>.}}
}}
Now, we repeat the same question but for error correction: given <math>n</math> and <math>k</math>, what is the upper bound on the length <math>\ell</math> of bursts that we can correct using any <math>(n, k)</math> code? The following theorem provides a preliminary answer to this question:

{{math theorem | name = Theorem (Burst error correction ability) | math_statement = The burst error correction ability of any <math>(n, k)</math> code satisfies <math>\ell \leqslant  n-k-\log_q (n-\ell)+2</math>

{{math proof | proof = First we observe that a code can correct all bursts of length <math>\leqslant  \ell</math> if and only if no two codewords differ by the sum of two bursts of length <math>\leqslant  \ell.</math> Suppose that two codewords <math>\mathbf{c}_1</math> and <math>\mathbf{c}_2</math> differ by bursts <math>\mathbf{b}_1</math> and <math>\mathbf{b}_2</math> of length <math>\leqslant  \ell</math> each. Upon receiving <math>\mathbf{c}_1</math> hit by a burst <math>\mathbf{b}_1</math>, we could interpret that as if it was <math>\mathbf{c}_2</math> hit by a burst <math>-\mathbf{b}_2</math>. We can not tell whether the transmitted word is <math>\mathbf{c}_1</math> or <math>\mathbf{c}_2</math>. Now, suppose that every two codewords differ by more than two bursts of length <math>\ell</math>. Even if the transmitted codeword <math>\mathbf{c}_1</math> is hit by a burst of length <math>\ell</math>, it is not going to look like another codeword that has been hit by another burst. For each codeword <math>\mathbf{c},</math> let <math>B(\mathbf{c})</math> denote the set of all words that differ from <math>\mathbf{c}</math> by a burst of length <math>\leqslant  \ell.</math> Notice that <math>B(\mathbf{c})</math> includes <math>\mathbf{c}</math> itself. By the above observation, we know that for two different codewords <math>\mathbf{c}_i</math> and <math>\mathbf{c}_j, B(\mathbf{c}_i)</math> and <math>B(\mathbf{c}_j)</math> are disjoint. We have <math>q^k</math> codewords. Therefore, we can say that <math>q^k |B( \mathbf{c})| \leqslant  q^n</math>. Moreover, we have <math>(n-\ell)q^{\ell-2} \leqslant  |B(\mathbf{c})|</math>. By plugging the latter inequality into the former, then taking the base <math>q</math> logarithm and rearranging, we get the above theorem.}}
}}
A stronger result is given by the Rieger bound:

{{math theorem | name = Theorem (Rieger bound) | math_statement = If <math>\ell</math> is the burst error correcting ability of an <math>(n, k)</math> linear block code, then <math>2\ell \leqslant n - k</math>.

{{math proof | proof = Any linear code that can correct any burst pattern of length <math>\leqslant \ell</math> cannot have a burst of length <math>\leqslant 2\ell</math> as a codeword. If it had a burst of length <math>\leqslant 2\ell</math> as a codeword, then a burst of length <math>\ell</math> could change the codeword to a burst pattern of length <math>\ell</math>, which also could be obtained by making a burst error of length <math>\ell</math> in all zero codeword. If vectors are non-zero in first <math>2\ell</math> symbols, then the vectors should be from different subsets of an array so that their difference is not a codeword of bursts of length <math>2\ell</math>. Ensuring this condition, the number of such subsets is at least equal to number of vectors. Thus, the number of subsets would be at least <math>q^{2\ell}</math>. Hence, we have at least <math>2\ell</math> distinct symbols, otherwise, the difference of two such polynomials would be a codeword that is a sum of two bursts of length <math>\leqslant  \ell.</math> Thus, this proves the Rieger Bound.}}}}

'''Definition.''' A linear burst-error-correcting code achieving the above Rieger bound is called an optimal burst-error-correcting code.

=== Further bounds on burst error correction ===

There is more than one upper bound on the achievable code rate of linear block codes for multiple phased-burst correction (MPBC). One such bound is constrained to a maximum correctable cyclic burst length within every subblock, or equivalently a constraint on the minimum error free length or gap within every phased-burst. This bound, when reduced to the special case of a bound for single burst correction, is the Abramson bound (a corollary of the Hamming bound for burst-error correction) when the cyclic burst length is less than half the block length.<ref name = "Ling, San">{{cite book |last1=Ling |first1=San |first2=Chaoping |last2=Xing |title=Coding Theory: A First Course |publisher=Cambridge University Press |date=2004 |isbn=978-0-521-52923-5 |pages= |url={{GBurl|N1jiL8V3ISwC|pg=PR7}}}}</ref>

{{math theorem | name = Theorem (number of bursts) | math_statement = For <math>1 \leqslant \ell \leqslant \tfrac{1}{2}(n+1),</math> over a binary alphabet, there are <math>n2^{\ell-1}+1</math> vectors of length <math>n</math> which are bursts of length <math>\leqslant \ell</math>.<ref name = "Coding Bounds" />

{{math proof | proof = Since the burst length is <math>\leqslant \tfrac{1}{2}(n+1),</math> there is a unique burst description associated with the burst. The burst can begin at any of the <math>n</math> positions of the pattern. Each pattern begins with <math>1</math> and contain a length of <math>\ell</math>. We can think of it as the set of all strings that begin with <math>1</math> and have length <math>\ell</math>. Thus, there are a total of <math>2^{\ell-1}</math> possible such patterns, and a total of <math>n2^{\ell-1}</math> bursts of length <math>\leqslant \ell.</math> If we include the all-zero burst, we have <math>n2^{\ell-1}+1</math> vectors representing bursts of length <math>\leqslant \ell.</math>}}}}

{{math theorem | name = Theorem (Bound on the number of codewords) | math_statement = If <math>1 \leqslant \ell \leqslant \tfrac{1}{2}(n+1),</math> a binary <math>\ell</math>-burst error correcting code has at most <math>2^n/(n2^{\ell-1}+1)</math> codewords.

{{math proof | proof = Since <math>\ell \leqslant \tfrac{1}{2}(n+1)</math>, we know that there are <math>n 2^{\ell-1}+1</math> bursts of length <math>\leqslant \ell</math>. Say the code has <math>M</math> codewords, then there are <math>Mn2^{\ell-1}</math> codewords that differ from a codeword by a burst of length <math>\leqslant \ell</math>. Each of the <math>M</math> words must be distinct, otherwise the code would have distance <math>< 1</math>. Therefore, <math>M(2^{\ell-1}+1) \leqslant 2^n</math> implies <math>M \leqslant 2^n/(n2^{\ell-1}+1).</math>}}}}

{{math theorem | name = Theorem (Abramson's bounds) | math_statement = If <math>1 \leqslant \ell \leqslant \tfrac{1}{2}(n+1)</math> is a binary [[Linear code|linear]] <math>(n,k), \ell</math>-burst error correcting code, its block-length must satisfy: <math display="block">n \leqslant 2^{n-k-\ell+1} -1.</math>

{{math proof | proof = For a linear <math>(n,k)</math> code, there are <math>2^k</math> codewords. By our previous result, we know that 
<math display="block">2^k \leqslant \frac{2^n}{n2^{\ell-1}+1}.</math>
Isolating <math>n</math>, we get <math>n \leqslant 2^{n-k-\ell+1}-2^{-\ell+1}</math>. Since <math>\ell \geqslant 1</math> and <math>n</math> must be an integer, we have <math>n \leqslant 2^{n-k-\ell+1}-1</math>.}}}}

''Remark.'' <math>r = n - k</math> is called the redundancy of the code and in an alternative formulation for the Abramson's bounds is <math>r \geqslant \lceil \log_2(n+1) \rceil + \ell -1.</math>

==Fire codes==
{{See also|Cyclic code#Fire codes as cyclic bounds}}

Sources:<ref name = "Ling, San" /><ref name = "Moon, Todd">{{cite book |first=Todd K. |last=Moon |title=Error Correction Coding: Mathematical Methods and Algorithms |publisher=Wiley |date=2005 |isbn=978-0-471-64800-0 |pages= |url={{GBurl|UMmbDwAAQBAJ|pg=PR13}}}}</ref><ref name = "Lin, Shu">{{cite book |last1=Lin |first1=Shu |first2=Daniel J. |last2=Costello  |title=Error Control Coding: Fundamentals and Applications |publisher=Pearson-Prentice Hall |date=2004 |edition=2nd |isbn=978-0-13-017973-9 |pages= |url=}}</ref>

While [[cyclic code]]s in general are powerful tools for detecting burst errors, we now consider a family of binary cyclic codes named Fire Codes, which possess good single burst error correction capabilities. By single burst, say of length <math>\ell</math>, we mean that all errors that a received codeword possess lie within a fixed span of <math>\ell</math> digits.

Let <math>p(x)</math> be an [[irreducible polynomial]] of degree <math>m</math> over <math>\mathbb{F}_2</math>, and let <math>p</math> be the period of <math>p(x)</math>. The period of <math>p(x)</math>, and indeed of any polynomial, is defined to be the least positive integer <math>r</math> such that <math>p(x) |x^r - 1.</math> Let <math>\ell</math> be a positive integer satisfying <math>\ell \leqslant m</math> and <math>2\ell -1</math> not divisible by <math>p</math>, where <math>m</math> and <math>p</math> are the degree and period of <math>p(x)</math>, respectively. Define the Fire Code <math>G</math> by the following [[generator polynomial]]:
<math display="block">g(x) = \left (x^{2\ell -1} + 1 \right )p(x).</math>

We will show that <math>G</math> is an <math>\ell</math>-burst-error correcting code.

{{math theorem | name = Lemma 1 | math_statement = <math>\gcd \left (p(x), x^{2\ell -1} + 1 \right ) = 1.</math>

{{math proof | proof = Let <math>d(x)</math> be the greatest common divisor of the two polynomials. Since <math>p(x)</math> is irreducible, <math>\deg(d(x)) = 0</math> or <math>\deg(p(x))</math>. Assume <math>\deg(d(x)) \neq 0,</math> then <math>p(x) = c d(x)</math> for some constant <math>c</math>. But, <math>(1/c)p(x)</math> is a divisor of <math>x^{2\ell -1}+1</math> since <math>d(x)</math> is a divisor of <math>x^{2\ell -1}+1</math>. But this contradicts our assumption that <math>p(x)</math> does not divide <math>x^{2\ell -1}+1.</math> Thus, <math>\deg(d(x)) = 0,</math> proving the lemma.}}}}

{{math theorem | name = Lemma 2 | math_statement = If <math>p(x)</math> is a polynomial of period <math>p</math>, then <math>p(x)|x^k-1</math> if and only if <math>p|k.</math>

{{math proof | proof = If <math>p|k</math>, then <math>x^k-1 = (x^p-1)(1 + x^p + x^{2p} + \dots + x^{k/p})</math>. Thus, <math>p(x)|x^k-1.</math>
{{pb}}
Now suppose <math>p(x) | x^k-1</math>. Then, <math>k \geqslant p</math>. We show that <math>k</math> is divisible by <math>p</math> by induction on  <math>k</math>. The base case <math>k = p</math> follows. Therefore, assume  <math>k > p</math>. We know that  <math>p(x)</math> divides both (since it has period <math>p</math>) <math display="block">x^p -1 = (x-1) \left (1 + x + \dots + x^{p-1} \right ) \quad  \text{and}  \quad x^k - 1 = (x-1) \left (1 + x + \dots + x^{k-1} \right ).</math> But  <math>p(x)</math> is irreducible, therefore it must divide both <math>(1 + x + \dots + x^{p-1})</math> and <math>(1 + x + \dots + x^{k-1})</math>; thus, it also divides the difference of the last two polynomials, <math>x^p(1 + x + \dots + x^{p-k-1})</math>. Then, it follows that  <math>p(x)</math> divides <math>(1 + x + \cdots + x^{p-k-1})</math>. Finally, it also divides: <math>x^{k-p}-1 = (x-1)(1 + x + \dots + x^{p-k-1})</math>. By the induction hypothesis, <math>p | k-p</math>, then  <math>p | k</math>.

A corollary to Lemma 2 is that since <math>p(x) = x^p - 1</math> has period <math>p</math>, then <math>p(x)</math> divides <math>x^k-1</math> if and only if <math>p|k</math>.}}}}

{{math theorem | math_statement = The Fire Code is <math>\ell</math>-burst error correcting<ref name = "Moon, Todd" /><ref name = "Lin, Shu" />}}

If we can show that all bursts of length <math>\ell</math> or less occur in different [[coset]]s, we can use them as [[coset leader]]s that form correctable error patterns. The reason is simple: we know that each coset has a unique [[syndrome decoding]] associated with it, and if all bursts of different lengths occur in different cosets, then all have unique syndromes, facilitating error correction.

===Proof of Theorem===
Let  <math>x^i a(x)</math> and <math>x^j b(x)</math> be polynomials with degrees  <math>\ell_1-1</math> and <math>\ell_2-1</math>, representing bursts of length <math>\ell_1</math> and  <math>\ell_2</math> respectively with <math>\ell_1, \ell_2 \leqslant \ell.</math> The integers <math>i, j</math> represent the starting positions of the bursts, and are less than the block length of the code. For contradiction sake, assume that  <math>x^i a(x)</math> and <math>x^j b(x)</math> are in the same coset. Then, <math>v(x) = x^i a(x) + x^j b(x)</math> is a valid codeword (since both terms are in the same coset). Without loss of generality, pick <math>i \leqslant j</math>. By the [[division theorem]] we can write: <math>j-i = g(2\ell -1)+r,</math> for integers <math>g</math> and <math>r, 0 \leqslant r < 2\ell -1</math>. We rewrite the polynomial <math>v(x)</math> as follows:
<math display="block">v(x) = x^ia(x) + x^{i + g(2\ell -1) + r} = x^ia(x) + x^{i + g(2\ell -1) + r} + 2x^{i+r}b(x) = x^i \left (a(x) + x^b b(x) \right ) + x^{i+r}b(x) \left (x^{g(2\ell -1)}+1 \right )</math>

Notice that at the second manipulation, we introduced the term  <math>2x^{i+r}b(x)</math>. We are allowed to do so, since Fire Codes operate on  <math>\mathbb{F}_2</math>. By our assumption, <math>v(x)</math> is a valid codeword, and thus, must be a multiple of <math>g(x)</math>. As mentioned earlier, since the factors of <math>g(x)</math> are relatively prime, <math>v(x)</math> has to be divisible by <math>x^{2\ell -1} + 1</math>. Looking closely at the last expression derived for <math>v(x)</math> we notice that <math>x^{g(2\ell -1)} + 1</math> is divisible by <math>x^{2\ell -1} + 1</math> (by the corollary of Lemma 2). Therefore, <math>a(x) + x^bb(x)</math> is either divisible by <math>x^{2\ell -1} + 1</math> or is <math>0</math>. Applying the division theorem again, we see that there exists a polynomial <math>d(x)</math> with degree <math>\delta</math> such that:
<math display="block">a(x) + x^b b(x) = d(x) (x^{2\ell -1}+1)</math>

Then we may write:
<math display="block">\begin{align}
\delta + 2\ell  -1 & = \deg \left (d(x) \left (x^{2\ell -1}+1 \right ) \right ) \\
&= \deg \left (a(x) + x^b b(x) \right ) \\
&= \deg \left (x^b b(x) \right ) && \deg(a(x)) = \ell_1 - 1 < 2\ell -1 \\
&= b + \ell_2 -1
\end{align}</math>

Equating the degree of both sides, gives us  <math>b = 2\ell  - \ell_2 + \delta.</math> Since  <math>\ell_1, \ell_2 \leqslant \ell</math> we can conclude  <math>b \geqslant \ell + \delta,</math> which implies  <math>b > \ell -1</math> and  <math>b > \delta</math>. Notice that in the expansion:
<math display="block">a(x) + x^bb(x) = 1 + a_1x + a_2x^2 + \dots + x^{\ell_1-1} + x^b \left (1 + b_1x + b_2x^2 + \dots + x^{\ell_2-1} \right ). </math>
The term  <math>x^b</math> appears, but since  <math>\delta < b < 2\ell -1</math>, the resulting expression <math>d(x)(x^{2\ell -1} + 1)</math> does not contain  <math>x^b</math>, therefore <math>d(x) = 0</math> and subsequently <math>a(x) + x^b b(x) = 0.</math> This requires that  <math>b = 0</math>, and  <math>a(x) = b(x)</math>. We can further revise our division of <math>j-i</math> by <math>g(2\ell -1)</math> to reflect <math>b = 0,</math> that is {{nowrap|<math>j-i = g(2\ell -1)</math>.}} Substituting back into <math>v(x)</math> gives us,
<math display="block">v(x) = x^i b(x) \left (x^{j-1}+1 \right ).</math>

Since <math>\deg(b(x)) = \ell_2-1 < \ell</math>, we have <math>\deg(b(x)) < \deg(p(x)) = m</math>. But <math>p(x)</math> is irreducible, therefore  <math>b(x)</math> and <math>p(x)</math> must be relatively prime. Since <math>v(x)</math> is a codeword, <math>x^{j-1}+1</math> must be divisible by <math>p(x)</math>, as it cannot be divisible by <math>x^{2\ell -1} + 1</math>. Therefore, <math>j-i</math> must be a multiple of <math>p</math>. But it must also be a multiple of <math>2\ell -1</math>, which implies it must be a multiple of <math>n = \text{lcm}(2\ell -1,p)</math> but that is precisely the block-length of the code. Therefore, <math>j-i</math> cannot be a multiple of <math>n</math> since they are both less than <math>n</math>. Thus, our assumption of <math>v(x)</math> being a codeword is incorrect, and therefore <math>x^i a(x)</math> and <math>x^j b(x)</math> are in different cosets, with unique syndromes, and therefore correctable.

===Example: 5-burst error correcting fire code===
With the theory presented in the above section, consider the construction of a  <math>5</math>-burst error correcting Fire Code. Remember that to construct a Fire Code, we need an irreducible polynomial <math>p(x)</math>, an integer <math>\ell</math>, representing the burst error correction capability of our code, and we need to satisfy the property that  
<math>2\ell -1</math> is not divisible by the period of <math>p(x)</math>. With these requirements in mind, consider the irreducible polynomial <math>p(x) = 1 + x^2 + x^5</math>, and let <math>\ell = 5</math>. Since <math>p(x)</math> is a primitive polynomial, its period is <math>2^5 - 1 = 31</math>. We confirm that <math>2\ell - 1 = 9</math> is not divisible by <math>31</math>. Thus,
<math display="block">g(x) = (x^9+1) \left (1 + x^2 + x^5 \right ) = 1 + x^2 + x^5 + x^9 + x^{11} + x^{14}</math>
is a Fire Code generator. We can calculate the block-length of the code by evaluating the least common multiple of <math>p</math> and <math>2\ell -1</math>. In other words, <math>n = \text{lcm}(9,31) = 279</math>. Thus, the Fire Code above is a cyclic code capable of correcting any burst of length <math>5</math> or less.

==Binary Reed–Solomon codes==
Certain families of codes, such as [[Reed–Solomon error correction|Reed–Solomon]], operate on alphabet sizes larger than binary. This property awards such codes powerful burst error correction capabilities. Consider a code operating on <math>\mathbb{F}_{2^m}</math>. Each symbol of the alphabet can be represented by <math>m</math> bits. If <math>C</math> is an <math>(n,k)</math> Reed–Solomon code over <math>\mathbb{F}_{2^m}</math>, we can think of <math>C</math> as an <math>[mn,mk]_2</math> code over <math>\mathbb{F}_{2}</math>.

The reason such codes are powerful for burst error correction is that each symbol is represented by <math>m</math> bits, and in general, it is irrelevant how many of those <math>m</math> bits are erroneous; whether a single bit, or all of the <math>m</math> bits contain errors, from a decoding perspective it is still a single symbol error. In other words, since burst errors tend to occur in clusters, there is a strong possibility of several binary errors contributing to a single symbol error.

Notice that a burst of <math>(m+1)</math> errors can affect at most <math>2</math> symbols, and a burst of <math>2m + 1</math> can affect at most <math>3</math> symbols. Then, a burst of <math>tm+1</math> can affect at most <math>t + 1</math> symbols; this implies that a <math>t</math>-symbols-error correcting code can correct a burst of length at most <math>(t-1)m+1</math>.

In general, a <math>t</math>-error correcting Reed–Solomon code over <math>\mathbb{F}_{2^m}</math> can correct any combination of
<math display="block">\frac{t}{1+\lfloor (l+m-2)/m \rfloor}</math>
or fewer bursts of length <math>l</math>, on top of being able to correct <math>t</math>-random worst case errors.

===An example of a binary RS code===
Let <math>G</math> be a <math>[255,223,33]</math> RS code over <math>\mathbb{F}_{2^8}</math>. This code was employed by [[NASA]] in their  [[Cassini-Huygens]] spacecraft.<ref>{{cite web |url=http://quest.arc.nasa.gov/saturn/qa/cassini/Error_correction.txt |title=Cassini: What kind of error correction? |date=1999 |website=quest.arc.nasa.gov |archive-url=https://web.archive.org/web/20120627022807/http://quest.arc.nasa.gov/saturn/qa/cassini/Error_correction.txt |archive-date=2012-06-27}}</ref> It is capable of correcting <math>\lfloor 33/2 \rfloor = 16</math> symbol errors. We now construct a Binary RS Code <math>G'</math> from <math>G</math>. Each symbol will be written using <math>\lceil \log_2(255) \rceil = 8</math> bits. Therefore, the Binary RS code will have <math>[2040,1784,33]_2</math> as its parameters. It is capable of correcting any single burst of length <math>l = 121</math>.

==Interleaved codes==
Interleaving is used to convert convolutional codes from random error correctors to burst error correctors. The basic idea behind the use of interleaved codes is to jumble symbols at the transmitter. This leads to randomization of bursts of received errors which are closely located and we can then apply the analysis for random channel. Thus, the main function performed by the interleaver at transmitter is to alter the input symbol sequence. At the receiver, the deinterleaver will alter the received sequence to get back the original unaltered sequence at the transmitter.

===Burst error correcting capacity of interleaver===
[[File:Row_and_column_major_order.svg|thumb|upright|Illustration of row- and column-major order]]
{{math theorem | math_statement = If the burst error correcting ability of some code is <math>\ell,</math> then the burst error correcting ability of its <math>\lambda</math>-way interleave is <math>\lambda \ell.</math>

{{math proof | proof = Suppose that we have an <math>(n, k)</math> code that can correct all bursts of length <math>\leqslant \ell.</math> [[Forward error correction#Interleaving|Interleaving]] can provide us with a <math>(\lambda n, \lambda k)</math> code that can correct all bursts of length <math>\leqslant \lambda \ell,</math> for any given <math>\lambda</math>. If we want to encode a message of an arbitrary length using interleaving, first we divide it into blocks of length <math>\lambda k</math>. We write the <math>\lambda k</math> entries of each block into a <math>\lambda \times k</math> matrix using row-major order. Then, we encode each row using the <math>(n, k)</math> code. What we will get is a <math>\lambda \times n</math> matrix. Now, this matrix is read out and transmitted in column-major order. The trick is that if there occurs a burst of length <math>h</math> in the transmitted word, then each row will contain approximately <math>\tfrac{h}{\lambda}</math> consecutive errors (More specifically, each row will contain a burst of length at least <math>\lfloor\tfrac{h}{\lambda}\rfloor</math> and at most <math>\lceil\tfrac{h}{\lambda}\rceil</math>). If <math>h \leqslant \lambda \ell,</math> then <math>\tfrac{h}{\lambda} \leqslant \ell</math> and the <math>(n, k)</math> code can correct each row. Therefore, the interleaved <math>(\lambda n, \lambda k)</math> code can correct the burst of length <math>h</math>. Conversely, if <math>h > \lambda \ell,</math> then at least one row will contain more than <math>\tfrac{h}{\lambda}</math> consecutive errors, and the <math>(n, k)</math> code might fail to correct them. Therefore, the error correcting ability of the interleaved <math>(\lambda n, \lambda k)</math> code is exactly <math>\lambda \ell.</math> The BEC efficiency of the interleaved code remains the same as the original <math>(n, k)</math> code. This is true because: 
<math display="block">\frac{2\lambda \ell}{\lambda n - \lambda k}= \frac{2\ell}{n-k}</math>}}
}}
===Block interleaver===
The figure below shows a 4 by 3 interleaver.
[[File:Block,Interleaver,Example.jpg|framed|center|An example of a block interleaver]]

The above interleaver is called as a '''block interleaver'''. Here, the input symbols are written sequentially in the rows and the output symbols are obtained by reading the columns sequentially. Thus, this is in the form of <math>M \times N</math> array. Generally, <math>N</math> is length of the codeword.

'''Capacity of block interleaver''': For an <math>M \times N</math> block interleaver and burst of length <math>\ell,</math> the upper limit on number of errors is <math>\tfrac{\ell}{M}.</math> This is obvious from the fact that we are reading the output column wise and the number of rows is <math>M</math>. By the theorem above for error correction capacity up to <math>t,</math> the maximum burst length allowed is <math>Mt.</math> For burst length of <math>Mt+1</math>, the decoder may fail.

'''Efficiency of block interleaver (<math>\gamma</math>):''' It is found by taking ratio of burst length where decoder may fail to the interleaver memory. Thus, we can formulate <math>\gamma</math> as
<math display="block">\gamma=\frac{Mt+1}{MN} \approx \frac{t}{N}.</math>

'''Drawbacks of block interleaver :''' As it is clear from the figure, the columns are read sequentially, the receiver can interpret single row only after it receives complete message and not before that. Also, the receiver requires a considerable amount of memory in order to store the received symbols and has to store the complete message. Thus, these factors give rise to two drawbacks, one is the latency and other is the storage (fairly large amount of memory). These drawbacks can be avoided by using the convolutional interleaver described below.

===Convolutional interleaver===
Cross interleaver is a kind of multiplexer-demultiplexer system. In this system, delay lines are used to progressively increase length. Delay line is basically an electronic circuit used to delay the signal by certain time duration. Let <math>n</math> be the number of delay lines and <math>d</math> be the number of symbols introduced by each delay line. Thus, the separation between consecutive inputs = <math>nd</math> symbols. Let the length of codeword <math>\leqslant n.</math> Thus, each symbol in the input codeword will be on distinct delay line. Let a burst error of length <math>\ell</math> occur. Since the separation between consecutive symbols is <math>nd,</math> the number of errors that the deinterleaved output may contain is <math>\tfrac{\ell}{nd+1}.</math> By the theorem above, for error correction capacity up to <math>t</math>, maximum burst length allowed is <math>(nd+1)(t-1).</math> For burst length of <math>(nd+1) (t-1) + 1,</math> decoder may fail.
[[File:Convolutional,Interleaver,Example.jpg|framed|center|An example of a convolutional interleaver]]
[[File:Deinterleaver,Example.jpg|framed|center|An example of a deinterleaver]]

'''Efficiency of cross interleaver (<math>\gamma</math>):''' It is found by taking the ratio of burst length where decoder may fail to the interleaver memory. In this case, the memory of interleaver can be calculated as
<math display="block">(0 + 1 + 2 + 3 + \cdots + (n-1))d = \frac{n(n-1)}{2} d.</math>

Thus, we can formulate <math>\gamma</math> as follows:
<math display="block">\gamma = \frac{(nd+1)(t-1)+1}{\frac{n(n-1)}{2}d}.</math>

'''Performance of cross interleaver :''' As shown in the above interleaver figure, the output is nothing but the diagonal symbols generated at the end of each delay line. In this case, when the input multiplexer switch completes around half switching, we can read first row at the receiver. Thus, we need to store maximum of around half message at receiver in order to read first row. This drastically brings down the storage requirement by half. Since just half message is now required to read first row, the latency is also reduced by half which is good improvement over the block interleaver. Thus, the total interleaver memory is split between transmitter and receiver.

==Applications==

=== Compact disc ===

Without error correcting codes, digital audio would not be technically feasible.<ref name = "Algebraic Error Control Codes">Algebraic Error Control Codes (Autumn 2012) – Handouts from Stanford University</ref> The [[Reed–Solomon error correction|Reed–Solomon codes]] can correct a corrupted symbol with a single bit error just as easily as it can correct a symbol with all bits wrong. This makes the RS codes particularly suitable for correcting burst errors.<ref name = "Lin, Shu" /> By far, the most common application of RS codes is in compact discs. In addition to basic error correction provided by RS codes, protection against burst errors due to scratches on the disc is provided by a cross interleaver.<ref name = "Ling, San" />

Current compact disc digital audio system was developed by N. V. Philips of The Netherlands and Sony Corporation of Japan (agreement signed in 1979).

A compact disc comprises a 120&nbsp;mm aluminized disc coated with a clear plastic coating, with spiral track, approximately 5&nbsp;km in length, which is optically scanned by a laser of wavelength ~0.8 μm, at a constant speed of ~1.25&nbsp;m/s. For achieving this constant speed, rotation of the disc is varied from ~8 rev/s while scanning at the inner portion of the track to ~3.5 rev/s at the outer portion. Pits and lands are the depressions (0.12 μm deep) and flat segments constituting the binary data along the track (0.6 μm width).<ref>{{cite book |last=McEliece |first=Robert J. |title=The Theory of Information and Coding: A Mathematical Framework for Communication |publisher=Addison-Wesley |series=Advanced Book Program |date=1977}}</ref>

The CD process can be abstracted as a sequence of the following sub-processes:
* Channel encoding of source of signals
* Mechanical sub-processes of preparing a master disc, producing user discs and sensing the signals embedded on user discs while playing – the channel
* Decoding the signals sensed from user discs

The process is subject to both burst errors and random errors.<ref name = "Algebraic Error Control Codes" /> Burst errors include those due to disc material (defects of aluminum reflecting film, poor reflective index of transparent disc material), disc production (faults during disc forming and disc cutting etc.), disc handling (scratches – generally thin, radial and orthogonal to direction of recording) and variations in play-back mechanism. Random errors include those due to jitter of reconstructed signal wave and interference in signal. CIRC ([[Cross-interleaved Reed–Solomon coding|Cross-Interleaved Reed–Solomon code]]) is the basis for error detection and correction in the CD process. It corrects error bursts up to 3,500 bits in sequence (2.4&nbsp;mm in length as seen on CD surface) and compensates for error bursts up to 12,000 bits (8.5&nbsp;mm) that may be caused by minor scratches.

'''Encoding:''' Sound-waves are sampled and converted to digital form by an A/D converter. The sound wave is sampled for amplitude (at 44.1&nbsp;kHz or 44,100 pairs, one each for the left and right channels of the stereo sound). The amplitude at an instance is assigned a binary string of length 16. Thus, each sample produces two binary vectors from <math>\mathbb{F}_2^{16}</math> or 4 <math>\mathbb{F}_2^{8}</math> bytes of data. Every second of sound recorded results in 44,100&nbsp;&times;&nbsp;32 = 1,411,200 bits (176,400 bytes) of data.<ref name = "Lin, Shu" /> The 1.41 Mbit/s sampled data stream passes through the error correction system eventually getting converted to a stream of 1.88&nbsp;Mbit/s.

Input for the encoder consists of input frames each of 24 8-bit symbols (12 16-bit samples from the A/D converter, 6 each from left and right data (sound) sources). A frame can be represented by  <math>L_1 R_1 L_2 R_2 \ldots L_6 R_6</math> where <math> L_i </math> and <math>R_i</math> are bytes from the left and right channels from the <math>i^{th}</math> sample of the frame.

Initially, the bytes are permuted to form new frames represented by <math>L_1 L_3 L_5 R_1 R_3 R_5 L_2 L_4 L_6 R_2 R_4 R_6</math> where <math>L_i,R_i </math> represent <math>i</math>-th left and right samples from the frame after 2 intervening frames.

Next, these 24 message symbols are encoded using C2 (28,24,5) Reed–Solomon code which is a shortened RS code over <math>\mathbb{F}_{256}</math>. This is two-error-correcting, being of minimum distance 5. This adds 4 bytes of redundancy, <math>P_1 P_2</math> forming a new frame: <math>L_1 L_3 L_5 R_1 R_3 R_5 P_1 P_2 L_2 L_4 L_6 R_2 R_4 R_6</math>. The resulting 28-symbol codeword is passed through a (28.4) cross interleaver leading to 28 interleaved symbols. These are then passed through C1 (32,28,5) RS code, resulting in codewords of 32 coded output symbols. Further regrouping of odd numbered symbols of a codeword with even numbered symbols of the next codeword is done to break up any short bursts that may still be present after the above 4-frame delay interleaving. Thus, for every 24 input symbols there will be 32 output symbols giving <math> R = 24/32 </math>. Finally one byte of control and display information is added.<ref name = "Lin, Shu" /> Each of the 33 bytes is then converted to 17 bits through EFM (eight to fourteen modulation) and addition of 3 merge bits. Therefore, the frame of six samples results in 33 bytes &times;&nbsp;17 bits (561 bits) to which are added 24 synchronization bits and 3 merging bits yielding a total of 588 bits.

'''Decoding:''' The CD player (CIRC decoder) receives the 32 output symbol data stream. This stream passes through the decoder D1 first. It is up to individual designers of CD systems to decide on decoding methods and optimize their product performance. Being of minimum distance 5 The D1, D2 decoders can each correct a combination of <math>e</math> errors and <math>f</math> erasures such that <math>2e+f<5</math>.<ref name = "Lin, Shu" /> In most decoding solutions, D1 is designed to correct single error. And in case of more than 1 error, this decoder outputs 28 erasures. The deinterleaver at the succeeding stage distributes these erasures across 28 D2 codewords. Again in most solutions, D2 is set to deal with erasures only (a simpler and less expensive solution). If more than 4 erasures were to be encountered, 24 erasures are output by D2. Thereafter, an error concealment system attempts to interpolate (from neighboring symbols) in case of uncorrectable symbols, failing which sounds corresponding to such erroneous symbols get muted.

'''Performance of CIRC:<ref name = "Algebraic Error Control Codes" />''' CIRC conceals long bust errors by simple linear interpolation. 2.5&nbsp;mm of track length (4000 bits) is the maximum completely correctable burst length. 7.7&nbsp;mm track length (12,300 bits) is the maximum burst length that can be interpolated. Sample interpolation rate is one every 10 hours at Bit Error Rate (BER) <math>= 10^{-4}</math> and 1000 samples per minute at BER = <math>10^{-3}</math> Undetectable error samples (clicks): less than one every 750 hours at BER = <math>10^{-3}</math> and negligible at BER = <math>10^{-4}</math>.

==See also==
* [[Error detection and correction]]
* [[Error-correcting codes with feedback]]
* [[Code rate]]
* [[Reed–Solomon error correction]]

==References==
{{Reflist|33em}}

[[Category:Coding theory]]
[[Category:Error detection and correction]]
[[Category:Computer errors]]