Extract properties from the options structure created with bvpset
Syntax
Description
val = bvpget(options,'name')
extracts the value of the named property from the structure options
, returning an empty matrix if the property value is not specified in options
. It is sufficient to type only the leading characters that uniquely identify the property. Case is ignored for property names. []
is a valid options
argument.
val = bvpget(options,'name',default)
extracts the named property as above, but returns val = default
if the named property is not specified in options
. For example,
returns val = 1e-4
if the RelTol
is not specified in opts
.
Form the initial guess for bvp4c
Syntax
solinit = bvpinit(x,v)
solinit = bvpinit(x,v,parameters)
solinit = bvpinit(sol,[anew bnew])
solinit = bvpinit(sol,[anew bnew],parameters)
Description
solinit = bvpinit(x,v)
forms the initial guess for bvp4c
in common circumstances.
x
is a vector that specifies an initial mesh. If you want to solve the boundary value problem (BVP) on , then specify x(1)
as and x(end)
as . The function bvp4c
adapts this mesh to the solution, so often a guess like x = linspace(a,b,10)
suffices. However, in difficult cases, you must place mesh points where the solution changes rapidly. The entries of x
must be ordered and distinct, so if , then x(1) <>, and similarly for .
v
is a guess for the solution. It can be either a vector, or a function:
- Vector - For each component of the solution,
bvpinit
replicates the corresponding element of the vector as a constant guess across all mesh points. That is,v(i)
is a constant guess for thei
th componenty(i,:)
of the solution at all the mesh points inx
. - Function - For a given mesh point, the function must return a vector whose elements are guesses for the corresponding components of the solution. The function must be of the form
- where
x
is a mesh point andy
is a vector whose length is the same as the number of components in the solution. For example, if you use@guess
,bvpinit
calls this function for each mesh pointy(:,j) = guess(x(j))
.
solinit = bvpinit(x,v,parameters)
indicates that the BVP involves unknown parameters. Use the vector parameters
to provide a guess for all unknown parameters.
solinit
is a structure with the following fields. The structure can have any name, but the fields must be named x
, y
, and parameters
.
solinit = bvpinit(sol,[anew bnew])
forms an initial guess on the interval [anew bnew]
from a solution sol
on an interval . The new interval must be larger than the previous one, so either anew <=
a <
b <= bnew
or anew >=
a >
b >= bnew
. The solution sol
is extrapolated to the new interval. If sol
contains parameters
, they are copied to solinit
.
solinit = bvpinit(sol,[anew bnew],parameters)
forms solinit
as described above, but uses parameters
as a guess for unknown parameters in solinit
.
Create/alter boundary value problem (BVP) options structure
Syntax
options = bvpset('name1',value1,'name2',value2,...)
options = bvpset(oldopts'name1',value1,...)
options = bvpset(oldopts,newopts)
bvpset
Description
options = bvpset('name1',value1,'name2',value2,...)
creates a structure options
in which the named properties have the specified values. Any unspecified properties have default values. It is sufficient to type only the leading characters that uniquely identify the property. Case is ignored for property names.
options = bvpset(oldopts,'name1',value1,...)
alters an existing options structure oldopts
.
options = bvpset(oldopts,newopts)
combines an existing options structure oldopts
with a new options structure newopts
. Any new properties overwrite corresponding old properties.
bvpset
with no input arguments displays all property names and their possible values.
BVP Properties
No comments:
Post a Comment