|
|
[交流]
请大家帮忙分析一下我这个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一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有7人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有7人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有19人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有8人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
» 抢金币啦！回帖就可以得到:
查看全部散金贴
|
回复此楼
+1/336
