dop54_dop853
所属分类:matlab编程
开发工具:matlab
文件大小:172KB
下载次数:141
上传日期:2011-11-08 19:43:35
上 传 者:
water628
说明: 求解时滞常微分方程的matlab程序代码
(codes for delay differential equation)
文件列表:
Arenstorf.m (1168, 2010-01-02)
Arenstorf_dop54_dop853.m (2164, 2010-03-18)
Ballode1_dop54_dop853.m (2519, 2010-03-08)
Ballode2_dop54_dop853.m (4229, 2010-03-08)
dop54.m (33153, 2010-03-08)
dop853.m (49195, 2010-03-08)
dopget.m (1945, 2010-03-05)
dopset.m (7209, 2010-03-05)
Evol_DonneesSoleilPlanetes.m (1677, 2010-03-18)
Evol_Ephemeride 2010 Helio.txt (6225, 2010-03-18)
Evol_HighDefMinimum_dop54.m (3249, 2010-02-17)
Evol_HighDefMinimum_dop853.m (3385, 2010-02-17)
Evol_Soleilplanete_dop54_dop853.m (12150, 2010-03-18)
Evol_VitAccSoleilPlanetes.m (1749, 2010-02-12)
MassOsci.m (455, 2010-02-07)
MasstFcn.m (67, 2010-03-07)
MasstOsci.m (526, 2010-03-07)
MasstyFcn.m (96, 2010-03-07)
MasstyOsci.m (551, 2010-03-07)
Mass_dop54_dop853.m (2187, 2010-03-17)
NonNegativFcn.m (57, 2010-02-13)
NonNegativ_dop54_dop853.m (1741, 2010-02-13)
Osci.m (422, 2010-02-07)
Osci_dop54_dop853.m (6136, 2010-02-12)
VanderPol.m (549, 2010-02-12)
VanderPol_dop54_dop853.m (2278, 2010-03-18)
VanderPol_t2_y2.m (1109, 2010-03-18)
VanderPol_t2_y2_dop54_dop853.m (1122, 2010-02-17)
% --------------------------------------------------------- %
% Readme.txt %
% %
% List of the examples used to test dop54 and dop853 %
% %
% Denis Bichsel %
% 58 Rue des Deurres %
% 2000 Neuchtel %
% Tel. 41 (0) 32 730 10 16 %
% email: dbichsel@infomaniak.ch %
% --------------------------------------------------------- %
The functions dop54 and dop853 are translation in Matlab from
the Fortran code written by :
E. Hairer and G. Wanner
Université de Genève.
% --------------------------------------------------------- %
The options for dop54 and dop853 can be set with dopset
and read with dopget.
These function are very similar to Matalb function ode45,
odeset and odeget.
The allowed options are essentially the same as in Matlab.
The options
BDF
Jacobian
JConstant
JPattern
Refine
Vectorized
MvPattern
MassSingular
InitialSlope
MaxOrder
are not implemented.
Remark : The following options are not implemented for ode45
in Matlab
BDF
Jacobian
JConstant
JPattern
Vectorized
MvPattern
InitialSlope
MassSingular
InitialSlope
MaxOrder
The option "Refine" is just useless with dop54 and dop 853,
because it's possible to choose all the output value for
the t variable without any effect on the step size calculation.
DenseOuputFcn, DenseOutputSel.
Like the "OutputFcn" and "OutputSel" options, "DenseOutputFcn"
and "DenseOutputSel" allow to get the t and y values after each
step calculation. More, with this option the user gets also the
coefficients which allow dense calculation on the intervalle
[t, t+h], h = step size, 5 coefficients for dop54 and
8 coefficients for dop853 (see example below).
These options allows for example, the calculation with high
precision of the extrema of a function
% --------------------------------------------------------- %
List of examples:
------------------
The name of the main function of the examples end always
with "_dop54_dop853.m". The other functions with a name
which begin on the same way, are used for the calculation.
For example : Arenstorf_dop54_dop853.m uses the function
"Arenstorf.m"
Arenstorf_dop54_dop853.m
This function tests the following options:
RelTol
AbsTol
MaxIter
MaxStep
and compare the execution time of dop54, dop853 and
Matlab ode45
Ballode1_dop54_dop853.m
Ballode2_dop54_dop853.m
These functions test the following options:
Events
OutputFcn
OutputSel
InitialStep
MaxStep
Evol_Soleilplanete_dop54_dop853.m
This function tests the following options:
RelTol
AbsTol
InitialStep
MaxStep
Stats
DenseOutputFcn
DenseOutputSel
This function shows very clearly the difference of
dop54 and dop853 in solving the Mercury's perihelion
problem. dop54 like ode45 don't succeed but dop853 result is
531.88 arc second (true value is 531.54 in the Newton's
theory)
Mass_dop54_dop853.m
This function tests the following options:
Mass
MStateDependence
NonNegativ_dop54_dop853.m
This function tests the following option:
NonNegative
Osci_dop54_dop853.m
This very simple function compare the time needed in solving
the same problem (harmonic oscillator) by dop54, dop853 and
the Matlab function ode45.
VanderPol_dop54_dop853.m
This function tests the following options:
OutputSel
DenseOutputFcn
DenseOutputSel
InitialStep
MaxStep
% --------------------------------------------------------- %
近期下载者:
相关文件:
收藏者: