Code covered by the BSD License  

Highlights from
Dynamical Systems with Applications using MATLAB

image thumbnail

Dynamical Systems with Applications using MATLAB

by

 

13 Sep 2002 (Updated )

Companion Software.

Editor's Notes:

This file was selected as MATLAB Central Pick of the Week

Programs_10c.m
% Chapter 10 - Limit Cycles.
% Programs_10c - Phase Portrait (Fig. 10.7).
% Copyright Birkhauser 2013. Stephen Lynch.

% Non-convex limit cycle of a Lienard system.
clear
hold on
a2=90;a4=-882;a6=2598.4;a8=-3359.997;a10=2133.34;a12=-651.638;a14=76.38;
sys1=@(t,x) [x(2);-x(1)-x(2)*(a2*(x(1))^2+a4*(x(1))^4+a6*(x(1))^6+a8*(x(1))^8+a10*(x(1))^10+a12*(x(1))^12+a14*(x(1))^14)]; 
[t,xt] = ode45(sys1,[0 30],[1.4 0]);
plot(xt(:,1),xt(:,2))      
hold off
axis([-3 3 -3 3])
fsize=15;
set(gca,'XTick',-3:1:3,'FontSize',fsize)
set(gca,'YTick',-3:1:3,'FontSize',fsize)
xlabel('x(t)','FontSize',fsize)
ylabel('y(t)','FontSize',fsize)
hold off

% End of Programs_10c.

Contact us