TP2 Matlab-combined

K=1.8e4;%raideur;
M=375;%maase;
C=1660;%amortissement;
w=15;%pulsation forcée;
F0=100;%amplitude de la force,
X0=5e-2;%solution;
V0=0;%vitesse initiale;
w0=sqrt(K/M);%pulsation propre;
B=(w/w0);
Z=C/(2*M*w0);
w1=w0*sqrt(1-(Z^2));
Te=1e-3;
T=0:1e-3:5;
F=F0*sin(w*T);
if Z<1
A1=X0;
B1=(Z*w0*X0)/w1;
for i=1:length(T)
xh(i)=exp(-Z*w0*T(i))*(A1*cos(w1*T(i))+B1*sin(w1*T(i)));
end
elseif Z==1
A1=X0;
B1=Z*w0*X0;
for i=1:length(T)
xh(i)=exp(-z*w0*T(i))*(A1+(B1*T(i)));
end
else Z>=1;
A1=X0;
B1=(Z*w0*X0)/w1;
for i=1:length(T)
xh(i)=exp(-Z*w0*T(i))*(A1*cosh(w1*T(i))+B1*sinh(w1*T(i))) ;
end
end
fi=atan(2*Z*B)/(1-B^2);
xp=F0*sin(w*T-fi)/(K*sqrt(((1-B^2)^2)+(2*Z*B)^2));
x=xh+xp;
plot(T,x)
figure %1
plot(T,xp)
figure %2
plot(T,x,T,xh,T,xp)
figure %3
A=[0,1;-K/M,-C/M];
y=[5e-2;0];
B=zeros(2,length(T));
for i=1:length(T)
B(2,i)=F(i)/M;
end
for i=2:length(T)
y(:,i)=y(:,i-1)+Te*(A*y(:,i-1)+B(:,i));
end
plot(T,y(1,:))
1

figure %4
plot(T,x,T,y(1,:))
figure %5
2

3

4

5

6

Published with MATLAB® R2022b
7

Maoag bmloei
l eWo e "(A,cin (w,t+ B, co(ut))
At o e (A, cone +b/Gm®)
A4
aNec cnO:1
m
O= o
lt)=- E w,e AACnWE) 5, m (w,b)e (-AASn()
+Bcou (waE}
At-o Vo = 0 = -'CW, An64
Atzon
- Ewst
ABt)
1o = A
H = - C e(A-6zt)-e 6
AE-o
Vo = - Ew A +Bz
, Az
Ag Xo