mathematica绘制一维波动方程曲线,
对于d2u(x,t)/dt2-c^2*d2u(x,t)/dx2=0,0
人气:194 ℃ 时间:2020-10-02 07:46:52
解答
你给的初始条件这个 u(x,0)=d sin(2 pi x) 好像不对吧,d难道是系数吗?是的,d是个系数!如果你是要画Plot3D的话,那么应该是求微分方程的数值解,因此要指定c,d的值,并用NDSolve求解。 c = 2; d = 1; sol = NDSolve[{Derivative[0, 2][u][x, t] - 4*Derivative[2, 0][u][x, t] == 0, Derivative[0, 1][u][x, 0] == 0, u[x, 0] == Sin[2*Pi*x], u[0, t] == 0, u[1, t] == 0},u, {x, 0, 1}, {t, 0, Pi}, MaxStepSize -> 0.001]求解出来的结果是一个插值形式的函数,你要用替换的方式才能转化成Plot3D能画的函数。 注意解方程和画图,运行时间都是很慢的,因为求解PDE真的很复杂。 Plot3D[Evaluate[u[x, t] /. sol], {x, 0, 1}, {t, 0, Pi}, PlotPoints -> 45, Mesh -> {25, 30}, PlotStyle -> Directive[Yellow, Specularity[White, 20], Opacity[0.8]]]
推荐
猜你喜欢
- 帮我翻译成英语谢谢【 我爱我的外公和外婆】
- 1,1,3,18,18,0,) 这是个数字推理题怎么做?
- 张老师床上的灯光无论风霜雨雪,暑往寒来,每天都到深夜这句话是什么意思
- 松柏冢累累:累累读音是什么
- 黑色混合物A高温下变成红色固体B和无色气体C
- The policeman was_____the matter,______to find out the truth as soon as possible.
- 死尸回光中,PH3和氧气反应生成什么
- nice to meet u,求翻译