На практиці час виконання чисельного вирішення IVP x ( t 0 ) = x 0 часто переважає тривалість оцінки правої частини (RHS) f . Отже, припустимо, що всі інші операції є миттєвими (тобто без обчислювальних витрат). Якщо загальний час виконання для вирішення IVP обмежений, то це еквівалентно обмеженню кількості оцінок f до деякого N ∈
Нас цікавить лише кінцеве значення .
Я шукаю теоретичні та практичні результати, які допомагають мені вибрати найкращий метод ODE в таких умовах.
Якщо, наприклад, ми могли б вирішити IVP, використовуючи два явні кроки Ейлера шириною ( t 1 - t 0 ) / 2 або один крок шириною t 1 - t 0, використовуючи метод середини. Мені не відразу зрозуміло, який із них кращий. Для більших N можна, звичайно, також думати про багатоступеневі методи, ітераційні схеми Runge-Kutta тощо.
Я шукаю результати, подібні до тих, що існують, наприклад, для правил квадратури: Ми можемо вибрати ваг { w i } і пов'язані з ними точки { x i } такі, що правило квадратури ∑ n i = 1 w i g ( x i ) є точним для всіх многочленів g таким, що d e g ( g ) ≤ 2 n - 1 .
Отже, я шукаю верхню або нижню межі щодо глобальної точності методів ODE, враховуючи обмежену кількість дозволених оцінок RHS . Це нормально, якщо межі дотримуються лише деяких класів RHS або створюють додаткові обмеження щодо рішення (так само, як результат для правила квадратури, який утримується лише для поліномів до певної міри).
РЕДАКТУВАННЯ: Деякі основні відомості: Це стосується важких додатків у режимі реального часу, тобто результат повинен бути доступний до відомого терміну. Звідси обмеження кількості оцінок РЗС як домінуючого коефіцієнта витрат. Зазвичай наші проблеми жорсткі та порівняно невеликі.
EDIT2: На жаль, у мене немає точних вимог до термінів, але можна припустити, що буде досить малим (безумовно, <100, можливо ближче до 10). Враховуючи вимогу в реальному часі, ми повинні знайти компроміс між точністю моделей (з кращими моделями, що призводять до більш тривалого часу виконання RHS і, отже, до нижнього N ), і точністю методу ODE (з кращими методами, що вимагають вищих значення N ).