66537

Численное решение задачи Коши для ОДУ

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

Информатика, кибернетика и программирование

Цель: Решить ОДУ методом Адамса(3-х шаговым неявным) и методом Рунге-Кутта 4-го порядка (классическим). Результаты предоставить графически. Задание 1 Условие: Решить методом Адамса (3-х шаговым неявным). Результат представить графически...

Русский

2014-08-22

160 KB

2 чел.

Выполнил: Марудо А.В., 2 курс, 3 группа

Проверил: Шапочкина Ирина Викторовна

Лабораторная работа #4(вариант #47)

Тема: Численное решение задачи Коши для ОДУ.

Цель: Решить ОДУ методом Адамса(3-х шаговым неявным) и методом Рунге-Кутта 4-го порядка (классическим). Результаты предоставить графически.

Задание 1

Условие: Решить методом Адамса (3-х шаговым неявным). Результат представить графически:

;          ;          ;

Ход работы: 

Теория метода Адамса-Моултона:

проинтегрировав по промежутку

интерполируя f(x, y(x))

делаем замену  в соответствии с которой

где

В 3-х шаговом методе интерполяционный многочлен будет иметь 2-ю степень.

Ординатная формула будет иметь следующий вид:

Код программы:

//вычислим решение ОДУ в 3-х начальных точках с помощью классического метода Рунге-Кутта 4-го порядка точности

//применим ординаторную формулу

while x<=b do

  begin

   y:=y+dx/24*(55*f(x,f0)-59*f(x-dx,f1)+37*f(x-2*dx,f2)-9*f(x-3*dx,f3));

   y:=y+dx/24*(9*f(x,y)+19*f(x-dx,f3)-5*f(x-2*dx,f2)+f(x-3*dx,f1));

   x:=x+dx;

   Series1.AddXY(x,y);

   f0:=f1;

   f1:=f2;

   f2:=y;

  end;

Полученные результаты:

Вывод:

Построен график y(x) решений ОДУ.

Порядок шаговости метода Адамса-Моултона на единицу ниже порядка его точности. Т.е. метод является методом 4-го порядка точности.

Задание 2

Условие: Решить методом Рунге-Кутта 4-го порядка (классическим). Результат представить графически:

;          ;          ;

Ход работы:

;

Теория метода:

1. Вычисляем значение коэффициента наклона касательной в начальной точке:

;

2. Вычисляем значение коэффициента наклона в середине интервала, используя предыдущее приближение:

;

3. Уточняем значение коэффициента наклона в центре интервала, используя полученное значение :

;

4. С помощью значения  находим приближенное значение функции в конце интервала и коэффициент наклона касательной в конечной точке:

;

5. «Усредняем» полученные коэффициенты и находим уточненное значение искомой функции

;

Код программы:

 begin

  x:=0;

  xx:=2;

  y0:=1;

  n:=100;

  dx:=(xx-x)/n;

  repeat

   Form1.Series1.AddXY(x,y);

   k[1]:=f(x,y);

   k[2]:=f(x+1/2*dx,y+1/2*dx*k[1]);

   k[3]:=f(x+1/2*dx,y+dx*k[2]/2);

   k[4]:=f(x+dx,y+dx*k[3]);

   y:=y+(k[1]+2*k[2]+2*k[3]+k[4])/6*dx;

   x:=x+dx;

  until x>=2;

 end;

Полученные результаты:

Вывод:

Построен график y(x) решений ОДУ.

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

БГУ

Физический факультет

2011/2012 учебный год

Минск

PAGE   \* MERGEFORMAT 2


 

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

26487. Особенности оформления писем, предаваемых электронной почтой 49 KB
  Особенности оформления писем предаваемых электронной почтой. Напомним вначале что электронным письмом называют документ передаваемый по каналам электронной почты. Адрес в системе электронной почты состоит из имени электронного почтового ящика которое обычно совпадает с регистрационным именем пользователя и домена который описывает место компьютер или локальную систему где этот электронный ящик на ходится. В целом требования к оформлению текста документов посылаемых электронной почтой аналогичны нормам изложенным в главе 3 п.
26488. Основные понятия делопроизводства 51 KB
  Организация работы с документами организация документооборота хранения и использования документов в текущей деятельности учреждения предприятия. Службой документационного обеспечения управления ДОУ называется структурное подразделение на которое возложены делопроизводственные операции регистрация контроль исполнения хранение использование документов и т. Структурными подразделениями службы ДОУ в зависимости от уровня организации и объема документов являются: управление делами; управление делопроизводством; канцелярия; отдел ДОУ;...
26489. Бланки документов и их оформление 49 KB
  Бланк стандартный лист бумаги на котором заранее воспроизводится информация об организации авторе от имени которого издается документ. Для организации ее структурного подразделения должностного лица устанавливают следующие виды бланков документов: общий бланк; бланк письма; бланк конкретного вида документа кроме письма. Реквизиты общего бланка документа: герб для организаций имеющих на это право; эмблема организации при наличии герба не проставляется; наименование вышестоящей организации если она имеется; наименование...
26490. ОТВЕТЫ К ПРОМЕЖУТОЧНОЙ АТТЕСТАЦИИ 10 КЛАССА ПРОФИЛЬНОГО КУРСА «ОФИСНЫЕ ТЕХНОЛОГИИ» 34.5 KB
  При адресовании документа должностному лицу инициалы указываются перед фамилией. При адресовании документа физическому лицу указывают фамилию получателя затем почтовый адрес. При подписании документа несколькими должностными лицами их подписи располагают одна под другой в последовательности соответствующей занимаемой должности. При подписании документа несколькими лицами равных должностей их подписи располагают на одном уровне.
26491. Многокритериальные задачи принятия решения 18.13 KB
  Смысл обоих подходов состоит в том что один из критериев оценки альтернатив переводится в ограничение. В ряде случаев можно использовать отношение двух указанных критериев. Третий подход к синтезу критериев стоимости и эффективности приводит к построению паретовского множества. Парето развивая исследования эджварда ввел в экономику понятия оптимальности для случая нескольких критериев.
26492. Принятие решений в задачах с детерминированными целочисленными параметрами 24.47 KB
  Первая категория задачи с неделимостями. Вторая категория комбинаторные задачи. задачи теории расписаний упорядочение планирование согласование. Третья категория задачи сводящиеся к задачам дискретного программирования.
26493. Основные понятия теории расписаний 29.8 KB
  Задачи теории расписаний делятся на детерминированные и стохастические. К детерминированным задачам теории расписаний относятся задачи упорядочения планирования и согласования. В этом случае задачи детерминированного календарного планирования сводятся к задачам упорядочения. В некоторых классификациях к задачам теории расписания могут быть отнесены например задачи распределения в которых множество работ с заданными временными характеристиками необходимо распределить по приборам у которых заранее установлены параметры производительности.
26494. Применение метода динамического программирования в задачах принятия решений 26.55 KB
  Концептуально динамическое программирование применяется для анализа систем которые характеризуются следующими признаками: процесс функционирования системы включает последовательные этапы текущие этапы i конечный этап m. предполагается что для системы выполняется принцип отсутствия последействия. Суть этого принципа заключается в том что состояние Si зависит только от состояния системы на предыдущем этапе то есть на Si1 а так же зависит от управляющего воздействия Ui. И не зависит от предыдущих состояний системы и предыдущих...
26495. Основные типы вероятностных задач и критериев оценки решения 30.14 KB
  Например допустим рассматривается детерминированная система на вход которой через равные промежутки времени Т1 поступают работы.ожидания времени простоя на стоимость 1ой единицы времени их работы зарплата отнесенная к суммарному фонду рабочего времени. 2 Математический аппарат используемый при разработке модели ПР Для конструирования вероятностных моделей ПР примем аппарат случайных процессов: Процесс называется случайным если для каждого момента времени его состояние представляет собой случайную величину. Если переходы между...