我正在尝试对 Bogoliubov-de Gennes Hamiltonian 进行对角化。问题是哈密顿量包含拉普拉斯算子。这可以通过使用离散的拉普拉斯算子来解决
如果我使用一个简单的哈密尔顿算子（只是拉普拉斯算子）并使用eig(l1dd(...,...))离散l1dd的拉普拉斯算子（2 和 -1's），如https://people.sc.fsu.edu/~jburkardt/m_src/laplacian/laplacian.html，那么我得到：
我期望线性或二次色散（$\omega \propto k 或 k^2$），但我没有得到。可能是什么问题呢？看起来eig做得很好，所以可能l1dd是问题所在？