Русская Википедия:Ортогональные многочлены
В математике последовательностью ортогональных многочленов называют бесконечную последовательность действительных многочленов
- <math>p_0(x),\ p_1(x),\ p_2(x),\ \ldots</math>,
где каждый многочлен <math>p_n(x)</math> имеет степень <math>n</math>, а также любые два различных многочлена этой последовательности ортогональны друг другу в смысле некоторого скалярного произведения, заданного в пространстве <math>L^2</math>.
Понятие ортогональных многочленов было введено в конце XIX в. в работах П. Л. Чебышёва по непрерывным дробям и позднее развито А. А. Марковым и Т. И. Стилтьесом и нашло различные применения во многих областях математики и физики.
Определение
Ортогональность с весом
Пусть <math>(a,b)</math> — промежуток на вещественной оси (конечный или бесконечный). Этот промежуток называется интервалом ортогональности. Пусть
- <math>w : ~(a,b) \to \mathbb{R}</math>
заданная непрерывная, строго положительная внутри промежутка <math>(a,b)</math> функция. Такая функция называется весовой или просто весом. Функция <math>w(x)</math> связана с пространством функций <math>L_{2}</math>, для которых сходится интеграл
- <math>\int_{a}^{b} \left[f(x)\right]^2 w(x) \; dx < \infty </math>.
В полученном пространстве можно ввести скалярное произведение по формуле
- <math>\langle f, g \rangle = \int_{a}^{b} f(x) g(x) w(x) \; dx</math> для вещественных функций,
- <math>\langle f, g \rangle = \int_{a}^{b} f(x) \overline{g(x)} w(x) \; dx</math> для комплекснозначных функций.
Если скалярное произведение двух функций равно нулю <math>\langle f, g \rangle = 0</math>, то такие функции называются ортогональными с весом <math>w(x)</math>. Как правило, среди ортогональных полиномов рассматриваются только вещественные функции.
Классическая формулировка
Систему многочленов
- <math>p_0(x),p_1(x),\cdots,p_n(x),\cdots</math>
называют ортогональной, если
- <math>p_n(x)</math> — многочлен степени <math>n</math>,
- <math>\langle p_m,p_n \rangle = \delta_{mn} h_{n}</math>, где <math>\delta_{mn}</math> — символ Кронекера, <math>h_{n}</math> — нормировочный множитель.
Ортогональный базис называется ортонормированным, если все его элементы имеют единичную норму <math>||p_n||=h_n=1</math>. Некоторые классические многочлены, представленные ниже, могут быть нормированы по какому-либо другому правилу. Для таких многочленов значения <math>h_n</math> отличаются от единицы и указаны в таблице внизу.
Общие свойства последовательностей ортогональных многочленов
Рекуррентные соотношения
Любые ортогональные полиномы удовлетворяют следующей рекуррентной формуле, связывающей три последовательных многочлена из системы:
- <math>{p_{n+1}(x)\ =\ (A_nx+B_n)\ p_n(x)\ -\ C_n\ p_{n-1}(x)},</math>
где
- <math>A_n=\frac{k_{n+1}}{k_n},\quad B_n=A_n \left(r_{n+1}-r_n \right), \quad C_n= \frac {A_n h_n} {A_{n-1} h_{n-1}}, </math>
- <math>r_n=\frac{k'_n}{k_n}, \quad h_n= \langle p_n(x),p_n(x) \rangle</math>,
- <math>k_n</math> и <math>k'_n</math> — коэффициенты при членах <math>x^n</math> и <math>x^{n-1}</math> в полиноме <math>p_n(x).</math>
Эта формула остаётся справедливой и для <math>n=0</math>, если положить <math>p_{-1}(x)=0</math>.
Формула Кристоффеля-Дарбу
<math>\sum_{k=0}^{n}\frac{p_k(x)p_k(y)}{h_k}=\frac{k_n}{k_{n+1}h_n}\frac{p_{n+1}(x)p_n(y)-p_{n+1}(y)p_n(x)}{x-y}</math>,
или при <math>y \to x</math>
<math>\sum_{k=0}^{n}h_k^{-1}\left[p_k(x)\right]^2=\frac{k_n}{k_{n+1}h_n}\left[p'_{n+1}(x)p_n(x)-p_{n+1}(x)p'_n(x)\right]</math>
Корни многочленов
Все корни многочлена <math>p_n(x)</math> являются простыми, вещественными и все расположены внутри интервала ортогональности <math>\left[a;b\right]</math>. Шаблон:Hider
Между двумя последовательными корнями многочлена <math>p_n(x)</math> расположен в точности один корень многочлена <math>p_{n+1}(x)</math> и, по крайней мере, один корень многочлена <math>p_m(x)</math>, при <math>m>n</math>.
Минимальность нормы
Каждый многочлен <math>p_n(x)</math> в ортогональной последовательности имеет минимальную норму среди всех многочленов <math>P_n(x)</math> такой же степени и с таким же первым коэффициентом.
Полнота системы
Система ортогональных многочленов <math>p_i(x)</math> является полной. Это значит, что любой многочлен <math>S(x)</math> степени n может быть представлен в виде ряда
- <math>S(x)=\sum_{i=0}^n {\alpha}_i\ p_i(x)</math>,
где <math>\alpha</math> коэффициенты разложения.
{{Hider|
title = Доказательство| hidden = 1 | title-style = text-align: left; | content-style = text-align: left; | content =
Доказывается с помощью математической индукции. Выберем <math>\ {\alpha}_n</math> так, чтобы <math>\ S(x)-{\alpha}_n P_n(x)</math> был многочленом степени меньше <math>n</math>. Далее по индукции. }}
Дифференциальные уравнения, приводящие к ортогональным многочленам
Очень важный класс ортогональных многочленов возникает при решении дифференциального уравнения следующего вида:
<math>{Q(x)}\,f + {L(x)}\,f' + {\lambda}f = 0,</math>
где <math>Q(x)</math> и <math>L(x)</math> заданные многочлены второго и первого порядка, соответственно, а <math>f(x)</math> и <math>\lambda</math> неизвестные функция и коэффициент. Это уравнение называется задачей Штурма — Лиувилля и может быть переписано в его более стандартной форме
<math>(R(x)y')' + W(x)\,\lambda\,y = 0,</math>
где <math>R(x) = e^{\int\frac{L(x)}{Q(x)}\,dx},\, W(x) =\frac{R(x)}{Q(x)}.</math> Решение этого уравнения приводит к множеству собственных чисел <math>{\lambda}_0, {\lambda}_1, {\lambda}_2, \dots</math> и множеству собственных функций <math>P_0, P_1, P_2, \dots</math>, обладающих следующими свойствами:
- <math>P_n(x)</math> — полином степени n, зависящий от <math>{\lambda}_n</math>
- последовательность <math>P_0, P_1, P_2, \dots</math> ортогональна с весовой функцией <math>W(x)</math>
- Промежуток ортогональности зависит от корней многочлена Q, причём корень L находится внутри промежутка ортогональности
- Числа <math>\lambda_n</math> и полиномы <math>P_n(x)</math> могут быть получены из формул
- <math>{\lambda}_n = - n \left( \frac{n-1}{2} Q + L' \right)</math>
- <math>P_n(x) = \frac{1}{e_n\,W(x)} \ \frac{d^n}{dx^n}\left(W(x)[Q(x)]^n\right)</math> формула Родрига.
Дифференциальное уравнение имеет нетривиальные решения только при выполнения одного из следующих условий. Во всех этих случаях при изменении масштаба или/и сдвига области определения и выбора способа нормировки многочлены решения сводятся к ограниченному набору классов, которые называются классическими ортогональными полиномами
- 1. Якобиподобные многочлены
- Q — многочлен второго порядка, L — первого. Корни Q различны и действительны, корень L лежит строго между корнями Q. Первые коэффициенты Q и L имеют один знак. При помощи линейного преобразования уравнение сводится к <math>Q(x)=1-x^2</math> с интервалом ортогональности <math>[-1,1]</math>. Решениями являются многочлены Якоби <math>P_n^{(\alpha, \beta)}(x)</math> или их частные случаи многочлены Гегенбауэра <math>C_n^{(\alpha)}(x)</math>, Лежандра <math>P_n(x)</math> или Чебышёва обоих типов <math>T_n(x)</math>, <math>U_n(x)</math>.
- 2. Лагерроподобные многочлены
- Q и L — многочлены первого порядка. Корни Q и L различны. Первые коэффициенты Q и L имеют один знак, если корень L меньше корня Q и наоборот. Сводится к <math>Q(x)=x</math> и интервалу ортогональности <math>[0,\infty)</math>. Решениями являются обобщённые многочлены Лагерра <math>L_n^{(\alpha)}(x)</math> или их частному случаю многочленам Лагерра <math>L_n(x)</math>.
- 3. Эрмитоподобные многочлены
- Q — ненулевая константа, L — многочлен первого порядка. Первые коэффициенты Q и L имеют противоположный знак. Сводится к <math>Q(x)=1</math> и интервалу ортогональности <math>(-\infty,\infty)</math>. Решениями являются многочлены Эрмита <math>H_n(x)</math>.
Производные ортогональных полиномов
Обозначим <math>P_n^{m}(x)</math> как m-ую производную полинома <math>P_n(x)</math>. Производная <math>P_n^{m}(x)</math> является полиномом степени <math>n-m</math> и обладает следующими свойствами:
- ортогональность
- Для заданного m последовательность полиномов <math>P_m^{m}, P_{m+1}^{m}, P_{m+2}^{m}, \dots</math> ортогональна с весовой функцией <math>W(x)[Q(x)]^m</math>
- <math>P_n^{m} = \frac{1}{e_nW(x)[Q(x)]^m} \ \frac{d^{n-m}}{dx^{n-m}}\left(W(x)[Q(x)]^m\right)</math>
- дифференциальное уравнение
- <math>{Q(x)}\,y + (m\,Q'(x)+L(x))\,y' + [{\lambda}_n-{\lambda}_m]\ y = 0</math>, где <math>y(x)=P_n^{m}(x)</math>
- дифференциальное уравнение второго вида
- <math>(R(x)[Q(x)]^m\ y')' + [\lambda_n-\lambda_m]W(x)[Q(x)]^{m}\ y = 0</math>, где <math>y(x)=P_n^{m}(x)</math>
- рекуррентные соотношения (для удобства у коэффициентов a, b и c опущены индексы n и m)
- <math>P_n^{m}(x) = aP_{n+1}^{m+1}(x) + bP_n^{m+1}(x) + cP_{n-1}^{m+1}(x),</math>
- <math>P_n^{m}(x) = (ax+b)P_n^{m+1}(x) + cP_{n-1}^{m+1}(x),</math>
- <math>Q(x)P_n^{m+1}(x) = (ax+b)P_n^{m}(x) + cP_{n-1}^{m}(x).</math>
Классические ортогональные многочлены
Классические ортогональные полиномы, которые происходят из дифференциального уравнения, описанного выше, имеют много важных приложений в таких областях как: математическая физика, численные методы, и многие другие. Ниже приводятся их определения и основные свойства.
Многочлены Якоби
Многочлены Якоби обозначаются <math>P_n^{(\alpha, \beta)}(x)</math>, где параметры <math>\alpha</math> и <math>\beta</math> вещественные числа больше −1. Если <math>\alpha</math> и <math>\beta</math> не равны, полиномы перестают быть симметричными относительно точки <math>x=0</math>.
- Весовая функция <math>W(x)=(1-x)^\alpha(1+x)^\beta</math> на промежутке ортогональности <math>[-1,1]</math>
- Дифференциальные уравнения
- <math>(1-x^2)\,y + (\beta-\alpha-[\alpha+\beta+2]\,x)\,y' + {\lambda}\,y = 0</math>
- Собственные числа
- <math>\lambda_n = n(n+1+\alpha+\beta)</math>
- Рекуррентная формула
- <math>P_{n+1}(x) = (A_n\,x+B_n)\,P_n(x) - C_n\,P_{n-1}(x),</math>
- где
- <math>A_n=\frac{(2n+1+\alpha+\beta)(2n+2+\alpha+\beta)}{2(n+1)(n+1+\alpha+\beta)},</math>
- <math>B_n=\frac{({\alpha}^2-{\beta}^2)(2n+1+\alpha+\beta)}{2(n+1)(2n+\alpha+\beta)(n+1+\alpha+\beta)},</math>
- <math>C_n=\frac{(n+\alpha)(n+\beta)(2n+2+\alpha+\beta)}{(n+1)(n+1+\alpha+\beta)(2n+\alpha+\beta)}</math>
- Нормировка
- <math>P_n^{(\alpha, \beta)}(1)=\frac{\Gamma(n+1+\alpha)}{n!\,\Gamma(1+\alpha)}, \qquad h_n=\frac{2^{\alpha+\beta+1}\,\Gamma(n\!+\!\alpha\!+\!1)\,\Gamma(n\!+\!\beta\!+\!1)}
{n!(2n\!+\!\alpha\!+\!\beta\!+\!1)\Gamma(n\!+\!\alpha\!+\!\beta\!+\!1)}, \qquad k_n=\frac{\Gamma(2n+1+\alpha+\beta)}{n!\,2^n\,\Gamma(n+1+\alpha+\beta)}, \qquad e_n=(-2)^n\,n!</math>
Многочлены Гегенбауэра
Многочлены Гегенбауэра обозначаются <math>C_n^{(\alpha)}(x)</math>, где параметр <math>\alpha</math> вещественное число больше −1/2. Он выводится из многочленов Якоби для равных параметров <math>\alpha</math> и <math>\beta</math>
<math>C_n^{(\alpha)}(x) = \frac{\Gamma(2\alpha\!+\!n)\,\Gamma(\alpha\!+\!1/2)} {\Gamma(2\alpha)\,\Gamma(\alpha\!+\!n\!+\!1/2)}\! \ P_n^{(\alpha-1/2, \alpha-1/2)}.</math>
Остальные Якобиподобные многочлены являются частным случаем полиномов Гегенбауэра с выбранным параметром <math>\alpha</math> и соответствующей нормализацией.
- Весовая функция <math>W(x)=(1-x^2)^{\alpha-1/2}</math> на промежутке ортогональности <math>[-1,1]</math>
- Дифференциальные уравнения
- <math>(1-x^2)\,y - (2\alpha+1)\,x\,\,y' + {\lambda}\,y = 0</math>
- Собственные числа
- <math>\lambda_n = n(n+2\alpha)</math>
- Рекуррентная формула
- <math>(n+1)\,C_{n+1}^{(\alpha)}(x) = 2(n+\alpha)x\,C_n^{(\alpha)}(x) - (n+2\alpha-1)\,C_{n-1}^{(\alpha)}(x)</math>
- Нормировка
- <math>C_n^{(\alpha)}(1)=\frac{\Gamma(n+2\alpha)}{n!\,\Gamma(2\alpha)}</math> если <math>\alpha\ne0,\qquad</math> <math> h_n=\frac{\pi\,2^{1-2\alpha}\Gamma(n+2\alpha)}{n!(n+\alpha)(\Gamma(\alpha))^2}, \qquad k_n=\frac{\Gamma(2n+2\alpha)\Gamma(\frac{1}{2}+\alpha)}{n!\,2^n\,\Gamma(2\alpha)\Gamma(n+\frac{1}{2}+\alpha)}, \qquad e_n = \frac{(-2)^n\,n!\,\Gamma(2\alpha)\,\Gamma(n+\frac{1}{2}+\alpha)}
{\Gamma(n+2\alpha)\Gamma(\alpha+\frac{1}{2})}</math>
- Прочие свойства
- <math>C_n^{(\alpha+1)}(x) = \frac{1}{2\alpha}\! \ \frac{d}{dx}C_{n+1}^{(\alpha)}(x)</math>
Многочлены Лежандра
Многочлены Лежандра обозначаются <math>P_n(x)</math> и являются частным случаем многочленов Гегенбауэра с параметром <math>\alpha=1/2</math>
<math>P_n(x) = C_n^{(1/2)}(x).</math>
- Весовая функция <math>W(x)=1</math> на промежутке ортогональности <math>[-1,1]</math>
- Дифференциальные уравнения
- <math>(1-x^2)\,y - 2x\,y' + {\lambda}\,y = 0, \qquad ([1-x^2]\,y')' + \lambda\,y = 0</math>
- Собственные числа
- <math>\lambda_n = n(n+1)</math>
- Рекуррентная формула
- <math>(n+1)\,P_{n+1}(x) = (2n+1)x\,P_n(x) - n\,P_{n-1}(x)</math>
- Нормировка
- <math>P_n(1)=1, \qquad h_n=\frac{2}{2n+1}, \qquad k_n=\frac{(2n)!}{2^n\,(n!)^2}, \qquad e_n=(-2)^n\,n!</math>
- Первые несколько многочленов
- <math>P_0(x)=1;</math>
- <math>P_1(x)=x;</math>
- <math>P_2(x)=(3x^2-1) / 2;</math>
- <math>P_3(x)=(5x^3-3x) / 2;</math>
- <math>P_4(x)=(35x^4-30x^2+3) / 8;</math>
Многочлены Чебышёва
Многочлен Чебышёва <math>T_n(x)</math> часто используется для аппроксимации функций как многочлен степени <math>n</math>, который меньше всего отклоняется от нуля на интервале <math>[-1,1]</math>
<math>T_n(x) = \cos(n\,arccos(x)).</math>
Является частным случаем нормированного многочлена Гегенбауэра для параметра <math>\alpha \to 0</math>
<math>T_n(x) = \lim_{\alpha \to 0}n\,\Gamma(\alpha)\,C_n^{(\alpha)}.</math>
- Весовая функция <math>W(x)=(1-x^2)^{-1/2}</math> на промежутке ортогональности <math>[-1,1]</math>
- Дифференциальное уравнение
- <math>(1-x^2)\,y - x\,y' + {\lambda}\,y = 0</math>
- Собственные числа
- <math>\lambda_n = n^2</math>
- Рекуррентная формула
- <math>T_{n+1}(x) = 2x\,T_n(x) - T_{n-1}(x)</math>
- Нормировка
- <math>T_n(1)=1, \qquad h_n=\left\{
\begin{matrix} \pi &:~n=0 \\ \pi/2 &:~n\ne 0 \end{matrix}\right. ,\qquad k_n = 2^{n-1}, \qquad e_n=(-2)^n\,\frac{\Gamma(n+1/2)}{\sqrt{\pi}} </math>
Многочлен Чебышёва второго рода <math>U_n(x)</math> характеризуются как многочлен, интеграл от абсолютной величины которого на интервале <math>[-1, +1]</math> меньше всего отклоняется от нуля
<math>U_n = \frac{1}{n+1}\,T_{n+1}'</math>
- Весовая функция <math>W(x)=(1-x^2)^{1/2}</math> на промежутке ортогональности <math>[-1,1]</math>
- Дифференциальное уравнение
- <math>(1-x^2)\,y - 3x\,y' + {\lambda}\,y = 0</math>
- Нормировка
- <math>U_n(1)=n+1, \qquad h_n=\pi / 2, \qquad k_n = 2^n, \qquad e_n=2(-2)^n\,\frac{\Gamma(n+3/2)}{(n+1)\,\sqrt{\pi}}</math>
Многочлены Лагерра
Ассоциированные или обобщённые многочлены Лагерра обозначаются <math>L_n^{(\alpha)}(x)</math>, где параметр <math>\alpha</math> вещественное число больше -1. Для <math>\alpha = 0</math> обобщённые многочлены сводятся к обычным многочленам Лагерра
<math>L_n(x) = L_n^{(0)}(x).</math>
- Весовая функция <math>W(x)=x^{\alpha}e^{-x}</math> на промежутке ортогональности <math>[0,\infty)</math>
- Дифференциальные уравнения
- <math>x\,y + (\alpha + 1-x)\,y' + {\lambda}\,y = 0\, \qquad (x^{\alpha+1}\,e^{-x}\, y')' + {\lambda}\,x^{\alpha}\,e^{-x}\,y = 0</math>
- Собственные числа
- <math>\lambda_n = n</math>
- Рекуррентная формула
- <math>(n+1)\,L_{n+1}^{(\alpha)}(x) = (2n+1+\alpha-x)\,L_n^{(\alpha)}(x) - (n+\alpha)\,L_{n-1}^{(\alpha)}(x)</math>
- Нормировка
- <math>k_n=\frac{(-1)^n}{n!}, \qquad h_n=\frac{\Gamma(n+\alpha+1)}{n!}, \qquad e_n=n!</math>
- Прочие свойства
- <math>L_n^{(\alpha+1)}(x) = - \frac{d}{dx}L_{n+1}^{(\alpha)}(x)</math>
Многочлены Эрмита
- Весовая функция <math>W(x)=e^{-x^2}</math> на промежутке ортогональности <math>[-\infty,\infty]</math>
- Дифференциальные уравнения
- <math>y - 2xy' + {\lambda}\,y = 0\, \qquad (e^{-x^2}\,y')' + e^{-x^2}\,\lambda\,y = 0</math>
- Собственные числа
- <math>\lambda_n = 2n</math>
- Рекуррентная формула
- <math>H_{n+1}(x) = 2x\,H_n(x) - 2n\,H_{n-1}(x)</math>
- Нормировка
- <math>k_n = 2^n, \qquad h_n=2^n\,n!\,\sqrt{\pi}, \qquad e_n=(-1)^n</math>
- Первые несколько многочленов
- <math>H_0(x) = 1</math>
- <math>H_1(x) = 2x</math>
- <math>H_2(x) = 4x^2-2</math>
- <math>H_3(x) = 8x^3-12x</math>
- <math>H_4(x) = 16x^4-48x^2+12</math>
Построение ортогональных многочленов
Процесс ортогонализации Грама-Шмидта
Система ортогональных многочленов <math>f_1, f_2, \ldots, f_k</math> может быть построена путём применения процесса Грама-Шмидта к системе многочленов <math>g_k(x)=x^k</math> следующим образом. Определим проектор как
- <math>\mathrm{proj}_{f}\,(g) = {\langle f, g\rangle\over\langle f, f\rangle}f = { \int_{x_1}^{x_2} f(x) g(x) W(x) \; dx \over \int_{x_1}^{x_2} (f(x))^2 W(x) \; dx} f(x)</math>,
тогда ортогональные полиномы последовательно вычисляются по схеме
- <math>
\begin{align} f_1 & = g_1, \\ f_2 & = g_2-\mathrm{proj}_{f_1}\,(g_2), \\ f_3 & = g_3-\mathrm{proj}_{f_1}\,(g_3)-\mathrm{proj}_{f_2}\,(g_3), \\ & {}\ \ \vdots \\ f_k & = g_k-\sum_{j=1}^{k-1}\mathrm{proj}_{f_j}\,(g_k). \end{align} </math>
Данный алгоритм относится к численно неустойчивым алгоритмам. При вычислении коэффициентов разложения ошибки округления и погрешности численного интегрирования накапливаются с увеличением номера полинома.
По моментам весовой функции
Весовая функция <math>w(x)</math>, заданная на промежутке <math>\left[a;b\right]</math>, однозначно определяет систему ортогональных многочленов <math>\{p_n(x)\}_{n=0}^{\infty}</math> с точностью до постоянного множителя. Обозначим через числа
- <math>\mu_n=\int_a^b{w(x)x^ndx}</math>
моменты весовой функции, тогда многочлен <math>p_n(x)</math> может быть представлен в виде:
- <math> p_n(x) = \det\left[
\begin{matrix} \mu_0 & \mu_1 & \mu_2 & \cdots & \mu_n \\ \mu_1 & \mu_2 & \mu_3 & \cdots & \mu_{n+1} \\ \mu_2 & \mu_3 & \mu_4 & \cdots & \mu_{n+2} \\ \vdots & \vdots & \vdots & & \vdots \\ \mu_{n-1} & \mu_n & \mu_{n+1} & \cdots & \mu_{2n-1} \\ 1 & x & x^2 & \cdots & x^n \end{matrix} \right] </math>.
Сложность вычисления ортогональных полиномов определяется сложностью вычисления определителя матрицы. Существующие алгоритмические реализации вычисления требуют минимум <math>O(n^3)</math> операций.
По рекуррентным формулам
Если выбрать нормировку многочлена <math>p_n(x)</math> таким образом, что коэффициент <math>k_n</math> при главном члене равен единице, рекуррентное соотношение может быть переписано в следующем виде:
- <math>{p_{n+1}(x)\ =\ (x-\alpha_n)\ p_n(x)\ -\ \gamma_n\ p_{n-1}(x)},</math>
где
- <math>\alpha_n = \frac{\langle x p_n, p_n\rangle}{\langle p_n, p_n\rangle}, \qquad \gamma_n = \frac{\langle x p_n, p_{n-1}\rangle}{\langle p_{n-1}, p_{n-1}\rangle}</math>.
Применение ортогональных многочленов
Ортогональные полиномы применяются для построения точных квадратурных формул
<math> \int\limits_{\Omega} f(x) w(x) dx \approx \sum\limits_{i=1}^{n} w_i f(x_i), </math>
где <math>x_i</math> и <math>w_i</math> являются узлами и весами квадратурной формулы. Квадратурная формула является точной для всех полиномов <math>f(x)</math> до степени <math>2n-1</math> включительно. При этом узлы <math>x_i</math> есть корни n-го полинома из последовательности полиномов <math>p_0(x),p_1(x), ...</math>, ортогональных с весовой функцией <math>w(x)</math>. Веса <math>w_i</math> вычисляются из формулы Кристоффеля-Дарбу.
Так же многочлены Чебышёва первого <math>T_n(x)</math> и второго <math>U_n(x)</math> типа часто используется для аппроксимации функций.
Примечания
Ссылки
Для дальнейшего чтения
Шаблон:Rq Шаблон:Ортогональные многочлены