Hubbard模型铁磁态平均场解法