问题是odeint解决了形式的一阶ODE

$$\dot u(t) = Au(t),$$ 但波动方程是（在空间离散化后）形式为的二阶ODE $$\ddot{u}(t) = Au(t).$$ （事实上​​，你的图对于抛物线方程的解看起来很合理。）要应用黑盒 ODE 求解器（顺便说一句，这真的不是一个好主意），你首先需要重写波动方程作为一阶方程组 $$\left\{\begin{aligned} \dot{u}(t) &= v(t),\\ \dot{v}(t) &= Au(t).\end{aligned}\right.$$

冒着听起来粗鲁的风险：在对此类代码进行细微修改之前，您应该首先阅读一些 PDE 以了解您要解决的问题。一个好的起点可以是Randy LeVeque的附录 E，普通和偏微分方程的有限差分方法，SIAM 2007。

由于波动方程涉及时间上的二阶导数，因此需要两个初始条件：初始位移$U(x,0)$和初速度$U_t(x,0)$. 我没有看到您的问题中指定的初始速度条件。这使我怀疑您链接到的教程有问题...

正如我所怀疑的，本教程中的代码是用于热方程，而不是波动方程。为了求解波动方程，您还需要完全使用不同的时间步进方案。教程中实现的那个不适用于波动方程。