|
Модераторы: Poseidon |
|
MrNull |
|
||||
Новичок Профиль Группа: Участник Сообщений: 3 Регистрация: 21.12.2007 Где: +7 906 059-05-56 Репутация: нет Всего: нет |
Так и не смог понять, разрешается ли на этом форуме помощь за деньги или нет, модеры, не ругайте сильно.
Задание: Реализовать программу на C++ поиска минимума функции методом Ньютона. Программу написать наиболее просто, без использования ссылок, указателей и т.д., так будет легче в ней разобраться. Задание и описание метода прикладываю к сообщению. Немного подробнее о функции штрафа:
Важно заметить, что в ходе аналитического решения выяснилось, что условный минимум (при наложении ограничений) и безусловный (минимум функции) совпадают, поэтому в принципе штрафные функции не нужны, так что вот так. Картинка функции с ограничениями: Срок исполнения: 23.12.2007 Оплата: от 1000 рублей и выше на твоё усмотрение Это сообщение отредактировал(а) MrNull - 21.12.2007, 01:34 Присоединённый файл ( Кол-во скачиваний: 16 ) task_and_method.rar 60,11 Kb |
||||
|
|||||
ST1 |
|
|||
Новичок Профиль Группа: Участник Сообщений: 20 Регистрация: 20.7.2005 Где: Уфа Репутация: 1 Всего: 1 |
Возьмусь за задание. Дай свой ответ
Это сообщение отредактировал(а) ST1 - 22.12.2007, 18:07 |
|||
|
||||
MrNull |
|
|||
Новичок Профиль Группа: Участник Сообщений: 3 Регистрация: 21.12.2007 Где: +7 906 059-05-56 Репутация: нет Всего: нет |
спасибо! стучись мне в аську 269-743-146 или звони +7 906 059-05-56.
ответ: глобальный минимум в точке x=2, y=3, на картинке это хорошо видно. Условный минимум (с учётом ограничений) совпадает с глобальным. |
|||
|
||||
MrNull |
|
|||
Новичок Профиль Группа: Участник Сообщений: 3 Регистрация: 21.12.2007 Где: +7 906 059-05-56 Репутация: нет Всего: нет |
Делюсь с остальными:
Программа поиска безусловного минимума функции f(x1, x2) = 2 * (x1 - 3)^2 + (x2 - 2)^2 Алгоритм расчета: методом Ньютона решается система |F1(x1, x2) = df/dx1 = 0| (1) | | |F2(x1, x2) = df/dx2 = 0|, состоящая из равенства нулю частных производных исходной функции (необходимое условие экстремума). Алгоритм метода Ньютона: Шаг 1. Задается начальное приближение для x1, x2,...xn (в нашем случае n = 2) Шаг 2. Вычисляется невязка fvec (вектор правых частей системы (1)) Проверяется сходимость (сумма модулей координат fvec меньше некоторого малого eps > 0) Если достигнута сходимость, цикл завершается. Иначе переход к следующему шагу. Шаг 3. Вычисляются элементы матрицы Якобиана системы (1): |dF1(x1, x2)/dx1 dF1(x1, x2)/dx2| J = | | |dF2(x1, x2)/dx1 dF2(x1, x2)/dx2| Шаг 4. Методом Гаусса решается система J * (dx) = - fvec, где dx = (dx1, dx2) - вектор смещений Шаг 5. Рассчитывается следующее приближение переменных: x1 = x1 + dx1, x2 = x2 + dx2 Шаг 6. Переход к следующей итерации (Шаг 1.) если не превышено максимальное число итераций Метод Гаусса решения линейных систем: Элементарными преобразованиями над строками матрицы она приводится к верхнетреугольному виду (прямой ход метода Гаусса); полученная система с треугольной матрицей решается в явном виде методом последовательного исключения неизвестных в порядке dxn,...,dx1 (обратный ход метода Гаусса)
|
|||
|
||||
iceberglife177 |
|
||||
Новичок Профиль Группа: Участник Сообщений: 1 Регистрация: 2.5.2017 Репутация: нет Всего: нет |
а функцию розенброка можно вашим методом считать? |
||||
|
|||||
Noubpoeno |
|
|||
Новичок Профиль Группа: Участник Сообщений: 32 Регистрация: 9.12.2022 Репутация: нет Всего: нет |
Модератор: Сообщение скрыто. |
|||
|
||||
Правила форума "Центр помощи" | |
ВНИМАНИЕ! Прежде чем создавать темы, или писать сообщения в данный раздел, ознакомьтесь, пожалуйста, с Правилами форума и конкретно этого раздела.
Более подробно с правилами данного раздела Вы можете ознакомится в этой теме. Если Вам помогли и атмосфера форума Вам понравилась, то заходите к нам чаще! С уважением, Poseidon, Rodman |
1 Пользователей читают эту тему (1 Гостей и 0 Скрытых Пользователей) | |
0 Пользователей: | |
« Предыдущая тема | Центр помощи | Следующая тема » |
|
По вопросам размещения рекламы пишите на vladimir(sobaka)vingrad.ru
Отказ от ответственности Powered by Invision Power Board(R) 1.3 © 2003 IPS, Inc. |