Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Математика / Шилкин Лекции.PDF
Скачиваний:
33
Добавлен:
11.05.2015
Размер:
7.25 Mб
Скачать

 

 

 

 

ЛЕКЦИЯ 23

 

 

 

 

Точка ( x0

, y0 , λ

0 ) , удовлетворяющая системе

f x

( x, y) + λgx ( x, y) =0

 

 

 

 

 

 

 

 

 

f y( x, y) + λgy ( 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 =( gy , gx ) xx

xy

 

y

< 0, то целевая функция

′′

′′

gx

 

 

 

 

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 . Далее Lxx′ = 0, Lyy′ = 0, Lxy′ =1, gy = 3, gx = 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

xy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

xy

 

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 axi2 + bxi xi yi ,

 

 

 

 

 

 

 

 

 

 

 

 

i=1

 

 

i=1

 

i=1

 

 

 

 

 

 

 

 

 

 

 

 

S

 

 

 

n

 

+ b yi )

2

 

 

n

 

 

 

 

n

 

n

 

 

b

=

(axi

 

b =2(axi + b yi )==2 axi + bn yi

 

 

i=1

 

 

 

 

 

 

 

 

i=1

 

 

 

i=1

 

i=1

 

 

 

 

Приравнивая к нулю, получим линейную систему двух уравнений с

двумя неизвестными a и b , которая называется нормальной системой.

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

n

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

axi2 + bxi =xi yi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i=n1

 

i=1

 

i=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

axi + 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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

axi4 + bxi3 + cxi2 =

xi2 yi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i=1

 

 

i=1

 

i=1

 

 

 

 

i=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

n

 

n

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

axi3 + bxi2 + cxi =

xi yi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i=1

 

 

 

i=1

 

i=1

 

n

 

i=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

axi2 + bxi + 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

2a +

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