PDE system - pdepe system example - error or incomplete description on Matlab website?

3 次查看（过去 30 天）
Leonid
Leonid 2019-12-19
评论： Bill Greene 2019-12-24
I am no longer sure how to use Matlab's pdepe properly for solving a system of PDEs after checking the example given on Matlab's help website.
If you look to an example here https://www.mathworks.com/help/matlab/math/solve-angio-pde-system.html under the 'Code Equation' you'll see a rewritten pde system in a matrix form. Comparing it to the original equation it is obvious that either the term " a*dndx*dcdx " is missing from the first (rewritten) equation, or Matlab implicity assumes that the derivative is acting only on the dudx terms when it deals with pde systems (this will be very weird). Could anyone with an extensive experiece in solving pde systems in Matlab clarify this?
Thanks

请先登录，再回答此问题。

采纳的回答

Bill Greene
Bill Greene 2019-12-20
Looks like you have found an error in this documentation page. If you look at equation 3 in the referenced paper you will see that this equation corresponds to what is actually coded in the pdepe example. The results from the pdepe analysis agree with figures 3 and 4 in the referenced paper. It is the first PDE in the documentation page that is incorrect.
  2 个评论
Leonid
Leonid 2019-12-23
Thanks for your answer. Based on the referenced paper it is the documentation page that contains error.
One additional thing that will be helpful to clarify:
If the coefficients d and a were functins of x [for example defined by some discrete data (and then interpolated within the function)], whould their derivatives through the flux term be taken into account by Matlab pdepe?

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Boundary Conditions 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by