clc;
clear all;
close all;

% dados de entrada
R=88.9e-3;
L=261.62e-3;
C=8.89e-3;

% 1o. passo - valor específico de q
q=pi/3;

% cálculo de variáveis secundárias
A=asin((R*sin(q)-C)/L);
X=(R*cos(q))+(L*cos(A));

% cálculo dos coeficientes de velocidade
Ka=(R*cos(q))/(L*cos(A));
Kx=-(C+X*tan(A));

% cálculo das derivadas dos coeficientes de velocidade
La=((-R*sin(q))/(L*cos(A)))+(Ka^2*tan(A));
Lx=-((Kx*tan(A))+(X*Ka/(cos(A)^2)));

disp(['Para q = ',num2str(q),' radianos, os valores desejados são:'])
disp(' ')
disp(['A',' = ',num2str(A),' radianos'])
disp(['X',' = ',num2str(X),' m'])
disp(['Ka',' = ',num2str(Ka)])
disp(['Kx',' = ',num2str(Kx),' m'])
disp(['La',' = ',num2str(La)])
disp(['Lx',' = ',num2str(Lx),' m'])
disp(' ')

% % 2o. passo - giro completo em q - posição
% qv=linspace(0,2*pi,73);
% 
% % cálculo e visualização das variáveis secundárias
% Av=asin((R*sin(qv)-C)/L);
% Xv=(R*cos(qv))+(L*cos(Av));
% figure(1)
% subplot(1,2,1)
% plot(qv*180/pi,Av*180/pi,'-r')
% axis([0 360 -25 25])
% grid on
% xlabel('Ângulo da manivela (graus)')
% ylabel('Ângulo da biela (graus)')
% title('Movimento da biela')
% subplot(1,2,2)
% plot(qv*180/pi,Xv*1000,'-g')
% axis([0 360 160 360])
% grid on
% xlabel('Ângulo da manivela (graus)')
% ylabel('Deslocamento do pistão (mm)')
% title('Movimento do pistão')

% % 3o. passo - giro completo em q - velocidade
% 
% % cálculo dos coeficientes de velocidade
% Kav=(R*cos(qv))./(L*cos(Av));
% Kxv=-(C+Xv.*tan(Av));
% figure(2)
% subplot(1,2,1)
% plot(qv,Kav,'-r')
% axis([0 2*pi -0.4 0.4])
% grid on
% xlabel('Ângulo da manivela (rad)')
% ylabel('Coef. de velocidade de A (rad/rad)')
% title('Movimento da biela')
% subplot(1,2,2)
% plot(qv,Kxv*1000,'-g')
% axis([0 2*pi -100 100])
% grid on
% xlabel('Ângulo da manivela (rad)')
% ylabel('Coef. de velocidade de X (mm/rad)')
% title('Movimento do pistão')

% % 4o. passo - giro completo em q - aceleração
% 
% % cálculo das derivadas dos coeficientes de velocidade
% Lav=((-R*sin(qv))./(L*cos(Av)))+((Kav.^2).*tan(Av));
% Lxv=-(Kxv.*tan(Av))+((Xv.*Kav)./(cos(Av).^2));
% figure(3)
% subplot(1,2,1)
% plot(qv,Lav,'-r')
% axis([0 2*pi -0.4 0.4])
% grid on
% xlabel('Ângulo da manivela (rad)')
% ylabel('Deriv de coef. de vel. de A (rad/rad^2)')
% title('Movimento da biela')
% subplot(1,2,2)
% plot(qv,Lxv*1000,'-g')
% axis([0 2*pi -60 140])
% grid on
% xlabel('Ângulo da manivela (rad)')
% ylabel('Deriv. de coef. de vel. de X (mm/rad^2)')
% title('Movimento do pistão')
