Programação

clear allclose allclc

% Parâmetros do sistemaparam.m = 1; % massaparam.c = 0.1; % coeficiente de amortecimentoparam.k = 1; % rigidezparam.A = 1; % Amplitude do deslocamento da base

wvector = logspace(-1,1,100);

% Parâmetros do tempo de integraçãodeltat = 0.01; % Passo de integração temporaltzero = 0; % Tempo de início da integraçãotend = 200; % Tempo de fim de integraçãotspan = [tzero:deltat:tend]; % Vetor de tempo de integração

% Condições iniciaisy0(1) = 0; % Condição inicial do deslocamento - Variável de estado y1y0(2) = 0; % Condição inicial da velocidade - Variável de estado y2

for c1 = 1:1:length(wvector)

param.w = wvector(c1); % Frequência de excitação externa

% Seção principal - Integração[T Y] = integrador(tspan, y0, param); % Chama a função integrador% T,Y parâmetros de saída || tspan, y0 e estrutura param, parâmetros % de entrada

a = 0.8*(length(tspan)-1);

end

% Resultados% gerar histórico no tempo do deslocamento% figure(1);% plot(T,Y(:,1),T,Y(:,2))% xlabel('time [s]');% ylabel('displacement [m], velocity [m/s]');% grid

% definir índice de início do regime permanente em tspan% a = 0.8*(length(tspan)-1);

% Plano de fase% figure(2);% plot(Y(1:a,1),Y(1:a,2),'b', 'linewidth', 0.5)

% hold on;% plot(Y(a:end,1),Y(a:end,2), 'g', 'linewidth', 0.5)% hold off;% xlabel('displacement [m]');% ylabel('velocity [m]');% grid

% figura em 3D - as duas variáveis de estado no tempo% deslocamento, velocidade, tempo %figure(3)%plot3(Y(:,1),Y(:,2),T(:,1),'b',Y(a:end,1),Y(a:end,2),T(a:end,1),'r')%xlabel('displacement [m]')%ylabel('velocity [m/s]')%zlabel('time [s]')%grid

%% Seção de energia% Ec0 = 0.5.*param.m*Y(1,2).^2;% Ep0 = 0.5.*param.k*Y(1,1).^2;% Et0 = Ec0+Ep0;% % Ec = 0.5.*param.m*Y(:,2).^2;% Ep = 0.5.*param.k*Y(:,1).^2;% Ed = Et0 - Ec - Ep;% % figure(4);% plot(T,Ep,T,Ec,T,Ed,T,(Ec+Ep),T,(Ed+Ec+Ep))% xlabel('time [s]');% ylabel('Energy [joule]');% grid% legend('Energia potencial','Energia cinética','Energia dissipada', 'Energia Total','Soma das energias')