0

我想计算朝北、朝南等方向的墙壁上的太阳辐射。我有该位置的太阳辐射,我尝试了以下等式:

S_module = S_incident * (math.cos(alpha)*math.sin(beta)*math.cos(psi-theta) + math.sin(alpha)*math.cos(beta))

alpha = 24.86 #sun elevation angle 
theta = 81.58 #sun azimuth angle
beta = 90 #wall tilt angle (vertical module beta=90)
psi = 135 #surface orientation: psi=0 == facing north, psi=180 == facing south; is the azimuth angle that the module faces 
S_incident = 663 [W]

这是 2019-6-21 6h timezone=utc

我尝试了 python 包 pysolar,上面的等式在这里找到:https ://www.pveducation.org/pvcdrom/properties-of-sunlight/arbitrary-orientation-and-tilt

我的问题是我没有得到正确的结果,例如对于上述属性,结果是-495 [W]。这不可能是真的,因为在早晨,辐射必须至少有一个小的正值。

感谢您的任何提示!

4

1 回答 1

1

阅读文档:三角函数以弧度工作;你喂它度数。你的 90 弧度墙与现实不符。:-)

import math

alpha = 24.86 #sun elevation angle 
theta = 81.58 #sun azimuth angle
beta = 90 #wall tilt angle (vertical module beta=90)
psi = 135 #surface orientation: psi=0 == facing north, psi=180 == facing south; is the azimuth angle that the module faces 
S_incident = 663 # Watts

S_module = S_incident * (math.cos(alpha)*math.sin(beta)*math.cos(psi-theta) +
                         math.sin(alpha)*math.cos(beta))

print(math.cos(alpha), math.sin(alpha))
print(math.cos(beta), math.sin(beta))

输出:

0.9630361043608977 -0.26937234768510715
-0.4480736161291701 0.8939966636005579

这些显然不是您期望的值。

用因子转换math.pi / 180

于 2019-06-24T17:09:02.743 回答