Search
Code Directory
 ASP
 ASP.NET
 C/C++
 CFML
 CGI/PERL
 Delphi
 Development
 Flash
 HTML
 Java
 JavaScript
 Pascal
 PHP
 Python
 SQL
 Tools
 Visual Basic & VB.NET
 XML
New Code
.Net Runtime Library for Delphi 6.0.4.0
Scimbo 1.64
AnyMap JS Maps 8.4.2
GetOrgChart 2.5.3
AnyChart JS Charts and Dashboards 8.4.2
OrgChart JS 3.8.0
dbForge Compare Bundle for MySQL 8.1
dbForge Search for SQL Server 2.3
Database Workbench Pro 5.5.0
Luxand FaceSDK 7.0
SSIS Data Flow Components 1.10
Entity Developer Professional 6.3
dbForge Index Manager for SQL Server 1.9
dbForge Data Generator For MySQL 2.2
Magento Australia Post eParcel Extension 1.0
Top Code
PHP Point of sale 10.0
MATLAB Support Package for Arduino (aka ArduinoIO Package) 1.0
Betting system 6.x-1.x-dev
Java-2-Pseudo 1.0
Faculty Evaluation System 1.1
TeeBI for RAD Studio Suite 2017
Cuckoo Search (CS) Algorithm 1.0
Student Information Management System 1.0
Java/RTR 1.0
JEDI Database Desktop 27012002
Mind Fighter 1.1
JAC (Java Asn.1 Compiler) 3.0
000-516 Free Test Exam Questions 10.0
ICDoctorAppointment - Doctor Appointment Script 1.2
CONTRAST CONTROLLER 1.0
Top Rated
Deals and Discounts Website Script 1.0.2
ADO.NET Provider for ExactTarget 1.0
Solid File System OS edition 5.1
Classified Ad Lister 1.0
Aglowsoft SQL Query Tools 8.2
ICPennyBid Penny Auction Script 4.0
PHP Review Script 1.0
ATN Resume Finder 2.0
ATN Site Builder 3.0
Availability Booking Calendar PHP 1.0
PHP GZ Blog Script 1.1
ATN Jobs Software 4.0
ATN Mall 2.0
WeBuilder 2015 13.3
PHP Digital Download Script 1.0.4
Calculation Lyapunov Exponents for ODE 1.0
File ID: 80884






Calculation Lyapunov Exponents for ODE 1.0
Download Calculation Lyapunov Exponents for ODE 1.0http://www.mathworks.comReport Error Link
License: Shareware
File Size: 10.0 KB
Downloads: 13
Submit Rating:
Calculation Lyapunov Exponents for ODE 1.0 Description
Description: Lyapunov exponent calcullation for ODE-system. The alogrithm employed in this m-file for determining Lyapunov exponents was proposed in A. Wolf, J. B. Swift, H. L. Swinney, and J. A. Vastano, "Determining Lyapunov Exponents from a Time Series," Physica D, Vol. 16, pp. 285-317, 1985.

For integrating ODE system can be used any MATLAB ODE-suite methods.

This function is a part of MATDS program - toolbox for dynamical system investigation
See: http://www.math.rsu.ru/mexmat/kvm/matds/

Input parameters:
n - number of equation
rhs_ext_fcn - handle of function with right hand side of extended ODE-system. This function must include RHS of ODE-system coupled with variational equation (n items of linearized systems, see Example).
fcn_integrator - handle of ODE integrator function, for example: @ode45
tstart - start values of independent value (time t)
stept - step on t-variable for Gram-Schmidt renormalization procedure.
tend - finish value of time
ystart - start point of trajectory of ODE system.
ioutp - step of print to MATLAB main window. ioutp==0 - no print, if ioutp>0 then each ioutp-th point will be print.

Output parameters:
Texp - time values
Lexp - Lyapunov exponents to each time value.

Users have to write their own ODE functions for their specified systems and use handle of this function as rhs_ext_fcn - parameter.

Example. Lorenz system:
dx/dt = sigma*(y - x)
dy/dt = r*x - y - x*z
dz/dt = x*y - b*z

The Jacobian of system:
| -sigma sigma 0 |
J = | r-z -1 -x |
| y x -b |

Then, the variational equation has a form:
F = J*Y
where Y is a square matrix with the same dimension as J.
Corresponding m-file:
function f=lorenz_ext(t,X)
SIGMA = 10; R = 28; BETA = 8/3;
x=X(1); y=X(2); z=X(3);
Y= [X(4), X(7), X(10);
X(5), X(8), X(11);
X(6), X(9), X(12)];
f=zeros(9,1);
f(1)=SIGMA*(y-x);
f(2)=-x*z+R*x-y;
f(3)=x*y-BETA*z;
Jac=[-SIGMA,SIGMA,0;
R-z,-1,-x; y, x,-BETA];
f(4:12)=Jac*Y;

% Run Lyapunov exponent calculation:

[T,Res]=lyapunov(3,@lorenz_ext,@ode45,0,0.5,200,[0 1 0],10);

See files: lyapunov.m - algorithm
lorenz_ext - Lorenz system extended rhs
run_lyap - example of calling and result visualization
Govorukhin V.N.
This file is intended for use with MATLAB and was produced for MATDS-program http://www.math.rsu.ru/mexmat/kvm/matds/ lyapunov.m is free software. lyapunov.m is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY.

License: Shareware

Related: jacobian, sigma, sigmay, lorenz, Parameter, Square, dlorenz exttx

O/S:BSD, Linux, Solaris, Mac OS X

File Size: 10.0 KB

Downloads: 13



More Similar Code

XVidCalc is a video bitrate calculator to provide quick calculation of bitrates for multimedia files for optimizing quality and filesize. PHP and X11/QT versions available.



The toolbox draws:
phase portrait,
vector field
Lyapunov exponents,
graphic solution
bifurcation diagram,
iteration graph
time series
Iterated function.

Are analyzed:
fixed points, periodic...



For longer operations on your computer, it is practical to inform the user that a calculation is running. For this reason, you can call the 'busyWindow' function. It displays your message to the user.
There are three modes: The busy window...



ode_solver provides a graphical user interface for solving ODE systems with additive noise influence.
Also it is possible to run solver in no-GUI (batch) mode. See required parameters in help-file.

example for Duffing's...



This contribution is made in order to calculate the maximum likelihood contour plots of any function. Just replace the <<fcnObjfun>> file with your function change some basic settings in the <<Example>> file and there you...



RKN1210 12th/10th order Runge-Kutta-Nystrom integrator

RKN1210() is a 12th/10th order numerical integrator for ordinary differential equations of the form

y'' = f(t, y) (1)

with initial conditions

...



A polynomial is a mathematical expression involving a sum of powers in one or more variables multiplied by coefficients. A general multivariate polynomial is captured with the syntax

p = sum( c_i * prod( x_j^p_ij ) ) + k



for derivatives calculation use, a chinese book----I wrote.



This is a GUI for Riccati formulation and H-infinity and H-2 Norms calculation for Robust Control.



for personal use only.
1: the interest rate structure calculation and rules.
2: the interest rate model, used frequently in real world.

User Review for Calculation Lyapunov Exponents for ODE
- required fields
     

Please enter text on the image