14641

Решить заданное дифференциальное уравнение методом Эйлера с применением «ручных» вычислений

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

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

Решить заданное дифференциальное уравнение методом Эйлера с применением ручных вычислений а также с помощью программы составленной на языке программирования Паскаль с шагом 2h и с шагом h. Свести результаты вычисления в одну таблицу и сопоставить точность полученных...

Русский

2013-06-08

123.98 KB

72 чел.

Решить заданное дифференциальное уравнение методом Эйлера с применением «ручных» вычислений, а также с помощью программы, составленной на языке программирования Паскаль с шагом 2h  и с шагом h. Свести результаты вычисления в одну таблицу и сопоставить точность полученных значений функции. Пользуясь таблицей, сделать прикидку графика интегральной кривой (ломаная Эйлера). С помощью прикладного программного средства (MathCAD) методом Эйлера обеспечить вывод полученных решений в виде таблиц и графиков.

1) Решение «вручную» с помощью MS Excel:

y=f(x,y)

a

b

y0

h

x+sin(y/3)

1.6

2.6

4.6

0.1

1.1) При h=0.1:

1.2) При  h=0.2:

i

Xi1

Yi1

h

i

Xi2

Yi2

h

 0

1,6

4,6

0,1

0

1,6

4,6

0,2

1

1,7

4,859930

1

1,8

5,119860

2

1,8

5,129809

2

2,0

5,678018

3

1,9

5,408842

3

2,2

6,267746

4

2,0

5,696160

4

2,4

6,881464

5

2,1

5,990831

5

2,6

7,511425

6

2,2

6,291888

6

-

8,150509

7

2,3

6,598345

8

2,4

6,909227

9

2,5

7,223592

10

2,6

7,540557

11

-

7,859315

Относительные погрешности метода:

Влияние шага на точность вычислений:

2) Решение в MathCad:

построим график:

Относительные погрешности метода:

3) Решение в Matlab:

3.1) При  h=0.1:

>> a=1.6;

>>  b=2.6;

>> h=0.1;

>> y=4.6;

>> for i=0:1:10

x=a+h*i;

f=x+sin(y/3);

y=y+h*f;

end

>> y

y =

   7.8593

3.2) При  h=0.2:

>> a=1.6;

>> b=2.6;

>> h=0.2;

>> y=4.6;

>> for i=0:1:5

x=a+h*i;

f=x+sin(y/3);

y=y+h*f;

end

>> y

y =

   8.1505

Относительные погрешности метода:

Построим график:

создадим документ Blank M-file, напечатаем

function f=fun2(t,y)

f=t+sin(y/3);

сохраним под именем fun2

>> [T,Y]=ode45('fun2',[1.6:0.1:2.6],[4.6],[])

T =

   1.6000

   1.7000

   1.8000

   1.9000

   2.0000

   2.1000

   2.2000

   2.3000

   2.4000

   2.5000

   2.6000

Y =

   4.6000

   4.8650

   5.1395

   5.4226

   5.7136

   6.0113

   6.3148

   6.6231

   6.9353

   7.2504

   7.5676

>> plot(T,Y)

>> [T1,Y1]=ode45('fun2',[1.6:0.2:2.6],[4.6],[])

T1 =

   1.6000

   1.8000

   2.0000

   2.2000

   2.4000

   2.6000

Y1 =

   4.6000

   5.1395

   5.7136

   6.3148

   6.9353

   7.5676

>> hold on

>> plot(T1,Y1)

Сводная таблица результатов:

Метод

Результат

δ(y), %

h

2h

h

2h

«вручную» с Excel

7,859315

8,150509

MathCad

7,859315

8,150509

Matlab

7.8593

8.1505

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

МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РТ

АЛЬМЕТЬЕВСКИЙ ГОСУДАРСТВЕННЫЙ

НЕФТЯНОЙ ИНСТИТУТ

Кафедра информатики

ЛАБОРАТОРНАЯ  РАБОТА №10

По дисциплине: «Прикладное программирование»

На тему: «Численное решение обыкновенного дифференциального

уравнения первого порядка методом Эйлера»

Вариант №42

Выполнил: студент группы 10-21   Хохлов Р.С. Проверил: доцент каф. информатики Амиров Д.Ф.

Альметьевск 2013


 

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

3552. Интерпретатор введенного кадра или УП 67 KB
  Программа интерпретатор Блок – схема алгоритма Kadr: ClearScreen белый фон PrintXY 0,0,WW Locate 0,2 CursorOn ввод кадра mov di, OFFSET BuferData mov...
3553. Реконструкция слесарно-механического участка в условиях ПЧ-4 1.1 MB
  Автомобили работают 365 дней в году, в две смены, 8 часов в смену. Климат умеренно холодный, рельеф слабохолмистый, дорожное покрытие преимущественно асфальтобетон. КУЭ – III. Осуществление такого рода услуг возможно благодаря универсальному подвижному составу и хорошей производственной базе
3554. Аерологія виробок 1.92 MB
  Модуль 1. Шахтне повітря, кліматичні умови, пилогазовий режим та його контроль. Вступ. Призначення та задачі шахтної вентиляції по створенню безпечних умов праці шахтарів. Мета та задачі дисципліни. Тема 1.1 Шахтне повітря. Поняття «шахтне повітря»....
3555. Сутність цивільного та земельного законодавств, їх компетенція у вирішенні питання набуття права приватної власності на земельні ділянки громадянами України 315.5 KB
  Актуальність даної теми обґрунтовується постійно зростаючою зацікавленістю встановлення режиму власності на землю, з метою здійснення господарської діяльності фізичними і юридичними особами на земельних ділянках. З утворенням України як незале...
3556. Історія економічних вчень 144.69 KB
  Предмет і завдання курсу Історія економічних вчень Економічне життя суспільства вивчається системою економічних наук - це науки про загальні закони экономічного розвитку, галузеві економічні науки: науки, шо відмічають конкретні процеси і яв...
3557. Поняття та особливості сільськогосподарського виробництва 36.41 KB
  Сільське господарство розвивається на основі різних форм власності і видів господарювання. Рівень господарювання і характер економічної відокремленості цих господарств визначають специфічні особливості їхніх взаємовідносин з державою і певні відмінності у способах використання механізму дії економічних законів.
3558. Адміністративне право 121.5 KB
  Вступ до адміністративного права. Адміністративне право — це одна з профільних, фундаментальних галузей правової системи України. Адміністративне право визначається як сукупність юридичних норм та правових інститутів, призначених для ре...
3559. Ремонт машин і обладнання підприємств виробництва будівельних матеріалів 304 KB
  Метою курсового проекту по ремонту механічного устаткування є привити студентам навички рішення інженерних питань ремонту, економічних обґрунтувань і планування різних видів ремонту та складання необхідної документації.
3560. Изучение внешнего фотоэффекта 67.5 KB
  Цель работы: изучение внешнего фотоэффекта. Задача: определение световой и вольт-амперной характеристики фотоэлемента. Техника безопасности: напряжение 220 В подается от сети на трансформатор и выпрямитель, поэтому соответствующие токоведущие...