Студопедия

КАТЕГОРИИ:


Архитектура-(3434)Астрономия-(809)Биология-(7483)Биотехнологии-(1457)Военное дело-(14632)Высокие технологии-(1363)География-(913)Геология-(1438)Государство-(451)Демография-(1065)Дом-(47672)Журналистика и СМИ-(912)Изобретательство-(14524)Иностранные языки-(4268)Информатика-(17799)Искусство-(1338)История-(13644)Компьютеры-(11121)Косметика-(55)Кулинария-(373)Культура-(8427)Лингвистика-(374)Литература-(1642)Маркетинг-(23702)Математика-(16968)Машиностроение-(1700)Медицина-(12668)Менеджмент-(24684)Механика-(15423)Науковедение-(506)Образование-(11852)Охрана труда-(3308)Педагогика-(5571)Полиграфия-(1312)Политика-(7869)Право-(5454)Приборостроение-(1369)Программирование-(2801)Производство-(97182)Промышленность-(8706)Психология-(18388)Религия-(3217)Связь-(10668)Сельское хозяйство-(299)Социология-(6455)Спорт-(42831)Строительство-(4793)Торговля-(5050)Транспорт-(2929)Туризм-(1568)Физика-(3942)Философия-(17015)Финансы-(26596)Химия-(22929)Экология-(12095)Экономика-(9961)Электроника-(8441)Электротехника-(4623)Энергетика-(12629)Юриспруденция-(1492)Ядерная техника-(1748)

Выбор шага и погрешность решения




Краткая характеристика методов прогноза и коррекции.

По сравнению с одношаговыми методами, методы прогноза и коррекции имеют ряд особенностей:

1) Для реализации методов прогноза и коррекции необходимо иметь информацию о нескольких предыдущих точках. Поэтому они не относятся к числу «самостартующихся» методов и начинать решение приходится с помощью какого-либо одношагового метода. По этой же причине в процессе решения дифференциальных уравнений нельзя изменять шаг интегрирования.

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

3) Применяя метод Рунге - Кутта четвёртого порядка, на каждом шаге приходится вычислять четыре значения функции, в то время как для обеспечения сходимости метода прогноза и коррекции того же порядка точности достаточно двух значений функции. Поэтому методы прогноза и коррекции требуют почти вдвое меньше машинного времени, чем методы Рунге - Кутта сравнимой точности.

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

отдельных шагах, складывающихся в суммарную ошибку, будет весьма велико. Если же наоборот - шаг выбран слишком большим, то значительной будет погрешность, обусловленная усечением рядов, и накопившаяся суммарная ошибка будет также недопустимо большой (рис.3.6).

Рис. 3.6. Выбор оптимального шага интегрирования

Кроме того, при больших значениях h возникает опасность появления цифровой неустойчивости решения. Важно различать между собой две меры погрешностей дискретизации [9],[10]:

а) локальная ошибка - это погрешность, вносимая в вычислительный процесс на каждом шаге вычислений.

б) глобальная ошибка - это разность между вычисленным и точным значениями величины на всем этапе реализации численного алгоритма, определяющая суммарную погрешность, накопившуюся с момента начала вычислений. Обычно, выбирая величину шага, стремятся, чтобы локальная ошибка на шаге была меньше некоторой заданной допустимой величины. Вообще говоря, если порядок точности метода Р, то локальная ошибка определяется выражением:

L(h) = c hp+1,

где c - некоторая постоянная; h - шаг интегрирования.

Указанное выражение с помощью θ-символики может быть записано более компактно:

L(h) = θ (hp+1).

Например, рассмотрим метод Эйлера:

Отсюда заключаем, что для метода Эйлера p =1, т.е. метод имеет первый порядок точности, а локальная ошибка определяется формулой:

L(h) = θ (h2).

Рассмотрим теперь глобальную ошибку дискретизации в конечной точке t=tN. По мере повышения требований к точности, длина шага h будет убывать, а общее из числа N, необходимое для достижения tN, будет возрастать:

Далее, глобальная ошибка E(h) может быть представлена как сумма N локальных ошибок и поэтому мы можем, округляя, записать:

E(h) = N θ(hp+1) = θ(hp).

Для метода Эйлера p =1, так что уменьшение длины шага в 2 раза уменьшает среднюю локальную ошибку примерно в 2p+1 = 4 раза. Но так как для достижения tN теперь потребуется вдвое больше шагов, то глобальная ошибка ум-еньшится лишь в 2p = 2 раза. Если используется один из методов прогноза и коррекции, то ошибка на шаге определяется величиной последнего члена в формуле коррекции. При использовании же одношаговых методов Рунге-Ку-тта, локальную ошибку не удается выразить в явной форме. Один из методов оценки этой ошибки основан на двойном счете. Если для вычисления значения искомой функции yn+1 в точке tn+1 используется шаг h, то разность между истинным и вычисленным значениями на данном шаге равна:

Если уменьшим шаг вдвое и выполнив два шага вычислим в точке tn+1, то получим:

(3.21)

Вычитая это выражение из предыдущего, найдем:

Отсюда можно найти локальную погрешность:

(3.22)

Данная формула называется правилом Рунге. Если ошибка на данном шаге слишком велика, то шаг делят пополам и вычисления повторяют вновь.

Запишем формулу (3.22) в виде:

Подставляя полученное выражение в (3.21) и полагая p=4, найдём приближенное значение ошибки на одном шаге:

(3.23)

Из (3.23) вытекает уточнённое решение:

Недостатком этого метода является то, что значение yn+1 приходится вычислять дважды, причем второй раз с помощью двух шагов. Тем не менее, эта процедура часто включается в вычислительный алгоритм для автоматического изменения шага в процессе вычислений и часто используется в методах Рунге - Кутта. Главное достоинство методов Рунге - Кутта - простота начала счета и возможность изменения величины шага в процессе вычисления. С

другой стороны, главным достоинством методов прогноза и коррекции является простота оценки ошибки на шаге.




Поделиться с друзьями:


Дата добавления: 2015-04-24; Просмотров: 1109; Нарушение авторских прав?; Мы поможем в написании вашей работы!


Нам важно ваше мнение! Был ли полезен опубликованный материал? Да | Нет



studopedia.su - Студопедия (2013 - 2024) год. Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав! Последнее добавление




Генерация страницы за: 0.008 сек.