实验5 常微分程的数值解.docVIP

实验5 常微分方程的数值解 概要：将装满放射性废物的圆桶扔到水深300ft的海底，圆桶体积55gal，装满废料的桶重为527.436lbf，在海中浮力为470.327lbf。此外，下沉时受到的阻力与速度成正比，比例系数为0.08lbf/s。实验发现当圆桶速度超过40ft/s时，就会因与海底冲撞而破裂。 要求：（1）建立解决上述问题的微分方程模型（2）用数值和解析两种方法求解微分方程，并回答谁赢得了官司。 模型建立 由牛顿第二定律可列出圆桶下沉速度的微分方程： 其中G为圆桶重量，F为浮力，b为下沉阻力与速度关系的比例系数。换算到国际单位制，dept=300*0.3048=91.4400 海深（m） ve=40*0.3048=12.1920 速度极限（超过就会使圆筒碰撞破裂）(m/s) G=527.436*0.4536*9.8＝2344.6 圆筒重量（N） F=470.327*0.4536*9.8＝2090.7 浮力(N) m=527.436*0.4536＝239.24 圆筒质量（kg） b=0.08*0.4536*9.8/0.3048＝1.1667 比例系数(Ns/m) 模型求解 一．求数值解 Matlab程序如下： s