我有三个偏微分方程（PDE）和一个洳图所示的变量的解析解 （x，yt），p（xy，t）C_ {a}（x，yt）和C_ {b}（x，yt）即在空间和时间方面。
我知道在Matlab中有一个函数pdepe()来解决一维抛物型椭圓偏微分方程的初边值问题我想知道这个函数或Matlab中的其他函数如何用于解决下面描述的二维和耦合问题。
以下两个等式分别表示用于两個物种a和b偏微分方程：
最后，最后一个方程给出如下：
域名总大小为10厘米x 5厘米y形子域的宽度为0.5厘米。这个子域的初始\ phi为0.50而在周围的矩阵\ phi = 0.26。 1Pa和0Pa的常数p分别保持在边界（1）和（2）对应于大约10^-3m m^-1的梯度。边界（3）和（4）上的p由边界（1）和（2）之间的线性梯度确定在边界（3）保持C_ {a} = 2mol m^-3和C_ {b} =
留下评论我的回答本来是不错的... –
嗨Rody：我把你的答案投下来，因为我认为（我也希望你）这是一个评论而不是答案。我将在下┅次通过留言来通知您对这次抱歉。 –
我不知道你在说什么建议'pdepe'可以解决您的问题，因此它是一个答案然后，你后来重新提出这个問题以包括你知道它存在的事实，但只是不知道如何使用它那么......这是非常重要的信息。所以：你用这个函数试过了什么 –
