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
Rapid PHP 2018 15.5
IBDAC 6.4
Online Course Booking Script 1.3.3
Database Workbench Pro 5.6.4
Job Portal Script 1.3.2
The C# PDF Library 5.2
PostgreSQL Data Access Components 6.0
Virtual Data Access Components 11.0
HTMLPad 2018 15.5
WeBuilder 2018 15.5
SentiMask SDK Trial 2.0.0
Track Order For Magento 2 1.0.0
Calendar 365 For Dynamics CRM 4.0
Scimbo 1.2
Odoo Furnito Theme 1.0
Top Code
Hotel Management System Restaurant Version 5.25
Azizi search engine script PHP 4.1.10
Travel Booking script 1.0.4
Chessgra 1.1
Keno Game 1.I
VChess 6.x-1.x-dev
Killer Guestbook Script 1.1
RLS Algorithm Demo 1.0
PHP MLM Software 2.0.1
College Management System Script 1.0.4
Spider(rader) Chart for Multi Data 1.0
Simple Radar_system_simulation 1.0
The Spherical K-means algorithm 1.0
Discrete Frechet Distance 1.0
AJ Matrix v2.6
Top Rated
VisualNEO Web 2018.12.15
Azizi search engine script PHP 4.1.10
Paste phpSoftPro 1.4.1
Extreme Injector 3.7
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
Invoice Manager by PHPJabbers 3.0
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
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