No BSD License  

Highlights from
Forced Duffing Oscillator

4.28571
4.3 | 8 ratings Rate this file 58 Downloads (last 30 days) File Size: 24.9 KB File ID: #16731 Version: 1.0
image thumbnail

Forced Duffing Oscillator

by

Housam Binous (view profile)

 

06 Oct 2007 (Updated )

Computes Poincaré section of the forced Duffing oscillator

| Watch this File

File Information
Description

The forced Duffing oscillator exhibits behavior ranging from limit cycles to chaos due to its nonlinear dynamics. When the periodic force that drives the system is large, chaotic behavior emerges and the phase space diagram is a strange attractor. In that case the behavior of the system is sensitive to the initial condition. In order to plot a Poincaré section, take one data point from phase space per period of the driving force. The Poincaré section is a complicated fractal curve when the phase diagram is a strange attractor. The Poincaré section is a single point when the phase space diagram is a limit cycle.

Here is a link to a Mathematica 6.0 Demonstration concerning the forced Duffing oscillator:

http://demonstrations.wolfram.com/ForcedDuffingOscillator/

MATLAB release MATLAB 7.0.1 (R14SP1)
Tags for This File   Please login to tag files.
Please login to add a comment or rating.
Comments and Ratings (9)
28 Apr 2015 Akash kumar

Why this code is not running ?
there is some problem which says ''Error using feval
Undefined function 'duffing' for input arguments of type 'double''.What should I do,any help will be appeciated.

08 Sep 2012 Woo Sang Pyo

thank you for your efforts
I have a question
in your m-file xdot=duffing(t,x), why omega^2*x(2) term is plus sign not negative sign?

Comment only
08 Sep 2012 Woo Sang Pyo  
17 Nov 2010 sabara thinam

i'm very happy to using this programme it's very nice

17 Nov 2010 sabara thinam

mr.housam binous done a wonderfull making of duffing forced oscillator.. i'm using this in this programme is very usfull my studies...

07 Oct 2010 Dalhousie Moghtadayi

thank you for the source codes;
they are really helpful;

05 May 2008 Przemyslaw Jasinski

very good, helped me a lot.

13 Apr 2008 Nicholas Fantuzzi

The m-file doesn't work with epsilon < 0

24 Jan 2008 R RAMLAN

very helpful and easy to understand

Contact us