| 查看: 962 | 回复: 4 | |||
[交流]
【求助】求解常微分方程出现警告，求助
|
|
求解如下微分方程时，出现警告 Warning: Failure at t=4.893789e-005. Unable to meet integration tolerances without reducing the step size below the smallest value allowed (1.738623e-019) at time t. 无法得到要求的解区间，请帮忙分析一下。 Y0=[0.0001; 15.2571]; tspan=[0.4894e-4 0]; [t,YY]=ode23s('DYdt',tspan,Y0); plot(YY(:,1),t) function Yd=DYdt(t,Y) Yd=[Y(2); (1+Y(2)^2)/Y(1)+ 9.4597e+008*(1+Y(2)^2)^1.5]; |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有16人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有14人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有15人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有17人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有21人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有21人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求解二阶非线性常微分方程！！！解决之后另有重谢！！！
已经有5人回复
- 哪位高手可以帮我解答，常微分方程，解的延拓课后习题求解
已经有11人回复
- 求解二阶变系数常微分方程
已经有7人回复
- 求助非线性常微分方程和常系数双曲偏微分方程求解
已经有7人回复
- 常微分方程求解公式
已经有7人回复
- 求解方程，在线等……求大神指点！
已经有7人回复
- Matlab数值求解二阶常微分方程
已经有9人回复
- Matlab数值求解非线性常微分方程
已经有5人回复
- 求助，matlab求解二元二阶的常微分方程组
已经有9人回复
- 求助ANSYS求解总是警告
已经有6人回复
- 一个二阶常微分方程的求解方法？
已经有5人回复
- matlab求解常微分方程
已经有6人回复
- 二阶常微分方程求解
已经有5人回复
- 1stopt4阶非线性常微分方程，帮小弟看一下，跪求。
已经有7人回复
- 高手进来用matlab求解常微分方程，十分感谢！
已经有17人回复
- 隐式欧拉法求解一阶常微分方程
已经有7人回复
- 1stOpt求解常微分方程边值问题
已经有5人回复
- matlab 求解隐式常微分方程时用ode15i 函数出现的问题求助
已经有3人回复
- 【求助】ansys求解时出现警告：有4个小支点调解方程求解，原因是？
已经有3人回复
- 【求助】出现以下警告一般是哪些设置出现了问题？
已经有4人回复
- 【求助】常微分方程组求解中系数与某变量值关联的问题
已经有12人回复
- 【求助】如何用Runge-Kutta迭代求解二阶常微分方程组【已解决】
已经有9人回复
2楼2011-03-11 15:55:53
3楼2011-03-14 09:07:47
★ ★
xiegangmai(金币+2): 谢谢应助！ 2011-03-14 12:37:04
micronano(金币+5): 2011-03-14 16:35:29
xiegangmai(金币+2): 谢谢应助！ 2011-03-14 12:37:04
micronano(金币+5): 2011-03-14 16:35:29
|
我也不知道你那常微分方程组写得对不对？再核实下。 按照你的方程组形式和初始条件帮你简化了下，两三行就可以解决问题了。 --------------------------------------------------------------------------------------------------- function main fun = @(t,Y)[Y(2);(1+Y(2)^2)/Y(1)+9.4597e+08*(1+Y(2)^2)^1.5]; [t,Y] = ode45(fun,[0.4894e-4 0],[0.0001 15.2571]); plot(Y(:,1),Y(:,2)) --------------------------------------------------------------------------------------------------- 但结果同样遇到问题Warning: Failure at t=4.893789e-005. Unable to meet integration tolerances without reducing the step size below the smallest value allowed (1.738623e-019) at time t. ------------------------------------------- 分析： 1. 按提示需要减小step size 2. 检查你的方程和初值条件，如果都没问题，那么即使warning提示，也不表示一定是错误 |
4楼2011-03-14 10:01:32
5楼2011-03-14 10:13:04
50