24小时热门版块排行榜

返回列表

wpynwpu

木虫 (正式写手)

应助: 4 (幼儿园)
金币: 1469.1
红花: 1
帖子: 550
在线: 262.5小时
虫号: 438036
注册: 2007-09-13
性别: GG
专业: 金属材料的凝固与结晶学

[求助] 采用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
请大家帮帮忙分析一下错误原因

回复此楼

» 猜你喜欢

什么是人一生最重要的？已经有6人回复
为什么中国大学工科教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？已经有11人回复
网上报道青年教师午睡中猝死、熬夜猝死的越来越多，主要哪些原因引起的？已经有9人回复
【博士招生】太原理工大学2026化工博士已经有5人回复
280求调剂已经有3人回复
面上可以超过30页吧？已经有11人回复
版面费该交吗已经有15人回复
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低已经有18人回复

» 本主题相关价值贴推荐，对您同样有帮助:

高手进来用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人回复

1楼 2012-05-26 10:17:58

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

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

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);