Пред.
След.
Макеты страниц
Распознанный текст, спецсимволы и формулы могут содержать ошибки, поэтому с корректным вариантом рекомендуем ознакомиться на отсканированных изображениях учебника выше Также, советуем воспользоваться поиском по сайту, мы уверены, что вы сможете найти больше информации по нужной Вам тематике ДЛЯ СТУДЕНТОВ И ШКОЛЬНИКОВ ЕСТЬ
ZADANIA.TO
11. Сгущение сетки.Как получить требуемую точность расчета? Априорные оценки точности для этого мало полезны. Во-первых, остаточные члены выражаются через производные решения, которое до начала расчета неизвестно. Во-вторых, априорные оценки обычно являются мажорантными и могут во много раз превосходить фактическую ошибку расчета. Имеются стандартные программы численного интегрирования дифференциальных уравнений с гак называемым «автоматическим выбором шага». В них каждый шаг выбирается так, чтобы вносимая на нем погрешность не превышала заданной величины. Но при этом не учитывается, что эта погрешность в ходе дальнейших расчетов умножается на величину типа экспоненты в (18), т. е. может сильно возрасти. Кроме того, общее число шагов заранее не определено. В результате фактическая точность расчета по подобным программам обычно неизвестна. Поэтому основным практическим приемом является апостериорная оценка точности. Для ее получения расчет проводят на двух или более сгущающихся сетках и применяют правило Рунге или Рунге — Ромберга (см. главу III, п. 3). Напомним, в чем оно заключается. Вспомним априорную оценку погрешности схемы ломаных (18). Запишем ее, опуская первое слагаемое, связанное с неточным заданием начальных данных:
здесь Если сетка равномерная,
Вместо оценки точности можно погрешность (43) прибавить к численному решению, уточнив его:
но тогда вопрос о погрешности уточненного решения остается открытым. Приведенное рассуждение справедливо и в том случае, если сетки с разным числом узлов не равномерны, но их можно описать функциями При выводе оценок типа (18) старшими членами формулы Тейлора (13) пренебрегают. Если их учесть (считая правую часть уравнения непрерывно дифференцируемой достаточное число раз), то погрешность выразится суммой, где последующие слагаемые содержат более высокие степени Правило Рунге применимо для сеток с любым отношением шагов При этом как для равномерных, так и для квазиравномерных сеток условие совпадения узлов выполняется. В тех узлах, которые являются общими для нескольких сеток, можно уточнить
Рис. 44. Однако можно уточнить функцию во всех узлах самой подробной сетки, если немного усложнить вычисления. Например, для двух нижних сеток на рис. 44 это делается так. Используем совпадающие узлы сеток для определения поправок к значениям функции
Значение поправок в остальных узлах найдем простейшей интерполяцией. Для равномерных или квазиравномерных сеток можно положить
Затем вычислим уточненные значения
Этот способ легко обобщается на произвольное число сеток. Такое уточнение выгодно для специальных схем третьего типа, имеющих невысокий порядок точности; выполнить уточнение обычдо проще, чем составить специальную схему высокого порядка точности. Примеры применения правила Рунге даны в таблице 18 (п. 5) и таблице 19 (п. 6), содержащих численное решение задачи
В таблице 18 интегрирование выполнено по схеме ломаных (15), и для уточнения использованы сетки с
|
1 |
Оглавление
|