La respuesta es la suma de varios pulsos que son equivalentes a la carga externa aplicada.

Ejemplo

Se determina la respuesta para la superposición de impulsos escalonados utilizando los datos del ejercicio base

Calculo de la respuesta ante impulsos expresados utilizando superposicion
clear all
clc

Datos
m=14700 %[kg];
k=17031767.414 %[N/m];
xi=0.10;
Po=0.20*m*9.81 % [N];
to=0.2; %[s]
tf=1.8; %[s]

wn=(k/m)^0.5;
Tn=2*pi()/wn

step=0.001;

t1=[0:step:to to+step:step:tf];
t2=[0:step:tf-to];

Calculo de la respuesta sin amortiguamiento
% Resolucion de la ecuacion para vibracion forzada para una carga
% escalonada (con y sin amortiguamiento) - Step function
%   Los datos requeridos son: [xt_a,xt_na] = SDOF_StepF(t,m,k,xi,Po)
%   [t] → intervalo de tiempo
%   m → masa del sistema
%   k → rigidez del sistema
%   xt_a → respuesta de desplazamiento en el tiempo sin amortiguamiento
%   xt_sa → respuesta de sesplazamiento en el tiempo con amortiguamiento

step=0.001;

t1=[0:step:to to+step:step:tf];
t2=[0:step:tf-to];

[y,xt_sa_1] = SDOF_StepF(t1,m,k,xi,Po);
[z,xt_sa_2] = SDOF_StepF(t2,m,k,xi,-Po);

% Suma de los dos vectores
xt_temp=zeros(1,length(0:step:to)-1);
xt_sa_2 = [xt_temp xt_sa_2];
t2=[xt_temp t2+to];

figure
plot(t1,xt_sa_1,'k')
hold on
plot(t2,xt_sa_2,'b')
hold on

figure
plot(t1,xt_sa_1+xt_sa_2,'k')

Calculo de la respuesta con amortiguamiento
% Resolucion de la ecuacion para vibracion forzada para una carga
% escalonada (con y sin amortiguamiento) - Step function
%   Los datos requeridos son: [xt_a,xt_na] = SDOF_StepF(t,m,k,xi,Po)
%   [t] → intervalo de tiempo
%   m → masa del sistema
%   k → rigidez del sistema
%   xt_a → respuesta de desplazamiento en el tiempo sin amortiguamiento
%   xt_sa → respuesta de sesplazamiento en el tiempo con amortiguamiento

step=0.001;

t1=[0:step:to to+step:step:tf];
t2=[0:step:tf-to];

[xt_a_1,xt_na_1] = SDOF_StepF(t1,m,k,xi,Po);
[xt_a_2,xt_na_2] = SDOF_StepF(t2,m,k,xi,-Po);

% Suma de los dos vectores
xt_temp=zeros(1,length(0:step:to)-1);
xt_a_2 = [xt_temp xt_a_2];
t2=[xt_temp t2+to];

figure
plot(t1,xt_a_1,'k')
hold on
plot(t2,xt_a_2,'b')
hold on

figure
plot(t1,xt_a_1+xt_a_2,'k')