帮我生成一段代码。该代码能随机生成一段脉动速度随时间的变化函数，周期在4到20秒之间，随机分成2-5个时间段，每个时间段不需要相等。在每一段中，速度是随机分配的线性、二次或正弦函数，每一段的结尾速度和下一段的开始速度相同。速度范围在1到4 m/s之间。

import random
import numpy as np
import matplotlib.pyplot as plt

# 随机生成一个周期
period = random.uniform(4, 20)
# 随机生成时间段数量
num_segments = random.randint(2, 5)
# 随机生成每个时间段的长度
segment_lengths = [period/num_segments]*num_segments
random.shuffle(segment_lengths)
# 随机生成每个时间段的函数类型
segment_types = random.choices(["linear", "quadratic", "sine"], k=num_segments)
# 随机生成每个时间段的速度变化范围
segment_speed_ranges = [random.uniform(1, 4) for _ in range(num_segments)]

# 初始速度
speed = random.uniform(1, 4)

# 定义 x 轴的时间点
time_points = np.linspace(0, period, 1000)
speeds = []

# 遍历每个时间段
for i in range(num_segments):
    segment_type = segment_types[i]
    segment_speed_range = segment_speed_ranges[i]
    segment_length = segment_lengths[i]
    x_points = np.linspace(0, segment_length, 1000)
    if segment_type == "linear":
        # 速度上升或下降
        if random.random() < 0.5:
            m = segment_speed_range / segment_length
            b = speed
            y_points = m*x_points + b
        else:
            m = -segment_speed_range / segment_length
            b = speed + segment_speed_range
            y_points = m*x_points + b
    elif segment_type == "quadratic":
        # 速度上升或下降
        if random.random() < 0.5:
            a = segment_speed_range / (segment_length**2)
            b = 0
            c = speed
            y_points = a*x_points**2 + b*x_points + c
        else:
            a = -segment_speed_range / (segment_length**2)
            b = 0
            c = speed + segment_speed_range
            y_points = a*x_points**2 + b*x_points + c
    else:
        # 正弦函数
        amplitude = segment_speed_range / 2
        angular_frequency = 2 * np.pi / segment_length
        phase = random.uniform(0, 2*np.pi)
        y_points = amplitude