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

清华大学《数学实验》实验报告  实验四 常微分方程数值解 姓名：芦琛璘 班级：化 33 学号：2013011934 实验目的： 1、 练习数值微分的计算；  2、 掌握用 MATLAB 软件求微分方程初值问题数值解的方法；  3、 通过实例学习用微分方程模型解决简化的实际问题；  4 、 了解欧拉方法和龙格‐库塔方法的基本思想和计算公式，及稳定性等概念。    实验内容： 【问题1】 放射性废物的处理：将反射性废物装入桶中，仍到水深为300ft 的大海中，桶重55gal，装满放射性物质后桶重527.436lbf，浮力470.327lbf。 阻力系数0.08lbf s/ft，当速度超过40ft/s时就会与海底碰撞而破裂。 ⑴ 建立解决上述问题的微分方程模型 ⑵ 用数值和解析的方法求解上述的微分方程模型 【模型建立】建立直角坐标系tOh，t=0时放射桶的位置在O(0,0)，质量 m=527.436*0.4536 kg；重力G=527.436*0.4536*9.8 N；浮力F=470.327*0.453 6*9.8 N；f为阻力，f=k*v；k为阻力系数，k=0.08l*0.4536*9.8s/0.3048m。 由牛顿第二定理知ma=G-F-f，所以列出微分方程:   初始条件为 0 0，x 0 0； 【模型求解】 （1）数值解 fangshe.m: function dx=fangshe(~,x) G=527.436*0.4536*9.8;%桶重 F=470.327*0.4536*9.8%浮力 m=527.436*0.4536;%质量 k=0.081*0.4536*9.8/0.3048;%阻力系数 dx=[x(2);G/m-F/m-k*x(2)/m];%微分方程 end fangsheode.m 化学工程系 化 33 芦琛璘  清华大学《数学实验》实验报告  function fangsheode ts=0:0.01:1000;%粗略估计到1000秒 x0=[0,0];%初值 opt=odeset(reltol,1e-3,abstol,1e-6); %相对误差1e-6，绝对误差1e-9  %使用5级4阶龙格—库塔公式计算 [t,x]=ode45(@fangshe,ts,x0); result=[t,x]% %输出计算结果 figure,plot(t,x(:,2)) %绘出速度-时间图 title(speed-t) grid; xlabel(speed),ylabel(t); figure,plot(t,x(:,1)) %绘出深度-时间图 title(deepth-t), grid  xlabel(t),ylabel(deepth); figure,plot(x(1:1500,1:1),x(1:1500,2:2),k)%绘出速度时间图 title(speed-deepth), grid% xlabel(deepth),ylabel(speed); hold on plot([91.440,91.440],[0,15], g); %绘出海底深度所对应的曲线 plot([0,120],[40*0.3048,40*0.3048], r); %绘出到达海底时所能承受的最大速  %找出盛有放射物的桶到达海底时的速度 度 [deepth,speed]=ginput end 【结果如下】