필터 지우기
필터 지우기

Matlab tool to plot 3d phase portrait

조회 수: 47 (최근 30일)
Shreshtha Chaturvedi
Shreshtha Chaturvedi 2024년 5월 11일
답변: Sam Chak 2024년 5월 13일
I have the following system of differential equations
x' = -alpha0 * x + beta0 * x * y;
y' = alpha1 * y - beta1 * x * y;
z' = -alpha2 * (z - 5.5)^3 - beta2 * x * y;
I wanted to take a look at the phase portrait of this system for various values of the coefficients. Is there an inbuilt matlab tool that can do that?

채택된 답변

Sam Chak
Sam Chak 2024년 5월 13일
If you want to use multiple slider widgets for alpha and beta to visualize the interactive 3D phase portrait, then you need to use the Graphical User Interface Design Environment or the App designer to create your desired graphical user interface (GUI).
help GUIDE
Graphical user interface design environment MATLAB Version 24.1 (R2024a) 19-Nov-2023 GUIDE functions. guide - Open the GUI Design Environment. Documentation for matlab/guide doc guide

추가 답변 (1개)

Steven Lord
Steven Lord 2024년 5월 12일
Use odeset to set the OutputFcn option to @odephas3. See the ballode or orbitode example files for a demonstration of how to use the OutputFcn option, though I believe those use @odeplot and @odephas2 instead of @odephas3.
  댓글 수: 3
Sam Chak
Sam Chak 2024년 5월 12일
@Steven Lord suggested calling '@odephas3' to generate the 3D phase plot. Is this what you were looking for? The 'quiver3' function is used for plotting a 3D vector field, which can be challenging for humans to visualize when projected on a static 2D plane. However, on your MATLAB machine, you can manually rotate the view using your mouse.
initX = [0 0 0
0 0 1
1 0 0
1 0 1
1 1 0
1 1 1];
tspan = [0 10];
X0 = initX(3,:); % test different initial values
opts = odeset('OutputFcn', @odephas3); % <-- insert this
[t, X] = ode45(@ode, tspan, X0, opts);
%% Differential Equations
function dX = ode(t, X)
dX = zeros(3, 1);
x = X(1);
y = X(2);
z = X(3);
alpha0 = 1;
beta0 = 1;
alpha1 = 1;
beta1 = 1;
alpha2 = 1;
beta2 = 1;
dX(1) = -alpha0 * x + beta0 * x * y;
dX(2) = alpha1 * y - beta1 * x * y;
dX(3) = -alpha2 * (z - 5.5)^3 - beta2 * x * y;
end
Shreshtha Chaturvedi
Shreshtha Chaturvedi 2024년 5월 12일
편집: Shreshtha Chaturvedi 2024년 5월 12일
I see. However, I just wanted to see how the solutions behave via the phase portrait, not plot the actual solutions by solving it using ode45. As far as plotting the solutions are concerned, I have separately used Euler's method to solve the ode, as one cannot specify a particular time step in ode45. I am just plotting that solution using plot3, but I will check this out too. My original doubt however, is whether there exists any interactive tool, a widget of sorts, through which we can change the coefficients and obtain various phase portraits.

댓글을 달려면 로그인하십시오.

카테고리

Help CenterFile Exchange에서 Surfaces, Volumes, and Polygons에 대해 자세히 알아보기

제품


릴리스

R2024a

Community Treasure Hunt

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

Start Hunting!

Translated by