phase portrait for system of ode

9 vues (au cours des 30 derniers jours)
nune pratyusha
nune pratyusha le 13 Fév 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
  1 commentaire
darova
darova le 13 Fév 2021
What is the question?

Connectez-vous pour commenter.

Réponse acceptée

SaiDileep Kola
SaiDileep Kola le 26 Fév 2021
Refer the answer to similar topic here

Plus de réponses (0)

Catégories

En savoir plus sur Programming dans Help Center et File Exchange

Produits

Community Treasure Hunt

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

Start Hunting!

Translated by