回 帖 发 新 帖 刷新版面

主题:[讨论]这样的二阶偏微分方程如何数值模拟?

我的数学和编程基础很差,勉强可以用matlab解微分方程,偏微分就找不到北了,好多概念性的东西搞不清。无奈目前做的一个课题碰到了这么一个方程,希望大家可以发表自己的看法、集思广益,真切的希望能得到各位高人的帮忙!

我遇到的问题是:对于下面的这个二阶偏微分方程,如何用Matlab(或者Fortran语言)去模拟A(x,y)随x、y的变化?

A''x(x,y)+A''y(x,y)+f[A'x(x,y), A'y(x,y)] A(x,y)=0

其中,
A''x(x,y)、A''y(x,y) 分别是A(x,y) 对x,y的二阶偏导数;
A'x(x,y)、A'y(x,y) 分别是A(x,y) 对x,y的一阶偏导数;
f[A'x(x,y), A'y(x,y)]表示f是与A'x(x,y)和A'y(x,y)都相关的函数。
只知道方程在x=0、y=0两个边界上的初始条件A(x=0,y)、A(x,y=0)、A’(x=0,y)、A’(x,y=0)。 

回复列表 (共1个回复)

沙发

看过Matlab的PDE工具箱的介绍,里面提供了一些最简单、经典的偏微分方程,如:椭圆型、双曲型、抛物型等少数的偏微分方程,并给出了求解方法,但仍然是一头雾水:

(1) 用Matlab的PDE工具箱解偏微分方程时都是先画出方程的求解区域再在此区域内求解,而我所要模拟的方程大致是一个衰减波的形式,其衰减深度(衰减的快慢)随着A’x(x,y)和A’y(x,y)前的系数的变化而变化,所以似乎不能先画出求解区域再求解;

(2) 用Matlab的PDE工具箱解偏微分方程的第二步是设置边界条件(Dirichlet和Neumann两种),而我所要模拟的方程其实只知道方程在x=0、y=0两个边界上的初始条件,而不知道其他的边界条件。

综上所述,我现在有以下疑问:
是否能用Matlab的PDE工具箱模拟上面这个只知道初始条件的二阶偏微分方程?如果能,边界条件和求解区域怎么处理?如果不能,用其他什么办法可以模拟?

盼望大家的指教! 

我来回复

您尚未登录,请登录后再回复。点此登录或注册