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

  • Вид работы:
    Доклад
  • Предмет:
    Биология
  • Язык:
    Русский
    ,
    Формат файла:
    MS Word
    4,64 kb
  • Опубликовано:
    2009-01-12
Вы можете узнать стоимость помощи в написании студенческой работы.
Помощь в написании работы, которую точно примут!

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

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

Отчет по лабораторной работе № 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.

 


    

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



















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





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

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











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















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

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

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



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





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

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


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



Видим что результаты расходятся, для построения графика и нахождения решения будем использовать матрицу 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] с начальными условиями , . Нарисуйте графики компонент и фазовый портрет решения.

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

Сформируем вектор начальных условий


Сформируем вектор правых частей

Найдем матрицу значений решения этой системы











Построим фазовый портрет решения графики компонент решения.


Фазовый портрет решения                        Графики компонент решения

Метод Рунге - Кутты Rkadapt

Сформируем вектор начальных условий


сформируем вектор правых частей

Найдем матрицу значений решения этой системы


Построим фазовый портрет решения графики компонент решения

Фазовый портрет решения              Графики компонент решения


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

Сформируем вектор начальных условий


Сформируем вектор правых частей

Найдем матрицу значений решения этой системы













Построим фазовый портрет решения графики компонент решения

Фазовый портрет решения                       Графики компонент решения



Odesolve

Введем в MathCAD










 


Построим фазовый портрет решения графики компонент решения

Фазовый портрет решения                  Графики компонент решения


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

Общий фазовый портрет решения


Общие графики компонент решения

Наиболее точным является метод Odesolve, так как из общих графиков видно что его функция более гладкая.

Z - rkfixed, R - Rkadapt, B - Bulstoer

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

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


Сформируем вектор начальных условий



Сформируем вектор правых частей

Найдем матрицу значений решения этой системы



















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


Метод Рунге - Кутты Rkadapt

Сформируем вектор

начальных условий


Сформируем вектор

правых частей

Найдем матрицу значений решения этой системы

дифференциальный уравнение MathCAD



















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


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

Сформируем вектор начальных условий


Сформируем вектор

правых частей

Найдем матрицу значений решения этой системы


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

Odesolve





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

Общий график

На графике видно что все методы совпадают по точности.Y - rkfixed, R - Rkadapt, B - Bulstoer

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

Похожие работы на - Анкилостомоз собак и кошек

 

Не нашел материал для своей работы?
Поможем написать качественную работу
Без плагиата!