How to solve the following System of first order differential equations using ode45?please help

2 views (last 30 days)
q1=[q11; q21; q13];
q10=[7.1892; 11.4189; 0.5826];
t1 = pi:0.1:2*pi;
dq11dt1= -k1*(12.7734-q11)*cos(q31);
dq21dt1=-k1*(12.7734-q11)*sin(q31);
dq31dt1= -k2*(0.3097+(13.6167-q21)*sin(t1));

Answers (1)

Aquatris
Aquatris on 7 Sep 2018
Edited: Aquatris on 7 Sep 2018
First you create your function that outputs derivative of q when time and q are given to it.
function qd = asd(t,q)
k1 = 1.2;
k2 = 1.3;
k3 = 1.5;
qd = [-k1*(12.7734-q(1))*cos(q(3));
-k1*(12.7734-q(1))*sin(q(3));
-k2*(0.3097+(13.6167-q(2))*sin(t))];
end
Then in the main script, you call
q10=[7.1892; 11.4189; 0.5826];
t1 = pi:0.1:2*pi;
[t,q] = ode45(@asd,t1,q10);
You did not specify k1 k2 k3 so I randomly selected them.

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!