ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 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 ) + в.
Общая теория наилучшего приближения функций с помощью
многочленов |
создана П. |
Л. Чебышевым *. |
|
Способ |
построения |
многочленов, |
равномерно приближа |
ющих данную непрерывную / (х), был |
предложен в частности |
С. Н. Бернштейном. Не нарушая общности, можно рассматривать