图不看它应该如何Python

Question

所以我试图根据一年中的某一天重建这一天的图表：

该图表遵循公式

该功能取决于一年中的某天t和纬度的小theta。其他信息给出：大THETA是地球轴线的倾斜度为23.4°，T_24是当天的24小时长度，T_365当年的长度（以天为单位）是365 * 24。 t0是春天的第一天，即79天* 24小时。

我们单独定义chi

因为这种依赖

到目前为止，这是我的代码：

import matplotlib.pyplot as plt
import numpy as np
from scipy import *

THETA=np.radians(23.4)
T24=24
T365=365*24
t0=79*24

ts = np.linspace(0, 365, 100)
thetas = np.linspace(-90, 90, 100)
t, theta = np.meshgrid(ts, thetas)

z=T24*(1-arccos(tan(np.radians(theta))*sin(THETA)*sin(2*pi*(t-t0)/T365)/sqrt(1-(sin(THETA))**2*(sin(2*pi*(t-t0)/T365))**2))/pi)

CS = plt.contour(t, theta, z)

我仍然没有包括chi的区别，因为我甚至更早遇到了问题。所以，问题是图表看起来不应该如何，我只得到7条水平线。

Answer 1

似乎等式中的时间单位是小时。因此ts必须以小时为单位：

ts = np.linspace(0, 365, 100)*24