How to make a 10mm impulse input

hello
i have this code, i want to know how can change the impulse to 10mm magnitude?
Thanks
Ur_M=53; m1=Ur_M/12.33; m2=Ur_M/12.33; m3=Ur_M/6.16; m4=Ur_M/1.48; k1=6; k2=6; k3=10; k4=10; k5=18; c1=0.3; c2=0.65; c4=1.9; ms=0.3;
ks=403; cs=2150;
ub=1;
KG = [100 350 850] ;
u1_disp_imp = cell(3,1) ;
for i = 1:length(KG)
kg = KG(i) ;
A= [0 1 0 0 0 0 0 0 0 0; (-kg-ks)/ms -cs/ms ks/ms cs/ms 0 0 0 0 0 0;
0 0 0 1 0 0 0 0 0 0; ks/m1 cs/m1 (-ks-k1-k2)/m1 (-cs-c1-c2)/m1 k2/m1 c2/m1 k1/m1 c1/m1 0 0;
0 0 0 0 0 1 0 0 0 0; 0 0 k2/m2 c2/m2 (-k3-k2)/m2 -c2/m2 k3/m2 0 0 0; 0 0 0 0 0 0 0 1 0 0;
0 0 k1/m3 c1/m3 k3/m3 0 (-k1-k3-k4-k5)/m3 (-c1-c4)/m3 (k4+k5)/m3 c4/m3;
0 0 0 0 0 0 0 0 0 1; 0 0 0 0 0 0 (k4+k5)/m4 c4/m4 (-k4-k5)/m4 -c4/m4];
B= [0; kg/ms; 0; 0; 0; 0; 0; 0; 0; 0];
C= [1 0 0 0 0 0 0 0 0 0; 0 0 1 0 0 0 0 0 0 0; 0 0 0 0 1 0 0 0 0 0; 0 0 0 0 0 0 1 0 0 0; 0 0 0 0 0 0 0 0 1 0];
D= [0; 0; 0; 0; 0];
My_sys=ss(A,B,C,D)
[U,T,X]= impulse(My_sys)
u1_disp_imp{i} = X(:,3) ;
end
figure
hold on
for i = 1:3
plot(u1_disp_imp{i})
end
legend('disp due to impulse input for soft soil','disp due to impulse input for hard soil','disp due to impulse input for very hard soil')

Respuestas (1)

Star Strider
Star Strider el 7 de Dic. de 2020

0 votos

The Dirac function is defined as having unit area, infinite amplitude, and infinitesimal width. None of those can be changed.
The step and stepplot functions allow changes in the amplitude of the step input with the StepAmplitude name-value pair in stepDataOptions.

2 comentarios

MN
MN el 8 de Dic. de 2020
Editada: MN el 8 de Dic. de 2020
is there a way to make this input using the lsim? I need to get my plots in the impulse format, so i can identify the amplitude and the time for each peak.
Star Strider
Star Strider el 8 de Dic. de 2020
Probably.
Give it a shot!

Iniciar sesión para comentar.

Categorías

Más información sobre Graphics en Centro de ayuda y File Exchange.

Etiquetas

Preguntada:

MN
el 7 de Dic. de 2020

Comentada:

el 8 de Dic. de 2020

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by