phase portrait for system of ode
11 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
nune pratyusha
el 13 de Feb. de 2021
Respondida: SaiDileep Kola
el 26 de Feb. de 2021
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
my output is coming fizzy but i want one single elliptical nature phase portrait
Respuesta aceptada
Más respuestas (0)
Ver también
Productos
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!