Введение 3
1 Численное решение задачи для математической модели системы 6
гарпактициды-диатомовы водоросли
1.1 Постановка задачи 6
1.2 Алгоритм решения задачи 7
1.3 Результаты вычислительных экспериментов 7
2 Численное решение обратной задачи 14
2.1 Постановка задачи 14
2.2 Алгоритм решения задачи 15
2.3 Тестовые расчёты для обратной задачи 16
Заключение 28
Список использованных источников 29
Приложения 30
Рассматривается система трофического сообщества гарпактициды- диатомовые водоросли. Это система основана на процессе трофотаксиса, т.е. перемещении микроорганизмов в направлении питательных веществ. В работе Тютюнова Ю.В.[1] была выведена математическая модель для системы такого рода. Полученная математическая модель является системой типа хищник жертва и в области Q=[0,T]x[0,1] имеет следующий вид:
N(0, x) = N0 (x), 0 < x < 1,
R(0, x) = R (x), 0 < x < 1,
v(0, x) = v0 (x), 0 < x < 1,
dN (t, x)
dx
v(t, x)| x=c., = °,
В этой система коэффициент k - коэффициент трофотаксиса. Этот коэффициент отображает поисковую активность хищника, другими словами, то на сколько целенаправленно хищник ищет свою жертву. При k=0 поиск жертв происходит случайным образом. Коэффициент a - коэффициент поисковой эффективности хищника, характеризующий интенсивность атак хищника. Аналитически в работе [1] было получено что при k>0, модель устойчива, пока коэффициент трофотаксиса не достигнет некоторого критического значения k*, после которого модель теряет устойчивость.
Условия определения этого критического значения не было получено. Известно только то, что при различных входных параметрах критическое значение также различно. Так же важно отметить, что на практике редко случается, когда коэффициент трофотаксиса равен нулевому значению.
Так как не всегда есть возможность экспериментально провести необходимые замеры значения параметров математической модели таких, как начальные условия, граничные условия, коэффициент смертности, коэффициенты диффузии, коэффициент трофотаксиса, коэффициент поисковой эффективности хищника, возникает проблема нахождения этих коэффициентов. Существуют различные методы определения значения различных параметров математической модели, которые по каким-либо причинам могут быть неизвестны. Один из методов определения значения этих параметров, это решение обратной задачи. Обратные задачи - это такой тип задач, в случае которых некоторые значения параметров задач могут быть неизвестны.
Первые работы по исследованию обратных задач появились в первой половине 20-ого века. С появлением мощных вычислительных ресурсов область приложений обратных задач стала охватывать всё большее и большее количество научных дисциплин. И проблема разработки методов численного решения обратных задач становится всё более актуальной.
При решении прямых задач математической физики, исследователи стремятся найти функции, которые описывают различные физические явления, например, распространения различных примесей в жидкостях, распространения тепла в различных веществах и так далее. В ходе решения подразумевается, что требующиеся входные данные для математической модели являются известными. На практике это редко достижимо и различные значения параметров модели неизвестны. В следствии чего возникает необходимость ставить и решать обратные задачи. В большинстве случаем эти задачи некорректны[2][3][4][5], что означает нарушение хотя бы одного из свойств корректности (условие существования, условие единственности, условие устойчивости решения по отношению к возмущениям во входных данных).
В проделанной работы численно решается задача для математической модели системы гарпактицид-диатомовых водорослей. Так же осуществляется постановка обратной задачи для рассматриваемой математической модели и её численное решение. И проводятся вычислительные эксперименты по исследованию поведения решения относительно возмущений во входных данных, а именно в условии переопределения.
Цель работы- численное решение прямой и обратной задачи для математической модели системы гарпактицид-диатомовых водорослей. Для выполнения поставленной цели необходимо выполнить следующие задачи:
- численное решение прямой и обратной задачи.
- разработка алгоритма решения обратной задачи.
- разработка программы для реализации численных методов при помощи пакета прикладных программ Matlab[6].
-проведение вычислительных экспериментов.
В работе получены следующие результаты:
1. Численно решена прямая и обратная задачи для математической модели системы гарпактицид-диатомовых водорослей.
2. Был разработан алгоритм решения обратной задачи.
3. Для реализации численных методов разработана программа при помощи использования пакета прикладных программ Matlab[6].
4. Проведены вычислительные эксперименты.
1. Тютюнов, Ю.В. Построение, исследование и приложения математических моделей пространственно временной динамики популяционных систем: автореф. дис. ... док. физ.-мат. наук : 03.00.02 / Тютюнов Юрий Викторович.
- Красноярск, 2009. -50 с.
2. Кучунова, Е.В. Численная идентификация коэффициентов параболических уравнений / Е.В. Кучунова, В.Е. Распопов - Вестник КрасГУ. Серия Физ.- мат. науки, 2004г. №5/2 с. 7-14.
3. Мандрик, Ю.В. Численная идентификация коэффициентов одномерного параболического уравнения/ Ю.В.Мандрик, В.Е. Распопов- Вестник КрасГУ, 2006г., №1.
4. Жак, Т.Ю. Численная идентификация свободного члена специального вида в параболическом уравнении/ Т.Ю. Жак, В.Е. Распопов// Международная конференция Алгебра и её приложения: Тезисы докладов. Красноярск, 200г.
- 186 с.
5. Самарский, А. А. Численные методы решения обратных задач математической физики/ А.А. Самарский, П.Н.Вабищевич-Москва, 2009г.
6. Дьяконов, В. П. MATLAB. Полный самоучитель: учеб. пособие / В.П. Дьяконов. - Москва: ДМК Пресс, 2012.-768 с.
7. Бахвалов, Н.С. Численные методы: учебное пособие / Н.С. Бахвалов, Н.П. Жидков, Г.М. Кобельков; Издание 7-е, Бином. Лаборатория знаний, 2012. - 636с.
8. Самарский, А.А. Численные методы: учебное пособие для студентов вузов, обучающихся по специальности «Прикладная математика»/ А.А. Самарский, А.В. Гулин; Москва: Наука, 1989. - 432 с.
9. Зализняк, В.Е. Теория и практика по вычислительной математике: учебное
пособие/ В.Е. Зализняк, Г.И. Щепановская.- Красноярск: Сибирский федеральный университет, 2012. - 174 с.