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

답변 (1개)

Star Strider
Star Strider 2020년 12월 7일

0 개 추천

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

MN
MN 2020년 12월 8일
편집: MN 2020년 12월 8일
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 2020년 12월 8일
Probably.
Give it a shot!

댓글을 달려면 로그인하십시오.

카테고리

도움말 센터File Exchange에서 Call Python from MATLAB에 대해 자세히 알아보기

질문:

MN
2020년 12월 7일

댓글:

2020년 12월 8일

Community Treasure Hunt

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

Start Hunting!

Translated by