Читать курсовая по математике: "Методы решения дифференциальных уравнений" Страница 1

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

КУРСОВАЯ РАБОТА

Методы решения дифференциальных уравнений Введение Дифференциа́льное уравне́ние -уравнение, связывающее значение некоторой неизвестной функции в некоторой точке и значение её производных различных порядков в той же точке. Дифференциальное уравнение содержит в своей записи неизвестную функцию, ее производные и независимые переменные. Решить дифференциальное уравнение не просто. Есть большое число способов их решения.1.Теоретическая часть .1 Модифицированный метод ЭйлераСуть этого метода состоит в следующем. Сначала вычисляются вспомогательные значения искомой функции y в точках t = ti +с помощью формулы:

y = yi +fi = yi +f(ti, yi).

Затем находится значение правой части уравнения в средней точке

= f(t, y) и затем полагается yi+1 = yi + h f, i = 0, 1, …, n - 1.

.2 Формула Ньютона

Интерполяционный многочлен легко определяется если его построить в виде:

(x) = С0 + С1(x - x0) + C2(x - x0) (x - x1) + ...+ Cn(x - x0)(x - x1) ... (x - xn-1)

Исходя из условия интерполяции , для коэффициентов Ci получим систему уравнений треугольного вида

(x0) = С0(x1) = С0 + С1(x1 - x0 )(x2) = С0 + С1(x1 - x0 ) + C2(x2 - x0 )(x2 - x1 )

. . . . . . . . . . . . . . . . . . . . . . . . . .. . . . . . . . . . . . . . . . . . . . . . . . . .(xn ) = С0 + С1(xn - x0) + C2(xn - x0)(xn - x1) + ...+ Cn(xn - x0)(xn - x1) ... (xn- xn-1)

Из этой системы легко находятся

и так далее.

.3 Дихотомия

Пусть задана функция .

Разобьём мысленно заданный отрезок пополам и возьмём две симметричные относительно центра точкиитак, что:

,

где- некоторое число в интервале

Отбросим тот из концов изначального интервала, к которому ближе оказалась одна из двух вновь поставленных точек с максимальным значением (напомним, мы ищем минимум ), то есть:

Если , то берётся отрезок , а отрезокотбрасывается.

Иначе берётся зеркальный относительно середины отрезок , а отбрасывается .

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

На каждой итерации приходится вычислять новые точки. Можно добиться того, чтобы на очередной итерации было необходимо высчитывать лишь одну новую точку, что заметно способствовало бы оптимизации процедуры. Это достигается путём зеркального деления отрезка в золотом сечении, в этом смысле метод золотого сечения можно рассматривать, как улучшение метода дихотомии с параметром , где- золотое сечение .

.4 Метод трапеций

Выведем формулу трапеций из геометрических соображений. Заменим график функции y = f(x) ломаной линией (рис.1.1), полученной следующим образом. Из точек a = x0, x1, x2,…, xn = b проведем ординаты до пересечения с кривой y = f(x). Концы ординат соединим прямолинейными отрезками.

Рис. 1.1

Тогда площадь криволинейной трапеции приближенно можно считать равной площади фигуры, составленной из трапеций. Так как площадь трапеции, построенной на отрезке [xi, xi+1] длины h =, равна h , то, пользуясь этой формулой для i = 0, 2, … , n - 1, получим квадратурную формулу трапеций:

=»Iтр =h=

2. Расчетная часть

kursovaya;q:real;zagacha1(q:real);x, y, yp, h, b: real;: integer;f (xn, yn: real): real;:= 2 * (xn * xn + yn);;:= 0;:= 0;:= 1;:= 0.1;:= 1;('Модифицированный метод Эйлера');('');('dy/dx = 2*(x^2+y), y(0) = 1');('');

writeln('| i | x | yp | y |');

writeln('-----------------------------');('|', i:2, ' |', x:4:1, ' |', yp:7:4, ' |', y:7:4, ' |');:= y + h/2 * f(x, y);:= y + h * (f(x+h/2, yp));:= i + 1;:= x + h;x > b;;;zagacha2(q:real);i, k0, k1,

Похожие работы