Я хотел бы изменить параметр (o) в моем ODE45 для механизма сцепления. Муфта включается, когда скорость и ускорение находятся в одном направлении (ускорение / уровень> 0), и отключается, когда (ускорение / уровень <0). Я не уверен, как реализовать l oop в файле fuction. </p>
*if dy(1)*dy(2)<0 %%% need the velocity and floor acceleration in the same direction
o=0
else
o=0.5;
end
dy(1)= y(2);
dy(2)=(1/u)*(-1*(1-v)*uu-2*e1*w1*y(2)+2*v*e2*w2*y(4)-w1*w1*y(1)+v*w2*w2*y(3));
dy(3) = y(4);
dy(4) = (1/u)*((-o*uu)+2*e1*w1*y(2)-2*(1+o)*e2*w2*y(4)+w1*w1*y(1)-(1+o)*w2*w2*y(3));*