Question

Would appreciate some help with this, please post the code as well as some plots, thank you.

FOUR- Matlab Assume the following equations of motion for a system: mLO cos 0 - m,L0 sin0+ kx = F, sin ot (mm m,L m2L cos0 + m,gL sin 0 L cos 0F sin ot For mm2 1 kg, k=1 N/m, L= 1m, o = 1 rad/s, fo-1 N, g=10m/s2, Solve the EOM using Matlab. Play with the initial conditions - try several of them. Plot q versus / and x versus t.

Answer 1

clc
clear all
close all
format long
F0=1;
w=1;
k=1;
m1=1;
m2=1;
g=10;
L=1;
f=@(t,y) [y(2); -(g*m2^2*sin(y(1)) - k*m2*cos(y(1))*y(3) + g*m1*m2*sin(y(1)) - F0*m1*cos(y(1))*sin(t*w) + L*m2^2*cos(y(1))*sin(y(1))*y(2)^2)/(L*m2*(m1 + m2 - m2*cos(y(1))^2));y(4);(F0*sin(t*w) - k*y(3) - F0*cos(y(1))^2*sin(t*w) + L*m2*sin(y(1))*y(2)^2 + g*m2*cos(y(1))*sin(y(1)))/(m1 + m2 - m2*cos(y(1))^2)];
[T,Y]=ode45(f,[0,10],[1;1;1;1]);
plot(T,Y(:,1));
figure;
plot(T,Y(:,3));

A MATLAB R2018 HOME PLOTS APPS EDITOR PUBLISH VIEW B E N IFTY 50 SENSEX 12,021.30 40,757.42 81.20 287.72 0.6996 0.72% Figure 1 - O X Run S Advar New Figure 2 File Edit 8 . Open View File Edit View Insert Tools Desktop Window Help SHOW @moko. @ IE 0 Insert Window Tools " Desktop A Help E O Kax ATLAB g1.m x + 0.5 Current Folder Name 0409 1029 1033 1036 1040 1045 1046 1049 1055 3082 Advanced af-ZA am-et AppLocke appraiser Details + L*m2^2 (v (1 ,1]); X Workspace -1.5 0 1 2 3 4 5 6 7 8 9 10 Name 0 1 2 3 4 5 6 7 8 9 10 FO Hk HL m1 m2 HT 245x1 double script . Ln 11 Col 5 ENG 10:27 AM IN 11/20/19 W Type here to search