Читать практическое задание по математике: "Решение задач методами Эйлера и Рунге-Кутта" Страница 1


назад (Назад)скачать (Cкачать работу)

Функция "чтения" служит для ознакомления с работой. Разметка, таблицы и картинки документа могут отображаться неверно или не в полном объёме!

1. Построить кубический сплайн, интерполирующий функцию у = (х) на [1,00; 1,20] для равномерного разбиения с шагом h = 0,04: (х) = ln x Найти значения в точках 1,05; 1,13; 1,17.

Решение

Построим таблицу значений функции на интервале [1,00; 1,20] с шагом

h = 0,04:

x

(х) = lnx

1

0

1,04

0,039221

1,08

0,076961

1,12

0,113329

1,16

0,14842

1,2

0,182322

Сплайн-интерполяция таблично заданной функции

    На отрезке [a, b] задать одномерную сетку

x = {xi / xi = xi –1 + hi, hi > 0, i = 1, 2, 3, …, n; x0 = a, xn = b} и значения yi = f(xi) в узлах сетки xi, i = 0, 1, 2, …, n.

Задать x*  (a, b).

    Положить ai = yj, i = 0, 1, 2, …, n. Составить и решить трех диагональную систему методом прогонки:

Определить значения коэффициентов ci, i = 0, 1, 2, …, n.

    Определить значения коэффициентов di и bi, i = 1, 2, 3, …, n, воспользовавшись формулами:

di = (ci – ci – 1) / hi, i = 1, 2, …

    Определить значение индекса 0 < k  n из условия x*  [xk – 1, xk]. Вычислить по формуле

S(x*) = Sk(x*) = ak + bk(x* – xk) + (ck / 2)(x* – xk)2 + (dk / 6)(x* – xk)3.

    Процесс завершен: S(x*) – результат интерполяции табличных данных в точку x*  (a, b).

Результаты вычислений удобнее представлять в виде таблицы:

ai

bi

ci

di

0,03922

0,96467

-1,188280

-29,70700

0,07696

0,92494

-0,798322

9,74897

0,11333

0,89366

-0,765997

0,80813

0,14842

0,85986

-0,92391

-3,94780

0,18232

0,84138

0,00000

23,09770

Значение функции в точке находится по формуле: S(x*) = Sk(x*) = ak + bk(x* – xk) + (ck / 2)(x* – xk)2 + (dk / 6)(x* – xk)3

2. Найти решение задачи Коши для дифференциального уравнения на равномерной сетке [a, b] с шагом 0,2 методом Эйлера и классическим методом Рунге-Кутта , , 0  х  1 Решение. Метод Эйлера - разностная аппроксимация Эйлера. Точность метода . Метод Рунге-Кутта

дифференциальный интерполирующий уравнение сплайн

Результаты вычислений удобнее представлять в виде таблиц:

Метод Эйлера

x

y

0

0

1

0,2

0,2

1

0,4

0,416

1.04

0,6

0,67392

1.1232

0,8

1,00639

1.25798

1

1,45926

1.45926

Метод Рунге-Кутта

i

=

0

0

1

0

0,02

0,0202

0,040808

1,0202

1

0,2

1,0202

0,0408081

0,0624363

0,0630852

0,0866629

1,08329

2

0,4

1,08329

0,086663

0,112662

0,113962

0,14367

1,19722

3

0,6

1,19722

0,143666

0,177667

0,180047

0,220362

1,37713

4

0,8

1,37713

0,22034

0,267713

0,271977

0,329821

1,64872

5

1

1,64872

0,329743

0,398989

0,406607

0,493278

2,05442



Интересная статья: Быстрое написание курсовой работы