我如何使用NxM矩阵作为我在`pdepe`中的初始条件
How can I use NxM matrix to be my initial condition in `pdepe`

根据pdepe文档 ，求解器的初始条件函数具有以下语法：

u = icFun(x);

在列向量u返回指定x的PDE初始值。 因此，只有当PDE是具有M空间网格点的N未知数的系统时，初始条件才会是N x M矩阵。

因此，可以使用N x M矩阵来填充初始条件，但是将需要一些映射将给定列与x的特定值相关联。 例如，在调用pdepe的main函数中，可能存在

% icData is the NxM matrix of data
% xMesh is an 1xM row vector that has the spatial value for each column of icData
icFun = @(x) icData(:,x==xMesh);

这种方法的唯一缺点是初始条件的网格以及因此的pdepe解都受到初始数据的约束。 这可以通过使用如下插值方案来克服：

% icData is the NxM matrix of data
% xMesh is an 1xM row vector that has the spatial value for each column of icData
icFun = @(x) interp1(xMesh,icData',x,'pchip')';

其中存在转置以符合interp1对数据的interp1 。

u使用“线法”样式在每个网格上定义不同条件比使用pdepe更容易

MOL还可以更灵活地用于诸如3D问题之类的不同情况下：)）