Systolic array implementation of recursive bispectrum estimation

Presents the systolic array implementation of recursive higher order spectrum in which the bispectrum estimate is updated every data sample. Forward and reverse sequence running Fourier transforms are first systolically realized. The results are then used to drive a second systolic array whose outputs represent the FT of the data third order moment. The proposed systolic array has no global communications with a number of PEs independent of the size of the employed two dimensional lag window.