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


 

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

50031. Инструментальные возможности программы Corel Draw 167 KB
  Это также наиболее известный из графических программных продуктов корпорации Corel которая наряду с dobe Corportion является ведущим производителем программных продуктов для компьютерной графики. Достоинствам продуктов Corel Corportion является разработка нескольких миллионов готовых изображений причем каждая линия в них поддается редактированию. В Corel Drw существуют не только мощные средства векторного редактирования но и средства верстки многостраничных документов а также подготовки их как в печатном так и в электронном виде.
50032. Измерение параметров индуктивности в цепи переменного тока 255 KB
  Цель работы: Определение импеданса сдвига фаз и измерение индуктивности на разных частотах в резистивно-индуктивной цепи. При работе на переменном токе с реактивными элементами в цепи индуктивность емкость следует обязательно учитывать их реактивный характер проводимости. Кроме того реактивные...
50033. Перевірка правил Кірхгофа 133.5 KB
  Мета роботи: перевірити правила Кірхгофа для кола постійного струму. Теоретичні пояснення правил Кірхгофа а також їх практичне використання для розрахунку розгалужених електричних кіл показані в розділі 3. Застосуємо перше правило Кірхгофа до вузла В...
50034. ИЗМЕРЕНИЕ РАЗРЕШАЮЩЕЙ СПОСОБНОСТИ ОБЪЕКТИВОВ 315 KB
  Как следствие фокусное расстояние объектива зависит от длины световой волны и если для одной длины волны изображение хорошо сфокусировано то для других длин волн хорошей фокусировки не наблюдается. Если как это обычно бывает оправа объектива круглая то изображение светящейся точки имеет вид круглого пятна окруженного концентрическими светлыми и темными кольцами рис. Способность объектива создавать раздельные изображения близко расположенных мелких деталей называется разрешающей способностью объектива. Чем меньше угол  тем ближе...
50035. Юридическая психология. Учебно-методический комплекс 677.5 KB
  Цель дисциплины – психологическая подготовка юриста к профессиональной деятельности, формирование эффективных приемов работы с людьми и овладение методами профессионально значимого самопознания и саморазвития личности.
50036. ОПРЕДЕЛЕНИЕ КОЛИЧЕСТВА ВРЕДНЫХ ГАЗОВ В ВОЗДУХЕ ПРОИЗВОДСТВЕННЫХ ПОМЕЩЕНИЙ 1.08 MB
  Приводятся величины предельно допустимых концентраций ПДК вредных газов выше которых следует применять средства защиты: респираторы с противогазовыми патронами если концентрация вредного газа не превышает 10 ПДК или противогазы если концентрация вредного газа в воздухе выше 10 ПДК. ИЗУЧЕНИЕ И ВЫБОР СРЕДСТВ ИНДИВИДУАЛЬНОЙ ЗАЩИТЫ ОРГАНОВ ДЫХАНИЯ Цель работы: Изучить назначение устройство и область применения средств индивидуальной защиты органов дыхания СИЗОД. Порядок выбора СИЗОД их классификация и требования к применению...
50037. ОПРЕДЕЛЕНИЕ РАДИУСА КРИВИЗНЫ ЛИНЗЫ С ПОМОЩЬЮ КОЛЕЦ НЬЮТОНА 140 KB
  Кольца Ньютона. Классическим примером полос равной толщины являются кольца Ньютона. 1а видно что толщина воздушного зазора δ связана с радиусами наблюдаемых колец rk и радиусом кривизны линзы R следующим образом: 2 Учтем что δ R2 1 и пренебрегая этим слагаемым в формуле 2 получим: 3 С учетом 3 выражение 1 будет иметь вид: 4 Интерференционные максимумы светлые кольца получаются при условии если на разности хода укладывается четное число полуволн целое число длин волн: 5 где k = 0 1 2. порядок интерференции или...