Asked by The Protoulis
on 11 Nov 2018

I have the following code trying to solve differential equation of second order:

[t,X] = ode45(@(t,X) odefun_4(t,X,u), [0 10], [0;0]);

function dx = odefun_4(t,x,u)

%xdisp(u);

m = 15;

b = 0.2;

k = 2;

%u = interp1(ut,u,t);

dx = [x(2); u/m - (b/m)*x(2) - (k/m)*x(1)];

end

Vector u is a 1x100 vector. I expect to get a 2x100 matrix as the return value. What am I doing wrong ?

Answer by madhan ravi
on 11 Nov 2018

Edited by madhan ravi
on 11 Nov 2018

Accepted Answer

EDITED

m = 15;

b = 0.2;

k = 2;

ut = linspace(0, 10);

u = 5 * sin(2 * ut) + 10.5; % input of our system - external forc

dx = @(t,x,u)[x(2); u/m - (b/m)*x(2) - (k/m)*x(1)]

for u=u

[t,X] = ode45(@(t,x)dx(t,x,u), [0 10], [0;0]); %function calling

plot(t,X)

hold on

end

madhan ravi
on 11 Nov 2018

The Protoulis
on 11 Nov 2018

Thank you very much!

madhan ravi
on 11 Nov 2018

Anytime :) make sure to accept the answer so that people know the question is solved

Sign in to comment.

Opportunities for recent engineering grads.

Apply Today
## 0 Comments

Sign in to comment.