Sitio web de programación diferencial

Definición de parámetros

m 1 = 1; m2 = 2; m3 = 3; 6;k4 = 7;

P0 = 8; w = 9;

Ecuación diferencial

a =[k 1 k2-k 1 0;-k2 k2 k3-k3; 0-k3 k3 k4];

dx = @(t,x)[x(4:6); (-A * x(1:3) [P0 * sin( w * t); 0; 0])./[m 1; m2; m3]];

Condiciones iniciales (x4~x6 son x1'~x3 ')

x0 = zero(6, 1);

Solución, dibujo directo

ode45(dx, [0 5], x0)

lstr = array fun( @ (I){ sprintf('Posición', 'S', 'Dirección', 'Nivel')