37840

Решение систем обыкновенных дифференциальных уравнений

Лабораторная работа

Архивоведение и делопроизводство

В классе неявных методов абсолютно устойчивыми являются неявный одношаговый метод Эйлера неявный одношаговый метод трапеций неявный двухшаговый метод Гира и его реализация с переменным шагом – метод Шихмана. В данной лабораторной работе изучаются следующие три наиболее часто используемые на практике численные метода: явный метод Эйлера неявный метод Эйлера неявный метод Шихмана. Явный метод Эйлера Формула интегрирования явного метода Эйлера имеет вид: 3.

Русский

2013-09-25

300 KB

64 чел.

Лабораторная работа № 3

 РЕШЕНИЕ  систем  обыкновенных

дифференциальных  уравнений

ЦЕЛЬ РАБОТЫ: изучить и программно реализовать на языке высокого уровня такие методы решения систем обыкновенных дифференциальных уравнений, как явный и неявный методы Эйлера, метод Шихмана, исследовать  их на тестовых задачах.

Элементы теории

Большая часть явлений и процессов в различных областях науки и техники описывается системами обыкновенных дифференциальных уравнений (СОДУ)

  (3.1)

с начальными условиями

   (3.2)

Такую математическую постановку задачи называют задачей Коши для системы обыкновенных дифференциальных уравнений.

Перепишем (3.1), (3.2) в векторном виде:

 (3.3)

, (3.4)

где , . Решить такую задачу Коши – значит определить траекторию , удовлетворяющую уравнению (3.3) и начальному условию (3.4).

Применение численных методов предполагает приближенное вычисление значений  в точках  отрезка . С этой целью задачу (3.3), (3.4)  заменяют разностной схемой

 (3.5)

из которой рекуррентно вычисляют приближенные значения , где через  обозначено приближенное значение .

Разностная схема связывает искомое решение  в текущий момент времени  с построенными к этому моменту времени решениями  в предыдущие моменты времени. Выбор функции  в разностной схеме определяет соответствующий метод численного решения, который носит название одношагового, если , многошагового при , явного, если   не зависит от , и неявного в противном случае.

В случае многошаговых численных методов на каждом шаге вычислительного процесса привлекается информация с нескольких предыдущих шагов. В свою очередь одношаговые методы основаны на информации, полученной лишь на одном предыдущем шаге, что, в зависимости от характера задачи, может быть как достоинством, так и недостатком. В частности, если траектория системы (3.3) меняется плавно, то более выигрышными оказываются многошаговые методы. Однако их использование требует предварительного  вычисления значений . Поэтому при разработке программы численного решения выбор расчетной формулы (3.5) осуществляется таким образом, чтобы .  Например, для нахождения  применяются одношаговые методы , для вычисления  – двухшаговые  и т. д.

Если величины  фиксированы, то схему (3.5) естественно рассматривать как систему уравнений

относительно неизвестного значения . Решение этой системы

,

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

Определение 1. Разностный метод называют абсолютно устойчивым, если малые вариации параметров СОДУ приводят к малым изменениям его решения при любых шагах интегрирования .  Если устойчивость обеспечивается лишь при некоторых ограничениях на , то такие методы называют ограниченно устойчивыми.

Для численных методов решения СОДУ важен также их порядок точности.

Определение 2. Если локальная погрешность численного метода (погрешность на шаге)

,

то порядок точности метода равен .

В научной литературе показано, что не существует абсолютно устойчивых явных линейных (как одношаговых, так и многошаговых) методов. В классе неявных методов абсолютно устойчивыми являются неявный одношаговый метод Эйлера, неявный одношаговый метод трапеций, неявный двухшаговый метод Гира и его реализация с переменным шагом – метод Шихмана.

В данной лабораторной работе изучаются следующие три наиболее часто используемые на практике численные метода:  явный метод Эйлера, неявный метод Эйлера, неявный метод Шихмана.

Явный метод Эйлера

Формула интегрирования явного метода Эйлера имеет вид:

 (3.6)

Его локальная погрешность

 (3.7)

пропорциональна , т. е. явный метод Эйлера имеет первый порядок точности.

Выбор величины шага интегрирования  в явном методе Эйлера необходимо делать, исходя из сохранения устойчивости и точности вычислений.

Условие устойчивости явного метода Эйлера

,

или

 (3.8)

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

Сформулируем условие устойчивости явного метода Эйлера в более грубой, но полезной для практического выбора шага, форме:

явный метод Эйлера устойчив, если

, (3.9)

где  – приращение переменной  на k-ом шаге интегрирования,  – максимально возможное абсолютное значение данной переменной.

Из (3.6), (3.7) следует, что

,

так как ошибка  имеет более высокий по сравнению с приращением

переменной порядок малости.

Зададимся локальной погрешностью (например, 1% от максимального значения):

.

Тогда условия соблюдения точности интегрирования, при которых  выполняются  условия устойчивости (3.9), будут иметь вид

. (3.10)

Из (3.6) следует, что

,

а, значит, неравенства (3.10) приводят к условиям

.

Недостатком этой формулы является тот факт, что , если . Поэтому на шаг интегрирования необходимо ввести ограничение сверху. Обозначим через  максимально допустимый шаг интегрирования. Тогда условия выбора шага примут вид:

.  (3.11)

Из (3.11) вытекает, что

.  (3.12)

Приведем возможную схему реализации явного метода Эйлера.

Алгоритм 3.1

  1.  Задать исходные данные: .
  2.  Задать начальные условия: .
  3.  Вычислить вектор .
  4.  Определить шаг по формулам (3.11), (3.12).
  5.  Выполнить шаг:

  1.  Вычислить .
  2.  Вывести на печать .
  3.  Если  то перейти к п. 3, иначе выйти из программы.

Неявный метод Эйлера

Разностная схема неявного метода Эйлера

  (3.13)

требует решения на каждом временном шаге алгебраической задачи

 (3.14)

относительно искомого вектора  (например, методом Ньютона).

Локальная погрешность неявного метода Эйлера

= (3.15)

по порядку величины  является такой же, как и локальная погрешность явного метода Эйлера, и противоположна по знаку.

Условие устойчивости неявного метода Эйлера

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

Для вычисления локальной погрешности  по формуле (3.15) необходимо оценить  в точке . Это можно сделать в разностной форме, используя значения

.

Следовательно

-.  (3.16)

Возможная стратегия выбора шага в неявном методе Эйлера сводится к следующему. Задается, как и в явном методе Эйлера, допустимая локальная погрешность, например, такая:

.

После выполнения шага интегрирования  вычисляются  по формуле (3.16) для всех переменных . Если хотя бы одно , то шаг  уменьшается вдвое, и вычисления повторяются с шагом  от точки . В случае , очередной шаг рассчитывается по формуле

,  (3.17)

вытекающей из соотношения для оценки локальной  погрешности.

При выборе шага интегрирования можно использовать также стратегию трех зон:

 (3.18)

В обоих стратегиях окончательный шаг интегрирования выбирается минимальным среди всех значений

.  (3.19)

Стратегия (3.17), (3.19) на первый взгляд является более предпочтительной, так как обеспечивает квазиоптимальный выбор шага. Однако при такой стратегии число возвратов (повторных расчетов с уменьшенным вдвое шагом интегрирования) будет существенно большим, чем в методе трех зон. Именно по этой причине метод трех зон является широко используемым на практике.

Возможная  схема реализации неявного метода Эйлера выглядит следующим образом.

Алгоритм 3.2

  1.  Задать исходные данные: .
  2.  Положить: .
  3.  Вычислить .
  4.  Решить методом Ньтона систему в общем случае нелинейных алгебраических уравнений (3.14) относительно .
  5.  Вычислить  по формуле (3.16).
  6.  Если , то принять  и перейти к п. 3.
  7.  Определить шаг  по формулам (3.18), (3.19) либо (3.17), (3.19).
  8.  Если , то положить .
  9.  Вывести на печать .
  10.  Выполнить сдвиг переменных и шагов интегрирования:

.

  1.  Если  то перейти к  п. 3, иначе выйти из программы.

Метод Шихмана

Среди линейных многошаговых методов абсолютной устойчивостью обладает метод Гира второго порядка

и его реализация с переменным шагом – метод Шихмана:

,  (3.20)

где

 (3.21)

Локальная погрешность метода Шихмана

 (3.22)

пропорциональна третьей производной, которую можно аппроксимировать разностной формулой

и использовать в (3.22) при оценке погрешности.

При реализации метода Шихмана необходимо учесть следующее обстоятельство. В начальный момент времени известен лишь вектор . Поэтому  первый шаг интегрирования следует выполнить неявным методом Эйлера, который также реализуется формулой (3.20) при . Второй шаг интегрирования можно было бы выполнить методом Шихмана, однако оценку погрешности по формуле (3.23) в этом случае провести еще невозможно. Значит и  второй шаг интегрирования следует выполнить неявным методом Эйлера. И лишь начиная с третьего шага интегрирования, вычисления проводятся по формулам (3.22), (3.23).

Задание

  1.  Написать, отладить и исследовать на задаче, предложенной преподавателем (перечень задач приведен ниже), программы  численного решения задачи Коши для системы обыкновенных дифференциальных уравнений (явным и неявным методами Эйлера, методом Шихмана). Вычисления выполнить для .

Задачи:

1.

2.

3.

4.

Числа  являются собственными числами симметричной матрицы . Если все , то решение системы стремится к точке покоя (-1; 1; -2). Значения  и  задаются преподавателем.

Содержание электронного  отчета

1. Тексты программ.

2. Задача, результаты  ее решения каждым методом  для .

22


 

А также другие работы, которые могут Вас заинтересовать

43264. Себестоимость продукции, ее учет и пути снижения 236.5 KB
  Затраты живого и овеществленного труда в процессе производства составляют издержки производства. В условиях товарно-денежных отношений и хозяйственной обособленности предприятия неизбежно сохраняются различия между общественными издержками производства и издержками предприятия. Общественные издержки производства - это совокупность живого и овеществленного труда, находящая выражение в стоимости продукции. Издержки предприятия состоят из всей суммы расходов предприятия на производство продукции и ее реализацию.
43265. Процесс сравнения двух четырёх разрядных чисел 132.5 KB
  Подсказанная Максимом схема. Схема сравнения или компаратор от англ compre сравнивать операционный узел ЦВМ предназначенный для сравнения двух чисел А и В. Схема.
43266. Проектирование и реализация информационно-поисковой системы Виртуальный Деканат с помощью программы MySQL-Front 458 KB
  База данных представляет собой структурированную совокупность данных. Эти данные могут быть любыми - от простого списка предстоящих покупок до перечня экспонатов картинной галереи или огромного количества информации в корпоративной сети. Для записи, выборки и обработки данных, хранящихся в компьютерной базе данных, необходима система управления базой данных, каковой и является ПО MySQL. Поскольку компьютеры замечательно справляются с обработкой больших объемов данных, управление базами данных играет центральную роль в вычислениях. Реализовано такое управление может быть по-разному - как в виде отдельных утилит, так и в виде кода, входящего в состав других приложений.
43267. СУЩНОСТЬ, ФОРМЫ И ВИДЫ ПРЕДПРИНИМАТЕЛЬСКОЙ ДЕЯТЕЛЬНОСТИ 176 KB
  Неотвратимо осуществляются экономические реформы. Результатом реформ является становление и развитие новых экономических, финансовых, социальных и других отношений, базирующихся на формировании рыночной экономики, в которой ведущим субъектом хозяйствования являются предприниматели (коллективные и индивидуальные).
43268. Методы и методики перевода 226 KB
  Современные технологии политика наука и конечно литература – без перевода все это было бы недоступно для многих наций. Когда в мире появилось как минимум два языка наверное тогда человек впервые столкнулся с проблемой адекватности перевода. Когда знания людей о каком-либо виде деятельности стали обобщаться превращаясь в науку перевод не стал исключением.
43269. Сестринский процесс при ревматоидном артрите 982.5 KB
  Ревматоидный артрит – аутоиммунное заболевание неизвестной этиологии, для которого характерны симметричный эрозивный артрит (синовиит) и широкий спектр внесуставных (системных) проявлений. Ревматоидный артрит является одним из весьма распространенных заболеваний, причем заболеваемость им за последние 50—60 лет увеличилась. Через 5 лет от начала заболевания, несмотря на лечение „базисными“ препаратами, 16% пациентов теряют трудоспособность, а через 20 лет – 90%, треть всех пациентов становятся полными инвалидами.
43270. Історія, види та характеристики веб-браузерів 579 KB
  Але доступ до нього неможливий без багатьох факторів одним з яких є Webбраузери. Webбраузери це програмне забезпечення для комп’ютера під’єднаного до Інтернет що дає можливість працювати з текстом малюнками або іншою інформацією на вебсторінці. Поняття веббраузера Браузер це програма що представляє в зручному для сприйняття вигляді інформацію що отримується з Інтернету.
43271. Решение задач оптимизации 730.5 KB
  Цель задания: Найти минимум функции методами прямого поиска и градиентными методами. Цель работы отработка навыков решения задач безусловной оптимизации функции нескольких переменных методами прямого поиска и отработка навыков решения задач безусловной оптимизации градиентными методами. Решена задача безусловной оптимизации функции нескольких переменных методами прямого поиска и градиентными методами. Нахождение стационарной точки Целевая функция: ; Частные производные f по x1 и x2: ∂f ∂x1=2x1 x2 –6; ∂f ∂x2=x1 2x2 –4; Приравниваем...
43272. Синтез линейной дискретной антенны с равномерным амплитудным распределением 892.5 KB
  Акустическая антенна обычно состоит из электроакустических преобразователей элементов антенны акустических экранов несущей конструкции акустических развязок амортизаторов и линий электрокоммуникаций. По способу создания пространственной избирательности антенны можно подразделять на интерференционные фокусирующие рупорные и параметрические. Интерференционные антенны можно подразделять на непрерывные и дискретны. И непрерывные и дискретные антенны подразделяют по конфигурации геометрического образования объединяющего активные...