Курсовая работа: Интеграционный метод Эйлера для решения линейных систем алгебраических уравнений

· жесткость системы затрудняет анализ результата решения

ЗАКЛЮЧЕНИЕ

В данной работе был исследован явный метод Эйлера для решения обычных и жестких систем ОДУ. Было проанализировано влияние величины шага интегрирования на ошибку аппроксимации, и ее влиянии на число итераций. Для этого была написана программа (Приложение 1), реализующая метод, и протестирована при различных исходных данных.

СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ

1 Ортега Дж., Рейнболдт В. Итерационные методы решения нелинейных систем уравнений со многими неизвестными.-М.: Мир, 1975.- 558 стр.

2 Самарский А.А., Гулин А.В. Численные методы: Учеб. Пособие для вузов.- М.: Наука,1989.- 432 стр.

3 Сарычева О.М. Численные методы в экономике / О.М. Сарычева.- Новосибирск, 1995.- 67 стр.

ПРИЛОЖЕНИЯ

ПРИЛОЖЕНИЕ 1

Текст головной программы:

h=0.01; % шаг интегрирования

t0=0; % начальное время интегрирования

x0=[0;0];

Edop=0.01;

[t_out,y_out]=RK1(t0,x0,h,Edop); % вызов RK1

ytoch=FunT(t_out); % точное решение

% построение графика решения методом Рунге-Кутта 1

plot(t_out,y_out);

grid;

title('Solution for x1 and x2 by method Runge-Kutta 1');

ylabel('x');

xlabel('t');

Текст программы для решения ОДУ методом Эйлера с постоянным шагом:

function [t_out,y_out]=RungeKutta1(t0,x0,h,Edop);

% функция решения методом Рунге-Кутта 1

t=t0;

x=x0;

t_out=t;

y_out=x0;