考虑两个情形: 其中, 保持粒子数守恒,模拟时请使用“好量子数”。 不保持粒子数守恒。 现在模拟方式确定,并通过验算了建模的正确性。 matlab建模完成 程序目录位于 e2019 2021thermalnonlocality 1Dchainsuperconductor