| 查看: 1336 | 回复: 3 | ||
wpynwpu木虫 (正式写手)
|
[求助]
采用matlab求解微分方程遇到问题，请帮忙分析一下错误-重谢
|
|
问题是：液滴运动的微分方程；液滴初始速度10m/s，求解时间时0.15s，采用龙格-库塔法求解，程序如下： function y=fly1(t,U) a=1.1616; %氮气密度 b=1.873e-5; %氮气黏度 d=3e-4; %熔滴直径 c=2385; %熔滴密度 g=9.8; %重力加速度 f=d*a*U/b %Re数 h=0.28+6/f.^0.5+21/f %拖曳系数 y=((1-a/c)*g-3*(0.28+6/f.^0.5+21/f)*a*U.^2/(4*d*c))/(1+0.5*a/c) options=odeset('outputfcn','odeplot') [t,U]=ode45('fly1',[0:0.001:0.15],[0;10],options) U1=U(:,1) plot(t,U1) title('Droplet Velocity VS Time') xlabel('Time/s') ylabel('Droplet Velocity/m/s') 出现的错误是： ??? Maximum recursion limit of 2000 reached. Use set(0,'RecursionLimit',N) to change the limit. Be aware that exceeding your available stack space can crash MATLAB and/or your computer. Error in ==> odeset 请大家帮帮忙分析一下错误原因 |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有6人回复
-
明天应该可查了！？
已经有5人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有7人回复
-
没有任何消息-是不是就凉了
已经有8人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
人气不行了
已经有10人回复
-
我面上完蛋了
已经有10人回复
-
能否退出参与的面上项目解除限项
已经有24人回复
-
filecode，4个jtjc了
已经有17人回复
-
2026国自然函评费到账
已经有19人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 高手进来用matlab求解常微分方程，十分感谢！
已经有17人回复
- matlab微分方程组参数拟合，以周为单位求解，汇总后以年为单位进行数值比较
已经有4人回复
- matlab 线性常微分方程求解
已经有5人回复
- matlab解微分方程
已经有10人回复
- Matlab求解偏微分方程组
已经有22人回复
- matlab数值求解边界条件微分方程组
已经有7人回复
- 高金求助matlab解微分方程组
已经有12人回复
- matlab解微分方程组
已经有15人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
wpynwpu: 金币+50, ★★★★★最佳答案, 真是高手，谢谢 2012-05-26 11:04:43
臭水沟: 金币+3, 谢谢交流~~ 2012-05-26 11:06:04
感谢参与，应助指数 +1
wpynwpu: 金币+50, ★★★★★最佳答案, 真是高手，谢谢 2012-05-26 11:04:43
臭水沟: 金币+3, 谢谢交流~~ 2012-05-26 11:06:04
|
function exercise11 options=odeset('outputfcn','odeplot'); [t,U]=ode45(@fly1,[0:0.001:0.15],10,options); U1=U(:,1); plot(t,U1) title('Droplet Velocity VS Time') xlabel('Time/s') ylabel('Droplet Velocity/m/s') function y=fly1(t,U) a=1.1616; %氮气密度 b=1.873e-5; %氮气黏度 d=3e-4; %熔滴直径 c=2385; %熔滴密度 g=9.8; %重力加速度 f=d*a*U/b; %Re数 h=0.28+6/f.^0.5+21/f; %拖曳系数 y=((1-a/c)*g-3*(0.28+6/f.^0.5+21/f)*a*U.^2/(4*d*c))/(1+0.5*a/c); |
2楼2012-05-26 10:51:19
wpynwpu
木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 1469.1
- 红花: 1
- 帖子: 550
- 在线: 262.5小时
- 虫号: 438036
- 注册: 2007-09-13
- 性别: GG
- 专业: 金属材料的凝固与结晶学
3楼2012-05-26 11:04:02
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
4楼2012-05-26 12:31:24
20