% Exemplo Capítulo 14 Pêndulo

clear all
close all
clc

m = 1; % massa do objeto
l = 1; % comprimento da haste
b = 1;  % constante da amortecimento (resistência do fluido)
g = 9.8; % aceleração da gravidade
a = 10*pi/180; % posição inicial (theta(0)=10º) em radianos

wn = sqrt(g/l);
xi = b/(2*m*wn);
wd = wn*sqrt(1-xi^2);

x = [];
y = [];
theta = [];
ind = 0;
for t = 0:0.1:12
    ind = ind+1;
    Theta = a*exp(-xi*wn*t)*(cos(wd*t)+(xi/sqrt(1-xi^2))*sin(wd*t));
    xp = l*sin(Theta);
    yp = -l*cos(Theta);
    theta = [theta, Theta];
    x = [x, xp];
    y = [y, yp];
    
    figure
    % desenhando a haste do pêndulo:
    plot([0 xp],[0 yp],'k')  
    hold on
    % desenhando a massa do pêndulo:
    plot(xp,yp,'ko','MarkerSize',15,'MarkerFaceColor','k') 
    axis([-0.3 0.3 -1.3 0])
    nome = strcat('Pend_b_',num2str(b),'_',num2str(ind));
    print( figure(1), '-djpeg', nome);
    close all

    
end

t = 0:0.1:12;
plot(t,theta)
xlabel('t')
ylabel('\theta(t)')