- •Лекция 1
- •Лекция 2
- •Лекция 3
- •Лекция 4
- •Лекция 5
- •Лекция 6
- •Лекция 7
- •Лекция 8
- •Лекция 9
- •Лекция 10
- •Лекция 11
- •Лекция 12
- •Лекция 13
- •Лекция 14
- •Лекция 15
- •Лекция 16
- •Лекция 17
- •Лекция 18
- •Лекция 19
- •Лекция 20
- •Лекция 21
- •Лекция 22
- •Лекция 23
- •Лекция 24
- •Лекция 25
- •Лекция 26
- •Лекция 27
- •Лекция 28
- •Лекция 29
- •Лекция 30
- •Лекция 31
- •Лекция 32
- •Лекция 33
- •Лекция 34
- •Лекция 35
- •Лекция 36
- •Лекция 37
- •Лекция 38
- •Лекция 39
- •Лекция 40
- •Лекция 41
- •Лекция 42
- •Лекция 43
- •Лекция 44
- •Лекция 45
- •Лекция 46
- •Лекция 47
|
|
|
|
ЛЕКЦИЯ 23 |
|
|
|
|
Точка ( x0 |
, y0 , λ |
0 ) , удовлетворяющая системе |
f x′ |
( x, y) + λg′x ( x, y) =0 |
||||
|
|
|
|
|||||
|
|
|
|
|
f y′( x, y) + λg′y ( x, y) =0 |
|||
называется стационарной точкой функции Лагранжа L( x, y, λ) . |
|
|||||||
Теорема (достаточные условия условного экстремума) |
|
|
||||||
Пусть функции |
f ( x, y), g( x, y) имеют непрерывные частные производные |
|||||||
второго порядка в окрестности стационарной точки ( x0 , y0 , λ0 ) |
функции Ла- |
|||||||
гранжа L( x, y, λ) = f ( x, y) + λg( x, y) . Тогда если |
в |
этой |
точке |
число |
||||
L′′ |
L′′ |
|
g′ |
|
|
z = f ( x, y) |
в точке |
|
Q =( g′y , g′x ) xx |
xy |
|
y |
< 0, то целевая функция |
||||
′′ |
′′ |
− g′x |
|
|
|
|
||
Lyx |
Lyy |
|
|
|
|
|||
( x0 , y0 ) имеет условный максимум, а при Q > 0 - условный минимум. Доказательство проводить не будем. Идея та же, что и в достаточных
условиях экстремума: разложить f ( x, y), g( x, y) по формуле Тейлора в точке ( x0 , y0 ) и сложить результаты (умножив предварительно g( x, y) на λ0 ), получим разложение функции Лагранжа по формуле Тейлора, которое нач-
нется с d 2 L( x0 , y0 , λ0 ) и знак приращения ∆z определится d 2 L. |
|
||
Пример 1. |
Найти экстремум |
функции z = xy при |
условии |
2 x + 3 y − 5 =0. |
|
|
|
Решение. |
Составим |
функцию |
Лагранжа: |
L( x, y, λ) = xy + λ(2 x + 3 y − 5) .
∂L =∂∂xL =
Составим систему ∂y
∂∂λL =
y + 2λ =0 |
|
x + 3λ =0 |
. Решаем ее и получаем |
|
|
2 x + 3 y − 5 =0 |
|
λ = − 5 2 , x = 5 4 , y = 5 6 . Далее L′xx′ = 0, L′yy′ = 0, L′xy′ =1, g′y = 3, g′x = 2.
|
|
|
|
0 |
|
1 |
3 |
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
||
(3, − 2) |
|
|
|
|
|
|
=− 6 |
− 6 =− 12 |
< 0 . |
Итак, |
|
в |
точке |
||||||||||
|
=(−2, 3) |
− 2 |
|
|
|||||||||||||||||||
|
|
|
|
1 0 − 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
5 |
|
5 |
|
|
|
|
|
|
|
|
|
|
5 |
|
5 |
|
|
25 |
|
|
||
M0 |
|
|
, |
|
|
|
|
функция имеет условный максимум |
zmax = z |
|
, |
|
|
= |
|
. |
|
||||||
|
6 |
|
4 |
6 |
24 |
|
|||||||||||||||||
|
4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Пример 2. При каких размерах прямоугольная банка объемом 32 см3 ,
открытая сверху, имеет наименьшую поверхность?
S = xy + 2 xz + 2 yz → min
Решение. xyz = 32 . Здесь 3 переменных, исключим
одну:
z = |
32 |
|
, S = xy + 2 |
|
32 |
( x + y) → min, |
∂S |
= y − |
64 |
|
=0, |
∂S |
= x − |
|
64 |
|
=0 y = |
64 |
, |
|||||||||||||||
xy |
|
xy |
∂x |
x2 |
∂y |
|
y2 |
x2 |
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
x − |
|
|
64 |
|
|
=0, x − |
|
x |
4 |
=0 x =4, |
y =4, |
z =2. |
∂2 S |
= |
128 |
, |
∂2 S |
= |
128 |
, |
∂2 S |
=1. |
|
|||||||||||
|
642 |
|
|
64 |
∂x2 |
|
x3 |
∂y2 |
|
y3 |
∂x∂y |
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
x4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
1 |
=4 − 3 =1 |
> 0 в точке (4, 4, 2) min S = 48 см2 . |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
1 |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Упражнение. При каких значениях диаметра основания d и высоты h цилиндрическая банка, объем которой равен 54π , имеет наименьшую поверхность?
Пример 3. Найти наибольшее и наименьшее значения функции
z = x2 + y2 в круге ( x −
2 )2 + ( y −
2 )2 ≤ 9 .
|
|
Решение. Находим стационарные точки: |
∂z |
= 2 x = 0, |
∂z |
= 2 y = 0 |
||||||||||||||||||||||||||||
|
|
∂x |
|
∂y |
|
|||||||||||||||||||||||||||||
|
∂2 z |
=2, ∂2 z2 = |
2, |
∂2 z |
=0 |
|
2 0 |
|
=4 > 0 O(0,0) − |
точка |
минимума. Мини- |
|||||||||||||||||||||||
|
|
|
||||||||||||||||||||||||||||||||
|
2 |
∂x∂y |
||||||||||||||||||||||||||||||||
|
∂x |
∂y |
|
|
|
|
|
|
|
|
0 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
мум глобальный, значит, |
zнаим = zmin = z(0,0) =0 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
Теперь |
|
|
проведем |
|
|
|
исследование |
|
на |
|
|
|
границе |
|
|
области: |
||||||||||||||||||
|
L( x, y, λ) = x 2 + y 2 + λ(( x − 2 )2 + ( y − 2 )2 − 9) . |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
∂L = 2 x + 2λ( x − 2 ) =0, |
∂L = 2 y + 2λ( y − 2 ) =0, |
∂L |
=( x − 2 )2 + ( y − 2 )2 =9 |
||||||||||||||||||||||||||||||
|
∂x |
|
|
|
|
|
|
|
|
|
|
|
∂y |
|
|
|
|
|
|
∂λ |
|
|
|
|
|
|
|
|
|
|||||
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Решая |
|
|
|
|
|
|
|
|
|
|
эту |
|
|
|
|
|
систему, |
|
|
|
|
|
|
|
|
получим |
||||||||
λ1 = − 5 , x1 = y1 = |
5 2 |
|
, λ2 |
= − 1 |
, x2 = y2 = − |
2 . |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
3 |
|
|
|
|
|
2 |
|
|
|
|
3 |
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
Смотрим |
|
|
|
|
|
далее |
|
|
|
|
|
по |
|
достаточным |
|
|
|
|
|
условиям: |
||||||||||||||
|
′′ |
|
|
|
|
′′ |
|
|
|
|
|
|
′′ |
′ |
− |
|
|
′ |
|
|
|
|
|
2 ) . Тогда в точ- |
||||||||||
|
Lxx =2 |
+ 2λ, Lyy =2 + 2λ, Lxy =0, gx =2( x |
2 ), g y =2( y − |
|||||||||||||||||||||||||||||||
|
|
|
5 2 |
|
5 2 |
|
|
5 |
|
имеем (3 2, − 3 |
|
− |
4 |
3 |
|
0 |
|
|
|
3 2 |
|
|
||||||||||||
|
|
|
, |
, − |
|
|
|
0 |
|
4 |
|
|
|
|
=− 48 < 0 . |
|||||||||||||||||||
ке M1 |
2 |
|
|
|
2 |
3 |
|
2 ) |
|
|
− |
3 |
|
− 3 |
2 |
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
5 2 5 2 |
|
|
|
|
|
|
|
|
|
2 |
|
2 |
|
|
1 |
|
|
|
|
||||||||
|
zmzx = z |
наиб |
|
|
|
|
|
, |
|
|
|
|
|
=25. |
|
|
|
|
|
|
, − |
|
|
, − |
|
|
имеем: |
|||||||
|
= z |
|
2 |
|
|
2 |
|
В точке M 2 − |
2 |
|
2 |
3 |
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
(− 3 2, 3 |
|
1 |
3 |
0 |
|
− 3 |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
2 ) |
|
|
|
1 |
|
|
|
|
|
= 12 > 0 - точка минимума. |
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
0 |
3 |
3 |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
Есть достаточные условия экстремума и в случае большего числа переменных и уравнений связи, но они, естественно, более сложны в написа-
нии. С экономической точки зрения множители Лагранжа λi (их столько, сколько условий связи) есть оценка единицы ресурса (в задачах линейного программирования c1 x1 +K+ cn xx → max, Ax ≤ b ).
Метод наименьших квадратов
Важное практическое значение имеет следующая задача. Пусть требуется установить функциональную зависимость между двумя переменными
величинами x, y по результатам |
n экспериментальных измерений, |
приве- |
||||||||
денных в таблице : |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
x1 |
x2 |
|
….. |
xi |
….. |
xn |
|
|
|
y |
y1 |
y2 |
|
….. |
yi |
….. |
yn |
|
|
Другими словами, требуется выразить зависимость между x, y |
анали- |
|||||||||
тически, то есть дать формулу, связывающую между собой соответствующие значения переменных. Формулы, служащие для аналитического представления опытных данных, принято называть эмпирическими формулами.
Следует заметить, что подбор эмпирической формулы не ставит задачу разгадать истинный смысл зависимости – эта задача математически неразрешима.
Ставится задача подобрать формулу, в каком-то смысле наилучшим образом отображающую полученные результаты. Применяются для этой цели различные методы. Можно построить многочлен, принимающий в заданных точках заданные значения. Достоинство этого способа – полученная формула точно воспроизводит заданные значения. Такого рода формулы называются интерполяционными. Но здесь по двум точкам строится прямая, по трем – парабола, по n точкам – многочлен степени (n − 1) , то есть с ростом n растет степень многочлена.
Можно разбить опытные данные на тройки точек и по каждой тройке строить параболу. В результате получится функция, “склеенная” из парабол. Достоинство, как и в предыдущем случае, недостаток – громоздкая запись, недифференцируемость в местах “склейки”. Это так называемая сплайн – интерполяция.
Третий способ заключается в следующем. Исходя из некоторых теоретических или практических соображений (можно нарисовать эти точки на плоскости) подбирается наиболее простая формула, которая дает наилучшее совпадение с опытными данными. Обычно в экономике достаточно линей-
ной, квадратичной, логарифмической, показательной, степенной и обратнопропорциональной зависимостей.
Слова “наилучшее совпадение” понимается здесь в том смысле, что из данного множества формул вида y = f ( x, a, b) наилучшим образом отображающим данные значения считается та, для которой сумма квадратов отклонений наблюдаемых значений от вычисленных по формуле, является наи-
n
меньшей, то есть ∑( yi − f ( xi , a, b))2 → min .
i=1
Итак, пусть дана таблица с n экспериментальными данными и из теоретических соображений или из характера расположения экспериментальных
точек ( xi , yi ), i =1, n на плоскости XOY известно, что эмпирическую формулу следует искать в виде y =ax + b , где параметры a, b подлежат определению. Следуя способу наименьших квадратов составим функцию двух пе-
ременных a и b : S ( a , b ) = ∑n |
( ax i + b − y i ) 2 |
, для которой |
i = 1 |
|
|
нужно найти минимум по a и b . Из необходимого условия экстремума имеем:
∂S |
|
|
|
∂S |
=0, т.е. |
∂S |
n |
|
|
|
|
2 |
′ |
+ b − y1 ) |
2 |
|
|
|
||||
|
=0, |
|
|
|
|
|
|
|
|
|
|
|
|
+ |
|
|
||||||
∂a |
∂b |
∂a |
= ∑(axi + b − yi ) |
|
a =[(ax1 |
|
|
|
||||||||||||||
|
|
|
|
|
|
i=1 |
|
|
|
|
|
|
|
|
|
|
||||||
(ax2 + b − y2 ) |
2 |
+K+ (axn + b |
− yn ) |
2 |
′ |
|
|
|
|
|
|
|
||||||||||
|
|
+]a = |
|
|
|
|
|
|
|
|||||||||||||
2(ax1 + b − y1 )x1 + 2(ax2 + b − y2 )x2 +K+ 2(axn + b − yn )xn = |
|
|
|
|
||||||||||||||||||
|
|
|
n |
|
|
n |
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
= 2 a∑xi2 + b∑xi − ∑xi yi , |
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
i=1 |
|
|
i=1 |
|
i=1 |
|
|
|
|
|
|
|
|
|
|
|
|
||
∂S |
|
|
|
n |
|
+ b − yi ) |
2 |
′ |
|
|
n |
|
|
|
|
n |
|
n |
|
|
||
∂b |
= |
∑(axi |
|
b =2∑(axi + b − yi )==2 a∑xi + bn − ∑ yi |
|
|||||||||||||||||
|
i=1 |
|
|
|
|
|
|
|
|
i=1 |
|
|
|
i=1 |
|
i=1 |
|
|
||||
|
|
Приравнивая к нулю, получим линейную систему двух уравнений с |
||||||||||||||||||||
двумя неизвестными a и b , которая называется нормальной системой. |
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
n |
|
n |
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a∑xi2 + b∑xi =∑xi yi |
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
i=n1 |
|
i=1 |
|
i=1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
a∑xi + bn=∑ yi |
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
i=1 |
|
i=1 |
|
|
|
|
|
|||
|
|
Легко показать, |
что на решении системы M0 (a0 , b0 ) матрица Гессе для |
|||||||||||||||||||
функции |
|
|
S(a, b) |
|
|
|
положительно |
определена, |
|
то |
есть |
|||||||||||
′′ |
|
|
|
|
|
′′ ′′ |
|
|
′′ |
2 |
> 0 и |
точка |
|
M0 (a0 , b0 ) есть |
точка |
минимума |
||||||
Saa |
> 0, ∆ = Saa Sbb |
− ( Sab ) |
|
|
||||||||||||||||||
функции S(a, b) .
Пример. Рост производительности труда на предприятии за 5 лет отражен в следующей таблице:
года |
|
x |
1 |
2 |
3 |
4 |
5 |
|
|
|
|
|
|
|
|
Среднее кол-во |
|
|
|
|
|
|
|
продукции, |
вы- |
y |
235 |
250 |
270 |
|
|
пускаемой |
ра- |
292 |
300 |
||||
бочим за смену |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Полагая эту зависимость линейной, установить ее аналитическую формулу. Решение. Уравнение искомой зависимости y =ax + b . Коэффициенты нормальной системы удобно искать с помощью следующей расчетной табли-
цы:
i |
xi |
yi |
xi2 |
xi yi |
1 |
1 |
835 |
1 |
235 |
2 |
2 |
250 |
4 |
500 |
3 |
3 |
270 |
9 |
810 |
|
|
|
|
1168 |
4 |
4 |
292 |
16 |
|
5 |
5 |
300 |
25 |
1500 |
∑15 1347 55 4213
Подставляем |
полученные |
в |
таблице |
данные |
и |
получаем: |
|
55a + 15b =4213 |
|
|
|
|
|
|
|
15a + 5b =1347 . Решая эту систему, получаем a |
=17,2, b = 217,8 . Искомая |
||||||
эмпирическая формула имеет вид: |
y =17,2 + 217,8 . |
|
|
|
|||
Предположим теперь, |
что |
f ( x) следует |
|
искать |
в виде |
параболы |
|
|
|
|
n |
|
|
|
|
f ( x) =ax2 + bx + c . Тогда S(a, b, c) =∑(axi2 + bxi |
+ c − yi )2 . Составим сис- |
||||||
i=1
тему для отыскания стационарной точки:
∂S = |
∑n [( ax i2 + bx i + c − y i ) 2 ]′a = ∑n |
2( ax i2 + bx i + c − y i )x i2 = |
|
|
|||||||||||||||||||||||||||||||||
∂a |
i =1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i =1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∑ 2( ax i4 + bx i3 + cx i2 − y i x i2 ) = |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
i =1 |
|
|
+ b ∑ x i3 |
+ c ∑ x i2 − ∑ x i2 |
y i |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
= 2 a ∑ x i4 |
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
|
n |
|
|
|
n |
|
|
|
|
n |
|
|
|
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
i =1 |
|
|
|
i =1 |
|
|
i =1 |
|
|
|
|
|
|
|
i =1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
∂S = |
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 ∑ ( ax i2 + bx i + c − y i )x i = |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
∂b |
i =1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
|
|
|
n |
|
|
|
|
|
n |
|
|
|
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
= 2 a ∑ x i3 + b ∑ x i2 + c ∑ x i − ∑ x i y i |
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
|
i =1 |
|
|
|
i =1 |
|
|
|
i =1 |
|
|
|
|
|
|
i =1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
∂S = 2 ∑ ( ax i2 + bx i |
+ c − y i ) = 2 a ∑ x i2 + b ∑ x i + cn − ∑ y i |
|
|
||||||||||||||||||||||||||||||||||
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
|
|
|
n |
|
|
|
|
n |
|
|
||
∂c |
i =1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i =1 |
|
|
|
i =1 |
|
|
i =1 |
|
|
||||||||
Получаем нормальную систему: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
n |
|
|
|
n |
|
n |
|
|
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
a∑xi4 + b∑xi3 + c∑xi2 = |
∑xi2 yi |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
i=1 |
|
|
i=1 |
|
i=1 |
|
|
|
|
i=1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
n |
|
|
|
n |
|
n |
|
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
a∑xi3 + b∑xi2 + c∑xi = |
∑xi yi |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
i=1 |
|
|
|
i=1 |
|
i=1 |
|
n |
|
i=1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
n |
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
a∑xi2 + b∑xi + cn = ∑ yi |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
i=1 |
|
|
|
i=1 |
|
|
|
|
i=1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
По гиперболе |
f ( x) =a + b |
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
n |
|
|
|
b |
|
2 |
|
∂S |
|
|
|
|
n |
|
|
b |
|
|
|
∂S |
|
n |
b |
|
|
1 |
|
|||||||||
S(a, b) =∑ |
|
|
|
|
|
|
= |
2∑ |
|
+ |
|
|
|
= |
|
|
|
|
. |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
a + |
|
xi |
− yi , |
|
∂a |
a |
xi |
− yi , |
∂b |
2∑ a + |
xi |
|
− yi |
|
|||||||||||||||||||||||
|
i=1 |
|
|
|
|
|
|
|
|
|
i=1 |
|
|
|
|
|
|
i=1 |
|
xi |
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
|
1 |
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
na + b∑ |
|
|
|
|
=∑ yi |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
xi |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
1 |
|
|
i=1 |
|
|
1 |
i=1 |
|
|
|
y |
. |
|
|
|
|
|
|
|
|
|
||||||
Нормальная система |
n |
|
|
|
|
|
|
n |
|
|
n |
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
a |
∑ |
|
|
|
+ b∑ |
|
|
|
=∑ |
|
|
i |
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
xi |
|
|
2 |
|
|
xi |
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
i=1 |
|
|
|
|
i=1 |
xi |
|
i=1 |
|
|
|
|
|
|
|
|
|
||||||||||||
Если |
выравнивание |
|
ведется |
по |
показательной |
функции |
y =abx , |
то |
|||||||||||||||||||||||||||||
можно предварительно прологарифмировать. Тогда имеем ln y =ln a + x ln b . Обозначая через Y =ln y, A =ln a, B =ln b , получим линейную функцию Y = A + Bx . Определив A и B , получим a =e A , b =e B .
Не существует общего правила для выбора подходящего вида эмпирической формулы. Для наиболее часто встречающихся зависимостей с двумя параметрами эмпирическую формулу можно выбирать с помощью следующей таблицы:
|
|
|
№п/п |
|
|
|
|
|
xs |
|
|
|
|
|
|
|
|
ys |
|
|
|
|
Вид эмпирической |
|
фор- |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
мулы |
|
|
|
|
|
|
|
|||
|
|
|
|
1 |
|
( x1 + xn ) |
|
|
|
( y1 + yn ) |
|
|
|
y ax + b |
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
x1 xn |
|
|
|
|
|
y1 yn |
|
|
|
|
|
|
|
y axb |
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
3 |
|
( x1 + x2 ) |
|
|
|
y1 yn |
|
|
|
|
|
y =ab x , y = aeβx , β =ln b |
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4 |
|
2x1 xn |
( x1 + xn ) |
|
( y1 + yn ) |
|
|
y =a |
+ b |
x |
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
5 |
|
( x1 + xn ) |
|
|
|
2 y1 yn |
( y1 + yn ) |
|
|
y = 1 |
(ax |
+ b) |
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
6 |
|
2x1 xn |
( x1 + xn ) |
|
2 y1 yn |
( y1 + yn ) |
|
|
y = x |
(ax |
+ b) |
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
7 |
|
|
x1 xn |
|
|
|
|
|
( y1 + yn ) |
|
|
y =a lg x + b |
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Пользоваться этой таблицей можно лишь в том случае, если в исход- |
|||||||||||||||||||||||||||||||||||
ных данных xi+1 − xi |
> 0, а |
yi+1 − yi обладает постоянным знаком. |
|
|
|
|||||||||||||||||||||||||||||||||
|
|
|
Для проверки пригодности выбранной эмпирической формулы, ис- |
|||||||||||||||||||||||||||||||||||
пользуя исходные данные находят значения |
|
, |
|
. Затем сравнивают ys , |
||||||||||||||||||||||||||||||||||
xs |
ys |
|||||||||||||||||||||||||||||||||||||
соответствующее |
|
|
|
|
в исходных данных, со значением |
|
. Если |
|
|
|
||||||||||||||||||||||||||||
xs |
ys |
xs |
|
не на- |
||||||||||||||||||||||||||||||||||
ходится среди исходных данных xi , то соответствующее значение |
yˆ s |
|
можно |
|||||||||||||||||||||||||||||||||||
определить с помощью линейной интерполяции |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
ˆ |
= |
|
+ |
yi+1 |
− yi |
|
|
|
|
|
|
− |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
yi |
( xs |
xi ) , |
где |
|
xi , xi+1 |
- промежуточные значения, |
|
между |
||||||||||||||||||||||||||||||
|
|
|
|
|
||||||||||||||||||||||||||||||||||
ys |
|
|
xi+1 |
− xi |
|
|
|
|
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
которыми находится |
|
xs ( xi |
< xs < xi+1 ) . |
Затем вычисляется |
|
|
ˆ |
− |
ys |
|
. Если |
|||||||||||||||||||||||||||
|
|
|
ys |
|
|
|||||||||||||||||||||||||||||||||
эта величина большая, то соответствующая эмпирическая формула непригодна.
Пример. Себестоимость y (в долларах) одного экземпляра книги в зависимости от тиража x (тысячи экземпляров) характеризуется данными, собранными издательством в течение ряда лет. Определить вид эмпирической формулы, отвечающей этой зависимости и найти ее параметры.
x |
1 |
2 |
3 |
5 |
10 |
20 |
30 |
50 |
100 |
200 |
|
|
|
|
|
|
|
|
|
|
|
y |
10,15 |
5,52 |
4,08 |
2,85 |
2,11 |
1,62 |
1,41 |
1,30 |
1,21 |
1,15 |
|
|
|
|
|
|
|
|
|
|
|
Составим таблицу:
№п/п |
|
|
|
|
xs |
|
|
|
|
ys |
|
|
|
yˆ s |
|
yˆ s − |
|
|
Вид эмпирической |
||||
|
|
|
|
|
|
|
|
|
|
|
ys |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
формулы |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
1 |
|
x1 + xn |
= 100, |
|
y1 + yn |
|
= 5,6 |
1,21 |
4,44 |
|
|
y =ax + b |
не подхо- |
||||||||||
|
|
|
2 |
|
2 |
|
|
|
|
|
|
|
|
дит |
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
2 |
|
x1 xn |
≈ 14 |
|
y1 yn |
≈ 3,42 |
1,91 |
1,51 |
|
|
y = axb |
не подхо- |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
дит |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
3 |
|
x1 + xn |
= 100, |
|
y1 yn |
≈ 3,42 |
1,21 |
1,21 |
|
|
y ab x |
не подходит |
|||||||||||
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
4 |
|
2x1 xn |
|
≈ 1,94 |
|
y1 + yn |
|
= 5,6 |
5,52 |
0,13 |
|
|
y = a + |
b |
|
подходит |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
x1 + xn |
|
2 |
|
|
|
|
|
|
|
|
x |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
лучше всего |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
5 |
|
x1 + xn |
|
|
|
2 y1 yn |
|
≈ 2,0 |
1,21 |
0,85 |
|
|
1 |
|
|
|
|||||||
|
|
|
|
|
= 100, |
|
|
|
|
|
|
|
|
y = |
|
не подхо- |
|||||||
|
|
|
2 |
|
|
y1 + yn |
|
|
|
|
|
|
ax + b |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
дит |
|
|||
6 |
|
2x1 xn |
|
≈ 1,94 |
|
2 y1 yn |
|
≈ 2,0 |
5,52 |
3,46 |
|
|
y = |
x |
не подхо- |
||||||||
|
|
x1 + xn |
|
|
y1 + yn |
|
|
|
|
|
|
ax + b |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
дит |
|
|||
7 |
|
x |
1 |
x |
n |
≈ 14 |
|
y1 + yn |
|
= 5,6 |
1,91 |
3,74 |
|
|
y =a lg x + b не подхо- |
||||||||
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
дит |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Для определения параметров a и b формулы |
y = a + |
b |
решим нормальную |
||||||
x |
|||||||||
|
|
|
|
|
|
|
|
||
систему, коэффициенты которой найдем по следующей таблице: |
|||||||||
i |
xi |
yi |
1 xi |
1 xi2 |
yi xi |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
1 |
1 |
10,15 |
1,000 |
1,000 |
10,150 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
2 |
5,52 |
0,500 |
0,250 |
2,760 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
3 |
4,08 |
0,333 |
0,111 |
1,360 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4 |
5 |
2,85 |
0,200 |
0,040 |
0,570 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
5 |
19 |
2,11 |
0,100 |
0,010 |
0,211 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
6 |
20 |
1,62 |
0,050 |
0,002 |
0,081 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
7 |
30 |
1,41 |
0,033 |
0,001 |
0,047 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
8 |
50 |
1,30 |
0,020 |
0,000 |
0,026 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
9 |
100 |
1,21 |
0,010 |
0,000 |
0,001 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
10 |
200 |
1,15 |
|
0,005 |
0,000 |
0,000 |
|
||
|
|
|
|
|
|
|
|
|
|
∑ |
|
|
31,40 |
|
2,25 |
1,414 |
15,206 |
|
|
|
10a + 2,251b = 31,4 |
|
|
|
|
|
|||
2,251a + 1,414b = 15,206 . Из этой системы находим a ≈ 1,112, b = 8,973 . |
|||||||||
Искомая формула имеет вид: |
y =1,112 + |
8,973 |
. |
||||||
|
|||||||||
|
|
|
|
|
|
|
|
x |
|
