zhendong_xt_matlab_响应时间_频响_nan_ZhenDong_

function dy = func_ode45(t,y) dy = zeros(6,1); m1 = 1; m2 = 2; m3 = 3; k1 = 1; k2 = 2; k3 = 3 c1 = 1; c2 = 2; c3 = 3; dy(1) = y(2); dy(2) = k1/m1*y(3) - k1/m1*y(1) + c1/m1*y(4) - c1/m1*y(2); dy(3) = y(4); dy(4) = k2/m2*y(5) + k1/m2*y(1) -(k1+k2)/m2*y(3) + c1/m2*y(2) - c1/m2*y(4); dy(5) = y(6); dy(6) = k2/m3*y(3) - (k2+k3)/m3*y(5) - c3/m3*y(6); end