How can I solve this matrix for y?

Question

Keith Grey 2020-6-15

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/548364-how-can-i-solve-this-matrix-for-y

编辑： Keith Grey 2020-6-15

M = [1 0 0 0; 0 1 0 0; 0 0 20 0; 0 0 0 10]; % Mass matrix
y = [x1; x2; v1; v2]; % Displacement & Velocity Matrix
B = [0 0 -1 0; 0 0 0 -1; -2 2 -2 2; 1 -1 1 -1]; % Spring Matrix
C = [0; 0; 0; -1]; % Force Matrix

I'm following an example:

Example

At the bottom of the page, Frequency vs Amplitude is achieved. That's what I'm trying to do.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Ameer Hamza 2020-6-15

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/548364-how-can-i-solve-this-matrix-for-y#answer_451335

See ode45()

tspan = [0; 10];
IC = [1 1 0 0]; % initial condition
[t, Y] = ode45(@odeFun, tspan, IC);
plot(t, Y);
legend({'x_1', 'x_2', 'v_1', 'v_2'});
function dYdt = odeFun(t, Y)
    M = [1 0 0 0; 0 1 0 0; 0 0 20 0; 0 0 0 10]; % Mass matrix
    B = [0 0 -1 0; 0 0 0 -1; -2 2 -2 2; 1 -1 1 -1]; % Spring Matrix
    C = [0; 0; 0; -1]; % Force Matrix
    
    dYdt = M\(C - B*Y);
end