$ $ ()

2002年，Berson等^[1]发现了哺乳动物视网膜的第三类感光细胞（ipRGC），打破了人眼只存在视觉通路的观念，非视觉通道逐渐被大众熟悉。ipRGC能参与调节人体光生物效应，包括激素分泌、昼夜节律和瞳孔大小等，光不仅影响着人体生理节律，在视觉工作、主观情绪等方面也存在影响。因此，照明质量的评价应由原来单一视觉评价过渡到视觉和非视觉的双重评价。目前对非视觉的评估尚未形成统一标准，Gall等^[2]通过非视觉光谱响应曲线和明视觉光谱响应曲线，提出用昼夜节律因子 $ {C}_{AF} $ （Circadian Action Factor）表示非视觉效应的大小；Rea等^[3]提出用昼夜节律刺激值CS （Circadian Stimulus）和昼夜节律照度 $ {CL}_{A} $ （Circadian Light）描述非视觉效应的效果，此模型考虑了S型视锥细胞光谱灵敏度曲线以及人眼黑视素光谱灵敏度曲线的影响，昼夜节律刺激的阈值和饱和值分别为0.1和0.7，昼夜节律照度是光源在2856 K黑体辐射体标准化之下的辐射照度值。

由于LED具有亮度大、可调节性强等特点，它的使用越来越广泛^[4]。为了设计满足生活工作需要的照明系统，需要兼顾视觉参数（显色性指数、辐射发光效率、蓝光危害效率）和非视觉参数（昼夜节律因子、昼夜节律刺激值）。郑莉莉等^[5]通过RGB照明系统，研究了昼夜节律因子和相关色温受各通道电流的影响，给出了昼夜节律因子、相关色温与电流之间的数学模型，并模拟了自然光在不同时段的光谱，不过该研究只探讨了电流对昼夜节律因子的影响，未涉及视觉参数。Jarboe等^[6-7]评估了几种LED在办公室的使用效果，研究表明，白天时至少维持2 h超过0.3的CS值，有利于提高工作效率和警觉性，该研究着重分析了昼夜节律刺激参数，但缺少昼夜节律因子的分析。

文中以LED芯片发光材料中光子能级分布规律的光谱模型为基础，通过改变LED所处的电流、温度来模拟光谱，确定RGBY各通道光源的绝对光谱和电流、温度之间的关系。通过遗传算法进行视觉和非视觉参数的多目标优化，对辐射发光效率和昼夜节律因子取Pareto最优解集，设置照度为300 lx，综合考量了昼夜节律刺激、蓝光危害效率、显色指数等指标，设计出8组适宜的照明系统，并探究8组方案非视觉参数与温度的关系。

以往光谱拟合主要通过高斯函数或洛伦兹函数来实现，但拟合效果都不理想。高斯模型拟合尾部收敛较快的光谱效果不错，而洛伦兹模型更适合尾部收敛较慢的曲线拟合^[8]。蒋腾等^[9]通过双高斯模型^[10-11]探究电流对RGB参数的影响，随着驱动电流增加，红光光谱向长波方向移动，绿光和蓝光光谱产生蓝移，拟合参数在0.9842以上。尽管以上研究综合了高斯和洛伦兹模型的优势，但拟合的效果依旧欠佳。徐广强等^[12]采用光子在二维空间内联合态密度函数作为单色LED的光谱模型，实现了对目标光谱的匹配。以上研究分析了单一的非视觉参数和视觉参数的关系，未分析两种非视觉参数之间的关系。

光谱产生本质上是光源在不同波长时能量分布不同^[13]。基于LED芯片发光材料中光子能级分布规律，采用如下数学模型表征单色LED光谱：

$ S(\lambda )=\frac{A}{\left[1+\mathrm{exp}\left(-\dfrac{\lambda -{\lambda }_{c1}}{{\omega }_{1}}\right)\right]\left[1+\mathrm{exp}\left(\dfrac{\lambda -{\lambda }_{c2}}{{\omega }_{2}}\right)\right]}\text{} $ (1)

式中： $ S\left(\lambda \right) $表示光谱辐射强度； $ \lambda $表示波长； $ A、 {\lambda }_{c1}、 {\lambda }_{c2}、 {\omega }_{1}、{\omega }_{2} $均为拟合参数。

利用RGBY四色光源来探究光源绝对光谱和电流、温度的关系。为防止外界光源对实验的影响，实验在密闭无光环境下进行，系统构建框图如图1所示。

光源驱动系统由微控制器、恒温加热台、驱动器、直流稳压电源（12 V）、直流降压模块（3.5 V）、RGBY灯珠组成。通过驱动模块DD313来驱动电路，驱动电流为350 mA。电流控制主要通过改变PWM通道占空比来实现，在每个通道串联一个0.1 Ω的电阻，通过蓝牙调节当前灯珠占空比进而改变电流，用万用表测量电阻的电压，这样就获得LED各通道电流值。建立光谱的电流模型时，将电流设置为100 mA、200 mA和300 mA获取数据。恒温加热台V-0505能灵活调节LED温度，范围设置在30°~120°内，以15°为间隔递增。使用SPIC-200光谱照度计测得不同电流、温度下的光谱。RGBY器件示意图以及数据测试实物图分别如图2、图3所示。

利用1 stopt软件对RGBY绝对光谱分别拟合， 以红灯在30°、300 mA下为例，此时在635 nm处光功率最大， $ {P}_{max}\left(635\right) $为0.0552 ${\rm{W}}\cdot {\left({\rm{m}}^{2}\cdot {\rm{nm}}^{-1}\right)}^{-1}$，光谱见图4。

根据以上数据，文中建立了光谱与电流、温度的二阶方程。拟合参数 $ {A、\lambda }_{c1}、{\lambda }_{c2}、{\omega }_{1}、{\omega }_{2} $与电流I和温度T之间的关系见公式（2）~（6）：

$ $$\begin{array}{r}\left[\begin{array}{l}{A}_r\\ A_g\\ A_b\\ A_y\end{array}\right]=\left[\begin{array}{l}6\times {10}^{-6}3\times {10}^{-5}{-10}^{-5}\\ -3\times {10}^{-4}-8\times {10}^{-4}-5\times {10}^{-4}\\ -0.00440.0280-0.0172\\ -2\times {10}^{-9}-7\times {10}^{-8}-2\times {10}^{-7}\end{array}\right]\cdot \left[\begin{array}{l}{I}^2\\ T^2\\ I\cdot T\end{array}\right]+\left[\begin{array}{l}0.0034-0.0072\\ 0.3227-0.0467\\ 4.9922-5.1855\\ 8\times {10}^{-5}{10}^{-5}\end{array}\right]\cdot \left[\begin{array}{l}I\\ T\end{array}\right]+\left[\begin{array}{l}0.4296\\ 7.7990\\ 135.5038\\ -0.0005\end{array}\right]\end{array}$$$ (2)

$ $$\begin{array}{r}\left[\begin{array}{l}{\lambda }_{c1}r\\ {\lambda }_{c1}g\\ {\lambda }_{c1}b\\ {\lambda }_{c1}y\end{array}\right]=\left[\begin{array}{l}-4\times {10}^{-5}7\times {10}^{-4}{10}^{-4}\\ 2\times {10}^{-5}-4\times {10}^{-4}2\times {10}^{-4}\\ -6\times {10}^{-5}{10}^{-4}6\times {10}^{-5}\\ 9\times {10}^{-6}4\times {10}^{-4}-3\times {10}^{-6}\end{array}\right]\cdot \left[\begin{array}{l}{I}^2\\ T^2\\ I\cdot T\end{array}\right]+\left[\begin{array}{l}0.03420.1280\\ 0.00880.3133\\ 0.03030.0314\\ 0.00110.0116\end{array}\right]\cdot \left[\begin{array}{l}I\\ T\end{array}\right]+\left[\begin{array}{l}650.4724\\ 575.4355\\ 448.0102\\ 514.9526\end{array}\right]\end{array}$$$ (3)

$ $$\begin{array}{r}\left[\begin{array}{l}{\lambda }_{c2}r\\ {\lambda }_{c2}g\\ {\lambda }_{c2}b\\ {\lambda }_{c2}y\end{array}\right]=\left[\begin{array}{l}-4\times {10}^{-5}3\times {10}^{-5}-7\times {10}^{-5}\\ 5\times {10}^{-6}9\times {10}^{-5}{10}^{-5}\\ 7\times {10}^{-6}2\times {10}^{-5}5\times {10}^{-6}\\ {10}^{-5}8\times {10}^{-6}-2\times {10}^{-5}\end{array}\right]\cdot \left[\begin{array}{l}{I}^2\\ T^2\\ I\cdot T\end{array}\right]+\left[\begin{array}{l}0.03090.1566\\ 0.00750.0548\\ 0.01640.0725\\ 0.00120.0569\end{array}\right]\cdot \left[\begin{array}{l}I\\ T\end{array}\right]+\left[\begin{array}{l}626.6856\\ 505.4527\\ 329.3168\\ 606.1363\end{array}\right]\end{array}$$ $ (4)

$ $$\begin{array}{r}\left[\begin{array}{l}{\omega }_{1}r\\ {\omega }_{1}g\\ {\omega }_{1}b\\ {\omega }_{1}y\end{array}\right]=\left[\begin{array}{l}{-10}^{-5}5\times {10}^{-5}-5\times {10}^{-6}\\ 7\times {10}^{-6}2\times {10}^{-5}5\times {10}^{-6}\\ 2\times {10}^{-7}-8\times {10}^{-6}4\times {10}^{-6}\\ 7\times {10}^{-6}4\times {10}^{-5}2\times {10}^{-5}\end{array}\right]\cdot \left[\begin{array}{l}{I}^2\\ T^2\\ I\cdot T\end{array}\right]+\left[\begin{array}{l}0.00700.0266\\ 0.00110.0311\\ 0.00090.0096\\ 0.00020.0257\end{array}\right]\cdot \left[\begin{array}{l}I\\ T\end{array}\right]+\left[\begin{array}{l}7.8166\\ 9.1606\\ 4.2916\\ 8.3980\end{array}\right]\end{array}$$ $ (5)

$ $$\begin{array}{r}\left[\begin{array}{l}{\omega }_{2}r\\ {\omega }_{2}g\\ {\omega }_{2}b\\ {\omega }_{2}y\end{array}\right]=\left[\begin{array}{l}-4\times {10}^{-6}{-10}^{-5}{-10}^{-5}\\ 3\times {10}^{-6}-2\times {10}^{-6}-2\times {10}^{-6}\\ {10}^{-6}-2\times {10}^{-5}-{10}^{-6}\\ {10}^{-5}5\times {10}^{-5}{10}^{-5}\end{array}\right]\cdot \left[\begin{array}{l}{I}^2\\ T^2\\ I\cdot T\end{array}\right]+\left[\begin{array}{l}0.00260.0103\\ 0.00040.0141\\ 0.00110.0118\\ 0.00040.0418\end{array}\right]\cdot \left[\begin{array}{l}I\\ T\end{array}\right]+\left[\begin{array}{l}3.4636\\ 6.0685\\ 10.1147\\ 37.0091\end{array}\right]\end{array}$$ $ (6)

以上拟合参数R-square均大于0.95，拟合效果较好。根据公式(6)可知，当给出单色LED所处的电流、温度值时，可以计算出公式(1)中的拟合参数，以红光为例，其电流和温度的综合影响模型如公式(7)所示，同理可以得出绿光、蓝光和黄光的光谱功率分布。这样就可以计算出当前光源系统的光谱功率分布，进而得出光源的视觉和非视觉参数。

$ $$\begin{array}{r}\{\begin{array}{l}{f}_{1r}=1+\exp (-{\displaystyle \frac{\lambda -[-4\times {10}^{-5}7\times {10}^{-4}{10}^{-4}]\cdot \left[\begin{array}{c}{I}^2\\ T^2\\ I\cdot T\end{array}\right]-\left[0.03420.1280\right]\cdot \left[\begin{array}{c}I\\ T\end{array}\right]-650.4724}{[-{10}^{-5}5\times {10}^{-5}-5\times {10}^{-6}]\cdot \left[\begin{array}{c}{I}^2\\ T^2\\ I\cdot T\end{array}\right]+\left[0.00700.0266\right]\cdot \left[\begin{array}{c}I\\ T\end{array}\right]+7.8166}})\\ A_r=[6\times {10}^{-6}3\times {10}^{-5}-{10}^{-5}]\cdot \left[\begin{array}{c}{I}^2\\ T^2\\ I\cdot T\end{array}\right]+\left[0.0034-0.0072\right]\cdot \left[\begin{array}{c}I\\ T\end{array}\right]+0.4296\\ f_{2r}=1+\exp \left({\displaystyle \frac{\lambda -[-4\times {10}^{-5}3\times {10}^{-5}-7\times {10}^{-5}]\cdot \left[\begin{array}{c}{I}^2\\ T^2\\ I\cdot T\end{array}\right]-\left[0.03090.1566\right]\cdot \left[\begin{array}{c}I\\ T\end{array}\right]-626.6856}{[-4\times {10}^{-6}-{10}^{-5}-{10}^{-5}]\cdot \left[\begin{array}{c}{I}^2\\ T^2\\ I\cdot T\end{array}\right]+\left[0.00260.0103\right]\cdot \left[\begin{array}{c}I\\ T\end{array}\right]+3.4636}}\right)\\ S_r(\lambda )={\displaystyle \frac{A_r}{f_{1r}\cdot f_{2r}}}\end{array}\end{array}$$ $ (7)

视觉参数包括^[14-15]辐射发光效率LER、蓝光危害效率 $ {\eta }_{B} $、显色指数 $ {R}_{a} $、相关色温CCT和照度E，非视觉参数包括昼夜节律因子 $ {C}_{AF} $、昼夜节律刺激CS。昼夜节律因子 $ {C}_{AF} $的计算见公式（8）：

$ {C}_{AF}=\frac{{\displaystyle {\int }_{380}^{780}P\left(\lambda ,I,T\right)\cdot C\left(\lambda \right)\text{d}\lambda }}{{\displaystyle {\int }_{380}^{780}P\left(\lambda ,I,T\right)\cdot V\left(\lambda \right)\text{d}\lambda }}$ (8)

式中： $ P（\lambda ,I,T） $为光源在电流I、温度T下光谱功率的叠加； $ C(\lambda ) $为非视觉光谱响应曲线； $ V(\lambda ) $为明视觉光谱响应曲线。昼夜节律刺激CS值可由公式（9）来计算：

$$$\begin{array}{rl}& C{L}_A=\{\begin{array}{c}1548[\int M_{C\lambda }{E}_{\lambda }\mathrm{d}\lambda +(a_{b-y}(\int \frac{S_{\lambda }}{m{p}_{\lambda }}{E}_{\lambda }\mathrm{d}\lambda -k\int \frac{V_{\lambda }}{m{p}_{\lambda }}\cdot \\ E_{\lambda }\mathrm{d}\lambda )-a_{rod}(1-\mathrm{e}{\displaystyle \frac{-{\displaystyle \int {V{}^{\prime }}_{\lambda }{E}_{\lambda }\mathrm{d}\lambda }}{RodSat}}))]l\\ \mathrm{i}\mathrm{f}\int \frac{S_{\lambda }}{m{p}_{\lambda }}{E}_{\lambda }\mathrm{d}\lambda -k\int \frac{V_{\lambda }}{m{p}_{\lambda }}{E}_{\lambda }\mathrm{d}\lambda ⩾0\\ 1548\int M_{C\lambda }{E}_{\lambda }\mathrm{d}\lambda \mathrm{i}\mathrm{f}\int \frac{S_{\lambda }}{m{p}_{\lambda }}{E}_{\lambda }\mathrm{d}\lambda -k\int \frac{V_{\lambda }}{m{p}_{\lambda }}{E}_{\lambda }\mathrm{d}\lambda <0\end{array}\\ & CS=0.7-{\displaystyle \frac{0.7}{1+{\left({\displaystyle \frac{C{L}_A}{355.7}}\right)}^{1.1026}}}\end{array}$$ $ (9)

式中： $ C{L_A} $为生理节律光照，1548是归一化常数； $ {E_\lambda } $为光谱辐照度分布； $ {S_\lambda } $为S型视锥细胞光谱灵敏度曲线； $ m{p_\lambda } $是黄斑色素透射率曲线； $ {V_\lambda } $和 $ {V_\lambda }^, $分别代表明视觉以及暗视觉的光谱响应曲线； $RodS at、 k、 {a}_{b-y}$和 $ {a}_{rod} $为常数，分别为6.5 W/ ${\rm{m}}^{2}$、0.2616、0.7和3.3。

辐射发光效率LER和照度E的计算如下：

$ LER=\frac{{K}_{m}{\displaystyle {\int }_{380}^{780}P(\lambda ,I,T)V(\lambda )\text{d}\lambda }}{{\displaystyle {\int }_{380}^{780}P(\lambda ,I,T)\text{d}\lambda }}$ (10)

$ E={K}_{m}{\displaystyle {\int }_{380}^{780}P(\lambda ,I,T)V(\lambda )\text{d}\lambda }$ (11)

$ {\mathrm{式}\mathrm{中}：K}_{m} $表示光源最大光通量，为683 lm/W。

为了设计既有较高 $ {C_{AF}} $又有较高LER的照明系统，需对两参数折衷处理。根据国际照明委员会CIE对办公室环境的参数要求，照度设为300 lx，显色指数 $ {R}_{a} $大于90。考虑到现实环境因素会阻碍CS传递到人眼，将CS值设置为0.3~0.4之间，蓝光危害效率 $ {\eta }_{B} $控制在0.3以下，在2500~8000 K色温下寻求满足上述条件的照明方案。采用遗传算法对 $ {C_{AF}} $和LER进行多目标优化，为了消除量纲的影响，把 $ {C_{AF}} $和LER归一化处理， 通过Matlab仿真得出 $ {C_{AF}} $在限制下最大值和最小值分别为1.1647和0.1963；LER最大值和最小值分别为368 lm/W和265 lm/W。优化的方程如下，以四路占空比大小作为未知量 $ {D_r},{D_g},{D_b},{D_y} $:

$ \left\{$$\begin{array}{rl}& f(1)=Max\left(\frac{C_{AF}(D_r,D_g,D_b,D_y)-0.1963}{1.1647-0.1963}\right)\\ & f(2)=Max\left(\frac{LER(D_r,D_g,D_b,D_y)-265}{368-265}\right)\\ & \mathrm{s}.\mathrm{t}.{\eta }_B(D_r,D_g,D_b,D_y)⩽0.3;\\ & R_a(D_r,D_g,D_b,D_y)⩾90;\\ & 0.3⩽\text{CS}(D_r,D_g,D_b,D_y)⩽0.4;\\ & E(D_r,D_g,D_b,D_y)=300\end{array}$$\right. $ (12)

求 $ {C_{AF}} $和LER的Pareto解集时，最优前端个体系数ParetoFraction取0.08，种群大小PopulationSize取100，最大遗传代数Generations设置为500代，停止代数StallGenLimit为500代，适应度函数偏差TolFun设为 $ 1{e}^{-100} $，得到的8个照明系统方案如表1所示。

由表1可知，照度为300 lx时，照明方案的色温分布在2855~7628 K，节律刺激值CS随色温的变化曲线如图5所示， CS在2800~3500 K、4000~8000 K时随色温增加而变大，在3500~4000 K呈现急剧下降的趋势，这是因为 $ {CL}_{A} $和色温有公式（9）的关系，而CS和 $ {CL}_{A} $呈正相关。要想维持CS值在0.3~0.4之间，色温的可选空间较大。

为验证该光谱模型的准确性，测试这8种方案下各路电流值（见表2），再改变灯珠温度，将各方案实测值和模型值进行比较。8种方案是按照色温的大小来排列的，即按照昼夜节律因子 $ {C_{AF}} $的大小排列，方案8的色温最高，即方案8对人体的非视觉效应影响最大。由于是在2500~8000 K色温内的寻优，所以这8种方案可以代表限制条件下的最佳方案。

由于目前常用的照明光源为D50和D65，而方案5和方案7的色温分别最接近D50和D65，以方案5和方案7为例，分析实测光谱和模型光谱的拟合情况，如图6所示。

从拟合效果看，两个方案均能准确反映实际光谱。为了更精确地体现模型和实测光谱的差距，表3计算了8种方案各参数指标的误差，并标明了最大误差时的温度。由于SPIC-200可以直接测出相关色温CCT、显色指数 $ {R}_{a} $、照度E，但辐射发光效率LER、昼夜节律因子 $ {C}_{AF} $等参数无法直接测出，对于无法直接测得的参数，通过SPIC-200测量出实际的光谱，再根据公式（10）计算得出。

表3中，模型和实测参数的绝对误差最大为7.06%，说明该模型能较好拟合实际光谱。初始照度为300 lx，但随着LED温度的升高，E逐渐衰减。图7给出了照度E和温度T的拟合曲线，两者呈线性负相关， $ {R}^{2} $为0.9869。不同方案参数的最大误差值有差别，影响参数误差的因素为LED的结温差、热耦合效应等。驱动LED芯片时，连续波电流会使LED因热效应进而结温上升，影响实验结果，在实验中采用脉冲电流驱动LED发光，并且仅在测量光谱的一瞬间闭合开关，最大限度减小LED的结温差对实验数据的影响。由于四路LED芯片之间最小的距离为19 mm左右，热耦合效应很微弱，对实验影响极小。

整个实验的温度调节是通过恒温加热平台V-0505实现的，尽管将温度从30°~120°进行了等分，但温度随时间的变化并非是均匀的，因此需要考虑梯度温度的影响。该实验的温度仅受时间影响，所以在一维状态下，梯度温度可看成是温度的变化率。经过测量可知，温度从30°升到120°需要18 min，以1 min为间隔，通过红外线测温仪来测量实验平台的实时温度，绘制出温度随时间的变化曲线见图8。

由图8可以看出，加热过程中，温度的变化率逐渐衰减，12 min之后，温度几乎平稳。说明在刚开始时，光源的非视觉参数变化速率较快，随着温度的逐渐平稳，非视觉参数也增加到最大，趋于平稳。

以方案1、5、8为例，两种非视觉参数随温度的变化见图9，可以看出：温度升高， $ {C}_{AF} $值逐渐增加，CS值逐渐减小。当从30 ℃上升到120 ℃时，方案1的 $ {C}_{AF} $从0.2223增至0.2459，CS从0.3008减至0.2597；方案5的 $ {C}_{AF} $由0.5783上升到0.6701，CS从0.3344减少到0.2998；方案8的 $ {C}_{AF} $值由1.0192上升到1.0902，CS从0.3986减少到0.3657。

随着温度变化，蓝光占比逐渐增加，非视觉效应增强，但 $ {C}_{AF} $和CS的计算受照度影响不同，说明二者在描述非视觉效应时存在一些差异。为消除照度的影响，对光谱进行照度补偿，使照度恒定在300 lx，通过电流、温度模型计算当前光谱分布，补偿公式如下：

$ \left\{$$\begin{array}{rl}& f=Min[k_1\cdot {\displaystyle \sum _{380}^{780}{(P_{I,T}(\lambda )-P_{I_0{\text{,T}}_0}(\lambda ))}^2}+k_2\cdot {(C_{{AF}_{I,T}}-C_{{AF}_{I_0{,T}_0}})}^2+\\ & k_3\cdot {({\text{CS}}_{I,T}-{\text{CS}}_{I_{0,}{T}_0})}^2+k_4\cdot {(x_{I,T}-x_{I_0{,T}_0})}^2+k_5\cdot {(y_{I,T}-y_{I_0{,T}_0})}^2]\\ & \mathrm{s}.\mathrm{t}.D_i\in [0,1]i=r,g,b,y\\ & E(D_r,D_g,D_b,D_y)=300\mathrm{l}\mathrm{x}\end{array}$$\right. $ (13)

式中： ${D_{{r}}},{D_{{g}}},{D_{{b}}},{D_{{y}}}$为温度T下补偿占空比； ${P_{{I,T}}}{\left( \lambda \right)}$、 ${P_{{{{I}}_0},{{{T}}_0}}}{\left( \lambda \right) }$指电流I、温度T和初始电流I₀；温度T₀下光源的光谱； ${C_{{{{A}}{{F}}}_{{I,T}}}}$、 ${C_{{{{A}}{{F}}}_{{{I}_0,{{T}}_0}}}}$为电流I、温度T和初始电流I₀，温度T₀下昼夜节律因子； ${\text{CS}_{{I,T}}}$、 ${\text{CS}_{{{{{I}}_0}},{{{T}_0}}}}$为电流I、温度T和初始电流 ${{{I}}_0}$，温度 ${{{T}}_0}$下昼夜节律刺激； $({x_{{I,T}}},{y_{{I,T}}})$、 $({x_{{{{I}}_0}{\text{,}}{{{T}}_0}}},{y_{{{{I}}_0},{{{T}}_0}}})$为电流I、温度T和初始电流I₀，温度T₀下的色坐标； $ {k_1},{k_2},{k_3},{k_4},{k_5} $为各项的权重。昼夜节律因子 $ {C}_{AF} $和昼夜节律刺激CS作为两种衡量光的非视觉效应的变量，理应随温度的变化呈现相同的趋势。但在补偿前发现二者随温度的变化趋势截然相反，是因温度升高后有物理量的变化，导致二者趋势不同。经分析发现照度受温度的影响很大，于是决定用PWM调光补偿照度，其采用高频率的周期数字脉冲驱动光源，通过改变周期脉冲内有效电平脉宽的宽度来调整输出的有效电流，从而调节光源参数。PWM调光补偿的响应速度更快，光源控制精度更高，可控性更强，是照明调光领域中的优选方式。补偿后照度始终维持在300 lx，其余的光源参数尽可能地接近补偿前，从而得到合理的补偿方案。补偿后的非视觉参数随温度的变化趋势如图10所示。

补偿后，随着温度升高，C_AF的增大趋势略微减缓，CS由下降变为上升。方案1的C_AF从0.2223逐渐增加到0.2405，CS的值从0.3008增到0.3204；方案5的C_AF从0.5783增加到0.6522，CS的值从0.3344增加到0.3637；方案8的C_AF从0.5783逐渐增加到0.6522，CS的值从0.3344增加到0.3637。C_AF和CS均随温度升高而增大，且二者表现出一定的正相关性。

文中基于LED芯片发光材料中光子能级分布的规律，提出了光谱的电流、温度模型。通过该模型可以根据LED当前的电流、温度参数，模拟光源的实际光谱。在综合考虑非视觉参数和视觉参数的条件下，设计出8种符合照明要求的照明方案，并针对各个方案验证了电流、温度模型的可靠性。同时比较各方案下，昼夜节律因子和昼夜节律刺激随温度的变化趋势，发现昼夜节律因子随温度的升高而增大，但昼夜节律刺激随温度增加而减小，照度随温度的增大而线性衰减。补偿照度后，两种非视觉参数的值均随温度的升高而增大，且二者表现出一定的正相关性，这表明随着温度的升高，人体的非视觉效应的确在逐渐增强。该研究为LED照明光源设计中非视觉效应的考量提供了一定参考依据。