Background of the Invention
1. Field of the Invention
The present invention relates to adaptive signal processors for signal classification, identification and enhancement.
2. Prior Art
U.S. Pat. No. 4,315,319, "Nonlinear Signal Processor" by inventors Stanley A. White and V. A. Vitols, issued Feb. 9, 1982 is of interest because it describes a nonlinear signal-processing technique to combat the problem of nonlinear distortion by mapping an input signal having any given amplitude-density function into a new signal with any given amplitude-density function, subject only to the constraint that the distorting mechanism be a single-valued nonlinear function, the slope of which is everywhere positive. The nonlinear signal-processing technique of that patent has been successfully applied to both speech and image-restoration problems.
U.S. Pat. No. 4,843,583 for "A Non-linear Adaptive Filter" filed Oct. 15, 1985, having a common inventor and assignee is of interest because it depicts an adaptive computation which is used in this invention.
Summary of Invention
A computing element for use in a neural network is described. A neural network uses an array of the invention computing elements arrayed into (i) rows and (j) columns. Each computing element has at least K (K>1) input terminals and at least one output terminal y.sub.ij.
The input terminals for each respective computing element are coupled to receive a sequence of concurrent input signals on K parallel input lines representing a parallel input-vector signal S.sub.ij =(s.sub.ij1, s.sub.ij2, s.sub.ij3, . . . , s.sub.ijK).sup.T.
Each computing element has a weighting function computational block or means that is responsive to the concurrent input signals for computing a K-dimensional weighting coefficient representing a K-dimensional weighting-coefficient vector, W.sub.ij =(w.sub.ij1, w.sub.ij2, . . . , w.sub.ijk).sup.T.
Each weighting function means also computes a scalar activation signal u.sub.ij by forming the vector inner product of the input vector S.sub.ij and the weighting-coefficient vector, W.sub.ij. The weighting function means also provides a backpropagated vector signal X.sub.ij =p.sub.ij times W.sub.ij, where p.sub.ij is a backpropagated signal characterized to provide a performance-error signal for the computing element.
Each computing element also has a nonlinear processor or nonlinear processor means that is coupled to receive and respond to the activation signal u.sub.ij for mapping successive values of u.sub.ij through a nonlinear mapping function, M.sub.ij, to provide a single-valued output signal, y.sub.ij. The nonlinear processor also responds to an m-dimensional backpropagated signal vector, X.sub.ij+1 =(x.sub.1,j+1,p ; x.sub.2,j+1,p ; x.sub.3,j+1,p ; . . . x.sub.m,j+1,p).sup.T from a subsequent layer of computing elements to provide a measure of the performance error of the computing element for computing the backpropagated signal, p.sub.ij. Each member of X.sub.i,j+1 is the pth member of the backpropagated error vector from each computing element in column j+1.
In a more particular alternative embodiment of the invention computing element, the nonlinear processor has a memory means that is responsive to the activation signal u.sub.ij for computing a single-valued monotonic activation function F(u.sub.ij). The activation function F(u.sub.ij) is activated by successive u.sub.ij signals to provide corresponding scalar signals e.sub.ij, for all (0.ltoreq.e.sub.ij .ltoreq.1).
Each e.sub.ij signal has a value corresponding to a point at u.sub.ij on the distribution function curve of F(u.sub.ij). The distribution function of a numerical-valued phenomenon has a value at its output which is equal to the probability that the sample value of the input phenomenon is less than or equal to its actual value.
In an even more particular embodiment, an inverse output distribution function means G(e.sub.ij) responds to each e.sub.ij signal as an address to a random access memory to provide a corresponding y.sub.ij signal at the output of the computing element. Each y.sub.ij signal delivered to the output has a value that is the inverse of the desired output-distribution function. Each e.sub.ij signal, which is a probability distribution function value, is used as an address to select a value corresponding to the inverse of the desired output distribution function.
A derivative output density function, g.sub.ij responds to each e.sub.ij signal to provide a corresponding .alpha..sub.ij signal. Each respective .alpha..sub.ij has a value that corresponds to the derivative of the inverse of the probability-distribution function of the desired output signal.
In another alternative embodiment, the memory means further comprises a computer means that responds to each successive value of u.sub.ij by computing the probability-density function f(u.sub.ij). The probability-density function f(u.sub.ij) is activated by successive u.sub.ij signals to provide corresponding scalar signals v.sub.ij. Each v.sub.ij signal has a value corresponding to a point on the probability-density function computed from all past values of u.sub.ij.
In yet another alternative of the computing element, each value of the activation signal u.sub.ij is used to select an address in an f(u.sub.ij) memory location having an address corresponding to the value of the activation signal u.sub.ij to read a value from the f(u.sub.ij) function corresponding to the probability of u.sub.ij at the same address location.
Each e.sub.ij signal selects successive addresses in the g.sub.ij function to readout .alpha..sub.ij signal values for each successive e.sub.ij signal. Each .alpha..sub.ij signal value selected corresponds to the derivative of the inverse of the desired output distribution function.
A backpropagation summer responds to a backpropagated signal X.sub.ij+1 through "m" input ports. The X.sub.ij+1 signal represents an m-dimensional vector; X.sub.ij+1 =(x.sub.1,j+1,p ; x.sub.2,j+1,p ; x.sub.3,j+1,p ; . . . x.sub.m,j+p,1).sup.T. The backpropagation summer sums the elements of each successive X.sub.ij signal to provide a scalar .beta..sub.ij signal.
A combiner multiplies the value of each successive f-mapped activation signal, v.sub.ij times each successive .alpha..sub.ij signal, times each successive .beta..sub.ij signal to provide a corresponding p.sub.ij signal.
Brief Description of the Drawings
FIG. 1 is a block diagram of a neural network shown as an interconnected array of computing elements.
FIG. 2 is a block diagram of a neural network computing element having a first block representing a weighting function and a second block representing a nonlinear function.
FIG. 3 is an expanded block diagram of a computing element.
FIG. 4 is a block diagram of a component of the weighting function.
FIG. 5 is a table of values illustrating the values present on particular signal lines within the block diagram of the component of the weighting function of FIG. 4.
FIG. 6 is a block diagram of the computer for calculating values for the f.sub.ij and F.sub.ij RAMs within the memory system in the block diagram of the nonlinear function of FIG. 4.
Preferred Embodiment
FIG. 1 is a schematic diagram 8 having a neural network of two layers or columns of computing elements. In general, a neural network is an interconnected array of computing elements such as those represented by circles 10, 11, 12, 13, 14, 15, 16. A computing element, such as computing element 10 combines a parallel input signal vector S.sub.ij =(s.sub.ij1, s.sub.ij2, s.sub.ij3, . . . , s.sub.ijK).sup.T on signal lines 17, 18, 19, 20 with a coefficient vector W.sub.ij =(w.sub.ij1, w.sub.ij2, w.sub.ij3, . . . , w.sub.ijK) developed within the respective computing element to generate an output scalar, y.sub.ij =M.sub.ij (u.sub.ij) where u.sub.ij is referred to as the activation signal and the values of u.sub.ij are obtained from the product of s.sub.ij.sup.T *w.sub.ij. The letter "T" is used to signify the transpose of the vector quantity. The components of W.sub.ij are adjustable according to the minimization of some performance criterion Q, such as the minimizing of the sum of the squares of the errors of all the outputs from the neural network. The nonlinear mapping function M.sub.ij, is a single-valued monotonically nondecreasing function, which may be fixed or variable.
The columns of computing elements forming the neural network are designated by the subscript j where j=1, 2, . . . , J. There are I(j) elements in each column in rows i=1, 2, . . . , I(j). The computing element at location (i,j) has inputs s.sub.ijk, where k extends from k=0, 1, 2, . . . , K(j-1). Inputs s.sub.ij0 =1 for all i and all j. Each computing element has an output y.sub.ij.
Typically, the parallel input signal vector, S.sub.0 =(s.sub.01, s.sub.02, . . . , s.sub.0K).sup.T, is coupled to the K(0) inputs of each computing element in the input, or first, layer of computing elements. The signal input to the kth input node of the ith computing element in the first column is s.sub.i1k =s.sub.0k. The output signal from the ith computing element in column 1 is y.sub.i1. The input signal to the kth input of the ith computing element in the second column is s.sub.i2k =y.sub.k1. In general, the input to the kth input of ith computing element in the ith column is s.sub.ijk =y.sub.kj-1.
For the purpose of indexing, the value of K(j)=I(j-1); J=0 for the source column of signals; j=1 for the input column of computing elements and j=J for the output column of computing elements. Each computing element is represented in FIG. 4 by a circle with a multiplicity of bidirectional inputs which connect to every output node of the previous layer and a bidirectional output which fans out to an input of every element in the next layer. Although a two-layer network is shown in FIG. 1, neural networks can be formed that include higher numbers of layers of computing elements.
In a backpropagation network, all signal paths are bilateral. The error in the performance of the neural network is measured, the contribution of each w.sub.ijk to that error is determined, and a correction signal is propagated back from the output to be used in the computation of the value of each w.sub.ijk.
FIG. 2 shows a computing element 10 to be composed of two sections or functions: a weighting function 21, which generates activation signal, u.sub.ij, on signal line 22 and a nonlinear function 24, which generates y.sub.ij, and outputs y.sub.ij on signal line 25. Parallel input signals, s.sub.ijk are supplied on signal lines 17, 18, 19, 20 and the backpropagated error signals are outputted on signal lines 26, 27, 28, 29. An m-dimensional backpropagated signal vector, X.sub.ij+1 =(x.sub.1,j+1,p ; x.sub.2,j+1,p ; x.sub.3,j+1,p ; . . . x.sub.m,j+1,p).sup.T is received by summer 30 on signal lines 32, 33, 34 and 35. Summer 30 is shown coupling the summed backpropagated signal .beta..sub.ij to the nonlinear function 24 on signal line 38. The backpropagated signal p.sub.ij is shown being coupled from the nonlinear function to the weighting function 21 by signal line 40.
Several neural network theorists have speculated that one of the functions of the nonlinearity within a computing element is to map the inner product of the input vector with its weighting vector to obtain the activation signal, u.sub.ij, on signal line 22, in a manner to control its amplitude-density function.
If the activation signal, u.sub.ij, has Gaussian statistics and is scaled such that the full dynamic range is used, then mapping u.sub.ij through the "conventional" nonlinearity (the logistic function, y=[(1+exp(-u)].sup.-1) yields a signal output with amplitude statistics that are nearly uniformly distributed.
Because of the greater entropy of such a signal, the network tends to be "livelier" and to respond faster. If a computing element receives many such inputs which are "somewhat" independent, then the next inner product will tend to be "somewhat" Gaussian, and the amplitude-density function controlling or shaping procedure will repeat. If the assumed amplitude statistics are violated, performance becomes comparatively more sluggish.
The above considerations lead to the embodiment of FIG. 3, which provides a more responsive solution. Phantom box 42 represents a computer, within a memory system 82, for measuring the amplitude statistics of the inner products, i.e. the activation signal u.sub.ij, on signal line 22, and actively maps them into the desired statistics through nonlinear function 24.
The embodiment of the computing element 10 of FIG. 3 is a "doubly-adaptive" backpropagation computing element. The configuration is termed "doubly adaptive" because both weighting function 21 and nonlinear function 24 are adaptive. The weighting function 21 provides a first level of adaptability by adjusting the gain of the parallel input signals, s.sub.ijk on signal lines 17, 18, 19, 20 and the amplitude of the backpropagated error signals on signal lines 26, 27, 28, 29.
The nonlinear function 24 provides a second level of adaptability by adjusting the gain of the bidirectional input signal comprising parallel input signals, s.sub.ijk on signal lines 17, 18, 19, 20 and the backpropagated error signal x.sub.ijk according to a dynamic performance measure which is transmitted to weighting function 21 by the backpropagated signal p.sub.ij. The nonlinear function 24 provides a third level of adaptability by adjusting the shape of the nonlinear mapping function F.sub.ij, provided by the data stored in RAM1 46.
The computer within phantom box 10 is a bidirectional version of the nonlinear adaptive filter described in U.S. Pat. No. 4,843,583, issued Jun. 27, 1989 for "A Nonlinear Adaptive Filter" to S. A. White and V. A. Vitols, and having a common assignee, which is incorporated herein by reference. Two modifications are incorporated into the nonlinear adaptive filter of the '583 patent.
The first modification to the nonlinear adaptive filter of the '583 patent is to replace all frequency-shaping functions within the adaptive loops with unity gain paths; i.e. by replacing each H(z) term with a unity gain path. The second modification is to remove the delay elements from the tapped-delay line and drive the tap points directly with the components of the signal vector which is the input to that neural element. These changes to the nonlinear adaptive filter of the '583 patent provide a unidirectional version of the computing element 10, but do not provide a computing element capable of backpropagation and therefore do not provide a computing element suitable for use in a neural network.
The mathematics of the modified nonlinear adaptive filter were combined with the established mathematics of backpropagation networks to assist in formulating a block diagram of the invention computing element 10.
FIG. 3 is a detailed block diagram of the backpropagation, doubly adaptive computing element. It is "all inclusive" in that any function which may be required at any level in a backpropagation network element has been included. Any known backpropagation neural network may be mechanized entirely using only subsets of the functions and topology of FIG. 3. For each parallel input signal, s.sub.ijk, to the weighting function 21, there is a backpropagated signal, x.sub.ijk =p.sub.ij *w.sub.ijk.
A bilateral scaling unit, such as the weighting function 21, has weighting function elements, such as those in FIG. 3, represented by circles 48, 50, 52, 54, 56 that generate three quantities: (1) a set of backpropagated error signals, x.sub.ijk, to all elements with which it is connected in the preceding stage; (2) a coefficient weighted signal, d.sub.ijk =s.sub.ijk *w.sub.ijk, to be summed with like feedforward signals to generate the activation signal, u.sub.ij, on signal line 22, to the local nonlinear function 24; (3) an updated concurrent weighting signal w.sub.ijk for use within the weighting function element.
There is a bilateral coupling between the weighting function and the nonlinear function. The forward signal is u.sub.ij, which is simply the sum of the coefficient weighted signals, the d.sub.ijk 's. The backpropagated signal is p.sub.ij on signal line 40 which is generated within the nonlinear function 24.
FIG. 3 shows an embodiment of the weighting function 21 that responds to the concurrent parallel input signals S.sub.ij =(s.sub.ij1, s.sub.ij2, . . . , s.sub.ijK).sup.T and computes a K-dimensional weighting coefficient W.sub.ij, computes the products d.sub.ijk =s.sub.ijk *w.sub.ijk, and computes u.sub.ij =d.sub.ij1 +d.sub.ij2 +. . . , +d.sub.ijK in weighting function summer 58 which is output on signal line 22. Each sequential value of u.sub.ij is obtained by computing a K-dimensional weighting-coefficient vector, W.sub.ij =(w.sub.ij1, w.sub.ij2, . . . , w.sub.ijK).sup.T in weighting function elements 48, 50, 52, 54, 56, and by forming the vector inner product of the input signal S.sub.ij vector elements and the weighting-coefficient vector, W.sub.ij by performing the scalar products to generate d.sub.ijk within each of the respective weighting function elements and summing these components within element 58.
The weighting function shown in phantom block 21 also provides a set of external backpropagated error signals x.sub.ijk =P.sub.ij times w.sub.ijk on signal lines 26, 27, 28 and 29 from the output backpropagated signal terminals. The backpropagated signal p.sub.ij, on signal line 40 is characterized to provide a performance error of all columns of computing elements subsequent to this computing element weighted by the instantaneous gain of the computing element.
FIG. 4 is a block diagram of a weighting function element such as that represented by circle 50 in FIG. 3. The parallel input signal s.sub.ij1 is sampled on signal line 17 and is multiplied by the concurrent weighting signal w.sub.ij1 on signal line 60 within the first multiplier 62 to generate the weighted coefficient signal, d.sub.ij1, on signal line 64. The backpropagated signal p.sub.ij on signal line 40 is multiplied by w.sub.ij1 in second multiplier 66 to generate the backpropagated error signal, x.sub.ij1, on signal line 26.
The weighting signal w.sub.ij1 is generated in the accumulator loop 68 comprising summer 70 and delay element 72. Summer 70 combines an update signal q.sub.ij1 on signal line 74 and the concurrent weighting signal, w.sub.ij1, on signal line 60 to produce the next sample weighted signal on signal line 76. The sample weighted signal drives the delay element 72. The update signal, q.sub.ij1, is the output of third multiplier 78. The third multiplier 78 is mechanized to multiply the parallel input signal, s.sub.ij, on signal line 17, times an internally stored convergence coefficient C1 provided on signal line 80, times the internally backpropagated signal, p.sub.ij, on signal line 40.
FIG. 5 provides an exemplary set of calculations to illustrate the performance of the weighting function element of FIG. 4 for an assumed value of C1=1.0, and an initial set of (GIVEN) initial conditions.
Referring again to FIG. 3, phantom block 24 is a block diagram of a nonlinear function 24 that maps successive values of the activation signal u.sub.ij through nonlinear mapping functions F.sub.ij, provided by the data in RAM1 46, and G.sub.ij provided by the data in ROM1, 81 to provide a single-valued output signal y.sub.ij at output signal line 25.
The computing element 10 and the nonlinear function 24, through the summer 30, is coupled to receive the m-dimensional backpropagated signal vector, X.sub.ij+1 =(x.sub.1,j+1,p ; x.sub.2,j+1,p ; x.sub.3,j+1,p ; . . . x.sub.m,j+1,p).sup.T on signal lines 32-35. Each component of the X.sub.ij+1 backpropagated signal vector originates from the pth weighting function of each computing element in column j+1. The X.sub.i,j+1 signal vector is characterized to provide a measure of the performance error of all computing elements subsequent to the computing element 10 such as computing elements 14, 15, 16.
Phantom block 82 represents a memory means connected to receive the activation signal u.sub.ij via signal line 22. The memory means 82 computes and stores a single-valued monotonic activation function F(u.sub.ij). The activation function F(u.sub.ij) is addressed by successive u.sub.ij signals to provide a sequence of corresponding scalar signals e.sub.ij, (0.ltoreq.e.sub.ij .ltoreq.1). Each e.sub.ij signal on signal line 84 has a value that is read out on signal line 84 from an address that corresponds to the value of u.sub.ij. The value of e.sub.ij that is obtained from the memory means 82 in response to each u.sub.ij address corresponds to the value or point on the distribution function curve compiled from all past values of the corresponding u.sub.ij signal on signal line 22.
ROM1, 81 in FIG. 3, is a first ROM that is preloaded with data that represents the inverse of the desired output distribution function of the output signal y.sub.ij on signal line 25. ROM1 is a read only memory that responds to each e.sub.ij signal as an address signal by providing a corresponding y.sub.ij output signal on output signal line 25. Each y.sub.ij output signal that is read out of the first ROM1 81 has an inverse probability value selected in response to the value of the e.sub.ij signal. The inverse probability values of G(e.sub.ij) are initialized before operation, G(e.sub.ij), by programming to correspond to the desired output distribution function of the output signal y.sub.ij.
ROM2, 86 in FIG. 3, is a second ROM, that is loaded with data that represents the derivative of the inverse of the output density function. The function g.sub.ij =dG.sub.ij (e.sub.ij)/de.sub.ij by definition. ROM2 is also addressed by values of e.sub.ij. Each e.sub.ij address signal commands the second ROM 86 to read out a corresponding .alpha..sub.ij signal on signal line 88. Each .alpha..sub.ij signal that e.sub.ij reads out has a value that corresponds to a point on the derivative of the inverse of the probability distribution function of the desired output signal y.sub.ij on signal line 25.
Dotted block 42 represents a computer that responds to each successive value of u.sub.ij on signal line 22 by computing and updating the corresponding stored value of the probability density function f(u.sub.ij) in RAM2, 90. The probability density function f(u.sub.ij) is activated by successive u.sub.ij signals which serve as address signals for RAM2 90, which responds by outputting a value from the corresponding address that represents a scalar signal v.sub.ij on signal line 91. Each v.sub.ij signal has a value that corresponds to a point on the probability density function of u.sub.ij computed from all past values of u.sub.ij on signal line 22.
Circle 30 in the lower right side of FIG. 3, represents a backpropagation summer connected to signal lines 32-35 to receive the vector components of the input backpropagated signal X.sub.ij+1. The X.sub.ij+1 signal represents an m-dimensional vector, X.sub.ij+1 =(x.sub.1,j+1,p ; x.sub.2,j+1,p ; x.sub.3,j+1,p, . . . x.sub.m,j+1,p).sup.T. The backpropagation summer 30 sums the values of the vector components of each successive X.sub.ij+1 signal to provide a summed backpropagated signal .beta..sub.ij on signal line 38.
In the design of a neural network, such as the network depicted in FIG. 1, and where j=J, we are referring to computing elements in the last column. The nonlinearity function characterized by block 92 in FIG. 3 in computing elements in the Jth columns is typically a function selected to control the performance criterion, Q, described above, to which the unit adapts.
The shape of the nonlinearity, such as a signum function, is predetermined to be the derivative of the performance criteria, Q, with respect to the error at the output. The error at the ith output is .beta..sub.ij where the inputs to the summer 30 are defined as:
If j=J, then y.sub.ij =.beta..sub.ij and element 92 is simply a unity gain function.
Block 94, in FIG. 3, represents a combiner that multiplies the value of each successive .gamma..sub.ij signal times each successive .alpha..sub.ij signal times each successive v.sub.ij signal to provide a corresponding backpropagated signal p.sub.ij.
The e.sub.ij signal on signal line 84 selects successive addresses in the ROM2, 86, storing the values generated by operation of the g.sub.ij function, to readout .alpha..sub.ij signal values on signal line 88. Each .alpha..sub.ij signal value corresponds to the derivative of the inverse of the desired output distribution function corresponding to the value of the e.sub.ij.
The combiner or multiplier, 94 scales signal .gamma..sub.ij on signal line 93 from the conditional nonlinearity 92 by the instantaneous feedforward gain, i.e. .alpha..sub.ij *v.sub.ij, of the nonlinear function 24 to provide the backpropagated signal p.sub.ij on signal line 40. The backpropagated signal p.sub.ij =.gamma..sub.ij (.alpha..sub.ij *v.sub.ij), where (.alpha..sub.ij *v.sub.ij) is the instantaneous feedforward gain of the nonlinear function 24.
The expression for the gain of the nonlinear function 24 is obtained from FIG. 3 by observing that:
The output of block 81 is:
The output of block 46 is:
In general, the gain of any function having an output y.sub.ij and an input signal u.sub.ij is obtained from the relationship:
Differentiating equations (1) and (2) provides:
Block 86 is defined as:
and block 90 is defined as:
By combining (4) and (6)
and by combining (5) and (7)
From (8) and (9)
dy.sub.ij =g.sub.ij (e.sub.ij)*f.sub.ij (u.sub.ij)*du.sub.ij(10)
But by definition:
and
By combining (10), (11) and (12):
From equations (3) and (13) therefore:
FIG. 6 is a block diagram of a particular embodiment of the computer shown in FIG. 3 as phantom box 42. The computer is characterized to respond to each successive value of u.sub.ij on signal line 22 by computing and updating the single-valued monotonic activation function F(u.sub.ij) via signal path 95.
The operational process of this embodiment is explained by examining the sequence of steps that occur with the selector switches 96a, 96b, 96c in each of two possible positions; i.e. position I for the input cycle and position U for the update cycle.
With the switches in the I position, the activation signal u.sub.ij passes along lines 22, through switch 96a to signal line 97 and then to the address bus of the LSB RAM 98 and the address bus for the MSB RAM 99. These RAMS are part of the same memory; the LSB RAM 98 being storage for the least significant bits and MSB RAM 99 being storage for the most significant bits of the same data word or value.
As each RAM address is selected by u.sub.ij, the value of the data word stored in the RAM MSB 99 at that address is read out on lines 100 and fed to the summer 102 and to the input of the -2.sup.-E block 104. With switch 96b in the I position, 1 is added to the value fed to the summer 102 from RAM MSB 99 via signal line 100. The number stored at the RAM MSB 99 and RAM LSB 98 address is reloaded back into the same address with a 1 bit increase.
The switches 96a, 96b, 96c are transferred to the U position for update. The counter, 106 is advanced one count to select the next address in sequence and the RAM MSB 99 and the RAM LSB 98 are addressed by this address via switch 96a in position U. The value of the number in the RAM MSB 99 at that sequential address is then fed via line 100 to the input of the -2.sup.-E block 104. This block, in conjunction with the direct feedback input to the adder 102, multiplies the value received on signal line 100 by 1-2.sup.-E, a number of less than 1, thereby diminishing the magnitude of the number received. The diminished value of the number that was in the RAM MSB 99 address location is then fed via switch 96b in the U position to the summer 102 and reloaded as a new value into the RAM MSB 99 and RAM LSB 98 at the selected address.
The value selected for E in the mechanization of the -2.sup.-E block 104 function is application dependent and establishes the time constant of the run down of the value in the RAM MSB 99 and RAM LSB 98 selected location addressed by the counter 106 via signal lines 107, 108, switch 96a and signal line 97 in the RAM MSB 99 and RAM LSB 98.
The 1-2.sup.-E multiplier, mechanized by block 104, line 100, and summer 102, operates as the feedback gain in a low pass filter and keeps the values of the data stored in locations of the RAM MSB 99 and RAM LSB 98 from overflowing. The value of the output of RAM MSB 99 is also fed to the accumulator 110, via switch 96c in the U position. The accumulator 110 operates with delay 111 to accumulate the amplitude density function F.sub.ij and output the value of the amplitude distribution function on signal line 95 to the appropriate address of RAM1.
The output of the RAM MSB 99 is also fed via signal line 112 to the same address of RAM2. Phantom block 114 characterizes the memory and control means function of the computer means 42, for receiving the activation signal, scaling, storing and controlling the process to provide update signals for the amplitude-density storage function 90 and the amplitude distribution storage function 46.