声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1654|回复: 3

[编程技巧] [求助]二维椭圆型偏微分方程求解的问题

[复制链接]
发表于 2009-8-11 21:23 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
-▽·(c▽u)+au=f 形式的方程,用函数assempde(b,p,e,t,c,a,f)求解
其中c a f是偏微分方程的参数
b是边界条件
p e t是网格划分
具体的不是很清楚,看了一些实例,还是不太明白

我想问的是
1.看到的几个实例中c a f 都是常数,
不过我需要求解的方程是-▽·(c▽u)+au=f(x,y)的形式,并且f是一个离散函数,这个要怎么办?
(假设f(x,y)和已知矩阵M的元素一一对应)

2.边界条件的设置,实例中看到有这些
g='circleg';     %  定义单位圆区域
b=’circleb2’;  % x^2 on the boundary
另一种
g='squareg';     %  定义单位方形区域
b='squareb3';    % 左右零边界条件,顶底零导数边界条件
我希望是矩形区域,边界四周函数值都是0,应该怎么定义b?
还有g定义区域,单位方形的时候是不是x,y都是-1~1的范围?

3.网格划分
[p,e,t] = initmesh(g);
这个完全不明白,p,e,t分别代表什么参数?
如果我要按照上面的M矩阵的规模(比如M是256*256的矩阵)来划分网格,该怎么定义?

4.最后我怎么调出u(x,y)的值啊?
貌似实例中都将pedsurf(p,t,u)画个图,不过我要具体数值

初次接触PDE问题,同学介绍来这个论坛,还望各位高手能帮小弟解答一下。万分感谢!
回复
分享到:

使用道具 举报

 楼主| 发表于 2009-8-12 16:58 | 显示全部楼层
额。。。不知道我有没有把问题讲清楚。。。
发表于 2009-8-12 21:43 | 显示全部楼层
个人水平专业有限, PDE Toolbox没使用过!
最近工作很忙, 实在没空看help文件学习! sorry
耐心待高人路过
发表于 2009-8-13 22:45 | 显示全部楼层
可以查一下武汉大学陆君安老师的书。
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-29 22:00 , Processed in 0.061318 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表