版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

>论坛更新日志 (5558)
>考研 (1102)
>导师招生 (1010)
>文献求助 (410)
>虫友互识 (264)
>基金申请 (202)
>休闲灌水 (180)
>论文投稿 (140)
>硕博家园 (123)
>考博 (120)
>招聘信息布告栏 (96)
>博后之家 (67)
>教师之家 (65)
>绿色求助(高悬赏) (40)
>公派出国 (37)
>留学生活 (32)

返回列表

iphysunba

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 5.5
帖子: 5
在线: 2.2小时
虫号: 1830507

[交流] 请大家帮忙分析一下我这个matlab解微分方程初值的例子,到底问题出在哪里

请大家帮忙分析一下我这个matlab解微分方程初值的例子,到底问题出在哪里：
m. file:
function free
options=odeset('Jacobian', @FJac)
[t,xx]=ode15s(@serr,[0 1000000],[0.01 0.8671],options);
save t t -ascii;
save xx xx -ascii;
function dx=serr(t,x)
dx=zeros(2,1);
dx(1)=exp(-(1/x(1)))*((cosh(x(2))-1)/120/x(1)-1/3);
dx(2)=4.8e-4-240*exp(-(1/x(1)))*sinh(x(2));
end
function dfdx=FJac(t,x)
dfdx=[(exp(-1/x(1))/x(1)^2*((cosh(x(2))-1)*(1/x(1)-1)/120-1/3))  (exp(-1/x(1))/x(1)/120*sinh(x(2)));
   (-exp(-1/x(1))/x(1)^2*240*sinh(x(2)))  (-240*exp(-1/x(1))*cosh(x(2)))]
end
end
在matlab 运行 free，结果出现下列提示：
Warning: Failure at t=1.729686e+005.  Unable to meet integration tolerances without
reducing the step size below the smallest value allowed (4.656613e-010) at time t.
> In ode15s at 753
  In free at 3

回复此楼

» 猜你喜欢

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

高手进来用matlab求解常微分方程，十分感谢！已经有17人回复
matlab微分方程组参数拟合，以周为单位求解，汇总后以年为单位进行数值比较已经有4人回复
matlab解微分方程已经有10人回复
matlab数值求解边界条件微分方程组已经有7人回复
高金求助matlab解微分方程组已经有12人回复
matlab解微分方程组已经有15人回复
【求助】matlab 求解微分方程中的未知参数已经有20人回复
【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题已经有21人回复
【求助】MATLAB 有限差分法（FDM）求解偏微分方程已经有22人回复
【求助】积分微分方程matlab求解已经有6人回复
【求助】求助一个一阶常微分方程的初值问题已经有16人回复
【求助】matlab怎么求解偏微分方程组啊，先谢谢了已经有13人回复
【求助】用MATLAB求解偏微分方程的数值解悬赏20 已经有8人回复
【求助】大家帮忙推荐一本微分方程数值解的经典教材已经有7人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2012-05-24 01:01:13

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

lijie169

铜虫 (著名写手)

应助: 159 (高中生)
金币: 161.5
帖子: 1338
在线: 491.5小时
虫号: 1180408

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

http://www.nd.edu/~nancy/Math306 ... Mlb8a/DEwMlb8a.html
减小区间大小

赞一下