- #1

- 719

- 1

## Main Question or Discussion Point

Hi Folks,

I have made up this code in an attempt to find the response time a 4 dof problem. However, I do not get an oscillatory response as I expect, instead I get a straight line which drops off towards the end. Not sure what is happening..any ideas..?

function f = Simple_Planetary(t,y)

f = zeros(8, 1);

f(1) = y(2);

f(2) = -y(7)+2*y(1);

f(3) = y(4);

f(4) = -y(7)+2*y(3);

f(5) = y(6);

f(6) = -y(7)+2*y(5);

f(7) = y(8);

f(8) = 3*y(7)-y(1)-y(2)-y(3)+10*cos(t);

y0 = [1;1;1;1;1;1;1;1];

[t,y] = ode45('Simple_Planetary', tspan, y0);

subplot (211)

plot (t,y (: , 1));

xlabel ('t');

ylabel ('x1 (t)');

subplot (212)

plot (t,y (: , 3));

xlabel ('t');

ylabel ('x2 (t)');

Regards

I have made up this code in an attempt to find the response time a 4 dof problem. However, I do not get an oscillatory response as I expect, instead I get a straight line which drops off towards the end. Not sure what is happening..any ideas..?

function f = Simple_Planetary(t,y)

f = zeros(8, 1);

f(1) = y(2);

f(2) = -y(7)+2*y(1);

f(3) = y(4);

f(4) = -y(7)+2*y(3);

f(5) = y(6);

f(6) = -y(7)+2*y(5);

f(7) = y(8);

f(8) = 3*y(7)-y(1)-y(2)-y(3)+10*cos(t);

y0 = [1;1;1;1;1;1;1;1];

[t,y] = ode45('Simple_Planetary', tspan, y0);

subplot (211)

plot (t,y (: , 1));

xlabel ('t');

ylabel ('x1 (t)');

subplot (212)

plot (t,y (: , 3));

xlabel ('t');

ylabel ('x2 (t)');

Regards