CA2725894C - Sine/cosine generator - Google Patents
Sine/cosine generator Download PDFInfo
- Publication number
- CA2725894C CA2725894C CA2725894A CA2725894A CA2725894C CA 2725894 C CA2725894 C CA 2725894C CA 2725894 A CA2725894 A CA 2725894A CA 2725894 A CA2725894 A CA 2725894A CA 2725894 C CA2725894 C CA 2725894C
- Authority
- CA
- Canada
- Prior art keywords
- sine
- angle
- cosine
- approximation
- coarse
- 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 - Fee Related
Links
- 238000000034 method Methods 0.000 abstract description 84
- 230000006870 function Effects 0.000 description 17
- 238000000354 decomposition reaction Methods 0.000 description 12
- 101100095718 Homo sapiens SH2D2A gene Proteins 0.000 description 6
- 102100031779 SH2 domain-containing protein 2A Human genes 0.000 description 6
- 238000006243 chemical reaction Methods 0.000 description 6
- 238000011156 evaluation Methods 0.000 description 6
- 230000014509 gene expression Effects 0.000 description 6
- 238000013507 mapping Methods 0.000 description 6
- 238000012545 processing Methods 0.000 description 6
- 238000000638 solvent extraction Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000004590 computer program Methods 0.000 description 4
- 238000007796 conventional method Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000007792 addition Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000004883 computer application Methods 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F7/00—Methods or arrangements for processing data by operating upon the order or content of the data handled
- G06F7/38—Methods or arrangements for performing computations using exclusively denominational number representation, e.g. using binary, ternary, decimal representation
- G06F7/48—Methods or arrangements for performing computations using exclusively denominational number representation, e.g. using binary, ternary, decimal representation using non-contact-making devices, e.g. tube, solid state device; using unspecified devices
- G06F7/544—Methods or arrangements for performing computations using exclusively denominational number representation, e.g. using binary, ternary, decimal representation using non-contact-making devices, e.g. tube, solid state device; using unspecified devices for evaluating functions by calculation
- G06F7/548—Trigonometric functions; Co-ordinate transformations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F1/00—Details not covered by groups G06F3/00 - G06F13/00 and G06F21/00
- G06F1/02—Digital function generators
- G06F1/03—Digital function generators working, at least partly, by table look-up
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F7/00—Methods or arrangements for processing data by operating upon the order or content of the data handled
- G06F7/38—Methods or arrangements for performing computations using exclusively denominational number representation, e.g. using binary, ternary, decimal representation
- G06F7/48—Methods or arrangements for performing computations using exclusively denominational number representation, e.g. using binary, ternary, decimal representation using non-contact-making devices, e.g. tube, solid state device; using unspecified devices
- G06F7/499—Denomination or exception handling, e.g. rounding or overflow
- G06F7/49942—Significance control
Abstract
A method (300) for generating a sine and cosine of an input angle (.SLZERO.102). The method involves decomposing .SLZERO.102 to an octant or quadrant, a coarse angle (A), and a fine angle (B), determining cos(A), and determining sin(A). The method also involves decomposing cos(A) and sin(A) to a most significant word (MSW) and a least significant word (LSW). The method further involves computing an approximation of 1-cos(B), an approximation of sin(B), and a plurality of products (P1,..., P4) using the MSWs and approximations. The method involves computing approximations of cos(.SLZERO.'102) and sin(.SLZERO.'102) using the values for cos(A), sin(A), and P1,..., P4. The method involves scaling the approximations of cos(.SLZERO.'102) and sin(.SLZERO.'102) to a desired resolution.
Description
SINE/COSINE GENERATOR
The invention concerns communications systems, and more particularly, the efficient digital generation of sine/cosine evaluations.
Conventional sine/cosine generators typically compute a sine and cosine of an input angle. According to one communications system application, the sine and cosine of the same input angle can be used together to form a quadrature form of a sinusoid, which is the equivalent to the complex-valued exponential evaluated at the input angle. Quadrature forms of sinusoids are well known to those having ordinary skill in the art, and therefore will not be described herein.
However, it should be understood that quadrature sinusoids are often used in digital up conversion applications, digital down conversion applications, numerically controlled oscillator applications, and Fourier transform applications.
One conventional method for generating a quadrature sinusoid employs look up tables (LUT). The LUT based method is well known to those having ordinary skill in the art, and therefore will not be described herein.
However, it should be understood that the LUT based method generally involves mapping a fixed-precision input angle to a pre-defined evaluation of a trigonometric function within some resulting accuracy. The LUT based method generally provides reasonably efficient computation for ten to twelve (10-12) bit precision sine and cosine trigonometric computations. If precisions greater than ten to twelve (10-12) bits is desired, then the LUT based method requires exponentially increasing amounts of hardware for linear increases in bits of precision. Therefore, the LUT based method is hardware inefficient for high-accuracy applications.
Other conventional methods for generating a quadrature sinusoid employ Coordinate Rotation Digital Computer (CORDIC) algorithms. CORDIC
algorithms generally use vector rotation to compute the sine and cosine of an input angle. CORDIC algorithms involve iteratively performing vector rotations by an arbitrary angle using shift and add techniques. As such, CORDIC based methods are computationally intensive. Further, the basic CORDIC method cannot take advantage of any natural symmetries in the underlying trigonometric operators. As such, the basic CORDIC method is less efficient than tailored hardware calculators.
Another conventional method of generating a quadrature sinusoid employs an out-of-phase algorithm. This out-of-phase based method is described in U.S. Patent No. 5,276,633. This out-of-phase based method more efficiently computes the sine and cosine of an input angle as compared to the LUT and CORDIC
based methods. In this regard, it should be understood that the out-of-phase based method provides ten to eighteen (10-18) bit precision sine and cosine trigonometric computations. However, such a ten to eighteen (10-18) bit precision is unsuitable for digital signal processing applications requiring a higher degree of accuracy.
In view of the forgoing, there is a need for an improved sine/cosine generator and method for efficient digital generation of a sine/cosine of an input angle. The improved method needs to be less computationally intensive than the conventional CORDIC based methods. The improved method also needs to provide a bit precision greater than the bit precisions of conventional LUT and out-of-phase based methods.
This Summary is provided to comply with 37 C.F.R. 1.73, requiring a summary of the invention briefly indicating the nature and substance of the invention.
It is submitted with the understanding that it will not be used to interpret or limit the scope or meaning of the claims.
The present invention concerns a method for efficient digital generation a sine and cosine of an input angle. The method involves decomposing the input angle into an octant or quadrant, a coarse angle, and a fine angle. The coarse angle can be in a quadrant, an octant, or a finer-resolution decomposition of a notional two dimensional (2D) complex-valued space. The fine angle can be bounded by a smallest coarse angle existing in the 2D space. The method also involves determining a cosine and a sine of the coarse angle. The cosine of the coarse angle is decomposed to a first most significant word and a first least significant word.
Similarly, the sine of the coarse angle is decomposed to a second most significant word and a second least significant word. The method also involves computing an approximation of one
The invention concerns communications systems, and more particularly, the efficient digital generation of sine/cosine evaluations.
Conventional sine/cosine generators typically compute a sine and cosine of an input angle. According to one communications system application, the sine and cosine of the same input angle can be used together to form a quadrature form of a sinusoid, which is the equivalent to the complex-valued exponential evaluated at the input angle. Quadrature forms of sinusoids are well known to those having ordinary skill in the art, and therefore will not be described herein.
However, it should be understood that quadrature sinusoids are often used in digital up conversion applications, digital down conversion applications, numerically controlled oscillator applications, and Fourier transform applications.
One conventional method for generating a quadrature sinusoid employs look up tables (LUT). The LUT based method is well known to those having ordinary skill in the art, and therefore will not be described herein.
However, it should be understood that the LUT based method generally involves mapping a fixed-precision input angle to a pre-defined evaluation of a trigonometric function within some resulting accuracy. The LUT based method generally provides reasonably efficient computation for ten to twelve (10-12) bit precision sine and cosine trigonometric computations. If precisions greater than ten to twelve (10-12) bits is desired, then the LUT based method requires exponentially increasing amounts of hardware for linear increases in bits of precision. Therefore, the LUT based method is hardware inefficient for high-accuracy applications.
Other conventional methods for generating a quadrature sinusoid employ Coordinate Rotation Digital Computer (CORDIC) algorithms. CORDIC
algorithms generally use vector rotation to compute the sine and cosine of an input angle. CORDIC algorithms involve iteratively performing vector rotations by an arbitrary angle using shift and add techniques. As such, CORDIC based methods are computationally intensive. Further, the basic CORDIC method cannot take advantage of any natural symmetries in the underlying trigonometric operators. As such, the basic CORDIC method is less efficient than tailored hardware calculators.
Another conventional method of generating a quadrature sinusoid employs an out-of-phase algorithm. This out-of-phase based method is described in U.S. Patent No. 5,276,633. This out-of-phase based method more efficiently computes the sine and cosine of an input angle as compared to the LUT and CORDIC
based methods. In this regard, it should be understood that the out-of-phase based method provides ten to eighteen (10-18) bit precision sine and cosine trigonometric computations. However, such a ten to eighteen (10-18) bit precision is unsuitable for digital signal processing applications requiring a higher degree of accuracy.
In view of the forgoing, there is a need for an improved sine/cosine generator and method for efficient digital generation of a sine/cosine of an input angle. The improved method needs to be less computationally intensive than the conventional CORDIC based methods. The improved method also needs to provide a bit precision greater than the bit precisions of conventional LUT and out-of-phase based methods.
This Summary is provided to comply with 37 C.F.R. 1.73, requiring a summary of the invention briefly indicating the nature and substance of the invention.
It is submitted with the understanding that it will not be used to interpret or limit the scope or meaning of the claims.
The present invention concerns a method for efficient digital generation a sine and cosine of an input angle. The method involves decomposing the input angle into an octant or quadrant, a coarse angle, and a fine angle. The coarse angle can be in a quadrant, an octant, or a finer-resolution decomposition of a notional two dimensional (2D) complex-valued space. The fine angle can be bounded by a smallest coarse angle existing in the 2D space. The method also involves determining a cosine and a sine of the coarse angle. The cosine of the coarse angle is decomposed to a first most significant word and a first least significant word.
Similarly, the sine of the coarse angle is decomposed to a second most significant word and a second least significant word. The method also involves computing an approximation of one
-2-minus a cosine of the fine angle and an approximation of a sine of the fine angle. The method further involves computing a plurality of products, generating an approximation of a cosine of the input angle modulo it/k, and generating an approximation of a sine of the input angle modulo it/k using the cosine of the coarse angle, the sine of the coarse angle, and each of the products. k is an integer value selected in accordance with a particular type of partitioning and reconstruction process employed. For example, if an octant partitioning and reconstruction process is employed, then k is selected to be an integer equal to four (4). If a quadrant partitioning and reconstruction process is employed, then k is selected to be an integer equal to two (2).
According to an aspect of the invention, a first product is computed by multiplying the first most significant word by the approximation of a sine of the fine angle. A second product is computed by multiplying the first most significant word by the approximation of one minus a cosine of the fine angle. A third product is computed by multiplying the second most significant word by the approximation of one minus a cosine of the fine angle. A fourth product is computed by multiplying the second most significant word by the approximation of a sine of the fine angle.
According to another aspect of the invention, the approximation of the cosine of the input angle modulo it/k is generated by computing a first difference and subtracting the fourth product from the first difference. The first difference is computed by subtracting the second product from the cosine of the coarse angle. The approximation of the sine of the input angle modulo it/k is generated by computing a second difference and combining the first product with the second difference.
The second difference is computed by subtracting the third product from the sine of the coarse angle.
According to another aspect of the invention, the method involves scaling the approximation of the sine of the input angle modulo it/k by a first scale factor. The scaled approximation of the sine of the input angle modulo it/k is combined with the approximation of the cosine of the input angle modulo it/k to yield the cosine of the input angle modulo it/k. The method also involves scaling the
According to an aspect of the invention, a first product is computed by multiplying the first most significant word by the approximation of a sine of the fine angle. A second product is computed by multiplying the first most significant word by the approximation of one minus a cosine of the fine angle. A third product is computed by multiplying the second most significant word by the approximation of one minus a cosine of the fine angle. A fourth product is computed by multiplying the second most significant word by the approximation of a sine of the fine angle.
According to another aspect of the invention, the approximation of the cosine of the input angle modulo it/k is generated by computing a first difference and subtracting the fourth product from the first difference. The first difference is computed by subtracting the second product from the cosine of the coarse angle. The approximation of the sine of the input angle modulo it/k is generated by computing a second difference and combining the first product with the second difference.
The second difference is computed by subtracting the third product from the sine of the coarse angle.
According to another aspect of the invention, the method involves scaling the approximation of the sine of the input angle modulo it/k by a first scale factor. The scaled approximation of the sine of the input angle modulo it/k is combined with the approximation of the cosine of the input angle modulo it/k to yield the cosine of the input angle modulo it/k. The method also involves scaling the
-3-
4 PCT/US2009/044653 approximation of the cosine of the input angle modulo it/k by a second scale factor.
The scaled approximation of the cosine of the input angle modulo 7L/k is combined with the approximation of the sine of the input angle modulo 7L/k to yield the sine of the input angle modulo it/k. The cosine of the input angle modulo it/k is mapped to a cosine of the input angle using N most significant bits from the input angle.
Similarly, the sine of the input angle it/k is mapped to a sine of the input angle using the N most significant bits from the input angle. N is an integer value selected in accordance with a particular type of partitioning and reconstruction process employed. For example, if an octant partitioning and reconstruction process is employed, then N is selected to be an integer equal to three (3). If a quadrant partitioning and reconstruction process is employed, then N is selected to be an integer equal to two (2).
A sine/cosine generator is also provided. The sine/cosine generator is comprised of a decomposition device, at least one coarse angle device, at least one fine angle device, a product computation device, and an output computation device.
The decomposition device is configured for decomposing the input angle into an octant or quadrant, a coarse angle, and a fine angle. Those with ordinary skill in the art will recognize that this decomposition technique and resulting processing may be extended to a coarse angle, a fine angle, and a hyper-fine angle using similar trigonometry identities. The coarse angle device is configured for determining a cosine of the coarse angle, determining a sine of the coarse angle, decomposing the cosine of the coarse angle to a first most significant word and a first least significant word, and decomposing the sine of the coarse angle to a second most significant word and a second least significant word. The fine angle device is configured for computing an approximation of one minus a cosine of the fine angle and an approximation of a sine of the fine angle. The product computation device is configured for computing a plurality of products.
The output computation device is coupled to the product computation device and the coarse angle device. The output computation device is configured for generating approximations of a cosine of the input angle modulo it/k and a sine of the input angle modulo it/k using the cosine of the coarse angle, the sine of the coarse angle, and each of the products. The output computation device is also configured for computing the cosine and sine of the input angle modulo it/k using scale factors, the approximation of the cosine of the input angle modulo it/k, and the approximation of the sine of the input angle modulo it/k. The output computation device is further configured for mapping the cosine and sine of the input angle modulo 7L/k to a cosine and sine of the input angle using the N most significant bits of the input angle.
Embodiments will be described with reference to the following drawing figures, in which like numerals represent like items throughout the figures, and in which:
FIG. 1 is a conceptual diagram of a process for generating a sine and cosine of an input angle that is useful for understanding the present invention.
FIG. 2 is a two-dimensional graph illustrating quadrants and octants of a two dimensional space that is useful for understanding the present invention.
FIGS. 3A-3B collectively provide a flow diagram of a method for generating a sine and cosine of an input angle that is useful for understanding the present invention.
FIG. 4 is a block diagram of an exemplary sine/cosine generator that is useful for understanding the present invention.
FIG. 5 is a block diagram of another exemplary sine/cosine generator that is useful for understanding the present invention.
The present invention concerns efficient, fully pipelined non-linear processor (NLP) based sine/cosine generators. The NLP based sine/cosine generator provides an improved output bit precision as compared to conventional sine/cosine generators. For example, the NLP based sine/cosine generator provides greater than eighteen (18) bits of output precision. In contrast, conventional sine/cosine generators often provide ten to eighteen (10-18) bits of output precision.
The invention will now be described more fully hereinafter with reference to accompanying drawings, in which illustrative embodiments of the invention are shown. This invention, may however, be embodied in many different
The scaled approximation of the cosine of the input angle modulo 7L/k is combined with the approximation of the sine of the input angle modulo 7L/k to yield the sine of the input angle modulo it/k. The cosine of the input angle modulo it/k is mapped to a cosine of the input angle using N most significant bits from the input angle.
Similarly, the sine of the input angle it/k is mapped to a sine of the input angle using the N most significant bits from the input angle. N is an integer value selected in accordance with a particular type of partitioning and reconstruction process employed. For example, if an octant partitioning and reconstruction process is employed, then N is selected to be an integer equal to three (3). If a quadrant partitioning and reconstruction process is employed, then N is selected to be an integer equal to two (2).
A sine/cosine generator is also provided. The sine/cosine generator is comprised of a decomposition device, at least one coarse angle device, at least one fine angle device, a product computation device, and an output computation device.
The decomposition device is configured for decomposing the input angle into an octant or quadrant, a coarse angle, and a fine angle. Those with ordinary skill in the art will recognize that this decomposition technique and resulting processing may be extended to a coarse angle, a fine angle, and a hyper-fine angle using similar trigonometry identities. The coarse angle device is configured for determining a cosine of the coarse angle, determining a sine of the coarse angle, decomposing the cosine of the coarse angle to a first most significant word and a first least significant word, and decomposing the sine of the coarse angle to a second most significant word and a second least significant word. The fine angle device is configured for computing an approximation of one minus a cosine of the fine angle and an approximation of a sine of the fine angle. The product computation device is configured for computing a plurality of products.
The output computation device is coupled to the product computation device and the coarse angle device. The output computation device is configured for generating approximations of a cosine of the input angle modulo it/k and a sine of the input angle modulo it/k using the cosine of the coarse angle, the sine of the coarse angle, and each of the products. The output computation device is also configured for computing the cosine and sine of the input angle modulo it/k using scale factors, the approximation of the cosine of the input angle modulo it/k, and the approximation of the sine of the input angle modulo it/k. The output computation device is further configured for mapping the cosine and sine of the input angle modulo 7L/k to a cosine and sine of the input angle using the N most significant bits of the input angle.
Embodiments will be described with reference to the following drawing figures, in which like numerals represent like items throughout the figures, and in which:
FIG. 1 is a conceptual diagram of a process for generating a sine and cosine of an input angle that is useful for understanding the present invention.
FIG. 2 is a two-dimensional graph illustrating quadrants and octants of a two dimensional space that is useful for understanding the present invention.
FIGS. 3A-3B collectively provide a flow diagram of a method for generating a sine and cosine of an input angle that is useful for understanding the present invention.
FIG. 4 is a block diagram of an exemplary sine/cosine generator that is useful for understanding the present invention.
FIG. 5 is a block diagram of another exemplary sine/cosine generator that is useful for understanding the present invention.
The present invention concerns efficient, fully pipelined non-linear processor (NLP) based sine/cosine generators. The NLP based sine/cosine generator provides an improved output bit precision as compared to conventional sine/cosine generators. For example, the NLP based sine/cosine generator provides greater than eighteen (18) bits of output precision. In contrast, conventional sine/cosine generators often provide ten to eighteen (10-18) bits of output precision.
The invention will now be described more fully hereinafter with reference to accompanying drawings, in which illustrative embodiments of the invention are shown. This invention, may however, be embodied in many different
-5-forms and should not be construed as limited to the embodiments set forth herein. For example, the present invention can be embodied as a method, a data processing system or a computer program product. Accordingly, the present invention can take the form as an entirely hardware embodiment, an entirely software embodiment, or a hardware/software embodiment.
Referring now to FIG. 1, there is provided a conceptual diagram of a process 100 for generating a sine and cosine of an input angle that is useful for understanding the present invention. As shown in FIG. 1, the process 100 can begin by receiving an input angle 0102 from an external device (not shown). The external device can include, but is not limited to, an angle generator that generates angles between zero and three hundred sixty degrees (0 -360 ), non-inclusive of the three hundred sixtieth degree (360 ) using an input signal. Angle generators and sine/cosine generators are often used in a variety of communications system applications. Such communications system applications include, but are not limited to, digital up conversion applications, digital down conversion applications, and Fourier transform applications After receiving an input angle 0102, the sine and cosine generation process 100 continues with a first, second, and third determination step 103, 104, 106.
The first determination step 103 extracts the three (3) most significant bits (MSBs) of the binary represented angular input and uses them to compute the octant of the input angle 0102. In this regard, it should be understood that step 103 can be provided by performing any method known to those having ordinary skill in the art for decomposing an angle into an octant designation. Those with ordinary skill in the art will recognize that step 103 may be implemented using simple digital logic and contributes three (3) effective bits of output accuracy to the final result.
The outputs of step 103 are used to adjust the octant-agnostic sine/cosine evaluation calculated through the remainder of the sine and cosine generation process 100. The outputs of step 103 are also used to pre-adjust the fine angle evaluations for any symmetries and/or reflections applied.
Referring now to FIG. 1, there is provided a conceptual diagram of a process 100 for generating a sine and cosine of an input angle that is useful for understanding the present invention. As shown in FIG. 1, the process 100 can begin by receiving an input angle 0102 from an external device (not shown). The external device can include, but is not limited to, an angle generator that generates angles between zero and three hundred sixty degrees (0 -360 ), non-inclusive of the three hundred sixtieth degree (360 ) using an input signal. Angle generators and sine/cosine generators are often used in a variety of communications system applications. Such communications system applications include, but are not limited to, digital up conversion applications, digital down conversion applications, and Fourier transform applications After receiving an input angle 0102, the sine and cosine generation process 100 continues with a first, second, and third determination step 103, 104, 106.
The first determination step 103 extracts the three (3) most significant bits (MSBs) of the binary represented angular input and uses them to compute the octant of the input angle 0102. In this regard, it should be understood that step 103 can be provided by performing any method known to those having ordinary skill in the art for decomposing an angle into an octant designation. Those with ordinary skill in the art will recognize that step 103 may be implemented using simple digital logic and contributes three (3) effective bits of output accuracy to the final result.
The outputs of step 103 are used to adjust the octant-agnostic sine/cosine evaluation calculated through the remainder of the sine and cosine generation process 100. The outputs of step 103 are also used to pre-adjust the fine angle evaluations for any symmetries and/or reflections applied.
-6-It should be noted that the invention is not limited to octant symmetry based applications as shown in FIG. 1. The present invention can also be used in quadrant symmetry based application. In such a scenario, the first determination step 103 can alternatively involve extracting the two (2) most significant bits (MSBs) of the binary represented angular input and use them to compute the quadrant of the input angle 0102.
Referring again to FIG. 1, step 104 can generally involve decomposing the input angle 0102 to a coarse angle A. In this regard, it should be understood that the second determination step 104 can be provided by performing any method known in the art for decomposing an angle into a coarse angle A. For example, the second determination step 104 can involve truncating or rounding a binary representation of the input angle 0102 after a fixed number of bits to obtain a coarse angle A.
In the case of octant symmetry (which is efficient in the concurrent computation of a cosine and sine), the coarse angle A is decomposed in the first octant exclusive of the three (3) most significant bits. If a rounding method is used, an adjustment (not shown) to the fine angle is required to recover the sign. In the case of quadrant symmetry, the coarse angle is decomposed in the first quadrant exclusive of the two (2) most significant bits. The invention is not limited in this regard.
According to an embodiment of the invention, the coarse angle A is an angle in a two dimensional (2D) space. For example, the coarse angle A is one of a plurality of angles in a first octant or quadrant. A graph 200 illustrating a 2D space divided into quadrants and octants is provided in FIG. 2. As shown in FIG. 2, each quadrant I, II, III, IV is comprised of a plurality of angles residing in a respective ninety degree (90 ) area which is non-inclusive of the ninetieth-degree (90 ).
The first quadrant I includes angles between zero degrees (0 ) and ninety degrees (90 ).
The second quadrant II includes angles between ninety degrees (90 ) and one hundred eight degrees (180 ). The third quadrant III includes angles between one hundred eight degrees (180 ) and two hundred seventy degrees (270 ). The fourth quadrant IV includes angles between two hundred seventy degrees (270 ) and three hundred sixty degrees (360 ). Each octant i, ii, iii, iv, v, vi, vii, viii is comprised of a plurality
Referring again to FIG. 1, step 104 can generally involve decomposing the input angle 0102 to a coarse angle A. In this regard, it should be understood that the second determination step 104 can be provided by performing any method known in the art for decomposing an angle into a coarse angle A. For example, the second determination step 104 can involve truncating or rounding a binary representation of the input angle 0102 after a fixed number of bits to obtain a coarse angle A.
In the case of octant symmetry (which is efficient in the concurrent computation of a cosine and sine), the coarse angle A is decomposed in the first octant exclusive of the three (3) most significant bits. If a rounding method is used, an adjustment (not shown) to the fine angle is required to recover the sign. In the case of quadrant symmetry, the coarse angle is decomposed in the first quadrant exclusive of the two (2) most significant bits. The invention is not limited in this regard.
According to an embodiment of the invention, the coarse angle A is an angle in a two dimensional (2D) space. For example, the coarse angle A is one of a plurality of angles in a first octant or quadrant. A graph 200 illustrating a 2D space divided into quadrants and octants is provided in FIG. 2. As shown in FIG. 2, each quadrant I, II, III, IV is comprised of a plurality of angles residing in a respective ninety degree (90 ) area which is non-inclusive of the ninetieth-degree (90 ).
The first quadrant I includes angles between zero degrees (0 ) and ninety degrees (90 ).
The second quadrant II includes angles between ninety degrees (90 ) and one hundred eight degrees (180 ). The third quadrant III includes angles between one hundred eight degrees (180 ) and two hundred seventy degrees (270 ). The fourth quadrant IV includes angles between two hundred seventy degrees (270 ) and three hundred sixty degrees (360 ). Each octant i, ii, iii, iv, v, vi, vii, viii is comprised of a plurality
-7-of angles residing in a respective forty-five degree (45 ) area. For example, the first octant i includes angles of zero degrees (0 ) to forty-five degrees (45 ) which is non-inclusive of the forty-fifth degree (45 ). The invention is not limited in this regard.
Referring again to FIG. 1, the third determination step 106 can generally involve decomposing the input angle 0102 to a fine angle B in the first octant. In this regard, it should understood that the third determination step 106 can be provided by performing any method known in the art for decomposing an input angle 0102 into a fine angle B. For example, the third determination step 106 can involve truncating or rounding a binary representation of an input phase angle after a fixed number of bits to obtain a remainder defining the fine angle B.
The invention is not limited in this regard.
According to an embodiment of the invention, the fine angle B can be bounded by a smallest coarse angle A. For example, if the coarse angle A is in a first octant of a 2D space, then the fine angle B can be bounded by an angle less than or equal to one degree (<1 ). Still, the invention is not limited in this regard. The fine angle B can be bounded by any angle selected in accordance with a particular sine and cosine generation application.
The following Example is provided in order to further illustrate the first, second, and third determination steps 103, 104, 106. The scope of the invention, however, is not to be considered limited in any way thereby.
Assume that an input angle 0102 is sixty five degrees (65 ) or 1.134464014 radians, which is approximately equivalent to 0.0010111000111000111 when scaled to a binary representation (where zero is zero degrees and one is three hundred sixty degrees non-inclusive). As such, only the representations 0.0000000000000000000 through 0.1111111111111111111 are allowed. In such a scenario, the first determination step 103 reduces the input angle 0102 to an angle between zero degrees (0 ) to forty-five degrees (45 ) through use of octant symmetry after the first octant results are computed. The angle between zero degrees (0 ) to
Referring again to FIG. 1, the third determination step 106 can generally involve decomposing the input angle 0102 to a fine angle B in the first octant. In this regard, it should understood that the third determination step 106 can be provided by performing any method known in the art for decomposing an input angle 0102 into a fine angle B. For example, the third determination step 106 can involve truncating or rounding a binary representation of an input phase angle after a fixed number of bits to obtain a remainder defining the fine angle B.
The invention is not limited in this regard.
According to an embodiment of the invention, the fine angle B can be bounded by a smallest coarse angle A. For example, if the coarse angle A is in a first octant of a 2D space, then the fine angle B can be bounded by an angle less than or equal to one degree (<1 ). Still, the invention is not limited in this regard. The fine angle B can be bounded by any angle selected in accordance with a particular sine and cosine generation application.
The following Example is provided in order to further illustrate the first, second, and third determination steps 103, 104, 106. The scope of the invention, however, is not to be considered limited in any way thereby.
Assume that an input angle 0102 is sixty five degrees (65 ) or 1.134464014 radians, which is approximately equivalent to 0.0010111000111000111 when scaled to a binary representation (where zero is zero degrees and one is three hundred sixty degrees non-inclusive). As such, only the representations 0.0000000000000000000 through 0.1111111111111111111 are allowed. In such a scenario, the first determination step 103 reduces the input angle 0102 to an angle between zero degrees (0 ) to forty-five degrees (45 ) through use of octant symmetry after the first octant results are computed. The angle between zero degrees (0 ) to
-8-forty-five degrees (45 ) is non-inclusive of the forty-fifth degree (45 ).
Octant symmetry is described in the following table.
Input angle d Cosine output as a function of Sin output as a function of the first quadrant angle ~' the first quadrant angle ~' (input angle modulo ir/4) (input angle modulo it /4) 0 < ~ < 7t/4 Cos(') Sin(') 7t/4 < ~ < 7t/2 Sin(7t/4-fl Cos(7t/4-fl 7t/2 < ~ < 37t/4 -Sin(') Cos(') 37t/4 < < zt -Cos(n/4-fl Sin(7t/4-fl 7t < ~ < 57t/4 -Cos(') -Sin(') 57t/4 < < 37t/2 -Sin(7t/4-fl - Cos(7t/4-fl 37t/2 < < 77t/4 Sin(') -Cos(') 77t/4 < < 27t Cos(7t/4-fl - Sin(7t/4-fl As shown in the table, octant symmetry refers to three (3) concurrent natural symmetries of sinusoids that may be exploited, wherein the three (3) concurrent natural symmetries include an even symmetry of a cosine function, an odd symmetry of a sine function, and a reflective symmetry of the sine and cosine functions about lines offset from the coordinate axes by ir/4 radians (45 degrees). The result of these symmetry properties is that a hardware implementation of any sinusoid may obtain an effective three (3) bit precision increase from reflecting the input angle to the first quadrant, calculate resulting sine and cosine values in the first octant, and reverse a reflection to obtain actual sine and cosine function evaluations. In
Octant symmetry is described in the following table.
Input angle d Cosine output as a function of Sin output as a function of the first quadrant angle ~' the first quadrant angle ~' (input angle modulo ir/4) (input angle modulo it /4) 0 < ~ < 7t/4 Cos(') Sin(') 7t/4 < ~ < 7t/2 Sin(7t/4-fl Cos(7t/4-fl 7t/2 < ~ < 37t/4 -Sin(') Cos(') 37t/4 < < zt -Cos(n/4-fl Sin(7t/4-fl 7t < ~ < 57t/4 -Cos(') -Sin(') 57t/4 < < 37t/2 -Sin(7t/4-fl - Cos(7t/4-fl 37t/2 < < 77t/4 Sin(') -Cos(') 77t/4 < < 27t Cos(7t/4-fl - Sin(7t/4-fl As shown in the table, octant symmetry refers to three (3) concurrent natural symmetries of sinusoids that may be exploited, wherein the three (3) concurrent natural symmetries include an even symmetry of a cosine function, an odd symmetry of a sine function, and a reflective symmetry of the sine and cosine functions about lines offset from the coordinate axes by ir/4 radians (45 degrees). The result of these symmetry properties is that a hardware implementation of any sinusoid may obtain an effective three (3) bit precision increase from reflecting the input angle to the first quadrant, calculate resulting sine and cosine values in the first octant, and reverse a reflection to obtain actual sine and cosine function evaluations. In
-9-numerical terms, the octant symmetry may be viewed as the signs and the reflection across ir/4 integer multiple boundaries.
After applying octant symmetry rules to the computed result for an angle in the first octant, the input angle OO 102 is reduced to twenty degrees (20 ) or 0.349065850 radians. In a binary representation, this angle is approximately 0.0000111000111000111. The coarse angle A can be selected as the first eight bits following the first three bits used to determine octant symmetry, i.e., the coarse angle A has a binary representation of 01110001. The fine angle B can be selected as the remaining bits to the right of the decimal point of the above listed binary sequence, i.e., the fine angle B has a binary representation of 11000111.
As shown in the above provided table, the third most significant bit (MSB) determines when the cosine and sine angles are inverted about a twenty-two and a half degree (22.5 ) reflected point. The exclusive ORing of the two (2) bits below the most significant bit (MSB) determines when the final computed first octant cosine and sine are transposed to the sine and cosine outputs. The exclusive ORing of the MSB and the bit below the MSB determines when to negate the output sine value while the MSB determines when to negate the output cosine value.
As noted above, the present invention is not limited to octant symmetry based application. For example, the present invention can also be used in quadrant symmetry application. Quadrant symmetry is described in the following table.
Input angle d Cosine output as a function of Sin output as a function of the first quadrant angle ~' the first quadrant angle ~' (input angle modulo ir/4) (input angle modulo ir/4) 0 < ~ < 7c/2 Cos(') Sin(') 7c/2 < ~ < 7c -Sin(') Cos(') 7c < ~ < 37c/2 -Cos(') -Sin(') 3ir/2 < ~ < 2t Sin(') -Cos(')
After applying octant symmetry rules to the computed result for an angle in the first octant, the input angle OO 102 is reduced to twenty degrees (20 ) or 0.349065850 radians. In a binary representation, this angle is approximately 0.0000111000111000111. The coarse angle A can be selected as the first eight bits following the first three bits used to determine octant symmetry, i.e., the coarse angle A has a binary representation of 01110001. The fine angle B can be selected as the remaining bits to the right of the decimal point of the above listed binary sequence, i.e., the fine angle B has a binary representation of 11000111.
As shown in the above provided table, the third most significant bit (MSB) determines when the cosine and sine angles are inverted about a twenty-two and a half degree (22.5 ) reflected point. The exclusive ORing of the two (2) bits below the most significant bit (MSB) determines when the final computed first octant cosine and sine are transposed to the sine and cosine outputs. The exclusive ORing of the MSB and the bit below the MSB determines when to negate the output sine value while the MSB determines when to negate the output cosine value.
As noted above, the present invention is not limited to octant symmetry based application. For example, the present invention can also be used in quadrant symmetry application. Quadrant symmetry is described in the following table.
Input angle d Cosine output as a function of Sin output as a function of the first quadrant angle ~' the first quadrant angle ~' (input angle modulo ir/4) (input angle modulo ir/4) 0 < ~ < 7c/2 Cos(') Sin(') 7c/2 < ~ < 7c -Sin(') Cos(') 7c < ~ < 37c/2 -Cos(') -Sin(') 3ir/2 < ~ < 2t Sin(') -Cos(')
-10-As shown in the table, quadrant symmetry refers to three (3) concurrent natural symmetries of sinusoids that may be exploited, wherein the three (3) concurrent natural symmetries include an even symmetry of a cosine function, an odd symmetry of a sine function, and a reflective symmetry of the sine and cosine functions about lines offset from the coordinate axes by it/2 radians (or 90 degrees).
Referring again to FIG. 1, the process 100 continues with a plurality of computational steps 108, 110, 112, 114, 116, 118. These computational steps 108, 110, 112, 114, 116, 118 are performed for determining a sine and cosine of the input angle 0102 modulo it /k using the previously determined coarse angle A and fine angle B. k is a value selected in accordance with a particular type of symmetry employed.
If the octant symmetry is employed, then k is selected to be four (4) so that it /k has a value equal to forty-five degrees (45 ). In contrast, k is selected to be two (2) when quadrant symmetry is employed so that it/k has a value equal to ninety-degrees (90 ).
The sine and cosine of the input angle 0102 modulo it/k can be defined by the following mathematical expressions (1) and (2). The input angle 0102 modulo k is defined as an input angle 0102.
sin(O' 102) = sin(A)cos(B) + cos(A)sin(B) (1) cos(0' 102) = cos(A)cos(B) - sin(A)sin(B) (2) where:
sin(O' 102) is the sine of the input angle 0102 modulo it/k;
cos(O' 102) is the cosine of the input angle 0102 modulo it/k;
sin(A) is a sine of a coarse angle A;
cos(A) is a cosine of coarse angle B;
cos(B) is a cosine of a fine angle B; and sin(B) is a sine of a fine angle B.
Referring again to FIG. 1, the process 100 continues with a plurality of computational steps 108, 110, 112, 114, 116, 118. These computational steps 108, 110, 112, 114, 116, 118 are performed for determining a sine and cosine of the input angle 0102 modulo it /k using the previously determined coarse angle A and fine angle B. k is a value selected in accordance with a particular type of symmetry employed.
If the octant symmetry is employed, then k is selected to be four (4) so that it /k has a value equal to forty-five degrees (45 ). In contrast, k is selected to be two (2) when quadrant symmetry is employed so that it/k has a value equal to ninety-degrees (90 ).
The sine and cosine of the input angle 0102 modulo it/k can be defined by the following mathematical expressions (1) and (2). The input angle 0102 modulo k is defined as an input angle 0102.
sin(O' 102) = sin(A)cos(B) + cos(A)sin(B) (1) cos(0' 102) = cos(A)cos(B) - sin(A)sin(B) (2) where:
sin(O' 102) is the sine of the input angle 0102 modulo it/k;
cos(O' 102) is the cosine of the input angle 0102 modulo it/k;
sin(A) is a sine of a coarse angle A;
cos(A) is a cosine of coarse angle B;
cos(B) is a cosine of a fine angle B; and sin(B) is a sine of a fine angle B.
-11-However, it should be understood that the mathematical expressions (1) and (2) only provide finite precision sine and cosine trigonometric computations, i.e., the equality statements in mathematical expressions (1) and (2) hold within some small bounded error. The desired output precision drives the precision required in intermediate computations as well as how precise the approximations of terms in mathematical expressions (1) and (2) must be, which in turn leads to more hardware processing. As such, the mathematical expressions (1) and (2) need to be modified to exploit controlled finite precision approximations of the sine and cosine trigonometric computations that are hardware efficient. Such modifications can generally involve computing approximations of the sine of the fine angle B, the cosine of the fine angle B, the sine of the input angle O' 102, and the cosine of the input angle O' 102. The particularities and advantages of these modifications will become more evident as the discussion progresses.
As shown in FIG. 1, the cosine of the coarse angle A is determined during the first computational step 108. The cosine determination can generally involve computing the cosine of the coarse angle A or performing look up table operations to select a pre-computed value for the cosine of the coarse angle A.
Thereafter, the cosine of the coarse angle A is decomposed to a most significant word MSW,os(A) and a least significant word LSW,os(A). The term "decompose" as used herein refers to a process of breaking a value (e.g., cos(A) or sin(A)) to its constituent or component parts. The decomposition process can generally involve computing values for the most significant word of cos(A) (MSW,os(A)) and least significant word of cos(A) (LSW,os(A)). Alternatively, the decomposition process can involve performing look up table operations to select pre-computed values for MSWCos(A) and/or LSWCos(A). It should be noted that the cosine determination and decomposition process can collectively involve performing look up table operations to select pre-computed values for MSWCos(A) and/or LSW,os(A). The phrase "most significant word"
as used herein refers to a set of left most Mi bits of a binary string representing the cosine or sine of the coarse angle A. The phrase "least significant word" as used herein refers to a set of right most M2 bits of a binary string representing the cosine or
As shown in FIG. 1, the cosine of the coarse angle A is determined during the first computational step 108. The cosine determination can generally involve computing the cosine of the coarse angle A or performing look up table operations to select a pre-computed value for the cosine of the coarse angle A.
Thereafter, the cosine of the coarse angle A is decomposed to a most significant word MSW,os(A) and a least significant word LSW,os(A). The term "decompose" as used herein refers to a process of breaking a value (e.g., cos(A) or sin(A)) to its constituent or component parts. The decomposition process can generally involve computing values for the most significant word of cos(A) (MSW,os(A)) and least significant word of cos(A) (LSW,os(A)). Alternatively, the decomposition process can involve performing look up table operations to select pre-computed values for MSWCos(A) and/or LSWCos(A). It should be noted that the cosine determination and decomposition process can collectively involve performing look up table operations to select pre-computed values for MSWCos(A) and/or LSW,os(A). The phrase "most significant word"
as used herein refers to a set of left most Mi bits of a binary string representing the cosine or sine of the coarse angle A. The phrase "least significant word" as used herein refers to a set of right most M2 bits of a binary string representing the cosine or
-12-sine of the coarse angle A. The value of Mi can be the same as or different from the value of M2.
Similarly, the sine of the coarse angle A is determined during the second computational step 110. The sine determination can generally involve computing the sine of the coarse angle A or performing look up table operations to select pre-computed values for the sine of the coarse angle A. Thereafter, the sine of the coarse angle A is decomposed to a most significant word MSWS,f(A) and a least significant word LSWS,f(A). The decomposition process can generally involve computing values for the most significant word of sin(A) (MSWs,f(A)) and least significant word of sin(A) (LSWS,f(A)). Alternatively, the decomposition process can involve performing look up table operations to select pre-computed values for MSWS,f(A) and/or LSWS,f(A). It should be noted that the sine determination and decomposition process can collectively involve performing look up table operations to select pre-computed values for MSWS,f(A) and/or LSWS,f(A).
Notably, the cosine and sine of the coarse angle A are decomposed to a most and least significant word for purposes of efficient hardware (e.g., look up table) implementations. In this regard, it should be understood that the address space of the coarse and fine angles A, B need not be approximately the same. However, the coarse angle A trigonometric values stored in one or more look up tables need to have an accuracy which allows for computing a desired output accuracy. Subsequent operations can be performed to apply angle sum trigonometric identities to the most significant words, thereby minimizing a hardware implementation of the present invention (i.e., less multipliers are needed for a hardware implementation as compared to conventional sine/cosine generators).
Referring again to FIG. 1, an approximation of one minus the cosine of the fine angle B is computed during the third computational step 112. It should be noted that the approximation of 1-cos(B) is computed instead of an approximation of cos(B) for a variety of reasons, such as providing an efficient hardware implementation.
Similarly, the sine of the coarse angle A is determined during the second computational step 110. The sine determination can generally involve computing the sine of the coarse angle A or performing look up table operations to select pre-computed values for the sine of the coarse angle A. Thereafter, the sine of the coarse angle A is decomposed to a most significant word MSWS,f(A) and a least significant word LSWS,f(A). The decomposition process can generally involve computing values for the most significant word of sin(A) (MSWs,f(A)) and least significant word of sin(A) (LSWS,f(A)). Alternatively, the decomposition process can involve performing look up table operations to select pre-computed values for MSWS,f(A) and/or LSWS,f(A). It should be noted that the sine determination and decomposition process can collectively involve performing look up table operations to select pre-computed values for MSWS,f(A) and/or LSWS,f(A).
Notably, the cosine and sine of the coarse angle A are decomposed to a most and least significant word for purposes of efficient hardware (e.g., look up table) implementations. In this regard, it should be understood that the address space of the coarse and fine angles A, B need not be approximately the same. However, the coarse angle A trigonometric values stored in one or more look up tables need to have an accuracy which allows for computing a desired output accuracy. Subsequent operations can be performed to apply angle sum trigonometric identities to the most significant words, thereby minimizing a hardware implementation of the present invention (i.e., less multipliers are needed for a hardware implementation as compared to conventional sine/cosine generators).
Referring again to FIG. 1, an approximation of one minus the cosine of the fine angle B is computed during the third computational step 112. It should be noted that the approximation of 1-cos(B) is computed instead of an approximation of cos(B) for a variety of reasons, such as providing an efficient hardware implementation.
-13-In this regard, it should be understood that if the fine angle B is a small angle then the cosine of the fine angle B is approximately equal to one (1).
Therefore, if the cosine of the fine angle B is greater than or equal to 1-2-16 and an output computation requires a fine angle accuracy of 2-23, then twenty-four (24) bits are required to store the cosine of the fine angle B in a look up table. However, only eight (8) least significant bits are required to store values for one minus the cosine of the fine angle B since the sixteen (16) most significant bits are a priori known to have values equal to zero (0). As a result, no storage is required for and no operations are required to combine the most significant bits of the words with another word.
According to an embodiment of the invention, the approximation computational step 112 can involve Taylor Series Approximations. Taylor Series Approximations are well known to those having ordinary skill in the art, and therefore will not be described in great detail herein. However, it should be understood that the Taylor Series Approximation employed in the third computational step 112 can be defined by the following mathematical equation (3).
1-cos(2nA) =1- [1- [(2nA)2i2!] + [(2nA)4i4!] - [(2nA)6i6!] + ...] (3) Still, the invention is not limited in this regard. Any known method known to those having ordinary skill in the art for computing an approximation of a cosine of an angle can be used.
Referring again to FIG. 1, the sine of the fine angle B is computed during the fourth computational step 114. According to an embodiment of the invention, the approximation computational step 114 can involve Taylor Series Approximations. Taylor Series Approximations are well known to those having ordinary skill in the art, and therefore will not be described in great detail herein.
However, it should be understood that the Taylor Series Approximation employed in the fourth computational step 114 can be defined by the following mathematical equation (4).
sin(2nB) = (2nB) - [(2nB)3/3!] + [(2nB)'i5 ! ] - [(2nB)7/7! ] + ... (4)
Therefore, if the cosine of the fine angle B is greater than or equal to 1-2-16 and an output computation requires a fine angle accuracy of 2-23, then twenty-four (24) bits are required to store the cosine of the fine angle B in a look up table. However, only eight (8) least significant bits are required to store values for one minus the cosine of the fine angle B since the sixteen (16) most significant bits are a priori known to have values equal to zero (0). As a result, no storage is required for and no operations are required to combine the most significant bits of the words with another word.
According to an embodiment of the invention, the approximation computational step 112 can involve Taylor Series Approximations. Taylor Series Approximations are well known to those having ordinary skill in the art, and therefore will not be described in great detail herein. However, it should be understood that the Taylor Series Approximation employed in the third computational step 112 can be defined by the following mathematical equation (3).
1-cos(2nA) =1- [1- [(2nA)2i2!] + [(2nA)4i4!] - [(2nA)6i6!] + ...] (3) Still, the invention is not limited in this regard. Any known method known to those having ordinary skill in the art for computing an approximation of a cosine of an angle can be used.
Referring again to FIG. 1, the sine of the fine angle B is computed during the fourth computational step 114. According to an embodiment of the invention, the approximation computational step 114 can involve Taylor Series Approximations. Taylor Series Approximations are well known to those having ordinary skill in the art, and therefore will not be described in great detail herein.
However, it should be understood that the Taylor Series Approximation employed in the fourth computational step 114 can be defined by the following mathematical equation (4).
sin(2nB) = (2nB) - [(2nB)3/3!] + [(2nB)'i5 ! ] - [(2nB)7/7! ] + ... (4)
-14-Still, the invention is not limited in this regard. Any known method known to those having ordinary skill in the art for computing an approximation of a sine of an angle can be used.
It should be noted that each of the mathematical equations (3) and (4) can be a Taylor Series of an Nth order selected for providing an n bit precision. The n bit precision can be any bit precision selected in accordance with a particular sine/cosine generation application. However, it should be appreciated that such a Taylor Series based sine and cosine generation process 100 overcomes the limitations of conventional sine/cosine generators. For example, the Taylor Series based sine and cosine generation process 100 can efficiently provide a sine and cosine generation with a bit precision greater than an eighteen (18) bit precision. Still, the invention is not limited in this regard. Notably, only two (2) expansion terms in mathematical equations (3) and (4) need to be computed to achieve a twenty-four bit precision, i.e., 1-cos(2riA) = 1- [1 - [(2riA)2/2!] and sin(2riB) = (2riB) - [(2riB)3/3!].
Referring again to FIG. 1, the process 100 continues with the fifth computational step 116. During the fifth computational step 116, a plurality of products P1, P2, P3, P4 are computed. The products Pi, P2, P3, P4 can be defined by the following mathematical equations (5)-(8).
Pi = MSWcos(A)'sin(B)appr (5) P2 = MSWcos(A)'[1-cos(B)appr] (6) P3 = MSWsin(A)'[1-cos(B)appr] (7) P4 = MSWsin(A)'sin(B)appr (8) where:
MSWcos(A) is the most significant word of the cosine of the coarse angle A;
sin(B)appr is the approximation of the sine of the fine angle B;
1-cos(B)appr is the approximation of one minus the cosine of the fine angle B;
and
It should be noted that each of the mathematical equations (3) and (4) can be a Taylor Series of an Nth order selected for providing an n bit precision. The n bit precision can be any bit precision selected in accordance with a particular sine/cosine generation application. However, it should be appreciated that such a Taylor Series based sine and cosine generation process 100 overcomes the limitations of conventional sine/cosine generators. For example, the Taylor Series based sine and cosine generation process 100 can efficiently provide a sine and cosine generation with a bit precision greater than an eighteen (18) bit precision. Still, the invention is not limited in this regard. Notably, only two (2) expansion terms in mathematical equations (3) and (4) need to be computed to achieve a twenty-four bit precision, i.e., 1-cos(2riA) = 1- [1 - [(2riA)2/2!] and sin(2riB) = (2riB) - [(2riB)3/3!].
Referring again to FIG. 1, the process 100 continues with the fifth computational step 116. During the fifth computational step 116, a plurality of products P1, P2, P3, P4 are computed. The products Pi, P2, P3, P4 can be defined by the following mathematical equations (5)-(8).
Pi = MSWcos(A)'sin(B)appr (5) P2 = MSWcos(A)'[1-cos(B)appr] (6) P3 = MSWsin(A)'[1-cos(B)appr] (7) P4 = MSWsin(A)'sin(B)appr (8) where:
MSWcos(A) is the most significant word of the cosine of the coarse angle A;
sin(B)appr is the approximation of the sine of the fine angle B;
1-cos(B)appr is the approximation of one minus the cosine of the fine angle B;
and
-15-MSWS,,,(A) is the most significant word of the sine of the coarse angle A.
After computing the products P1, ..., P4, the process 100 continues with a sixth computational step 118. During the sixth computational step 118, the cosine and sine of the input angle 0102 are computed. In this regard, it should be understood that the sixth computational step 118 can generally involve computing approximations to the cosine and sine of the input angle 0102 using the product terms P1, ..., P4, the coarse angle sinusoid evaluation, cos(A), sin(A), and the selected octant from determination step 103. The approximations in the first octant can be defined by the following mathematical equations (9) and (10).
COS(0' 102)appr = [cos(A)-P2] - P4 =
[cos(A) - [MSWCOS(Aj1-cos(B)appr]]] - [MSWS;,,(A)=sin(B)appr] (9) sin(0' 102)appr = P1 + [sin(A)-P3] =
[MSWcos(A)=sin(B)appr] + [sin(A) - [MSWS;,,(Aj1-cos(B)appr]]] (10) where:
cos(0' 1o2)appr is an approximation to a cosine of the input angle 0102 modulo it/k;
sin(0' 102)appr is an approximation to a sine of the input angle 0102 modulo 7L/k. The first octant approximations are mapped to the required octant in the sixth computational step 118.
The approximations to the cosine and sine of the input angle 0102 (computed in the sixth computational step 118) given in mathematical equations (9) and (10) are not of the desired resolution. The term "resolution" as used herein refers to a numerical accuracy of a result of a computation or series of computations. As such, additional computations must be performed for obtaining the cosine and sine of the input angle 0'102 that are of the desired resolution. The additional computations can be defined by the following mathematical equations (11) and (12).
After computing the products P1, ..., P4, the process 100 continues with a sixth computational step 118. During the sixth computational step 118, the cosine and sine of the input angle 0102 are computed. In this regard, it should be understood that the sixth computational step 118 can generally involve computing approximations to the cosine and sine of the input angle 0102 using the product terms P1, ..., P4, the coarse angle sinusoid evaluation, cos(A), sin(A), and the selected octant from determination step 103. The approximations in the first octant can be defined by the following mathematical equations (9) and (10).
COS(0' 102)appr = [cos(A)-P2] - P4 =
[cos(A) - [MSWCOS(Aj1-cos(B)appr]]] - [MSWS;,,(A)=sin(B)appr] (9) sin(0' 102)appr = P1 + [sin(A)-P3] =
[MSWcos(A)=sin(B)appr] + [sin(A) - [MSWS;,,(Aj1-cos(B)appr]]] (10) where:
cos(0' 1o2)appr is an approximation to a cosine of the input angle 0102 modulo it/k;
sin(0' 102)appr is an approximation to a sine of the input angle 0102 modulo 7L/k. The first octant approximations are mapped to the required octant in the sixth computational step 118.
The approximations to the cosine and sine of the input angle 0102 (computed in the sixth computational step 118) given in mathematical equations (9) and (10) are not of the desired resolution. The term "resolution" as used herein refers to a numerical accuracy of a result of a computation or series of computations. As such, additional computations must be performed for obtaining the cosine and sine of the input angle 0'102 that are of the desired resolution. The additional computations can be defined by the following mathematical equations (11) and (12).
-16-COS(0102) = [Sffi(O102)appr'SF1] + COS(O102)appr (11) sin(0102) = [COS(O102)appr'SF2] + Sm(O102)appr (12) where:
cos(0102) is the sine of the input angle 0102;
sin(0102) is the sine of the input angle 0102;
SF1 is a first scale factor; and SF2 is a second scale factor.
According to an embodiment of the invention, the first and second scale factors SF1, SF2 are a function of the weighting of the most significant word of the cosine of the coarse angle MSW,os(A) and the most significant word of the sine of the coarse angle MSWS,,,(A). The term "function" as used here refers to a mathematical relation in which each element of one word is associated with at least one element of another word. The term "weighting" as used herein refers to a process of aligning binary points of two (2) most significant words on their locations in a combined word and not in separate words. For example, a combined word is defined as a02 + a12-1 + a22-2 + a32-3 + a42-4 + a52-5 + a62-6 + a72-7 + a82-8 + a92-9 + a102-10 +
a112-11 + a122-12 + a132-13 + a142-14 + a152-15 + a162-16 + a172-17 + a182-18 + a192-19 + a202-20 + a212-21 + a222-22 + a232-23. The combined word is broken into the following two (2) words having the same implied binary point positions: a 2 + a12-1 + a22-2 +
a32-3 +
a42-4 + a52-5 + a62-6 + a72-7 + a82-8 + a92-9 + a102-10 + a112-11; and a122 +
a132-1 + a142-2 + a152-3 + a162-4 + a172-5 + a182-6 + a192-7 + a202-8 + a212-9 + a222-10 +
a232-11. In such a scenario, the least significant word must be weighted by 2-12 before the least significant word can be combined with the most significant word to obtain the combined word. Still, the invention is not limited in this regard.
FIGS. 3A-3B and accompanying text described below illustrate an exemplary method 300 according to an embodiment of the invention. The method
cos(0102) is the sine of the input angle 0102;
sin(0102) is the sine of the input angle 0102;
SF1 is a first scale factor; and SF2 is a second scale factor.
According to an embodiment of the invention, the first and second scale factors SF1, SF2 are a function of the weighting of the most significant word of the cosine of the coarse angle MSW,os(A) and the most significant word of the sine of the coarse angle MSWS,,,(A). The term "function" as used here refers to a mathematical relation in which each element of one word is associated with at least one element of another word. The term "weighting" as used herein refers to a process of aligning binary points of two (2) most significant words on their locations in a combined word and not in separate words. For example, a combined word is defined as a02 + a12-1 + a22-2 + a32-3 + a42-4 + a52-5 + a62-6 + a72-7 + a82-8 + a92-9 + a102-10 +
a112-11 + a122-12 + a132-13 + a142-14 + a152-15 + a162-16 + a172-17 + a182-18 + a192-19 + a202-20 + a212-21 + a222-22 + a232-23. The combined word is broken into the following two (2) words having the same implied binary point positions: a 2 + a12-1 + a22-2 +
a32-3 +
a42-4 + a52-5 + a62-6 + a72-7 + a82-8 + a92-9 + a102-10 + a112-11; and a122 +
a132-1 + a142-2 + a152-3 + a162-4 + a172-5 + a182-6 + a192-7 + a202-8 + a212-9 + a222-10 +
a232-11. In such a scenario, the least significant word must be weighted by 2-12 before the least significant word can be combined with the most significant word to obtain the combined word. Still, the invention is not limited in this regard.
FIGS. 3A-3B and accompanying text described below illustrate an exemplary method 300 according to an embodiment of the invention. The method
-17-300 is provided for generating a sine and cosine of an input angle 0102. It should be appreciated, however, that the method 300 disclosed herein is provided for purposes of illustration only and that the present invention is not limited solely to the particular method shown.
Referring now to FIG. 3A, there is provided a method for generating a sine and cosine of an input angle 0102 that is useful for understanding the present invention. As shown in FIG. 3A, the method begins at step 302 and continues with step 304. In step 304, an input angle 0102 is decomposed into an octant (or quadrant), a coarse angle A, and a fine angle B. As noted above, the fine angle B can be bounded by a smallest coarse angle A.
After step 304, the method 300 continues with step 306. In step 306, a cosine of the coarse angle A is computed or determined using look-up table operations. Look-up table operations are well known to those having ordinary skill in the art, and therefore will not be described in herein. Subsequently, step 308 is performed where the sine of the coarse angle A is computed or determined using look-up table operations. Thereafter, steps 310 and 312 are performed for decomposing cos(A) and sin(A) to a most significant word and a least significant word.
In step 314, an approximation of one minus the cosine of the fine angle B is computed. As noted above, this approximation can be obtained using a Taylor Series Approximation. Similarly, an approximation of the sine of the fine angle B is computed in step 316. This approximation can also be obtained using a Taylor Series Approximation.
After completing step 316, step 318 is performed where a plurality of products P1, ..., P4 are computed using the previously computed or determined values for the most significant word of the cosine of the coarse angle A, the most significant word of the sine of the coarse angle A, the approximation of the sine of the fine angle B, and the approximation of one minus the cosine of the fine angle B. Each of the products P1, ..., P4 is described above in relation to FIG. 1 and mathematical equations (5) through (8). The description above is sufficient for understanding the
Referring now to FIG. 3A, there is provided a method for generating a sine and cosine of an input angle 0102 that is useful for understanding the present invention. As shown in FIG. 3A, the method begins at step 302 and continues with step 304. In step 304, an input angle 0102 is decomposed into an octant (or quadrant), a coarse angle A, and a fine angle B. As noted above, the fine angle B can be bounded by a smallest coarse angle A.
After step 304, the method 300 continues with step 306. In step 306, a cosine of the coarse angle A is computed or determined using look-up table operations. Look-up table operations are well known to those having ordinary skill in the art, and therefore will not be described in herein. Subsequently, step 308 is performed where the sine of the coarse angle A is computed or determined using look-up table operations. Thereafter, steps 310 and 312 are performed for decomposing cos(A) and sin(A) to a most significant word and a least significant word.
In step 314, an approximation of one minus the cosine of the fine angle B is computed. As noted above, this approximation can be obtained using a Taylor Series Approximation. Similarly, an approximation of the sine of the fine angle B is computed in step 316. This approximation can also be obtained using a Taylor Series Approximation.
After completing step 316, step 318 is performed where a plurality of products P1, ..., P4 are computed using the previously computed or determined values for the most significant word of the cosine of the coarse angle A, the most significant word of the sine of the coarse angle A, the approximation of the sine of the fine angle B, and the approximation of one minus the cosine of the fine angle B. Each of the products P1, ..., P4 is described above in relation to FIG. 1 and mathematical equations (5) through (8). The description above is sufficient for understanding the
-18-product P1, ..., P4 computations. Notably, the level of precision in calculating the product terms varies P1, ..., P4 and is chosen based on the overall desired output precision. Subsequent to computing the products Pi, ..., P4, the method 300 continues with a step 320 of FIG. 3B.
Referring now to FIG. 3B, a first difference D1 is computed in step 320. The first difference is obtained by subtracting the second product P2 from the cosine of the coarse angle A. Accordingly, the first difference can be defined by the following mathematical equation (13).
Di = cos(A) - P2 (13) Similarly, a second difference D2 is computed in step 322. The second difference D2 is obtained by subtracting the third product P3 from the sine of the coarse angle A. As such, the second difference D2 can be defined by the following mathematical equation (14).
D2 = sin(A) - P3 (14) In step 324, an approximation of the cosine of the input angle 0'102 is generated by computing a third difference D3. The third difference D3 is computed by subtracting the fourth product P4 from the first difference D1. Accordingly, the approximation of the cosine of the input angle cos(D' 1o2)appr can be defined by the above recited mathematical equation (9) and the following mathematical equation (15).
COS(O'102)appr = D3 = D1 - P4 (15) After completing step 324, the method 300 continues with step 326.
In step 326, an approximation to the sine of the input angle 0'102 is generated by computing a first sum S1. The first sum S1 is computed by adding together the first product P1 and the second difference D2. As such, the approximation of the sine of the input angle sin(D' 102)appr can be defined by the above recited mathematical equation (10) and the following mathematical equation (16).
Referring now to FIG. 3B, a first difference D1 is computed in step 320. The first difference is obtained by subtracting the second product P2 from the cosine of the coarse angle A. Accordingly, the first difference can be defined by the following mathematical equation (13).
Di = cos(A) - P2 (13) Similarly, a second difference D2 is computed in step 322. The second difference D2 is obtained by subtracting the third product P3 from the sine of the coarse angle A. As such, the second difference D2 can be defined by the following mathematical equation (14).
D2 = sin(A) - P3 (14) In step 324, an approximation of the cosine of the input angle 0'102 is generated by computing a third difference D3. The third difference D3 is computed by subtracting the fourth product P4 from the first difference D1. Accordingly, the approximation of the cosine of the input angle cos(D' 1o2)appr can be defined by the above recited mathematical equation (9) and the following mathematical equation (15).
COS(O'102)appr = D3 = D1 - P4 (15) After completing step 324, the method 300 continues with step 326.
In step 326, an approximation to the sine of the input angle 0'102 is generated by computing a first sum S1. The first sum S1 is computed by adding together the first product P1 and the second difference D2. As such, the approximation of the sine of the input angle sin(D' 102)appr can be defined by the above recited mathematical equation (10) and the following mathematical equation (16).
-19-sin(O'102)appr = D2 + P1 (16) Subsequently, the method 300 continues with step 328.
In step 328, the approximation of cos(O'1o2) is scaled by a first scale factor SF1 to place the same in a desired resolution. Step 328 also involves scaling the approximation of sin(O' 102) by a second scale factor SF2 to place the same in a desired resolution. These scaling operations can be defined by the following mathematical equations (17) and (18).
COS(O'102)appr' = COS(O'102)appr SF1 = D3 - SF1 (17) sin(O' 102)appr' = sin(O' 102)appr SF2 = S 1 = SF2 (18) After scaling the approximation of the cosine and sine of the input angle 0'102, the method 300 continues with step 330. In step 330, the cosine of the input angle 0'102 is computed (or determined). More particularly, the scaled approximation of the sine of the input angle sin(O' 102)appr' is combined with the approximation of the cosine of the input angle cos(O' 1o2)appr to yield the cosine of the input angle 0'102. This computation (or determination) can be defined by the above recited mathematical equation (11) and the following mathematical equation (19).
COS(O' 102) = sin(O' 102)appr' + COS(O' 102)appr (19) Similarly, the sine of the input angle 0'102 is computed or determined in step 332. Step 332 can generally involve combining the scaled approximation of the cosine of the input angle cos(O' 102) appr' with the approximation of the sine of the input angle sin(O' 102)appr to yield the sine of the input angle 0'102. This computation or determination can be defined by the above recited mathematical equation (12) and the following mathematical equation (20).
sm(O' 102) = COS(O' 102)appr' + sin(O' 102)appr (20)
In step 328, the approximation of cos(O'1o2) is scaled by a first scale factor SF1 to place the same in a desired resolution. Step 328 also involves scaling the approximation of sin(O' 102) by a second scale factor SF2 to place the same in a desired resolution. These scaling operations can be defined by the following mathematical equations (17) and (18).
COS(O'102)appr' = COS(O'102)appr SF1 = D3 - SF1 (17) sin(O' 102)appr' = sin(O' 102)appr SF2 = S 1 = SF2 (18) After scaling the approximation of the cosine and sine of the input angle 0'102, the method 300 continues with step 330. In step 330, the cosine of the input angle 0'102 is computed (or determined). More particularly, the scaled approximation of the sine of the input angle sin(O' 102)appr' is combined with the approximation of the cosine of the input angle cos(O' 1o2)appr to yield the cosine of the input angle 0'102. This computation (or determination) can be defined by the above recited mathematical equation (11) and the following mathematical equation (19).
COS(O' 102) = sin(O' 102)appr' + COS(O' 102)appr (19) Similarly, the sine of the input angle 0'102 is computed or determined in step 332. Step 332 can generally involve combining the scaled approximation of the cosine of the input angle cos(O' 102) appr' with the approximation of the sine of the input angle sin(O' 102)appr to yield the sine of the input angle 0'102. This computation or determination can be defined by the above recited mathematical equation (12) and the following mathematical equation (20).
sm(O' 102) = COS(O' 102)appr' + sin(O' 102)appr (20)
-20-Notably, the scaled additions are provided as an adjustment that makes the hardware implementation more efficient and may be replaced by similar approximations without loss of generality. After completing step 332, step 334 is performed.
In step 334, N MSBs from the input angle 0102 are used to control the mapping from sin(D'102) and cos(D'102) values to the sin(D102) and cos(D102) values described previously in the mapping table. N is selected in accordance with a particular type of symmetry employed. For example, if octant symmetry is employed, then N is selected to be three (3). If quadrant symmetry is employed, then N is selected to be two (2). After completing step 334, step 336 is performed where the method 300 ends.
Referring now to FIG. 4, there is provided an exemplary sine/cosine generator 400 that is useful for understanding the present invention. It should be noted that the sine/cosine generator 400 implements the method 300 described above in relation to FIGS. 3A and 3B. As shown in FIG. 4, the sine/cosine generator 400 is comprised of an octant segmenting device 403, a coarse angle device (CAD) 404, a fine angle device (FAD) 406, a cosine computation and decomposition device (CCDD) 408, a sine computation and decomposition device (SCDD) 410, and Taylor Series Approximation devices (TSADs) 412, 414. The sine/cosine generator 400 is also comprised of a product device 416 and a cosine and sine computation device (CSCD) 418.
The octant segmenting device 403 is comprised of hardware and/or software configured for receiving an input angle 0102 and determining the octant of the input angle 0102. It should be understood that that octant segmenting device 403 can employ any method known to those having ordinary skill in the art for decomposing an angle into an octant designation. The octant segmenting device is also comprised of hardware and/or software configured for communicating the determined octant to the CSCD 418. The invention is not limited in this regard. For example, the sine/cosine generator 400 can alternatively include a quadrant segmenting device (not shown).
In step 334, N MSBs from the input angle 0102 are used to control the mapping from sin(D'102) and cos(D'102) values to the sin(D102) and cos(D102) values described previously in the mapping table. N is selected in accordance with a particular type of symmetry employed. For example, if octant symmetry is employed, then N is selected to be three (3). If quadrant symmetry is employed, then N is selected to be two (2). After completing step 334, step 336 is performed where the method 300 ends.
Referring now to FIG. 4, there is provided an exemplary sine/cosine generator 400 that is useful for understanding the present invention. It should be noted that the sine/cosine generator 400 implements the method 300 described above in relation to FIGS. 3A and 3B. As shown in FIG. 4, the sine/cosine generator 400 is comprised of an octant segmenting device 403, a coarse angle device (CAD) 404, a fine angle device (FAD) 406, a cosine computation and decomposition device (CCDD) 408, a sine computation and decomposition device (SCDD) 410, and Taylor Series Approximation devices (TSADs) 412, 414. The sine/cosine generator 400 is also comprised of a product device 416 and a cosine and sine computation device (CSCD) 418.
The octant segmenting device 403 is comprised of hardware and/or software configured for receiving an input angle 0102 and determining the octant of the input angle 0102. It should be understood that that octant segmenting device 403 can employ any method known to those having ordinary skill in the art for decomposing an angle into an octant designation. The octant segmenting device is also comprised of hardware and/or software configured for communicating the determined octant to the CSCD 418. The invention is not limited in this regard. For example, the sine/cosine generator 400 can alternatively include a quadrant segmenting device (not shown).
-21-Referring again to FIG. 4, the CAD 404 is comprised of hardware and/or software configured for receiving an input angle 0102 and determining the coarse angle A of the same. As noted above, the coarse angle A can be one of a plurality of angles selected from the first octant (or quadrant if quadrant symmetry is employed). The CAD 404 is also comprised of hardware and/or software configured for communicating the coarse angle A to the CCDD 408 and the SCDD 410.
The CCDD 408 is comprised of hardware and/or software configured for computing or determining the cosine of the coarse angle A and decomposing the same to a most significant word and a least significant word. The CCDD 408 is also comprised of hardware and/or software configured for communicating the most significant word of the cosine of the coarse angle A to the product device 416. The CCDD 408 is further comprised of hardware and/or software configured for communicating the cosine of the coarse angle A to the CSCD 418.
It should be understood that the CCDD 408 can be comprised of one or more look up tables containing pre-computed values for the cosine of the coarse angle A, the most significant word, and/or the least significant word. The pre-computed values can be expressed as a binary number system representation. Stated differently, each of the pre-computed values can be expressed as a string of bits, wherein each bit has a one (1) value or a zero (0) value. Still, the invention is not limited in this regard.
The pre-computed values can be expressed in any number system representation selected in accordance with a particular sine/cosine generator 400 application.
According to an embodiment of the invention, the least significant word is scaled such that all bits stored in the look up table(s) are significant. Also, the results of the arithmetic operations based on the content of the look up table(s) are scaled as the inverse of the look up table scaling. For example, a combined word is defined as a02 + a12-1 + a22-2 + a32-3 + a42-4 + a52-5 + a62-6 + a72-7 + a82-8 + a92-9 +
a102-10 + a112-11 + a122-12 + a132-13 + a142-14 + a152-15 + a162-16 + a172-17 + a182-18 + a192-19 + a202-20 + a212-21 + a222-22 + a232-23. The combined word is broken into the following two (2) words having the same implied binary point positions: a02 + a12-1 +
a22-2 +
a32-3 + a42-4 + a52-5 + a62-6 + a72-7 + a82-8 + a92-9 + a102-10 + a112-11 +
0X2-12 + 0X2-13 +
The CCDD 408 is comprised of hardware and/or software configured for computing or determining the cosine of the coarse angle A and decomposing the same to a most significant word and a least significant word. The CCDD 408 is also comprised of hardware and/or software configured for communicating the most significant word of the cosine of the coarse angle A to the product device 416. The CCDD 408 is further comprised of hardware and/or software configured for communicating the cosine of the coarse angle A to the CSCD 418.
It should be understood that the CCDD 408 can be comprised of one or more look up tables containing pre-computed values for the cosine of the coarse angle A, the most significant word, and/or the least significant word. The pre-computed values can be expressed as a binary number system representation. Stated differently, each of the pre-computed values can be expressed as a string of bits, wherein each bit has a one (1) value or a zero (0) value. Still, the invention is not limited in this regard.
The pre-computed values can be expressed in any number system representation selected in accordance with a particular sine/cosine generator 400 application.
According to an embodiment of the invention, the least significant word is scaled such that all bits stored in the look up table(s) are significant. Also, the results of the arithmetic operations based on the content of the look up table(s) are scaled as the inverse of the look up table scaling. For example, a combined word is defined as a02 + a12-1 + a22-2 + a32-3 + a42-4 + a52-5 + a62-6 + a72-7 + a82-8 + a92-9 +
a102-10 + a112-11 + a122-12 + a132-13 + a142-14 + a152-15 + a162-16 + a172-17 + a182-18 + a192-19 + a202-20 + a212-21 + a222-22 + a232-23. The combined word is broken into the following two (2) words having the same implied binary point positions: a02 + a12-1 +
a22-2 +
a32-3 + a42-4 + a52-5 + a62-6 + a72-7 + a82-8 + a92-9 + a102-10 + a112-11 +
0X2-12 + 0X2-13 +
-22-0X2-14 + 0X2-15 + 0X2-16 + 0X2-17 + 0X2-18 + 0X2-19 + 0X2-20 + 0X2-21 + 0X2-22 +
0X2-23; and 0X2 + 0X2-1 + 0X2-2 + 0X2-3 + 0X2-4 + 0X2-5 + 0X2-6 + 0X2-7 + 0X2-8 +
0X2-9 + 0X2-10 + 0X2-11 + a122-12 + a132-13 + a142-14 + a152-15 + a162-16 +
a172-17 + a182-18 + a192-19 + a202-2 + a212-21 + a222-22 + a232-23. As should be understood, the digits that are identically zero (0) have no significance. Therefore, only the significant digits of the two (2) words need to be stored. More particularly, the last twelve (12) bits of the first word are discarded and the first twelve (12) bits are stored as the most significant word. Similarly, the first twelve (12) bits of the second word are discarded, the bit weights of the last twelve (12) bits of the second word are scaled by 212, and the scaled last twelve (12) bits are stored as the least significant word. Still, the invention is not limited in this regard.
Referring again to FIG. 4, the SCDD 410 is comprised of hardware and/or software configured for computing or determining the sine of the coarse angle A and decomposing the same to a most significant word and a least significant word.
The SCDD 410 is also comprised of hardware and/or software configured for communicating the most significant word of the sine of the coarse angle A to the product device 416. The SCDD 410 is further comprised of hardware and/or software configured for communicating the sine of the coarse angle A to the CSCD 418.
It should be understood that the SCDD 410 can be comprised of one or more look up tables containing pre-computed values for the sine of the coarse angle A, the most significant word, and/or the least significant word. The pre-computed values can be expressed as a binary number system representation. Stated differently, each of the pre-computed values can be expressed as a string of bits, each bit having a one (1) value or a zero (0) value. Still, the invention is not limited in this regard. The pre-computed values can be expressed in any number system representation selected in accordance with a particular sine/cosine generator 400 application.
According to an embodiment of the invention, the least significant word is scaled such that all bits stored in the look up table(s) are significant. Also, the results of the arithmetic operations based on the content of the look up table(s) are
0X2-23; and 0X2 + 0X2-1 + 0X2-2 + 0X2-3 + 0X2-4 + 0X2-5 + 0X2-6 + 0X2-7 + 0X2-8 +
0X2-9 + 0X2-10 + 0X2-11 + a122-12 + a132-13 + a142-14 + a152-15 + a162-16 +
a172-17 + a182-18 + a192-19 + a202-2 + a212-21 + a222-22 + a232-23. As should be understood, the digits that are identically zero (0) have no significance. Therefore, only the significant digits of the two (2) words need to be stored. More particularly, the last twelve (12) bits of the first word are discarded and the first twelve (12) bits are stored as the most significant word. Similarly, the first twelve (12) bits of the second word are discarded, the bit weights of the last twelve (12) bits of the second word are scaled by 212, and the scaled last twelve (12) bits are stored as the least significant word. Still, the invention is not limited in this regard.
Referring again to FIG. 4, the SCDD 410 is comprised of hardware and/or software configured for computing or determining the sine of the coarse angle A and decomposing the same to a most significant word and a least significant word.
The SCDD 410 is also comprised of hardware and/or software configured for communicating the most significant word of the sine of the coarse angle A to the product device 416. The SCDD 410 is further comprised of hardware and/or software configured for communicating the sine of the coarse angle A to the CSCD 418.
It should be understood that the SCDD 410 can be comprised of one or more look up tables containing pre-computed values for the sine of the coarse angle A, the most significant word, and/or the least significant word. The pre-computed values can be expressed as a binary number system representation. Stated differently, each of the pre-computed values can be expressed as a string of bits, each bit having a one (1) value or a zero (0) value. Still, the invention is not limited in this regard. The pre-computed values can be expressed in any number system representation selected in accordance with a particular sine/cosine generator 400 application.
According to an embodiment of the invention, the least significant word is scaled such that all bits stored in the look up table(s) are significant. Also, the results of the arithmetic operations based on the content of the look up table(s) are
-23-scaled as the inverse of the look up table scaling. Still, the invention is not limited in this regard.
Referring again to FIG. 4, the FAD 406 is comprised of hardware and/or software configured for receiving an input angle 0102 and determining the fine angle B of the same. As noted above, the fine angle B can be bounded by the smallest coarse angle A. The FAD 406 is also comprised of hardware and/or software configured for communicating the fine angle B to the TSADs 412, 414.
The TSAD 412 is comprised of hardware and/or software configured for computing or determining an approximation of one minus the cosine of the fine angle B. This computation or determination can involve a Taylor Series Approximation. The Taylor Series Approximation can involve using a Taylor Series of an Nth order selected for providing at least an n bit precision output. The TSAD
412 is also comprised of hardware and/or software configured for communicating the approximation of one minus the cosine of the fine angle B to the product device 416.
It should be noted that the TSAD 412 can be comprised of one or more look up tables containing pre-computed values for the approximation of one minus the cosine of the fine angle B or other fine angle approximations.
The TSAD 414 is comprised of hardware and/or software configured for computing or determining an approximation of sine of the fine angle B.
This computation or determination can involve a Taylor Series Approximation. The Taylor Series Approximation can involve using a Taylor Series of an Nth order selected for providing an n bit precision. The TSAD 414 is also comprised of hardware and/or software configured for communicating the approximation of sine of the fine angle B to the product device 416. It should be noted that the TSAD
414 can be comprised of one or more look up tables containing pre-computed values for the approximation of sine of the fine angle B or other fine angle approximations.
The product device 416 is comprised of hardware and/or software configured for computing a plurality of products P1, ..., P4. The products P1, ..., P4 are computed using the pre-computed or pre-determined values for the most significant word of the cosine of the coarse angle A, the most significant word of the
Referring again to FIG. 4, the FAD 406 is comprised of hardware and/or software configured for receiving an input angle 0102 and determining the fine angle B of the same. As noted above, the fine angle B can be bounded by the smallest coarse angle A. The FAD 406 is also comprised of hardware and/or software configured for communicating the fine angle B to the TSADs 412, 414.
The TSAD 412 is comprised of hardware and/or software configured for computing or determining an approximation of one minus the cosine of the fine angle B. This computation or determination can involve a Taylor Series Approximation. The Taylor Series Approximation can involve using a Taylor Series of an Nth order selected for providing at least an n bit precision output. The TSAD
412 is also comprised of hardware and/or software configured for communicating the approximation of one minus the cosine of the fine angle B to the product device 416.
It should be noted that the TSAD 412 can be comprised of one or more look up tables containing pre-computed values for the approximation of one minus the cosine of the fine angle B or other fine angle approximations.
The TSAD 414 is comprised of hardware and/or software configured for computing or determining an approximation of sine of the fine angle B.
This computation or determination can involve a Taylor Series Approximation. The Taylor Series Approximation can involve using a Taylor Series of an Nth order selected for providing an n bit precision. The TSAD 414 is also comprised of hardware and/or software configured for communicating the approximation of sine of the fine angle B to the product device 416. It should be noted that the TSAD
414 can be comprised of one or more look up tables containing pre-computed values for the approximation of sine of the fine angle B or other fine angle approximations.
The product device 416 is comprised of hardware and/or software configured for computing a plurality of products P1, ..., P4. The products P1, ..., P4 are computed using the pre-computed or pre-determined values for the most significant word of the cosine of the coarse angle A, the most significant word of the
-24-sine of the coarse angle A, the approximation of one minus the sine of the fine angle B, and the approximation of the sine of the fine angle B. The products Pi, ..., P4 are described in detail above in relation to FIG. 1. The description provided above is sufficient for understanding the product P1, ..., P4 computations. It should be noted that the product device 416 can be comprised of a plurality of multipliers.
Multipliers are well known to those having ordinary skill in the art, and therefore will not be described herein.
The product device 416 is also comprised of hardware and/or software configured for communicating the products Pi,. .., P4 to the CSCD 418. The CSCD
418 is comprised of hardware and/or software configured for generating two (2) outputs. A first output represents the cosine of the angle 0'102. A second output represents the sine of the angle 0'102. The outputs can be generated by performing a plurality of computations using the cosine of the coarse angle A, the sine of the coarse angle A, the products P1, ..., P4, scale factors SF1, SF2 (described above in relation to FIG. 1), and octant mappings. The computations are described above in relation to FIGS. 1 and 3A-3B. The description provided above is sufficient for understanding the computations performed for generating the cosine of the input angle 0102 and the sine of the input angle 0102.
Referring now to FIG. 5, there is provided another exemplary sine/cosine generator 500 that is useful for understanding the present invention. It should be noted that the sine/cosine generator 500 implements the method 300 described above in relation to FIGS. 3A and 3B minus the octant mapping of step 334. As shown in FIG. 5, the sine/cosine generator 500 is comprised of a coarse address/address compliment device 502, a coarse cosine look up table (LUT) 504, a coarse sine LUT 506, a cosine address device 518, a single sided address to angle conversion device 522, a symmetric 1-cos(B) LUT 516, a sine/sine complement device 518, and multipliers 508, 510, 512, 514.
It should be noted that the coarse cosine LUT 504 contains pre-computed values for the cosine of the coarse angle A and the most significant words for the cosine of the coarse angle A. By taking the calculated number of most
Multipliers are well known to those having ordinary skill in the art, and therefore will not be described herein.
The product device 416 is also comprised of hardware and/or software configured for communicating the products Pi,. .., P4 to the CSCD 418. The CSCD
418 is comprised of hardware and/or software configured for generating two (2) outputs. A first output represents the cosine of the angle 0'102. A second output represents the sine of the angle 0'102. The outputs can be generated by performing a plurality of computations using the cosine of the coarse angle A, the sine of the coarse angle A, the products P1, ..., P4, scale factors SF1, SF2 (described above in relation to FIG. 1), and octant mappings. The computations are described above in relation to FIGS. 1 and 3A-3B. The description provided above is sufficient for understanding the computations performed for generating the cosine of the input angle 0102 and the sine of the input angle 0102.
Referring now to FIG. 5, there is provided another exemplary sine/cosine generator 500 that is useful for understanding the present invention. It should be noted that the sine/cosine generator 500 implements the method 300 described above in relation to FIGS. 3A and 3B minus the octant mapping of step 334. As shown in FIG. 5, the sine/cosine generator 500 is comprised of a coarse address/address compliment device 502, a coarse cosine look up table (LUT) 504, a coarse sine LUT 506, a cosine address device 518, a single sided address to angle conversion device 522, a symmetric 1-cos(B) LUT 516, a sine/sine complement device 518, and multipliers 508, 510, 512, 514.
It should be noted that the coarse cosine LUT 504 contains pre-computed values for the cosine of the coarse angle A and the most significant words for the cosine of the coarse angle A. By taking the calculated number of most
-25-significant bits (exclusive of the three most significant bits) at this point of the sin(input angle) and cos(input angle) computations, the size of the multiplications required to compute the sine and cosine of an input angle is reduced. It should also be noted that the number of bits for each most significant word is a function of the desired output accuracy of the sine/cosine generator 500. It should further be noted that each of the multipliers 508, 510, 512, 514 is provided for computing a respective product P1, ..., P4 (described above in relation to FIGS. 1 and 3A-3B).
Referring again to FIG. 5, the sine/cosine generator 500 is also comprised of scalers 524, 528, 550, 552, multipliers 526, 530, 546, 548, and multiplexers (MUXs) 532, 542. The sine/cosine generator 500 is further comprised of scalers 536, 540, compliment devices 534, 544, and adders 534, 538. It should be noted that the respective outputs of the compliment devices 534, 544 are the approximations of the cosine and sine of the first octant equivalent of the input angle.
However, the approximations are not to the desired output accuracy. As a result, small compensation terms must be added by the scalers 536, 540 to supply the desired output accuracy. The compensation terms are a function of the input angle.
In light of the forgoing description of the invention, it should be recognized that the present invention can be realized in hardware, software, or a combination of hardware and software. A method for generating a sine and cosine of an input angle according to the present invention can be realized in a centralized fashion in one processing system, or in a distributed fashion where different elements are spread across several interconnected processing systems. Any kind of computer system, or other apparatus adapted for carrying out the methods described herein, is suited. A typical combination of hardware and software could be a general purpose computer processor, with a computer program that, when being loaded and executed, controls the computer processor such that it carries out the methods described herein.
Of course, an application specific integrated circuit (ASIC), and/or a field programmable gate array (FPGA) could also be used to achieve a similar result.
The present invention can also be embedded in a computer program product, which comprises all the features enabling the implementation of the methods
Referring again to FIG. 5, the sine/cosine generator 500 is also comprised of scalers 524, 528, 550, 552, multipliers 526, 530, 546, 548, and multiplexers (MUXs) 532, 542. The sine/cosine generator 500 is further comprised of scalers 536, 540, compliment devices 534, 544, and adders 534, 538. It should be noted that the respective outputs of the compliment devices 534, 544 are the approximations of the cosine and sine of the first octant equivalent of the input angle.
However, the approximations are not to the desired output accuracy. As a result, small compensation terms must be added by the scalers 536, 540 to supply the desired output accuracy. The compensation terms are a function of the input angle.
In light of the forgoing description of the invention, it should be recognized that the present invention can be realized in hardware, software, or a combination of hardware and software. A method for generating a sine and cosine of an input angle according to the present invention can be realized in a centralized fashion in one processing system, or in a distributed fashion where different elements are spread across several interconnected processing systems. Any kind of computer system, or other apparatus adapted for carrying out the methods described herein, is suited. A typical combination of hardware and software could be a general purpose computer processor, with a computer program that, when being loaded and executed, controls the computer processor such that it carries out the methods described herein.
Of course, an application specific integrated circuit (ASIC), and/or a field programmable gate array (FPGA) could also be used to achieve a similar result.
The present invention can also be embedded in a computer program product, which comprises all the features enabling the implementation of the methods
-26-described herein, and which, when loaded in a computer system, is able to carry out these methods. Computer program or application in the present context means any expression, in any language, code or notation, of a set of instructions intended to cause a system having an information processing capability to perform a particular function either directly or after either or both of the following: (a) conversion to another language, code or notation; (b) reproduction in a different material form.
Additionally, the description above is intended by way of example only and is not intended to limit the present invention in any way, except as set forth in the following claims.
Additionally, the description above is intended by way of example only and is not intended to limit the present invention in any way, except as set forth in the following claims.
-27-
Claims (10)
1. A sine/cosine generator operable to generate a sine and cosine of an input angle utilizing digital circuitry, or one or more processors executing computer instructions, or a combination thereof, said sine/cosine generator configured to:
decompose said input angle into an octant or quadrant, a coarse angle, and a fine angle;
determine a cosine and a sine of said coarse angle;
decompose said cosine of said coarse angle to a first most significant word and a first least significant word and said sine of said coarse angle to a second most significant word and a second least significant word;
compute an approximation of one minus a cosine of said fine angle and an approximation of a sine of said fine angle;
compute a first product by multiplying said first most significant word by said approximation of a sine of said fine angle, a second product by multiplying said first most significant word by said approximation of one minus a cosine of said fine angle, a third product by multiplying said second most significant word by said approximation of one minus a cosine of said fine angle, and a fourth product by multiplying said second most significant word by said approximation of a sine of said fine angle;
generate approximations of a cosine of said input angle modulo .pi./k and a sine of said input angle modulo .pi./k using said cosine of said coarse angle, said sine of said coarse angle, and each of said products, wherein k is an integer selected in accordance with octant or quadrant symmetry; and using the sine and cosine of the input angle in a communication system to digitally up convert a signal, digitally down convert the signal, or transform the signal.
decompose said input angle into an octant or quadrant, a coarse angle, and a fine angle;
determine a cosine and a sine of said coarse angle;
decompose said cosine of said coarse angle to a first most significant word and a first least significant word and said sine of said coarse angle to a second most significant word and a second least significant word;
compute an approximation of one minus a cosine of said fine angle and an approximation of a sine of said fine angle;
compute a first product by multiplying said first most significant word by said approximation of a sine of said fine angle, a second product by multiplying said first most significant word by said approximation of one minus a cosine of said fine angle, a third product by multiplying said second most significant word by said approximation of one minus a cosine of said fine angle, and a fourth product by multiplying said second most significant word by said approximation of a sine of said fine angle;
generate approximations of a cosine of said input angle modulo .pi./k and a sine of said input angle modulo .pi./k using said cosine of said coarse angle, said sine of said coarse angle, and each of said products, wherein k is an integer selected in accordance with octant or quadrant symmetry; and using the sine and cosine of the input angle in a communication system to digitally up convert a signal, digitally down convert the signal, or transform the signal.
2. The sine/cosine generator according to claim 1, wherein said generate feature further includes:
computing a first difference by subtracting said second product from said cosine of said coarse angle and a second difference by subtracting said third product from said sine of said coarse angle, subtracting said fourth product from said first difference to yield said approximation of said cosine of said input angle modulo .pi./k, and combining said first product and said second difference to yield said approximation of said sine of said input angle modulo .pi./k.
computing a first difference by subtracting said second product from said cosine of said coarse angle and a second difference by subtracting said third product from said sine of said coarse angle, subtracting said fourth product from said first difference to yield said approximation of said cosine of said input angle modulo .pi./k, and combining said first product and said second difference to yield said approximation of said sine of said input angle modulo .pi./k.
3. The sine/cosine generator according to claim 1, wherein said sine/cosine generator is further configured to scale said approximation of said sine of said input angle by a first scale factor, and combine said scaled approximation of said sine of said input angle with said approximation of said cosine of said input angle to yield said cosine of said input angle.
4. The sine/cosine generator according to claim 3, wherein said sine/cosine generator is further configured to scale said approximation of said cosine of said input angle modulo .pi./k by a second scale factor, and combine said scaled approximation of said cosine of said input angle modulo .pi./k with said approximation of said sine of said input angle modulo .pi./k to yield said sine of said input angle modulo .pi./k.
5. The sine/cosine generator according to claim 1, wherein said cosine and sine of said coarse angle are pre-computed values selected from at least one look up table.
6. The sine/cosine generator according to claim 1, wherein said approximations of one minus a cosine of said fine angle and said sine of said fine angle are computed using Taylor Series Approximations.
7. The sine/cosine generator according to claim 1, wherein said first and second most and least significant words are pre-computed values selected from at least one look up table.
8. The sine/cosine generator according to claim 1, wherein said coarse angle is in a quadrant or octant of a two dimensional space.
9. The sine/cosine generator according claim 8, wherein said coarse angle is one of a plurality of angles in said quadrant or octant of a two dimensional space.
10. The sine/cosine generator according to claim 1, wherein said fine angle is bounded by a smallest coarse angle existing in a two dimensional space.
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US12/129,654 US8200728B2 (en) | 2008-05-29 | 2008-05-29 | Sine/cosine generator |
| US12/129,654 | 2008-05-29 | ||
| PCT/US2009/044653 WO2009146284A1 (en) | 2008-05-29 | 2009-05-20 | Sine/cosine generator |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CA2725894A1 CA2725894A1 (en) | 2009-12-03 |
| CA2725894C true CA2725894C (en) | 2013-07-09 |
Family
ID=40974430
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CA2725894A Expired - Fee Related CA2725894C (en) | 2008-05-29 | 2009-05-20 | Sine/cosine generator |
Country Status (7)
| Country | Link |
|---|---|
| US (1) | US8200728B2 (en) |
| EP (1) | EP2304546B1 (en) |
| JP (1) | JP5236804B2 (en) |
| KR (1) | KR101226402B1 (en) |
| CA (1) | CA2725894C (en) |
| TW (1) | TW200951674A (en) |
| WO (1) | WO2009146284A1 (en) |
Families Citing this family (87)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8620980B1 (en) | 2005-09-27 | 2013-12-31 | Altera Corporation | Programmable device with specialized multiplier blocks |
| US8266198B2 (en) | 2006-02-09 | 2012-09-11 | Altera Corporation | Specialized processing block for programmable logic device |
| US8266199B2 (en) | 2006-02-09 | 2012-09-11 | Altera Corporation | Specialized processing block for programmable logic device |
| US8301681B1 (en) | 2006-02-09 | 2012-10-30 | Altera Corporation | Specialized processing block for programmable logic device |
| US8386550B1 (en) | 2006-09-20 | 2013-02-26 | Altera Corporation | Method for configuring a finite impulse response filter in a programmable logic device |
| US8386553B1 (en) | 2006-12-05 | 2013-02-26 | Altera Corporation | Large multiplier for programmable logic device |
| US7930336B2 (en) | 2006-12-05 | 2011-04-19 | Altera Corporation | Large multiplier for programmable logic device |
| US8650231B1 (en) | 2007-01-22 | 2014-02-11 | Altera Corporation | Configuring floating point operations in a programmable device |
| US8312551B2 (en) | 2007-02-15 | 2012-11-13 | Harris Corporation | Low level sequence as an anti-tamper Mechanism |
| US8645450B1 (en) | 2007-03-02 | 2014-02-04 | Altera Corporation | Multiplier-accumulator circuitry and methods |
| US7937427B2 (en) * | 2007-04-19 | 2011-05-03 | Harris Corporation | Digital generation of a chaotic numerical sequence |
| US8611530B2 (en) | 2007-05-22 | 2013-12-17 | Harris Corporation | Encryption via induced unweighted errors |
| US7921145B2 (en) * | 2007-05-22 | 2011-04-05 | Harris Corporation | Extending a repetition period of a random sequence |
| US7995757B2 (en) | 2007-05-31 | 2011-08-09 | Harris Corporation | Closed galois field combination |
| US7970809B2 (en) * | 2007-06-07 | 2011-06-28 | Harris Corporation | Mixed radix conversion with a priori defined statistical artifacts |
| US7962540B2 (en) * | 2007-06-07 | 2011-06-14 | Harris Corporation | Mixed radix number generator with chosen statistical artifacts |
| US7974413B2 (en) * | 2007-06-07 | 2011-07-05 | Harris Corporation | Spread spectrum communications system and method utilizing chaotic sequence |
| US8005221B2 (en) * | 2007-08-01 | 2011-08-23 | Harris Corporation | Chaotic spread spectrum communications system receiver |
| US7995749B2 (en) * | 2007-10-30 | 2011-08-09 | Harris Corporation | Cryptographic system configured for extending a repetition period of a random sequence |
| US8180055B2 (en) | 2008-02-05 | 2012-05-15 | Harris Corporation | Cryptographic system incorporating a digitally generated chaotic numerical sequence |
| US8363830B2 (en) | 2008-02-07 | 2013-01-29 | Harris Corporation | Cryptographic system configured to perform a mixed radix conversion with a priori defined statistical artifacts |
| US8959137B1 (en) | 2008-02-20 | 2015-02-17 | Altera Corporation | Implementing large multipliers in a programmable integrated circuit device |
| US8040937B2 (en) * | 2008-03-26 | 2011-10-18 | Harris Corporation | Selective noise cancellation of a spread spectrum signal |
| US8139764B2 (en) | 2008-05-06 | 2012-03-20 | Harris Corporation | Closed galois field cryptographic system |
| US8320557B2 (en) | 2008-05-08 | 2012-11-27 | Harris Corporation | Cryptographic system including a mixed radix number generator with chosen statistical artifacts |
| US8145692B2 (en) | 2008-05-29 | 2012-03-27 | Harris Corporation | Digital generation of an accelerated or decelerated chaotic numerical sequence |
| US8064552B2 (en) | 2008-06-02 | 2011-11-22 | Harris Corporation | Adaptive correlation |
| US8068571B2 (en) | 2008-06-12 | 2011-11-29 | Harris Corporation | Featureless coherent chaotic amplitude modulation |
| US8594250B2 (en) * | 2008-07-25 | 2013-11-26 | Qualcomm Incorporated | Apparatus and methods for computing constant amplitude zero auto-correlation sequences |
| US8325702B2 (en) * | 2008-08-29 | 2012-12-04 | Harris Corporation | Multi-tier ad-hoc network in which at least two types of non-interfering waveforms are communicated during a timeslot |
| US8165065B2 (en) * | 2008-10-09 | 2012-04-24 | Harris Corporation | Ad-hoc network acquisition using chaotic sequence spread waveform |
| US8307023B1 (en) | 2008-10-10 | 2012-11-06 | Altera Corporation | DSP block for implementing large multiplier on a programmable integrated circuit device |
| US8351484B2 (en) | 2008-12-29 | 2013-01-08 | Harris Corporation | Communications system employing chaotic spreading codes with static offsets |
| US8406276B2 (en) | 2008-12-29 | 2013-03-26 | Harris Corporation | Communications system employing orthogonal chaotic spreading codes |
| US8392492B2 (en) * | 2009-01-19 | 2013-03-05 | Himax Media Solutions, Inc. | Sine/cosine value generation apparatus |
| US8706790B1 (en) | 2009-03-03 | 2014-04-22 | Altera Corporation | Implementing mixed-precision floating-point operations in a programmable integrated circuit device |
| US8645449B1 (en) | 2009-03-03 | 2014-02-04 | Altera Corporation | Combined floating point adder and subtractor |
| US8468192B1 (en) | 2009-03-03 | 2013-06-18 | Altera Corporation | Implementing multipliers in a programmable integrated circuit device |
| US8457077B2 (en) | 2009-03-03 | 2013-06-04 | Harris Corporation | Communications system employing orthogonal chaotic spreading codes |
| US8428102B2 (en) | 2009-06-08 | 2013-04-23 | Harris Corporation | Continuous time chaos dithering |
| US8509284B2 (en) | 2009-06-08 | 2013-08-13 | Harris Corporation | Symbol duration dithering for secured chaotic communications |
| US8428103B2 (en) | 2009-06-10 | 2013-04-23 | Harris Corporation | Discrete time chaos dithering |
| US8340295B2 (en) * | 2009-07-01 | 2012-12-25 | Harris Corporation | High-speed cryptographic system using chaotic sequences |
| US8363700B2 (en) | 2009-07-01 | 2013-01-29 | Harris Corporation | Rake receiver for spread spectrum chaotic communications systems |
| US8428104B2 (en) | 2009-07-01 | 2013-04-23 | Harris Corporation | Permission-based multiple access communications systems |
| US8406352B2 (en) | 2009-07-01 | 2013-03-26 | Harris Corporation | Symbol estimation for chaotic spread spectrum signal |
| US8379689B2 (en) | 2009-07-01 | 2013-02-19 | Harris Corporation | Anti-jam communications having selectively variable peak-to-average power ratio including a chaotic constant amplitude zero autocorrelation waveform |
| US8385385B2 (en) | 2009-07-01 | 2013-02-26 | Harris Corporation | Permission-based secure multiple access communication systems |
| US8369376B2 (en) | 2009-07-01 | 2013-02-05 | Harris Corporation | Bit error rate reduction in chaotic communications |
| US8848909B2 (en) * | 2009-07-22 | 2014-09-30 | Harris Corporation | Permission-based TDMA chaotic communication systems |
| US8369377B2 (en) | 2009-07-22 | 2013-02-05 | Harris Corporation | Adaptive link communications using adaptive chaotic spread waveform |
| US8650236B1 (en) | 2009-08-04 | 2014-02-11 | Altera Corporation | High-rate interpolation or decimation filter in integrated circuit device |
| US8412756B1 (en) | 2009-09-11 | 2013-04-02 | Altera Corporation | Multi-operand floating point operations in a programmable integrated circuit device |
| US8396914B1 (en) | 2009-09-11 | 2013-03-12 | Altera Corporation | Matrix decomposition in an integrated circuit device |
| US8539016B1 (en) | 2010-02-09 | 2013-09-17 | Altera Corporation | QR decomposition in an integrated circuit device |
| US8601044B2 (en) | 2010-03-02 | 2013-12-03 | Altera Corporation | Discrete Fourier Transform in an integrated circuit device |
| US8484265B1 (en) | 2010-03-04 | 2013-07-09 | Altera Corporation | Angular range reduction in an integrated circuit device |
| US8345725B2 (en) | 2010-03-11 | 2013-01-01 | Harris Corporation | Hidden Markov Model detection for spread spectrum waveforms |
| US8510354B1 (en) | 2010-03-12 | 2013-08-13 | Altera Corporation | Calculation of trigonometric functions in an integrated circuit device |
| US8539014B2 (en) | 2010-03-25 | 2013-09-17 | Altera Corporation | Solving linear matrices in an integrated circuit device |
| US8862650B2 (en) | 2010-06-25 | 2014-10-14 | Altera Corporation | Calculation of trigonometric functions in an integrated circuit device |
| US8589463B2 (en) | 2010-06-25 | 2013-11-19 | Altera Corporation | Calculation of trigonometric functions in an integrated circuit device |
| US8577951B1 (en) | 2010-08-19 | 2013-11-05 | Altera Corporation | Matrix operations in an integrated circuit device |
| US8645451B2 (en) | 2011-03-10 | 2014-02-04 | Altera Corporation | Double-clocked specialized processing block in an integrated circuit device |
| CN102156626B (en) * | 2011-04-27 | 2012-10-10 | 哈尔滨工业大学 | Sine-cosine function IP core capable of reconfiguring spaceborne computer |
| US9600278B1 (en) | 2011-05-09 | 2017-03-21 | Altera Corporation | Programmable device using fixed and configurable logic to implement recursive trees |
| US8654819B2 (en) | 2011-06-22 | 2014-02-18 | Harris Corporation | Systems and methods for pulse rotation modulation encoding and decoding |
| US8812576B1 (en) | 2011-09-12 | 2014-08-19 | Altera Corporation | QR decomposition in an integrated circuit device |
| US8949298B1 (en) | 2011-09-16 | 2015-02-03 | Altera Corporation | Computing floating-point polynomials in an integrated circuit device |
| US9053045B1 (en) | 2011-09-16 | 2015-06-09 | Altera Corporation | Computing floating-point polynomials in an integrated circuit device |
| US8762443B1 (en) | 2011-11-15 | 2014-06-24 | Altera Corporation | Matrix operations in an integrated circuit device |
| US8543634B1 (en) | 2012-03-30 | 2013-09-24 | Altera Corporation | Specialized processing block for programmable integrated circuit device |
| US9098332B1 (en) | 2012-06-01 | 2015-08-04 | Altera Corporation | Specialized processing block with fixed- and floating-point structures |
| US8996600B1 (en) | 2012-08-03 | 2015-03-31 | Altera Corporation | Specialized processing block for implementing floating-point multiplier with subnormal operation support |
| US9207909B1 (en) | 2012-11-26 | 2015-12-08 | Altera Corporation | Polynomial calculations optimized for programmable integrated circuit device structures |
| US9189200B1 (en) | 2013-03-14 | 2015-11-17 | Altera Corporation | Multiple-precision processing block in a programmable integrated circuit device |
| US9348795B1 (en) | 2013-07-03 | 2016-05-24 | Altera Corporation | Programmable device using fixed and configurable logic to implement floating-point rounding |
| US9811503B1 (en) * | 2015-01-28 | 2017-11-07 | Altera Corporation | Methods for implementing arithmetic functions with user-defined input and output formats |
| US9684488B2 (en) | 2015-03-26 | 2017-06-20 | Altera Corporation | Combined adder and pre-adder for high-radix multiplier circuit |
| US9875084B2 (en) | 2016-04-28 | 2018-01-23 | Vivante Corporation | Calculating trigonometric functions using a four input dot product circuit |
| US10942706B2 (en) | 2017-05-05 | 2021-03-09 | Intel Corporation | Implementation of floating-point trigonometric functions in an integrated circuit device |
| US10997510B1 (en) * | 2018-05-22 | 2021-05-04 | Marvell Asia Pte, Ltd. | Architecture to support tanh and sigmoid operations for inference acceleration in machine learning |
| US10929779B1 (en) | 2018-05-22 | 2021-02-23 | Marvell Asia Pte, Ltd. | Architecture to support synchronization between core and inference engine for machine learning |
| US10929778B1 (en) | 2018-05-22 | 2021-02-23 | Marvell Asia Pte, Ltd. | Address interleaving for machine learning |
| US10740432B1 (en) * | 2018-12-13 | 2020-08-11 | Amazon Technologies, Inc. | Hardware implementation of mathematical functions |
| JP7262028B2 (en) * | 2019-07-01 | 2023-04-21 | ラピスセミコンダクタ株式会社 | trigonometric function unit |
| US20220317972A1 (en) * | 2021-03-30 | 2022-10-06 | Micron Technology, Inc. | Hardware for concurrent sine and cosine determination |
Family Cites Families (58)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR1501059A (en) | 1966-09-26 | 1967-11-10 | Csf | New key generator |
| GB1198263A (en) * | 1967-06-06 | 1970-07-08 | Nat Res Dev | Improvements in or relating to Digital Differential Analysers |
| US4077063A (en) * | 1976-08-11 | 1978-02-28 | The Singer Company | Apparatus for rapidly determining the trigonometric functions of an input angle |
| US4486846A (en) * | 1982-06-28 | 1984-12-04 | Motorola Inc. | Numerically controlled oscillator using quadrant replication and function decomposition |
| US4646326A (en) * | 1983-10-20 | 1987-02-24 | Motorola Inc. | QAM modulator circuit |
| US4703507A (en) * | 1984-04-05 | 1987-10-27 | Holden Thomas W | Noise reduction system |
| US4855946A (en) * | 1987-10-22 | 1989-08-08 | Motorola, Inc. | Reduced size phase-to-amplitude converter in a numerically controlled oscillator |
| US5077793A (en) * | 1989-09-29 | 1991-12-31 | The Boeing Company | Residue number encryption and decryption system |
| US5113361A (en) * | 1990-08-02 | 1992-05-12 | Harris Corporation | Sin/cos generator implementation |
| US5319735A (en) * | 1991-12-17 | 1994-06-07 | Bolt Beranek And Newman Inc. | Embedded signalling |
| US5276633A (en) * | 1992-08-14 | 1994-01-04 | Harris Corporation | Sine/cosine generator and method |
| JPH06132731A (en) * | 1992-10-20 | 1994-05-13 | Fujitsu Ltd | Cosine wave generation circuit |
| US6614914B1 (en) * | 1995-05-08 | 2003-09-02 | Digimarc Corporation | Watermark embedder and reader |
| JPH07175636A (en) * | 1993-12-16 | 1995-07-14 | Fuji Film Micro Device Kk | Digital trigonometric function generating circuit |
| US5467294A (en) * | 1994-03-09 | 1995-11-14 | Hu; Vince | High speed, low power direct digital synthesizer |
| US5937000A (en) * | 1995-09-06 | 1999-08-10 | Solana Technology Development Corporation | Method and apparatus for embedding auxiliary data in a primary data signal |
| WO1997044935A1 (en) * | 1996-05-20 | 1997-11-27 | Philips Electronics N.V. | Cryptographic method and apparatus for non-linearly merging a data block and a key |
| AU3567997A (en) * | 1996-07-05 | 1998-02-02 | Paulo Correa | Controller-based radio frequency amplifier module and method |
| US5963460A (en) | 1996-12-17 | 1999-10-05 | Metaflow Technologies, Inc. | Apparatus for computing transcendental functions quickly |
| WO1998059458A1 (en) * | 1997-06-23 | 1998-12-30 | The Regents Of The University Of California | Chaotic digital code-division multiple access (cdma) communication systems |
| US6633226B1 (en) * | 1997-08-18 | 2003-10-14 | X-Cyte, Inc. | Frequency hopping spread spectrum passive acoustic wave identification device |
| US6285761B1 (en) | 1998-03-04 | 2001-09-04 | Lucent Technologies, Inc. | Method for generating pseudo-random numbers |
| US5924980A (en) * | 1998-03-11 | 1999-07-20 | Siemens Corporate Research, Inc. | Method and apparatus for adaptively reducing the level of noise in an acquired signal |
| US6363104B1 (en) * | 1998-10-02 | 2002-03-26 | Ericsson Inc. | Method and apparatus for interference cancellation in a rake receiver |
| US6570909B1 (en) * | 1999-07-09 | 2003-05-27 | Nokia Mobile Phones | Interference suppression in a CDMA receiver |
| US20020099746A1 (en) * | 1999-07-26 | 2002-07-25 | Tie Teck Sing | T-sequence apparatus and method for general deterministic polynomial-time primality testing and composite factoring |
| US6909785B1 (en) | 1999-11-11 | 2005-06-21 | Qualcomm, Inc. | Method and apparatus for efficient irregular synchronization of a stream cipher |
| US7010559B2 (en) * | 2000-11-14 | 2006-03-07 | Parkervision, Inc. | Method and apparatus for a parallel correlator and applications thereof |
| JP4188571B2 (en) * | 2001-03-30 | 2008-11-26 | 株式会社日立製作所 | Arithmetic method of information processing apparatus and tamper resistant arithmetic disturbance implementation method |
| US7076065B2 (en) * | 2001-05-11 | 2006-07-11 | Lockheed Martin Corporation | Chaotic privacy system and method |
| US7269198B1 (en) * | 2001-11-19 | 2007-09-11 | Bbn Technologies Corp. | Systems and methods for beaconing in wireless networks with low probability of detection |
| JP2003218835A (en) * | 2002-01-18 | 2003-07-31 | Mitsubishi Electric Corp | Spread spectrum transmitter and spread spectrum receiver |
| FR2837331B1 (en) * | 2002-03-13 | 2004-06-18 | Canon Kk | METHOD OF INTERLEAVING A BINARY SEQUENCE |
| KR20030029064A (en) * | 2003-02-10 | 2003-04-11 | 구시경 | SINE/COSINE waveform generator without LUT(Look Up Table) by periodical data shifting |
| US7328228B2 (en) * | 2003-09-02 | 2008-02-05 | Sap Aktiengesellschaft | Mapping pseudo-random numbers to predefined number ranges |
| US7512645B2 (en) * | 2004-03-19 | 2009-03-31 | Texas Instruments Incorporated | System and method for generating pseudorandom numbers |
| US7078981B2 (en) * | 2004-07-27 | 2006-07-18 | Lucent Technologies Inc. | 16 QAM modulator and method of 16 QAM modulation |
| WO2006110954A1 (en) | 2005-04-20 | 2006-10-26 | Synaptic Laboratories Limited | Process of and apparatus for counting |
| ITVA20050027A1 (en) * | 2005-05-03 | 2006-11-04 | St Microelectronics Srl | METHOD OF GENERATION OF SUCCESSIONS OF NUMBERS OR BIT PSEUDO CASUALI |
| US8312551B2 (en) * | 2007-02-15 | 2012-11-13 | Harris Corporation | Low level sequence as an anti-tamper Mechanism |
| WO2008099367A2 (en) | 2007-02-15 | 2008-08-21 | Koninklijke Philips Electronics N.V. | Coordination in wireless networks having devices with different physical layer transmission schemes |
| US7937427B2 (en) * | 2007-04-19 | 2011-05-03 | Harris Corporation | Digital generation of a chaotic numerical sequence |
| US7921145B2 (en) * | 2007-05-22 | 2011-04-05 | Harris Corporation | Extending a repetition period of a random sequence |
| US8611530B2 (en) * | 2007-05-22 | 2013-12-17 | Harris Corporation | Encryption via induced unweighted errors |
| US7995757B2 (en) * | 2007-05-31 | 2011-08-09 | Harris Corporation | Closed galois field combination |
| US7970809B2 (en) * | 2007-06-07 | 2011-06-28 | Harris Corporation | Mixed radix conversion with a priori defined statistical artifacts |
| US7974413B2 (en) * | 2007-06-07 | 2011-07-05 | Harris Corporation | Spread spectrum communications system and method utilizing chaotic sequence |
| US7962540B2 (en) * | 2007-06-07 | 2011-06-14 | Harris Corporation | Mixed radix number generator with chosen statistical artifacts |
| US8005221B2 (en) * | 2007-08-01 | 2011-08-23 | Harris Corporation | Chaotic spread spectrum communications system receiver |
| US7995749B2 (en) * | 2007-10-30 | 2011-08-09 | Harris Corporation | Cryptographic system configured for extending a repetition period of a random sequence |
| US8180055B2 (en) * | 2008-02-05 | 2012-05-15 | Harris Corporation | Cryptographic system incorporating a digitally generated chaotic numerical sequence |
| US8363830B2 (en) * | 2008-02-07 | 2013-01-29 | Harris Corporation | Cryptographic system configured to perform a mixed radix conversion with a priori defined statistical artifacts |
| US8040937B2 (en) * | 2008-03-26 | 2011-10-18 | Harris Corporation | Selective noise cancellation of a spread spectrum signal |
| US8139764B2 (en) * | 2008-05-06 | 2012-03-20 | Harris Corporation | Closed galois field cryptographic system |
| US8320557B2 (en) * | 2008-05-08 | 2012-11-27 | Harris Corporation | Cryptographic system including a mixed radix number generator with chosen statistical artifacts |
| US8145692B2 (en) * | 2008-05-29 | 2012-03-27 | Harris Corporation | Digital generation of an accelerated or decelerated chaotic numerical sequence |
| US8064552B2 (en) * | 2008-06-02 | 2011-11-22 | Harris Corporation | Adaptive correlation |
| US8068571B2 (en) * | 2008-06-12 | 2011-11-29 | Harris Corporation | Featureless coherent chaotic amplitude modulation |
-
2008
- 2008-05-29 US US12/129,654 patent/US8200728B2/en active Active
-
2009
- 2009-05-20 JP JP2011511719A patent/JP5236804B2/en not_active Expired - Fee Related
- 2009-05-20 KR KR1020107029516A patent/KR101226402B1/en active IP Right Grant
- 2009-05-20 WO PCT/US2009/044653 patent/WO2009146284A1/en active Application Filing
- 2009-05-20 CA CA2725894A patent/CA2725894C/en not_active Expired - Fee Related
- 2009-05-20 EP EP09755646.8A patent/EP2304546B1/en active Active
- 2009-05-22 TW TW098117128A patent/TW200951674A/en unknown
Also Published As
| Publication number | Publication date |
|---|---|
| US8200728B2 (en) | 2012-06-12 |
| TW200951674A (en) | 2009-12-16 |
| CA2725894A1 (en) | 2009-12-03 |
| KR20110015026A (en) | 2011-02-14 |
| JP5236804B2 (en) | 2013-07-17 |
| US20090300088A1 (en) | 2009-12-03 |
| EP2304546B1 (en) | 2014-04-02 |
| JP2011522330A (en) | 2011-07-28 |
| KR101226402B1 (en) | 2013-01-24 |
| EP2304546A1 (en) | 2011-04-06 |
| WO2009146284A1 (en) | 2009-12-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CA2725894C (en) | Sine/cosine generator | |
| US5737253A (en) | Method and apparatus for direct digital frequency synthesizer | |
| JPH03204720A (en) | Elementary function arithmetic unit | |
| US4747067A (en) | Apparatus and method for approximating the magnitude of a complex number | |
| US9151842B2 (en) | Method and apparatus for time of flight sensor 2-dimensional and 3-dimensional map generation | |
| EP0288720B1 (en) | Curve generation in a display system | |
| Fu et al. | A high-speed processor for digital sine/cosine generation and angle rotation | |
| Deng et al. | High-speed parameterisable Hough transform using reconfigurable hardware | |
| KR100959215B1 (en) | Speed detection device | |
| Mokhtar et al. | Implementation of trigonometric function using CORDIC algorithms | |
| Zyuzina et al. | Monotone approximation of a scalar conservation law based on the CABARET scheme in the case of a sign-changing characteristic field | |
| CN108182050A (en) | Utilize four input dot product circuit counting trigonometric functions | |
| JP2004172794A (en) | Composite trigonometric function computing element and phase rotation computing element | |
| KR20220170751A (en) | Cordic computation of sin/cos using combined approach in associative memory | |
| Alnafutchy | Calculation of trigonometric functions using cordic algorithm | |
| Fousse | Accurate multiple-precision Gauss-Legendre quadrature | |
| JPH08241187A (en) | Booth multiplier for trigonometric function | |
| Sharmila | VLSI Implementation of Coordinate Rotation Based Design Methodology using Verilog HDL | |
| CN114063701A (en) | Low-delay-based accurate frequency output method and system | |
| Akhil | Cordic Algorithm based Design and Implementation of Trigonometric Co-Processor | |
| Duhan | Design and Implementation of Floating Point Sine and Cosine Functions Using Xilinx ISE for Analog Circuits. | |
| Singh et al. | Analysis of simple CORDIC algorithm using MATLAB | |
| Sanchez et al. | Architecture implementation of an improved decimal CORDIC method | |
| JP6308845B2 (en) | Arithmetic apparatus, arithmetic method, and program | |
| Ijaz | A Digital Integrated Inertial Navigation System for Aerial Vehicles |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| EEER | Examination request | ||
| MKLA | Lapsed |
Effective date: 20170523 |