Zolotarev polynomials

In mathematics, Zolotarev polynomials are polynomials used in approximation theory. They are sometimes used as an alternative to the Chebyshev polynomials where accuracy of approximation near the origin is of less importance. Zolotarev polynomials differ from the Chebyshev polynomials in that two of the coefficients are fixed in advance rather than allowed to take on any value. The Chebyshev polynomials of the first kind are a special case of Zolotarev polynomials. These polynomials were introduced by Russian mathematician Yegor Ivanovich Zolotarev in 1868.

Zolotarev polynomials of degree ${\displaystyle n}$ in ${\displaystyle x}$ are of the form

${\displaystyle Z_{n}(x,\sigma )=x^{n}-\sigma x^{n-1}+\cdots +a_{k}x^{k}+\cdots +a_{0}\ ,}$

where ${\displaystyle \sigma }$ is a prescribed value for ${\displaystyle a_{n-1}}$ and the ${\displaystyle a_{k}\in \mathbb {R} }$ are otherwise chosen such that the deviation of ${\displaystyle Z_{n}(x)}$ from zero is minimum in the interval ${\displaystyle [-1,1]}$.^[1]

A subset of Zolotarev polynomials can be expressed in terms of Chebyshev polynomials of the first kind, ${\displaystyle T_{n}(x)}$. For

${\displaystyle 0\leq \sigma \leq {\dfrac {1}{n}}\tan ^{2}{\dfrac {\pi }{2n}}}$

then

${\displaystyle Z_{n}(x,\sigma )=(1+\sigma )^{n}T_{n}\left({\frac {x-\sigma }{1+\sigma }}\right)\ .}$

For values of ${\displaystyle \sigma }$ greater than the maximum of this range, Zolotarev polynomials can be expressed in terms of elliptic functions. For ${\displaystyle \sigma =0}$, the Zolotarev polynomial is identical to the equivalent Chebyshev polynomial. For negative values of ${\displaystyle \sigma }$, the polynomial can be found from the polynomial of the positive value,^[2]

${\displaystyle Z_{n}(x,-\sigma )=(-1)^{n}Z_{n}(-x,\sigma )\ .}$

The Zolotarev polynomial can be expanded into a sum of Chebyshev polynomials using the relationship^[3]

${\displaystyle Z_{n}(x)=\sum _{k=0}^{n}a_{k}T_{k}(x)\ .}$

The original solution to the approximation problem given by Zolotarev was in terms of Jacobi elliptic functions. Zolotarev gave the general solution where the number of zeroes to the left of the peak value (${\displaystyle q}$) in the interval ${\displaystyle [-1,1]}$ is not equal to the number of zeroes to the right of this peak (${\displaystyle p}$). The degree of the polynomial is ${\displaystyle n=p+q}$. For many applications, ${\displaystyle p=q}$ is used and then only ${\displaystyle n}$ need be considered. The general Zolotarev polynomials are defined as^[5]

${\displaystyle Z_{n}(x|\kappa )={\frac {(-1)^{p}}{2}}\left[\left({\dfrac {H(u-v)}{H(u+v)}}\right)^{n}+\left({\dfrac {H(u+v)}{H(u-v)}}\right)^{n}\right]}$

where

${\displaystyle u=F\left(\left.\operatorname {sn} \left(\left.v\right|\kappa \right){\sqrt {\dfrac {1+x}{x+2\operatorname {sn} ^{2}\left(\left.v\right|\kappa \right)-1}}}\right|\kappa \right)}$

${\displaystyle v={\dfrac {p}{n}}K(\kappa )}$

${\displaystyle H(\varphi )}$ is the Jacobi eta function

${\displaystyle F(\varphi |\kappa )}$ is the incomplete elliptic integral of the first kind

${\displaystyle K(\kappa )}$ is the quarter-wave complete elliptic integral of the first kind. That is, ${\displaystyle K(\kappa )=F\left(\left.{\frac {\pi }{2}}\right|\kappa \right)}$^[6]

${\displaystyle \kappa }$ is the Jacobi elliptic modulus

${\displaystyle \operatorname {sn} (\varphi |\kappa )}$ is the Jacobi elliptic sine.

The variation of the function within the interval [−1,1] is equiripple except for one peak which is larger than the rest. The position and width of this peak can be set independently. The position of the peak is given by^[7]

${\displaystyle x_{\text{max}}=1-2\operatorname {sn} ^{2}(v|\kappa )+2{\dfrac {\operatorname {sn} (v|\kappa )\operatorname {cn} (v|\kappa )}{\operatorname {dn} (v|\kappa )}}Z(v|\kappa )}$

where

${\displaystyle \operatorname {cn} (\varphi |\kappa )}$ is the Jacobi elliptic cosine

${\displaystyle \operatorname {dn} (\varphi |\kappa )}$ is the Jacobi delta amplitude

${\displaystyle Z(\varphi |\kappa )}$ is the Jacobi zeta function

${\displaystyle v}$ is as defined above.

The height of the peak is given by^[8]

${\displaystyle Z_{n}(x_{\text{max}}|\kappa )=\cosh 2n{\bigl (}\sigma _{\text{max}}Z(v|\kappa )-\varPi (\sigma _{\text{max}},v|\kappa ){\bigr )}}$

where

${\displaystyle \varPi (\phi _{1},\phi _{2}|\kappa )}$ is the incomplete elliptic integral of the third kind

${\displaystyle \sigma _{\text{max}}=F\left(\left.\sin ^{-1}\left({\dfrac {1}{\kappa \operatorname {sn} (v|\kappa )}}{\sqrt {\dfrac {x_{\text{max}}-x_{\mathrm {L} }}{x_{\text{max}}+1}}}\right)\right|\kappa \right)}$

${\displaystyle x_{\mathrm {L} }}$ is the position on the left limb of the peak which is the same height as the equiripple peaks.

The Jacobi eta function can be defined in terms of a Jacobi auxiliary theta function,^[9]

${\displaystyle H(\varphi |\kappa )=\theta _{1}(a|b)}$

where,

${\displaystyle a={\frac {\pi \varphi }{2K'(\kappa )}}}$

${\displaystyle b=\exp \left(-{\frac {\pi K'(\kappa )}{K(\kappa )}}\right)}$

${\displaystyle K'(\kappa )=K({\sqrt {1-\kappa ^{2}}})\ .}$^[10]

The polynomials were introduced by Yegor Ivanovich Zolotarev in 1868 as a means of uniformly approximating polynomials of degree ${\displaystyle x^{n+1}}$ on the interval [−1,1]. Pafnuty Chebyshev had shown in 1858 that ${\displaystyle x^{n+1}}$ could be approximated in this interval with a polynomial of degree at most ${\displaystyle n}$ with an error of ${\displaystyle 2^{-n}}$. In 1868, Zolotarev showed that ${\displaystyle x^{n+1}-\sigma x^{n}}$ could be approximated with a polynomial of degree at most ${\displaystyle n-1}$, two degrees lower. The error in Zolotarev's method is given by,^[11]

${\displaystyle 2^{-n}\left({\dfrac {1+\sigma }{1+n}}\right)^{n+1}\ .}$

The procedure was further developed by Naum Achieser in 1956.^[12]

Zolotarev polynomials are used in the design of Achieser-Zolotarev filters. They were first used in this role in 1970 by Ralph Levy in the design of microwave waveguide filters.^[13] Achieser-Zolotarev filters are similar to Chebyshev filters in that they have an equal ripple attenuation through the passband, except that the attenuation exceeds the preset ripple for the peak closest to the origin.^[14]

Zolotarev polynomials can be used to synthesise the radiation patterns of linear antenna arrays, first suggested by D.A. McNamara in 1985. The work was based on the filter application with beam angle used as the variable instead of frequency. The Zolotarev beam pattern has equal-level sidelobes.^[15]

• Achieser, Naum, Hymnan, C.J. (trans), Theory of Approximation, New York: Frederick Ungar Publishing, 1956. Dover reprint 2013 ISBN 0486495434.
• Beebe, Nelson H.F., The Mathematical-Function Computation Handbook, Springer, 2017 ISBN 978-3-319-64110-2.
• Cameron, Richard J.; Kudsia, Chandra M.; Mansour, Raafat R., Microwave Filters for Communication Systems, John Wiley & Sons, 2018 ISBN 1118274342.
• Hansen, Robert C., Phased Array Antennas, Wiley, 2009 ISBN 0470529172.
• McNamara, D.A., "Optimum monopulse linear array excitations using Zolotarev Polynomials", Electron, vol. 21, iss. 16, pp. 681–682, August 1985.
• Newman, D.J., Reddy, A.R., "Rational approximations to ${\displaystyle x^{n}}$ II", Canadian Journal of Mathematics, vol. 32, no. 2, pp. 310–316, April 1980.
• Pinkus, Allan, "Zolotarev polynomials", in, Hazewinkel, Michiel (ed), Encyclopaedia of Mathematics, Supplement III, Springer Science & Business Media, 2001 ISBN 1402001983.
• Vlček, Miroslav, Unbehauen, Rolf, "Zolotarev polynomials and optimal FIR filters", IEEE Transactions on Signal Processing, vol. 47, iss. 3, pp. 717–730, March 1999 (corrections July 2000).
• Zahradnik, Pavel; Vlček, Miroslav, "Analytical design of 2-D narrow bandstop FIR filters", pp. 56–63 in, Computational Science — ICCS 2004: Proceedings of the 4th International Conference, Bubak, Marian; van Albada, Geert D.; Sloot, Peter M.A.; Dongarra, Jack (eds), Springer Science & Business Media, 2004 ISBN 3540221298.