我想 用。在使用 Galerkin 近似时，我得到 其中是质量矩阵，是刚度矩阵。我注意到我可以写成

$$\alpha u_{tt} = -\nabla^2u$$ $\frac{\partial u}{\partial n} = 1$ $$M\ddot{c}=\frac{1}{\alpha}(Dc+b)$$ $M$$D$ $$\dot{U} = ZU+B$$

其中和在右上角有对角线块：和左下角和在给出的向量。$U = [c, \dot{c}]^T$$Z$$I$$M^{-1}D$$B$$M^{-1}b$

那么我们可以继续在这个问题上使用隐式有限差分方案。（？）

我知道我的解释很少，但我希望以前可能解决过这个问题的人可以确认这是否是一种可行的方法。我的解决方案看起来非常“不稳定”。

也许我使用这种方法真的很天真。有没有简单的方法来解决这个问题？