我想在数值格式的每一步提取状态变量的值。
deSolve R包中的ode()函数,使用实现的ODE求解器之一对一个常微分方程组进行数值求解。为了做到这一点,它使用了一个动态调整的集成步骤,该步骤基于每个集成末尾的本地截断错误的值。用户基本上指定了所需的输出时间,它定义了时间步长网格,它可以有一个等于或小于输出时间的步骤。如果我们以Lotka Vo
我试图用隐式有限差分法来模拟木材圆柱体内的热传导。对于圆柱形和球形形状,我使用的一般热方程是:其中p为形状因子,p=1为圆柱,p=2为球面。边界条件包括表面对流。有关该模型的更多细节,请参见下面Matlab代码中的评论。of char production from wood (rho/s)
运行上述代码可以得到木材圆柱体中心和表面的温度分布:从这幅图中可以看出,由于某种原因,中心温度和表面温度在2分钟内迅速会聚,这是不正确的