How to solve the following System of first order differential equations using ode45?please help
1 view (last 30 days)
Show older comments
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));
0 Comments
Answers (1)
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.
1 Comment
sangita kamat
on 8 Sep 2018
Edited: sangita kamat
on 8 Sep 2018
@Aquatris, thank you very much for your help. Shall do.its perfectly working . thanks once again
See Also
Categories
Find more on Ordinary Differential Equations in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!