实验4 常微分方程数值解.docxVIP

实验 4 常微分方程数值解 实验目的 ： 1、 用 MATLAB软件掌握求微分方程数值解的方法，并对结果作初步分析； 2、 通过实例学习用微分方程模型解决简化的实际问题。 【题目 1】放射性废物的处理：有一段时间，美国原子能委员会（现为核管理 委员会）处理浓缩放射性废物时， 把它们装入密封性能很好的圆桶中， 然后扔到水深 300ft 的大海中。这种做法是否会造成放射性污染，自然引起生态学家及社 会各界的关注。 原子能委员会一再保证圆桶非常坚固， 绝不会破漏， 这种做法是绝对安全的。然而一些工程师们却对此表示怀疑， 认为圆桶在海底相撞时有可能发生破裂。于是双方展开了一场笔墨官司。 究竟谁的意见正确呢？原子能委员会使用的是 55gal 的圆桶，装满放射性废物时的圆桶重量为 527.436lbf ，在海水中受到的浮力为 470.327lbf 。此外，下沉时圆桶还要受到水的阻力， 阻力与下沉速度成正比， 工程师们做了大量实验， 测得其比例系数为 0.08lbf s/ft 。同时，大量破坏性实验发现当圆桶速度超过 40ft/s 时，就会因为与海底冲撞而发生破裂。 1） 建立解决上述问题的微分方程模型。 2） 用数值和解析两种方法求解微分方程，并回答谁赢了这场官司。 【问题分析】 本题应该求出圆桶下沉到海底时的速度，比较其与 40ft/s 的大 小，来决定谁最后赢了这场官司。 已知水深