Лекция 16

. Функции, ортогональные на точечном множестве.

Функции φ (x) и ψ (x) называют ортогональными на множестве

точек x = {x₀, x₁, x₂,…, x_n}, если

_n

Σ φ (x_і) × ψ (x_і) = 0.

^і=0

Пример: φ (x) = 3x² - 15x + 10; ψ (x) = 2x – 5

ортогональны на множестве x_і = і (і = 0,1,2,3,4,5).

Имеем: φ (0) = 10;

x
φ		-2	-8	-8	-2
ψ	-5	-3	-1
φ × ψ	-50			-8	-6		Σ = 0

Система функций {φ_k (x)} называется ортогональной на данном множестве X, если функции

системы попарно ортогональны между собой на этом множестве X.

Функция φ(x)≡0, естественно, ортогональна любой функции. Поэтому будем рассматривать только такие функции,

_n

Σ φ² (x_і) > 0

^і=0

т.е. хотя бы одно значение φ (x_і) ≠ 0

Система ортогональных функций {φ_k (x)} называется ортонормированной, если для всех k

выполнено равенство

║φ_k (x)║ = 1, где - норма функции φ(х)

Если {φ_k (x)} – система функций ортогональна на множестве Х, то система функций

{φ_k (x) ⁄ ║φ_k║_х} – ортогональная на Х.

Функции f_k (x) (k = 0, 1,…,m) называются линейно независимыми на множестве Х, если они

определены на этом множестве и из равенства

λ₀f₀ (x_і) + λ₁f₁ (x_і) +…+ λ_mf_m (x_і) = 0 (і = 1, n)

Следует, что все постоянные λ_k = 0 (k = 0,m)

В противном случае функции f_k (x) – линейно зависимые на Х.

Если множество Х не точечное, а континумум, т.е. х принадлежит отрезку a < х < b, то условие линейной независимости то же самое, только рассматриваются не точки x_і, а множество х € (a, b), т.е. рассматриваем условие

_m

Σ λ_kf_k (x) = 0 х € (a, b)

^k=0

Легко можно доказать лемму:

Функции φ_k (x) (k = 0, 1,…,m), ортогональные на множестве Х = {x₀, x₁, x₂,…,x_m} и имеющие ненулевые нормы, линейно независимы на этом множестве.

Рассмотрим систему полиномов

P₀ (x), P₁ (x),…, P_m (x), (1)

ортогональны на точечном множестве Х = {x₀, x₁,…,x_n}, т.е.

_n

Σ P_і (x_і) P_k (x_і) = 0 при j ≠ k (2)

^і=0 _{2 n 2}

и таких, что S_j = ║P_j║ = Σ P_j (x_і) > 0 – квадрат нормы

^{x і=0}

Пусть степень полинома P_j = j

Т.к. полином P_і (x) (j = 0, 1, 2,…,m) линейно независимые на Х поскольку они ортогональны, то любой полином Q_m (x) степени не выше m может быть представлен в виде линейной комбинации полиномов (2), т.е.

Q_m (x) = b₀ P₀ (x) + b₁P₁ (x) +…+ b_mP_m (x) (3),

где b_і (i = 0,m) – некоторые постоянные числа.

Выражение (3) называется разложением полинома Q_m (x) по системе (1)

Если полином P_j (x) ортогональны, то коэффициенты b_k равны

k = 0, 1,…,m (4)

Действительно, умножим (3) на полином P_k (x) (k ≤ m) и просуммируем результат по всем x_i

(i = 0,n)

_{n n n n 2 n}

Σ Q_m (x_i) × P_k (x_i) = b₀ Σ P₀ (x_i) × P_k (x_i) + b₁ Σ P₁ (x_i) × P_k (x_i) +...+ b_k Σ P_k (x_i) +...+ b_m Σ P_m (x_i) ×

^{і=0 і=0 і=0 і=0 і=0}

× P_k (x_i) (5)

В силу условия ортогональности из (5) следует (4), т.к. все Σ P_j (x) × P_k (x) = 0 для j≠k

Коэффициенты (4) называют коэффициентами Фурье полинома Q_m (x) относительно данной системы функций P_k (x) (k = 0,m), ортогональных на Х.

_n

Если система (1) ортогональна на Х, т.е. Σ Q_m (x_i) × P_k (x_i) (6)

^і=0

Теперь будем аппроксимировать заданную функцию y = f (х) полиномом Q_m (x) на множестве Х.

Q_m (x) = C₀P₀ (x) + C₁P₁ (x) +…+ C_mP_m (x),

где {P_k (x)} (k = 0,m), - ортогональные полиномы.

_{n 2}

Минимизируя S _m = Σ [C₀P₀ (x_i) + C₁P₁ (x_i) +…+ C_mP_m (x_i) - f (х_i)] ­по C₁, C₂,…, C_m, т.е. приравнивая

^і=0

и разрешая полученную систему уравнений, получим:

k= (7)

C_k – коэффициенты Фурье функции f (x) относительно ортогональной системы {P_k (x)} на Х

Беря вторую производную , можно убедиться, что и, следовательно

C_m (C_k) – минимальна!

Т.о.

Обобщенный полином фиксированного порядка m с коэффициентами Фурье данной функции f (х) на множестве Х = {x₀, x₁,…,x_n} – полином Фурье - обладает наименьшим квадратичным отклонением от этой функции на Х по сравнению со всеми полиномами того же порядка m.

Можно показать, что для полинома Фурье

_{2 m 2 2}

S_{m =} ║f (х)║ - Σ C_k × ║P_k║

^{k=0 x}

_{2 2 m n}

Если система {P_k (x)} ортонормирована, то ║P_k║ = 1 и S_{m =} ║f (х)║ - Σ C_{k ,} C_{k =} Σ f (х_i) × P_k (х_i)

^{х k=0 і=0}

III. Полиномы Чебышева, ортогональные на системе равноотстоящих точек.

Пусть дана система n + 1 равноотстоящих точек х = {x₀, x₁,…,x_n} с шагом n. С помощью

линейного преобразования переведем эти точки в t = 0, 1, 2,…, n.

Полиномы P_{0, n} (t), P_{1, n} (t),…, P_{m, n} (t) (m ≤ n) степеней 0, 1,…, m, ортогональные на множестве {0, 1, 2,…,n} и отличные от нуля на этом множестве, называются ортогональными полиномами Чебышева.

(P_{k, n} (t): k – степень полинома, n – число точек, уменьшенное на 1)

Полиномы Чебышева можно задать формулой

k s s s t^[s]____ (1)

P_{k, n} (t) = Σ (-1) C C ×

^{S = 0 k k k + s} n^[s]

где k = 0, 1,…, m;

s

C - число сочетаний из k по s

^k

t^[s] = t × (t - 1) × (t - 2) ×...× (t-s +1)

обобщенные степени t и n

n ^[s] = n × (n - 1) × (n - 2) ×...× (n – s+1)

Т.о. четыре первых ортогональных полинома равны:

P_{0, n} (t) = 1

; (n≥1) , k=1, s=1

; (n≥2) (2)

; (n≥3)

Возвращаясь к прежней переменной х, получим систему полиномов, ортогональных на множестве Х

, (k=0,1,…,m; m≤n)

Можно показать, что квадрат нормы полинома равен

^{(3) [n+k+1], [k+1] – обобщенные степени}

Разделив многочлены на их нормы, мы получим нормированную систему ортогональных полиномов Чебышева

(k=0,1,2,…,m; m≤n) (4)

Пример. Получить систему полиномов до третьей степени включительно, ортонормированных на системе точек ; ; ; ; ;

Решение. Полагая , переведем точки в целочисленные точки t=0,1,2,3,4,5. Теперь в формулах (1) принимаем n=5.

Имеем ;

(k=0, n=s)

k=1, n=5 s=0 1 2 переходим к Х

1-0,4t=1-0,4*

И т.д.

По формуле (3) вычисляем нормы по формуле (3) получаем квадрат нормы, затем вычисляем корень, т.е. норму!

Делим полиномы на их нормы и переходим от переменной t к переменной х, получим нормированную систему ортогональных полиномов Чебышева:

Если функция Y=f(x) задана на множестве узлов с шагом h, то наилучший аппроксимирующий ее на Х полином степени m будет иметь вид:

(5)

(k=0,1,2,…,m) (6)
-коэффициенты Фурье ф-ии f(x) относительно системы ортогональных полиномов Чебышева

Из (5) и (6) следует, что полином не изменится, если ортогональные полиномы Чебышева умножить на постоянные множители, отличные от нуля. Т.к.

Поэтому часто вместо полиномов пользуются полиномами подобраны так, чтобы для целочисленных значений аргумента t значения тоже были целочисленными. Имеются (в справочной литературе) таблицы полиномов , что значительно упрощает процедуры построения полиномов .

Дата добавления: 2014-01-04