Group Abstract

Message Boards

WOLFRAM COMMUNITY

7.9K Views

4 Replies

2 Total Likes

View groups...

Follow this post

Share this post:

GROUPS:

Mathematics Equation Solving Wolfram Language

Solve a 3x3 ODE system?

Ermanno Citraro

Posted 6 years ago

POSTED BY: Ermanno Citraro

4 Replies

Sort By:

Alexander Trounev

Alexander Trounev, A&E Trounev IT Consulting

Posted 6 years ago

This can be used as a regular analytical solution, for example s = DSolve[{i[t] == Cr vc'[t], vc[t] == vco[t] (1 + m) - Lr i'[t], -(1 + m) i[t] - vco[t]/RL == Co vco'[t], vc[0] == Vcrmax, i[0] == 0, vco[0] == Vo}, {i, vc, vco}, t]; {Plot[Evaluate[ i[t] /. s[[1, 1]] /. {Cr -> 1, Lr -> 1, m -> 1, Co -> 1, RL -> 1, Vcrmax -> 1, Vo -> 1}], {t, 0, 2 Pi}, AxesLabel -> {"t", "i"}], Plot[Evaluate[ vc[t] /. s[[1, 2]] /. {Cr -> 1, Lr -> 1, m -> 1, Co -> 1, RL -> 1, Vcrmax -> 1, Vo -> 1}], {t, 0, 2 Pi}, AxesLabel -> {"t", "vc"}], Plot[Evaluate[ vco[t] /. s[[1, 3]] /. {Cr -> 1, Lr -> 1, m -> 1, Co -> 1, RL -> 1, Vcrmax -> 1, Vo -> 1}], {t, 0, 2 Pi}, AxesLabel -> {"t", "vco"}]}

This can be used as a regular analytical solution, for example

s = DSolve[{i[t] == Cr vc'[t], 
    vc[t] == vco[t] (1 + m) - Lr i'[t], -(1 + m) i[t] - vco[t]/RL == 
     Co vco'[t], vc[0] == Vcrmax, i[0] == 0, vco[0] == Vo}, {i, vc, 
    vco}, t];

{Plot[Evaluate[
   i[t] /. s[[1, 1]] /. {Cr -> 1, Lr -> 1, m -> 1, Co -> 1, RL -> 1, 
     Vcrmax -> 1, Vo -> 1}], {t, 0, 2 Pi}, AxesLabel -> {"t", "i"}], 
 Plot[Evaluate[
   vc[t] /. s[[1, 2]] /. {Cr -> 1, Lr -> 1, m -> 1, Co -> 1, RL -> 1, 
     Vcrmax -> 1, Vo -> 1}], {t, 0, 2 Pi}, AxesLabel -> {"t", "vc"}], 
 Plot[Evaluate[
   vco[t] /. s[[1, 3]] /. {Cr -> 1, Lr -> 1, m -> 1, Co -> 1, RL -> 1,
      Vcrmax -> 1, Vo -> 1}], {t, 0, 2 Pi}, 
  AxesLabel -> {"t", "vco"}]}

POSTED BY: Alexander Trounev

Gianluca Gorni

Gianluca Gorni, University of Udine

Posted 6 years ago

They are algebraic objects that are comlicated, because they involve equations of degree 3. Try with numerical values for the parameters: With[{Cr = 1, Lr = 1, Co = 1, m = 1, RL = 1, Vcrmax = 1, Vo = 1}, DSolveValue[{i[t] == Cr vc'[t], vc[t] == vco[t] (1 + m) - Lr i'[t], -(1 + m) i[t] - vco[t]/RL == Co vco'[t], vc[0] == Vcrmax, i[0] == 0, vco[0] == Vo}, {i[t], vc[t], vco[t]}, t]] // N

They are algebraic objects that are comlicated, because they involve equations of degree 3. Try with numerical values for the parameters:

With[{Cr = 1, Lr = 1, Co = 1, m = 1, RL = 1, Vcrmax = 1, Vo = 1}, 
  DSolveValue[{i[t] == Cr vc'[t], 
    vc[t] == vco[t] (1 + m) - Lr i'[t], -(1 + m) i[t] - vco[t]/RL == 
     Co vco'[t], vc[0] == Vcrmax, i[0] == 0, vco[0] == Vo}, {i[t], 
    vc[t], vco[t]}, t]] // N

POSTED BY: Gianluca Gorni

Ermanno Citraro

Posted 6 years ago

Thanks sir.

POSTED BY: Ermanno Citraro

Ermanno Citraro

Posted 6 years ago

Grazie Gianluca.

POSTED BY: Ermanno Citraro

Reply to this discussion

Reply Preview

Attachments

Remove Add a file to this post

Follow this discussion

or Discard

Feedback