f



Problem in plotting phase portraits between x1 and x2.

This is the code which i made but  is not functioning. please help to sort the problem.
clc
clear all;
syms x1(x) x2(x)
q=0.03;
omega=0.1;
u=1;
l=1/(sqrt(3));
m=1/(sqrt(3));
n=1/(sqrt(3));
mu=0.2;
alpha=0.01;
v=sqrt((2*(1+alpha))/(q+1));
beta1=(q+1)/2;
beta2=((q+1)*(3-q))/8;
A=(v^3/2)*((3/v^4)-((2*beta2)/(1+alpha)));
B=v^3/(2*(1+alpha));
C=(v^3/2)*((1/(1+alpha))+(1/(omega^2)));
D=mu/2;
a=u/((B*l^3)+(l*(m^2+n^2)*C));
b=(A*l)/(2*((B*l^3)+(l*(m^2+n^2)*C)));
c=D/((B*l^3)+(l*(m^2+n^2)*C));
d=a/b;

x1(x) = dsolve(diff(x1,x)== x2, x1(0) == 0);
 
x2(x) = dsolve(diff(x2,x)== a*x1-b*(x1^2)+c*x2, x2(0) == d);
% ezplot(x2)

% plot(x1)
%  hold on
% plot(x2)
%     grid on
%  legend('x1Sol','x2Sol')
% plot(x1)

% hold on
% plot(x2)
  % grid on
0
Barjinder
12/24/2016 7:34:03 AM
comp.soft-sys.matlab 211266 articles. 18 followers. lunamoonmoon (258) is leader. Post Follow

0 Replies
169 Views

Similar Articles

[PageSpeed] 47

Reply: