马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
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问题,同学介绍来这个论坛,还望各位高手能帮小弟解答一下。万分感谢! |