中易网

用mathematica解一个简单的二阶偏微分方程

答案:1  悬赏:60  
解决时间 2021-01-17 15:54
  • 提问者网友:孤山下
  • 2021-01-17 01:20
用mathematica解一个简单的二阶偏微分方程
最佳答案
  • 二级知识专家网友:渡鹤影
  • 2021-01-17 01:35
1 你的代码里混了中文标点。
2 你的方程是热传导方程,它的解析解一般是级数解。Mathematica截止目前,是不用级数来表示方程的解的。(软件的这种处理方法可能和级数的收敛判定困难有关——Mathematica是个非常严谨的数学软件。)所以DSolve无法求解你的方程。
3 退一步讲,即使你想补上a的具体数值,使用NDSolve来求解这个方程的数值解,在你所给的条件下,这也是做不到的。如果你学习过偏微分方程的相关知识,或者你手头有《数学物理方程》之类的课本,你就会知道,你所给的限制条件,不属于教科书里通常会给出的限制条件的任何一种。如果你具备更深入的有限差分方面的知识,你就会知道,仅仅给出三个孤立的点上的函数值,也是根本无法求得这个方程的定解的。你的限制条件是你随手给的?还是你只是单纯地写错了条件?总之你再检查检查吧。
追问:多谢你详细的解答。我问了其他人,确实是我的初始条件有问题。如果把初始条件近似成
p(x,0) = \delta(x) 方程就有了一个通解,是一个高斯函数,指数为something*x^2/t . 我之前参考国内数学物理方法分离变量法解的话,那个常数不好确定,有人通过对那个常数积分也能得到同样的高斯函数,可是意义不明确。对x进行傅立叶变换的话,求解方程变的相对容易很多。
我要举报
如以上回答内容为低俗、色情、不良、暴力、侵权、涉及违法等信息,可以点下面链接进行举报!
点此我要举报以上问答信息!
大家都在看
推荐信息