phase plane in MATLAB?

2 views (last 30 days)
JIE YUANG LU
JIE YUANG LU on 10 Nov 2016
Edited: JIE YUANG LU on 11 Nov 2016
Hi I have a equations here, and I wonder it to do phase plane.How to do it?

Accepted Answer

Mischa Kim
Mischa Kim on 10 Nov 2016
Hello, check out the code below:
function my_ODE()
th0 = 0;
Dth0 = 1;
tspan = linspace(0,2,200);
options = odeset('RelTol',1e-8,'AbsTol',1e-10);
[~,X] = ode45(@DE, tspan,[th0; Dth0],options);
plot(X(:,1),X(:,2))
grid
end
function dX = DE(~,x)
dX = [x(2); u(x(1))];
end
function uval = u(x1)
if (x1 < 0)
uval = +5;
else
uval = -5;
end
end
  1 Comment
JIE YUANG LU
JIE YUANG LU on 11 Nov 2016
Edited: JIE YUANG LU on 11 Nov 2016
so when I change U condition which code to change? and what is DE?

Sign in to comment.

More Answers (0)

Categories

Find more on MATLAB in Help Center and File Exchange

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!