定义一个函数 ContrastStretch 实现该函数s=1/(1+math.pow(m/r,E))的计算，输入向量r，给出参数m=0.5,输出向量s。编写测试代码，计算该函数在不同E值下，区间0到1的1000个采样点的输出，并绘制相应的函数曲线。

首先，我们来定义一个名为`ContrastStretch`的Python函数，它接受一个向量`r`作为输入，并应用给定的公式来计算输出向量`s`。我们将使用numpy库来进行数学运算。

import numpy as np

def ContrastStretch(r, m=0.5, E=1):
    s = 1 / (1 + np.power((r - m), E))
    return s

接下来，我们将编写测试代码，生成区间0到1之间的1000个等间距采样点，并分别对不同的`E`值应用这个函数，最后将结果绘制成图表。

# 生成1000个0到1之间的均匀采样点
x = np.linspace(0, 1, 1000)

# 测试不同E值
Es = [0.5, 1, 2]  # 可选的E值列表

fig, axs = plt.subplots(len(Es), sharex=True)
for i, E in enumerate(Es):
    s_values = ContrastStretch(x, m=0.5, E=E)
    axs[i].plot(x, s_values, label=f'E={E}')
    axs[i].set_title(f'Contrast Stretch for E = {E}')
    axs[i].legend()

plt.xlabel('Input (r)')
plt.ylabel('Output (s)')
plt.show()

上述代码首先创建了一个包含三个不同E值的子图，然后对于每个子图，计算对应的`s`值并绘制出来。这会生成多条曲线，每一条代表一种E值下的对比度拉伸效果。