matlab解耦合方程,如何在Matlab中求解這些耦合微分方程？

我有三個偏微分方程(PDE)和一個如圖所示的變量的解析解。 (x，y，t)，p(x，y，t)，C_ {a}(x，y，t)和C_ {b}(x，y，t)即在空間和時間方面。如何在Matlab中求解這些耦合微分方程？

我知道在Matlab中有一個函數pdepe()來解決一維拋物型橢圓偏微分方程的初邊值問題。我想知道這個函數或Matlab中的其他函數如何用於解決下面描述的二維和耦合問題。

問題：

以下兩個等式分別表示用於兩個物種a和b，偏微分方程：

凡D_ {H}和q給出爲：

這裏，R_ {A} = R_ {B} = R，其中R由下式給出：

最後，最後一個方程給出如下：

初始和邊界條件：

域名總大小爲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} = 0.2302mol m^-3的常數C，而邊界(4)處的濃度被設定爲C_ {a} = 1mol m^-3和C_ {b} = 0.4603mol m^-3。邊界(1)處的濃度由邊界(3)和(4)之間的恆定梯度確定，而對流通量邊界條件$$(\ frac {\ partial C} {\ partial x} = 0)$$在(2)的出口處。

2012-10-24

Pupil

+0

留下評論我的回答本來是不錯的... –

+0

嗨Rody：我把你的答案投下來，因爲我認爲(我也希望你)，這是一個評論而不是答案。我將在下一次通過留言來通知您。對這次抱歉。 –

+0

我不知道你在說什麼。建議'pdepe'可以解決您的問題，因此它是一個答案。然後，你後來重新提出這個問題，以包括你知道它存在的事實，但只是不知道如何使用它，那麼......這是非常重要的信息。所以：你用這個函數試過了什麼？ –