[PDF] Efficient Approximations for the Arctangent Function T



Previous PDF Next PDF
















[PDF] fonction circulaire réciproque cours

[PDF] limite de arctan

[PDF] limite arctan en 0

[PDF] le pouvoir du peuple par le peuple pour le peuple

[PDF] fonctions trigonométriques réciproques pdf

[PDF] shlomo sand livres

[PDF] le peuple est il souverain dissertation

[PDF] exercices corrigés fonction arctangente

[PDF] fonction circulatoire définition

[PDF] comment la terre d'israël fut inventée pdf

[PDF] origine des juifs d'israel

[PDF] appareil circulatoire cours

[PDF] système circulatoire

[PDF] comment la terre d'israël fut inventée

[PDF] appareil circulatoire schéma

Efficient Approximations for the Arctangent Function T

Efficient Approximations

for the Arctangent Function

IEEE SIGNAL PROCESSING MAGAZINE[108]MAY 2006

dspTIPS&TRICKS

Sreeraman Rajan, Sichun Wang,

Robert Inkol, and Alain Joyal

1053-5888/06/$20.00©2006IEEE

T his article provides several efficient approximations for the arctangent function using Lagrange interpolation and minimax optimization techniques. These approximations are particularly useful when processing power, memory, and power consumption are important issues. In addition to com- paring the errors and the computational workload of these approximations, we also extend them to all four quadrants.

ARCTANGENT APPROXIMATIONS

The evaluation of the arctangent func-

tion is commonly encountered in real- time signal processing applications.

Numerous algorithms are available to

implement the arctangent function when computational cost is unimpor- tant. The most direct solution is based on the Taylor series expansion on [

Š1,

1]. However, this series converges

slowly for arguments close to one and hence is inefficient.Iterative algorithms, such as the coordinated rotation digital computer (CORDIC) algorithm (requiring only shifts and add operations), have been successfully used to compute trigono- metric functions [1], [2]. However, their sequential nature makes them less attractive when speed is a major concern; attempts to increase speed have been at the expense of additional hardware [3]. Lookup-table-based approaches to the computation of inverse trigonometric functions are very fast but require considerable memory [4], [5].

Polynomial and rational function

approximations that have been pro- posed in the literature [3] are more suitable for numerical coprocessors.

Approximations using polynomials of

large degrees are computationally expensive. Rational approximations are in principle more accurate than polyno- mial approximations for the same num- ber of coefficients. However, therequired division operations are rela- tively complex to implement; iterative techniques, such as those based on

NewtonÕs method, are often used.

The aforementioned approaches are

best suited for applications where pro- cessing power and/or memory are read- ily available. However, for many applications, simpler and more effi- cient ways of evaluating arctan( x) are desirable. This article proposes simple approximations for evaluating the arct- angent function that may be easily implemented in hardware with limited memory and processing power.

PROPOSED METHODOLOGY

AND APPROXIMATIONS

Approximations to the arctangent

function can be obtained using second- and third-order polynomials and simple rational functions. Lagrange interpolation-based and minimax criterion-basedapproaches are used to obtain the polynomial coefficients. The following is our initial derivation of an arctangent approximation using

Lagrange interpolation.

Consider the three points

x0 =Š1, x 1 =0, and x 2 =1. Let (x)=arctan (1+x )/(1Šx), Š1x1. According to the Lagrange interpolation formula, we have (x)(xŠx1 )(xŠx 2 (x 0 Šx 1 )(x 0 Šx 2 )(x 0 (xŠx 0 )(xŠx2 (x 1 Šx 0 )(x 1 Šx 2 )(x 1 (xŠx 0 )(xŠx 1 (x 2 Šx 0 )(x 2 Šx 1 )(x 2)

4(x+1),Š1x1.(1)

It can be verified that

(x)=/4+ arctan( x); hence, arctan(x) can be approximated by the first-order formula [FIG1] Approximation errors using (2).0

Error (rad)

x

1-1-8-4048

× 10

-2

IEEE SIGNAL PROCESSING MAGAZINE[109]MAY 2006

arctan(x)π

4x,Š1x1.(2)

This linear approximation has been used

in [6] for FM demodulation due to its minimal complexity. It requires only a scaling operation by a fixed constant and can be computed in one cycle in most processors. Figure 1 shows the deviation of this linear (first-order) approximation from the arctangent function on the interval [

Š1, 1]. This error function,

given in (3), is antisymmetric with respect to x=0and attains a maximum of about

0.07 rad (4¼) at

x max =±(4/πŠ1) 1/2 (x)=arctan(x)Šπ 4x,

Š1x1.(3)

HereÕs a useful trick to improve the

accuracy of the first-order approximation (2). It can be seen from Figure 1 that the error curve is approximately quadratic on the interval [0, 1] and hence can be approximated by a second-order polyno- mial. Following this strategy and apply- ing LagrangeÕs interpolation formula to (x)using x0 =0, x 1 =x max , and x 2 =1, we obtain the following formula: (x)0.285x(1Š|x|),

Š1x1,(4)

where the odd symmetry of (x)has been applied. A second-order, and more accurate, approximation for arctan( x)quotesdbs_dbs2.pdfusesText_2