Excel批量计算月相求助：现有公式结果存在错误

问题根源 #

你之前使用固定平均朔望月长度（27.3或29.53）的公式会产生误差，因为实际朔望月长度并非固定值，它在29.27到29.83天之间波动，长期累积的偏差会导致月相判断错位。

解决方案一：改进的Excel公式（无需VBA） #

使用更精确的朔望月平均长度（29.530588853天，国际天文联合会认可的长期平均值），结合相位区间判断来匹配月相：

步骤1：计算相位比例 #

在空白单元格输入公式（假设目标日期在B2）：

=MOD((B2-DATE(2000,1,6))/29.530588853,1)

这个公式返回0到1之间的数值，代表目标日期相对于最近一次新月的相位比例：

• 0 = 新月
• 0.25 = 上弦月
• 0.5 = 满月
• 0.75 = 下弦月

步骤2：匹配月相编号 #

基于相位比例，用LOOKUP函数映射到你定义的月相编号：

=LOOKUP(MOD((B2-DATE(2000,1,6))/29.530588853,1),
{0,0.125,0.375,0.625,0.875},
{1,2,3,4,1})

区间设置说明：

• 0~0.125：新月（1）
• 0.125~0.375：上弦月（2）
• 0.375~0.625：满月（3）
• 0.625~0.875：下弦月（4）
• 0.875~1：回到新月（1）

这种区间缓冲可以抵消平均长度带来的小幅误差，适合大多数日常场景。

解决方案二：VBA自定义函数（高精度） #

如果需要更高精度（比如专业天文级匹配），可以用VBA实现Meeus天文算法，直接计算日期的精确月相：

步骤1：添加VBA代码 #

1. 打开Excel，按Alt+F11打开VBA编辑器
2. 右键点击工作簿名称 → 插入 → 模块
3. 粘贴以下代码：

Function GetMoonPhase(DateVal As Date) As Integer
    Dim JD As Double
    Dim T As Double
    Dim M As Double, Mprime As Double, F As Double
    Dim Phase As Double
    
    ' 转换日期为儒略日
    JD = DateVal + 2415018.5 + (TimeValue(DateVal) / 24) - 0.5
    
    ' 计算J2000.0以来的儒略世纪数
    T = (JD - 2451545#) / 36525#
    
    ' 计算太阳平近点角
    M = WorksheetFunction.Radians((357.52911 + 35999.05029 * T - 0.0001537 * T ^ 2))
    
    ' 计算月球平近点角
    Mprime = WorksheetFunction.Radians((134.96340251 + 477198.8675055 * T + 0.0087414 * T ^ 2 + 0.00000186 * T ^ 3))
    
    ' 计算日月平距角
    F = WorksheetFunction.Radians((93.2720950 + 483202.0175233 * T - 0.0034029 * T ^ 2 - 0.00000317 * T ^ 3))
    
    ' 计算月相比例（0=新月，0.25=上弦，0.5=满月，0.75=下弦）
    Phase = (F + 2 * WorksheetFunction.Sin(M) * (0.1644 + 0.000304 * T) - 2 * WorksheetFunction.Sin(Mprime) * (0.0251 + 0.000014 * T) + 2 * WorksheetFunction.Sin(2 * F - Mprime) * (0.0118 + 0.000004 * T) - WorksheetFunction.Sin(2 * F) * 0.0041 - WorksheetFunction.Sin(M - Mprime) * 0.0032) / (2 * 3.1415926535)
    Phase = Phase - Int(Phase)
    If Phase < 0 Then Phase = Phase + 1
    
    ' 映射到月相编号
    Select Case Phase
        Case Is < 0.125: GetMoonPhase = 1 ' 新月
        Case Is < 0.375: GetMoonPhase = 2 ' 上弦月
        Case Is < 0.625: GetMoonPhase = 3 ' 满月
        Case Is < 0.875: GetMoonPhase = 4 ' 下弦月
        Case Else: GetMoonPhase = 1 ' 新月（循环）
    End Select
End Function

步骤2：使用自定义函数 #

在单元格中输入公式即可获取月相编号：

=GetMoonPhase(B2)

这个方法计算精度极高，适合专业项目的大量数据处理。

内容的提问来源于stack exchange，提问作者Muller