Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

MatLab help using ode23/ode45

  1. Mar 8, 2004 #1
    I want to use matlab ode23 or ode45 to solve the following problem,
    (ignore the periods '.' as they are there to line up the vector spaces)

    .... Ids ....0.5Ids-0.3Iqs-0.4Idr-0.1IdrIqr ........ 0.2Vd
    d ..Iqs = -0.3Ids+0.5Iqs+0.8Idr-0.2Iqr^2 . . + . 0.2Vds
    dt .Idr ...-0.5Ids+0.1Iqs+0.7Idr+0.6IdrIqr ......... 0
    .. . Iqr ... 0.9Iqs-0.5Iqs-0.4Idr+0.3Iqr^2 .......... 0
    .... Wm .. 0.27(IdrIqs-IqrIds)-0.1Wm-0.85 ........ 0

    Any help explaining the correct input syntax would be greatly appreciated (as I'm not that familiar with MatLab).

  2. jcsd
  3. Mar 8, 2004 #2

    correction original post....should say -0.1Iqr not -0.1IdrIqr,
    -0.2Iqr not -0.2Iqr^2, 0.6Iqr not 0.6IdrIqr, 0.3Iqr not 0.3 Iqr^2
  4. Mar 8, 2004 #3


    User Avatar
    Staff Emeritus
    Science Advisor
    Gold Member

    Hi louie,

    I'm not sure I understand what you're trying to do.

    For ode45, you need two files. One which is your differential updater, and one which is the main function which calls it.

    the differential function should return a vector which is the derivatives and read in time steps you want to look at and the original state vector.

    The main file should set up the state vector and pass it.

    for example, the differential function:

    Code (Text):

    function dW=differential(t,W)

    for the calling function:
    Code (Text):

    function [t,W]=odesolver("input variables")

    W_0= [X_0;Y_0];


    This will return the values for W at the points designated by t
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook