15.11.2014 Views

MATLAB Mathematics - SERC - Index of

MATLAB Mathematics - SERC - Index of

MATLAB Mathematics - SERC - Index of

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

5 Differential Equations<br />

odephas3<br />

odeprint<br />

Three-dimensional phase plane plot<br />

Print to command window<br />

Introduction to Initial Value ODE Problems<br />

What Is an Ordinary Differential Equation?<br />

The ODE solvers are designed to handle ordinary differential equations. An<br />

ordinary differential equation contains one or more derivatives <strong>of</strong> a dependent<br />

variable y with respect to a single independent variable t, usually referred to<br />

as time. The derivative <strong>of</strong> y with respect to t is denoted as y′ , the second<br />

derivative as y′′ , and so on. Often yt () is a vector, having elements<br />

y 1 , y 2 ,…,<br />

y n .<br />

Types <strong>of</strong> Problems Handled by the ODE Solvers<br />

The ODE solvers handle the following types <strong>of</strong> first-order ODEs:<br />

• Explicit ODEs <strong>of</strong> the form y′ = fty ( , )<br />

• Linearly implicit ODEs <strong>of</strong> the form<br />

matrix<br />

• Fully implicit ODEs <strong>of</strong> the form<br />

Mt ( , y) ⋅ y′<br />

ftyy′ ( , , ) = 0<br />

= fty ( , ), where M(t,y) is a<br />

(ode15i only)<br />

Using Initial Conditions to Specify the Solution <strong>of</strong> Interest<br />

Generally there are many functions yt () that satisfy a given ODE, and<br />

additional information is necessary to specify the solution <strong>of</strong> interest. In an<br />

initial value problem, the solution <strong>of</strong> interest satisfies a specific initial<br />

condition, that is, y is equal to y 0 at a given initial time t 0 . An initial value<br />

problem for an ODE is then<br />

y′ = fty ( , )<br />

yt ( 0 ) = y 0<br />

(5-1)<br />

If the function ft ( , y)<br />

is sufficiently smooth, this problem has one and only one<br />

solution. Generally there is no analytic expression for the solution, so it is<br />

necessary to approximate yt () by numerical means, such as using one <strong>of</strong> the<br />

ODE solvers.<br />

5-4

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!