问题描述
我正在尝试求解二维热方程,我正在使用 ode15 求解,有人指示我必须调整 dT/dt 方程。是否可以获得一些有关我应该如何处理的指导?
我遵循了 Matlab 中的这个示例:https://uk.mathworks.com/matlabcentral/fileexchange/59916-simple-heat-equation-solver?s_tid=srchtitle
怎么可能
2
示例中的 getRHS 函数的一部分要更改以限制其升高的温度?
解决方法
只需将其添加在 ode 之前作为 thr
的阈值 u0
,例如:
thr=0.5; % just some # I picked,u0(u0>thr)=thr;