% solução numérica de sistema com 1 gdl
% rigidez cúbica, sem excitação e com excitação harmônica

clc
clearvars
close all

% vetor de tempo
  % sem excitação
% ts=[0 100];
% ts=linspace(0,100,500);
  % com excitação harmônica
    % alteração na frequência
% ts=linspace(0,100,500); % periódica
% ts=0:(2*pi/1.29):100*(2*pi/1.29); % caótica
% ts=0:(2*pi/1.4):100*(2*pi/1.4); % periódica (de onde até onde?)
    % alteração na amplitude (exercício)
% ts=linspace(0,100,500); 

% condições iniciais
  % sem excitação
% ci=[1.0 0.52]; % 0.52 0.54 0.56
  % com excitação
% ci=[0 0]; % alteração na frequência
ci=[1 0]; % alteração na amplitude (exercício)

[t,x]=ode45('rceh',ts,ci);

% -b: sem excitação; .b: com excitação (freq.); -b: com excitação (amp.)
plot(x(:,1),x(:,2),'-b') 
ylabel('velocidade (m/s)')
xlabel('deslocamento (m)')
title('Plano de fase')
axis([-2 2 -1 1])
grid on