I have a vibration problem which asks to find and plot the t

function tt, yy model02 x0 0.01 0.01 100 3.5 [t, x] ode45 f ts,x0) plot (t,x(: 1 t,x 2) end M 10 0 0 5 K [1000 400 400 1200 B 1 300 B2-800 2 50 A zeros (2) eye (2) inv (M) *K zeros (2) fl inv (M) *B1; 2 inv (M) *B2; v1 A*x+ 0 0 fl] cos (w1*t); end

Solution

solution:

for given forced undamped system equation of motion in matrix form can be writtena s

[K-M*w^2]X=F

hecne programme in matlab as follows

functionX=forced_vibration(K,M,F,omega)

K=[1000 -400;-400 1200];

M=[10 0;0 5];

F=[300;800];

omega=[5;50];

for i=1:2

omega=omega i

Xi=(K-M*(omega i)^2)\\F

fprintf(\'value of amplitude for first system\');

end

end