phase portrait for system of ode

Question

nune pratyusha 2021-2-13

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/744347-phase-portrait-for-system-of-ode

回答： SaiDileep Kola 2021-2-26

采纳的回答： SaiDileep Kola

Untitled.png

在 MATLAB Online 中打开

dx/dt=y

dy/dt=-(1/r+g+a+b*|x|)y/c-x/(lc)

for this system i am finding phase portraits ,my code is

clc
clear all
%data:
 r = 1430; 
    a = -0.0683; 
    b = 0.0676;
    c = 36*10^-9;
    g = -0.0676;
    l=27
[x,y]=meshgrid(-2:0.01:2.05, -2:0.01:2.05);
dx = y;
dy = -(1/c)*((1/r)+g+a+b*abs(x)).*y-x/(l*c)
streamslice(x, y, dx, dy,'filled');
% hold on
% figure
quiver(x, y, dx, dy);
% plot(dx,dy)
% figure
title('Phase portrait')
% axis tight