狄利克雷边界条件

计算科学 matlab
2021-12-11 05:53:18

我正在尝试使用 ode15s 在 matlab 中解决 ODE。而不是以格式指定 ODE

  M * dC/dt = f(C,t)  where C is a function of x and t.

我想用

  M * dC/dt = J*C

由于向量 C 非常大,我预计计算时间会大大减少。我试图弄清楚如何对 Dirichlet 边界条件 C(x=0,t) = a 进行编码。谢谢

1个回答

通常,您会将狄利克雷条件作为初始条件。如果Ci对应于每个索引的初始条件i在某些集合中,您将设置其中的每一个Cis 为一个固定值(狄利克雷初始条件),然后确保fi设置为零,并且i第 行M是个i单位矩阵的第 th 行。