Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Высшая математика - 2.Уч.-практ.пособие.pdf
Скачиваний:
93
Добавлен:
20.02.2016
Размер:
1.66 Mб
Скачать

Пример

 

1.14.

Исследовать

 

 

на

экстремум

функцию

z = x3 + y3 3xy .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение. Находим частные производные первого и второго

порядков:

zx

= fx(x, y) = (x3 + y3 3xy)x = 3x2 3y ;

 

 

 

 

 

zy = fy(x, y) = (x3 + y3 3xy)y = 3y2 3x ;

 

 

 

 

 

 

z′′xx = (3x2 3y)x = 6x ;

 

 

 

 

 

 

 

 

 

 

 

 

 

z′′xy = (3x2 3y)y = −3 ;

z′′yy = (3x2 3y)y = 6y .

 

 

 

 

 

 

 

 

 

 

2

3y =

0,

 

 

= x

2

,

 

Составляем систему (1.11):

3x

 

 

y

 

 

из ко-

 

 

 

 

 

 

 

= y2

 

 

 

 

 

 

3y2 3x =

0

 

x

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

торой находим две стационарные точки M0 (0,0)

и M1(1,1) .

 

 

 

 

 

 

 

 

 

 

В точке M0

имеем A = 6x

 

= 0,

C = 6y

 

= 0 , B = −3,

 

 

 

 

(0;0)

 

 

 

 

 

(0;0)

 

 

 

 

 

 

Δ= AC B2

= 0 (3)2 = −9 . Так как

< 0 ,

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

точке M0 нет. В точке M1 имеем:

A = 6x

 

 

= 6 , C = 6y

 

 

= 6 ,

 

 

 

 

 

 

 

 

 

(1;1)

 

 

 

(1;1)

B = −3, Δ= AC B2 = 36 (3)2 = 27 .

 

 

 

 

 

 

 

 

 

 

 

Так как

< 0 и

A > 0 , то в точке

M1

функция имеет мини-

мум, который равен

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

zmin = f (1,1) =13 +13 3 1 1 = −1. Ответ: zmin = f (1,1) = −1.

1.4. Метод наименьших квадратов

В различных экономических и других практических задачах часто возникает необходимость установления аналитической зависимости между интересующими переменными, которые заданы, например, в виде статистических данных за определенный период. Одним из распростра-

26

ненных способов решения подобных задач является метод наименьших квадратов, который основывается, по сути, на нахождении экстремума функций нескольких переменных.

1.4.1. Понятие эмпирической формулы

Важное значение имеет следующая задача: требуется установить вид функциональной зависимости между двумя переменными

величинами x и

y по результатам n экспериментальных измере-

ний, приведенных в таблице 1.1:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

x1

x2

 

xi

xn

 

 

y

y1

y2

 

yi

yn

Табл.1.1.

Иначе говоря, требуется выразить зависимость между x и y

аналитически, т.е. указать формулу, связывающую между собой соответствующие значения переменных. Формулы, служащие для аналитического представления опытных данных, принято называть

эмпирическими формулами.

Следует заметить, что подбор эмпирической формулы не ставит задачу разгадать истинный вид зависимости эта задача ма-

тематически неразрешима. Ставится задача подобрать формулу, в каком-то смысле наилучшим образом отображающую полученные результаты. Для этой цели применяются различные методы. Мож-

но построить многочлен, принимающий в данных точках xi значения yi , приведенные в табл. 1.1. Достоинство этого метода в том,

что полученная формула точно воспроизводит заданные значения. Такого рода формулы называются интерполяционными. Здесь по

двум точкам строится прямая, по трем точкам

парабола, по n

точкам

многочлен степени (n 1) , т.е. с ростом n степень мно-

гочлена растет.

 

Можно разбить опытные данные на тройки точек, а затем по каждой тройке точек строить параболу. Получится функция, «склеенная» из парабол. Достоинства такого способа состоят в том, что, во-первых, данные табл.1.1 отображаются точно, во-вторых, степень многочлена невысокая (вторая). Недостатком является громоздкая запись и недифферен-

27

цируемость полученной функции в местах «склейки». Такой подход называется сплайн-интерполяция.

Третий способ заключается в следующем. Исходя из некоторых теоретических или практических соображений (например, конфигурации расположения точек на координатной плоскости) подбирается наиболее простая формула, которая дает наилучшее совпадение с опытными данными. Наиболее типичными в экономических исследованиях являются формы зависимостей в виде ли-

нейной

(y = a x +a ) , квадратичной ( y = a x2

+a x +a ) , степен-

 

0

1

0

 

 

1

2

ной (y = a xa1 ) ,

логарифмической ( y = a

log

a

x) ,

показательной

 

0

0

 

 

 

 

 

 

 

 

1

 

 

(y = a

a x ) и обратно-пропорциональной

(y = a

+ a1 ) функций.

0

1

 

 

 

 

0

x

 

 

 

 

 

 

 

Слова «наилучшее совпадение» понимаются здесь в том смысле, что из данного множества формул вида y = f (x, a0 , a1,..., am )

наилучшей считается та, для которой сумма квадратов отклонений наблюдаемых значений yi от вычисленных по формуле значений

y = f (xi , a0 , a1,..., am ) является наименьшей. Описанный третий

способ построения эмпирической формулы называется методом наименьших квадратов, а вычисленные путем решения задачи

n

S(a0 , a1,...am ) = (yi f (xi , a0 , a1,...am ))2 min

i=1

значения параметров a

0 , a 0

,..., a

0

задают наилучшую в смысле ме-

0

1

 

m

 

 

 

 

тода наименьших квадратов формулу

y = f (x, a

0 , a 0

,...a

0 ) . Исходя

 

 

 

 

0

1

 

m

из необходимых условий экстремума функций многих переменных, минимум функции S(a0 , a1,..., am ) будет в тех точках, где частные произ-

водные

 

S

,

 

S

, …

S

обращаются в нуль. Таким образом, полу-

 

 

 

 

a

 

 

a

a

 

 

 

 

0

 

1

 

m

 

чается система

(m +1)

уравнений с (m +1) неизвестными

 

S = 2[ f (xi , a0 ,..., am ) yi ] f (xi , a0 ,..., am ) = 0, j = 0,1,2,..., m,

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a j

 

a j

 

i=1

 

 

 

 

 

именуемую нормальной системой метода наименьших квадратов.

28

1.4.2. Выравнивание экспериментальных данных по прямой

Пусть для данных табл. 1.1 из теоретических или практических соображений известно, что эмпирическую функцию следует искать

в виде y = a1x + a0 . Тогда наилучшие значения параметров a1 и a0

являются решением нормальной системы метода наименьших квадратов, которая в данном случае имеет вид:

 

n

2

 

n

n

 

a1

xi

+ a0

xi

= xi yi,

 

 

i=1

 

 

i=1

i=1

(1.12)

 

n

 

 

n

 

a1

xi

+ a0n = yi .

 

 

i=1

 

 

i=1

 

 

Пример 1.15. Рост среднемесячной зарплаты сотрудников фирмы за 5 лет отражен в следующей таблице:

Годы

1

2

3

4

5

 

 

 

 

 

 

Средняя зарплата (у.е.)

235

250

270

292

300

 

 

 

 

 

 

Полагая эту зависимость линейной, установить ее аналитическую форму.

Решение. Уравнение искомой зависимости будем искать в виде y = a1x + a0 . Коэффициенты системы (1.11) удобно вычислять с

помощью расчетной таблицы, в столбцы которой заносим исходные данные, а также требуемые для коэффициентов системы (1.12) промежуточные вычисления:

i

x

y

x 2

x y

 

 

i

i

i

i i

 

1

1

235

1

235

 

 

 

 

 

 

 

2

2

250

4

500

 

 

 

 

 

 

 

3

3

270

9

810

 

 

 

 

 

 

 

4

4

292

16

1168

 

 

 

 

 

 

 

5

5

300

25

1500

 

 

 

 

 

 

 

Итого:

15

1347

55

4213

Табл.1.2.

 

 

 

 

 

29