phase portrait for system of ode

조회 수: 11 (최근 30일)
nune pratyusha
nune pratyusha 2021년 2월 13일
답변: SaiDileep Kola 2021년 2월 26일
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
darova
darova 2021년 2월 13일
What is the question?

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

채택된 답변

SaiDileep Kola
SaiDileep Kola 2021년 2월 26일
Refer the answer to similar topic here

추가 답변 (0개)

카테고리

Help CenterFile Exchange에서 Programming에 대해 자세히 알아보기

제품

Community Treasure Hunt

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

Start Hunting!

Translated by