Method for calculating phase shift coefficients of an M sequence

Abstract: An n-bit binary value corresponding to an amount of phase shift d is assigned to a shift register (302), and a shift operation is performed. An n-bit vector value corresponding to a decimal value "1" is assigned to a latch (305) as an initial value. Thereafter, the input from a switch (303) is sequentially stored. A multiplier (307) performs a square operation within a Galois field GF (2 n ) for the output of the latch. A doubler (309) performs a double operation within the Galois field GF (2 n ) for the output of the multiplier (307). The switch (303) selects either of the outputs of the multiplier (307) or the doubler (309) according to the output value from the MSB side of the shift register (302). After the shift operation and the latch operation are performed a number of times n, the n-bit output of the latch (305) is output as respective phase shift coefficients b 0 through b n-1 .

Description: Background of the Invention Field of the Invention The present invention relates to a method for generating a pseudorandom noise sequence with an arbitrarily designated phase, which is applied to a communication system using a spread-modulation method. Description of the Related Art In a communications system using spread spectrum modulation, the spectra of many signals can be spread over a broadband, multiplexed, and transmitted by using a Code Division Multiple Access (CDMA). Fig. 1 shows one principle of the configuration of a CDMA communication system. On a CDMA transmitter 901 side, a spread-modulation unit 905 spread-spectrum-modulates the transmission signal which is output from a transmission signal source 903, for example, which is frequency-modulated or phase-modulated by using the spread code generated by a spread code generating unit 904, and the resultant transmission signal is transmitted to a transmission line 906. On a CDMA receiver 902 side, a despread-demodulation unit 908 must despread (demodulate) a reception signal by using the despread code which has the same sequence and phase as those of the spread code on the transmitter side and which is output from a despread code generating unit 907 in synchronization with the timing of the transmitter side. Accordingly, the despread code generating unit 907 must have the capability for generating a sequence code having an arbitrary phase according to a timing synchronization signal (normally, this signal is autonomously generated from a reception signal within the CDMA receiver 902). In a CDMA communication, the spread code (and the despread code) for spreading a spectrum must satisfy the following conditions in addition to the condition that the spread code must be a broadband signal: (1) the number of types of spread codes must be large in order to allow codes to be assigned to many users; (2) a cross-correlation must be small in order to allow the spread code to be identified from a different user code; (3) a self-correlation must be strictly identified in order to ensure the synchronization with the signal addressed to a local station; and (4) the spread code must be as random as possible, have a long cycle, and be difficult to be decoded in order to improve the confidentiality of a communication signal. Conventionally, a PN (Pseudorandom Noise) sequence is known as the code for satisfying such conditions. Since the PN sequence can be generated by using a shift register, its generation process is not really random but deterministic. However, the PN sequence is a code having the following properties of randomness. Therefore, this is suitable as the spread code of the CDMA communication, which requires the above described conditions. Property 1: Balance Property The numbers of times that "1" and "0" respectively appear in one cycle of the sequence are different only by 1. Property 2: Run Property In the runs of "1s" and runs of "0s" which are included in one cycle, the length of each run is "1" when a classification number of that run is 1/2, "2" when a classification number of that run is 1/4, "3" when a classification number of that run is 1/8, ···. That is, there is a number of {(classification number of run) × (1/2k)} of the runs which have the run number k. Note that, for the runs of which this number is less than 1, these runs become meaningless runs. Property 3: Correlation Property If sequences are made cyclic and a comparison is made between the respective code values of the corresponding digits of the two sequences in every state, the number of the digits whose code values match and that of the digits whose code values do not match are different only by 1. An M sequence (maximum length sequence) is known as a typical PN sequence satisfying such properties. The M sequence is generated using the circuit including an n-stage shift register, which is shown in Fig. 2. In Fig. 2, the respective outputs of the stages of the n-stage shift register are multiplied with a coefficient f i (0 or 1), and the multiplication results are fed back to the input side of the shift register via exclusive-OR circuits (+ signs encircled in this figure). If the coefficient f 1 satisfies a particular condition when the values of all of the stages of the shift register are not "0" in an initial state, the cycle of the sequence a i output from the shift register will become the maximum cycle (2 n -1) that the n-stage shift register can generate. Such a sequence is referred to as an M sequence. The circuit shown in Fig. 2 can be represented by the following equation. If f n =1 is assigned to this equation, the following equation can be obtained. The above described equations (1) and (2) are referred to as linear recurring equations. Here, if a delay operator x which satisfies a i+1 =x j a i is assigned, the equation (2) becomes as follows. The polynomial f(x) of the following equation, which is represented by the term on the left side of the above equation (3), is referred to as a characteristic polynomial. If the coefficient f j included in this equation (4) belongs to a Galois field GF (2), and if f(x) is the minimal polynomial possessed by a primitive element α of the Galois field GF (2 n ), it is known that the circuit shown in Fig. 2, which includes the n-stage shift register, can generate the M sequence having the maximum cycle (2 n -1). This minimal polynomial is referred to as a primitive polynomial of degree k. Its details are referred to, for example, in the document "Sensing/Recognition Series Vol. 8, M Sequence and its Application, pp. 16-", written by J. Kashiwagi and published by Shokodo. The primitive polynomial can be calculated as described on pp. 171 to 191 of this document, and many types of primitive polynomials were previously obtained in some of the papers cited in this document. For example, the coefficient f j , which is included in the equation (4) and corresponds to a primitive polynomial f(x)=x 4 +x+1 of the Galois field GF (2 4 ), becomes f 0 =1, f 1 =1, f 2 =f 3 =0, and f 4 =1. As a result, the M sequence generating circuit shown in Fig. 3 can be configured based on the circuit shown in Fig. 2. Here, an M sequence x d a i whose phase is shifted by d bits from the output a i of the M sequence will be obtained. If a predetermined initial state of the n-stage (4 stages in Fig. 3) shift register is provided, all of the states of the M sequence, which succeed the initial state, are determined. Therefore, the M sequence having an arbitrary phase is proved to be obtained with the linear combination of the outputs of the respective stages of the shift register, as represented by the following equation (5). x d a i = b 0 x 0 a i + b 1 x 1 a i + b 2 x 2 a i + ··· + b n -1 x n -1 a i Consequently, the circuit for generating the M sequence having an arbitrary phase can be configured from the M sequence generating circuit shown in Fig. 3, which includes the 4-stage shift register, as shown in Fig. 4. In Fig. 4, an initial value is assigned to each of the stages of the 4-stage shift register (SR) 1203 in a PN generator (PNG) 1201. The feedback equivalent to that shown in Fig. 3 is provided by a TAP 1204. TAP information (TAPINFO) 1205, which corresponds to the respective coefficients b 0 through b 3 included in the equation (5), is provided to 4 AND circuits (ANDs) 1206 in a variable tap (ATAP) 1202. Consequently, the output selected according to the TAPINFO 1205 from among the outputs of the respective stages of the SR 1203 is added to another output by the corresponding AND circuit (AND) 1206 and exclusive-OR circuit (EXOR) 1207, whereby the M sequence x d a i having an arbitrary phase d is output as the result of the addition. The means for providing the initial value to the SR 1203, the means for providing the clock for implementing a shift operation, etc. are omitted and are not shown in Fig. 4. Provided next is the explanation about the principle for calculating the coefficients b 0 through b 3 , which are included in the equation (5) and form the TAPINFO 1205. In the PNG 1201 shown in Fig. 4, the M sequences x 1 a i to x 3 a i , whose phases are respectively shifted by 1 to 3 bits from the output a i of the M sequence are the outputs themselves of the shift register stages SR1 through SR4 in the second to fourth stages, as shown in Fig. 5. Next, the M sequence x 4 a i , whose phase is shifted by 4 bits from the output a i of the M sequence will be considered. In this case, if the primitive polynomial f(x)=x 4 +x+1 is set f ( x ) = x 4 + x +1 = 0 the following equations are satisfied with the operation in the Galois field. x 4 = 1+ x x 4 a i = x 0 a i + x 1 a i According to the above described equation (9), the M sequence x 4 a i whose phase is shifted by 4 bits from the output a i of the M sequence, is proved to be represented by the exclusive-OR operation performed between the output a i and the output x 1 a i whose phase is shifted by 1 bit. That is, the M sequence x 4 a i whose phase is shifted by 4 bits from the output a i can be obtained with the exclusive-OR operation performed between the output of the first shift register stage SRO and that of the second shift register stage SR1, as shown in Fig. 6. Namely Then, the M sequence x 5 a i whose phase is shifted by 5 bits from the output a i of the M sequence can be obtained with the exclusive-OR operation performed between the output of the second shift register stage SR1 and that of the third shift register stage SR2 according to following equation (11) obtained by multiplying both sides of the equation (9) by x, as shown in Fig. 6. x 5 a i = x 1 a i + x 2 a i Namely, Furthermore, the M sequence x 6 a i whose phase is shifted by 6 bits from the output a i of the M sequence can be obtained with the exclusive-OR operation performed between the output of the third shift register stage SR2 and that of the fourth shift register stage SR3 according to the following equation (13) obtained by multiplying both sides of the equation (11) by x. x 6 a i = x 2 a i + x 3 a i (13) Namely, The configuration shown in Fig. 7 is conventionally known as the configuration for sequentially outputting the coefficients b 0 through b 3 , which form the TAPINFO 1205, according to the above described rules. In this configuration, a 4-stage shift register (SR) 1501 is used in correspondence with the primitive polynomial f(x)=x 4 +x+1, and an exclusive-OR circuit EXOR 1502 is inserted between the output side of the shift register stage a 0 and the input side of the shift register stage a 1 , which respectively correspond to the terms 1=x 0 and x=x 1 forming the right side of the equation (7). The output of the output stage a 3 of the shift register (corresponding to the term x 4 on the left side of the equation (7)) is fed back to the EXOR 1502. Since the stage preceding the shift register stage a 0 does not exist, the output of the output stage a 3 is directly fed back to the input side of the shift register stage a 0 . More typically, a shift register with the number of stages corresponding to the degree of the primitive polynomial f(x) is used. Similar to the equation (7), an equation is formed by setting f(x) to "0" so that the term of the highest degree is included on the left side and the terms of other degrees are included on the right side. Then, an exclusive-OR circuit is inserted into the input side of the shift register stage corresponding to each of the terms forming the right side of the formed equation, and the output (corresponding to each of the terms on the left side of the equation) of an output stage of the shift register is fed back to the exclusive-OR circuit. Next, "1" is assigned to the first shift register stage a 0 as the initial value of the SR 1501, while "0" is assigned to the remaining shift register stages a 1 through a 3 shown in Fig. 7. The shift operation is then performed the number of times which corresponds to a desired amount of phase shift, so that the respective coefficients b 0 through b 3 forming the TAPINFO 1205 shown in Fig. 4 are determined as the respective outputs of the stages a 0 through a 3 of the SR 1501. With the conventional method shown in Fig. 7, however, the shift operation must be performed the number of times which corresponds to a desired amount of phase shift in order to calculate the TAPINFO 1205 corresponding to the desired amount of phase shift. Therefore, if the cycle of the M sequence becomes longer (for example, 10 minutes or so), it requires a huge amount of time to calculate the TAPINFO -1205. Summary of the Invention The present invention was developed in the above described background, and aims at quick calculation of each piece of tap information corresponding to each amount of phase shift. The present invention assumes the technique for calculating the respective phase shift coefficients b 0 through b n-1 intended for obtaining the M sequence x d a i whose phase is shifted by d bits from the output a i of the M sequence generated with a primitive polynomial f(x) of degree n, by using a linear combination b 0 x 0 a i + b 1 x 1 a i + b 2 x 2 a i + ··· + b n-1 x n-1 a i of the respective M sequences x 0 a i through x n-1 a i whose phases are shifted by 0 to n-1 bits from the output a i . In a first step of the present invention, the binary value of the amount of phase shift d is input. In a second step, the n-bit vector value corresponding to a decimal value "1" is set as an initial input vector value. In a third step, a target bit is set as the most significant bit of the binary value of the amount of phase shift d input in the first step. In a fourth step, the vector value of the primitive element α is multiplied with the input vector value within the Galois field GF (2 n ) and the result of the multiplication is defined as an output vector value if the value of the target bit is "1", and the input vector value is defined as the output vector value unchanged if the value of the target bit is "0". In a fifth step, a square operation is performed for the output vector value obtained in the fourth step within the Galois field GF (2 n ). In a sixth step, the position of the target bit is shifted to the least significant bit side by 1 bit in the binary value of the amount of phase shift d input in the first step, and the result of the square operation performed in the fifth step is used as a newly input vector value, with which the fourth and fifth steps are performed. In a seventh step, the respective elements of the operation results obtained so far are output as the phase shift coefficients b 0 through b n-1 , when the fourth and fifth steps have been performed for all of the bits structuring the binary value of the amount of phase shift d input in the first step. According to the present invention, phase shift coefficients can be calculated simply by repeatedly performing the operation for multiplying a current vector value with the vector representation value of the primitive element α and squaring the result of the vector multiplication, or the operation for squaring the unchanged current vector value, an approximate number of {log(d)/log(2)} times. These two types of operations within the Galois field GF (2 n ) can be implemented by hardware of a compact size. In this way, as the cycle of an M sequence becomes longer, the number of repetition times that the phase shift coefficients of the M sequence are calculated can be significantly reduced compared with that of the conventional method. Brief Description of the Drawings One skilled in the art can easily understand additional objects and features of the present invention from the description of the preferred embodiments and the attached drawings. In the drawings: Fig. 1 shows one principle of the configuration of a CDMA communication system; Fig. 2 shows the configuration of the circuit for generating an M sequence; Fig. 3 shows the configuration of the circuit for generating a 4-stage M sequence; Fig. 4 shows the configuration of the circuit for generating an M sequence having an arbitrary phase; Fig. 5 is a configuration diagram for explaining the principle of the generation of an M sequence having an arbitrary phase (Nol. 1); Fig. 6 is a configuration diagram for explaining the principle of the generation of an M sequence having an arbitrary phase (No. 2); Fig. 7 is a configuration diagram showing the configuration of a conventional circuit for generating tap information TAPINFO in order to generate an M sequence having an arbitrary phase; Fig. 8 shows the table including the vector representation of the primitive element α of the Galois field GF (2 4 ), which corresponds to the primitive polynomial f(x)=x 4 +x+1, and its power; Fig. 9 shows the principle of the operations according to the present invention; Fig. 10 shows the configuration of the circuit according to a preferred embodiment of the present invention; Fig. 11 is a timing chart showing the operations of the circuit according to the preferred embodiment of the present invention; Fig. 12 shows the first configuration of a square operation unit (MUL); Fig. 13 shows the second configuration of the square operation unit (MUL); Fig. 14 shows the first configuration of a double operation unit (DBL); and Fig. 15 shows the second configuration of the double operation unit (DBL). Description of the Preferred Embodiments Hereinafter, the details of a preferred embodiment according to the present invention will be described by referring to the drawings. <Principle of the preferred embodiment according to the present invention> Provided first is the explanation about the principle of the preferred embodiment according to the present invention. The preferred embodiment according to the present invention assumes the configuration shown in Fig. 4, which includes a 4-stage shift register (SR) 1203 and can generate an M sequence having an arbitrary phase, as an example. Assume that the primitive element of the Galois field GF (2 n ) is α. In this case, each of the coefficients b 0 through b n-1 included in the above described equation (5), which is the TAPINFO 1205 shown in Fig. 4, is equal to the value of each bit obtained by calculating α d within the Galois field GF (2 n ) for the n-bit vector which represents the primitive element α. Fig. 8 shows the vector representation of the primitive element a of the Galois field GF (2 4 ), which corresponds to the primitive polynomial f(x)=x 4 +x+1, and the vector representation of the power of the primitive element α. In this case, the vector representation of the primitive element a is (0(MSB), 0, 1, 0(LSB)) (a decimal value "2"). For example, the respective bit values of the vector representations α 1 through α 6 match the respective coefficient values b 0 through b 3 represented by the above described equations (6), (10), (12), and (14), which are used in the conventional method. Here, α d can be decomposed as represented by the following equation (15). By repeating the operation for replacing s with d in the term α s obtained by the decomposition operation of the equation (15) and for further performing the decomposition operation of the equation (15), for example, α 6 can be decomposed as follows. As a result, the vector representation value of α d can be calculated only by repeatedly performing the operation for multiplying the current vector value with the vector representation value of a and squaring the vector multiplication result, or the operation for squaring the unchanged current vector value sequentially from the innermost term of the structure which is parenthesized and decomposed as the right side of the equation (16) within the Galois field GF (2 n ) an approximate number of {log(d)/log(2)} times. If the operation for multiplying the current vector value with the vector representation value of α and squaring the vector multiplication result or the operation for squaring the unchanged current vector value within the Galois field GF (2 n ), can be implemented by hardware, it becomes possible to calculate the vector representation value of the element α d , that is, the above described coefficients b 0 through b n-1 , by the number of repetition times, which becomes much fewer than the conventionally required number of times d that the shift operation is performed, as the number of stages n is larger. Fig. 9 shows the principle of the algorithm for identifying the above described repetition operation for calculating the coefficients b 0 through b n-1 of the equation (5), which is the TAPINFO 1205 shown in Fig. 4, if an arbitrary amount of phase shift d is provided. Step 1: First of all, the binary value of the amount of phase shift d is input. Assuming that the amount of phase shift d equals 6, the corresponding 4 bits binary value will become (0(MSB), 1, 1, 0(LSB)) as shown in Fig. 9(a). Step 2: The vector value of n bits (0(MSB), 0, 0, 1(LSB)) (a decimal value "1") is defined as an initial input vector value. Step 3: A target bit is defined as an MSB (Most Significant Bit) of the binary value of the amount of phase shift d input in Step 1. Step 4: If the target bit is "1", the input vector value is multiplied with the vector value of the primitive element α within the Galois field GF (2 n ), and its result is used as an output vector value in Step 4 (operation I). For example, the vector value of the primitive element a of the Galois field GF (2 4 ), which corresponds to the primitive polynomial f(x)=x 4 +x+1, is (0(MSB), 0, 1, 0(LSB)), that is, a decimal value 2 as shown in Fig. 8. Accordingly, a 1-bit-shift operation is performed for the input vector value in this case. If the target bit is "0", the input vector value is used as the output vector value in Step 2 as it is. Step 5: The output vector value in Step 4 is squared within the Galois field GF (2 n ), and its operation result is used as the output vector value of Step 5 (operation II). Specifically, if an "i"th element a i (0≤i≤n-1) of the output vector value of Step 4 is set, the element a i is cleared and "1" resulting from modulo 2 is added to an element a 2i in modulo 2. The result of the addition is used as the output vector value of Step 5. Note that, however, if 2 i is equal to or larger than n, "1" resulting from modulo 2 is added to one or more elements a j which are equivalent to α 2i determined by the primitive polynomial f(x) within the Galois field GF (2 n ), and the result of the addition is defined as the output vector value of Step 5. Step 6: The position of the target bit is shifted to the LSB side by 1 bit in the binary value of the amount of phase shift d input in Step 1, and the output vector value of Step 5 is used as a newly input vector value. Then. Steps 4 and 5 are again performed. Upon termination of performing Steps 4 and 5 for the LSB, each element of the output vector value obtained in Step 5 is used as the TAPINFO 1205 shown in Fig. 4. The algorithm used in Steps 1 through 6 is implemented by hardware, thereby quickly calculating the TAPINFO 1205 shown in Fig. 4 for an arbitrary amount of phase shift d. Fig. 9(b) shows, the procedure for performing the above described Steps 1 through 6 for 4 bits of phase shi ft (0(MSB), 1, 1, 0(LSB)), which corresponds to the amount of phase shift d=6 and is shown in Fig. 9(a). Fig. 9(c) shows the results of the respective arithmetic operations shown in Fig. 9(b). The final operation result α 6 =α 3 +α 2 in Fig . 9(c) matches the conventional operation result represented by the equation (13). <The Detailed Configuration of the Preferred Embodiment According to the Present Invention> Fig. 10 assumes the circuit shown in Fig. 4, which includes the 4-stage shift register (SR) 1203 and can generate an M sequence having an arbitrary phase. This figure shows the configuration of the circuit according to the preferred embodiment according of the present invention, which is intended for generating the TAPINFO 1205 shown in Fig. 4. Fig. 11 is a timing chart showing the operations of the circuit shown in Fig. 10. If the primitive polynomial is assumed to be f(x)=x 4 +x+1, the signals such as the amount of phase shift (SFTVAL) 301, a latch output (LATOUT) 306, a square operation unit output (MULOUT) 308, and a double operation output (DBLOUT) 310, etc. will become the 4 bits vector signals. The SFTVAL 301 is assigned to a shift amount register (SREG) 302 beforehand. This operation corresponds to Step 1 of the operational algorithm referred to in the above described principle of the preferred embodiment according to the present invention. In the example shown in Fig. 11, for example, the amount of phase shift d equals 6, and its 4-bit binary value (0(MSB), 1, 1, 0(LSB)) is assigned to the SREG 302. Additionally, the 4-bit vector value (0(MSB), 0, 0, 1(LSB)) is assigned to a latch unit (LAT) 305 beforehand as an initial input vector value (LATINI) 304. This operation corresponds to Step 2 of the operational algorithm referred to in the above described principle of the preferred embodiment according to the present invention. As a result, the LATOUT 306 indicates the initial vector value (0(MSB), 0, 0, 1(LSB)) (a decimal value "1"). Accordingly, also the MULOUT 308 output from the square operation unit (MUL) 307 will become the same vector value. Next, the fourth bit value (the value of the MSB) of the binary value (0(MSB), 1, 1, 0(LSB)) of the SFTVAL 301 is output from the SREG 302 to a switch (SW) 303 as a select control signal (SEL) 311 at timing T0 shown in Fig. 11. This operation corresponds to Step 3 of the operational algorithm referred to in the above described principle of the preferred embodiment according to the present invention. The SW 303 selects the MULOUT 308 if the value of the SEL 311 is "0". The vector value (0(MSB), 0, 0, 1(LSB)) of the MULOUT 308, which is selected by the SW 303, is latched by the LAT 305 at timing T1 as shown in Fig. 11. This operation corresponds to the process performed when the target bit is "0" in Step 4 of the operational algorithm referred to in the above described principle of the preferred embodiment according to the present invention. Consequently, the LATOUT 306 indicates the vector value (0(MSB), 0, 0, 1(LSB)) (a decimal value "1"). Accordingly, also the MULOUT 308 output from the MUL 307 will become the same vector value. This operation corresponds to Step 5 of the operational algorithm referred to in the above described principle of the preferred embodiment according to the present invention. Next, the contents of the SREG 302 are shifted, and the third bit value "1" of the binary value (0(MSB), 1, 1, 0(LSB)) of the SFTVAL 301 is output from the SREG 302 to the SW 303 as the SEL 311 at timing T2 as shown in Fig. 11. This operation corresponds to Step 6 of the operational algorithm referred to in the above described principle of the preferred embodiment according to the present invention. If the value of the SEL 311 is "1", the SW 303 selects the DBLOUT 310. Here, the double operation unit (DBL) 309 outputs the vector value (0(MSB)), 0, 1, 0(LSB)) as the DBLOUT 310 by performing the 1-bit-left shift operation for the vector value (0(MSB), 0, 0, 1(LSB)) of the MULOUT 308. Then, the above described vector value of the DBLOUT 310, which is selected by the SW 303, is latched by the LAT 305 at timing T3 as shown in Fig. 11. This operation corresponds to the process performed when the value of the target bit is "1" in Step 4 of the operational algorithm referred to in the above described principle of the preferred embodiment according to the present invention. As a result, the LATOUT 306 indicates the vector value (0(MSB), 0, 1, 0(LSB)). Accordingly, the MULOUT 308 output from the MUL 307 indicates the vector value (0(MSB), 1, 0, 0(LSB)). This operation corresponds to Step 5 of the operational algorithm referred to in the above described principle of the preferred embodiment according to the present invention. Next, the contents of the SREG 302 are shifted, and the second bit value "1" of the binary value (0(MSB), 1, 1, 0(LSB)) of the SFTVAL 301 is output from the SREG 302 to the SW 303 as the SEL 311 at timing T4 as shown in Fig. 11. This operation corresponds to Step 6 of the operational algorithm referred to in the above described principle of the preferred embodiment according to the present invention. If the value of the SEL 311 is "1", the SW 303 selects the DBLOUT 310. Here, the double operation unit (DBL) 309 outputs the vector value (1(MSB), 0, 0, 0(LSB)) as the DBLOUT 310 by performing a 1-bit-left shift operation for the vector value (0(MSB), 1, 0, 0(LSB)) of the MULOUT 308. Then, the above described vector value of the DBLOUT 310, which is selected by the SW 303, is latched by the LAT 305 at timing T5. This operation corresponds to the process performed when the target bit is "1" in Step 4 of the operational algorithm referred to in the above described principle of the preferred embodiment according to the present invention. Consequently, the LATOUT 306 indicates the vector value (1(MSB), 0, 0, 0(LSB)). Accordingly, the MULOUT 308 output from the MUL 307 indicates the vector value (1(MSB), 1, 0, 0(LSB)). This operation corresponds to Step 5 of the operational algorithm referred to in the above described principle of the preferred embodiment according to the present invention. Then, the contents of the SREG 302 are shifted, and the first bit (LSB) value "0" of the binary value (0(MSB), 1, 1, 0(LSB)) of the SFTVAL 301 is output from the SREG 302 to the SW 303 as the SEL 311 at timing T6 as shown in Fig. 11. This operation corresponds to Step 6 of the operational algorithm referred to in the above described principle of the preferred embodiment according to the present invention. If the value of the SEL 311 is "0", the SW 303 selects the vector value (1(MSB), 1, 0, 0(LSB)) of the MULOUT 308. Then, the above described vector value of the MULOUT 308, which is selected by the SW 303, is latched by the LAT 305 at timing T7 as shown in Fig. 11. This operation corresponds to the process performed when the target bit is "0" in Step 4 of the operational algorithm referred to in the above described principle of the preferred embodiment according to the present invention. As a result, the LATOUT 306 indicates the vector value (1(MSB), 1, 0, 0(LSB)). Then, the vector value (1(MSB), 1, 0, 0(LSB)) of the LATOUT 306 is obtained as the TAPINFO 1205 shown in Fig. 4, which corresponds to the amount of phase shift d equal to 6, in synchronization with the termination of the eighth timing, that is, the timing T7. This vector value matches the result of the operation α 6 =α 3 +α 2 shown in Fig. 9(c), or the result of the conventional operation represented by the equation (13). <The First Configuration of the MUL 307 According to the Preferred Embodiment of the Present Invention> Fig. 12 shows the first configuration of the MUL 307 shown in Fig. 10. This configuration implements the following connections between the input signal lines a 0 through a 3 (the LATOUT 306 shown in Fig. 10) and the output signal lines b 0 through b 3 (the MULOUT 308 shown in Fig. 10). (1) Normally, if i≤[n/2] (note that [x] is an integer which does not exceed x) for an input signal line a i (0≤i≤n-1), the input signal line a i is connected- to an output signal line b 2i . Specifically, the input signal line a 0 is connected to the output signal line b 0 via one of exclusive-OR circuits (EXORs) 501, while the input signal line a 1 is connected to the output signal line b 2 via an EXOR 501. (2) Normally, if i>[n/2] for the input signal line a i (0≤i≤n-1), the input signal line a i is connected to one or more signal output lines b j (note that o≤j≤n-1), which are equivalent to α 2 i determined by the primitive polynomial f(x)=x 4 +x+1 within the Galois field GF (2 n ). Specifically, the input signal line a 2 is connected to the output signal line b 0 via the output signal line b 1 and the EXOR 501, while the input signal line a 3 is connected to the output signal line b 2 via the output signal line b 3 and an EXOR 501. If a plurality of input signal lines are connected to one output signal line as described above, the plurality of input signal lines are connected to the one output signal line via an EXOR 501. With the above described configuration, the MUL 307 for performing the square operation within the Galois field GF (2 n ) can be implemented as the hardware running with 1 clock. <The Second Configuration of the MUL 307 According to the Preferred Embodiment of the Present Invention> Fig. 13 shows the second configuration of the MUL 307 shown in Fig. 10. This configuration comprises an input register 601 for holding the respective signal values of input signal lines a 0 through a n-1 (n=4 in Fig. 13); a tap table 602 for storing the coefficient value indicating each of the operation values α 0 , α 2 , ···, α 2n , which corresponds to the primitive element α within the Galois field GF (2 n ); a latch 603 for holding the respective n-bit signal values b 0 through b n-1 ; and n exclusive-OR circuits (EXORs) 604 for performing an exclusive-OR operation between each output value of the latch 603 and that of the tap table 602, and for feeding back the result of the operation to the input side of the latch 603. When operations are started, the entire contents of the latch 603 are reset to 0. Next, the following operations are repeatedly performed sequentially from the value of the input signal line a 0 to that of the signal line a n-1 , among the values of the input signal lines a 0 through a n-1 (the LATOUT 306 shown in Fig. 10), which are assigned to the input register 601. (1) The coefficient value group stored in the tap table 602, which corresponds to the value of the input signal line a i currently being processed, is read out, and the respective coefficient values are input to the respective n EXORs 604. The respective outputs of the latch 603 are also input to the respective EXORs 604. (2) The respective outputs of the EXORs 604 are fed back to the respective inputs of the latch 603. (3) The next input signal line a i+1 is selected, and the operations of the above described (1) and (2) are again performed. (4) Upon termination of the operations of (1) and (2) for the input signal line a n-1 , the outputs of the respective EXORs 604 are output as the output signals b 0 through b 3 (the MULOUT 308 shown in Fig. 10). With the above described configuration, the MUL 307 which performs the square operation within the Galois field GF (2 n ) shown in Fig. 10 can be implemented as the hardware running with a maximum of n clocks. <The First Configuration of the DBL 309 According to the Preferred Embodiment of the Present Invention> Fig. 14 shows the first configuration of the DBL 309 shown in Fig. 10. This configuration implements the following connections between the input signal lines a 0 through a 3 (the MULOUT 308 shown in Fig. 10) and the output signal lines b 0 through b 3 (the DBLOUT 310 shown in Fig. 10). (1) Normally, if i<n-1 for an input signal line a i (0≤i<-1), the input signal line a i is connected to an output signal line b i+1 . Specifically, the input signal lines a 0 , a 1 a 2 are respectively connected to the output signal lines b 1 , b 2 , b 3 . (2) Normally, if i=n-1 for the input signal line a i (0≤i≤n-1), the input signal line a i is connected to one or more output signal lines b j determined by the primitive polynomial f(x)=x 4 +x+1 within the Galois field GF (2 n ). Specifically, the input signal line a 3 is connected to the output signal line b 1 via an EXOR 701 and the output signal line b 0 . If a plurality of input signal lines are connected to one output signal line as described above, the plurality of input signal lines are connected to the output signal line via the EXOR 701 as described above. With the above described configuration, the DBL 309 for performing the double operation within the Galois field GF (2 n ) shown in Fig. 10 can be implemented as the hardware running with one clock. <The Second Configuration of the DBL 309 According to the Preferred Embodiment of the Present Invention> Fig. 15 shows the second configuration of the DBL 309 shown in Fig. 10. This configuration comprises one or more exclusive-OR circuits (EXORs) 802, each of which performs an exclusive-OR operation between the outputs of one or more predetermined stages, which correspond to the positions of one or more coefficients equivalent to α n determined by the primitive polynomial f(x)=x 4 +x+1 within the Galois field GF (2 n ), and the output of the last stage of the SR 801, and for outputting each operation result to the stage next to each predetermined stage as each input signal, in addition to the n-stage (n=4 in Fig. 15) shift register SR 801. In the example shown in Fig. 15, an EXOR 802 is inserted between the first and second stages a 0 and a 1 of the SR 801. On starting operation, an input signal (the MULOUT 308 shown in Fig. 10) is input to each of the stages of the SR 801, and the content of each of the stages of the SR 801 is output as an output signal (the DBLOUT 310 shown in Fig. 10) after the shift operation is performed for one stage. With the above described configuration, the DBL 309 for performing a double operation within the Galois field GF (2 n ) shown in Fig. 10 can be implemented as the hardware running with one clock.

Claims: An M sequence phase shift coefficient calculating method for calculating respective phase shift coefficients b 0 through b n-1 intended for obtaining an M sequence x d a i whose phase is shifted by d bits from an output a i of an M sequence generated by a primitive polynomial f(x) of degree n by using a linear combination of respective M sequences x 0 a i through x n-1 a i , whose phases are shifted by 0 to n-1 bits from the output a i , comprising: a first step of inputting a binary value of an amount of phase shift d; a second step of setting an n-bit vector value, which corresponds to a decimal value "1", as an initial input vector value; a third step of setting a target bit as a most significant bit of a binary value of the amount of phase shift d input in said first step; a fourth step of multiplying an input vector value with a vector value of a primitive element α within a Galois field GF (2 n ) and using a multiplication result as an output vector value if the target bit is "1", and of defining the input vector value as the unchanged output vector value if the target bit is "0"; a fifth step of performing a square operation within the Galois field GF (2 n ) for the output vector value obtained in said fourth step; a sixth step of shifting a position of the target bit to a least significant bit side by 1 bit in the binary value of the amount of phase shift d, which is input in said first step, and of performing said fourth and fifth steps by using a result of the square operation performed in said fifth step as a newly input vector value; and a seventh step of outputting respective elements of operation results which are obtained so far as the respective phase shift coefficients b 0 through b n-1 , when said fourth and fifth steps has been performed for all of bits structuring the binary value of the amount of phase shift d input in said first step. An M sequence phase shift coefficient calculating device for calculating respective phase shift coefficients b 0 through b n-1 intended for obtaining an M sequence x d a i whose phase is shifted by d bits from an output a i of an M sequence generated by a primitive polynomial f(x) of degree n by using a linear combination of respective M sequences x 0 a i through x n-1 a i , whose phases are shifted by 0 to n-1 bits from the output a i , comprising: an n-stage shift register circuit (302), to which an n-bit binary value corresponding to an amount of phase shift d is assigned, for performing a shift operation in synchronization with a clock; and outputting a most significant bit side output as a switch control signal; a latch circuit (305), to which an n-bit vector value corresponding to a decimal value "1" is assigned as an initial vector value, for thereafter sequentially holding an n-bit input in synchronization with the clock; a square operation circuit (307) for performing a square operation within a Galois field GF (2 n ) for an n-bit output of said latch circuit; a double operation circuit (309) for performing a double operation within the Galois field GF (2 n ) for an n-bit output of said square operation circuit; and a switch circuit (303) for selecting either of the n-bit output of said square operation circuit and the n-bit output of said double operation circuit according to the switch control signal, and for outputting a selected output to said latch circuit, wherein the n-bit output of said latch circuit is output as the respective phase shift coefficients b 0 through b n-1 after a shift operation and a latch operation are performed a number of times n. The M sequence phase shift coefficient calculating device according to claim 2, wherein said square operation circuit (307) is configured so that: between n-bit input signal lines a 0 through a n-1 and n-bit output signal lines b 0 through b n-1 , an input signal line a i is connected to an output signal line b 2i if i is equal to or smaller than a maximum integer which does not exceed n/2; the input signal line a i is connected to one or more output signal lines b j (0≤j≤n-1), which are equivalent to a value α 2i determined by a primitive polynomial f(x) possessed by a primitive element a within the Galois field GF (2 n ) if i is larger than the maximum integer which does not exceed n/2; and a plurality of input signal lines are connected to a single output signal line via an exclusive-OR circuit (501) if the plurality of input signal lines are connected to the single output signal line. The M sequence phase shift coefficient calculating device according to claim 2, wherein said square operation circuit (307) comprises: an input register circuit (601) for holding respective values of n-bit input signal lines; a tap table storing circuit (602) for storing coefficient values indicating respective operation values α 0 , α 2 , ···, α 2n , which are determined by the primitive polynomial f(x) possessed by a primitive element α within the Galois field GF (2 n ); a latch circuit (603) for holding respective signal values of n bits; and n exclusive-OR circuits (604), each of which performs an exclusive-OR operation between each output value of said latch circuit (603) and each output value of said tap table storing circuit (602), and feeds back each operation result to an input side of said latch circuit (603), wherein: after entire contents of said latch circuit (603) are reset to "0" when operations are activated, a coefficient value group stored in said tap table storing circuit (602), which corresponds to signal values of input signal lines currently being processed, is read out; respective coefficient values are input to said n exclusive-OR circuits (604); and an operation for making said latch circuit (603) perform the latch operation is sequentially and repeatedly performed for respective signal values of input signal lines, which are held in said input register circuit, starting from a signal value of an input signal line of a least significant bit to a signal value of an input signal line of a most significant bit. The M sequence phase shift coefficient calculating device according to claim 2, wherein said double operation circuit (309) is configured so that: between n-bit input signal lines a 0 through a n-1 and output signal lines b 0 through b n-1 , an input signal line a i is connected to an output signal line b i+1 if i is smaller than n-1; the input signal line a i is connected to one or more output signal lines b j (0≤j≤n-1) equivalent to a value α n determined by a primitive polynomial f(x) possessed by a primitive element α within the Galois field GF (2 n ) if i equals n-1; and a plurality of input signal lines are connected to a single output signal line b j (0≤j≤n-1) via an exclusive-OR circuit (701) if the plurality of input signal lines are connected to the single output signal line. The M sequence phase shift coefficient calculating device according to claim 2, wherein said double operation circuit (309) comprises: an n-stage shift register circuit (801); and one or more exclusive-OR circuits (802) for performing an exclusive-OR operation between outputs of one or more predetermined stages of said shift register circuit, which correspond to positions of one or more coefficients equivalent to a value α n determined by a primitive polynomial f(x) possessed by a primitive element α within the Galois field GF (2 n ), and an output of a last stage of said shift register circuit, and for outputting each operation result as each input signal to a stage next to each of the predetermined stages, wherein when operations are started, an n-bit input signal is input to the respective stages of said shift register circuit, and contents of the respective stages of said shift register circuit are output as an n-bit output signal.

Inventors: KAWABATA KAZUO, NAKAMURA TAKAHARU, OHBUCHI KAZUHISA, KAWABATA, KAZUO, NAKAMURA, TAKAHARU, OHBUCHI, KAZUHISA, Kawabata, Kazuo c/o Fujitsu Limited, Nakamura, Takaharu c/o Fujitsu Limited, Ohbuchi, Kazuhisa c/o Fujitsu Limited
Applicants: FUJITSU LTD, FUJITSU LIMITED
Application year: 1998
Main classification: H04B1
IPC codes: H03K3, H04L9, G06F7, H04J13, G06F17, H04B1