scispace - formally typeset
Search or ask a question

Showing papers on "Fast Fourier transform published in 1973"


Journal ArticleDOI
TL;DR: A lower bound for the number of additions necessary to compute a family of linear functions by a linear algorithm is given when an upper bound c can be assigned to the modulus of the complex numbers involved in the computation.
Abstract: A lower bound for the number of additions necessary to compute a family of linear functions by a linear algorithm is given when an upper bound c can be assigned to the modulus of the complex numbers involved in the computation. In the case of the fast Fourier transform, the lower bound is (n/2) log2n when c = 1.

115 citations


Journal ArticleDOI
TL;DR: A diagrammatic representation of mixed radix and highest radix FFT algorithms is derived, and two broad classes of FFT hard-ware are explored, from the point of view of speed, parallelism, radix number, and type of memory.
Abstract: The fast Fourier transform algorithm is derived by means of successive fracturing of one-dimensional data strings into two-dimensional arrays. Using this formulation, a diagrammatic representation of mixed radix and highest radix FFT algorithms is derived. Using this representation, two broad classes of FFT hard-ware are explored, from the point of view of speed, parallelism, radix number, and type of memory.

106 citations


Journal ArticleDOI
TL;DR: There are many kinds of digital calculation methods for transient phenomena on a electrical system and the most highly developed technique among them is the iterative method of travelling waves which has been given as the O'Alembert's solution of the wave equation as mentioned in this paper.
Abstract: There are many kinds of digital calculation methods for transient phenomena on a electrical system. The most highly developed technique among them is the iterative method of travelling waves which has been given as the O'Alembert's solution of the wave equation. The technique was first devised by Schnyder! and Bergeron\" and developed by Frey3 , Barthold , McElroy 5, Ametani\" and others. The other famous technique is the Fourier transform which is based on the Bernoulli's solution of the wave equation. This technique has been developed by Leg0 , Mullineux\", Wedepohl? and others. This, however, has a severe disadvantage, that is, computer programs are using hundreds of hours of computer time with procedures requiring N 2 operations to compute Fourier transforms of N data points.

85 citations


Journal ArticleDOI
TL;DR: A technique for the digital simulation of multicorrelated Gaussian random processes is described, based upon generating discrete frequency functions which correspond to the Fourier transform of the random processes and then using the fast Fourier Transform algorithm to obtain the actual random processes.
Abstract: A technique for the digital simulation of multicorrelated Gaussian random processes is described. This technique is based upon generating discrete frequency functions which correspond to the Fourier transform of the random processes and then using the fast Fourier transform (FFT) algorithm to obtain the actual random processes. The main advantage of this method over other methods is computation time; it appears to be more than an order of magnitude faster than present methods of simulation. One of the main uses of multicorrelated simulated random processes is in solving nonlinear random vibration problems by numerical integration of the governing differential equations. [This research is supported in part by NASA.]

78 citations


Journal ArticleDOI
TL;DR: In this article, a mode-superposition method by the use of damped free vibration mode shapes is developed for such structural systems and the Fast Fourier Transform (FFT) procedures are used for the dynamic response analysis of linear structural systems.
Abstract: In the dynamic response analysis of extremely complex structural systems in which the damping characteristics of each element are independent, the damping matrix is not always diagonalized by the use of undamped free vibration mode shapes In the present paper, a mode-superposition method by the use of damped free vibration mode shapes is developed for such structural systems It is also shown that the Fast Fourier Transform (FFT) procedures, that are available for the dynamic response analysis of linear structural systems, are used effectively in this mode-superposition method with good accuracy

59 citations


Patent
13 Dec 1973
TL;DR: In this article, the stored data pattern of N complex data points or record samples in the same predetermined order each time after each iteration or pass of the data through an FFT arithmetic unit.
Abstract: Method and apparatus for reducing the hardware requirements and complex control of an in-place Fast Fourier Transform (FFT) processor by a unique repetitive arrangement of the stored data pattern of N complex data points or record samples in the same predetermined order each time after each iteration or pass of the data through an FFT arithmetic unit. This results in the final output sequence of the calculated data points automatically having the same format as the initial input sequence. Four data memories are utilized in pairs in conjunction with the arithmetic unit wherein the calculated samples are stored in successive iterations. For example, the first pair of memory units are selectively addressed for receiving the calculated results in a predetermined order and then used as the input means for consecutively feeding its contents to the arithmetic unit in the next pass at which time the second pair of memories are selectively addressed in the same predetermined manner for storing the calculated data. Thus the memories alternate as input and output memory devices always having the data stored therein in the same predetermined order while having a subsequent serial readout.

28 citations


Journal ArticleDOI
TL;DR: The advantages and drawbacks of Hadamard transform spectroscopy are examined in detail and existing theory, as reported in the literature, is not consistent with the claims made for the technique.
Abstract: Recent articles have claimed a significant S/N advantage of Hadamard transform spectroscopy over Fourier transform spectroscopy. The scanty published data does not support this assertion, and the possibility that the claim is valid in theory is examined. Existing theory, as reported in the literature, is not consistent with the claims made for the technique. The advantages and drawbacks of Hadamard transform spectroscopy are examined in detail.

23 citations


Patent
18 May 1973
TL;DR: In this paper, a special purpose data handling equipment is utilized with a fast Fourier transform algorithm computer, which buffers input data such that each block of input data includes data from the immediately preceding and following data blocks.
Abstract: Special purpose data handling equipment is utilized with a fast Fourier transform algorithm computer. The equipment buffers input data such that each block of input data includes data from the immediately preceding and following data blocks. This input redundancy is useful when input data is processed through a data window which attenuates data at the beginning and end of the data block, since, when this data is transformed and subsequently inversely transformed, simple addition of the output data points removes the effect of the attenuating window.

16 citations


Journal ArticleDOI
TL;DR: It is shown that the Haar transform can be computed using a Cooley-Tukey-type algorithm that is implemented in 2(N−1) additions/subtractions.
Abstract: It is shown that the Haar transform can be computed using a Cooley-Tukey-type algorithm that is implemented in 2(N−1) additions/subtractions. This algorithm is derived by relating the Haar transform to the modified Walsh-Hadamard transform using a simple bit-reversal scheme.

13 citations


Journal ArticleDOI
TL;DR: In this article, the initial value problem due to a transient distributed transverse force on the plate is solved using Laplace and Fourier transforms, and stress contours for various stresses and times after application of load are obtained for a graphite fiber-epoxy matrix composite plate.

12 citations


Journal ArticleDOI
TL;DR: In this paper, a scattering-matrix theory for the measurement of electroacoustic transducers in fluids is presented, and deconvolution of transverse scanning data is used to obtain complete effective directivity functions, corrected for the effects of the measuring transducers.
Abstract: Recently developed analytical techniques for the measurement of microwave antennas at reduced distances are “translated” into corresponding techniques for the measurement of electroacoustic transducers in fluids. The basic theory is formulated in scattering‐matrix form and emphasizes the use of plane‐wave spectra for the representation of sound fields. This theory, in contrast to those based on asymptotic description of transducer characteristics, is suitable for the formulation and solution of problems involving interactions at arbitrary distances. Two new techniques (in particular) are described: One, utilizing deconvolution of transverse scanning data, which may be taken at distances d much less than the Rayleigh distance dR(≡ a2/2λ), provides a means of obtaining complete effective directivity functions, corrected for the effects of the measuring transducer. Applicability of a (two‐dimensional, spatial) sampling theorem and the “fast Fourier transform” algorithm, which greatly facilitate the necessary...

ReportDOI
01 Mar 1973
TL;DR: In this article, it was shown that the classes of modular rings of integers in which the FFT may be computed is completely characterized by the prime decomposition of the modulus.
Abstract: : The paper presents several results relating modular arithmetic schemes and the Fast Fourier transform. In particular, the classes of modular rings of integers in which the FFT may be computed is completely characterized by the prime decomposition of the modulus. Also, an extension of this result for computation structures similar to modular rings of integers yields a sufficiency hypothesis for the computation of FFT.

Journal ArticleDOI
TL;DR: In this paper, a new technique for the numerical reconstruction of images from a long-wavelength hologram is proposed, which is to calculate images by the fast Fourier transform (FFT) algorithm with the hologram data sampled by a straight line approximation.
Abstract: A new technique for the numerical reconstruction of images from a long-wavelength hologram is proposed. The technique is to calculate images by the fast Fourier transform (FFT) algorithm with the hologram data sampled by a straight line approximation. This sampling technique is convenient to reduce the number of the sampling points, and the image can be reconstructed with less data than the conventional equally spaced sampling method. A one-dimensional hologram is constructed at S band and an image is reconstructed by the proposed method. These results are discussed and compared with the results of the conventional numerical reconstruction.

Journal ArticleDOI
TL;DR: This FFT method, in conjunction with the optimization technique, is found to be an efficient, general procedure to facilitate convolutions of phase time distributions using the FFT.
Abstract: This paper provides a numerical technique for the solution of the fraction of labelled mitoses curve. The essence of the method is to facilitate convolutions of phase time distributions using the FFT. This FFT method, in conjunction with the optimization technique, is found to be an efficient, general procedure.

Journal ArticleDOI
TL;DR: Concepts of group theory are used to explain the difference between two major fast Fourier transform algorithms and the relation between FFT and fast Walsh transforms.
Abstract: Concepts of group theory are used to explain the difference between two major fast Fourier transform (FFT) algorithms and the relation between FFT and fast Walsh transforms.

Proceedings ArticleDOI
24 Sep 1973
TL;DR: The SPS-81 is a very high speed multiprocessor computer designed for signal processing applications which provides the capability to perform signal processing functions with about a 100 fold increase in speed over a conventional computer.
Abstract: The SPS-41 and SPS-81 are very high speed multiprocessor computers designed for signal processing applications. They can either function as stand-alone processors or operate in conjunction with a conventional computer. The SPS-41 provides the capability to perform signal processing functions with about a 100 fold increase in speed over a conventional computer—typically 8.3 milliseconds for a 1024 complex Fast Fourier Transform. The SPS-81 is roughly twice as fast as the SPS-41. This speed is achieved through a highly parallel architecture which is optimized for signal processing operations.


Journal ArticleDOI
TL;DR: The Signal Analysis Program is described, a computer utility which employs the Fast Fourier Transform (FFT) algorithm and is designed to relieve the user from time consuming and tedious calculations.
Abstract: One of the problems in effectively teaching a course in signal and system theory is the difficulty of performing the mathematical operations involved in the application of theory to meaningful problems. In theory, the Fourier transform, convolution and correlation operations usually are readily grasped by the student but the application often proves to be difficult. Performing the integral operations requires painstaking care to get proper limits on the integrals, displacing functions properly, etc. It seems that the student gains little insight to the problem under analysis by these manipulations. Typically, only the simplest type of functions are handled by this analytic approach. This paper describes the Signal Analysis Program, a computer utility which employs the Fast Fourier Transform (FFT) algorithm and is designed to relieve the user from time consuming and tedious calculations. Classes of problems and signal functions, previously deemed too complicated to treat in an undergraduate course, require little more concern than the simplest of problems. The effectiveness and applicability of this utility has been demonstrated by students applying it to many other areas where it is not a course requirement.

A. K. Sinha1
01 Jul 1973
TL;DR: In this article, two methods are described to simulate, on a digital computer, a set of correlated, stationary, and Gaussian time series with zero mean from the given matrix of power spectral densities and cross spectral density.
Abstract: Two methods are described to simulate, on a digital computer, a set of correlated, stationary, and Gaussian time series with zero mean from the given matrix of power spectral densities and cross spectral densities. The first method is based upon trigonometric series with random amplitudes and deterministic phase angles. The random amplitudes are generated by using a standard random number generator subroutine. An example is given which corresponds to three components of wind velocities at two different spatial locations for a total of six correlated time series. In the second method, the whole process is carried out using the Fast Fourier Transform approach. This method gives more accurate results and works about twenty times faster for a set of six correlated time series.

Journal ArticleDOI
TL;DR: An introduction to the generation and use of new transform techniques which have important applications in binary control and processing methods and some applications of the transform are discussed, which include spectral analysis, filtering, non‐linear control and communications uses.
Abstract: An introduction is given to the generation and use of new transform techniques which have important applications in binary control and processing methods. A comparison is made between the fast Fourier transform and the equivalent fast Walsh transform together with the steps required to produce a transform algorithm and computer program. Some applications of the transform are then discussed and which include spectral analysis, filtering, non‐linear control and communications uses. 18 references to current work in these applications areas are included.

11 Apr 1973
TL;DR: It is found that the Hanning window is a reasonable compromise for achieving minimum bias, because in addition to being the optimum for one bandwidth constraint, it is very close to the optima for two other bandwidth constraints.
Abstract: : The time-limited nonnegative data windows that minimize the bias in auto- and cross-spectral estimation of stationary random processes by means of overlapped Fast Fourier Transform (FFT) processing are derived for a variety of constraints. The author constrains the time duration of the data window, or constrain the bandwidth in various ways, such as the equivalent-noise bandwidth, the half-power bandwidth, or the root-mean-square bandwidth. In the three bandwidth-constrained cases, the window duration is adjusted to meet the constraint. We find that the Hanning window is a reasonable compromise for achieving minimum bias, because in addition to being the optimum for one bandwidth constraint, it is very close to the optima for two other bandwidth constraints. The relative merits of the spectral characteristics of the windows are also discussed.

Journal ArticleDOI
T. Walzman1, M. Schwartz
TL;DR: This work develops another projection method to optimize the discrete frequency parameters and shows that for the channels considered, convergence is equivalent to comparable time domain equalizers.
Abstract: In a recent paper [1] we reported on a new mean-square error automatic equalizer utilizing Rosen's gradient projection theorem to optimize parameters in the discrete frequency domain. Here we develop another projection method to optimize the discrete frequency parameters. The algorithm converges (in the mean) for any channel, even in the presence of noise. It is shown that for the channels considered, convergence is equivalent to comparable time domain equalizers. The method makes use of fast Fourier transform (FFT) algorithms for computation of the iteration matrix, the gradient, and the projection operation. Use of the FFT for parameter iterations reduces the necessary computations per parameter to a number proportional to log_{2} M compared to M for a time domain equalizer, where M is the number of equalizer parameters. The method results in fewer computations per parameter for each iteration, but a somewhat slower rate of convergence than the method employing Rosen's gradient projection.

Journal ArticleDOI
TL;DR: In this article, the deconvolution of the response by fast Fourier transform and matrix inversion techniques is dealt with by using matrix inversus-Fast Fourier Transform (FFT) transform.
Abstract: Deals with the deconvolution of the response by fast Fourier transform and matrix inversion techniques.


Journal ArticleDOI
TL;DR: In this article, the analysis and stability properties of a special class of R-L-C networks, where the energy storage elements vary periodically with time and the dissipative elements are held constant, are discussed.
Abstract: The analysis and stability properties of a special class of R-L-C network, where the energy storage elements vary periodically with time and the dissipative elements are held constant, is discussed. In the analysis, fast Fourier transform technique has been utilized to obtain rapidly an accurate response of the network when excited by any arbitrary periodic signal. In the stability results, it has been shown that the network's stability can be ensured in a very simple manner if the connection matrix of the network is in Jordan normal form.

01 Jan 1973
TL;DR: This thesis aims to develop the tree graph idea, presented by E. Brigham and R. Morrow in the article "The fast Fourier transform," into another computer program that will calculate both the forward and inverse Fourier transforms using nonsymmetrical periodic functions.
Abstract: Many computer programs have been developed to calculate the Fast Fourier Transform (FFT). It is the intention of this thesis to develop the tree graph idea, presented by E. 0. Brigham and R. E. Morrow in the article "The fast Fourier transform," into another computer program that will calculate both the forward and inverse Fourier transforms using nonsymmetrical periodic functions. Most of the information obtained was from periodicals with relatively little available in text books. General matrix equations and calculations presented in proofs were developed by this thesis. The tree graph idea was also expanded to a general form. First the matrix equations for the Discrete Fourier Transform (DFT) were developed. Starting with the Fourier transform pair for an aperiodic continuous function and the matrix equations for the DFT, steps were taken to show how these equations were tied together. Two methods for the FFT wereeevaluated and proven equal to the DFT matrix equations. The more efficient method was written into a computer program which was explained in detail. Finally examples using this computer program were presented. The forward FFT was approximated using input functions such as a rectangular wave, a step function, a sawtooth wave, and a constant. A test for accuracy was also presented using a rectangular wave input time function.

Journal ArticleDOI
TL;DR: In this paper, a spectral analysis procedure for model identification and parameter evaluation of dynamic systems is given, assuming the process can be represented as a linear system, and an algorithm is given for converting the frequency response to spectra from which the poles and zeros of the processes can be identified.

Journal ArticleDOI
TL;DR: Fast Fourier Transforms (FFT) are widely known as useful tools for the evaluation of spectral data as mentioned in this paper, and the applicability of FFT methods to crystallographic problems is discussed.
Abstract: Fast Fourier Transforms (FFT) are widely known as useful tools for the evaluation of spectral data. This article discusses the applicability of FFT methods to crystallographic problems. Formulae are derived which make it possible to use fast Fourier transforms for the general Fourier summation of crystallographic data in all space groups and for the computation of slant planes at arbitrary positions in the unit cell. As even moderate resolutions would produce arrays of data too large to fit within internal computer memory, they must be kept to an external storage device. The organization of data and the resulting time requirements are thoroughly discussed. An ALGOL 60 program has been developed with as many redundancies eliminated as possible. An example shows that for a problem of moderate size this algorithm is faster by an order of magnitude than those which have traditionally been used.

01 Jun 1973
TL;DR: The performance of a modified version of NASA's multivariate spectrum analysis program is discussed and improvements were documented and compared to the standard fast Fourier transform by Monte Carlo techniques.
Abstract: The performance of a modified version of NASA's multivariate spectrum analysis program is discussed. A multiple regression model was used to make the revisions. Performance improvements were documented and compared to the standard fast Fourier transform by Monte Carlo techniques.

01 Apr 1973
TL;DR: It is shown that spatial domain recursive implementation of filters is much faster than nonrecursive implementations using fast fourier transforms (FFT) for the cases of interest in this work.
Abstract: Some techniques are presented and the software documentation for the digital enhancement of radiographs. Both image handling and image processing operations are considered. The image handling operations dealt with are: (1) conversion of format of data from packed to unpacked and vice versa; (2) automatic extraction of image data arrays; (3) transposition and 90 deg rotations of large data arrays; (4) translation of data arrays for registration; and (5) reduction of the dimensions of data arrays by integral factors. Both the frequency and the spatial domain approaches are presented for the design and implementation of the image processing operation. It is shown that spatial domain recursive implementation of filters is much faster than nonrecursive implementations using fast fourier transforms (FFT) for the cases of interest in this work. The recursive implementation of a class of matched filters for enhancing image signal to noise ratio is described. Test patterns are used to illustrate the filtering operations. The application of the techniques to radiographic images of metallic structures is demonstrated through several examples.