【讨论】这个偏微分方程用matlab怎么解？

exploiter

这个偏微分方程用matlab怎么解？

taohe

安装pde toolbox了吗？如果有的话，先在matlab里demo一下以获得大概了解，然后再参考这个工具库的文档。不是太难。
  
:I

exploiter

看了
还没看明白
  
有没有人见过关于传热、传质、吸附、解吸、流体运动等方面matlab参考书？
急需，请推荐。

exploiter

没人回答
  
先不说用matlab求解，pde tools 估计能解。
关键是要把这个微分方程用simulink仿真，simulink里没有偏微分模块
拉普拉斯变换把偏微分变成全微分？变完后是对位移x的全微分，不是对时间t,
还是没法仿真。
  
问了matlab大观园薛老师，也没有好办法。
谁做过类似的？郁闷啊

承辕

Matlab偏微分方程的解法参考书，我到是见过一本，不知道能不能给你一帮助！
《偏微分方程的matlab解法》武汉大学出版社 陆群安等编

taohe

有没有看到pdedemo5?那是个关于传热的例子。
  
一般来说在matlab里doc pde，然后看看get started，在用我上面说的demo运行一下关于传热的例子，不理解就把pdedemo5.m打印出来研究研究。
  
再matlab里用pde toolbox价偏微分方程一般来说有两种方法：GUI方式和命令行方式。GUI方法是使用pdetool，其用法参考"help pdetool" 或者"doc pdetool"
  
你的问题属一维传热问题，在pde toolbox里解应该是小菜一蝶。
看看：
help pdetool，看看怎么使用它来定义你要解的问题（偏微分方程和边界条件），定义完问题后就可以解了。
  
:I

exploiter

跟传热那个不一样
这个是类似对流问题
  
看看matlab网站上这么说的：
  
Solution Number: 4382  
Date Last Modified: 2002-02-15     
Product:    Partial Differential Equation 1.0 ==> Current Version  
Platform:    All Builds  
  
Problem Description
  
Can the Partial Differential Equation (PDE) Toolbox handle convection coefficients (terms with first order space derivatives) in heat transfer problems?
  
Solution:
  
There is no way of entering convection terms in the PDE Toolbox. Convection terms are first order space derivatives (such as du/dx); these are different than first order time derivatives, which can be handled in the parabolic equation.  
  
For static problems, it is often possible to formulate a nonlinear problem that contains the convection terms by putting the first order space derivative in the 'f' coefficient of the standard equation. The first order space derivatives are available as ux and uy in the expression.  
  
If the problem is not too strongly hyperbolic, the nonlinear solver will probably converge. If you turn on full Jacobian, convergence is even more likely.
  
看来对空间变量一阶比二阶的还不好办啊

taohe

原来这样，是我搞错啦。:I
  
现在看来好像只是没法告诉pde toolbox我们要解的pde。如果要解的话可以这样：
  
1。使用相关的方法自己解，比如用差分法。
  
2。试试femlab，是第三方公司出的有限元工具库，用在matlab。看过演示，好像比matlab的pde toolbox强。
  
:I

exploiter

是的femlab应该能解
但是我不想再利用别的工具了，就利用matlab和simulink，最终的目标是要把一个模型做成类似matlab的工具箱

tylangzi

我试了，可以解，1D和2D的
到Femlab里面选择coefficient形式，设置方程前面的系数，就可以表达出这个方程，然后画出模型，就可以求解了，我做了一个最简单的，u/r=10的，可以解