ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 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