2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки




 
 Использование функции ode45 в Matlab
Сообщение22.04.2007, 20:56 
Помогите пожалуйста!!! У меня имеется 3 уравнения, их необходимо решить с помощью функции ode45 в Matlab...
Вот и сама система уравнений:
m1*(d^2 x1/dt^2)= -k1x1-k2(x1-x2)
m2*(d^2 x2/dt^2)= -k2(x2-x1)-k3(x2-x3)
m3*(d^2 x3/dt^2)= -k4x3-k3(x3-x2)
Необходимо найти 3 значения: x1, x2 и x3.
Кто знает, как пользоваться этой функцией, подскажите!!!

 
 
 
 
Сообщение24.04.2007, 16:56 
Аватара пользователя
Читаем help, там естъ разобранный пример для:

$y_1'=y_2$
$y_2'=\mu(1-y_1^2)y_2-y_1$

Цитата:
The code below represents the van der Pol system in the function, vdp1. The vdp1 function assumes that $\mu=1$. The variables $y_1$ and $y_2$ are the entries $y(1)$ and $y(2)$ of a two-element vector.

Код:
function dydt = vdp1(t,y)
dydt = [y(2); (1-y(1)^2)*y(2)-y(1)];


For the van der Pol system, you can use ode45 on time interval [0 20] with initial values $y(1) = 2$ and $y(2) = 0$.
Код:
[t,y] = ode45(@vdp1,[0 20],[2; 0]);


а дальше делаете просто по аналогии для своей системы

 
 
 [ Сообщений: 2 ] 


Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group