| 查看: 546 | 回复: 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 |
回复此楼» 猜你喜欢
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
面上没中，邀请各位路过的虫友分析一下分数
已经有13人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有8人回复
-
两块石头
已经有8人回复
-
初秋的晨风
已经有6人回复
-
学科评审组评审是指会评吗？
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有14人回复
-
广西大学-广州大学招聘博士后 欢迎广大优秀人才！！！
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 高手进来用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