src/d_spektra.py


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20

import numpy as np
from shared import *
import crl, crl.tables

plot_setup(width=0.7)

for T in reversed([5000, 5500, 6500, 7500]):
	wvl = np.linspace(380, 780, 100)
	n = crl.CIEDaylightSpectrum(T).Yi(wvl)

	plt.plot(wvl, n,
	         label="$D_{%d}$" % (T // 100))

plt.xlabel("$\\lambda$ [nm]")
plt.ylabel("$n_\\lambda(\\lambda)$")
plt.xlim(380, 780)
plt.ylim(0, 200)
plt.legend()
plot_grid()
plot_export()