CN105099468A - High-speed QC-LDPC (quasi-cyclic-low-density parity-check) encoder based on four-level flow lines and used for deep space communication - Google Patents

Abstract

The invention provides a high-speed QC-LDPC (quasi-cyclic-low-density parity-check) encoder based on four-level flow lines and used for deep space communication. The encoder comprises a sparse matrix and vector multiplier, an I-type backward iteration circuit, a high-density matrix and vector multiplier and a II-type backward iteration circuit, wherein the sparse matrix and vector multiplier performs multiplication of a sparse matrix and a vector; the high-density matrix and vector multiplier performs multiplication of a high-density matrix and a vector; the I-type and II-type backward iteration circuits perform backward iterative operation. The whole encoding process is divided into four-level flow lines. The high-speed QC-LDPC encoder adopting a 1/2 code rate and used for a deep space communication system has the advantages of simple structure, low cost, high throughput capacity and the like.

Description

Based on the high speed QC-LDPC encoder of four level production lines in deep space communication

Technical field

The present invention relates to field of channel coding, particularly in a kind of CCSDS deep space communication system based on the high speed QC-LDPC encoder of four level production lines.

Background technology

Low-density checksum (Low-DensityParity-Check, LDPC) code is one of efficient channel coding technology, and quasi-cyclic LDPC (Quasi-CyclicLDPC, QC-LDPC) code is a kind of special LDPC code.The generator matrix G of QC-LDPC code and check matrix H are all the arrays be made up of circular matrix, have the feature of stages cycle, therefore are called as QC-LDPC code.The first trip of circular matrix is the result of footline ring shift right 1, and all the other each provisional capitals are results of its lastrow ring shift right 1, and therefore, circular matrix is characterized by its first trip completely.Usually, the first trip of circular matrix is called as its generator polynomial.

Deep space communication standard adopts the QC-LDPC code of system form, and the left-half of its generator matrix G is a unit matrix, and right half part is by e × c b × b rank circular matrix G _i,jthe array that (0≤i<e, e≤j<t, t=e+c) is formed, as follows:

Wherein, I is b × b rank unit matrixs, and 0 is the full null matrix in b × b rank.The continuous b of G capable and b row are called as the capable and block row of block respectively.From formula (1), G has e block capable and t block row.Deep space communication standard have employed a kind of QC-LDPC code of code check η=1/2, for this code, and t=20, e=8, c=12, b=2048.

In deep space communication standard, the existing solution of 1/2 code check QC-LDPC encoder is the serial encoder adding accumulator (Type-IShift-Register-Adder-Accumulator, SRAA-I) circuit based on 12 I type shift registers.The serial encoder be made up of 12 SRAA-I circuit, completes coding within 16384 clock cycle.The program needs 49152 registers, 24576 two inputs input XOR gate with door and 24576 two, also needs 196608 bit ROM to store the generator polynomial of circular matrix.The program has two shortcomings: one is need a large amount of memory, causes circuit cost high; Two is serial input information bits, and coding rate is slow.

Summary of the invention

In deep space communication system there is the shortcoming that cost is high, coding rate is slow in the existing implementation of 1/2 code check QC-LDPC encoder, for these technical problems, the invention provides a kind of high speed QC-LDPC encoder based on four level production lines.

As shown in Figure 2, the high speed QC-LDPC encoder based on four level production lines in deep space communication system forms primarily of 4 parts: sparse matrix and the multiplier of vector, after I type after iterative circuit, high-density matrix and vectorial multiplier and II type to iterative circuit.Cataloged procedure divides 4 steps to complete: the 1st step, uses sparse matrix and vectorial multiplier compute vector f and w; 2nd step, to iterative circuit compute vector q and x after use I type; 3rd step, uses high-density matrix to verify vectorial p with the multiplier calculating section of vector _x; 4th step, obtains part to iterative circuit compute vector y, y and vectorial q XOR after using II type and verifies vectorial p _y, thus obtain verifying vectorial p=(p _x, p _y).

In deep space communication system provided by the invention, 1/2 code check high speed QC-LDPC coder structure is simple, under the condition significantly improving coding rate, can reduce memory, thus reduces costs, improve throughput.

Be further understood by detailed description and accompanying drawings below about advantage of the present invention and method.

Accompanying drawing explanation

Fig. 1 is the structural representation of near lower triangular check matrix after ranks exchange;

Fig. 2 is the QC-LDPC cataloged procedure based on four level production lines;

Fig. 3 is the functional block diagram of ring shift left accumulator RLA circuit;

Fig. 4 is the multiplier of a kind of high-density matrix and the vector be made up of 4 RLA circuit;

Fig. 5 is sparse matrix and vectorial multiplier;

Fig. 6 gives the annexation of each multi input XOR gate and register in the multiplier of sparse matrix and vector;

Fig. 7 is to iterative circuit after I type;

Fig. 8 gives block position and the ring shift right figure place thereof at nonzero circle matrix place in matrix Q;

Fig. 9 is to iterative circuit after II type;

Figure 10 gives block position and the ring shift right figure place thereof at nonzero circle matrix place in matrix Y;

Figure 11 summarizes each coding step of encoder and the hardware resource needed for whole cataloged procedure and processing time.

Embodiment

Below in conjunction with accompanying drawing, preferred embodiment of the present invention is elaborated, can be easier to make advantages and features of the invention be readily appreciated by one skilled in the art, thus more explicit defining is made to protection scope of the present invention.

The row of circular matrix is heavy identical with column weight, is denoted as w.If w=0, so this circular matrix is full null matrix.If w=1, so this circular matrix is replaceable, is called permutation matrix, and it is by obtaining the some positions of unit matrix I ring shift right.The check matrix H of QC-LDPC code is by c × t b × b rank circular matrix H _j,kthe following array that (1≤j≤c, 1≤k≤t, t=e+c) is formed:

Under normal circumstances, the arbitrary circular matrix in check matrix H is full null matrix (w=0) or is permutation matrix (w=1).Make circular matrix H _j,kfirst trip g _j,k=(g _{j, k, 1}, g _{j, k, 2}..., g _{j, k, b}) be its generator polynomial, wherein g _{j, k, m}=0 or 1 (1≤m≤b).Because H is sparse, so g _j,konly have 1 ' 1 ', even there is no ' 1 '.

For the QC-LDPC code of 1/2 code check in deep space communication system, that front 8 pieces of row of H are corresponding is information vector a, and that rear 12 pieces of row are corresponding is the vectorial p of verification.Be one section with b bit, information vector a is divided into 8 sections, i.e. a=(a ₁, a ₂..., a ₈); Verify vectorial p and be divided into 12 sections, be i.e. p=(p ₁, p ₂..., p ₁₂).

Row swap operation is carried out to check matrix H, is converted near lower triangular shape H _aLT, as shown in Figure 1.The process that row exchange is as follows: front 8 pieces of row keep motionless, and 9th ~ 16 pieces of row arrange exchange with latter 4 pieces.

In FIG, the unit of all matrixes is all b=2048 bit instead of 1 bit.A is made up of 8 × 8 b × b rank circular matrixes, B is made up of 8 × 4 b × b rank circular matrixes, T is made up of 8 × 8 b × b rank circular matrixes, C is made up of 4 × 8 b × b rank circular matrixes, D is made up of 4 × 4 b × b rank circular matrixes, and E is made up of 4 × 8 b × b rank circular matrixes.T is lower triangular matrix, and u=4 reflects check matrix H _aLTwith the degree of closeness of lower triangular matrix.In FIG, matrix A and C corresponding informance vector a, matrix B and the vectorial p of the corresponding part verification of D _x=(p ₁, p ₂..., p ₄), matrix T and E be corresponding remaining verification vector p then _y=(p ₅, p ₆..., p ₁₂).p＝(p _x,p _y)。Above-mentioned matrix and vector meet following relation:

p _x ^Τ＝Φ(ET ^-1Aa ^Τ+Ca ^Τ)(3)

p _y ^Τ＝T ^-1(Aa ^Τ+Bp _x ^Τ)(4)

Wherein, Φ=(ET ^-1b+D) ^-1, subscript ^Τwith ^-1represent transposition and inverse respectively.As everyone knows, circular matrix inverse, product and remain circular matrix.Therefore, Φ is also the array be made up of circular matrix.Although matrix E, T, B and D are sparse matrixes, Φ is no longer sparse but highdensity under normal circumstances.

Make f ^t=Aa ^t, q ^t=T ^{– 1}f ^t, w ^t=Ca ^t, x ^t=Eq ^t+ w ^t, p _x ^t=Φ x ^t, y ^t=T ^{– 1}bp _x ^tand p _y ^t=q ^t+ y ^t.Vector f and w can be calculated by following formula:

${\left[\begin{array}{cc}f& w\end{array}\right]}^T=\left[\begin{array}{c}A\\ C\end{array}\right]a^T={\mathrm{Fa}}^T---\left(5\right)$

Wherein,

$F=\left[\begin{array}{c}A\\ C\end{array}\right]---\left(6\right)$

Q ^t=T ^{– 1}f ^tand x ^t=Eq ^t+ w ^tcan matrix equality be constructed as follows:

$\left[\begin{array}{cc}T& 0\\ E& I\end{array}\right]{\left[\begin{array}{cc}q& x\end{array}\right]}^T=Q{\left[\begin{array}{cc}q& x\end{array}\right]}^T={\left[\begin{array}{cc}f& w\end{array}\right]}^T---\left(7\right)$

Wherein,

$Q=\left[\begin{array}{cc}T& 0\\ E& I\end{array}\right]---\left(8\right)$

Once calculate p _x, y ^t=T ^{– 1}bp _x ^tcan be rewritten as:

[BT][p _xy] ^Τ＝Y[p _xy] ^Τ＝0(9)

Wherein,

Y＝[BT](10)

Because Q is the same with Y and T is all lower triangular matrix, so [qx] in formula (7) and the y in formula (9) can adopt the account form of backward iteration.

Φ relates to high-density matrix and vectorial multiplication, and F relates to sparse matrix and vectorial multiplication, and Q and Y relates to backward iterative computation.Based on the above discussion, a kind of QC-LDPC cataloged procedure based on four level production lines can be provided, as shown in Figure 2.

P _x ^t=Φ x ^tbe equivalent to p _x=x Φ ^t.Make x=(x ₁, x ₂..., x _{u × b}).Definition u bit vectors s _n=(x _n, x _n+b..., x _{n+ (u-1) × b}), wherein 1≤n≤b.Make Φ _j(1≤j≤u) is by Φ ^tjth block row in u × b rank matrix of forming of all circular matrix generator polynomials.Then have

p _j＝(…((0+s ₁Φ _j) ^ls(1)+s ₂Φ _j) ^ls(1)+…+s _bΦ _j) ^ls(1)(11)

Wherein, subscript ^{ls (1)}represent ring shift left 1.

A kind of ring shift left accumulator (Rotate-Left-Accumulator, RLA) circuit can be obtained, as shown in Figure 3 by formula (11).The index of look-up table is u bit vectors s _n, look-up table L _jthe u bit vectors that prior storage is variable and fixing Φ _jinstitute's likely product, therefore need 2 ^uthe read-only memory (Read-OnlyMemory, ROM) of b bit.B bit register R ₁, R ₂..., R _ube respectively used to the array section x cushioning vector x ₁, x ₂..., x _u, b bit register R _u+jfor storing p _xverification section p _j.1 RLA circuit counting vector p _jneed b clock cycle.

For deep space communication system, use 4 RLA circuit to calculate p simultaneously _x=(p ₁, p ₂..., p ₄) be a kind of reasonable plan, high-density matrix as shown in Figure 4 and vectorial multiplier.High-density matrix and vectorial multiplier are by 4 look-up table L ₁, L ₂..., L ₄, 8 b bit register R _3,1, R _3,2..., R _3,8xOR gate X is inputted with 4 b positions two _3,1, X _3,2..., X _3,4composition.Look-up table L ₁, L ₂..., L ₄store 4 variable bit vectors and fixing matrix Φ respectively ₁, Φ ₂..., Φ ₄institute's likely product, register R _3,1, R _3,2..., R _3,4be respectively used to the array section x cushioning vector x ₁, x ₂..., x ₄, register R _3,5, R _3,6..., R _3,8be respectively used to store p _xverification section p ₁, p ₂..., p ₄.4 RLA circuit need use 8192 two to input XOR gate, ROM and 16384 register of 131072 bits.4 RLA circuit counting vector p _xneed 2048 clock cycle.Use high-density matrix and vectorial multiplier compute vector p _xstep as follows:

1st step, resets register R _3,5, R _3,6..., R _3,8, input vector section x ₁, x ₂..., x ₄, by them respectively stored in register R _3,1, R _3,2..., R _3,4in;

2nd step, register R _3,1, R _3,2..., R _3,4ring shift left 1 time simultaneously, XOR gate X _3,1, X _3,2..., X _3,4respectively to look-up table L ₁, L ₂..., L ₄output and register R _3,5, R _3,6..., R _3,8content carry out XOR, XOR result is recycled to move to left after 1 time deposits back register R respectively _3,5, R _3,6..., R _3,8;

3rd step, repeats the 2nd step 2047 times, after completing, and register R _3,5, R _3,6..., R _3,8the content stored is verification section p respectively ₁, p ₂..., p ₄, they constitute part and verify vectorial p _x.

Make f=(f ₁, f ₂..., f ₈) and w=(f ₉, f ₁₀..., f ₁₂), then [fw]=(f ₁, f ₂..., f ₁₂).From formula (5), f _jthe capable and a of the jth block of matrix F ^tproduct, namely

$f_j=H_{j,1}{a}_1^T+H_{j,2}{a}_2^T+...+H_{j,i}{a}_i^T+...+H_{j,9}{a}_9^T---\left(12\right)$

Wherein, 1≤i≤8,1≤j≤12.F _jthe n-th bit f _j,n(1≤n≤b) is

$\begin{array}{c}{f}_{j,n}=g_{j,1}^{rs\left(n-1\right)}{a}_1+g_{j,2}^{rs\left(n-1\right)}{a}_2+\mathrm{...}+g_{j,i}^{rs\left(n-1\right)}{a}_i+\mathrm{...}+g_{j,9}^{rs\left(n-1\right)}{a}_9\\ =g_{j,1}{a}_1^{ls\left(n-1\right)}+g_{j,2}{a}_2^{ls\left(n-1\right)}+\mathrm{...}+g_{j,i}{a}_i^{ls\left(n-1\right)}+\mathrm{...}+g_{j,9}{a}_9^{ls\left(n-1\right)}\end{array}---\left(13\right)$

Wherein, subscript ^{rs (n – 1)}with ^{ls (n – 1)}represent ring shift right n – 1 and ring shift left n – 1 respectively.Since arbitrary circular matrix generator polynomial g _j,ionly have a small amount of ' 1 ' or even complete zero, the inner product so in formula (13) realizes by suing for peace to the tap of ring shift left register, sparse matrix as shown in Figure 5 and vectorial multiplier.Sparse matrix and vectorial multiplier are by 20 2048 bit register R _1,1, R _1,2..., R _1,20with 12 multi input XOR gate X _1,1, X _1,2..., X _1,12composition.Register R _1,1, R _1,2..., R _1,8for loading and ring shift left message segment a ₁, a ₂..., a ₈, register R _1,9, R _1,10..., R _1,20for storing the array section f of [fw] ₁, f ₂..., f ₁₂.Partially connected in Fig. 5 depends on all circular matrix generator polynomials in matrix F.If g _{j, i, m}=1 (1≤m≤b), so message segment a _im bit be connected to XOR gate X _{1, j}.Therefore, register R _{1, i}all taps depend on the nonzero element position of all circular matrix generator polynomials in matrix F i-th piece row, and multi input XOR gate X _{1, j}input depend on matrix F jth block capable in the nonzero element position of all circular matrix generator polynomials.Fig. 6 gives the annexation of each multi input XOR gate and register in the multiplier of sparse matrix and vector.Since all circular matrix generator polynomials in F have α=20 ' 1 ', so sparse matrix needs to use (α – c with the multiplier of vector)=8 two input XOR gate and calculate f simultaneously _{1, n}, f _{2, n}..., f _{12, n}.F and w can calculate complete within 2048 clock cycle.Use sparse matrix as follows with the step of vectorial multiplier compute vector f and w:

1st step, input message segment a ₁, a ₂..., a ₈, by them respectively stored in register R _1,1, R _1,2..., R _1,8in;

2nd step, register R _1,1, R _1,2..., R _1,8ring shift left 1 time simultaneously, XOR gate X _1,1, X _1,2..., X _1,12respectively XOR result is moved to left into register R _1,9, R _1,10..., R _1,20in;

3rd step, repeats the 2nd step 2047 times, after completing, and register R _1,9, R _1,10..., R _1,20the content stored is array section f respectively ₁, f ₂..., f ₁₂, they constitute vector f and w.

Formula (7) implies backward iterative operation, must solve vectorial q and x piecemeal.Definition [qx]=(q ₁, q ₂..., q ₁₂), and be initialized as complete zero.First, q ₁just f is equaled ₁.Secondly, q ₂the 2nd piece of row of matrix Q and vector [qx] ^tlong-pending and f ₂mould 2 He.Then, q ₃the 3rd piece of row of matrix Q and vector [qx] ^tlong-pending and f ₃mould 2 He.Repeat said process, until calculated q ₁₂till, to iterative circuit after I type as shown in Figure 7.After I type to iterative circuit by 12 2048 bit register R _2,1, R _2,2..., R _2,12with 11 multi input modulo 2 adder A _2,2, A _2,3..., A _2,12composition.

To calculate q _j(1≤j≤12) are example.The ring shift right version of the normally unit matrix of the nonzero circle matrix in check matrix H.Have N number of nonzero circle matrix during the jth block of hypothesis matrix Q is capable, their ring shift right figure place is s respectively _{j, k1}, s _{j, k2}..., s _{j, kN}(1≤k1, k2 ..., kN<j).Then,

$\begin{array}{c}{q}_j=f_j+I^{rs\left(s_{j,k1}\right)}{q}_{k1}+I^{rs\left(s_{j,k2}\right)}{q}_{k2}+\mathrm{...}+I^{rs\left(s_{j,kN}\right)}{q}_{kN}\\ =f_j+q_{k1}^{ls\left(s_{j,k1}\right)}+q_{k2}^{ls\left(s_{j,k2}\right)}+\mathrm{...}+q_{kN}^{ls\left(s_{j,kN}\right)}\end{array}---\left(14\right)$

Because N is very little, so formula (14) can calculate complete to the multi input modulo 2 adder of input ring shift left by one within 1 clock cycle.Therefore, compute vector [qx] needs 12 clock cycle altogether.Since a total β=20 nonzero circle matrix in matrix Q, so need to use (β – c) b=16384 two input XOR gate to iterative circuit after I type.

Matrix Q is by 12 × 12 b × b rank circular matrix Q _j,kthe array that (1≤j≤12,1≤k≤12) are formed.Nonzero circle matrix Q _j,ks relative to the ring shift right figure place of b × b rank unit matrix _j,k, 0≤s _j,k<b.For ease of describing, complete zero circular matrix is denoted as s relative to the ring shift right figure place of b × b rank circular matrix _j,k='-'.In the figure 7, nonzero circle matrix Q _j,kcorresponding array section q _kbe recycled the s that moves to left _j,kmulti input modulo 2 adder A is sent into behind position _{2, j}in with array section f _jcarry out XOR, the array section that complete zero circular matrix is corresponding does not participate in XOR, A _{2, j}result of calculation be q _j, stored in register R _{2, j}in.Fig. 8 gives block position and the ring shift right figure place thereof at nonzero circle matrix place in matrix Q.Step to iterative circuit compute vector q and x after use I type is as follows:

1st step, input vector section f ₁, by array section q ₁=f ₁stored in register R _2,1in;

2nd step, input vector section f _j, nonzero circle matrix Q _j,kcorresponding array section q _kbe recycled the s that moves to left _j,kmulti input modulo 2 adder A is sent into behind position _{2, j}in with array section f _jcarry out XOR, XOR result q _jbe stored into register R _{2, j}in, wherein, 2≤j≤12,1≤k<j, 0≤s _j,k<2048;

3rd step, with 1 for step-length increases progressively the value changing j, repeats the 2nd step 10 times, finally, and register R _2,1, R _2,2..., R _2,12that store is array section q respectively ₁, q ₂..., q ₁₂, they constitute vectorial q and x.

Formula (9) also implies backward iterative operation, must solve vectorial y piecemeal.Definition y=(y ₁, y ₂..., y ₈), and be initialized as complete zero.First, y ₁the 1st piece of row of matrix Y and vector [p _xy] ^tlong-pending.Secondly, y ₂the 2nd piece of row of matrix Y and vector [p _xy] ^tlong-pending.Repeat said process, until calculated y ₈till, to iterative circuit after II type as shown in Figure 9.After II type to iterative circuit by 12 2048 bit register R _4,1, R _4,2..., R _4,12with 8 multi input modulo 2 adder A _4,1, A _4,2..., A _4,8composition.Compute vector y needs 8 clock cycle altogether.Since a total ξ=28 nonzero circle matrix in matrix Y, so need to use (ξ – 2c+2u) b=24576 two input XOR gate to iterative circuit after II type.Matrix Y is by 8 × 12 b × b rank circular matrix Y _j,kthe array that (1≤j≤8,1≤k≤12) are formed.Nonzero circle matrix Y _j,ks relative to the ring shift right figure place of b × b rank unit matrix _j,k, 0≤s _j,k<b.Figure 10 gives block position and the ring shift right figure place thereof at nonzero circle matrix place in matrix Y.Step to iterative circuit compute vector y after use II type is as follows:

1st step, input validation section p ₁, p ₂..., p ₄, by them respectively stored in register R _4,9, R _4,10..., R _4,12in;

2nd step, nonzero circle matrix Y _j,kcorresponding array section p _kor y _kbe recycled the s that moves to left _j,kmulti input modulo 2 adder A is sent into behind position _{4, j}in carry out XOR, XOR result y _jbe stored into register R _{4, j}in, wherein, 1≤j≤8,1≤k<4+j, 0≤s _j,k<b;

3rd step, with 1 for step-length increases progressively the value changing j, repeats the 2nd step 7 times, finally, and register R _4,1, R _4,2..., R _4,8that store is array section y respectively ₁, y ₂..., y ₈, they constitute vectorial y.

The invention provides a kind of high speed QC-LDPC coding method based on four level production lines, be applicable to 1/2 code check QC-LDPC code in deep space communication system, its coding step is described below:

1st step, uses sparse matrix and vectorial multiplier compute vector f and w;

2nd step, to iterative circuit compute vector q and x after use I type;

3rd step, uses high-density matrix to verify vectorial p with the multiplier calculating section of vector _x;

4th step, obtains part to iterative circuit compute vector y, y and vectorial q XOR after using II type and verifies vectorial p _y, thus obtain verifying vectorial p=(p _x, p _y).

Figure 11 summarizes each coding step of encoder and the hardware resource consumption needed for whole cataloged procedure and processing time.

Be not difficult to find out from Figure 11, when streamline is full of, whole cataloged procedure needs max (t – c+b, c, u+b)=2056 clock cycle altogether, is less than based on 16384 clock cycle needed for the serial encoding method of 12 SRAA-I circuit.The former coding rate is 8 times of the latter.

In deep space communication standard, the existing solution of 1/2 code check QC-LDPC encoder needs 49152 registers, 24576 two inputs inputs XOR gate with door and 24576 two, also needs the generator polynomial of 196608 bit ROM storage circular matrixes.And the present invention needs 106496 registers, 0 two input inputs XOR gate with door and 49160 two, only need 131072 bit ROM.

As fully visible, compared with traditional serial SRAA method, the present invention has the advantages such as coding rate is fast, memory consumption is few.

The above; be only one of the specific embodiment of the present invention; but protection scope of the present invention is not limited thereto; any those of ordinary skill in the art are in the technical scope disclosed by the present invention; the change can expected without creative work or replacement, all should be encompassed within protection scope of the present invention.Therefore, the protection range that protection scope of the present invention should limit with claims is as the criterion.

Claims (7)

1. in a deep space communication based on the high speed QC-LDPC encoder of four level production lines, the check matrix H of 1/2 code check QC-LDPC code is the array be made up of c × t b × b rank circular matrix, wherein, c=12, t=20, b=2048, e=t-c=8, check matrix H is exchanged by ranks and is transformed near lower triangular shape, can be divided into 6 submatrixs $H=\left[\begin{array}{ccc}A& B& T\\ C& D& E\end{array}\right],$ A is made up of 8 × 8 b × b rank circular matrixes, B is made up of 8 × 4 b × b rank circular matrixes, lower triangular matrix T is made up of 8 × 8 b × b rank circular matrixes, C is made up of 4 × 8 b × b rank circular matrixes, D is made up of 4 × 4 b × b rank circular matrixes, E is made up of 4 × 8 b × b rank circular matrixes, Φ=(ET ^-1b+D) ^-1be made up of 4 × 4 b × b rank circular matrixes, Φ _jby Φ ^tjth block row in 4 × b rank matrix of forming of all circular matrix generator polynomials, wherein, subscript ^Τwith ^-1represent transposition and inverse respectively, 1≤j≤4, $Q=\left[\begin{array}{cc}T& 0\\ E& I\end{array}\right]$ By 12 × 12 b × b rank circular matrix Q _j,kform, wherein, I is unit matrix, and 0 is full null matrix, 1≤j≤12,1≤k≤12, nonzero circle matrix Q _j,ks relative to the ring shift right figure place of b × b rank unit matrix _j,k, wherein, 0≤s _j,k<b, Y=[BT] are by 8 × 12 b × b rank circular matrix Y _j,kform, wherein, 1≤j≤8,1≤k≤12, nonzero circle matrix Y _j,ks relative to the ring shift right figure place of b × b rank unit matrix _j,k, wherein, 0≤s _j,k<b, A and C corresponding informance vector a, matrix B and the vectorial p of the corresponding part verification of D _x, matrix T and E be corresponding remaining verification vector p then _y, verify vectorial p=(p _x, p _y), be one section with b bit, information vector a is divided into 8 sections, i.e. a=(a ₁, a ₂..., a ₈), verify vectorial p and be divided into 12 sections, be i.e. p=(p ₁, p ₂..., p ₁₂), p _x=(p ₁, p ₂..., p ₄), p _y=(p ₅, p ₆..., p ₁₂), vector f is divided into 8 sections, i.e. f=(f ₁, f ₂..., f ₈), vectorial w is divided into 4 sections, i.e. w=(f ₉, f ₁₀..., f ₁₂), [fw]=(f ₁, f ₂..., f ₁₂), vectorial q is divided into 8 sections, i.e. q=(q ₁, q ₂..., q ₈), vector x is divided into 4 sections, i.e. x=(q ₉, q ₁₀..., q ₁₂), [qx]=(q ₁, q ₂..., q ₁₂), vectorial y is divided into 8 sections, i.e. y=(y ₁, y ₂..., y ₈), it is characterized in that, described encoder comprises following parts:

Sparse matrix and vectorial multiplier, by 20 2048 bit register R _1,1, R _1,2..., R _1,20with 12 multi input XOR gate X _1,1, X _1,2..., X _1,12composition, for compute vector f and w;

To iterative circuit after I type, by 12 2048 bit register R _2,1, R _2,2..., R _2,12with 11 multi input modulo 2 adder A _2,2, A _2,3..., A _2,12composition, for compute vector q and x;

High-density matrix and vectorial multiplier, by 4 look-up table L ₁, L ₂..., L ₄, 8 b bit register R _3,1, R _3,2..., R _3,8xOR gate X is inputted with 4 b positions two _3,1, X _3,2..., X _3,4composition, verifies vectorial p for calculating section _x, look-up table L ₁, L ₂..., L ₄store 4 variable bit vectors and fixing matrix Φ respectively ₁, Φ ₂..., Φ ₄institute's likely product;

To iterative circuit after II type, by 12 2048 bit register R _4,1, R _4,2..., R _4,12with 8 multi input modulo 2 adder A _4,1, A _4,2..., A _4,8composition, obtains part for compute vector y, y and vectorial q XOR and verifies vectorial p _y, thus obtain verifying vectorial p=(p _x, p _y).

2. in a kind of deep space communication according to claim 1 based on the high speed QC-LDPC encoder of four level production lines, it is characterized in that, the process that described ranks exchange is as follows: front 8 pieces of row keep motionless, and 9th ~ 16 pieces of row arrange exchange with latter 4 pieces.

3. in a kind of deep space communication according to claim 1 based on the high speed QC-LDPC encoder of four level production lines, it is characterized in that, multiplier compute vector f and the step of w of described sparse matrix and vector are as follows:

1st step, input message segment a ₁, a ₂..., a ₈, by them respectively stored in register R _1,1, R _1,2..., R _1,8in;

2nd step, register R _1,1, R _1,2..., R _1,8ring shift left 1 time simultaneously, XOR gate X _1,1, X _1,2..., X _1,12respectively XOR result is moved to left into register R _1,9, R _1,10..., R _1,20in;

3rd step, repeats the 2nd step 2047 times, after completing, and register R _1,9, R _1,10..., R _1,20the content stored is array section f respectively ₁, f ₂..., f ₁₂, they constitute vector f and w.

4. in a kind of deep space communication according to claim 1 based on the high speed QC-LDPC encoder of four level production lines, it is characterized in that, the step to iterative circuit compute vector q and x after described I type is as follows:

1st step, input vector section f ₁, by array section q ₁=f ₁stored in register R _2,1in;

2nd step, input vector section f _j, nonzero circle matrix Q _j,kcorresponding array section q _kbe recycled the s that moves to left _j,kmulti input modulo 2 adder A is sent into behind position _{2, j}in with array section f _jcarry out XOR, XOR result q _jbe stored into register R _{2, j}in, wherein, 2≤j≤12,1≤k<j, 0≤s _j,k<2048;

3rd step, with 1 for step-length increases progressively the value changing j, repeats the 2nd step 10 times, finally, and register R _2,1, R _2,2..., R _2,12that store is array section q respectively ₁, q ₂..., q ₁₂, they constitute vectorial q and x.

5. in a kind of deep space communication according to claim 1 based on the high speed QC-LDPC encoder of four level production lines, it is characterized in that, described high-density matrix with vector multiplier compute vector p _xstep as follows:

1st step, resets register R _3,5, R _3,6..., R _3,8, input vector section x ₁, x ₂..., x ₄, by them respectively stored in register R _3,1, R _3,2..., R _3,4in;

2nd step, register R _3,1, R _3,2..., R _3,4ring shift left 1 time simultaneously, XOR gate X _3,1, X _3,2..., X _3,4respectively to look-up table L ₁, L ₂..., L ₄output and register R _3,5, R _3,6..., R _3,8content carry out XOR, XOR result is recycled to move to left after 1 time deposits back register R respectively _3,5, R _3,6..., R _3,8;

3rd step, repeats the 2nd step 2047 times, after completing, and register R _3,5, R _3,6..., R _3,8the content stored is verification section p respectively ₁, p ₂..., p ₄, they constitute part and verify vectorial p _x.

6. in a kind of deep space communication according to claim 1 based on the high speed QC-LDPC encoder of four level production lines, it is characterized in that, the step to iterative circuit compute vector y after described II type is as follows:

1st step, input validation section p ₁, p ₂..., p ₄, by them respectively stored in register R _4,9, R _4,10..., R _4,12in;

2nd step, nonzero circle matrix Y _j,kcorresponding array section p _kor y _kbe recycled the s that moves to left _j,kmulti input modulo 2 adder A is sent into behind position _{4, j}in carry out XOR, XOR result y _jbe stored into register R _{4, j}in, wherein, 1≤j≤8,1≤k<4+j, 0≤s _j,k<b;

3rd step, with 1 for step-length increases progressively the value changing j, repeats the 2nd step 7 times, finally, and register R _4,1, R _4,2..., R _4,8that store is array section y respectively ₁, y ₂..., y ₈, they constitute vectorial y.

7. in a deep space communication based on the high speed QC-LDPC coding method of four level production lines, the check matrix H of 1/2 code check QC-LDPC code is the array be made up of c × t b × b rank circular matrix, wherein, c=12, t=20, b=2048, e=t-c=8, check matrix H is exchanged by ranks and is transformed near lower triangular shape, can be divided into 6 submatrixs $H=\left[\begin{array}{ccc}A& B& T\\ C& D& E\end{array}\right],$ A is made up of 8 × 8 b × b rank circular matrixes, B is made up of 8 × 4 b × b rank circular matrixes, lower triangular matrix T is made up of 8 × 8 b × b rank circular matrixes, C is made up of 4 × 8 b × b rank circular matrixes, D is made up of 4 × 4 b × b rank circular matrixes, E is made up of 4 × 8 b × b rank circular matrixes, Φ=(ET ^-1b+D) ^-1be made up of 4 × 4 b × b rank circular matrixes, Φ _jby Φ ^tjth block row in 4 × b rank matrix of forming of all circular matrix generator polynomials, wherein, subscript ^Τwith ^-1represent transposition and inverse respectively, 1≤j≤4, $Q=\left[\begin{array}{cc}T& 0\\ E& I\end{array}\right]$ By 12 × 12 b × b rank circular matrix Q _j,kform, wherein, I is unit matrix, and 0 is full null matrix, 1≤j≤12,1≤k≤12, nonzero circle matrix Q _j,ks relative to the ring shift right figure place of b × b rank unit matrix _j,k, wherein, 0≤s _j,k<b, Y=[BT] are by 8 × 12 b × b rank circular matrix Y _j,kform, wherein, 1≤j≤8,1≤k≤12, nonzero circle matrix Y _j,ks relative to the ring shift right figure place of b × b rank unit matrix _j,k, wherein, 0≤s _j,k<b, A and C corresponding informance vector a, matrix B and the vectorial p of the corresponding part verification of D _x, matrix T and E be corresponding remaining verification vector p then _y, verify vectorial p=(p _x, p _y), be one section with b bit, information vector a is divided into 8 sections, i.e. a=(a ₁, a ₂..., a ₈), verify vectorial p and be divided into 12 sections, be i.e. p=(p ₁, p ₂..., p ₁₂), p _x=(p ₁, p ₂..., p ₄), p _y=(p ₅, p ₆..., p ₁₂), vector f is divided into 8 sections, i.e. f=(f ₁, f ₂..., f ₈), vectorial w is divided into 4 sections, i.e. w=(f ₉, f ₁₀..., f ₁₂), [fw]=(f ₁, f ₂..., f ₁₂), vectorial q is divided into 8 sections, i.e. q=(q ₁, q ₂..., q ₈), vector x is divided into 4 sections, i.e. x=(q ₉, q ₁₀..., q ₁₂), [qx]=(q ₁, q ₂..., q ₁₂), vectorial y is divided into 8 sections, i.e. y=(y ₁, y ₂..., y ₈), it is characterized in that, described coding method comprises the following steps:

1st step, uses sparse matrix and vectorial multiplier compute vector f and w;

2nd step, to iterative circuit compute vector q and x after use I type;

3rd step, uses high-density matrix to verify vectorial p with the multiplier calculating section of vector _x;

4th step, obtains part to iterative circuit compute vector y, y and vectorial q XOR after using II type and verifies vectorial p _y, thus obtain verifying vectorial p=(p _x, p _y).