24小时热门版块排行榜    

查看: 1146  |  回复: 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
请大家帮帮忙分析一下错误原因
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与,应助指数 +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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wpynwpu

木虫 (正式写手)

你真是高手呀,大侠,我的程序错在什么地方了?
3楼2012-05-26 11:04:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师

引用回帖:
3楼: Originally posted by wpynwpu at 2012-05-26 11:04:02
你真是高手呀,大侠,我的程序错在什么地方了?

主程序、子程序不分;
ode45函数不会用;
4楼2012-05-26 12:31:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wpynwpu 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 354求调剂 +3 Tyoumou 2026-03-18 5/250 2026-03-18 16:41 by Tyoumou
[考研] 材料与化工一志愿南昌大学327求调剂推荐 +8 Ncdx123456 2026-03-13 9/450 2026-03-18 14:40 by haxia
[考研] 一志愿武理材料305分求调剂 +3 想上岸的鲤鱼 2026-03-18 4/200 2026-03-18 14:33 by laoshidan
[考研] 材料专硕306英一数二 +10 z1z2z3879 2026-03-16 13/650 2026-03-18 14:20 by 007_lilei
[考研] 材料专硕274一志愿陕西师范大学求调剂 +6 薛云鹏 2026-03-13 6/300 2026-03-18 14:14 by 脱颖而出
[考研] 0703化学求调剂 总分331 +3 ZY-05 2026-03-13 3/150 2026-03-18 10:58 by macy2011
[考研] 考研化学学硕调剂,一志愿985 +4 张vvvv 2026-03-15 6/300 2026-03-17 17:15 by ruiyingmiao
[考研] 302求调剂 +9 负心者当诛 2026-03-11 9/450 2026-03-17 17:13 by ruiyingmiao
[考研] 材料专硕326求调剂 +6 墨煜姒莘 2026-03-15 7/350 2026-03-17 17:10 by ruiyingmiao
[考研] 有没有道铁/土木的想调剂南林,给自己招师弟中~ +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
[考研] 267一志愿南京工业大学0817化工求调剂 +6 SUICHILD 2026-03-12 6/300 2026-03-17 09:24 by 雾散后相遇lc
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 0703化学调剂 290分有科研经历,论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 5/250 2026-03-14 16:30 by JourneyLucky
[考研] 328求调剂 +3 5201314Lsy! 2026-03-13 6/300 2026-03-14 15:31 by hyswxzs
[基金申请] 有必要更换申报口吗 20+3 fannyamoy 2026-03-11 3/150 2026-03-14 00:52 by zhanghaozhu
[考研] 材料工程调剂 +9 咪咪空空 2026-03-12 9/450 2026-03-13 22:05 by 星空星月
[考研] 085600材料与化工 309分请求调剂 +7 dtdxzxx 2026-03-12 8/400 2026-03-13 14:43 by jxchenghu
[考研] 材料301分求调剂 +5 Liyouyumairs 2026-03-12 5/250 2026-03-13 14:42 by JourneyLucky
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生,2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
信息提示
请填处理意见