Boyarshinov_ChM_T1

Принимая во внимание выражение (4.30), окончательно получаем

F ~c v Av,v F ~c .

Сучетом положительной определенности матрицы А последнее выражение приводит

кнеравенству F ~c v F ~c v Rn 1 , но это и означает минимальность функционала

аргумента .

В силу допущения о минимальности функционала (4.29) имеем F ~c u F ~c , то есть g g 0 . Это, в свою очередь означает, что =0 доставляет минимум функции

g( ), откуда следует

Поскольку последнее равенство справедливо u Rn 1 , получаем Ac~ f 0, что и требовалось доказать.

В компонентной записи система линейных алгебраических уравнений выглядит следующим образом:

Алгоритм определения элемента наилучшего приближения:

b

- вычисление коэффициентов матрицы akp jk ,jp H jk x jp x dx, k,p 0,n;

a

119

b

- вычисление значений правых частей fp f,jp H f x jp x dx, p 0,n ;

a

- решение системы уравнений (4.31);

В этом случае погрешность приближения может быть оценена следующим образом:

f2H 2 ~ck f, k H ~ck~cp k , p H f2H 2 ~ck2 ~ck~cp kp f2H ~ck2.

k 0 k,p 0 k 0 k,p 0 k 0

n

коэффициенты Фурье.

На рис. 4.6 приведены графики приближения | x | c помощью степенных рядов.

30 Фурье Жан Батист Жозеф [21.3.1768 - 16.5.1830] - французский математик. Окончил военную школу в Осере, там же работал преподавателем. В 1796 - 1798 годах преподавал в Политехнической школе в Париже. В 1798 году принимал участие в Египетской экспедиции Наполеона Бонапарта. В 1802 -1815 годах был префектом департамента Изер. С 1817 года избран членом Парижской академии наук.В 1829 году стал иностранным почетным членом Петербургской академии наук.

120

В практических исследования часто возникает ситуация, когда необходимо аппроксимировать табличные значения fk f xk , k 0,n с помощью приближения j(x),

содержащего определяемые коэффициенты ck , k 0,m в количестве, меньшем чем число узловых точек, m < n. По этой причине, в отличие от рассмотренных ранее способов аппроксимации функции полиномами Ньютона, Лагранжа, сплайнами, не используется условие f xk j xk , k 0,n - равенство значений функции f(x) и ее приближения j(x) для заданного числа значений аргумента. Так, в рассматриваемом методе наименьших квадратов, “близость” аппроксимирующего многочлена к самой функции оценивается с помощью какой-либо нормы, то есть “в среднем” для всего отрезка, на котором строится аппроксимация. Для получения алгоритма построения приближения воспользуемся полученными соотношениями (4.26) - (4.31).

Пусть известен набор значений fk f xk , k 0,n функции для ряда значений ее аргумента. Для рассматриваемого случая положим H = Rn 1. В линейном пространстве размерности (n+1) скалярное произведение и норма определяются известным образом:

n

u u,u ukuk .

k 0

121

Пусть отыскиваемое приближение j x;c0,c1, ,cm зависит от известного числа m

параметров c0,c1, ,cm .

Степень отклонения функции f(x) от ее приближения j(x) определяется соотношением

f 2 f2 2 f, 2

fk2 2 fk xk ;c0,c1, ,cm 2 xk;c0,c1, , cm .

Для определения наименьшего отклонения воспользуемся необходимыми условиями минимума функции нескольких переменных:

Иными словами, речь идет о решении системы алгебраических уравнений, в общем случае - нелинейных:

В частном случае, когда приближение j(x) представимо в виде

m

j x;c0,c1, ,cm jk x ,

k 0

оценку (4.32) отклонения функции от ее приближения можно записать в форме:

122

f 2 fk2 2 fk xk ;c0,c1, ,cm 2 xk ;c0,c1, ,cm

Условие минимальности отклонения приближения от функции записывается аналогично представленному выше:

В итоге получена система линейных алгебраических уравнений относительно коэффициентов разложения cp, p 0,m:

Пример приближения функции | x | на отрезке [-1, 1] с использованием метода

наименьших квадратов приведен на рис. 4.7.

123

Контрольные вопросы и задания

Сформулируйте задачу аппроксимации. Каковы условия разрешимости этой задачи?

Укажите требования к аппроксимирующим функциям. В каком случае интерполяция называется линейной?

Что представляют собой разделенные разности? Поясните их геометрический смысл.

Укажите порядок построение интерполяционного многочлена Ньютона. Что представляет собой схема Горнера?

Укажите порядок построение интерполяционного многочлена Лагранжа.

Покажите, что полиномы Ньютона и Лагранжа, построенные на одном множестве табличных значений функции, тождественны.

Приведите оценку погрешности аппроксимации функции, заданной таблично,

полиномом Ньютона (Лагранжа).

Укажите условия сходимости процесса интерполяции полиномами. Приведите примеры.

Обоснуйте преимущества метода интерполяции для решения нелинейного уравнения.

Опишите способ решение нелинейного уравнения с помощью обратной интерполяции.

Укажите порядок построение интерполяционного многочлена Эрмита.

Опишите идею сплайн-аппроксимации функции и порядок построения кубического сплайна.

124

Что понимается под сходимостью процесса интерполяции кубическими сплайнами?

Сформулируйте и докажите лемму об оценке сходимости по "сеточной" норме.

Сформулируйте и докажите теорему о сходимости процесса интерполяции функции кубическими сплайнами.

Опишите порядок построения наилучшего приближения функции с использованием теории гильбертовых пространств.

Вчем заключается метод наименьших квадратов для аппроксимации функции,

заданной таблично?

125

степени n от переменной . Это, в свою очередь, означает, что существует n корней

126

Определим первый собственный вектор, соответствующий 1 6,

Здесь обозначено: X11, X12 - компоненты вектора X1 .

Очевидно, что последняя система содержит линейно зависимые уравнения (как это и следовало ожидать при det A 1E 0). Используя любое из уравнений этой системы,

получим X11 2 X12 . 5

Примем для однозначности определения собственных векторов условие нормирования

Xk 1, k 1,n , то есть

Теперь найдем второй собственный вектор, соответствующий 2 1,

Отсюда получаем связь компонент второго собственного вектора: X12 X22 .

Из условия нормирования следует:

Устойчивость собственных значений и векторов

Рассмотрим две матрицы, A и AT . Собственные числа обеих матриц, в чем нетрудно

убедиться, одинаковы. Собственные векторы матриц обозначим соответственно

Xk, Yk, k 1,n :

AXi iXi , ATY j jY j.

Для скалярных произведений Y j,AXi и ATYj,Xi вычислим разность

Y j,AXi ATYj,Xi Yj, iXi jY j,Xi i j Yj,Xi .

127

Поскольку

Для оценки устойчивости собственных значений в формуле (5.4) положим i = j :

A Xi ,Yi A Xi ,Yi i Xi ,Yi i Xi ,Yi .

Поскольку

A Xi ,Yi Xi ,ATYi Xi , iYi i Xi ,Yi ,

получаем

A Xi ,Yi i Xi ,Yi ,

i Xi ,Yi A Xi ,Yi A Xi Yi AXi Yi.

называется коэффициентом перекоса; i - угол между векторами Xi, Yi .

Далее будем предполагать, что матрица A симметрична31, то есть A AT , и все ее собственные числа различны. В этом случае собственные векторы матрицы А образуют в

31 В этом случае, согласно [7], все собственные числа матрицы А вещественны.

128