Problem #10

For the system

(a) Find the eigenvalues and eigenvectors.

(b) Classify the critical point (0,0) as to type and determine whether it is stable, asymptotically stable, or unstable.

(c) Sketch several trajectories in the phase plane and also sketch some typical graphs of x1 versus t.

(d) Use a computer to plot accurately the curves requested in part (c).

Part (a)

MATLAB can easily find the eigenvectors and eigenvalues for this system:

  >> A=[1 2;-5 -1];
  >> [V,D]=eig(A)
  V =
    -0.1690 - 0.5071i  -0.1690 + 0.5071i
     0.8452             0.8452          
  D =
     0.0000 + 3.0000i        0          
          0             0.0000 - 3.0000i
  >>

Thus, we see that the eigenvectors &xi1 and &xi2 and eigenvalues v1 and v2 are:

Part (b)

The eigenvalues are complex and purely imaginary. Thus, according to Table 9.1.1 in Section 9.1 of the text, the origin is a center and it is stable.

Part (d)

The general solution to this system is

The MATLAB code below will plot 25 curves (not necessarily distinct) for values of c1 and c2 between -2 and 2.

  >> c = [-2 -1 0 1 2];       % These are the coefficients for the solutions 
  >> tpts = [-50:0.05:10];
  >> % Open figure, add labels and title 
  >> figure; hold on;
  >> xlabel('x_{1}(t)'); ylabel('x_{2}(t)');
  >> title('Phase Portrait for Problem #10 in Section 9.1');
  >> % Loop to use c-values in each slot of general solution 
  >> for i = 1:5
  for j = 1:5
  strx1=sprintf('%d*(-0.1690*cos(3*t) + 0.5071*sin(3*t)) + %d*(-0.1690*sin(3*t)-0.5071*cos(3*t))',c(i),c(j));
  strx2=sprintf('%d*(0.8452*cos(3*t)) + %d*(0.8452*sin(3*t))',c(i),c(j));
  x1=inline(strx1,'t');
  x2=inline(strx2,'t');
  plot(x1(tpts),x2(tpts));
  end
  end
  >> axis([-3 3 -3 3]);        % Re-size window 
  >>

We can also plot a few representative graphs of x1 versus t.

  >> % Open figure, add labels and title 
  >> figure; hold on;
  >> xlabel('t'); ylabel('x_{1}(t)');
  >> title('x_{1}(t) versus t for Problem #10 in Section 9.1');
  >> % These are the coefficients for the 8 curves to graph 
  >> c1=[1 1 -1 -1  3 3 -3 -3];
  >> c2=[1 -1 -1 1 1 -1 1 -1];
  >> % Loop to graph the 8 curves
  >> for i = 1:8
  strx1=sprintf('%d*(-0.1690*cos(3*t) + 0.5071*sin(3*t)) + %d*(-0.1690*sin(3*t)-0.5071*cos(3*t))',c1(i),c2(j));
  x1=inline(strx1);
  plot(tpts, x1(tpts));
  end
  >> axis([-2*pi 2*pi -2 2]);             % Re-size window 
  >>

As we expected, the curves do not decay as t approaches infinity, but the amplitudes remain constant.