加盖温度-带ode15的二维热方程

Question

我正在尝试求解二维热方程，我正在使用ode15进行求解，我被指示必须调整dT/dT方程。有可能得到一些关于我应该如何做的指导吗？

我从Matlab中学习了这个例子：https://uk.mathworks.com/matlabcentral/fileexchange/59916-simple-heat-equation-solver?s_tid=srchtitle

怎么可能

du = alpha/dx^2*(u(1:end-2,2:end-1)-2*u(2:end-1,2:end-1)+u(2:end-1,3:end)...
+ u(2:end-1,1:end-2)-2*u(2:end-1,2:end-1)+u(3:end,2:end-1));
du = du(:);

更改示例中getRHS函数的一部分，以限制其上升到的温度？

Answer 1

只需将其添加到ode之前，作为u0的阈值thr，例如：

thr=0.5; % just some # I picked, 
u0(u0>thr)=thr;