24小时热门版块排行榜    

查看: 88  |  回复: 1
当前主题已经存档。
【悬赏金币】回答本帖问题,作者xinkaifish将赠送您 10 个金币

xinkaifish

木虫 (职业作家)

[求助] 这样的二阶偏微分方程如何数值模拟?

我的数学和编程基础很差,勉强可以用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)。

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xinkaifish

木虫 (职业作家)

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

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

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

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

盼望大家的指教!
2楼2007-08-22 10:54:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xinkaifish 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
信息提示
请填处理意见