科学技术当中的微分方程模型建模总结.pptVIP

①　求解问题（１），Matlab命令为 用Matlab软件求解。 求解出函数y(t)后，求出圆桶下落到水深为91米海底的时间，即求满足y(t)=91的时间t1,然后求出在时间t1的速度即可以得到问题的解答。这里的速度关系由直接求函数y(t)对t的导数得到。 syms t; syms y; y=dsolve(m*D2y+k*Dy-w+b=0,y(0)=0,Dy(0)=0) 运行结果为 y =-m/k^2*exp(-k/m*t)*(-w+b)-(-w+b)/k*t+m*(-w+b)/k^2 即 运行结果如图 b=1025.94*0.208; k=0.12; w=239.456; m=w/9.8; % 这里重力加速度取为9.8米/秒2 t=0?:0.1?:20?; y=-m/k^2*exp(-k/m*t)*(-w+b)-(-w+b)/k*t+m*(-w+b)/k^2-91?; plot(t,y) ②　求解方程y(t)=91，先画出图形，命令为 可见方程的根在t=13附近， 为求解方程，建立M文件： function y=my0(t) b=1025.94*0.208; k=0.12; w=239.456; m=w/9.8; y=-m/k^2*exp(-k/m*t)*(-w+b)-(-w+b)/k*t+m*(-w+b)/k^2-91?;