differential equation 2DOF

1 次查看（过去 30 天）

alsgud qor 2020-3-20

评论： Birdman 2020-3-23

------------------------------------------------------

function homework2 = main2(t,X1,X2)

global m1 c1 k1 m2 c2 k2

M=[m1 0;0 m2];

C=[c1+c2 -c2; -c2 c2];

K=[k1+k2 -k2; -k2 k2];

F=[3*cos(t);10*cos(3*t)];

X=[X1;X2];

homework2=M*X(3)+C*X(2)+K*X(1)-F;

end

-------------------------------------------------------

tspan= 0:0.1:50 ;

X0 = [0 0;0 0] ;

[t, X] = ode45(@main2, tspan, X0);

plot(X1,X2)

I do it here and it 's ain't working....

Birdman 2020-3-20

You may try Symbolic Toolbox to solve the problem. Run the attached script to see the result.

alsgud qor 2020-3-21

um sorry.because of time lag, I just recieved it. let me ask some Q

dydt = zeros(4,1); what does it mean?
dydt(2) = (-k1/m1)*y(1)+(-c1/m1)*y(2)+F1; this part should be changed into 'F1/m1'? same with dydt(4)

Birdman 2020-3-23

MATLAB Programming

在 Help Center 和 File Exchange 中查找有关 Programming 的更多信息

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

Translated by