How to solve the following System of first order differential equations using ode45?please help
    3 views (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 Numerical Integration and 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!

