| 查看: 561 | 回复: 1 | |||
[交流]
请大家帮忙分析一下我这个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 |
回复此楼» 猜你喜欢
40+中年人失业后，还能找到好工作吗？
已经有16人回复
-
人工智能虚构史料、伪造文献频现，多家期刊声明抵制
已经有4人回复
-
CsPbI2Br全无机钙钛矿SEM图像异常分析
已经有8人回复
-
请问社科规划办官网上的“项目查询模块”与“项目申报模块”的账号密码一样吗？
已经有3人回复
-
用NBS溴化时溶剂除了CCl4外还有什么
已经有9人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有5人回复
-
申请科研项目纯靠实力能成功吗？
已经有12人回复
-
极限拉扯，如何抉择
已经有15人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 高手进来用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人回复
2楼2012-05-26 15:43:57