An n stage shift register can be operated as a digital pulse sequence generator by feeding back to the first stage the modulo 2 sum of the signals produced by at least two stages of the register. By properly selecting the stages from which the modulo 2 sum feedback signal is derived, the register can be made to produce a pseudo random sequence of 2n-1 digits (the maximum length sequence possible for n stages). This application relates to a machine implemented process for calculating from a specific feedback connection producing a maximum length sequence all other feedback connections that will produce a maximum length sequence.

1. A machine implemented method for selecting which of the n stages of a shift register generator should be summed modulo-2 (in addition to the last stage) to produce a feedback signal to the first stage of the shift register to produce an output pulse sequence of maximum length of 2n-1 bits when one such feedback combination of stages is known, comprising the steps of: deriving a characteristic equation, a polynomial of degree n-1, the variable of which is the response matrix of the known feedback combination and the coefficients of whch are 1 or 0 depending on whether the stage represented by the exponent of the variable is summed or not, the zero power of the variable being an identity matrix with a coefficient of one; generating a first matrix having n columns and 2n-1 rows, the columns representing the exponents of a polynomial of degree n1 and the rows representing successive powers of the response matrix of the known combination, the elements of which are 1 or 0 depending on whether the powers of the variable corresponding to the column number is summed or not to equal the power of the variable corresponding to the row number, where the characteristic equation is substituted for the n-th power of the variable where it appears; generating a second matrix of n rows and n-1 columns, the successive rows of which are every k mod (2n-1)-th row of the first matrix with the column corresponding to the zero power of the variable eliminated, where k is an odd prime number that is also prime to 2n-1; the known solution, all the solutions have been found. This is also known because which is Euler''s Function. reducing the second matrix by adding modulo-2 one column to another so that each row except the last has only one element with a value of one and each such element in a differeent column for each row; and selecting as the feedback stages, in addition to the last stage, those stages corresponding to the row number having an element with a value of one in a column in which the last row has an element with a value of one.
2. A machine implemented method for finding which of n stages in a shift register should be summed modulo-2 to the last stage to provide a feedback signal to the first stage of the shift register to generate an output from one stage of a maximum length sequence of 2n-1 bits when one such feedback combination of stages is known, comprising the steps of: deriving a characteristic equation vector having elements aj such that

United States Patent [72] Inventors Edward Charles Farnett Cinnaminson; Lee Oliver Upton, Jr., Edgewater Park, both of NJ. 21 Appl. No. 880,028 [22] Filed Nov. 26, 1969 [45] Patented Oct. 19,1971 [73] Assignee RCA Corporation [54] MAXIMUM LENGTH PULSE SEQUENCE GENERATORS 2 Claims, 5 Drawing Figs.
[52] US. Cl 235/152, 340/146. 1, 340/348 [51] Int. Cl G06f1/02 [50] Field ofSearch 235/152, 153; 340/146.1, 348; 331/78 [56] References Cited UNITED STATES PATENTS 3,l 19,097 l/1964 Tullos 340/l46.1 X 3,155,818 11/1964 Goetz 340/146.1X
SET vi moor s msrs p PRIME NUMBEP ALSO PRIME T0 2M (ii CHARACTERlSTIC EQUATION VECTOR m NO OF SOLUTIONS PRESET Primary Examiner-Charles E. Atkinson Attorney-H. Christoffersen ABSTRACT: An n stage shift register can be operated as a digital pulse sequence generator by feeding back to the first stage the modulo 2 sum of the signals produced by at least two stages of the register. By properly selecting the stages from which the modulo 2 sum feedback signal is derived, the register can be made to produce a pseudo random sequence of 2"l digits (the maximum length sequence possible for n stages). This application relates to a machine implemented process for calculating from a specific feedback connection producing a maximum length sequence all other feedback connections that will produce a maximum length sequence.
STAGE STAGE STAGE STAGE 1 OUTPUT PATENTEDUCT 19 l97| SHEET 1. UF 3 SET= REPEAT Fig.1.
Edward C. Fame and Lee 0. Upf0n,Jr.'
ATTORNEY PATENTEUHET 19 l97l 3 6 14,400
SHEET 20F 3 5| PROGRAM CONTROL UNIT 50 52 is E PROGRAMMABLE TAPE GENERAL PURPOSE STATIONS COMPUTER 54 53 /PR|NTER E CONTROL CONSOLE CARD READER Flg. 3. PRESET IO 30 40 4 STAGE STAGE E STAGE STAGE 1 2 3 4 OUTPUT Fig. 4. PRESET STAGE STAGE E STAGE STAGE 1 2 3 4 OUTPUT INVIiN'IuRi Flg. Edward C. Farnett and Lee 0. Upton, Jr. BKJJ ATTORNEY PAIENTEDum 19 IHYI SHEET 3 [1F 3 kkJFI KkJi' kkJ FINISHED Fig.2.
l REPEAT 1mm).
NVIZN'IYJRS Edward C. Farrnetl and L66 0. Upton, Jr.
A TTORNE Y MAXIMUM LENGTH PULSE SEQUENCE GENERATORS BACKGROUND OF THE INVENTION The invention described in this application was made in the course of a contract with the Department of the Navy.
Pseudorandom sequences of digital pulses are useful in communications and radar. For example, modulating the output of a radar with such a sequence of pulses improves target resolution without increasing power output. Target resolution (the ability to separate targets) is determined by the width of the individual pulses, whereas the power output is determined by the total duration of the pulse widths transmitted.
A pulse sequence is a serially occurring train of pulses having a pseudorandom nature in that the time between pulses varies in an almost random way. The pulse sequence can be represented by a series of 1's and Os, where the digit 1 indicates a pulse at the time interval represented by the position of the digit and a zero indicates the absence of a pulse.
One of the simplest circuits for generating a pulse sequence uses a number of shift register stages connected serially so that the output of each stage is the input for the following stage and the output of the last stage is the input for the first. The number of stages required is the total number of digits (1s and (TS) in the pulse sequence. The stages are preset to represent the sequence desired and the shift register is shifted continuously. The output sequence can be taken from the output of any stage. This method is limited by space and cost consideration as the number of pulses in the sequence increases.
A modification of the simple shift register generator permits a sequence length of 2"-l pulses using only n stages, the sequence of all zeroes being excluded. For example, a 20- stage shift register properly connected can repeat a given sequence of more than one million pulses.
The shift register is modified by disconnecting the output of the last stage from the input of the first and providing an input value of l to the first stage if an odd number of certain selected stages are storing a value of l. The output pulse sequence is taken from the output of the last stage.
The output of the certain selected stages (feedback taps) are used as inputs to a modulo-2 adder, the output of which is the input to the first stage. The output of the last stage must always provide an input to the modulo-2 adder; otherwise, the effect is to shorten the length of the register.
Only certain feedback tap combinations will produce a pulse sequence of maximum length. Other combinations will produce shorter sequences. Because the last stage must always be used as a feedback tap, there is a total of 2" different possible combinations of feedback taps.
At least one combination of feedbacks is known in the art for n stages up to n==l that will produce a maximum length sequence. Other maximum length combinations can be derived by trial and error methods from a known combination.
One such method is based on the fact that taking every k-th pulse produces a different sequence, where k is an odd number that is prime to 2"l. Some of the 2'"! combinations can be eliminated and the remaining combinations are tested to find one that will produce a maximum length sequence.
It is an object of this invention to describe a method whereby all possible combinations of feedback taps producing a maximum length sequence can be directly determined, given one such combination.
BRIEF DESCRIPTION OF THE INVENTION A machine implemented method is described for selecting which stages of an n stage shift register generator, in addition to the last stage, should be summed modulo-2 to produce a feedback signal to the first stage of the shift register so that the output of the last stage will be a pulse sequence of maximum length (2"1 bits) if one such feedback combination is known.
The first step is the derivation of a characteristic equation, a polynomial of degree n--l, the variable of which is the response matrix of the known combination. The coefficients are one or zero depending on whether the stage represented by the exponent of the. variable is summed or not; the zero power of the variable is an identity matrix and has a coeffcient of one.
The next step is the generation of a first matrix having n columns and 2"-l rows. The columns represent theexponents of a polynomial of degree n-l and the rows represent successive powers of the response matrix of the known combination. The elements of the first matrix are one or zero depending on whether the powers of the variable corresponding to the column number, are summed or not to equal the power of the variable corresponding to the row number. If the n-th power appears, it is replaced by the characteristic equation.
Next, a second matrix is generated having n rows and n-l columns. The successive rows of the second matrix are every k mod (2"l)-th row of the first matrix with the column corresponding to the zero 'pow'er ofthe variable eliminated.
The second matrix is then reduced by adding columns together so that each row execpt the last has only one element with a value of one and each such element in each row is in a different column. The last step is the selection as a feedback stage, in addition to the last stage, those stages corresponding to a row number .having an element with a vlaue of one in a column in which the last row has an element with a value of one.
BRIEF DESCRIPTION OF THE DRAWINGS FIG. 1 is a flow chart of the steps to be followed in the execution of the first part of a program which is the preferred embodiment of the invention.
FIG. 2 is a flow chart of the second part of the program which is the preferred embodiment of the invention.
FIG. 3 is a block diagram of a general purpose programmable computer comprising the apparatus for performing the invention.
FIG. 4 is a logic diagram of a four-stage shift register generator providing a maximum length pulse sequence.
FIG. 5 is a logic diagram of another four-stage shift register generator providing a maximum length pulse sequence.
DETAILED DESCRIPTION OF THE INVENTION The method of the invention uses Boolean matrices which are manipulated by modulo-2 operations. A Boolean matrix is a matrix in which the elements are zero or one. Boolean matrices can be manipulated using the logical operators AND, OR and NOT. Modulo-2 addition is an EXCLUSIVE 0R function. That is, (a+b) mod 2 a@b= aB-l-Eb. The same result can be reached using conventional arithmetic by The right-hand member of the equation above has the value of 1 if, and only if, a or b has a value of l, but not both.
Addition (mod-2) of Boolean matrices is possible only when the matrices have the same number of rows and the same number of columns. The usual convention of denoting a matrix by an upper case letter, e.g., A, and its elements by a subscripted lower case letter, e.g., a will be followed in this description. The subscripts i and j denote the row and column numbers, respectively. Therefore, A638 a @b for i l,2,..., m andj= l,2,3,...,n.
A property of modulo-2 matrix addition is A Q A 0.
Also, the values of matrices are unchanged when columns are added to columns or rows are added to rows.
Multiplication (mod-2 of Boolean matrices is similar to the usual method of multiplying matrices, the summation of products for each element being modulo-2. The number of columns in the first matrix must be equal the number of rows in the second, and the number of rows in the first must equal the number of columns in the second. The elements of the product matrix C that result from A 93? B are 1n the method taught by this invention, a vector is used which is a matrix having m column and row one.
The first step in practicing this invention is the derivation of a characteristic equation vector from the known feedback combination. For a shift register generator having it stages, the characteristic equation vector has n 1 columns. The first and last columns have a value of 1. The second column of the vector corresponds to the first stage of the shift register, the third column corresponds to the second stage of the shift register, and so on up to the n 1 column corresponding to the n stage of the shift register. if the output of a stage is used as an input to the feedback network in the known shift register generator, a value of 1 is inserted in the characteristic equation vector column corresponding to that stage. if the output of a stage is not summed, the value of the element in the corresponding column is 0.
The characteristic equation vector is a function of the sum of matrices derived from the known feedback combination. The known shift register generator can be used to derive A matrix. The rows of the A matrix correspond to the inputs of each stage and the columns, to the outputs. if the j-i-th stages output signal is the input signal for the i-th stage, a, ,is given the value of 1. Otherwise, it is given a vlaue of 0. Such a matrix is called a response matrix.
if various powers of A are taken, i.e., A (which is A @A), A, and so on, it will be found that the n-th power is equal to the sum of the powers corresponding to the feedback. stages plus the identity matrix, A". The identity matrix has the elements when vie-j am 1 when 1 2 in other words, the elements of the main diagonal (upper left to lower right) have a value of 1 and the rest, 0.
The first column of the characteristic equation vector corresponds to the A" matrix; column 2, to A; and so on.
As an illustration, if, in a 15 stage shift register generator, feedback taps are taken from stages 3, 8, l2, and 15, the characteristic equation is A =A EB A 61 A 62 A and the characteristic equation vector is 1001000010001001.
it should be noted that the vector is used for programming purposes and has no significance, per se.
A B mitrix is constructed based on the characteristic equation vector and powers of the A matrix. The B matrix in the described embodiment has n 1. columns and 2"-l rows. The B matrix can be considered as having an infinite number of rows in which every group of 2"l rows are repeated. it is, therefore, necessary to calculate only the first 2"1 rows. The elements b in the first n rows of the B matrix are 0 except in the columns where =j-l which are given a value of 1. Each of the remaining rows of the B matrix is derived from its preceding row by shifting each column one column to the right and setting the value of the element in the first column to 0. if the element in the n 1 column of a row has a value of l, the characteristic equation vector is added (mod-2) to that row.
The B matrix elements are defined as foloows. The row number i, of the B matrix corresponds to the i-th power of the A matrix. The column numbers of the B matrix correspond to the powers of the A matrix (including the zero power), the summation of which equal the power denoted by the row number.
For example, in the first row of the B matrix, which corresponds to the A (A) matrix, b =l and the other elements are 0 because A is a function of A only. Similarly, in row 2, b l and the other elements are 0. in the n-th row, b,, ,=l.
The rule stated above requires that the characteristic equation vector be added to that row. in effect, this substitutes the summation representing the characteristic equation of A".
In the example above, where n=l5, the th row of the B matrix would be represents lll A15 @A11@ A8 The element in the n 1 column has the value of 1 so the characteristic equation vector is added to give 100000101001 1000 as the 18th row. This in effect substitutes the characteristic equation for A so that l8 12@ HG 9A3@ 0 QAIIQAO @A8 As stated above, .4 Q3 A =0 so that row 18 represents AIK=A12@AU 8 46 @A0 The 19th row represents A=A Q A" Q A Q .4 Q3 A. That is, each power is increased by multiplication by A. Therefore, the 19th row is 0100000101001100.
This process of generating rows is continued through row 2"-1, which completes the B martix.
A C matrix is then derived from the B matrix by selecting as each row of the Cmatrix every k-th row of the B matrix, where k is a prime number other than 2 and is also prime to 2"1. The C matrix has n rows and 011-1 columns; the first and last columns of the B matrix rows are dropped when deriving the C matrix.
The C matrix is then reduced by adding columns so that each row of the C matrix except the last has a single element with a value of 1, each in a different column.
A straightforward method of reducing the C martix is to add (mod-2) the first column having an element with a value of l in the first row to all other columns having a value of 'l in the first row. All the elements in the first row will then have a value of 0 except that in the column which was added to the others.
Next, the first column having an element with a value of l in the second row and a value of 0 in the first row is added to all other columns having an element with a value of 1 in the second row. This process is repeated for each row in succession except the last. The new solution uses as feedback stages those stages corresponding to row numbers having an element with a value of l in a column in which the last row has an element with a value of 1.
Another solution can be found by repeating the steps described above using a characteristic equation vector derived from the solution just obtained and using the same value of k. When the original solution is repeated, all the possible combinations of feedback taps producing a maximum length pulse sequence have been found.
HO. 1 shows a flow chart for the first part of a program for performing the method of the invention. The first step as shown in the block it) sets up the initial parameters of the problem. The value of n is stored to indicate the number of stages, p is selected from the prime numbers (except 21) which are also prime to 2"1, the characteristic equation vector, a,, is determined by the known solution, and the number of solutions is denoted by m.
The next step as shown in the block sets an indexing variable i to the value of l and the second step as shwon in the block 16 sets another indexing variable j to the value of 1.
During the running of the program, the variables may be addressed using any one of several conventions. ln this embodiment, the values of b are stored in a sector of the memory in serial fashion so that the first subscript varies faster than the second. For example, b, is the first element stored in the assigned sector and is followed in succession by b b and so on. Referencing a particular element, denoted by the subscripts i and j, causes the computer to calculate an address based on the values of i and j. if L is the memory address of the location preceding that of the element bu. the memory address of b is found by L i rU-l where r is the number of rows. For purpose of clarity, this step has not been shown, but it is implied by denoting a subscripted variable.
Next, each element in the first n-l rows of the B matrix is set to 0 as is shown in the block 20 in FIG. 1.
A decision point in the program, indicated by a diamond shape, is shown by the block 23. It indicates that the value ofi is compared to the value of j-l. If the values are equal, the program takes the branch indicated by the arrow emerging from the bottom of the block and labelled with an sign. If i is greater or less than j-l, the program branches to the other path.
In FIG. 1, if the value of i is equal to j-l, then the element b, is given the value of I. Ifi is not equal to j-l b remains 0 by skipping the b ,=l step.
Next, it is determined whether the value of j exceeds the quantity n I. If it does not, j is replaced by the value of j 1 and the program returns to the step indicated by the block 20. The steps are repeated as described above. When the value of j is greater than the quantity n l, the value of i is examined as is shown in the block 27 to determine whether it is greater than the quantity n I. If it is not, i is replaced by the value of i+ l and the steps of the program are repeated starting at the point where the indexing variable j is set to a value of l as is shown in the block 16. The program continues in this manner until the value ofi exceeds the quantity n-l.
The effect of manipulating the values of the indexing variables i and j as shown in the flow chart of FIG. 1 from the block to the block 27 is that i and j are stepped through all the (n-l) (n+1) combinations. That is to say, the elements of each successive row of the matrix are addressed by incrementing the value of i by l and each successive column of each successive row is addressed by incrementing the value of j and 1. Each time a new row is designated by incrementing the value of i, the value of j is initialized to the value of l and incremented through its range before the value of i is again incremented by 1. Therefore, all the elements in the n+1 columns of the first n-l rows (i.e., n -l elements) are set to 0 except for those in which the value of i is equal to j-l Next, the value of the indexing variable i is set to the quantity n as is shown in the block 31. The first element of the n-th row is set to 0. The index variable j is set to the value of 2 and each element of the n-th row is then set to the value of the element in the preceding column of the preceding row. This is shown in the block 33 by subtracting 1 from each of the indicies on the right-hand side of the equation.
This is done for each column of the n-th row as indicated by incrementing the indexing variable j from a value of 2 to n l.
The element in the last column (b is then examined to determine whether it is equal to 0 as is shown in the block 37.
[f the element is equal to 0, then the indexing variable i is examined to see whether the last row (2"l) has been completed. If not, the value of i is incremented by 1 and the next row is processed.
If the element is equal to l, the characteristic equation vector is added to that row as is shown by the block 38, using the indexing variable j to process all the columns.
The elements of the characteristic equation vector are denoted by a, where j is the column number. The elements of the characteristic equation vector are added (mod'2) to the elements of the i-th row by taking the absolute value of their difference. The indexing variable j is set to l and incremented through its range to the value n 1 so that each column of the i-th row is added (mod-2) in turn to the corresponding column of the characteristic equation vector.
The value of the indexing variable i is checked to determine whether it exceeds the value of 2"l. If not, its value is incremented by 1 and the next row is processed.
When the (Z -1 )th row has been processed, the B matrix is complete.
The C matrix is derived by taking every k-th row of the B matrix (excluding the first and last column) as successive rows in the C matrix. The indexing variables iand j are set to values of l and 2, respectively, and the value of k is determined by multiplying the value of i, which now denotes the row of the C matrix being processed, by the value of p which was set as an input parameter.
The value of k is next compared to the value 2"l. If k is greater than 2-l the quantity 2"-l is subtracted from k and the comparison repeated. This permits the program to treat the B matrix as if it had an infinite number of rows by keeping the value of k between 1 and 2"l, inclusive, i.e., k mod (2" I), letting 0 mod (2"-1 =2"-I.
When k is less than or equal to 2"l, each element of the B martix in the k-th row from columns 2 through it is transferred to the corresponding columns of the i-th row of the C matrix. (Columns 1 and n l of the C martix are ignored. Corresponding column numbers are used to simplify the program.)
When i=n, i.e., when the n-th row of the C matrix has been derived, the C matrix is complete.
FIG. 2 shows the flow chart for the final part of the program. The indexing variables i and j are set to l and 2, respectively, and each element in the i-th row of the C matrix is compared to 1. In the first row, a column is located in which the element has the value of l and this column is then added to all other columns having an element with a value of l in the same row. In the following rows, an element with a value of l is located in a column which has all elements in all of the preceding rows equal to zero. This column is added to all other columns having an element with a value of l in the same row. This procedure is repeated for all except the last row. When every row except the last has one element having a value of l in a difierent column than the other rows, the matrix C is reduced.
The C matrix is reduced in the described embodiment of the invention in the following way. Each element of the i-th row starting with the first row is examined to see if it has the value of l as is shown in the block 41.
If c, ,=l, a check is made to determine if the first row is being processed. The first row is a special case because the elements in the same column of the preceding rows are not checked. If the first row is not being processed, an indexing variable ii is set to l as is shown in the block 43. Each element in the proceeding rows of the column is checked to determine whether it has a value of 0.
If a preceding row is found that does not have the value of 0, the value of j is checked to determine whether all columns of the i-th row have been processed. If not, j is incremented by l and the next column examined.
If all of the preceding rows have a value of zero, the indexing variables ii and jj are set to the values of i and 2, respectively.
Using the indexing variables ii and 1 the j-th column is added (mod-2) to all other columns having a l in the same row. The elements in the preceding rows (1 to i-lhave been found to have a value of 0 so they need not be added.
It must be determined, however, during this addition process that the column which is being added to the other columns is not inadvertently added to itself. This is done by skipping the addition process when the indexing variable jj is equal to j, the later representing the column being added to the other columns. This is performed as is shown in the block 44. The addition process is skipped when jj=j.
The element in the ii-th row and ji-th column is checked to determine if it has a value of I. If not, the addition process is skipped. if c ,;=l, then each element of the j-th column is added (mod-2) 2) to the element in the corresponding row of the jj-th column.
The indexing variable kk is used as a row index to process all the elements of the columns involved.
When the addition has been completed for each row, the
value of i is checked to see whether it exceeds n-l. If it does shown in the block 41$. Then, the element in each column of the C matrix in the last (n-th) row is examined to determine whether it has a value of 1. if it does not, the value of j is incremented by 1 if it does not exceed the value of 12-1, and the process is repeated.
When an element in a column of the last row of the C matrix has a value of 1, each row above it is examined to determine the number of the row in which the value of the element is equal to 1. This is done by setting the indexing variable i to a value of 1 and examining the value of the element of each row.
When an element having a value of l is found in the i-th row, a value of l is inserted in the solution vector in the i-th column. This is continued until the value of the indexing variable j exceeds the value of n.
The solution vector is then written as the solution as is shown in die bloclr on. Only the values of the columns 2 through n are wirtten. The columns written correspond to each stage of the shift register generator except the last. The last stage is always used as a feedback tap. The solution vector indicates by l in the j-th column that a feedback tap is to be taken from the j-th stage of the shift register.
The value of m, which denotes the number of solutions, is decremented by l and checked to determine whether it is equal to 1. if m=0, the program is finished. if m is greater than 0, the characteristic equation is replaced by the solution vector, the first and last columns being set with a value of 1, respectively. The program is then repeated beginning at block l5 on 1F 1G. 11 to find another solution.
instead of setting a value of m, the program could store the original characteristic equation vector and compare it to each subsequently derived characteristic equation vector. When they were equal, all the solutions would have been calculated. However, as the value of it increases, the number of solutions increases. Therefore, it is preferable to set m as the desired number of solutions.
The characteristic equation itself can be derived by the program. For example, instead of using the characteristic equation vector as an input parameter, the number of the stages from which feedback taps are taken in the lcnown solution can be used as an input parameter. The program can generate the vector by setting all columns to 0 except the first and last which will be set to 1. Then the input (stage numbers) can be used as column indicies to set the vector element to a value of 1.
When the value of n is large, the 113 matrix will have a large number of rows. For example, if n=20, the matrix will have over a million rows. When n==l00, the value of 2-1 exceeds 1.25Xl0, which is a prohibitive number of rows to store.
There are techniques well known in the art for alleviating this problem. One is to use a single bit to store a matrix element value because such values are restricted to l and 0. However, for n=100, the memory requirements would still be prohibitive.
Another technique would be to store each row on magnetic tape as a block of information. The indexing variables would then be used to find the row being processed by any one of several ways. The simplest is to label each bloclr with a number identifying the row the block represents. The tape can then be operated to read each block starting with the first until the block identification number matches the index variable indicating the number of the row being processed.
A better technique, based on the fact that only n rows must be saved for the C matrix, would be to retain in memory or on tape only every k-th row. ln effect, the C matrix would be developed as the rows of the B matrix were processed. This technique can be modified for the case where !c is varied.
FIG. 3 illustrates the apparatus on which the method of the invention can be practiced. it is a block diagram of a computer system comprising a program control unit 50 to which are connected tape stations 511, a printer 52, a card reader 53, and a control console $4. A program based on the invention written in a computer language suitable for performance on the particular computer for which the program is written, is stored in a memory which may be an integral part of the program control unit 50. The program can be stored on a magnetic tape such as a program library tape and read into the program control unit 50 of the computer from tape stations 51. Alternatively, the program can be stored on punched cards and read into the computer program control unit 50 by card reader 53.
The initial data and parameters required during the execution of the program may be punched on cards and read from the card reader 53 or from the tape stations 51 as required during the execution of the program. Alternatively, the parameters might be generated by a preceding program and stored in specific locations in the memory within the process control unit 50.
As the various solutions are derived, results can be printed on the printer 52.
As an example of how this method derives a solution when one solution is known, consider the shift register generator shown in 1? l6. 4. It consists of four shift register stages and the contents of the first stage it) is shifted into the second stage and so on by means not shown. The output of the generator is taken from the output of the fourth stage 40 and the input to the first stage it is the output of an EXCLUSIVE-OR gate 61. The input signals of the EXCLUSIVE-0R gate 61 are the output signals of the third stage 30 and the fourth stage 40.
Therefore, stage 1 will have an input when a signal representing a l is stored in stage 3 and a signal representing a 0 is stored in stage 4, or vice versa.
The first stage is set to the value of 1 by means of a preset signal so that the contents of the shift register before the first shift pulse occurs can be represented by 1000. After the first shift pulse, the contents of the register can be represented by 0100. A 0 is shifted into the first stage because stage 3-30 and stage 4-40 both contained a value of 0 before the shift. Table 1 shows the contents of the register after each of the first 16 shift pulses.
TABLE] Shift Register Contents Shift Pulse Number Contents of Stage 2 o o 1 o 14 o o o 1 1s 1 0 o o 16 o 1 o o The contents of the shift register shown in FIG. 4 after the 16th shifi are seen in table 1 to be the same as the contents after the first shift. Therefore, shift pulses 16 through 30 will repeat the sequence shown in table i from pulses 1 through 15.
The output of stage 4 is the output pulse sequence of the generator, 001001101011 110, which will be repeated as long as the contents of the register are shifted.
It can be determined by means of Euler's phi function that there are two combinations of feedback taps for a four-stage shift register that will provide the maximum length pulse sequence. The second solution will now be derived from the known solution by the method of the invention.
The characteristic equation vector of the solution shown in 1 16.415 10011.
The first three rows of the B matrix are as follows:
The next rows are found by shifting each value of the preceding row to the right so that the fourth row of the B matrix would be 00001. The last column, however, contains a 1 so it is necessary to add (modulo-2) the characteristic equation vector. Therefore, the fourth row of the B matrix is 10010.
The next row of the B matrix would be 01001, obtained by shifting the preceding row one column to the right and inserting a zero in the first column. The last column, however, contains a value of 1 so that it is necessary to add the characteristic equation vector. The fifth row is 1 1010.
The complete 8 matrix is l 1000 01100 00110 10000 The next row of the 13 matrix will be the same as the first row and the matrix shown above would be repeated every 15 rows. The C matrix is derived by taking every k-th row of the 13 matrix, where k is a prime number other than 2 and also prime to 2 -1 (15 The smallest prime other than 2 also prime to 15 is 7. The C matrix is, therefore,
1 10 O1 1 1 l l 110 because the first and last columns of the rows from the B matrix are eliminated.
The C matrix is reduced first by adding (mod-2) the first column to the second column. The first step gives 100 01 l 101 100 The second step, adding (mod-2) the second column to the third column, yields 100 010 101 100 Finally, the third column is added (mod-2) to the first column to give the reduced C 'matrix The last row has a value of 1 only in the first column. The one in the first column above the last row is in the first row so the solution vector is 100.
The shift register shown in FIG. 5 is the shift register generator that corresponds to the solution vector. There is a feedback tap from stage 1, corresponding to the one in the first column of the solution vector.
A new characteristic equation vector can be derived by adding a column containing a 1 at each end of the solution vector, viz, 11001 If this new characteristic equation vector is used to find another solution in the manner just described, the resulting solution will be 001 and the resulting characteristic equation, 1001 1. Because this is the same as What is claimed is:
l. A machine implemented method for selecting which of the n stages of a shift register generator should be summed modulo-2 (in addition to the last stage) to produce a feedback signal to the first stage of the shift register to produce an output pulse sequence of maximum length of 2"-l bits when one such feedback combination of stages is known, comprising the steps of:
deriving a characteristic equation, a polynomial of degree n-1, the variable of which is the response matrix of the known feedback combination and the coefficients of whch are 1 or 0 depending on whether the stage represented by the exponent of the variable is summed or not, the zero power of the variable being an identity matrix with a coefficient of one;
generating a first matrix having I: columns and 2"l rows,
the columns representing the exponents of a polynomial of degree n-l and the rows representing successive powers of the response matrix of the known combination, the elements of which are 1 or 0 depending on whether the powers of the variable corresponding to the column number is summed or not to equal the power of the variable corresponding to the row number, where the characteristic equation is substituted for the n-th power of the variable where it appears;
generating a second matrix of n rows and n1 columns, the
successive rows of which are every k mod (2"-1)-th row of the first matrix with the column corresponding to the zero power of the variable eliminated, where k is an odd prime number that is also prime to 2"l; the known solution, all the solutions have been found. This is also known because which is Eulers Function.
reducing the second matrix by adding modulo-2 one column to another so that each row except the last has only one element with a value of one and each such element in a differeent column for each row; and selecting as the feedback stages, in addition to the last stage, those stages corresponding to the row number having an element with a value of one in a column in which the last row has an element with a value of one.
2. A machine implemented method for finding which of n stages in a shift register should be summed modulo-2 to the last stage to provide a feedback signal to the first stage of the shift register to generate an output from one stage of a maximum length sequence of 2"1 bits when one such feedback combination of stages is known, comprising the steps of:
deriving a characteristic equation vector having elements a,
such that 1 if the jth stage is summed 0 if the jth stage is not summed and a =l;
generating a first matrix having n+1 columns and 2"l rows with elements b such that for the first row b ==0 except b fl and the elements of each successive row are generating a second matrix having n-l columns and n rows, the elements c of each successive row corresponding to the elements of every k mod(2"1)-th row of the first matrix with the first and last columns deleted, where k is an odd prime number that is also prime to 2"1;
UNITED STATES PATENT OFFICE CERTIFICATE OF CORRECTION Patent N D t d October 19,
Inventor(s) Edward Charles Farnett G Lee Oliver Upton, Jr.
It is certified that error appears in the above-identified patent and that said Letters Patent are hereby corrected as shown below:
Col. 9, line 41, delete --llOOO--; line 49, "2 -1" should be -2 l-- Col. 1, line 48, "2 should be ---2 -line 56, "2 1" should be --2 Col. 2, line 18, "execpt" should be --except-- line 57, "a should be --a j-; line 59, "a b should be ---a j b j- Col. 3, line 2, "column" should be --columns---; line 23, "vlaue" should be -value-; line 47, "mitrix" should be --matrix;
line 54, j l" should be -i j-lline 60, "foloows" should be -foll0ws-- Col. 4, line 20, "martix" should be -matrix--; line 31, "martix" should be -matrix-- line 56, "(except 21)" should be -(except 2)--; line 61, "shwon" should be --shown---. Col. 6, line 10, "martix" should be --matrix--; line 12, "martix" should be --matrix---' line 60, "later" should be --latter-. Col. 7, line 21, "wirtten" should be --writtenline 28 "equal to 1" should be --equal to 0-"; line 48, "indicies" should be --indices-; line 53, "1.25Xl0" should be -1.25xlO *RM PC3-1050 {10-69) USCOMM-DC 60376-P69 V U 5 GOVERNMENT PRINTING OFFICE |9690-366-334 UNITED STATES PATENT OFFICE CERTIFICATE OF CORRECTION Patent No. 3,614,400 Dated October 19 1971 Inventor(s) Edward Charles Farnett G Lee Oliver Upton, Jr.
It is certified that error appears in the above-identified patent and that said Letters Patent are hereby corrected as shown below:
Col. 10, after line 14 insert,
-the known solution, all the solutions have been found.
This is also known because w -l) M 8 2 which is Euler's Function.-
Col. 10, line 26, "whch" should be ---which-- lines 45-51,
delete -the known solution. .which is Euler's Function.--
n Col. 11, line 1, b j b jr l 1 1,2,...,2 should be -b. -=b
Signed and sealed this 2nd day of May 1 972.
(SEAL) Attest:
EDWARD M.FLETCHER,JR. ROBERT GOTTSCHALK Attesting Officer Commissioner of Patents ORM powso (10459) USCOMM-DC 60376-F'69 U 5 GO ERNMENT PRINHNG OFFICE I969 0365334

Publication number | Publication date | Assignee | Title |
---|---|---|---|

US-3119097-A | January 21, 1964 | Jersey Prod Res Co | Electrical signal generator |

US-3155818-A | November 03, 1964 | Bell Telephone Labor Inc | Error-correcting systems |

US-3159810-A | December 01, 1964 | Sylvania Electric Prod | Data transmission systems with error detection and correction capabilities |

US-3162837-A | December 22, 1964 | Ibm | Error correcting code device with modulo-2 adder and feedback means |

Title |
---|

Chow and Davies, The Synthesis of Cyclic Code Generators, Electronic Engineering, April 1964, pp. 253 259. |

W. W. Peterson Error Correcting Codes, MIT Press and John Wiley & Sons, 1961, pp. 147, 148. |

Publication number | Publication date | Assignee | Title |
---|---|---|---|

US-6557129-B1 | April 29, 2003 | Janusz Rajski, Jerzy Tyszer, Mark Kassab, Nilanjan Mukherjee | Method and apparatus for selectively compacting test responses |

US-6539409-B2 | March 25, 2003 | Janusz Rajski, Mark Kassab, Nilanjan Mukherjee, Jerzy Tyszer | Method for synthesizing linear finite state machines |

US-5761265-A | June 02, 1998 | Board Of Regents, The University Of Texas System | Parallel architecture for generating pseudo-random sequences |

US-4890252-A | December 26, 1989 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Long period pseudo random number sequence generator |

US-8533547-B2 | September 10, 2013 | Mentor Graphics Corporation | Continuous application and decompression of test patterns and selective compaction of test responses |

US-9134370-B2 | September 15, 2015 | Mentor Graphics Corporation | Continuous application and decompression of test patterns and selective compaction of test responses |

US-7509546-B2 | March 24, 2009 | Janusz Rajski, Mark Kassab, Nilanjan Mukherjee, Jerzy Tyszer | Test pattern compression for an integrated circuit test environment |

WO-0138955-A1 | May 31, 2001 | Mentor Graphics Corporation | Procede de synthese de machines d'etats finis lineaires |

US-5412587-A | May 02, 1995 | The Boeing Company | Pseudorandom stochastic data processing |

US-7865794-B2 | January 04, 2011 | Mentor Graphics Corporation | Decompressor/PRPG for applying pseudo-random and deterministic test patterns |

US-3751648-A | August 07, 1973 | Communications Satellite Corp | Generalized shift register pulse sequence generator |

US-4125898-A | November 14, 1978 | The Singer Company | Digitally shaped noise generating system |

US-7500163-B2 | March 03, 2009 | Janusz Rajski, Mark Kassab, Nilanjan Mukherjee, Jerzy Tyszer | Method and apparatus for selectively compacting test responses |

US-6684358-B1 | January 27, 2004 | Janusz Rajski, Jerzy Tyszer, Mark Kassab, Nilanjan Mukherjee | Decompressor/PRPG for applying pseudo-random and deterministic test patterns |

US-7523372-B2 | April 21, 2009 | Janusz Rajski, Jerzy Tyszer, Nagesh Tamarapalli | Phase shifter with reduced linear dependency |

US-6708192-B2 | March 16, 2004 | Janusz Rajski, Mark Kassab, Nilanjan Mukherjee, Jerzy Tyszer | Method for synthesizing linear finite state machines |

US-7263641-B2 | August 28, 2007 | Janusz Rajski, Jerzy Tyszer, Nagesh Tamarapalli | Phase shifter with reduced linear dependency |

US-6327687-B1 | December 04, 2001 | Janusz Rajski, Jerzy Tyszer, Mark Kassab, Nilanjan Mukherjee | Test pattern compression for an integrated circuit test environment |

US-7493540-B1 | February 17, 2009 | Jansuz Rajski, Jerzy Tyszer, Mark Kassab, Nilanjan Mukherjee | Continuous application and decompression of test patterns to a circuit-under-test |

US-4188583-A | February 12, 1980 | Rca Corporation | Sampling method and apparatuses |

US-7877656-B2 | January 25, 2011 | Mentor Graphics Corporation | Continuous application and decompression of test patterns to a circuit-under-test |

US-6543020-B2 | April 01, 2003 | Janusz Rajski, Mark Kassab, Nilanjan Mukherjee, Jerzy Tyszer | Test pattern compression for an integrated circuit test environment |

US-6829740-B2 | December 07, 2004 | Janusz Rajski, Mark Kassab, Nilanjan Mukherjee, Jerzy Tyszer | Method and apparatus for selectively compacting test responses |

US-6353842-B1 | March 05, 2002 | Janusz Rajski, Jerzy Tyszer, Mark Kassab, Nilanjan Mukherjee | Method for synthesizing linear finite state machines |

US-7478296-B2 | January 13, 2009 | Janusz Rajski, Mark Kassab, Nilanjan Mukherjee, Jerzy Tyszer | Continuous application and decompression of test patterns to a circuit-under-test |

US-4769777-A | September 06, 1988 | General Electric Company | Unpredictable bit stream generator |

US-5574673-A | November 12, 1996 | Board Of Regents, The University Of Texas System | Parallel architecture for generating pseudo-random sequences |

US-9664739-B2 | May 30, 2017 | Mentor Graphics Corporation | Continuous application and decompression of test patterns and selective compaction of test responses |

US-3746847-A | July 17, 1973 | D Maritsas | Generating pseudo-random sequences |

US-8024387-B2 | September 20, 2011 | Mentor Graphics Corporation | Method for synthesizing linear finite state machines |

US-8108743-B2 | January 31, 2012 | Mentor Graphics Corporation | Method and apparatus for selectively compacting test responses |

US-2005015688-A1 | January 20, 2005 | Janusz Rajski, Jerzy Tyszer, Nagesh Tamarapalli | Phase shifter with reduced linear dependency |

US-7900104-B2 | March 01, 2011 | Mentor Graphics Corporation | Test pattern compression for an integrated circuit test environment |

US-7805649-B2 | September 28, 2010 | Mentor Graphics Corporation | Method and apparatus for selectively compacting test responses |

US-7111209-B2 | September 19, 2006 | Janusz Rajski, Mark Kassab, Nilanjan Mukherjee, Jerzy Tyszer | Test pattern compression for an integrated circuit test environment |

US-6874109-B1 | March 29, 2005 | Janusz Rajski, Jerzy Tyszer, Nagesh Tamarapalli | Phase shifter with reduced linear dependency |

US-7506232-B2 | March 17, 2009 | Janusz Rajski, Jerzy Tyszer, Mark Kassab, Nilanjan Mukherjee | Decompressor/PRPG for applying pseudo-random and deterministic test patterns |

US-7260591-B2 | August 21, 2007 | Janusz Rajski, Jerzy Tyszer, Mark Kassab, Nilanjan Mukherjee | Method for synthesizing linear finite state machines |

US-4023026-A | May 10, 1977 | International Telephone And Telegraph Corporation | Pseudo-random coder with improved near range rejection |

US-7093175-B2 | August 15, 2006 | Janusz Rajski, Jerzy Tyszer, Mark Kassab, Nilanjan Mukherjee | Decompressor/PRPG for applying pseudo-random and deterministic test patterns |

US-7653851-B2 | January 26, 2010 | Janusz Rajski, Jerzy Tyszer, Nagesh Tamarapalli | Phase shifter with reduced linear dependency |

EP-2144134-A1 | January 13, 2010 | Mentor Graphics Corporation | Method for synthesizing linear finite state machines |

US-2005097419-A1 | May 05, 2005 | Mentor Graphics Corporation | Method and apparatus for selectively compacting test reponses |

US-2007300110-A1 | December 27, 2007 | Janusz Rajski, Jerzy Tyszer, Nagesh Tamarapalli | Phase shifter with reduced linear dependency |

US-5446741-A | August 29, 1995 | International Business Machines Corporation | Fast memory power-on diagnostics using DMA |

US-2003131298-A1 | July 10, 2003 | Mentor Graphics Corporation | Test pattern compression for an integrated circuit test environment |

US-2004172431-A1 | September 02, 2004 | Mentor Graphics Corporation | Method for synthesizing linear finite state machines |

US-2003120988-A1 | June 26, 2003 | Janusz Rajski, Mark Kassab, Nilanjan Mukherjee, Jerzy Tyszer | Continuous application and decompression of test patterns to a circuit-under-test |

US-2004128599-A1 | July 01, 2004 | Mentor Graphics Corporation | Decompressor/PRPG for applying pseudo-random and deterministic test patterns |

US-7805651-B2 | September 28, 2010 | Mentor Graphics Corporation | Phase shifter with reduced linear dependency |

US-2007011530-A1 | January 11, 2007 | Janusz Rajski, Jerzy Tyszer, Mark Kassab, Nilanjan Mukherjee | Decompressor/PRPG for applying pseudo-random and deterministic test patterns |