Abstract:
A method for fast Fourier transforming data in large quantities stored in an external storage at a high speed is disclosed. In the method, data stored in the external storage is supposed that the data consists of data elements of N1 rows and N2 columns. And the data is read out from the external storage along the column direction, transferred to the processor, fast Fourier transformed, and returned to the external storage. Then, the data is read out from the external storage along the row direction, transferred to the processor, fast Fourier transformed, and returned to the external storage. Bit reversal of data elements and transposition are done in the appropriate step of the procedure.

Description:
FIELD OF THE INVENTION 
     The present invention relates to a method operating an FFT processor in communication with an external storage device for fast Fourier transforming data in large quantities at a high speed. 
     BACKGROUND OF THE INVENTION 
     When data is fast Fourier transformed using a computer, complex coefficients, known as twiddle factors, are used. A twiddle factor is given by 
     
         W.sub.n.sup.k =cos (2π/n)k-sin (2π/n)k 
    
     That is, it is composed of the real part of a cosine function and the imaginary part of a sine function. In the above equation, n is the number of data elements to be fast Fourier transformed, k=0, . . . , n-1. FIG. 1 illustrates the algorithm of an FFT performed on 16 points. A butterfly operation is a unit operation of an FFT. The butterfly operation is conducted from left to right in stage 1 of FIG. 1. Then, stages 2 and 3 are handled. Correspondingly, the twiddle factor is quoted in the same order. 
     An FFT hardware processor may not handle whole data to be Fourier transformed in one FFT processing. In this context, Fourier transform processing involves a fast Fourier transform and transposition of data, or bit reversal. One reason why whole data cannot be Fourier transformed in one FFT processing is that the data cannot be placed into the memory of the processor simultaneously. 
     In order to take the Fourier transform of such a large quantity of data, it is necessary to provide access to the data in a complicated manner. 
     FIG. 2 shows the conventional method by which a memory is used to take the Fourier transform of data in large quantities. FIG. 2 also shows the processing. Data exists outside an FFT hardware processor. Where the processor is a host computer, data is stored in an external storage such as a magnetic disk. Where the processor is different from the host computer, data is stored in the memory of the host computer or in an external storage such as a magnetic disk. 
     If the capacity of the memory of the FFT processor is small, it cannot implement an FFT in one operation. In this case, the processor fetches as much data as possible to the memory as shown in FIG. 2 and performs butterfly operations on the data. The FFT implementation is divided into N (=log 2  n) stages. When each stage is separated into M parts, M times of transfers (fetch and store) of data and M times of butterfly operations must be executed. 
     When the processing is complex as shown in FIG. 2, fast Fourier transform (FFT) cannot be implemented at a high speed. Furthermore, the twiddle factor table must be extended up to a requisite size. Hence, the cost performance deteriorates. 
     In the case of FIG. 2, the data is transferred from the external storage to the processor and returned to the external storage from the processor, at each stage. Therefore, the total number of data transferred during the whole FFT implementation is given by 
     
         2×total amount of data×N stages 
    
     For this reason, a reduction in the speed cannot be avoided even if double buffering techniques are exploited. 
     Where data exists in a host computer, bit reversal may be done by the computer and so little problems take place. However, where data exists on an external storage such as a magnetic disk, the speed is reduced considerably for the following reason. Two data elements to be interchanged are placed on a position (x) and on a bit reversal position of the position (x), therefore, some data element pairs cannot be exchanged at a time. That is, one data element must be moved when one transfer instruction is issued. This reduces the transfer speed down to a minimum, where a minimum amount of data, e.g., Y data elements, must be transferred from an external storage such as magnetic disk, a wasteful transfer given by (Y-1)×4 takes place to interchange the elements of one pair. Since this transfer is conducted for all the data, all the wasteful transfer is given by 
     
         (Y-1)×4×(n/2) 
    
     In this way, wasteful transfer slows down the processing. 
     As described above, the method illustrated in FIG. 2 enables Fourier transform of data in large quantities, but the processing speed is quite low. 
     Other conventional methods for FFT of data in large quantities are described in an article entitled &#34;Fast Fourier Transform of Externally Stored Data&#34; by Norman M. Brenner in IEEE Transactions on Audio and Electroacoustics, Vol. AU-17, No. 2, pp. 128-132, June 1969. 
     In the article, two methods for fast Fourier transform of data in large quantities, named Ryder&#39;s method and Granger&#39;s method, are disclosed. Both methods have the following problems regarding processing speed. 
     Ryder&#39;s method contains twice reading accesses from the external storage and twice writing accesses to the external storage. When data to be fast Fourier transformed is supposed to consist of data elements of C1 rows and C2 columns, the number of data transferred during whole FFT implementation is estimated as follows: ##EQU1## Ordinarily, since C1&gt;&gt;1, the number of data read from the external storage in the first reading access is quite large. Therefore, the number of data transferred during whole FFT implementation is quite large. Then, it cannot be avoided to lower the processing speed. 
     On the other hand, Granger&#39;s method contains five times of reading accesses from the external storage and five times of writing accesses to the external storage. In the method, the number of data transferred during whole FFT implementation is estimated to 10×n. It can be understood that 10 times of reading/writing accesses cause to lower the processing speed. 
     SUMMARY OF THE INVENTION 
     Accordingly, it is an object of the present invention to provide a process of operating a FFT hardware processor in communication with an external data storage device for fast Fourier transform of data in large quantities stored in said external storage device at a high speed. 
     According to the present invention, data stored in the external storage is considered to consist of data elements organized in N1 rows and N2 columns such that the data in external storage may be considered to consist of data elements organized in rows and columns defined by index bits. The data is read out from the external storage device along the column or row direction, transferred to the FFT processor, fast Fourier transformed, and returned to the external storage. Then, the data is again read out from the external storage along the row or column direction, transferred to the processor, fast Fourier transformed, and returned to the external storage. Index bit reversal of data elements and transposition are performed in the appropriate step(s) of the procedure. 
     According to a preferred method, the data read from the external storage device along the column direction is transferred to the processor. It is fast Fourier transformed by the processor. Index bits of data elements are reversed in the processor. The data is returned to the external storage. Next, the data is read from the external storage device along the row direction and transferred to the processor. It is multiplied by compensating coefficients and fast Fourier transformed by the processor. Index bits of data elements are reversed in the processor. The data is returned to the external storage. Finally, the data stored in the external storage is transposed. 
    
    
     In a further preferred method, only part of the data is transferred to the FFT processor at one time. The steps following transfer of data are the same and are repeated until all the data stored in the external storage is transformed and returned. 
     BRIEF DESCRIPTION OF THE DRAWINGS 
     FIG. 1 is a diagram illustrating an algorithm by which FFT operations are performed to handle 16 data elements; 
     FIG. 2 is a diagram illustrating the usage of a memory and the processing in the prior art Fourier transform handling data in large quantities; 
     FIG. 3 is a diagram illustrating an algorithm by which data is expanded in two dimensions; 
     FIG. 4 is a flowchart illustrating an FFT method of handling data in large quantities; 
     FIG. 5 is a conceptual view illustrating the use of a memory and the transfer of data in the first step of an FFT method according to the invention: 
     FIG. 6 is a conceptual view illustrating the use of a memory and the transfer of data in the second step of an FFT method according to the invention; 
     FIG. 7 is a diagram illustrating the correspondence between compensating coefficients and data elements; 
     FIG. 8 is a diagram illustrating the manner in which data is moved when data expanded in two dimensions is transposed; 
     FIGS. 9a-b show the relation between the arrangement of 16 twiddle factors and the arrangement of 4 rotation factors; 
     FIG. 10 is a diagram illustrating an algorithm by which butterfly operations are performed to handle 16 data elements as 4×4 data elements, i.e., in two dimensions; 
     FIG. 11 is a diagram illustrating the manner in which vectors are created; and 
     FIGS. 12a-c are diagrams of an FFT hardware system handling data in large quantities, the system being built in accordance with the invention. 
    
    
     DETAILED DESCRIPTION OF THE INVENTION 
     The present invention is applied to a combination of any type of FFT hardware processor, e.g., a general-purpose computer, array processor, or the like, and an associated external storage device such as a magnetic disk or the memory of a host computer. The present invention is intended to reduce the total number of data transfers to a minimum for preventing wasteful transfer of data. To realize this, data elements arranged in a line are considered rearranged in rows and columns and then the data is so treated. This method is referred to as the two-dimensional expansion method herein. 
     Referring to FIG. 3, one-dimensional data consisting of n data elements arranged in a line is stored in an external storage. We now assume that a processor is available by which up to m (n&gt;m) data elements can be Fourier transformed. As shown in FIG. 3, the n data elements are interpreted as consisting of N1 rows and N2 columns, i.e., n=N1×N2, where N1≦m and N2≦m. Note that this is merely an interpretation. Therefore, the actual addresses at which the data elements exist on the external storage need not be modified. 
     On this assumption, data is fast Fourier transformed. The number of stages of the butterfly operations is given by 
     
         log.sub.2 n=log.sub.2 N1+log.sub.2 N2 
    
     Intrinsically, all the stages equal to log 2  n should be fast Fourier transformed in one FFT processing. However, the present invention is characterized in that all the stages to be fast Fourier transformed in one FFT processing are divided into stages equal to log 2  N1 and stages equal to log 2  N2. They are referred to as the first step FFT implementation and the second step FFT implementation, respectively. 
     The flow of the FFT implementation is illustrated in FIG. 4. The flow is divided into three processing stages. The first processing stage is the first step FFT implementation. The second processing stage comprises compensation of data and the second step FFT implementation. The third processing stage is to transpose the data expanded in two dimensions, for accomplishing the index bit reversal of the data. 
     In the process illustrated in FIG. 4, N1 data elements for the first step are fast Fourier transformed, and their bits are reversed. N2 data elements for the second step are fast Fourier transformed, and they are bit reversed. Finally, data expanded in two dimensions is transposed to accomplish the bit reversal of the data. 
     For bit reversal of all the data elements, the following two modifications are available. 
     (1) The data about the first step is fast Fourier transformed without bit reversal. The data is fast Fourier transformed in the second step without bit reversal. Then, the data is bit reversed longitudinally and transversely. The data expanded in two dimensions is transposed. 
     (2) Where the memory of a host computer is used as an external storage, no bit reversal is executed at the time of FFT implementation. Finally, the bit reversal of the data is executed by the host computer. 
     FIG. 5 conceptually illustrates the usage of the memory and the transfer of data in the first step FFT implementation. 
     In the first step FFT implementation, data is subjected to cycle (=N2/x) division as indicated by the hatching in FIG. 5. Each divided portion is treated as follows. A set of data S1 of the first column consisting of N1 data elements and stored in an external storage is transferred to the memory of the FFT processor. The remaining sets of data (S2, S3, . . . Sx) are subsequently transferred in turn. Thus, one divided portion is placed into the memory of the FFT processor. If the plural sets of data spaced from each other by N2 data elements in the external storage are made continuous with each other in the memory, then FFT of the divided portion can be easily implemented. After FFT, x bit reversals are executed for each N1 data elements in the FFT processor are reversed. This bit reversal is carried out in the same memory of the FFT processor. The divided portion of the data which has been fast Fourier transformed and bit reversed is returned to its original position or to a location corresponding to another data region in the external storage. The above-described procedure is repeated N2/x times until the whole divided portions are fast Fourier transformed and bit reversed. 
     FIG. 6 conceptually illustrates the usage of a memory and the transfer of data in the second step FFT implementation. FIG. 7 shows the correspondence between compensating factors and data elements. 
     In the second step FFT implementation, data stored in the external storage after the first step FFT is subjected to cycle (=N1/y) division as indicated by the hatching in FIG. 6. Each divided portion is treated in the manner described below. First divided portion consisting of a succession of data, i.e., y by N2 data elements is transferred into the memory of the FFT processor from the external storage. Y requisite compensating coefficients are created for every set of N2 data elements and stored in the coefficient memory. Then, the respective N2 data elements constituting the first divided portion are multiplied by the corresponding coefficients and fast Fourier transformed. The bit reversal of the transformed data of the divided portion is executed. The fully processed data of the divided portion is returned to its original positions in the external storage or to locations corresponding to a different data region. 
     FIG. 7 shows by what compensating coefficients are data elements multiplied prior to the second step FFT implementation. The compensating coefficients are used for the data which has been bit reversed in the first step. If the bit reversal is not executed in the first step, the order is different. One compensating coefficient is needed for each one data element. The compensating coefficients shown in FIG. 7 do not exist at all times, but they are created as the need arises as described already. In FIG. 7, W n   i  · j  means ##EQU2## In the case of inverse Fourier transform, W n   i  · j  takes the conjugate complex number. 
     FIG. 8 shows the manner in which data is shifted when data expanded in two dimensions is transposed. When N2=N1, data can be transposed inside the same data region. When N2≠N1, data cannot be transposed quickly unless data is moved into a different data region. Where the external storage for the FFT processor is a magnetic disk or the like, data is once transferred into the memory of the host computer which plays key roles in the processing and then the data is returned to the external storage. At this time, the data expanded in two dimensions can be transposed, making use of the memory of the computer. If the data cannot fully enter the host computer, the data is divided into plural portions and they are separately handled, in the same way as in the first and second steps of FFT implementation. If the external storage is the memory of the host computer, data may be transposed within the memory. 
     Data has been interpreted as being expanded in two dimensions. Here the data is reinterpreted as one-dimensional data. Thus, the fast Fourier transform of data in large quantities is completed. Obviously, it is not necessary to modify the addresses of the data. 
     We next discuss the case in which 16 data elements are fast Fourier transformed, using a processor capable of fast Fourier transforming only 4 data elements. 
     FIG. 9 shows the relation between the arrangement of 16 twiddle factors and the arrangement of 4 twiddle factors. FIG. 10 is a diagram illustrating the algorithm by which data elements are treated as 4×4 data elements, i.e., as two-dimensional data, in butterfly operations. FIG. 11 is a diagram illustrating how to create compensating coefficients vector. The algorithm permitting 16 data elements to be fast Fourier transformed by the use of an FFT processor capable of fast Fourier transforming only 4 data elements depends much on how efficiently is FFT reduced to discrete Fourier transform (DFT). In principle, 16 data elements are treated as two-dimensional data composed of x×y=16. Fast Fourier transform of 4 data elements is implemented in the x-direction and in the y-direction separately. The obtained data is coupled by DFT to implement Fourier transform of 16 data elements. 
     We now notice the relation between the 16 twiddle factors W N   k  and the 4 twiddle factors W N   k . As can be seen from FIG. 9, 
     
         W.sub.16.sup.0 =W.sub.4.sup.0 
    
     
         W.sub.16.sup.4 =W.sub.4.sup.1 
    
     
         W.sub.16.sup.8 =W.sub.4.sup.2 
    
     
         W.sub.16.sup.12 =W.sub.4.sup.3 
    
     The butterfly operation in which 16 data elements are treated as 4×4 two-dimensional data elements, is illustrated in FIG. 10. Comparing the butterfly operation on the 16 data elements as shown in FIG. 10 with the butterfly operation shown in FIG. 1 reveals that the results E&#39;s of the calculations in the portion a (FFT in the y-direction) are exactly the same, since W 16   4  =W 4   1  is held. Therefore, we should take notice of portion b and can determine the compensating coefficients vector R of DFT acting on the E&#39;s. 
     Noting a of portion b, we resolve the butterfly. From the 16 data elements, we have ##EQU3## From 4×4 data elements, we get ##EQU4## Comparing equations (1) with equations (1)&#39; while considering the relation between W 16   k  and W 4   k , we see that they are exactly the same. Therefore, we have the following relationships: 
     
         E.sub.0 =E.sub.0 &#39;, E.sub.4 =E.sub.1 &#39; 
    
     
         E.sub.8 =E.sub.2 &#39;, E.sub.12 =E.sub.3 &#39; 
    
     Noting β of portion b, we resolve the butterflies. From the 16 data elements, we have ##EQU5## From 4×4 data elements, we have ##EQU6## 
     We compare equations (2) with equations (2)&#39;, considering the relation between W 16   k  and W 4   k . In order that the results of equations (2) be equal to the results of equations (2)&#39;, the following relations must hold: ##EQU7## Similarly, we respect to γ and δ of portion b, the following relations must be satisfied: ##EQU8## 
     In this way, the vector R is given by (1, 1, 1, 1, 1, W 16   2 , W 16   4 , W 16   6 , 1, W 16   1 , W 16   2 , W 16   3 , 1, W 16   3 , W 16   6 , W 16   9 ). 
     This vector is obtained by replacing W 4   k  of the matrix of the twiddle factors acting during DFT of four data elements with W 16   k , causing the bit reversal of 4 (2 bits) to act on the row of the matrix shown in FIG. 11. The use of this method enables a large amount of data to be Fourier transformed. 
     FIG. 12 shows examples of the fast Fourier transform hardware system according to the invention. FIG. 12(a) shows an example in which an FFT processor 1 is combined with an external storage 2 storing data. FIG. 12(b) shows an example in which a host computer 3 and an auxiliary processor 4 together form an FFT processor. This processor is combined with an external storage 5 in which data is stored. FIG. 12(c) shows an example in which data is stored in the memory of a host computer 6 and in an external storage. An auxiliary processor 7 other than the host computer 6 is employed as an FFT processor. 
     As can be seen from the description made thus far, in accordance with the present invention, the amount of transferred data including the amount of data transferred for FFT and bit reversal is given by total number of data transferred=2×total amount of data×3.  Therefore, the total number of data transferred is much smaller than that of the prior art method described in connection with FIG. 2. Hence, the processing can be carried out at a high speed. If double buffering techniques are employed in the processing, the effect of this processing method can be enhanced further. Additionally, slowdown of Fourier transform of data in large quantities can be prevented even if the capacity of the FFT hardware processor is insufficient. Since a large amount of data can be Fourier transformed at a high speed without adding any memory to the FFT processor, the system is excellent in cost performance. 
     Having thus described our invention with the detail and particularity required by the Patent Laws, what is claimed and desired to be protected by Letters Patent is set forth in the following claims.