自由电荷体密度计算方法