规定奇怪的边界条件

时间:2010-10-27 07:06:17

标签: wolfram-mathematica

有没有人知道如何使用NDSolve在Mathematica中规定类似u [t,0,y] == u [t,1,1-y]的边界条件......它总是抱怨因变量的参数应该字面上匹配自变量。

提前致谢。

1 个答案:

答案 0 :(得分:0)

这种对称条件可能以Derivative [0,1] [u] [x,1/2] == 0的形式重铸。当然,关于这个问题的更多信息会有所帮助。

编辑以响应rcollyer: (0,1)中所有x的代数同一性f(x)= f(1-x)意味着几何对称性:f的图将关于线x = 1/2对称。现在绘制这样一个函数的图形;如果它是可微分的,你会发现f'(1/2)= 0。

现在,我不确定OP的问题可以通过这种方式重铸;它取决于问题的具体细节。当处理磁盘上的PDE时,经常出现这种情况,其中函数u是极坐标r和θ的函数。如果磁盘代表一个夹紧的鼓,那么也许你的u(1,t)= 0。但是,你(0,t)呢?如果函数是对称且平滑的,则u_x(0,t)= 0是合理的条件。