Форум программистов, компьютерный форум, киберфорум
Методы оптимизации
Войти
Регистрация
Восстановить пароль
Карта форума Темы раздела Блоги Сообщество Поиск Заказать работу  
 
Рейтинг 5.00/11: Рейтинг темы: голосов - 11, средняя оценка - 5.00
23 / 16 / 7
Регистрация: 27.10.2013
Сообщений: 95
1

Странный Гессиан

03.02.2015, 13:04. Показов 1935. Ответов 8
Метки нет (Все метки)

Author24 — интернет-сервис помощи студентам
Добрый день! Решал следующую задачу:
Нужно найти экстремумы функции https://www.cyberforum.ru/cgi-bin/latex.cgi?e^{x+y}-2e^{x}-e^{y} при ограничении https://www.cyberforum.ru/cgi-bin/latex.cgi?e^{x}+e^{y}-5=0.

Нашел функцию Лагранжа https://www.cyberforum.ru/cgi-bin/latex.cgi?e^{x+y}-2e^{x}-e^{y}+\lambda(e^{x}+e^{y}-5). По системе частных производных функции Лагранжа приравненных к нулю нашел, что https://www.cyberforum.ru/cgi-bin/latex.cgi?e^{x}=2; e^{y}=3; \lambda=-1. Составил матрицу квадратичной формы из вторых дифференциалов функции Лагранжа (её Гессиан). Получилась следующая картина:

https://www.cyberforum.ru/cgi-bin/latex.cgi?\begin{pmatrix}<br />
e^{x+y}-2e^{x}+\lambda e^{x}&e^{x+y}\\<br />
e^{x+y}&e^{x+y}-e^{y}+\lambda e^{y}\\<br />
\end{pmatrix}

Соответственно, в найденной, подозрительной на экстремум точке получил:

https://www.cyberforum.ru/cgi-bin/latex.cgi?\begin{pmatrix}<br />
0&6\	<br />
6&0\\<br />
\end{pmatrix}

По критерию Сильвестра, поскольку второй главный минор отрицательный (0*0-6*6=-36), то данная форма не является отрицательно знакоорпделенной и следовательно эта точка не может быть максимумом. Но, в ответе, а также при проверке через wolframalpha (команда maximize e^(x+y)-2e^x-e^y on e^x+e^y-5=0) оказывается, что это точка локального максимума. Подскажите, пожалуйста, в чем может быть загвоздка. Может ли в данном случае критерий Сильвестра дать сбой? Что в таком случае предпринимать?
0
Programming
Эксперт
94731 / 64177 / 26122
Регистрация: 12.04.2006
Сообщений: 116,782
03.02.2015, 13:04
Ответы с готовыми решениями:

Странный doc(rtf) файл (совсем странный)
Добрый день, сегодня на работе столкнулся с очень странным .doc файлом. 3 страницы - вес 45,9 МБ....

Странный float
Здравствуйте, при повторном изучении переменных я смотрел их диапазоны и у float 3.4000000000... до...

Странный факт
Здравствуйте. У меня 2 приложения, одно весит 6,5 МВ, а другое – 12. Способ отрисовки графики...

Странный глюк
Перестал работать один простой запрос, все вроде правильно сделал и раньше работало как положено....

8
Эксперт по математике/физике
10743 / 7124 / 3875
Регистрация: 14.01.2014
Сообщений: 16,305
03.02.2015, 14:20 2
Не то написал

Добавлено через 10 минут
Ошибок не нашел. Сам по себе критерий Сильвестра не может давать сбой, так как очевидно, что матрица Гессе в этой точке является знаконеопределенной с собственными значениями https://www.cyberforum.ru/cgi-bin/latex.cgi?\pm 6.
1
23 / 16 / 7
Регистрация: 27.10.2013
Сообщений: 95
03.02.2015, 14:27  [ТС] 3
Однако, если проверить в wolframalpha (можно и любой иной программой) то все-таки локальный максимум существует и располагается именно в этой точке. Следовательно, должна быть какая-то подоплека, позволяющая вопреки отсутствию знакоопределенности матрицы в этой точки установить то что последняя все-таки является локальным максимумом. Вопрос в том, как это сделать.
0
Эксперт по математике/физике
10743 / 7124 / 3875
Регистрация: 14.01.2014
Сообщений: 16,305
03.02.2015, 14:34 4
У Вас задача на условный экстремум, в этом случае критерии наличия экстремума усложняются - надо выписывать матрицу Гессе уже третьего порядка с производными по трем переменным!
1
23 / 16 / 7
Регистрация: 27.10.2013
Сообщений: 95
03.02.2015, 14:41  [ТС] 5
Я пробовал в качестве третьей переменной брать производную по лямбде (множитель лагранжа) и уже получал матрицу 3 на 3, но там все равно критерий Сильвестра свидетельствовал в пользу отсутствия знакоопределенности. Или я что-то напутал?
0
Эксперт по математике/физике
10743 / 7124 / 3875
Регистрация: 14.01.2014
Сообщений: 16,305
03.02.2015, 14:47 6
Лучший ответ Сообщение было отмечено Potanin как решение

Решение

Достаточно взять определитель этой новой матрицы Гессе со знаком минус. Если это число меньше 0, то это максимум. Сейчас проверил, он равен -60 в этой точке, т.е. критерий максимума выполняется!
1
23 / 16 / 7
Регистрация: 27.10.2013
Сообщений: 95
03.02.2015, 14:50  [ТС] 7
Большое спасибо, не знал о таком способе, крайне признателен! Не могли бы вы, пожалуйста, скинуть ссылку где бы можно было почитать про это подробнее? А то везде где я смотрел описывался не верный способ, которым я раньше пользовался, и нигде не было внятного доказательства или пояснения этого метода.
0
Эксперт по математике/физике
10743 / 7124 / 3875
Регистрация: 14.01.2014
Сообщений: 16,305
03.02.2015, 14:52 8
Лучший ответ Сообщение было отмечено Potanin как решение

Решение

Цитата Сообщение от Potanin Посмотреть сообщение
Я пробовал в качестве третьей переменной брать производную по лямбде (множитель лагранжа) и уже получал матрицу 3 на 3, но там все равно критерий Сильвестра свидетельствовал в пользу отсутствия знакоопределенности. Или я что-то напутал?
В случае условного экстремума критерии несколько меняются - матрица Гессе может оказаться знаконеопределенной. Посмотрите лекцию Волченко: http://math.volchenko.com/Lectures/ExtremN.pdf
1
23 / 16 / 7
Регистрация: 27.10.2013
Сообщений: 95
03.02.2015, 20:04  [ТС] 9
Еще раз большое спасибо! Пятая теорема оказалась именно тем, что нужно. Теперь, дополнительно, стало ясно как решать задачи с произвольным числом множителей Лагранжа.
0
03.02.2015, 20:04
IT_Exp
Эксперт
87844 / 49110 / 22898
Регистрация: 17.06.2006
Сообщений: 92,604
03.02.2015, 20:04
Помогаю со студенческими работами здесь

Странный конструктор
class Int { private: int intez; public: ...

Странный процессор
Скажите, пожалуйста, почему у процессора i5 3450 температура 3-го ядра всегда больше минимум на 10...

Странный проект
У меня есть проект на C#, с этим проектом работаю уже несколько месяцев и до сегодняшнего дня всё...

Странный суп
Рецепт странного супа выглядит следующим образом: Возьмем четырехзначные десятичные натуральные...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
9
Ответ Создать тему
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin
Copyright ©2000 - 2024, CyberForum.ru