Рабочие листы
к вашим урокам
Скачать
1 слайд
Моделирование реальных процессов в пакете расширений
«dynamics» системы Maxima
2 слайд
Однако в ряде случаев отыскать символьное решение ОДУ в достаточно компактном виде невозможно. В этом случае целесообразно использовать численные методы. Maxima включает пакет расширения dynamics, позволяющий проинтегрировать систему ОДУ методом Рунге-Кутта.
3 слайд
Пример:
Решить ОДУ
при .
Используем пакет dynamics.
(%i1) load(''dynamics'')$
Выбираем шаг интегрирования 0,02.
(%i2) sol:rk([4*x^2-4*y^2,y^2-x^2+1],[x,y],
[-1.25,0.75],[t,0,4,0.02]);
В результате решения получаем список значений в формате [[ ]].
(%i1) load("dynamics")$
(%i2) rp1:4*x^2-4*y^2;
(%i3) rp2:y^2-x^2+1;
4 слайд
Изменение концентраций при моделировании автоколебательной химической реакции (брюсселятора)
Пример
Описание модели брюсселятора в Maxima приведено в следующем командном файле:
load("dynamics");
load("draw");
B:0.5;
eq1:-(B+1)*y0+y0^2*y1+1;
eq2:B*y0-y0^2+1;
t_range:[t,0,10,0.1];
sol: rk([eq1,eq2],[y0,y1],[1,1],t_range)$
len:length(sol);
t:makelist(sol[k][1],k,1,len)$
y0:makelist(sol[k][2],k,1,len)$
y1:makelist(sol[k][3],k,1,len)$
draw2d(title="Brusselator",xlabel="t",ylabel="y0,y1",
grid=true,points_joined = true,
points(t,y0),points(t,y1),terminal=eps);
5 слайд
система Лотка-Вольтерра
6 слайд
Пример формирования динамического хаоса (аттрактор Лоренца)
Рабочие листы
к вашим урокам
Скачать
6 672 333 материала в базе
Настоящий материал опубликован пользователем Тоболев Геннадий Юрьевич. Инфоурок является информационным посредником и предоставляет пользователям возможность размещать на сайте методические материалы. Всю ответственность за опубликованные материалы, содержащиеся в них сведения, а также за соблюдение авторских прав несут пользователи, загрузившие материал на сайт
Если Вы считаете, что материал нарушает авторские права либо по каким-то другим причинам должен быть удален с сайта, Вы можете оставить жалобу на материал.
Удалить материал
Ваша скидка на курсы
40%
Курс профессиональной переподготовки
500/1000 ч.
Курс профессиональной переподготовки
600 ч.
Курс профессиональной переподготовки
300/600 ч.
Курс профессиональной переподготовки
300/600 ч.
Мини-курс
4 ч.
Оставьте свой комментарий
Авторизуйтесь, чтобы задавать вопросы.