我正在尝试使用 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 进行编码。谢谢
我正在尝试使用 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 进行编码。谢谢
通常,您会将狄利克雷条件作为初始条件。如果对应于每个索引的初始条件在某些集合中,您将设置其中的每一个s 为一个固定值(狄利克雷初始条件),然后确保设置为零,并且第 行是个单位矩阵的第 th 行。