24小时热门版块排行榜    

Znn3bq.jpeg
查看: 362  |  回复: 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

木虫 (著名写手)

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

xinkaifish

木虫 (职业作家)

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

确实是光波传播的过程。
“或者是采用吸收边界”,能解释一下什么是“吸收边界”吗?
4楼2007-08-23 10:40:56
已阅   回复此楼   关注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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 284求调剂 +4 让我上岸吧阿西 2026-04-09 4/200 2026-04-09 11:52 by 5268321
[考博] 博士自荐 +3 可可小胖 2026-04-08 3/150 2026-04-09 11:25 by juanzi_88
[考研] 085801 总分275 本科新能源 求调剂 +7 bradoner 2026-04-08 8/400 2026-04-09 11:17 by upczlm1989
[考研] 349学科化学045106求调剂,化学类都可以 +7 保好懂懂 2026-04-08 7/350 2026-04-09 07:37 by jp9609
[考研] 一志愿211,化学学硕,310分,本科重点双非,求调剂 +15 努力奋斗112 2026-04-06 17/850 2026-04-08 22:27 by 凯凯要变帅
[考研] 求机械专硕297第二批调剂 +3 拾柒12。 2026-04-08 3/150 2026-04-08 22:25 by laoshidan
[考研] 本科211,293分请求调剂 +12 莲菜就是藕吧 2026-04-03 13/650 2026-04-08 20:30 by 背对大海出发
[考研] 266求调剂 +5 08电气工程 2026-04-03 5/250 2026-04-08 20:22 by 逆水乘风
[考研] 320分人工智能调剂 +9 振—TZ 2026-04-03 10/500 2026-04-08 19:56 by 振—TZ
[考研] 求调剂 +14 张zic 2026-04-05 15/750 2026-04-08 16:14 by 一只好果子?
[考研] 307求调剂 +14 超级伊昂大王 2026-04-06 14/700 2026-04-08 07:03 by 无际的草原
[考研] 080500求调剂 +12 黄宇博 2026-04-06 12/600 2026-04-07 12:41 by upczlm1989
[考研] 一志愿太原理工大学计算机技术专硕348,求调剂指导 +3 nexious 2026-04-05 3/150 2026-04-07 08:19 by jp9609
[考研] 一志愿安徽某211 0703化学总分339求调剂 +7 晚风不晚 2026-04-04 7/350 2026-04-06 14:06 by houyaoxu
[考研] 313求调剂 +5 海日海日 2026-04-04 7/350 2026-04-05 13:58 by imissbao
[考研] 考研调剂 +5 四川王涛 2026-04-04 5/250 2026-04-04 22:18 by 啵啵啵0119
[考研] 341求调剂 +3 洛多罗 2026-04-02 4/200 2026-04-04 21:36 by 智能智慧
[考研] 085701求调剂 +7 龚禹铭 2026-04-04 8/400 2026-04-04 13:49 by 小小树2024
[考研] 266求调剂 +8 学员97LZgn 2026-04-03 8/400 2026-04-04 09:02 by 20021109
[考研] 317分 一志愿江南大学 化学工程学硕 求调剂 +6 YinTai 2026-04-03 6/300 2026-04-03 22:30 by 无际的草原
信息提示
请填处理意见