Attention! Your ePaper is waiting for publication!
By publishing your document, the content will be optimally indexed by Google via AI and sorted into the right category for over 500 million ePaper readers on YUMPU.
This will ensure high visibility and many readers!
4.8. Глобальное решение уравнений и систем 483(магазин) или куча и т. п. (см. [4]) В целом же алгоритм глобальногодоказательного решения системы уравнений организуем в виде повторяющейсяпоследовательности следующих действий:• извлечение некоторого бруса из списка L,• дробление этого бруса на потомки,• проверка существования решений в каждомиз подбрусов-потомков, по результатам которой мы− либо выдаём этот подбрус в качестве ответак решаемой задаче,− либо заносим его в рабочий список Lдля последующей обработки,− либо исключаем из дальнейшего рассмотрения,как не содержащий решений рассматриваемой системы.Кроме того, чтобы обеспечить ограниченность времени работы алгоритма,на практике имеет смысл задаться некоторым порогом мелкости(малости размеров) брусов δ, при достижении которого дальшедробить брус уже не имеет смысла. В Табл. 4.2 приведён псевдокодполучающегося алгоритма, который называется методом ветвленийи отсечений: ветвления соответствуют разбиениям исходного бруса наподбрусы (фактически, разбиениям исходной задачи на подзадачи), аотсечения — это отбрасывание бесперспективных подбрусов исходнойобласти поиска. 2Отметим, что неизбежные ограничения на вычислительные ресурсыЭВМ могут воспрепятствовать решению этим алгоритмом задачи(4.19) «до конца», поскольку могут возникнуть ситуации, когда1) размеры обрабатываемого бруса уже меньше δ, нонам ещё не удаётся ни доказать существование в нёмрешений, ни показать их отсутствие;2) размеры обрабатываемого бруса всё ещё больше δ, новычислительные ресурсы уже не позволяют производитьего обработку дальше: исчерпались выделенное время,память и т. п.2 <strong>С</strong>тандартный английский термин для обозначения подобного типа алгоритмов— «branch-and-prune». <strong>С</strong> ними тесно связаны методы ветвей и границ, широкоприменяемые в вычислительной оптимизации.
4.8. Глобальное решение уравнений и систем 483(магазин) или куча и т. п. (см. [4]) В целом же алгоритм глобальногодоказательного решения системы уравнений организуем в виде повторяющейсяпоследовательности следующих действий:• извлечение некоторого бруса из списка L,• дробление этого бруса на потомки,• проверка существования решений в каждомиз подбрусов-потомков, по результатам которой мы− либо выдаём этот подбрус в качестве ответак решаемой задаче,− либо заносим его в рабочий список Lдля последующей обработки,− либо исключаем из дальнейшего рассмотрения,как не содержащий решений рассматриваемой системы.Кроме того, чтобы обеспечить ограниченность времени работы алгоритма,на практике имеет смысл задаться некоторым порогом мелкости(малости размеров) брусов δ, при достижении которого дальшедробить брус уже не имеет смысла. В Табл. 4.2 приведён псевдокодполучающегося алгоритма, который называется методом ветвленийи отсечений: ветвления соответствуют разбиениям исходного бруса наподбрусы (фактически, разбиениям исходной задачи на подзадачи), аотсечения — это отбрасывание бесперспективных подбрусов исходнойобласти поиска. 2Отметим, что неизбежные ограничения на вычислительные ресурсыЭВМ могут воспрепятствовать решению этим алгоритмом задачи(4.19) «до конца», поскольку могут возникнуть ситуации, когда1) размеры обрабатываемого бруса уже меньше δ, нонам ещё не удаётся ни доказать существование в нёмрешений, ни показать их отсутствие;2) размеры обрабатываемого бруса всё ещё больше δ, новычислительные ресурсы уже не позволяют производитьего обработку дальше: исчерпались выделенное время,память и т. п.2 <strong>С</strong>тандартный английский термин для обозначения подобного типа алгоритмов— «branch-and-prune». <strong>С</strong> ними тесно связаны методы ветвей и границ, широкоприменяемые в вычислительной оптимизации.
484 4. Решение нелинейных уравнений и их системТаблица 4.2. Интервальный метод ветвлений и отсеченийдля глобального доказательного решения уравненийВход<strong>С</strong>истема уравнений F(x) = 0. Брус X ∈ IR n .Интервальное расширение F : IX → IR n функции F.Заданная точность δ > 0 локализации решений системы.Выход<strong>С</strong>писок НавернякаРешения из брусов размера менее δ, которыегарантированно содержат решения системы уравнений в X.<strong>С</strong>писок ВозможноРешения из брусов размера менее δ, которыемогут содержать решения системы уравнений в X.<strong>С</strong>писок Недообработанные из брусов размера более δ, которыемогут содержать решения системы уравнений в X.Алгоритминициализируем рабочий список L исходным брусом X ;DO WHILE ( (L ≠ ∅) и (не исчерпаны ресурсы ЭВМ ) )извлекаем из рабочего списка L брус Y ;применяем к Y тест существования решения,его результат обозначаем также через Y ;IF ( в Y доказано отсутствие решений ) THENудаляем брус Y из рассмотренияELSEIF ( (размер бруса Y ) < δ ) THENзаносим Y в соответствующий из списковНавернякаРешения или ВозможноРешенияELSEрассекаем Y на потомки Y ′ и Y ′′и заносим их в рабочий список LEND IFEND IFEND DOвсе брусы из L перемещаем в список Недообработанные;