如何在'pdepe'中使用NxM矩阵作为初始条件

Question

我正在使用pdepe求解器同时解决4个PDE。 现在我想使用该解决方案作为下一步的初始边界条件，例如u0（1，：）= sol（end，：，1）这很有效，只要我只使用单个值作为初始值条件，例如u0（1）= sol（end，end，1），但是一旦我尝试整个向量就失败了。

这似乎实际上非常相似，但在我的情况下不起作用： 如何在`pdepe中使用NxM矩阵作为我的初始条件？

sol = pdepe(m,...,@(x)icfun(x,sol)...
...

%icfun
u0=icfun(x,sol)

u0(1,:)=sol(end,:,1)
u0(2,:)=sol(end,:,2)
u0(3,:)=sol(end,:,3)
u0(4,:)=sol(end,:,4)

end

什么不起作用，是我解释它的方式从上面的链接提示。

A=sol;
a=sol(end,:,1);

sol = pdepe(m,...,@(x)interp1(a,A',x,'pchip')'....

我希望能够将矢量用作仅一个pde的初始边界条件，这为4个PDE提供了NxM矩阵。 相反，会发生错误，说明初始边界条件需要是列向量。

'使用pdepe时出错（第231行）ICFUN的输出无效。 ICFUN必须返回一个列向量。

如何交出整个向量而不是单个值？