Chebyshev rational functions

For the Chebyshev rational functions used in the design of elliptic filters, see Elliptic rational functions.

File:ChebychevRational1.png

Plot of the Chebyshev rational functions for n = 0, 1, 2, 3, 4 for 0.01 ≤ x ≤ 100, log scale.

In mathematics, the Chebyshev rational functions are a sequence of functions which are both rational and orthogonal. They are named after Pafnuty Chebyshev. A rational Chebyshev function of degree n is defined as:

$R_n(x)\stackrel{\mathrm{d}\mathrm{e}\mathrm{f}}{=}{T}_n\left(\frac{x-1}{x+1}\right)$

where T_n(x) is a Chebyshev polynomial of the first kind.

Properties

Many properties can be derived from the properties of the Chebyshev polynomials of the first kind. Other properties are unique to the functions themselves.

Recursion

$R_{n+1}(x)=2\left(\frac{x-1}{x+1}\right)R_n(x)-R_{n-1}(x)\text{for}n\ge 1$

Differential equations

$(x+1{)}^2{R}_n(x)=\frac{1}{n+1}\frac{\mathrm{d}}{\mathrm{d}x}{R}_{n+1}(x)-\frac{1}{n-1}\frac{\mathrm{d}}{\mathrm{d}x}{R}_{n-1}(x)\text{for }n\ge 2$

$(x+1{)}^{2}x\frac{{\mathrm{d}}^2}{\mathrm{d}{x}^2}{R}_n(x)+\frac{(3x+1)(x+1)}{2}\frac{\mathrm{d}}{\mathrm{d}x}{R}_n(x)+n^2{R}_n(x)=0$

Orthogonality

File:ChebychevRational2.png

Plot of the absolute value of the seventh-order (n = 7) Chebyshev rational function for 0.01 ≤ x ≤ 100. Note that there are n zeroes arranged symmetrically about x = 1 and if x₀ is a zero, then ⁠1/x₀⁠ is a zero as well. The maximum value between the zeros is unity. These properties hold for all orders.

Defining:

$\omega (x)\stackrel{\mathrm{d}\mathrm{e}\mathrm{f}}{=}\frac{1}{(x+1)\sqrt{x}}$

The orthogonality of the Chebyshev rational functions may be written:

$\underset{0}{\overset{\mathrm{\infty }}{\int }}{R}_m(x)R_n(x)\omega (x)\mathrm{d}x=\frac{\pi c_n}{2}{\delta }_{nm}$

where c_n = 2 for n = 0 and c_n = 1 for n ≥ 1; δ_nm is the Kronecker delta function.

Expansion of an arbitrary function

For an arbitrary function f(x) ∈ L²
_ω the orthogonality relationship can be used to expand f(x):

$f(x)=\sum _{n=0}^{\mathrm{\infty }}{F}_n{R}_n(x)$

where

$F_n=\frac{2}{c_n\pi }\underset{0}{\overset{\mathrm{\infty }}{\int }}f(x)R_n(x)\omega (x)\mathrm{d}x.$

Particular values

$\begin{array}{cc}{R}_0(x)& =1\\ R_1(x)& =\frac{x-1}{x+1}\\ R_2(x)& =\frac{x^2-6x+1}{(x+1{)}^2}\\ R_3(x)& =\frac{x^3-15x^2+15x-1}{(x+1{)}^3}\\ R_4(x)& =\frac{x^4-28x^3+70x^2-28x+1}{(x+1{)}^4}\\ R_n(x)& =(x+1{)}^{-n}\sum _{m=0}^n(-1{)}^m(\genfrac{}{}{0ex}{}{2n}{2m})x^{n-m}\end{array}$

Partial fraction expansion

$R_n(x)=\sum _{m=0}^n\frac{(m!{)}^2}{(2m)!}(\genfrac{}{}{0ex}{}{n+m-1}{m})(\genfrac{}{}{0ex}{}{n}{m})\frac{(-4{)}^m}{(x+1{)}^m}$