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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 计算数学 » 请大家帮忙分析一下我这个matlab解微分方程初值的例子,到底问题出在哪里
21/1返回列表
查看: 494  |  回复: 1
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

iphysunba

新虫 (初入文坛)

[交流] 请大家帮忙分析一下我这个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
回复此楼

» 猜你喜欢

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

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

查看全部散金贴
1楼 2012-05-24 01:01:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijie169

铜虫 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
http://www.nd.edu/~nancy/Math306 ... Mlb8a/DEwMlb8a.html
减小区间大小
赞一下 回复此楼
2楼2012-05-26 15:43:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 iphysunba 的主题更新
21/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭