¡Master resuelve ecuaciones diferenciales dinámicas usando la programación ode45 en matlab! Gracias

Código de referencia:

Definición del parámetro

m1 = 1; m2 = 2; = 5; k3 = 6; k4 = 7;

P0 = 8; w = 9;

Ecuación diferencial

A = [k1 k2 -k1 0 ; -k2 k2 k3 -k3; 0 -k3 k3 k4];

dx = @(t,x)[x(4:6); *sin(w*t);0;0])./[m1;m2;m3]];

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

x0 = zeros(6, 1);

Resolver, dibujar directamente

ode45(dx, [0 5], x0)

lstr = arrayfun( @(i){sprintf('xi',i)},1:6);

legend(lstr{:})

legend(lstr{:},'Ubicación ', 'S', 'Orientación', 'horizontal')