科学，这个看似高深莫测的领域，其实充满了趣味和惊喜。今天，就让我们一起来揭开那些神奇现象的神秘面纱，感受科学的魅力吧！

神奇现象一：彩虹

提到彩虹，大家都会想到它那绚丽的七彩光芒。其实，彩虹的形成原理非常简单。当太阳光穿过雨滴时，光线会发生折射、反射和色散，最终形成我们看到的彩虹。这个过程可以用以下代码来模拟：

import matplotlib.pyplot as plt
import numpy as np

# 定义光的折射和反射角度
def refract_angle(angle):
    return angle - 42.4  # 42.4是光的折射率

# 模拟光线在雨滴中的传播
def simulate_light_path(angle):
    angles = [angle]
    for _ in range(20):  # 模拟光线在雨滴中传播20次
        angle = refract_angle(angle)
        angles.append(angle)
    return angles

# 生成彩虹图
def generate_rainbow():
    angles = [i for i in range(-180, 180, 1)]
    colors = ['red', 'orange', 'yellow', 'green', 'blue', 'indigo', 'violet']
    angles = [simulate_light_path(angle) for angle in angles]
    fig, ax = plt.subplots()
    for angle, color in zip(angles, colors):
        ax.plot(angle, [0]*len(angle), color=color)
    ax.set_xlim(-180, 180)
    ax.set_ylim(-1, 1)
    plt.show()

generate_rainbow()

神奇现象二：海市蜃楼

海市蜃楼是一种自然现象，当光线穿过不同密度的空气层时，会发生折射，形成虚像。这个现象可以用以下代码来模拟：

import matplotlib.pyplot as plt
import numpy as np

# 定义折射率函数
def refractive_index(density):
    return 1 + 0.00025 * density

# 模拟光线在空气中的传播
def simulate_light_path(density):
    angles = [0]
    for _ in range(100):  # 模拟光线在空气中传播100次
        density += 0.1  # 假设密度每传播1单位距离增加0.1
        angle = np.arcsin(np.sin(angles[-1]) / refractive_index(density))
        angles.append(angle)
    return angles

# 生成海市蜃楼图
def generate_halo():
    densities = np.linspace(0, 100, 1000)
    angles = [simulate_light_path(density) for density in densities]
    fig, ax = plt.subplots()
    for angle in angles:
        ax.plot(angle, [0]*len(angle), color='blue')
    ax.set_xlim(-180, 180)
    ax.set_ylim(-1, 1)
    plt.show()

generate_halo()

神奇现象三：雪花

雪花是自然界中最美丽的晶体之一。雪花之所以形态各异，是因为水分子在结冰过程中会按照一定的规律排列。这个现象可以用以下代码来模拟：

import numpy as np
import matplotlib.pyplot as plt

# 定义雪花晶体的生长过程
def grow_snowflake():
    positions = []
    for _ in range(1000):  # 生长1000个雪花晶体
        position = np.random.rand(2) * 2 - 1  # 在[-1, 1]范围内随机生成位置
        positions.append(position)
        # 根据晶体的规则生长
        for i in range(len(positions) - 1):
            dx = positions[-1][0] - positions[i][0]
            dy = positions[-1][1] - positions[i][1]
            distance = np.sqrt(dx**2 + dy**2)
            if distance < 0.1:
                positions[i] = positions[i] + (dx / distance, dy / distance)
    return positions

# 生成雪花图
def generate_snowflake():
    positions = grow_snowflake()
    fig, ax = plt.subplots()
    ax.scatter(*zip(*positions), c='blue', s=10)
    ax.set_xlim(-1.5, 1.5)
    ax.set_ylim(-1.5, 1.5)
    plt.show()

generate_snowflake()

通过这些有趣的例子，我们可以更加直观地了解科学现象的原理。希望这篇文章能让你在轻松愉快的氛围中，对科学产生更深的兴趣！