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')

Réponses (1)

Star Strider
Star Strider le 7 Déc 2020

0 votes

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 commentaires

MN
MN le 8 Déc 2020
Modifié(e) : MN le 8 Déc 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 le 8 Déc 2020
Probably.
Give it a shot!

Connectez-vous pour commenter.

Catégories

En savoir plus sur Graphics dans Centre d'aide et File Exchange

Question posée :

MN
le 7 Déc 2020

Commenté :

le 8 Déc 2020

Community Treasure Hunt

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

Start Hunting!

Translated by