Немонотонная сходимость в задаче с фиксированной точкой


13

Фон

Я решаю вариант уравнения Орнштейна-Цернике из теории жидкости. Абстрактно, задачу можно представить как решение задачи с неподвижной точкой , где A - интегроалгебраический оператор, а c ( r ) - функция решения (функция прямой корреляции OZ). Я решаю с помощью итерации Пикара, где я даю исходное пробное решение c 0 ( r ) и генерирую новые пробные решения по схеме c j + 1 = α (Aс(р)знак равнос(р)Aс(р)с0(р) где α - это настраиваемый параметр, который управляет сочетанием c и A c, используемым в следующем пробном решении. Для этого обсуждения предположим, что значение α неважно. Повторяюпока итерации сходится с точностью до желаемой толерантности, ε : Δ J + 1 d г | c j + 1 ( r ) - c

сJ+1знак равноα(AсJ)+(1-α)сJ ,
αсAсαε В моем варианте задачи A зависит от параметра λ , и мой вопрос заключается в том, как сходимость A c = c зависит от этого параметра.
ΔJ+1dр|сJ+1(р)-сJ(р)|<ε ,
AλAсзнак равнос

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

Ключевые вопросы

Имеет ли немонотонная сходимость в итерационных решениях задач с фиксированной точкой какое-либо особое значение? Означает ли это, что моя итерационная схема находится на грани нестабильности? Самое главное , должна ли немонотонная конвергенция вызывать у меня подозрение, что «конвергентное» решение не является хорошим решением проблемы с фиксированной точкой?

Ответы:


1

ИксИксзнак равное(Икс)Икс*еИкс(Икс*)αα<1Икс*

  1. λ

  2. Если ваше решение сошлось в рамках правильно установленного относительного допуска, который также учитывает небольшие числа, то оно имеет.


Можете ли вы уточнить свой второй пункт?
Endulum

|ИксJ+1-ИксJ||ИксJ|εε
Используя наш сайт, вы подтверждаете, что прочитали и поняли нашу Политику в отношении файлов cookie и Политику конфиденциальности.
Licensed under cc by-sa 3.0 with attribution required.