Приложение к разделу I
ФОРМУЛА ТЕЙЛОРА. ВЫВОД ФОРМУЛЫ ТЕЙЛОРА

Если функция y=f(x) имеет в точке Х₀ производную, то ее приращение можно представить в виде

y=x+0(x), ∆x

где , y₀= и , т. е. =y₀+A(x-x₀)+0(x-x₀).

. Иначе говоря, существует линейная функция

P₁(x)=y₀+A(x-x₀) (13.1)

такая, что

f(x)=P₁(x)+o(x-x₀), х x₀,

причем

p₁ (x₀) =y₀=f (x₀), Р’₁ (x₀)=А=f'(х₀).

Поставим более общую задачу. Пусть функция f имеет в точ­ке Х₀ n производных. Требуется выяснить, существует ли мно­гочлен Р_n(х) степени не выше п, такой, что

f(x)=P_n(x)+o((x-x₀)ⁿ), xx₀, (13.2)

f (x₀)= P_n(x₀), f'(х₀)= Р’_n (x₀),....,f⁽ⁿ⁾(x₀)= Р⁽ⁿ⁾_n (x₀) . (13.3)

Будем искать этот многочлен, по аналогии с формулой (13.1), в виде

P_n(x)=A₀₊A₁(x-x₀)+A₂(x-x₀)²+...+A_n(x-x₀)ⁿ.

Замечая, что Р_n (x₀)=A₀, из nеpвoгo условия (13.3), т. е. f(x₀)=Р_n(х₀),имем A₀=f(x₀). Далее,

Р’_n (x)= A₁+2A₂(х-х₀)+... +nA_n(х-х₀)^n-1,

отсюда P^’_n(x₀)=A₁, и так как P^’_n(x₀)=f^’(x₀),то A₁= f^’(x₀). Затем найдем вторую производную многочлена Р_n(x):

Р^''_n (x) = 2 • 1 • A₂,+... +n (n - 1)-A_n(х – х₀)^n-2. Отсюда и из условия f" (x₀)= Р^''_n (x₀) получим вообще k=0, 1, 2,..., n.

В силу самого построения, для многочлена Р_n(х)=

=f(x₀)+f^’(x₀)(x+x₀)+...+ (x-x₀₎+...+

выполнены все соотношения (13.3). Проверим, удовлетворяет ли он условию (13.2). Пусть

r_n(x) =f(x)-P_n(x)

Из условия (13.3) следует, что -

r_n(x₀)= r^’_n(x₀)=0 (13.4)

Поэтому, применяя п раз правило Лопиталя для раскрытия нео­пределенностипри х->хо, а именно сначала п—1 раз

теорему 2 из § 12, а затем теорему 1 того же параграфа, получаем:


т. е. действительно. r_n(x)=0 ((x-x₀)ⁿ),

Итак, доказана следующая очень важная теорема. Теорема 1. Пусть функция f (х) определенная на интервале (а, b), имеет в точке Х₀ (а, b) производные до порядка п вклю­чительно. Тогда при

f(x)=f(x₀)+ (13.5)

или f(x)=

Эта теорема остается справедливой, вмсте с ее доказательством, и для функции /, определенной на отрезке[а, b], при х₀ [а, b]
если для x₀=a и x₀=b под производные понимать соответствующие односторонние производные.

Формула (13.5), называется формулой Тейлора п-гопорядка с остаточным членом в форме Псчно. Многочлен

f(x)=f(x₀)+ (13.6) называется многочленом Тейлора функция

r_n(x)=f(x)-P_n(x) (13.7)

—остаточным членом п-го порядка форму Тейлора. Как пока­зано, остаточный член r_n(x) является бесконечно малой, при х, более высокого порядка, чем все члены многочлена Тейлора (13.6).

Укажем другой вид записи формулы (13.5). Полагая

Применение аппроксимирующих дифференциальных уравнений для решения обыкновенных дифференциальных уравнений.

Для численного решения обыкновенных дифференциальных уравнений (ОДУ) можно использовать вспомогательные, аппроксимирующие дифференциальные уравнения (АДУ) со свободными параметрами.

Рассмотрим пример. Пусть требуется решить ОДУ:

(1.1.1) y'=f(x,y), у(хо)=уд, [Xo,B].

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

(1.1.2) Yj'-KjyrO, j-=0,l,...,n.

(1.1.3) y(Xj+hj)=y(Xj) exp(Kj(Xj+hj)) —решение (1.1.2). Величины Kj определяются на каждом j-м шаге из (1.1.2);

(1.1.4) Kf=y'(Xj+hj)/y(Xj+hj).

Формула (1.1.3) является экстраполяционной.

Решение (1.1.3) можно уточнить с помощью расчётной формулы интерполяционного типа:

(1.1.5) у(х, +h,)=y(x,)-exp(K„(x,+h,)).

Здесь

(1.1.6) Ko,=f(Xf+hj,y(Xj+hj))/y(Xj+hj). Пример/ Для задачи Коши

(1.1.7) у'=у-2х/у, у(хо)=1, хо=0.

Вычислим величины y(Xg+h) и yf-Xg+h) при h=0,4 по расчётным

формулам (1.1.3) и (1.1.5):

(1.1.За) у(0+0,4)=у(0)-ехр(К,(0+0,4))=1,492; у(0)= 1; К, = 1;

6 =((1,492 -1,342) /1,342) -100% ==11%. (1.1.5а) у(0 + 0,4) = у(0) • ехр(К„ (0 + 0,4)) = 1,285; у(0) = 1; К„ = 0,638;

8 = ((1,285 -1,342)/1,342) • 100% =. 3,8%.

При решении (1.1.1) с использованием МВС вычисления проводятся с h=ho, ho/I, 2ho для контроля ^выбора шага.

Формулы (1.1.3) и (1.1.5) характеризуются тем, что в качестве аппроксимирующих функций используются экспоненциальные функции.

Применение ОМНК при регуляризации таблично заданных функций своими значениями и значениями своих производных

Пусть аппроксимируемая функция задана своими значениями и значениями своих производных/^ (Xkm), k=0, 1, .... г, т=1, 2, ..., mk.

При необходимости сглаживания аппроксимируемой функции Fn(C„ x) коэффициенты С, определяются по А.Н.Тихонову из условия минимума выражения

(1) Ф(F_n)+λ_iΨ(F_n)=min

В случае применения ОМНК

(2) Ф(F_n)=
(3)
В общем случае функционал Ψ (F_n) выбирается из условия, что, если значение этого функционала невелико, то функция Fn(C_i, x) обладает определённой гладкостью.

При λ_k=0 метод регуляризации переходит в ОМНК.

При r=0 и λ_k=0 (=о метод регуляризации А.Н.Тихонова переходит в классический метод наименьших квадратов Лежандра-Гаусса.

УДК517.91+518(075.8)+519.6

6Л2 Гольцов НА. Лекции о применении теории вероятностей для решения задач прикладного численного анализа.— М.: МГУ Л, 2001.

Автор: Гольцов Николай Алексеевич Издается в авторской редакции

| Московский государственный университет леса, 2002

Лицензия ЛР №020718 от 2.02.1998 г. Лицензия ПД №00326 от 14.02.2000 г.

Подписано к печати 17.5.2002 Формат 60х88/16 Бумага 80 г/м² "Снегурочка" Ризография Объем____п.л.______________Тираж 100 экз._______Заказ №_____

Издательство Московского государственного университета леса 141005, Московская обл., Мытищи, ул. 1-я Институтская, 1, МГУЛ. Телефон: (095)588-57-62 E-mail: [email protected]

Литература:

1. 
   Gauss C.F., Theoria interpolationist methodo nova tractada.—Werke, 1805. 13.3, p. 265-330.
   
2. 
   Legendre A.M., Nouvelles methodes pour la determination des orlites comets.—Paris, 1806
   
3. 
   Гаус К.Ф., Способ наименьших квадратов // Избранные геодезические сочинения / Под. ред. Г.В.Багратуни—М.: Геодезиздат, 1957.—144 стр.
   
4. 
   Гаус К.Ф., Способ наименьших квадратов.—М.: ГИФМЛ, 1959.—160 стр.
   
5. 
   Гмурман В.Е., Теория вероятностей и математическая статистика.—М.: Высшая школа, 1998.—478 стр.
   
6. 
   Гмурман В.Е., Руководство к решению задач по теории вероятности и математической статистике.—М.: Высшая школа, 1998.—400 стр.
   
7. 
   Венцель Е.С., Теория вероятностей.—М.: Наука 1968.—576 стр.
   
8. 
   Гольцов Н.А., Конструирование алгоритмов прикладного численного анализа для многопроцессорных вычислительных комплексов.—М.: Изд.МГУЛ, 2001.—55 стр.
   
9. 
   Соболь И.М., Численные методы Монте-Карло.—М.: Наука, 1972.—312 стр.
   
10. 
    Сурков А.В. Использование MS Excel для решения экономических задач.—М.: Изд.МГУЛ, 2000.—123 стр.
    

11. Колмогоров А.Н. К обоснованию метода наименьших квадратов // Успехи математических наук.-1946. Вып. 1 (11). Т. 1.-С. 57-70
2 Дополнительная литература

1. 
   Митропольский А.К., Элементы математической статистики.—Л.: Издат. ВЗЛТИ, 1969.—205 стр.
   
2. 
   Математическая энциклопедия. Авторы и предметный указатель даны в пятом томе.
   

Лекции подготовлены к печати студентами В.Б.Вакулинским, М.В.Максимовым, В.К.Бабаяном.