Файл: Сахарников Н.А. Высшая математика учебник.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 09.04.2024

Просмотров: 223

Скачиваний: 2

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

^

а тт

Следовательно, экстремум может быть лишь в точке х0 — у 0 =

—. По смыслу

а2

потому что при

задачи имеем в этой точке относительный максимум и = - ,

4

 

 

, т. е.

и(х„+ Ах, у 0 + Ау) < и (ж0, уо).

 

§ 22. ВОПРОСЫ ТЕОРИИ ПРИБЛИЖЕНИЯ ФУНКЦИИ

Пусть в рассматриваемом явлении две величины х и у связаны функциональной зависимостью у = / (х), которая нам неизве­ стна. Однако известна таблица экспериментальных данных

 

 

 

 

 

 

 

 

 

 

 

(1)

содержащая приближенные значения величин х и у.

 

 

Для нахождения функции у = f

(х)

табличных

данных (1)

недостаточно, но можно поставить вопрос о нахождении функ­

ции,

приближенно

представляющей

/ (х),

причем

это

при­

ближение должно

быть наилучшим

в каком-то смысле,

т. е.

должно удовлетворять

некоторому

условию.

Возможна различ­

ная постановка вопроса о наилучшем приближении

функции.

Рассмотрим

некоторые

из них.

 

 

 

 

 

 

 

138.

Интерполяционная формула Лагранжа. Дана таблица (1).

Требуется найти многочлен Ln (х) степени п, принимающий в точ­

ках

х 0, х х,

. . ., хп соответственно

значения

у 0,

у 2, . . ., уп\

 

 

L n{xk) = yk при к = 0,

1, ... ,

п.

 

 

(2)

Геометрическая постановка вопроса

такова — требуется

найти

многочлен, график которого проходит через точки,

координаты

которых заданы таблицей (1).

в

виде

 

 

 

 

Запишем

искомый

многочлен

 

 

 

 

L n{x) = c0(x —x1) .. . [х — хп) + сх (х — х0)( х ~ х2) . . . (х — хп) +

 

 

+ . . . + С п(х — х0)(х — хг) . . . — хп_г),

 

 

3

 

 

 

 

( )

где

с0, сх,

. . ., сп — неопределенные

коэффициенты.

Полагая

в (3) последовательно

х = х0, х

=

х г, . . .,

х = хп,

получим

согласно (2) систему

равенств

 

 

 

 

 

 

 


из которой

легко

найти

все ck. Обозначим

 

со (х) — (х — х0)...

...(х хп).

Систему (4)

можно

записать

короче:

 

 

 

 

г/0=

с0ю' (ж0),

ÿx= c 1tû'(a:1),

. .

 

г/„ =

с„со' {хп).

 

 

Отсюда следует, что ck =

(Хк)

где

к =

0, 1, . .

 

п.

дает

При этих значениях коэффициентов ск

равенство

(3)

интерполяционную

формулу

Лагранжа

 

 

 

 

 

 

 

 

 

 

 

 

 

 

П

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h=0

 

 

 

 

 

 

 

 

 

(5)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

График многочлена (5) проходит через

 

точки

Mk (xk,

yk).

Однако между этими точками он может уклоняться

от графика

функции / (х).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

П р и м е р . * Определим плотность 26%-го раствора фосфорной кислоты

Н3РО4 при 20 °С,

пользуясь следующими данными:

 

1,3350

 

 

 

 

у (плотности)

. . . 1,0764

1,1134

1,2160

 

 

 

 

а: (%Н3Р 0 4)

. . .

 

14

20

 

35

 

50

 

 

 

 

По формуле Лагранжа (5) при п = 4 получим

 

 

 

 

 

 

 

 

 

2

 

24

 

 

64

 

 

 

1

0,3350 = 1,1528.

 

у = 1 — J 0,0764 + -^- 0,1134+

0,2160 — ~

 

Действительное значение плотности 26%-го Н3РО4 составляет 1,1529.

І39.

О

методе

наименьших

квадратов.

В

естествознании

пользуются

так

называемыми эмпирическими формулами, соста­

вленными с помощью экспериментальных данных таблицы (1).

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

таких

формул

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

Ниже

изложена

идея

этого способа.

 

 

 

 

 

 

 

 

 

 

 

 

 

прибли­

Пусть дана таблица (1). Предположим, что в первом

жении величина

у

линейно

зависит

от х,

т. е.

 

 

 

 

 

 

 

 

 

у — ах-^Ъ.

 

 

 

 

 

 

 

(6)

Назовем уклонением (или отклонением) разность между точ­ ным значением функции (6) в точке xk и соответствующим зна­ чением yk из таблицы (1): гк = axk + Ъ ук.

Сумма квадратов уклонений есть функция величин а и Ъ:

и (а, Ь) = е\f + ej-f . . . +е £ .

(7)

В методе наименьших квадратов на коэффициенты а и b функ­ ции (7) накладывается условие — они должны доставлять мини­ мум сумме квадратов уклонений и (а, Ъ). Требуется найти а и Ъ,

удовлетворяющие этому условию.

* Этот пример взят из работы Л. М. БатунераиМ. Е. Позина «Математи­ ческие методы в химической технике», изд. 5-е. М., Химиздат, 1968, с. 740.


Для нахождения точки экстремума функции

71

и (а, 6) = 2 (axk + b - y k)2 к=о

вычислим ее частные производные по а и по b и приравняем их нулю :

иа =

 

=

+

 

 

=

(9)

и'ъ = 2

2 (axé+ 6ук)= 2а 2 х* -f 2 (га

1)6 2

*/* = О,

где индекс суммирования А изменяется от 0 до п.

 

системы

Введем обозначения для коэффициентов полученной

уравнений 2 Х1 =

Iх , x h 2 х* =

Ы , 2 z/fe == Ы , 2 хкѴп =

tx>ÿl.

Эти постоянные могут быть найдены непосредственно по

 

исход­

ным данным (1).

 

 

 

 

системе

Необходимые условия экстремума (9) приводят нас к

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

коэф­

фициентами

 

 

 

 

 

 

 

[X,

х]а -г {х)Ъ=[х, у],

[х]а + (і-{-іі)Ъ=-[у},

 

(10)

которая называется нормальной

системой

метода

наименьших

квадратов.

а0,

b = Ъ0 — решение системы (10). Можно дока­

Пусть а =

зать, что найденные а0 и Ь0 доставляют величине и (а, Ь) именно минимум. Функция (6) при значениях параметров а0 и Ь0 дает

эмпирическую формулу у = а0 х + Ь0.

Уточним эмпирическую формулу. Эмпирическая формула у = а0х + Ь0 характеризуется значением суммы квадратов укло­

нений и (а0, Ь0). Эту эмпирическую

формулу

можно

улучшить

так, чтобы уменьшилась сумма квадратов уклонений.

Для этого

положим

 

 

 

 

у = а0х + Ь 0 + сц>(х),

 

 

(И)

где с — искомый коэффициент, а ф (х)

— любая

данная функция

(например, х 2 или sin х), удовлетворяющая

условию 2 ф2 (хк) ^

Ф 0. Составим сумму квадратов уклонений функции (11)

П

П

 

 

 

V (с) = 2 К х* + &о+ сф Ы —г/*]2= 2

[ссры + 8*]2-

h=o

h-o

 

 

 

Найдем значение с, при котором эта сумма минимальна. Для

этого составим необходимое

условие экстремума

П

 

= 2 2 и

ы + гк\ ф ы = о,

h

 


из которого следует, что с = с0 = —

При этом

 

Z.ф2(**)

значении с равенство (11) дает уточненную эмпирическую формулу

у =

а 0х +

&о +

с0 ф (ж).

сумма квадратов

уклонений

меньше

Если

с0 ф. О,

то новая

прежней:

ѵ (с0) •< и (а0, 60). Действительно,

ѵ" =

2 ^

фа (^) > 0 ,

поэтому

V (с)

в точке с0 имеет именно минимум.

Других

точек

экстремума функция

ѵ (с) не имеет, при с =

0 она имеет значе­

ние

V (0)

=

и (а0,

b0).

Отсюда следует

наше утверждение.

 

Процесс улучшения эмпирической формулы может быть про­

должен

по

рассмотренной

схеме путем

добавления

новых

сла­

гаемых.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

140.

 

О

равномерном

приближе­

 

 

 

 

 

 

нии

функции многочленом. Иногда

 

 

 

 

 

 

целесообразна

постановка вопроса о

 

 

 

 

 

 

наилучшем

приближении в

смысле

 

 

 

 

 

 

равномерного

приближения.

 

 

 

 

 

 

 

Пусть дана функция

/ (х). Гово­

 

 

 

 

 

 

рят,

что

функция

ф (х)

в

проме­

 

 

 

 

 

 

жутке

[а,

Ъ\

 

равномерно

прибли­

 

 

 

 

 

 

жает

функцию /

(X) с точностью до

 

 

 

 

 

 

е, если

для

всех

х из

[a,

b]

выпол­

 

 

 

 

 

 

нено

условие |/ (х)

ф (х) I

<Г е.

 

 

Тейлора

Тп (х)

Можно

доказать (см. п. 216), что многочлен

функции / (х) при достаточно большом

п

равномерно

прибли­

жает / (х) в промежутке

х 0| <; г, если

/ (х)

имеет производ­

ные всех порядков, которые ограничены одним и тем же числом.

Ниже сформулирована более сильная теорема

(налагающая

меньшие требования на функцию) о равномерном

приближении

функции многочленом.

в промежутке

Теорема Вейерштрасса. Если f (х) непрерывна

[а, 6], то для любого е

> 0 существует соответствующий много­

член Рп (х) такой,

что для всех х

из [а,

б]

выполняется условие

 

 

I f ( x ) - P n(x )\< z .

 

(12)

Геометрически

условие (12)

означает,

что график функции

Рп (х) не выходит

из

пределов

полосы

(рис. 96)

f{x) — e < P n (х) <f ( x ) + в.

Общая теория наилучшего приближения функций с помощью

многочленов

создана П.

Л. Чебышевым *.

Способ

построения

многочленов,

равномерно приближа­

ющих данную непрерывную / (х), был

предложен в частности

С. Н. Бернштейном. Не нарушая общности, можно рассматривать