Анкилостомоз собак и кошек

Вид работы:

Доклад
Предмет:

Биология
Язык:

Русский
,
Формат файла:
MS Word

4,64 kb
Опубликовано:

2009-01-12

Все доклады по биологии

Скачать доклад Читать текст online Посмотреть все доклады

Вы можете узнать стоимость помощи в написании студенческой работы.

Анкилостомоз собак и кошек

Кафедра высшей математики

Отчет по лабораторной работе № 4

Тема: Решение дифференциальных уравнений в системе MathCAD.

Выполнил студент ИГЭУ

Теплоэнергетического Факультета

Кафедры тепловых электростанций

Группы I-2x

Атаманчук М. С.

Проверил: Крутов А.О.

Лабораторная работа № 4.

Тема: Решение дифференциальных уравнений

Цель: Научиться решать дифференциальные уравнения в системе MathCAD методами Рунге - Кутты (rkfixed, rkadapt), Булирша - Штера (Bulstoer) и Odesolve

Описание метода Рунге - Кутты

Метод Рунге-Кутта четвертого порядка для решения уравнения первого порядка.

Идея Рунге-Кута состоит в том, чтобы использовать метод неопределённых коэффициентов. Наиболее употребительным методом Рунге-Кутта решения уравнения первого порядка y' = F(x,y) (2.1.1) является метод четвертого порядка, в котором вычисления производятся по формуле:

yk+1 = yk +(k1 +2k2 +2k3 +k4 )/6, (2.1.2)

где = Fk h = F(xk , yk )h = F(xk +h/2, yk +k1 /2)h = F(xk +h/2, yk +k2 /2)h = F(xk +h, yk +k3 )h, = 0, ..., n-1

h = (xf -x0 )/n (2.1.3)

Описание метода Булирша - Штера (Bulstoer)

Метод Булирша-Штера - это метод решения системы обыкновенных дифференциальных уравнений первого порядка с гладкими правыми частями. Гладкость правых частей является необходимой для работы метода. Если правые части вашей системы не являются гладкими или содержат разрывы, то лучше использовать метод Рунге-Кутта. В случае же гладкой системы метод Булирша-Штера позволяет добиться существенно большей точности, чем метод Рунге-Кутта.

Принцип работы метода

Основной идеей метода является вычисление состояния системы в точке x+h, как результата двух шагов длины h/2, четырех шагов длины h/4, восьми шагов длины h/8 и так далее с последующей экстраполяцией результатов. Метод строит рациональную интерполирующую функцию, которая в точке h/2 проходит через состояние системы после двух таких шагов, в точке h/4 проходит через состояние системы после четырех таких шагов, и т.д., а затем вычисляет значение этой функции в точке h = 0, проводя экстраполяцию.

Гладкость правых частей приводит к тому, что вычисленное при помощи экстраполяции состояние системы оказывается очень близко к действительному, а использование рациональной экстраполяции вместо полиномиальной позволяет ещё больше повысить точность.

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

Задача №1

Найдите два решения дифференциального уравнения на отрезке [1; 5] с начальными условиями соответственно , . Проконтролируйте достижение точности 0.001. Нарисуйте их графики (на одном чертеже). Определите значения этих решений в точке x = 4.85.

Метод Рунге - Кутты (rkfixed)

Создадим функцию

Далее применим к ней функцию rkfixed

Для достижения точности 0.001 ниже полученных матриц Z, Z1 выполним оператор TOL :=0.001.

После этого снова найдем решение, но обозначим его иначе:

Возьмем два вектора и , являющиеся столбцами значений первого и второго решений. Найдем модуль разности этих векторов.

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

Построим графики этих функций.

Для нахождения значения этих решений в точке x = 4.85 применим к функции функцию rkfixed, но правое ограничение выстави в виде x = 4.85 и построим матрицу решений.

Решения дифференциального уравнения на отрезке [1; 5] с начальными условиями, в точке x = 4.85 будет y(4.85)=1.239 и y(4.85)=1.233 соответственно.

Метод Рунге - Кутты (rkadapt)

Создадим функцию

Далее применим к ней функцию rkadapt

Для достижения точности 0.001 ниже полученных матриц Z, Z1 выполним оператор TOL :=0.001.

После этого снова найдем решение, но обозначим его иначе:

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

Построим графики этих функций.

Для нахождения значения этих решений в точке x = 4.85 применим к функции функцию rkadapt, но правое ограничение выставим в виде x = 4.85 и построим матрицу решений.

Решения дифференциального уравнения на отрезке [1; 5] с начальными условиями, в точке x = 4.85 будет y(4.85)=-0.127 и y(4.85)=-0.142 соответственно.

Метод Булирша - Штера (Bulstoer)

Создадим функцию

Далее применим к ней функцию Bulstoer

Для достижения точности 0.001 ниже полученных матриц Z, Z1 выполним оператор TOL :=0.001.

После этого снова найдем решение, но обозначим его иначе:

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

Построим графики этих функций.

Для нахождения значения этих решений в точке x = 4.85 применим к функции функцию Bulstoer но правое

ограничение выставим в виде x = 4.85 и построим матрицу решений.

Решения дифференциального уравнения на отрезке [1; 5] с начальными условиями, в точке x = 4.85 будет y(4.85)=-0.132 и y(4.85)=-0.146 соответственно.

Odesolve

Набираем в MathCad

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

Решения дифференциального уравнения на отрезке [1; 5] с начальными условиями, в точке x = 4.85 будет y(4.85)=-0.129 и y(4.85)=-0.143 соответственно.

Общий график решения (Рунге - Кутты (rkfixed, rkadapt), Булирша

Штера (Bulstoer) и Odesolve)

Таблица результатов с разными методами

	rkfixed	Rkadapt	Bulstoer	Odesolve
y(4.85)=1.239y(4.85)=-0.127y(4.85)=-0.132y(4.85)=-0.129
y(4.85)=1.233y(4.85)=-0.142y(4.85)=-0.146y(4.85)=-0.143

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

Z - rkfixed, R - Rkadapt, B - Bulstoer

Задача №2

Найдите решение дифференциального уравнения на отрезке [2; 4] с начальным условием . Нарисуйте график этого решения. Найдите значение этого решения в точке x = 3.75 с точностью .

Метод Рунге - Кутты (rkfixed)

Создадим функцию

Далее применим к ней функцию rkfixed.

Для достижения точности 10^-5 ниже полученной матрицы Z выполним оператор TOL := 10^-5.

После этого снова найдем решение, но обозначим его иначе:

Возьмем два вектора Z^<1> и Z1^<1>, являющиеся столбцами значений первого и второго решений. Найдем модуль разности этих векторов.

Видим что результаты совпадают в пределах заданной погрешности что и было необходимо для нашей задачи.

Построим графики этих функций.

Для нахождения значения этих решений в точке x = 3.75 применим к функции функцию rkfixed но правое ограничение выставим в виде x = 3.75 и построим матрицу решений.

Решением дифференциального уравнения на отрезке [2; 4] с начальными условиями, в точке x = 3.75 будет y(3.75)=0.45

Метод Рунге - Кутты (rkadapt)

Создадим функцию

Далее применим к ней функцию rkadapt

Для достижения точности 10^-5 ниже полученной матрицы Z выполним оператор TOL := 10^-5.

После этого снова найдем решение, но обозначим его иначе:

Видим что результаты совпадают в пределах заданной погрешности что и было необходимо для нашей задачи.

Построим графики этих функций.

Решением дифференциального уравнения на отрезке [2; 4] с начальными условиями, в точке x = 3.75 будет y(3.75)=0.45

Метод Булирша - Штера (Bulstoer)

Создадим функцию

Далее применим к ней функцию Bulstoer

Для достижения точности 10^-5 ниже полученной матрицы Z выполним оператор TOL := 10^-5.

После этого снова найдем решение, но обозначим его иначе

Видим что результаты расходятся, для построения графика и нахождения решения будем использовать матрицу Z1, так как она более точная.

Построим графики этой функции.

Для нахождения значения этих решений в точке x = 3.75 применим к функции функцию Bulstoer но правое ограничение выставим в виде x = 3.75 и построим матрицу решений.

Решением дифференциального уравнения на отрезке [2; 4] с начальными условиями , в точке x = 3.75 будет y(3.75)=0.45

Odesolve

Набираем в MathCad

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

Общий график решения (Рунге - Кутты (rkfixed, rkadapt), Булирша - Штера (Bulstoer) и Odesolve)

Таблица результатов с разными методами

	Rkfixed	Rkadapt	Bulstoer	Odesolve
	y(3.75)=0.45	y(3.75)=0.45	y(3.75)=0.45	y(3.75)=0.45

Результаты различных методов одинаковы но по общему графику видно что наиболее гладкая функция у метода Odesolve. Z - rkfixed, R - Rkadapt, B - Bulstoer

Задача №3

Найдите решение системы дифференциальных уравнений на отрезке [0; 4] с начальными условиями , . Нарисуйте графики компонент и фазовый портрет решения.