视频

拓展学习

（1）求最大安全体重

for m=100:-0.5:20

    [t,x,y_w]=sim('bengji',0:0.01:100);

    if min(y_w)>1.5

        break;

    end

end

disp(['最大安全体重是',num2str(m)])

dis=min(y_w);

disp(['最小的安全距离是',num2str(dis)])

plot(t,y_w)

grid on

（2）求最小弹性系数

m=65; 

for k=10:0.1:50

    [t,x,y_w]=sim('bengji',0:0.01:100);

    if min(y_w)>1.5

        break;

    end 

end

disp(['最小弹性系数k是',num2str(k)])

dis=min(y_w);

disp(['最小的安全距离是',num2str(dis)])

plot(t,y_w)

grid on