Scribd Logo
Upload

Welcome to Scribd!

  • 71 views

    Integration of Torus.m Plot The Poincare Section Torus.m File - Two Uncoupled Oscilators With Omega1 2 and Omega2 3

    Uploaded by

    Poornima Sridaran
    The document describes calculating and plotting the Poincare section of two uncoupled oscillators where the ratio of their frequencies is rational, resulting in closed orbits and periodic behavior, versus irrational where it is quasiperiodic. It integrates the system, calculates the Poincare section points where one variable crosses a multiple of 2pi, and plots the section points.

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd

    Integration of Torus.m Plot The Poincare Section Torus.m File - Two Uncoupled Oscilators With Omega1 2 and Omega2 3

    Uploaded by

    Poornima Sridaran
    71 views3 pages
    The document describes calculating and plotting the Poincare section of two uncoupled oscillators where the ratio of their frequencies is rational, resulting in closed orbits and periodic behavior, versus irrational where it is quasiperiodic. It integrates the system, calculates the Poincare section points where one variable crosses a multiple of 2pi, and plots the section points.

    Original Description:

    how to plot poin care map

    Original Title

    Poincare Plot

    Copyright

    © © All Rights Reserved

    Available Formats

    PDF, TXT or read online from Scribd

    Did you find this document useful?

    Is this content inappropriate?

    The document describes calculating and plotting the Poincare section of two uncoupled oscillators where the ratio of their frequencies is rational, resulting in closed orbits and periodic behavior, versus irrational where it is quasiperiodic. It integrates the system, calculates the Poincare section points where one variable crosses a multiple of 2pi, and plots the section points.

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd
    Download as pdf or txt
    71 views3 pages

    Integration of Torus.m Plot The Poincare Section Torus.m File - Two Uncoupled Oscilators With Omega1 2 and Omega2 3

    Uploaded by

    Poornima Sridaran
    The document describes calculating and plotting the Poincare section of two uncoupled oscillators where the ratio of their frequencies is rational, resulting in closed orbits and periodic behavior, versus irrational where it is quasiperiodic. It integrates the system, calculates the Poincare section points where one variable crosses a multiple of 2pi, and plots the section points.

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd
    Download as pdf or txt
    You are on page 1of 3

    1/25/2018 Calculation and plot of the poincare section of two uncoupled osciators with omega1 and omega2 : omega1/omega2=p/q=rational

    -> closed orbi…

    Calculation and plot of the poincare section of two


    uncoupled osciators with omega1 and omega2 :
    omega1/omega2=p/q=rational -> closed orbit, periodic
    behaviour; omega1/omega2=irational -> quasiperiodic
    behaviour
    Contents

    integration of torus.m
    calculation of poincare section
    plot the poincare section
    torus.m file - two uncoupled oscilators with omega1=2 and omega2=3

    close;
    clear;
    ps=zeros(2,1);

    integration of torus.m

    [t,y]=ode45(@torus,[0:0.001:100],[0.5 0.5 0.5]);


    %plot projection of the phas space orbit to y1-y2 plane
    plot(y(:,1),y(:,2));
    hold;

    Current plot held

    http://lab.fs.uni-lj.si/lasin/www/teaching/chaotic/Poincare/poincare1.html#3 1/3
    1/25/2018 Calculation and plot of the poincare section of two uncoupled osciators with omega1 and omega2 : omega1/omega2=p/q=rational -> closed orbi…

    calculation of poincare section

    calculate number of integration steps

    n=size(t);
    %set the index of poincare points to 1
    np=1;
    for i=1:n(1)
    % detect the cros-section of the trajectory with the plane y1-y2
    if(y(i,3)>=(2*pi)*np)
    % store detected cross-section point y1,y2 to ps1,ps2
    ps(np,1)=y(i,1);
    ps(np,2)=y(i,2);
    % increase the index of poincare point
    np=np+1;
    end
    end

    plot the poincare section

    for i=1:np-1
    plot(ps(i,1),ps(i,2),'r+')
    % use pause to folow the plot of the poincare section
    %pause(2);
    end
    %plot(ps(:,1),ps(:,2),'r+');

    torus.m file - two uncoupled oscilators with omega1=2 and omega2=3

    %function dy=torus(t,y)
    %dy=zeros(3,1);
    %dy(1)=y(2);

    http://lab.fs.uni-lj.si/lasin/www/teaching/chaotic/Poincare/poincare1.html#3 2/3
    1/25/2018 Calculation and plot of the poincare section of two uncoupled osciators with omega1 and omega2 : omega1/omega2=p/q=rational -> closed orbi…
    %dy(2)=-(omega1)*y(1);
    %dy(3)=omega2;

    Published with MATLAB® 7.5

    http://lab.fs.uni-lj.si/lasin/www/teaching/chaotic/Poincare/poincare1.html#3 3/3

    You might also like