10662

Интегрирование дифференциальных уравнений второго порядка методом Рунге-Кутта

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

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

Лабораторная работа 12 Интегрирование дифференциальных уравнений второго порядка методом РунгеКутта. Цель работы. Научиться решать дифференциальное уравнение второго порядка путем преобразования его к системе двух уравнений первого порядка с последующ

Русский

2013-03-30

310 KB

23 чел.

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

Интегрирование дифференциальных уравнений

второго порядка методом Рунге-Кутта.

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

 Теоретические положения .  Пусть дано дифференциальное уравнение второго порядка

                          (1)

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

,                                 

.                                (2)

Требуется найти его решение    на отрезке  ,  где

,  ,                  (3)

.

Если ввести  вспомогательную функцию   такую, что

 

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

то уравнение (1) можно преобразовать к системе  уравнений вида:

                          (4)

Согласно методу Рунге-Кутта система (4) решается следующим образом:  cчитая известными величины    и  ,

- найдем две группы чисел

,

,

,                                (5)

,

и

,

,             

,            (6)

.

- на основании (5)  и  (6) имеем

,                             (7)

- последующие значения переменных найдем как

                                      (8)

Подставив в (5) и (6) вместо  соответственно  , по формулам (7)  найдем  , после чего  и т.д., пока не будет пройден  весь отрезок   и найдена искомая  функция  .

Порядок выполнения работы.   

- переписать  исходные данные к работе:

     - дифференциальное уравнение,

    ,   -  начальные  условия,

      - отрезок  для поиска решения,

- количество точек решения  ,

- точность расчетов принять равной  .

- ввести вспомогательную функцию , найти к ней начальное условие  ,

- преобразовать ДУ второго порядка в систему уравнений первого порядка,

- выполнить в  MathCad  расчеты  по формулам (5) – (8):

а) задать функцию ,

б) ввести начальные условия,

в) вычислить  ,

г) задать выражения из формул (5) и (6)  в следующем порядке:

,

д) по формулам (7) найти   и  ,

е) по формулам (8) найти последующие значения переменных,

ж) вычисления по пунктам б) – е) повторить еще 9 раз,

- составить итоговую таблицу расчетов ( колонки - ),

- выполнить проверку результатов:

а) сделать в Excel таблицу разностей по функции  ,

б) вычислить  по формуле численного дифференцирования,

в) подставить  в правую часть (1) величины   и найти  ,

г) сравнить   и   (теоретически они должны совпасть).

д) найти абсолютную и относительную погрешности метода,

- сделать выводы по работе.

Варианты исходных данных. Задано дифференциальное уравнение второго порядка,  начальные условия,  и интервал   для поиска его решения методом Рунге-Кутта.

Дифференциальное уравнение

Начальные

условия

Интервал

решения

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

24

25

  

         Пример расчета.   

  1.  Цель работы: проинтегрировать дифференциальное уравнение второго порядка методом Рунге-Кутта.
  2.  Исходные данные:

                                 n=10    h=(b-a)/10   

  1.   Преобразуем дифференциальное уравнение второго порядка в систему уравнений первого порядка.

          

  1.  Реализация в Mathcad , все расчеты выполним по следующей схеме:

                                      Продолжение                    Продолжение

                                        

 

  1.  Составляем таблицу:

I

X

dy

dz

yi

zi

0

0,000

0,013

0,251

0,000

0,000

1

0,105

0,039

0,250

0,013

0,251

2

0,210

0,066

0,249

0,052

0,502

3

0,315

0,066

0,247

0,118

0,751

4

0,420

0,091

0,244

0,184

0,749

5

0,525

0,091

0,241

0,275

0,995

6

0,630

0,117

0,237

0,366

0,990

7

0,735

0,116

0,232

0,483

1,232

8

0,840

0,141

0,227

0,624

1,464

9

0,945

0,165

0,221

0,765

1,459

10

1,050

0,164

0,214

0,930

1,685

  1.  Проверка:

x1=0.105

y1=0.013

y1’=z1=0.251

y’’1=z1                   найдем из последней колонки, создав таблицу разностей

                        dz,  d2z,  d3z,  d4z

dz

d2z

d3z

d4z

0

0,251

0

-0,002

-0,247

0,251

0,251

-0,002

-0,249

 

0,502

0,249

-0,251

 

 

0,751

-0,002

 

 

 

0,749

 

 

 

 

Для  определения точного значения, берем первые две разности  из

таблицы:

посчитаем погрешность при расчете методом Рунге-Кутта

7.Выводы по работе:  Выполняются студентом самостоятельно.

PAGE  5


 

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

137. Польско-Шведская интервенция 1609-1612 годов 44.22 KB
  В истории государства нашего были такие периоды, в которые под угрозу была поставлена его независимость и самобытность народа, если хотите. Так как объектом моего реферата является Польско-Шведская интервенция 1609-1612 гг., то для начала я определилась с одним из центральных понятий реферата – интервенция.
138. Теория и методология социологических исследований 1.54 MB
  Социология представляет собой профессиональный взгляд на мир человеческих отношений. Три типа социологического дискурса. Логика гипотетического рассуждения. Соотношение переменных и единиц исследования.
139. Разработка программного продукта для быстрого создания заметок и организации разнородных данных 1.41 MB
  Деление программ и технических средств на классы по типу предоставляемого функционала. Исследование возможных опасных и вредных факторов при эксплуатации ЭВМ и их влияния на пользователей.
140. Аналіз маркетингової діяльності ПП Діоніс 1.23 MB
  Теоретичні засади маркетингової діяльності торгівельного підприємства. Нормативно-правова база маркетингової діяльності в Україні. Удосконалення маркетингової діяльності ПП Діоніс. Організування розподілу продукції.
141. Феномен безработицы в России 38.65 KB
  В системе рыночных отношений важное место занимает рынок. Официальными безработными считаются трудоспособные граждане в трудоспособном возрасте (определяется законодательством). Исходя из необходимости учета безработных и принятия соответствующих государственных мер по обеспечению работой всех желающих.
142. Системный анализ московского метрополитена 273.5 KB
  Московский метрополитен, открытый 15 мая 1935 года, - это основа транспортной системы столицы. Он надежно связывает центр города с промышленными районами и жилыми массивами. На сегодняшний день доля Московского метрополитена в перевозке пассажиров
143. Історія української державності 40.71 KB
  Трипільська культура ,як зародки державотворчих процесів на сучасних Українських землях. Українська революція,або як ми прогавили незалежність. Суспільний нігілізм, тотальне розкрадання і не бажання влади змінити життя народу (а не своє) на краще.
144. Эскизный расчет радиолокационных систем 373.5 KB
  Расчет дальности действия РЛС при самоприкрытии цели активной шумовой помехой (АШП). Дальность действия и зона обнаружения в условиях применения маскирующих пассивных помех. Выбор и обоснование передающего устройства.
145. Аппараты с перемешивающими устройствами 240 KB
  В аппарате применяется лопастная мешалка. Основная задача переме-шивания – равномерное распределение вещества или температуры в перемешивающем объёме. Привод выбирается по ОСТ 26-01-1205-75 в соответствии с частотой вращения мешалки, номинальным давлением в корпусе аппарата.