Файл: КПрМА Практикум 7.pdf

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

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

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

Добавлен: 15.03.2025

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

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

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

Практикум 7. Экстремумы функции двух переменных

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

Продолжительность работы - 2 часа.

Оборудование – работа выполняется в компьютерном классе с использованием пакета MATLAB.

Порядок выполнения

1.Знакомство со справочным материалом по математике

2.Знакомство со справочным материалом по пакету MATLAB.

3.Изучение примеров.

4.Самостоятельное выполнение упражнений. При выполнении упражнений в случае сообщения системы об ошибке рекомендуется найти и исправить ошибку самостоятельно; однако, если после многократных попыток сделать

это не удается, то можно и нужно проконсультироваться с преподавателем.

P.S. Отчитываться перед преподавателем о выполнении упражнений не нужно. Однако, следует учесть, что их выполнение – залог успешного написания контрольной работы по модулю, поскольку контрольная работа составлена из аналогов упражнений.

Справочный материал по математике

1. Экстремумы функции. Функция z = f (x, y) имеет максимум (минимум) в точке M 0 ( x0 , y0 ) , если существует такая окрестность точки M0 , для всех точек M которой, отличных от M0 , выполняется неравенство f (M0 ) f (M ) ( f (M0 ) f (M ) ).

Необходимое условие экстремума. Если дифференцируемая функция z = f (x, y)

достигает экстремума в точке M0 , то в этой точке fx (M 0 ) = 0 , f y (M 0 ) = 0 .

Равенства fx = 0 , f y = 0 могут выполняться не только в точках экстремумов, но и в точках, которые не являются точками экстремума. Точки, в которых выполняются условия fx = 0 , f y = 0 называют стационарными.

Отбор точек экстремума среди стационарных точек связан с исследованием знака второго дифференциала функции. Этот способ применим только в том случае, если функция z = f (x, y) дважды дифференцируема в некоторой окрестности стационарной

точки и все ее производные непрерывны в этой точке.

 

 

 

Достаточное условие экстремума. Пусть S - стационарная точка.

 

 

(1) Если d 2 f (S, x, y) 0

при любых x, y одновременно не равных нулю, то S

- точка

минимума. (Второй

дифференциал

удовлетворяет

этому

условию,

если

 

2

 

 

 

 

 

 

 

AC B

 

0 , A 0 . Здесь A = fxx (S) ,

B = fxy (S ) , C = f yy (S ) ).

 

 

 

(2) Если d 2 f (S, x, y) 0

при любых x, y одновременно не равных нулю, то S

- точка максимума. (Второй

дифференциал

удовлетворяет

этому

условию,

если

AC B2 0 , A 0 ) .

 

 

 

 

 

 

1


(3) Если d 2 f (S, x, y) при одних

x, y принимает положительные значения, а

при других - отрицательные, то S

не является точкой экстремума. (Второй

дифференциал удовлетворяет этому условию, если AC B2 0 ).

(4) Если d 2 f (S, x, y) 0 при любых x, y или d 2 f (S, x, y) 0 при любых

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

(Второй дифференциал удовлетворяет этим условиям, если AC B2 = 0 ).

2. Линии уровня функции. Рассмотрим способ описания функции посредством

линий

уровня. Построим

сечения поверхности z = f (x, y)

горизонтальными

плоскостями z = c

(здесь c

- числа). Эти сечения представляют собой некоторые

кривые

– каждая

кривая

соответствует определенному числу.

Затем нанесем

(спроектируем) полученные кривые на плоскость Oxy . Получим линии уровня функции. Заметим, что линии уровня применяют на географических картах для изображения

рельефа местности.

Линии уровня могут помочь в локализации точек максимума и минимума функции.

Справочный материал по пакету MATLAB

Построение линий уровня функции

Для построения линий уровня служит команда contour. Базовый формат вызова функции:

contour(X,Y,Z)

Здесь X,Y – массивы координат узлов сетки, в которых определены аппликаты Z точек поверхности.

Другой вариант вызова:

contour(X,Y,Z,k)

Аргумент k задает число уровней (высот по оси z ), для которых должны быть построены линии уровня.

Еще один вариант вызова:

contour(X,Y,Z,v)

Аргумент v – вектор, который задает набор уровней (высот по оси z ), для которых должны быть построены линии уровня. Если какие-либо компоненты этого вектора не принадлежат диапазону [zmin , zmax ], то соответствующие линии уровня не существуют и

не строятся.

Чтобы рисунок стал более информативным, можно сделать оцифровку линий – записать около каждой линии уровня соответствующее значение z . Для автоматической оцифровки линий служит функция clabel(C,h). Ее входные параметры – матрица изолиний и вектор указателей на линии уровня. Их нужно предварительно сформировать, добавив к любому из форматов вызова функции contour выходные аргументы (например, [C,h]=contour(X,Y,Z,k)). Оцифровка линий уровня позволит оценить значения функции в точках экстремума.

u=-4:0.05:4; v=-2:0.05:2; [X,Y]=meshgrid(u,v); Z=X.^2+4*Y.^2; [C,h]=contour(X,Y,Z,10); clabel(C,h)

2


2

 

 

 

 

 

 

 

 

 

 

 

14.

 

 

17.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.4545

 

 

 

4545

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5455

 

 

23

 

 

 

 

 

 

 

 

 

 

17

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.5455

11

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

14

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

20

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

6364

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

1.5

 

 

2727

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

23

 

 

 

 

 

.6364

 

 

 

 

8

 

 

 

 

2727

 

 

 

 

 

 

 

 

11

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

20

3636

 

 

 

 

.72727

 

 

 

 

72727

 

3636

 

 

 

 

 

 

 

 

 

8

 

5.

 

 

 

14

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

81818

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

11

 

 

 

 

 

 

 

5455

 

 

 

 

 

 

 

 

 

 

 

 

 

5455

 

 

 

 

 

 

 

 

.81818

 

 

 

2

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6364

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

14

.

 

 

.90909

 

 

90909

 

5

 

 

 

 

 

0.5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

 

81818

8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

72727

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

11

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.5455

 

 

 

 

 

 

. .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

636472727

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

-0.5

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

90909

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

90909

 

 

. 6364

 

14

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

11

 

 

 

 

 

 

 

 

 

 

 

5.

 

 

 

 

 

 

 

 

 

 

 

-1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

14

 

 

 

81818

 

 

 

.81818

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

 

20

 

 

5455

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

11.

8.

 

 

 

 

8

 

 

.3636

 

-1.5

23

 

 

 

 

 

 

72727

 

 

 

 

 

 

 

 

 

 

 

.3636

 

 

 

6364

 

 

 

.6364

 

 

 

20

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

11

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5455

 

 

 

 

2727

 

17.

 

 

 

 

 

 

 

 

 

14

 

23

2727

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4545

 

14.

 

 

 

 

 

 

.4545

 

 

 

 

-2

 

 

 

 

 

 

 

5455

 

 

 

 

 

 

17

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-4

 

 

 

 

-3

-2

 

 

-1

0

 

 

1

 

 

2

 

3

 

 

4

Упражнения

Упражнение 1

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

а) z = sin xsin y sin(x + y) ( 0 x ; 0 y );

б) z = (5x + 7 y 25)e( x2 +xy+ y2 ) .

Упражнение 2

Если мы хотим выяснить, является ли стационарная точка точкой экстремума, нам нужно посмотреть, какие по знаку значения принимает второй дифференциал функции (п. 1 Справочного материала по математике). С другой стороны, представление об экстремумах функции можно получить, построив ее линии уровня. Зададимся вопросом: «Есть ли взаимосвязь между знаком второго дифференциала функции и видом ее линий уровня вблизи стационарной точки?». Исследуем вопрос экспериментально (на конкретных примерах). Результаты занесем в таблицу, проанализируем, сформулируем подмеченные закономерности.

S -

Знак

Изображение графика

Изображение линий

стационарная

d 2 f (S, x, y)

функции z = f (x, y)

уровня функции

точка

вблизи

 

z = f (x, y)

 

 

 

 

стационарной

 

 

 

точки

 

 

 

 

 

 

z = 0, 25x2 + 0, 09 y2

 

 

 

 

 

 

 

 

 

 

 

 

 

3


z = x2 ( y 2)2

z = (x y +1)2

z = x2 2xy2 + y4 y5

z = x2 y3 (6 x y)

z = 2x4 + y4 x2 2 y2

 

Список литературы и информационных ресурсов

1.

Сборник задач по математике для втузов [Текст]: Учеб. пособие

для втузов: В 4-х ч. Ч. 2: [Введение в анализ; Дифференциальное и интегральное исчисление функций одной переменной; Дифференциальное исчисление функций нескольких переменных; Кратные интегралы; Дифференциальные уравнения] / Под ред. А.В. Ефимова, А.С. Поспелова. - 5-е изд., перераб. и доп. - М.: Физматлит, 2009.

2.

В.Г.Потемкин "Введение в Matlab" (v 5.3)

http://matlab.exponenta.ru

3.

Мещеряков В.В. Задачи по математике с MATLAB&SIMULINK

– М.: ДИАЛОГ-МИФИ, 2007

4.

Амос Гилат. MATLAB. Теория и практика. 5-е изд./ Пер. с англ.

Смоленцев Н.К. – М.:ДМК Пресс, 2016.

5.

http://matlab.exponenta.ru

4