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

登录注册

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

24小时热门版块排行榜

返回列表

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

回复此楼

» 猜你喜欢

今天系统多次维护，明天很可能放榜！已经有5人回复
欢迎发来filecode的Mz6后的代码验证其规律已经有74人回复
岩土工程学报什么时候才能终审完呐已经有4人回复
时间戳又变了8-15 已经有25人回复
2027广东省杰青已经有6人回复
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。已经有3人回复
哪位老哥知道今年的国自然具体哪一天放榜？已经有9人回复
咱们一起用铁证分析2026国家社科基金中标与否已经有26人回复
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人已经有16人回复
有时候，自然基金真的不能太认真（我的申报经验）已经有12人回复

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

高手进来用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
减小区间大小

赞一下