phase portrait for system of ode
9 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
nune pratyusha
le 13 Fév 2021
Réponse apportée : SaiDileep Kola
le 26 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
Réponse acceptée
Plus de réponses (0)
Voir également
Produits
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!