24小时热门版块排行榜    

Znn3bq.jpeg
查看: 365  |  回复: 4
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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)。

[ Last edited by cuplgz on 2008-2-25 at 19:22 ]
回复此楼

» 猜你喜欢

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

redcrag

木虫 (著名写手)

找找时域有限差分法(FDTD)的资料看看其中的吸收边界条件(absorbing boundary condition )(现在流行perfectly matched layer边界条件),简单说就是用有限区域模拟波的传播,因此边界上要保证没有或仅有很小的反射。
引用回帖:
Originally posted by xinkaifish at 2007-8-23 10:40 AM:


确实是光波传播的过程。
“或者是采用吸收边界”,能解释一下什么是“吸收边界”吗?

http://redcrag.spaces.live.com/
5楼2007-08-24 08:07:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

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 15:58:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

redcrag

木虫 (著名写手)

Matlab的PDE工具箱能求解若干种典型的方程,如果你的不在其内,估计是不能用的。另外,你的边界条件如果不全的话,理论上是不是能保证解的唯一性?不过看你的问题似乎是波传播的过程,边界一般是零值或者是采用吸收边界。
http://redcrag.spaces.live.com/
3楼2007-08-23 08:08:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xinkaifish

木虫 (职业作家)

引用回帖:
Originally posted by redcrag at 2007-8-23 08:08:
Matlab的PDE工具箱能求解若干种典型的方程,如果你的不在其内,估计是不能用的。另外,你的边界条件如果不全的话,理论上是不是能保证解的唯一性?不过看你的问题似乎是波传播的过程,边界一般是零值或者是采用吸 ...

确实是光波传播的过程。
“或者是采用吸收边界”,能解释一下什么是“吸收边界”吗?
4楼2007-08-23 10:40:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料与化工371求调剂 +17 陪琳看海 2026-04-04 18/900 2026-04-10 07:51 by 314126402
[考研] 080100力学316求调剂 +5 L_Hairui 2026-04-07 5/250 2026-04-10 05:48 by 小小树2024
[考研] 一志愿2110,化学学硕310分,本科重点双非求调剂 +18 努力奋斗112 2026-04-08 18/900 2026-04-09 23:28 by wolf97
[考研] 297求调剂 +22 ljy20040718! 2026-04-03 24/1200 2026-04-09 20:48 by yanenwang
[考研] 求调剂希望还是希望在山河四省附近 +3 快乐的小白鸽 2026-04-05 3/150 2026-04-09 17:36 by wp06
[考研] 337求调剂 +4 Gky09300550, 2026-04-09 4/200 2026-04-09 17:18 by 帕尔马拉特
[考研] 083200 初试305分 求调剂 暂不考虑跨专业 +15 Claireyyyy 2026-04-09 15/750 2026-04-09 16:11 by zhuimr
[考研] 求调剂 +8 吃口冰激凌 2026-04-07 8/400 2026-04-09 08:03 by 5268321
[考研] 353求调剂 +8 晴空万里air 2026-04-07 8/400 2026-04-09 00:18 by GouQ
[考研] 软件工程求调剂22软工296分求调剂,接受跨调 +4 yangchen2017 2026-04-08 5/250 2026-04-08 21:56 by 土木硕士招生
[考研] 277、学硕,求调剂 数一104, +11 瓶子PZ 2026-04-07 12/600 2026-04-07 23:30 by 一只好果子?
[考研] 288环境专硕,求调材料方向 +35 lllllos 2026-04-04 39/1950 2026-04-07 23:24 by 一只好果子?
[考研] 307求调剂 +3 Youth@@ 2026-04-07 3/150 2026-04-07 22:00 by hemengdong
[考研] 085602调剂 初试总分335 +10 19123253302 2026-04-05 10/500 2026-04-07 15:23 by 小乔同学ya
[硕博家园] 0856材料化工求调剂,一志愿211,初试成绩349 +4 江淮北月 2026-04-05 4/200 2026-04-06 22:44 by chenzhimin
[考研] 071000生物学调剂 +7 拉提桃 2026-04-06 7/350 2026-04-06 18:55 by 52305043001
[考研] 求调剂到材料 +5 程9915 2026-04-06 5/250 2026-04-06 15:21 by yulian1987
[考研] 材料工程310专硕调剂 +14 捞捞我…. 2026-04-04 15/750 2026-04-06 14:18 by lqwchd
[考研] 316求调剂 +9 墨辰_Orion926 2026-04-04 9/450 2026-04-04 21:35 by lbsjt
[考研] 复试调剂 +6 范根培 2026-04-04 6/300 2026-04-04 14:27 by 土木硕士招生
信息提示
请填处理意见