Embedded Software and Motor Control Libraries for PXR40xx

[- π, π ), the input f16In must be multiplied by π. So, the fixed point fractional
implementation of the GFLIB_Sin_F16 function, using 9th order Taylor approximation,
is given as follows:
Equation GFLIB_Cos_Eq3
The 7th order polynomial approximation of the sine function has a good accuracy in the
range [- π/2, π/2) of the argument, but in wider ranges the calculation error quickly
increases. To minimize the error without having to use a higher order polynomial, the
symmetry of the sine function sin(x) = sin( π - x) is utilized. Therefore, the input
argument is transferred to be always in the range [- π/2, π/2) and the Taylor polynomial
is calculated only in the range of the argument [- π/2, π/2).
To make calculations more precise, the given argument value f16In (that is to be
transferred into the range [-0.5,0.5) due to the sine function symmetry) is shifted by 1 bit
to the left (multiplied by 2). Then, the value of f16In 2 , used in the calculations, is in the
range [-1, 1) instead of [-0.25, 0.25]. Shifting the input value by 1 bit to the left will
increase the accuracy of the calculated sin( π * f16In) function. Implementing such a
scale on the approximation function described by equation GFLIB_Cos_Eq2, results in
the following:
Equation GFLIB_Cos_Eq4
Equation GFLIB_Cos_Eq3 can be further rewritten into the following form:
Equation GFLIB_Cos_Eq5
Chapter 4 API References
where a 1...a 5 are coefficients of the approximation polynomial, which are calculated as
follows (represented as 16-bit signed fractional numbers):
Embedded Software and Motor Control Libraries for PXR40xx, Rev. 1.0
Freescale Semiconductor, Inc. 323