US3514585A - Multichannel correlator system - Google Patents

Multichannel correlator system Download PDF

Info

Publication number
US3514585A
US3514585A US552994A US3514585DA US3514585A US 3514585 A US3514585 A US 3514585A US 552994 A US552994 A US 552994A US 3514585D A US3514585D A US 3514585DA US 3514585 A US3514585 A US 3514585A
Authority
US
United States
Prior art keywords
signal
time
sampling
input
samples
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Lifetime
Application number
US552994A
Other languages
English (en)
Inventor
Keith H Norsworthy
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Boeing Co
Original Assignee
Boeing Co
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Boeing Co filed Critical Boeing Co
Application granted granted Critical
Publication of US3514585A publication Critical patent/US3514585A/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06G—ANALOGUE COMPUTERS
    • G06G7/00—Devices in which the computing operation is performed by varying electric or magnetic quantities
    • G06G7/12—Arrangements for performing computing operations, e.g. operational amplifiers specially adapted therefor
    • G06G7/19—Arrangements for performing computing operations, e.g. operational amplifiers specially adapted therefor for forming integrals of products, e.g. Fourier integrals, Laplace integrals or correlation integrals; for analysis or synthesis of functions using orthogonal functions
    • G06G7/1928—Arrangements for performing computing operations, e.g. operational amplifiers specially adapted therefor for forming integrals of products, e.g. Fourier integrals, Laplace integrals or correlation integrals; for analysis or synthesis of functions using orthogonal functions for forming correlation integrals; for forming convolution integrals
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06G—ANALOGUE COMPUTERS
    • G06G7/00—Devices in which the computing operation is performed by varying electric or magnetic quantities
    • G06G7/12—Arrangements for performing computing operations, e.g. operational amplifiers specially adapted therefor
    • G06G7/19—Arrangements for performing computing operations, e.g. operational amplifiers specially adapted therefor for forming integrals of products, e.g. Fourier integrals, Laplace integrals or correlation integrals; for analysis or synthesis of functions using orthogonal functions
    • G06G7/1921—Arrangements for performing computing operations, e.g. operational amplifiers specially adapted therefor for forming integrals of products, e.g. Fourier integrals, Laplace integrals or correlation integrals; for analysis or synthesis of functions using orthogonal functions for forming Fourier integrals, harmonic analysis and synthesis
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06G—ANALOGUE COMPUTERS
    • G06G7/00—Devices in which the computing operation is performed by varying electric or magnetic quantities
    • G06G7/12—Arrangements for performing computing operations, e.g. operational amplifiers specially adapted therefor
    • G06G7/19—Arrangements for performing computing operations, e.g. operational amplifiers specially adapted therefor for forming integrals of products, e.g. Fourier integrals, Laplace integrals or correlation integrals; for analysis or synthesis of functions using orthogonal functions
    • G06G7/1928—Arrangements for performing computing operations, e.g. operational amplifiers specially adapted therefor for forming integrals of products, e.g. Fourier integrals, Laplace integrals or correlation integrals; for analysis or synthesis of functions using orthogonal functions for forming correlation integrals; for forming convolution integrals
    • G06G7/1935—Arrangements for performing computing operations, e.g. operational amplifiers specially adapted therefor for forming integrals of products, e.g. Fourier integrals, Laplace integrals or correlation integrals; for analysis or synthesis of functions using orthogonal functions for forming correlation integrals; for forming convolution integrals by converting at least one the input signals into a two level signal, e.g. polarity correlators
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06J—HYBRID COMPUTING ARRANGEMENTS
    • G06J1/00—Hybrid computing arrangements
    • G06J1/005—Hybrid computing arrangements for correlation; for convolution; for Z or Fourier Transform

Definitions

  • a multichannel time correlation computer which includes means for efficiently establishing an array of time delay differentials between signal representations to be multiplied in deriving correlation coefficients.
  • the disclosed combinations of circuit features enable the real time presentation of a large number of simultaneous, continuously produced correlation coetlicient signals which collectively represent the correlation function, so that it can ⁇ be observed as it changes with time.
  • the array of time delay differentials enables conservation of multiplier circuitry..Specic improvements are disclosed in arrangements for establishing the time delay differential arrays and for multiplying signal values in analog and digital form.
  • This invention relates to systems for computing correlation functions, and rnore particularly to a multichannel instrument for producing continuous representations of time correlation functions by simultaneously computing time correlation coefficients for a multiplicity of different delay times. It further relates to systems for cornputation of Fourier transformations and similar functions, thereby prviding the capability of obtaining power spectra and cross power spectra directly from a continuously presented correlation function.
  • the invention encompasses certain specific improvements in time correlation instruments and the like, including improved techniques and devices for sampling input signals for both high and low frequencies.
  • correlation studies produce information concerning the ⁇ relationship between two functions or variables, such as two electrical signals (cross-correlation), or the relationship of a function to itself when displaced or delayed -by known amounts (autocorrelation).
  • cross-correlation or the relationship of a function to itself when displaced or delayed -by known amounts
  • autocorrelation In electronics, for example, correlation techniques are .useful in determining what effect noise produced at a selected point in a system has upon the output of that system, i.e. determining the correlation between output noise and noise produced at an internal point.
  • Correlation studies, ⁇ are by no means restricted to electronics, however.
  • Correlation techniques are used to detect common characteristics of signals obscured -by presence of uncorrelated noise, ⁇ and the ⁇ like, which statistically averages out to zero or a constant when averaged over a sufficient period of time.
  • time correlation function of any two time-varying signals S1(t) and S2(t) is where 1- is a time delay parameter, i.e. a selected amount of time delay occurring after the time instant t.
  • the time correlation function is a function Vof time delay parameter 1- and corresponds to the long term average of the product of two time signals S1(t) and S2(t) multiplied together after a time delay -r has been applied to one.
  • the valve of this function for any particular value of r is termed the correlation coeicient.
  • Such a function is useful because it contains information on the amount of signal which is common to both S1(t) and S20). If signals S1(t) and S2(t) are completed unrelated (in the time domain) the correlation function will be zero for all values of T; but if they are related by common frequency components the time correlation function will be dependent on a composite of the correlation components. In such applications the function is referred to as a crosscorrelation function.
  • Instruments have heretofore been constructed for computing time correlation functions, operating on the principle of implementing the above formulae. That is, they have included means for delaying one input signal with respect to another, multiplying the signals together, and averaging the product over a period of time.
  • Such instruments are of two general classes: the continuous input type and the sampling type.
  • the means for establishing delay between continuous input signals usually comprises means for recording or storing one signal on a magnetic tape or drum and repeating the same back after a selected length of time of transit between recording and pickup heads.
  • Acoustical delay devices have also been used in such instruments. While this type of instrument has advantages for some applications, it involves high costs attending complicated recording and pickup devices and necessary compensations because of distortions due to frequency response characteristics of magnetic tape recorders or other delay instruments used, and for other reasons.
  • Correlator instruments of the sampling type generally implement the following statistical forms of the above integration formulae for the cross-correlation and autocorrelation functions, respectively:
  • One method for obtaining power spectra or crosspower spectra has been to compute the Fourier transform of a previously obtained correlation function by digital computing techniques, or the like, involving rather extensive and time-consuming treatment of information.
  • Other well-known means exist for obtaining power spectra and cross-power spectra independently of correlation functions, namely spectral analyzers and the like.
  • Another important object hereof is to provide in such an instrument the capability of computing a large number of correlation coeicients, with minimum apparatus requirements, thus permitting achievement of economical compact and easily maintained systems for generating continuous representations of correlation functions.
  • a further object is to provide an instrument capable of computing correlation functions for time delays ranging from negative to positive values.
  • Another object hereof is to provide a. time correlation computer in which, by appropriate selection of values of delay time 'r, the computations may be centered about any selected portion of a correlation curve.
  • a related object is to provide such an instrument, which permits varying the extent of the range over which a correlation function may -be computed.
  • a further object is to provide such an instrument which is adaptable to permit computing simultaneously therein the respective correlation functions of different pairs of input signals.
  • An additional highly important object is to provide a time correlation instrument with means permitting sampling a signal or signals at high sampling rates, with computation delay times 1- as long as may be desired, and capable of handling a range of input signal frequencies as broad as may be desired.
  • a related object hereof is to provide in one embodiment of the invention a multichannel instrument capable of handling much higher input signal frequencies, and in another embodiment much lower input signal frequencies, than was possible heretofore in a time correlation computer.
  • a further object of the present invention is to provide a multichannel instrument which makes use of digital components such as binary shift registers.
  • An additional object of the present invention is to provide an improved multichannel correlator making use of digital components and wherein the product rate is readily increased for those applications wherein A1- is fixed by application requirements.
  • a further object of the present invention is to provide improved multiplication circuits for receiving a pair of multi-bit binary signals and providing an analog output signal which is proportional to the product of the two binary input signals.
  • Another object of the present invention is to provide a multichannel correlator system using digital techniques which are readily implemented through the use of multichannel shift registers, magnetic core storage units, or a magnetic drum for producing Cartesian product patterns of delays.
  • Still another object of the invention is to provide improved electronic time correlation instrument circuit arrangements which are well suited to transistorization and miniaturization.
  • sampling type correlators primarily in terms of its application to sampling type correlators, but it should be understood that most of its features also apply to continuous input type correlators.
  • sampling type correlator and continuous input type correlator will be obvious to those skilled in the art.
  • a system intermittently samples each input signal a selected number of times, delays the set of samples of one signal by predetermined amounts with relation ⁇ to those of the other signal, and applies both the delayed and undelayed samples to multiplying means adapted to produce the Cartesian product thereof in terms of circuit values (Le. voltages or currents).
  • the Cartesian product of the two sets of signal samples as thus derived is herein defined as the total set of products of all samplesin one set multiplied individually by all samples of the other set.
  • the multiplying means thus produces a plurality of output products equal in number to the .product of the numbers of samples in the respective sets.
  • Additional features of the invention reside in certain disclosed means for sampling input signals and for delaying and ⁇ pairing for multiplication the samples thus obtained, such means being suited especially for signals comprising high frequency components.
  • This portion of the system in its preferred form employs two sets of sampling gates of improved construction, one set for each input signal.
  • the gates in one set are individually actuated to sample one input signal at successively delayed points of time during each operating cycle.
  • the individual gates in the other set are successively actuated to sample the other signal.
  • the individual samples from the second signal are respectively paired with and multiplied by the immediately preceding samples from the rst signal, thereby to produce a plurality of products, respectively corresponding to different delay times r, equal in number to the product of the numbers of gates in therespective sets.
  • the cyclically produced products are individually time-averaged to produce continuously a multiplicity of output currents or voltages representative of computed correlation coefficients and which collectively represent the correlation function.
  • sampling and delaying means ⁇ especially adapted for input signals having relatively low frequency components comprises means for sampling each input signal at least once during each system cycle, and a plurality of means for storing the individual samples over time periods spanning at least several cycles. Selection means pair these samples, one sample of each signal in each pair, in the different possible combinations for multiplication.
  • the members of each pair are of different ages, having been stored for different lengths of time, and the relative ages vary from pair to pair to establish a suitable" ⁇ array of time delays r as in the embodiment previously described. A maximum number of such pairs are simultaneously presented to the multiplying means, the number being determined by the number of different time delays r obtainable from the available number of stored samples.
  • the two input signals are first applied to analog-to-digital converters for periodic sampling and conversion to digital information.
  • the digitizing rate of the first converter for one signal is fixed by the frequency of a clock pulse generator, while the digitizing rate of the other is set at a frequency of f/N where N is equal to the number of words or stages in a shift register coupled with the output of the rst converter.
  • the output of the second converter is applied to an M stage shift register having a shift rate of f/N. On each Nth clock pulse the contents of the N stage shift register are transferred to an N stage temporary lstorage unit.
  • the contents of the temporary storage unit are then matched with the various stages of the M stage register and applied to a plurality of multiplying circuits disposed intermediate the temporary storage unit and the output averaging circuits.
  • the system includes improved multiplying circuit arrangements providing pulse width-pulse height modulation, binary weighting of a width modulated pulse, and digital combination of paired bits words from two with the resulting signals being added by binary weighted summing resistors.
  • the increment of time delays between correlation coefficients is fixed and yet it would be desirable to increase the product rate to a value higher than that normally obtainable in a system producing an MXN delay pattern, even with N equal to unity.
  • a system which makes use of a single stage N shift register and a kM stage M shift register. Every kth stage of the M shift register has its output circuits paired with adjacent or intermittent output circuits of the N shift register.
  • the clock pulse frequency which controls the sample rate and shift rate of the M and N registers is then increased by a factor of k times the frequency f which would be used in an MXN delay pattern correlator operating on the same fixed increment of time delays between correlation coeiiicient.
  • the results is that the product rate is increased with A1- being held fixed in accordance with application requirements.
  • the invention also includes means for varying the timing control frequency so that increments between successive time delays (T) may be uniformly and simultaneously increased or decreased.
  • the system includes variable delay means operative upon the sampling and delay means to permit changing the relative time delays between the two sets of samples taken from signals S1(t) and 82(1), respectively.
  • FIG. 1 is a block diagram of the improved system for computing correlation functions
  • FIG. 2 is a schematic diagram of a preferred sampling gate
  • FIGS. 3a, 3by and 3c constitute a timing diagram illustrating the operating theory of the system.
  • FIG. 4 is a schematic diagram of sampling means for low frequency signals.
  • FIG. 5 is a block diagram of a digital embodiment of the present invention.
  • FIGS. 6, 7, and 8 are improved multiplying arrangements having particular advantage in the present systems and,
  • FIG. 9 is a block diagram of another form of digital system for providing increased product rates.
  • the system illustrated in FIG. 1 includes a sampling means 10 having inputs 10a and 10b to which the respective signals S1(t) and S2(t) are applied. In the case of cross correlation these are different signals, whereas they are the same signal S1(t) in the case of auto-correlation.
  • Sampling means 10 controlled by a timer 20 samples each input signal at timed intervals. The timing is such that a succession of samples of signal S2(t) are taken following a sucession of samples of signal S1(t). This alternate sampling of the two signals is repeated cyclically under control of timer 20.
  • Multiplying means 30 controlled by timer 20 has a series of inputs r, s, t, u, etc. to which the respective samples are fed (as voltages or currents).
  • the multiplier 30 includes pulse-width, pulse-height modulation means to produce at its outputs 31 a plurality of separate signals each representative of a product of two samples, one of each input signal. These output signals correspond to different values of computation time delay. During recurrent cycles of operation these output product signals (pulses) vary in magnitude with variation in magnitudes of samples of the input signals. Conversion of these product signals into continuous voltages or currents is accomplished by individual averaging circuits 41.
  • each averaging circuit 41 is a continuous analog representation of one correlation coefficient of input signals S1(t) and S2(t) for a particular delay time -r.
  • the collective averaging circuit outputs 42 carry a continuous analog representation of the complete correlation function.
  • the utilization means 50 is responsive to all the averaging circuit outputs 42 and can be display or recording apparatus for the correlation function, or when used for computing power spectra the utilization means 50 comprises means to continuously transform the computed correlation function information into power spectrum information.
  • control pulse generator 25 which may constitute any suitable oscillator of high stability.
  • a crystal-controlled blocking oscillator operating at one megacycle per second is suitable.
  • the basic control pulses from generator 25 are fed into a timing control selector 26, preferably comprising a binary-base frequency divider.
  • Timing control intervals established by selector 26 may be adjusted or varied by factor increments of two, such as from one microsecond to two microseconds, four microseconds 1,024 microseconds.
  • the resultant control pulses of selectively reduced frequency are then fed into a master timing control circuit 32 whereby they are distributed to various parts of the correlator system.
  • Timing control circuit 32 preferably comprising binary-base digital components, includes a series of switching circuits 34 (individually numbered l through 6) for controlling the respective sampling gates 11, 12, 21, 22,- 23 and 24 through connections 35a, 35b 3 51.
  • the switching circuits 34 preferably comprises D-C controlled bistable (scale-of-two) switching circuits or flip-flops of any suitable or known type. When triggered recurringly each produces two square waves of opposite phasing.
  • the function of the variable delay means 38 in the timing control circuit will be explained hereinafter.
  • signal S1(t) is continuously applied to the inputs of gates 11 and 12 and S20) similarly to the inputs of gates 21, 22, 23 and 24.
  • the timing of gate actuation for momentary signal sampling normally takes place in the following order: Gate 11 followed by gates 12 and 21 actuated simultaneously, thence followed successively by gates 22, 23 and 24. The total cycle is repeated recurringly by the timer 20.
  • these gates preferably are of the sample and hold type, each retaining its sample voltage or current for subsequent use by multiplier 30.
  • Such a gate circuit, illustrated in FIG. 2 consists mainly of a diode bridge switched alternately to on and o states by voltage waves from timer 20 applied through control transistors T1 and T2.
  • the collectors of these transistors are connected to DC voltage sources (not shown) typically of +12 and -12 volts and their emitters are connected directly to bridge terminals 104 and 106, respectively. Their emitters are also connected through resistances R1 and R2 to the DC voltages sources of -12 and +12 volts, respectively.
  • the base voltage of each control transistor controls its emitter voltage, and therefore the bias voltage applied to the bridge.
  • the timing control voltage applied to the bases of the control transistors T1 and T2 are -ll volts and +11 volts, respectively, the bridge is switched on, whereas a reversal of these voltages switches it off.
  • a double transistor emitter-follower circuit comprising transistors T3 and T4 accepts the input signal S(t) and ensures low drive impedance to the bridge input terminal 108.
  • the gate circuit is capable of following both positive and negative signal excursions. Moreover, because the time period from on to oi condition can be made extremely short, and because the capacitor voltage does not have change in this time period, very high-frequency signals S) can be sampled accurately to represent the Value of the input signal at any instant.
  • FIGS. 3a, 3b and 3c illustrates two complete cycles of operation. In each part of ⁇ the diagram both the sampling phase and the multiplying phase are shown.
  • wave-forms 54 depict the on and folf conditions of the sampling gates. ⁇
  • the consecutively numbered time markers 52 represent the output pulses from timing control selector 2'6(FIG.V 1) applied Vto the timing control circuit 32.
  • the sampling phase begins when the first gate 11 is switched on Capacitor C rapidly charges to a voltage corresponding to instantaneous input signal voltage and thereafter follows the latter for the duration of time the gate is on. In the ⁇ example this follow time is three time units 52.
  • the gate 11 is switched of ⁇ r ⁇ , the condenser voltage of that instant represents the desired sample of the signal, and is held duing the remainder of the cycle, extending through the multiplying phase. It is held until the gate is again switched on for resampling the signal. All the sampling gates operate in this way beginning at different times during the cycle.
  • gate 11 takes its sample of input signals S10) at point of time l
  • gate 12 takes its sample ⁇ at point 2.
  • Gates 21, 22, 23 and 24 sample signal S20) sequentially beginning with operation of ⁇ gate 21. This may begin, as illustrated, at the same ⁇ time (point 2) as sampling by gate 12.
  • the samples taken ⁇ by gates 11, 12, 21, 22, 23 and 24 are respectively designated a, b, c, d, e, and f. It will be observed that the time delays which exist between either of the samples a or b and different individual samples c, d, e and f can be arranged in a linear succession, as follows:
  • the multiplying portion of the cycle begins.
  • the ⁇ function of multiplying is to produce a number of output voltages each of which is the product of two different signal samples (from S10) and 820)) which are cross-correlated for a particular time delay. In the example, eight such output voltages are provided, corresponding to the eight time delays.
  • the multiplying means comprises a triggered saw-tooth generator 36 responsively connected through lead j' to the last stage (6) of time circuit 32,
  • the saw-tooth output waveform 55 appearing in FIGS. 3a and 3b passes from a negative voltage -vs to a positive peak voltage ml-vS, this signal excursion exceeding negative-to-positive excursions of input signal voltage S10). It is applied to one input of each of two amplitude comparison trigger circuits 44a and ⁇ 4417 which include input addition stages 44a' and 44b, respectively. Circuit 44a has a second input (r) to which the sample voltage from gate 11 is applied, and circuit 44b has a second input(s) to 'which the sample voltage from gate 12 is applied.
  • Each of these trigger circuits adds its two inputs in the addition stage, compares this sum to a reference voltage which in this case is zero, and produces fannoutput trigger pulse when the inputvolt-V age sum becomes equal to that reference voltage.
  • sample a gate 11 output
  • sample b gate 12 output
  • Trigger circuit 44b thus produces an output pulse when the saw-tooth output 55 reaches a value equal to the negative of the sample b voltage (-vb)
  • trigger circuit 44a produces an output pulse later in time when the sawtooth wave-form reaches a value equal to the negative of the sample a voltage va).
  • Two variable-length pulse generators 45a and 45h are respectively connected to be initiated by these output trigger pulses from trigger circuits 44a and 44h. These generators produce pulses 56a and 56h (FIG. 3b), respectively, of constant heights and of lengths determined by times at which they are triggered by the respective comparison circuits. Both generator circuits 45a and 45h are also connected to the saw-tooth generator 36 to be triggered olf by termination of the sawtooth.
  • variable-length pulses from generators 45a and 45b therefore begin at varying intervals after the instant the last sample f is taken depending upon the magnitudes of the respective samples from gates 11 and 12, and terminate at a constant interval after sample f is taken, because the saw-tooth wave is of xed duration.
  • variable length pulse generator outputs 56a and 56h variations are seen from one cycle to the next, successive pulses varying in length in linear proportion to successive sample values of signal S10) according to the following formula:
  • Pulse length L4-A181012) (7) where L and A1 are constants.
  • the saw-tooth, coincidence and variable-length pulse generative circuits may be of any suitable or well-known form so as to require no detailed description herein.
  • the multiplying means 30 further includes amplitude modulators 48, corresponding in number to the number of available time delays r (see Table I), each having two input terminals, one connected to a variable-length pulse generator 45a or 4511 and the other to one of the S20) gate circuit output terminals t, u, v or w.
  • the function of each amplitude modulator circuit 48 is to modulate the constant-magnitude pulse received from one of the variable-length pulse generators (45a or 4517) in linear proportion to the magnitude of one of the samples of signal S20) according to the following formula:
  • Equation 9 contains a first term which corresponds to the products of the respective pairs of original signal samples, and a second term proportional to sample values of S20). When averaged over many cycles this second term reduces to zero or a constant.
  • FIG. 3c the set of output pulses 57 (respectively designated b-c, a-c, b-d, etc., according to the signal sample products represented) from the multiplying modulators ⁇ 48 is depicted during two operating cycles of the system. It will be noted that some of these output pulses are positive and some negative, that each may change from one polarity to another and that their magnitudeduration characteristics are subject to variation, all in accordance with the multiplied values of the paired signal samples. It will also be seen that because of the assigned time delays in the order of sampling signals S1(t) and S2(t) each of the pulses 57 corresponds to a product S1(tk)S2(tk-i-T) for a different time delay fr.
  • each series ⁇ of combined variable-length, variable-magnitude, variable-polarity modulator output pulses 57 appearing at outputs 31 is next resolved into an output voltage, current or other measurement quantity such as a delay interval, a digital representation, etc.
  • the product information carried by pulses 57 is converted directly by averaging means 40 into a continuous, variable voltage or current in a process which suitably performs the time integration operation in the crosscorrelation formula (Equation 1).
  • the individual outputs 31 are applied to separate averaging circuits 41, each preferably consisting of a simple RC filter having a time constant many times longer than the length of a cycle of operation of the system.
  • a cross-correlation coeflicient (rn) is produced as each series of pulses 57 for a different time delay -rn is averaged individually during system operation.
  • variations in the average of a series of pulses 57 represent variations in the individual cross-correlation coeflicient thus computed.
  • the correlation function (1) itself comprising the collective correlation coefficient outputs, may change in shape and magnitude with time because of time variations in individual correlation coeicients.
  • the averaging circuits 41 perform the integration and division, collectively producing at outputs 42 the correlation function q (-r), made up of eight uniformly distributed correlation coefficients (vn) produced in analog form continuously during system operation. These may be displayed directly on an oscilloscope or recording apparatus (not shown), or may be otherwise employed in utilization means 50 as later described.
  • the timing control selector 26 permits varying the length of the basic opera-V tional time intervals 52, and thereby the intervals between signal samples, and in fact the length of the entire system operating cycle. Thus the spread between successive delays vyby which the cross-correlation function is represented may be varied at will.
  • the timing control apparatus 20 also includes means permitting production of the correlation function of the signals in regions displaced from zero.
  • a variable delay means 38 of any suitable or well-known form is provided in the timer circuit interposed between switching circuits 2 and 3. Depending upon the setting of this delay means 38, the samplings of signal S2(t) may be delayed by any desired amount from the samplings of signal S2(t). This is illustrated in FIG.
  • any change in the setting of timeY delay means 38 produces an appropriate shift in the starting time for the saw-tooth pulse, as shown by dotted line 55.
  • the multiplying phase of the cycle is always made to follow the sampling phase.
  • the saw-tooth generator 36 may have an adjustable saw-tooth slope control in order to permit varying the length of the multiplying phase of the cycle to compensate for changes in gain due to control frequency or time delay changes.
  • this variable delay capability permits real time display (or other utilization) of the correlation function in any region thereon including or removed from zero delay
  • the frequency selection capability permits expanding or contracting that region in order to vary the degree of definition of the cross-correlation function with the available number of delay values r. Further advantages of the variable delay capability will be discussed hereinafter.
  • the sampling means functionally illustrated in FIG. 4 for low-frequency signals comprises two stepping switches 58 and 59, ,having ganged sets of wipers W and each having nine contacts.
  • Contacts X1 to X9 of stepping switch 58 are connected to individual signal sample storage condenser c', and the contacts Y1 to YQ of switch 59 are individually connected to storage condenser C.
  • Switch 58 has a condenser charging Wiper W1 and condenser voltage sampling wipers Wa and Wb ganged there- 'with in respective positions angularly offset from wiper W1 by oneand two-contact intervals.
  • S1(t) is applied to wiper W1 through appropriate connections to input terminal 10a.
  • switch 59 has a condenser charging wiper W2, electrically connected to receive input signal S2(t) applied at input 10b, and four condenser voltage sampling wipers Wc, Wd, We and Wf ganged to wiper W2 in respective positions angularly offset therefrom at two, four, sixand eight-contact intervals.
  • the wipers of both switches are rotated clock-wise in stepwise fashion by stepping mechanism 61 actfuated periodically by a timer 60. Intervals between successive samplings and increments (T) between correlation coefficients are determined by the stepping periodof the switch, which may be varied by timer 60.
  • a typical range of variation is from 0.025 second to 10.0 seconds.
  • Storage condensers C' and C function in the same manner as output condensers C (FIG. 2) in the diode sampling gates in the first embodiment.
  • input signals S1(t) and S2(t) are here applied to a pair of condensers ⁇ C and C, that is, when charging wipers W1 and W2 are positioned on a pair of X-Y contacts during intervals between switching movements, the condenser voltages follow the input signal voltages until the wipers are stepped to the next set of X-Y contacts. When this occurs the two charged condensers remain charged to voltages corresponding to the respective input signal values at the instant of switching.
  • the digitized sample is available at the output of each gate and converter unit immediately following sampling of the signal.
  • the signal S1 is sampled at a rate corresponding to the frequency f of the clock pulse generator.
  • a count-down or dividing circuit 116 operates to divide the frequency f by the number of stages in the shift register 112, N.
  • the output of the divider 116 is applied to the sample gate of unit 111 so that the second input signal S2 is sampled at a rate equal to ⁇ f/N and the contents of register 213 are shifted at a rate of f/N.
  • the divider 116 also controls a transfer circuit 117 which is disposed between the shift register 112 and a temporary storage unit 118 having N stages therein with each stage capable of storing a word having the same number of bits as;the number of bits in each word in the shift register 112. On each Nth clock pulse the contents of the shift register 112 are transferred in parallel to the storage unit 118. ⁇ The output circuits 119 and 120 of the storage unit 118;are each paired with each of the output circuits 12111122, and 123 of the shift register 113. Each pair of circuits will be seen to be connected to one of the multiplying circuits 13D-135. The multiplying circuits are indicated as being of the analog output type with each having its output circuit coupled with an associated averaging circuit 136.
  • Each averaging circuit 136 has its output coupled with the utilization means 137.
  • the circuits 130-135 in one embodiment serve to provide rectangular output pulses 138 ⁇ which are pulse width-pulse lheight modulated in accordance with the binary values of the two input signals applied to the multiplier circuit.
  • the ⁇ signal S1 is digitized at the frequency of the clock pulse generator with the digitized information being applied to the shift register 112 and shifted in parallel therein at the frequency of the clock pulse generator.
  • the contents of the shift register 112 are transferred to the storage unit 118.
  • the signal S2 is sampled and the contents of the shift register 113 advanced to the right by one stage.
  • the multiplying means 130-135 provide M XN products where M and N are the number of words in the shift registers 113 and 112.
  • the delay value related to each multiplier output is fixed, and in the embodiment of FIG.
  • the individual ⁇ delays take on the values (LAT, ZAT, 3A ⁇ r-(MN-1)A1, Where is the reciprocal of the clock frequency f.
  • the sample product rate an any multiplier output channel is equal to the clock frequency divided by N, or f/N.
  • the sample product rate in the system of FIG. 5 is greater than the reciprocal of the maximum delay value (MN-DA1- as long as M is greater than unity.
  • each word in register 112 (and storage unit 118) is cornposed of six binary bits and that each word in register 113 is composed of four binary bits.
  • the sixl blnary signals identified as A0--A5 of the output circuit .119 are applied to a parallel comparator 140.
  • the comparator also receives six level binary information from the binary counter 141 which is repeatedly cycled from zero through (ZL-1) by pulses received from the pulse generator 142.
  • the pulse generator 142 operates at a frequency nf which is sufficiently high to permit complete cycling of the counter 141 between each shift of register 113.
  • nf a signal is provided on stop circuit 143 to stop a variable length pulse generator 144 which is started in operation by a start signal on start circuit 141A when the first pulse is applied to the counter 141.
  • the pulse generator 141 is operated in a manner such that a pulse is applied therefrom to the pulse height modulator 145 with the length of such pulse being proportional to the binary value of the A stage of the storage unit 118.
  • the four bits from the C stage of the shift register 113 are applied to a conventional digital-to-analog converter 146 having its output circuit 147 coupled for control of the pulse height modulator 145.
  • the result is a rectangular pulse 138 having a width corresponding to the binary value of the A stage of the storage unit 118 and an amplitude proportional to the value of the C stage of the shift register 113.
  • the multiplying circuit receives binary input signals and provides an output pulse having an area proportional to the product AXC. Accordingly, the multiplier circuit 130 is referred to as an analog output multiplier.
  • FIG. 7 it will be seen that the various components 140-144 of the system in FIG. 6 are utilized.
  • the variable length pulse from the generator 144 is applied simultaneously to each of the four binary AND gates 150-153.
  • the second input for each of the gates 150-153 is respectively one of the binary value C0-C3 corresponding to one of the four bits of the C word in the C stage of shift register 113.
  • Output signals from the gates 150-153 are combined by means of binary weighted resistors 154-1S7 so that an output signal proportional to the product AXC is provided to the averaging cir- 150-153 applies zero voltage (or a reference voltage V0) to the resistor associated therewith when the gate is closed When the gate is open it passes a voltage which is a function of the state of the binary bit preceding it. Thus if the input C0 to gate 150 is O the gate 150 passes zero voltage (or a reference voltage V0) whereas if the input C0 is l the gate passes a different reference voltage V1.
  • the binary input signals Co-Ca therefore serve to open the associated gate when the input is a l -so that the reference voltage V1 is applied to the weighting resistor for a length of time determined by the length of the pulse from the variable pulse generator 144.
  • Correlation coefiicient outputs could be derived by placing an averaging capacitor between the resistor summation terminals and ground or by feeding the summation terminal to an integrator.
  • FIG. 8 a circuit arrangement is illustrated for providing the required digital multiplication in a simplified manner by first combining the individual -words digitally and then producing an analog output by means of a digital-to-analog conversion.
  • the individual -words digitally and then producing an analog output by means of a digital-to-analog conversion.
  • Voltage sampling wipers Wa, Wb, Wc, Wd, We and Wf are electrically connected to apply the stored samples to respective inputs r, s, t, u, v and w of multiplier 30 to be by pairing samples stored diierent lengths of time over previous cycles, greatly increases the possible number of product pairs available for presentation to each output channel per averaging circuit time constant, over that previously possible in analyzing low frequency signals.
  • this procedure is Y*wiper pairs formed -by virtue of wiper connections to the implemented by a switch 82 selectively operable to intermultiplier inputs land further by virtue of appropriate pose a single inverter 80 in the signal (S2(t)) inputs of the cross-multiplication connections in the multiplier 30 (FIG. two sampling gates 21 and 22. So that both inverted and 1) itself.
  • Each succeeding row illustrates for succeeding uninverted values of the correlation coeiiicients 1p (Tn) cycles corresponding pairing of respective X-Y contacts will be available in all output channels of the correlator, by virtue of the above-mentioned connections and the re- Switches 84 in the switching circuit leads for gates 23 and spective wiper positions on the contacts.
  • Each column in r 24 are ganged with inverter switch 82 and Iare operable the table illustrates successive pairings of different con- 30 therewith to connect these two leads respectively to those tacts (and corresponding sampling condensers) to profor sampling gates 21 and 22.
  • this relation coeflicient for a particular time delay 'r. Comparichanges the relative timing of the last two sampling gates son of Table III with FIG. 4 shows that a complete set of 40 23 and 24 in the series to sample at times e and f" as six samples applied to the multiplier during each cycle illustrated by dot-dash lines.
  • circuitry will be aiected by this input inversion technique AS llldleated PreVlQnSlyg var la'flons 1n leng'fbs 0f u me 50 for producing inverted correlation coeicients, so that it delays me made POSSlble m t1115 embodiment by Vafylng is normally preferred to invert the signal input involving mei Period betnfeen (llePPlng Inovfnentts wthes s; the fewer number of input (sampling) circuits. 31.1. 59' Morflme e ays T may e es a 1S e y pr. In the embodiment of the invention illustrated in FIG.
  • s S1 and S2 are respectively applied to ers, or their equivalents, a small number having been the sam 1in ate and anale t, di it l o t .t chosen herein for simplicity of illustration only-
  • the con- 110 andpu lhich riodic 0S i all c Ivl: 6.? nils cept disclosed is of course not limited to use of stepping in ut Si Hals Th Se 1t. ady .t al p an lgldlze t e switches provided with storage condensers.
  • the stability (accuracy) of the measurement is re- Shift re gister 112 is shown as having two stages, or word lated to the number of samples or sample pairs received b4 each out ut Channel durin the integrating time of the 65 storage sections, each word being of six binary bits.
  • the y p g shift register 113 is shown as having three stages of four output averaging filters. When high-frequency signals are bits each The number of bits er Word in th t analyzed, increments (Af) between time delays can be ters have'be h b p.
  • c ook Pulse generator 114 @Peratmg a frequency has its signal output circuit 115, connected directly with the first gate analog-to-digital converter unit 110.
  • the unit and 111 preferably each include a sampling gate of the type shown in FIG.
  • A3(23)+A4(24) +A5(25) C C(2)+C1(21)+C2(22)+C3(23) where A0, A1 A5 and C0 C3 take on values of O or 1, the product A XC can be broken into a number of components and a Weight assigned to each as follows:
  • each product such as A5C2
  • A5C2 is unity if both A5 and C2 are "1 and zero if not.
  • the product term may be sutliciently ⁇ defined if only the -iirst four components in the above chart are accounted for.
  • These product terms are readily obtained through the use of simple diode AND gates and weighting resistors.
  • a number of diode AND gates 160 are arranged for receipt of the Various binary bits from the two words to be multiplied with the outputs of each individual AND gate being applied through the indicated resistor networks to a common summing junction 161 for application to an averaging circuit 136.
  • the ⁇ values corresponding to 28 have the greatest weight by ⁇ being applied to the summing junction 161 through a resistor 162 having a resistance R, the 2z components are ⁇ summed by a pair of resistors 163 each having a value of 2R, the 26 values are summed through the resistors 164each having a resistance of 4R, and the 25 values are summed through the resistors 165 each having a value of SR.
  • the sample product rate is given by the reciprocal of the quantity (nAr) Where N is the number of Words in the shift register 112 (and also the number of words in the storage unit 118).
  • Ar is fixed by the requirements of the system and accordingly it will be seen that the product rate is maximized when N equals unity. In that case the system then becomes one having a l XM delay pattern capability rather than a Cartesian N XM delay pattern.
  • A1- must remain fixed it is often desirable to obtain a product rate even greater than that obtained by having N equal to unity in a system such as that of FIG. 5.
  • the system shown in FIG. 9 makes possible the achievement of such higher product rates.
  • the input signals S1 and S2 are respectively applied to the sampling gate and analog-todigital converter units 210 and 211 which are similar to the units and 111 in FIG. 5. However in the system of FIG. 9 these units serve to sample both of the input signals at the rate of the occurrence of pulses from the clock pulse generator 214.
  • the digitized signals from the unit 210 are applied to a single word storage unit 212 (which can also be referred to as a single stage shift register) and the digitized signals from the converter 211 are applied to the stage shift register 213.
  • the clock pulse generator 214 is shown as operating at a frequency of K times the frequency f of FIG.
  • each Kth stage of the register 213 has an output circuit.
  • the frequency of the clock pulse generator is 3 times that of the frequency of the clock pulse generator in the system of FIG. 5 and accordingly every third stage of the shift register 213 has its signal output circuit paired with the output circuit 215 of the single stage register 212.
  • the input signals are sampled at the rate Kf of clock pulse generator 214, and that the product rate is also Kf. Since only every Kth stage of shift register 213 provides output signals AT remains at the desired constant value. It will therefore be seen that the system of FIG. 9 provides a convenient method of obtaining an increased product rate even though A1- must remain fixed for a given application. Even though the number of stages in the register 213 must be increased by a factor of K it should be noted that since only every Kth stage must be capable of providing output signals the intermediate stages can be low cost magnetic shift register stages.
  • addressable digital memories such as magnetic cores or magnetic drums could be used.
  • a core memory could readily be used with the digital words corresponding to signal samples being entered into the core memory by sequentially in creasing the write address registers.
  • Several read addresses would be selected in parallel with their addresses sequencing around at an appropriate iixed number of increments behind the write addresses.
  • the words read out of the memory would enter registers such as flip-hops connected to the multiplying means.
  • a system for computing a correlation function of electrical signals appearing at first and second circuit conductors comprising: signal input means connected to said conductors and providing first and second pluralities of time-related time successions of signal values at said first and second conductors, multiplication circuit means connected to said input means and recurringly cross multiplying each of the plurality of time-related time successions of signal values at the 'first conductor by each of the plurality of time-related time successions of signal values at the second conductor, the latter time successions of values bearing a predetermined time relationship to the former, to produce recurring sets of circuit responses representing recurringly the Cartesian product of the cross-multiplied signal values, and signal averaging means coupled with said multiplication means for separately time averaging the respective recurring circuit responses substantially continuously to represent respective correlation coefficients collectively defining said correlation function.
  • Correlation computer apparatus comprising rst and second signal inputs, first and second groups of signal sampling devices to which the inputs are respectively applied, timing means operable through recurring cycles and connected to the sampling devices to operate those in the iirst group at successively delayed times and those in the second group at successively delayed times in timed relation with the sampling times of the first group, each such device having a storage element therein which tracks the successive sample values of signal and which, between sampling times, retains the signal sample therein, means forming a plurality of output channels each with two inputs operatively associated respectively with sampling device storage elements of the first and second groups, whereby said output channels receive respectively different combinations of two signal samples including one from each sampling device group, means in each output channel for multiplying the recurringly received signal samples, and separate means coupled with each output channel for averaging recurring products therein to derive respective correlation coeiicient outputs.
  • the output channel inputs operatively associated with the sampling device storage elements in the first group comprise means operated cyclically by the timing means to generate in each such latter input recurring pulses each having an amplitude-duration pulse envelope product proportional to the currently stored signal sample in the associated sampling device, and wherein the inputs of all output channels comprise pulse modulators receiving such pulses and operable to vary the amplitude-duration products thereof in proportion to the stored signal samples in signal sampling devices of the second group associated with such inputs, thereby to multiply together the recurringly received signal samples operatively associated with each output channel, and means in each output channel responsive to the modulated pulses to produce an output response proportional to the time average of varying amplitude-duration recurrent products.
  • the first-mentioned pulse generator comprises in each channel a triggerable variable-length pulse generator having a starting trigger input and a stopping trigger input, means for applying the stopping trigger simultaneously to corresponding inputs of such pulse generators in response to operation of the timing means, a sawtooth wave generator cyclically triggered by said timing means, amplitude comparison circuits having outputs connected to the respective starting trigger inputs and each having a irst input connected to the sawtooth wave generator and a second input connected to the respectively associated sampling device storage element, said comparison circuits further including means for adding said first and second inputs, means for comparing the sum of said inputs to a reference circuit value, and means operable to produce an output trigger pulse in response to amplitude coincidence between the instantaneous input sum and said reference value.
  • a system for computing the correlation function of a pair of input signals comprising a first input means to which one signal is applied, a second input means to which the other signal is applied, each of said input means having a plurality of outputs and being operable to produce in said outputs representations of the applied input signal bearing a predetermined set of successively delayed relationships to one another, the signal representations from one input means also bearing predetermined time relationships with the signal representations from the other input means; multiplying circuit means having a plurality of inputs; circuit means operatively associated with said input means and with said multiplying circuit means and operable to apply said signal representations to said multiplying circuit means in pairs of one signal representation from each of said input means with the signal representations in each pair bearing to each other a time-delay relationship differing from time-delay relationships of signal representatons comprising all other pairs, said multiplying means having a plurality of output channels each responsive to a different one of said pairs of signal representations and being operable to produce in said output channels respective output responses each representative of the product
  • each delay means comprises a set of sampling devices and cyclically operable timing control means actuating said devices at successively delayed times to derive successive sample values of the applied input signal as said output signal representations.
  • timing control means is operable to actuate the sampling devices of each set in successive Order during each cycle and is connected to said multiplying means to operate the same cyclically for multiplying said pairs of signal sample representations derived during each such cycle.
  • timing control means is operable to actuate the sampling devices for the first signal at equal intervals and those sampling the second signal at intervals of a length equal to the product of the interval between samplings of the rst signal multiplied by the number of samples taken of said first signal in each cycle, thereby to establish a uniform array of equal time-delay differentials between pairs of signal sample representations.
  • At least one such set of sampling devices comprises a set of storage elements and storing means sequentially cooperable therewith to derive a single signal sample during each cycle in successive storage elements
  • said multiplying means includes a plurality of separate multiplier units each adapted to receive a different pair of signals from said delay means, and sample selection means cyclically operable in synchronous relationship with said storing means to apply to said multiplying means in each cycle a plurality of stored samples which bear the same storage time relationship to each other as the plurality of samples applied in each preceding cycle and wherein each said stored sample is applied to a different multiplier unit.
  • said storing means includes a series of terminals electrically connected respectively to the individual storage elements, first switching means connected to receive said applied input signal and operable to engage said terminals sequentially, second and third switching means comprising said sample selection means and being also operable to engage said terminals sequentially, said second and third switching means being respectively connected to first and second outputs to which said multiplying means is responsive, and actuating means connected to said iirst switching means to engage and disengage the same with successive terminals in step-by-step manner whereby to derive successive signal samples in the respective storage elements, said actuating means being similarly connected to the second and third switching means to engage and disengage the second and third switching means with the respective terminals a predetermined whole number of steps after engagement therewith of the first switching means, and to engage and disengage the third switching means with the respective terminals a predetermined whole number of steps after the second switching means, thereby to apply the stored samples successively to each of said iirst and second outputs, said second delays means comprising
  • a first switch means including a series of terminals connected respectively to the individual condensers and including first, second and third switching elements oper- -able to engage said terminals sequentially, said first switching element being connected to the first signal input, first and second outputs to which the second and third switching elements are respectively connected, actuating means connected to said first switching element to operate ⁇ the same in step-by-step manner through recurring cycles, thereby to store in the respective condensers successive samples of voltage in the rst signal input during each cycle, said actuating means being similarly operatively connected to said second and third switching elements, with the second switching element engaging the respective terminals a predetermined whole number of steps after the first switching element and with the third switching element engaging such terminals a predetermined whole number of steps after the second element, thereby to apply the stored samples successively to ⁇ each of said first and second outputs, a second series of ⁇ storage condensers
  • means for establishing time-delay differentials between signal values for multiplication comprising in combination: a first series of storage condensers, a first switch means including a series of terminals connected respectively to the individu-al condensers and in cluding a plurality of switching elements 1, 2, 3 p, operable to engage said terminals sequentially, the first switching element being connected to the first signal input, a first plurality of outputs to which switching elements 2, 3 p are respectively connected, actuating means connected to switching element 1 to operate the same in ⁇ step-by-step manner through recurring cycles, thereby to store in the respective condensers successive samples of voltage in the first signal input during each cycle, said actuating means being similarly operatively connected to said switching elements 2, 3 p, with switching element 2 ⁇ engaging the respective terminals a predetermined whole number of steps after switching element 1 and with each succeeding switching element engaging such terminals a predetermined whole number of steps after its predecessor, thereby to apply the stored samples
  • a system for computing a correlation function of a pair of input signals comprising separate means for deriving a set of samples of each input signal, with the individual samples of one set occurring in predetermined time-spaced relationship with individual samples of the other set, multiplying means having a plurality of outputs and operatively associated with said separate means and operable to multiply each individual sample in one set by each individual sa-mple in the other set, thereby to produce a plurality of output product circuit responses in the respective outputs, means in such system for cyclically repeating such sampling and multiplying, whereby to produce such output product responses recurringly, and separate means operatively associated with the respective multiplying means outputs for continuously averaging the responses therein whereby to derive separate correlation coefficient outputs which collectively represent the correlation function of said input signals.
  • samplederiving means comprises first and second sets of sampling devices to which the input signals are respectively applied
  • means for cyclically repeating sampling and multiplying comprises timing control means operable cyclically to actuate said devices at successively delayed times during each cycle whereby to derive successive samples of the applied input signals and being further operable to actuate said multiplying means cyclically for multiplying the samples derived during each such cycle.
  • sampling devices comprise gate circuits each having a sample storage element therein operable to retain a signal sample value between successive operations thereof, and wherein said multiplying means outputs are responsive to respectively different combinations of two of such currently retained signal samplesincluding one of each input signal.
  • said multiplying means comprises a plurality of input circuits each responsive to one of the gate circuits in the first set to generate recurring pulses each having an amplitude-duration pulse envelope product proportional to the sample currently stored in such circuit, and a plurality of output circuits responsive to respectively different combinations of one of said input circuits and one of the gate circuits in the second set, each output circuit being operable to modulate the amplitude-duration product pulse in the associated input circuit in proportion to the sample currently stored in the associated gate circuit.
  • At least one of said separate signal sample deriving means comprises a plurality of storage elements and storing means cooperable cyclically therewith to store the successively derived signal samples in said elements, and sample selection means cyclically operable in synchronous relationship with said storing means to apply to said multiplying means in each cycle a plurality of stored samples which bear the same time relationship to each other as the plurality of samples applied in each preceding cycle.
  • a system for computing a correlation function of a pair of input signals comprising a separate means for successively sampling each input signal, with the individual samples of one input signal derived in time-spaced relationship with individual samples of the other input signal, circuit means forming pairs of such samples with each pair having one sample from each input signal and members of the respective pairs having different time relationships whereby to correspond to different predetermined amounts of time delay, said circuits means including means operable to multiply together substantially simultaneously the members of such sample pairs, thereby to produce respective output product responses corresponding to said predetermined amounts of time delay, means in such system for cyclically repeating such sampling, pairing and multiplying to produce respective series of such output product responses, and means operative through recurring cycles to average each such series individually.
  • a system for computing the correlation function of first and second input signals comprising in combination: first signal input means to Which the first input signal is applied having a plurality N of outputs and operable to produce in said outputs N simultaneous representations of the applied signal having equal time-delay differentials AT therebetween; second signal input means to which the second signal is applied having a plurality M of outputs and being operable to produce in its said outputs M simultaneous representations of the applied input signal having equal time-delay differentals NAT therebetween; multiplying means having a ⁇ first plurality N of input circuits operatively associated with said first signal input means and operable to receive said first signal representations and a second plurality M of input circuits operatively associated with said second input means and operable to receive said second signal representations, said multiplying means having a plurality MN of output channels each providing an output signal representative of the product of a different pair of said input signal representations with each such pair including one input signal representation from each of said signal input means.
  • said -first signal input means comprises means for deriving samples of said first input signal, including first signal storage means having N storage locations for said samples; said second input means comprises means for deriving samples of said second input signal, including second signal storage means having M storage locations for said samples; and wherein said system further includes cyclically operable timing control means for actuating each of said sample deriving means at successively delayed times to store successive sample values of the applied input signal as said input signal representations.
  • timing control means includes means for actuating the first sample deriving means to store successive first signal samples separated by intervals AT and means for actuating the second sample deriving means to store second signal samples separated by time intervals NAT.
  • said first signal input means includes a first signal sampling gate, an analog to digital converter unit and an N stage shift register
  • said second signal input means includes a second signal sampling gate, a second analog to digital converter, and an M stage shift register
  • said system further including cyclically operable timing control means coupled with said first and second signal input means causing said first input signal to be sampled periodically at a frequency f and said second input signal to be sampled at a frequency of f/N
  • said first converter unit 1 being coupled with said N stage shift register and operable to digitize samples of said first signal and apply the same sequentially to said N stage shift register
  • said second converter unit being coupled with said M stage shift register and operable to digitize samples of said second signal and apply the same sequentially to said M stage shift register.
  • the system defined in claim 24 including means applying N samples of said first signal derived from said N stage register and M samples of said second signal derived from said M stage register simultaneously to said multiplying means.
  • said input means include respective pluralities of N and M readout means coupled respectively with said N and M outputs and cooperable in synchronous relationship with the respective sample deriving means to apply to said multiplying means in each time interval N/f respective pluralities of N and M stored samples of said signals bearing the same storage-time relationship to each other as the plurality of samples applied in each preceding cycle.
  • a signal sampling and multiplying system for a correlation computing apparatus comprising in combination: first signal sampling means including a first signal sampling gate adapted to sample a first signal periodically at time intervals AT apart; an N stage signal storage unit; first circuit means coupled with said first sampling means and with said storage unit adapted to apply N output signals representing N samples from said sampling means to said storage unit, said storage unit having N output circuits for providing simultaneously N output signals representing N signal samples; second signal sampling means including a second signal sampling gate adapted to sample a second signal periodically at time intervals NAT apart; a multistage signal storage unit having M stages; second circuit means coupled with said second sampling means and with said M stage storage unit adapted to :fill said M stage storage unit with M signals representing M signal samples which are NAT time intervals apart; and a plurality of multiplication circuit means coupled with said M stage storage unit and with said N stage storage unit for simultaneously multiplying the signal in each stage of said N stage unit times a plurality of the stages in said M stage unit.
  • said first circuit means includes an N stage shift register for receiving samples from said first sampling means sequentially and shifting the same therein, signal transfer means coupling each stage of said shift register with a stage of said N stage storage unit, and means operating said signal transfer means at time intervals NAT apart.
  • each of said signal sampling means includes an analog to digital converter, wherein the signals stored in said storage units are ⁇ digital signals representing the input signal samples, and wherein said multiplication circuit means include MN digital input-analog output multiplication units.
  • said multiplication circuit means includes a number MN of simultaneously operable multiplication units each adapted to multiply a different pair of input signals with each said pair including one signal from one stage of said M stage storage unit and one signal from said N stage storage unit, to thereby provide MN products simultaneously.
  • N equals VVone
  • ⁇ only every Kth-stage of said M stage unit is coupled with said multiplication circuit means and wherein the output of said N stage unit is paired with each said Kth stage Where K is a number greater than one.
  • each of said signal sampling means includes an analog to digital converter, N equals one, said N stage storage unit comprises a single stage 4multibit register, said M stage storage unit comprises a multistage multibit shift register, and each Kth stage of said M stage unit has a signal output circuit paired with the output circuit of said N stage unit as an input for said multiplication means, where K is a number greater than one.
  • each stage of said M stage unit other than each said Kth stage is a magnetic signal storage stage.
  • the said first input means includes N outputs and means for establishing equal time-delay differentials Ar between the signal representations in said N outputs
  • the second input means includes M outputs and means for establishing equal timedelay differentials NA-r between the signal representations in said M outputs
  • said multiplying circuit means includes a plurality of MN of said output channels providing output signals representative of respective products of signal representations corresponding to time-delay differentials separated by increments A-r.
  • a system for computing a correlation function of electrical signals appearing at first and second circuit conductors comprising: sampling means coupled with said conductors for selectively sampling the respective signals at the first and second conductors at time-related successive intervals, signal storage means coupled t0 said sampling means for storing such samples individually to represent said signal values, multiplication circuit means coupled with said storage means for cross multiplying each of the plurality of time-related time successions of signal values in the first storage means by each of the plurality of time-related time successions of signal values in the second storage means, the latter time successions of values bearing a predetermined time relationship to the former, to produce a set of circuit responses representing the Cartesian product of the cross-multiplied signal values, and signal averaging circuit means coupled with said multiplication circuit means for electrically time averaging the respective responses substantially continuously to represent respective correlation coefficients defining said correlation function.

Landscapes

  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Software Systems (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Automation & Control Theory (AREA)
  • Evolutionary Computation (AREA)
  • Fuzzy Systems (AREA)
  • Complex Calculations (AREA)
  • Electronic Switches (AREA)
  • Amplifiers (AREA)
  • Analogue/Digital Conversion (AREA)
  • Measuring Volume Flow (AREA)
US552994A 1963-01-17 1966-04-04 Multichannel correlator system Expired - Lifetime US3514585A (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US25618763A 1963-01-17 1963-01-17
US55299466A 1966-04-04 1966-04-04

Publications (1)

Publication Number Publication Date
US3514585A true US3514585A (en) 1970-05-26

Family

ID=26945190

Family Applications (1)

Application Number Title Priority Date Filing Date
US552994A Expired - Lifetime US3514585A (en) 1963-01-17 1966-04-04 Multichannel correlator system

Country Status (5)

Country Link
US (1) US3514585A (fr)
JP (1) JPS4216148B1 (fr)
DE (1) DE1474101A1 (fr)
FR (1) FR1401210A (fr)
GB (2) GB1065248A (fr)

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3621223A (en) * 1969-09-26 1971-11-16 Atomic Energy Commission High-frequency signal correlator
US3626168A (en) * 1969-07-07 1971-12-07 Boeing Co Measuring correlation, amplitude probability and probability density distributions, and signal response averaging
US3660647A (en) * 1969-12-24 1972-05-02 Us Navy Automatic signal delay tracking system
US3689750A (en) * 1969-12-13 1972-09-05 Philips Corp Phase-independent digital correlator for use in radar systems
US3701894A (en) * 1970-09-11 1972-10-31 Nasa Apparatus for deriving synchronizing pulses from pulses in a single channel pcm communications system
US3717756A (en) * 1970-10-30 1973-02-20 Electronic Communications High precision circulating digital correlator
US3906213A (en) * 1973-03-27 1975-09-16 Thomson Csf Correlation system for delay measurement
US4312040A (en) * 1970-09-09 1982-01-19 Schlumberger Limited Well log depth aligning
US4422165A (en) * 1981-02-11 1983-12-20 Mobil Oil Corporation Maximum likelihood estimation of the ratio of the velocities of compressional and shear waves
US4905209A (en) * 1988-05-27 1990-02-27 Honeywell Inc. Correlation match filter for passive signal detection
US4954828A (en) * 1989-10-13 1990-09-04 Cincinnati Microwave, Inc. Long range police radar warning receiver
US5049885A (en) * 1989-10-13 1991-09-17 Cincinnati Microwave, Inc. Long range police radar warning receiver
US5079553A (en) * 1989-10-13 1992-01-07 Cincinnati Microwave, Inc. Long range police radar warning receiver
US5134406A (en) * 1989-10-13 1992-07-28 Cincinnati Microwave, Inc. Long range police radar warning receiver with multiple array averaging
RU2177637C2 (ru) * 1999-01-18 2001-12-27 Самарский государственный технический университет Многоканальный знаковый коррелометр

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2100550B2 (fr) * 1970-05-26 1973-05-25 Commissariat Energie Atomique
DE2062589C3 (de) * 1970-12-18 1981-03-12 Siemens AG, 1000 Berlin und 8000 München Verfahren zur Ermittlung der Grundfrequenze eines wenigstens zeitweise periodischen Signales
US6509423B1 (en) 2001-08-21 2003-01-21 Dow Corning Corporation Silicone composition and cured silicone product

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US2908889A (en) * 1955-12-16 1959-10-13 Phillips Petroleum Co Computer
US3217289A (en) * 1961-08-18 1965-11-09 Phillips Petroleum Co Signal analysis

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US2908889A (en) * 1955-12-16 1959-10-13 Phillips Petroleum Co Computer
US3217289A (en) * 1961-08-18 1965-11-09 Phillips Petroleum Co Signal analysis

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3626168A (en) * 1969-07-07 1971-12-07 Boeing Co Measuring correlation, amplitude probability and probability density distributions, and signal response averaging
US3621223A (en) * 1969-09-26 1971-11-16 Atomic Energy Commission High-frequency signal correlator
US3689750A (en) * 1969-12-13 1972-09-05 Philips Corp Phase-independent digital correlator for use in radar systems
US3660647A (en) * 1969-12-24 1972-05-02 Us Navy Automatic signal delay tracking system
US4312040A (en) * 1970-09-09 1982-01-19 Schlumberger Limited Well log depth aligning
US3701894A (en) * 1970-09-11 1972-10-31 Nasa Apparatus for deriving synchronizing pulses from pulses in a single channel pcm communications system
US3717756A (en) * 1970-10-30 1973-02-20 Electronic Communications High precision circulating digital correlator
US3906213A (en) * 1973-03-27 1975-09-16 Thomson Csf Correlation system for delay measurement
US4422165A (en) * 1981-02-11 1983-12-20 Mobil Oil Corporation Maximum likelihood estimation of the ratio of the velocities of compressional and shear waves
US4905209A (en) * 1988-05-27 1990-02-27 Honeywell Inc. Correlation match filter for passive signal detection
US4954828A (en) * 1989-10-13 1990-09-04 Cincinnati Microwave, Inc. Long range police radar warning receiver
US5049885A (en) * 1989-10-13 1991-09-17 Cincinnati Microwave, Inc. Long range police radar warning receiver
US5079553A (en) * 1989-10-13 1992-01-07 Cincinnati Microwave, Inc. Long range police radar warning receiver
US5134406A (en) * 1989-10-13 1992-07-28 Cincinnati Microwave, Inc. Long range police radar warning receiver with multiple array averaging
RU2177637C2 (ru) * 1999-01-18 2001-12-27 Самарский государственный технический университет Многоканальный знаковый коррелометр

Also Published As

Publication number Publication date
FR1401210A (fr) 1965-06-04
DE1474101A1 (de) 1970-03-05
JPS4216148B1 (fr) 1967-09-02
GB1065248A (en) 1967-04-12
GB1065249A (en) 1967-04-12

Similar Documents

Publication Publication Date Title
US3514585A (en) Multichannel correlator system
US3626168A (en) Measuring correlation, amplitude probability and probability density distributions, and signal response averaging
US3544775A (en) Digital processor for calculating fourier coefficients
EP0132807B1 (fr) Générateur de grandeurs représentatives d'un signal
Johnson Constructions of particular random processes
US2643819A (en) Apparatus for computing correlation functions
US3404261A (en) Correlation apparatus for computing time averages of functions
US3466431A (en) D.c. power spectrum and fourier transform analyzer
US4071906A (en) CTD programmable convolver
US3697703A (en) Signal processing utilizing basic functions
US3555258A (en) Multicorrelator for analogue signals employing pulse width-amplitude multiplication and operating in real time
US3359409A (en) Correlation function apparatus
US3639695A (en) Systems for processing frequency modulated signals
US3745317A (en) System for generating the fourier transform of a function
US4296374A (en) Wideband digital spectrometer
US3548107A (en) Signal processing apparatus for multiplex transmission
WO1991004496A1 (fr) Methode et appareil pour le mesurage de frequence
RU2174706C1 (ru) Устройство для определения плотности распределения вероятностей случайного процесса
MacKay High-speed electronic-analogue computing techniques
SU813458A1 (ru) Устройство дл вычислени корре-л циОННыХ фуНКций
US3373359A (en) One-sample real time cross-correlator
SU859950A1 (ru) Цифровой спектроанализатор
SU917119A1 (ru) Анализатор комплексного спектра
SU789866A1 (ru) Спектральный анализатор
RU2007856C1 (ru) Демодулятор