python模拟月相亮面变化

要模拟月相亮面变化，需要考虑到月亮的公转和自转。我们可以根据当前日期计算出月亮的相位（即亮面的比例），然后根据相位绘制出亮面和阴面。

以下是一个简单的 Python 示例程序，可以根据输入的日期计算月相亮面变化：

import math

# 计算月相亮面比例
def moon_phase(year, month, day):
    # 基准日期：2000年1月6日12时
    base_year = 2000
    base_month = 1
    base_day = 6
    base_time = 12
    # 计算天数差
    days = (365.242 * (year - base_year) +
            30.6 * (month - base_month) +
            day - base_day)
    # 计算相位
    phase = days / 29.53
    return phase - math.floor(phase)

# 绘制月亮
def draw_moon(phase):
    # 设置画布大小
    canvas_width = 40
    canvas_height = 20
    # 计算亮面和阴面大小
    light_size = int(canvas_width * (0.5 + 0.5 * math.sin(phase * 2 * math.pi)))
    dark_size = canvas_width - light_size
    # 绘制亮面
    print('*' * light_size + ' ' * dark_size)
    # 绘制阴面
    print(' ' * light_size + '*' * dark_size)

# 输入日期
year = int(input('请输入年份：'))
month = int(input('请输入月份：'))
day = int(input('请输入日期：'))

# 计算相位并绘制月亮
phase = moon_phase(year, month, day)
draw_moon(phase)

运行程序后，会提示输入年份、月份和日期，然后输出对应的月相亮面变化。例如，输入2022年1月1日，输出如下：

请输入年份：2022
请输入月份：1
请输入日期：1
                    *                    
                ********                
             **************             
           ******************           
          ********************          
         **********************         
         ************************        
        **************************      
        ****************************    
       ******************************   
       *******************************  
      ********************************* 
     ***********************************
     ***********************************
    *************************************
    *************************************
   ***************************************
   ****************************************
  *****************************************
 *******************************************
 *******************************************
*********************************************
          阴影区域

可以看到，程序根据输入日期计算出了月相亮面比例，然后绘制出了对应的月亮。