一种多色led实现光源最佳显色性的光谱优化方法
【专利摘要】本发明公开了一种多色LED实现光源最佳显色性的光谱优化方法。在本方法中，首先，通过采样和测量获取各色单颗LED的光谱样本数据；然后，将目标光源色品和照度数据进行转换以获取其绝对三刺激值；最后，采用一种约束型差分进化算法迭代优化获得对应LED的驱动值，同时通过色温误差百分比、色品差、照度误差百分比以及显色性评价指标对优化测量结果进行评价。本发明解决了当前光谱优化方法存在的不足，能够同时满足实际应用对光源色品、照度以及显色性等照明品质方面的需求，也为实现光源色温和照度的连续可调提供了一种切实有效的方法，对促进智能照明相关技术的发展具有十分重要的意义。
【专利说明】
一种多色LED实现光源最佳显色性的光谱优化方法
技术领域
[0001] 本发明涉及LED照明领域，尤其涉及一种多色LED实现光源最佳显色性的光谱优化 方法。
【背景技术】
[0002] LED光源显著的节能、长寿等特性越来越引起人们的广泛关注。当前主要有两种方 法用于产生照明白光:一种为蓝光LED激发黄色荧光粉，即荧光粉转换型白光LED，这种白光 LED由于存在Stokes漂移故效率较低，且难以获得较高的显色性;另一种为借助LED丰富的 颜色分布及色光混合原理组合多色LED产生白光，该方式通过调整各色LED的混合比例或强 度能够产生想要的白光，而这一特性也为实现光源光谱优化提供了无限可能。显色性作为 评价照明光源最为重要的品质参数，在一般照明应用中，如办公、家居、商场、博物馆以及产 品颜色检测等场合，往往均对其有特定的要求，因此以光源最佳显色性为目标优化多色LED 混合光谱具有十分重要的实际意义。
[0003] 目前，采用多色LED对其混合光谱进行优化多集中于通过调整LED的峰值波长和带 宽来实现特定显色性评价指标或辐射效率的最大化，这种优化方式旨在为促进新型LED的 发展和LED类型的选择提供指导，并不是对实际已有的LED进行优化。另外，已有的优化方法 也不能很好地针对特定的光源色品和照度来优化多色LED的混合光谱。为此，本发明公开了 一种多色LED实现光源最佳显色性的光谱优化方法，该方法能够在匹配给定光源色品和照 度的条件下实现其显色性的最大化。
【发明内容】
[0004] 为了克服已有光谱优化方法的不足，本发明提供了一种多色LED实现光源最佳显 色性的光谱优化方法。
[0005] 为了实现上述目的，本发明采用的技术方案如下：
[0006] -种多色LED实现光源最佳显色性的光谱优化方法，包括如下步骤：
[0007] SI: LED光谱样本数据的获取；
[0008] S2:光源目标色品和照度数据的输入及处理，得到其对应的目标绝对三刺激值；
[0009] S3:采用约束型差分进化算法进行迭代优化，获得匹配光源目标绝对三刺激值且 实现其最佳显色性的各类LED的驱动值；
[0010] S4:对上述驱动值下测量得到的光谱进行优化评价。
[0011] 作为优选，所述步骤S1中LED光谱样本数据的获取过程具体包括以下步骤：
[0012] S101:在各类LED的驱动值动态范围内以一定采样间隔建立LED驱动数据样本，记 为d^，表示第i类LED的第j个驱动数据样本，用行向量将每类LED的驱动数据样本表示为cU = [(^▲,2,. . .，di,max]，其中di,max表示第i类LED的最大采样驱动值，并将各类LED最大采 样驱动值构成的行向量记为0 111￡?=[(11,111￡?，(12,111 ￡?，...，(111,111￡?]，其中11为]^1)的种类数；
[0013] S102:将步骤S101中的驱动数据样本分别依次驱动点亮各类单颗LED，并等待其发 光稳定；
[0014] S103:在各类单颗LED垂直正下方放置标准板，所述的标准板为标准灰板或标准白 板，采用光谱辐射计依次测量获得各类单颗LED点亮稳定时经标准板反射的光谱辐射亮度 分布数据，记为S^(A)，表示第i类单颗LED在驱动值点亮稳定后经标准板反射的光谱辐 射亮度分布，单位为W/(sr · m2 · ηπι)，λ为可见光波长。
[0015]作为优选，所述步骤S2中光源目标色品和照度数据的输入及处理具体包括以下步 骤：
[0016] S201:输入光源的目标色品数据，其中目标色品数据可以为光源的目标相关色温 Τ、光源的目标相关色温Τ和目标Duv或者光源的任意一种目标色品坐标，然后根据色度学原 理将目标色品输入数据转换为目标CIE1931色品坐标(x t，yt)和目标CIE1976UCS色品坐标 (u/，ν/ );
[0017] S202:输入待匹配光源照度Et，单位为lx，并结合标准板的反射比，根据光度学和 色度学原理计算待匹配光源照度E t和步骤S201中目标色品坐标(xt，yt)对应的目标绝对三 刺激值(1，1，匕），标准板的反射比可通过分光光度计测量计算获得。
[0018] 作为优选，所述步骤S3中光谱优化采用约束型差分进化算法进行迭代优化，其具 体步骤为：
[0019] S301:根据光照度叠加原理计算多色LED光源中各类LED的等效数量，并用行向量 表示为Nmax = [Nmax, 1，Nmax, 2，…，Nmax, η]，其中η为LE：D的种类数；
[0020] S302:按照下式计算各类LED在采样驱动值均为du打开时的混合光谱辐射亮度分 布Si，j，M(X):
[0021] Si;j;M(A)=Nmax>iSi>j(A) (9)
[0022] 式中Nmax,i为步骤S301中第i类LED的等效数量，Si,j(A)为步骤S103中单颗LED在采 样驱动值为di,j时测量的光谱福射亮度分布；
[0023] S303:通过三次样条插值算法预测各类LED在其驱动值动态范围内任一驱动值时 的光谱福射亮度分布Si,M(X)，gp
[0024] Si，M(X) =spline(di，Si(A)，Di) (10)
[0025] 式中spline表示执行三次样条插值，di为步骤S101中第i类LED采样驱动值构成的 行向量为步骤S302中第i类LED在波长λ处不同采样驱动值的光谱辐射亮度构成的行 向量，Di为第i类LED在其驱动值动态范围内的任一驱动值；
[0026] S304:定义目标优化函数，将给定光源目标色品和照度参数的光谱优化转换为在 满足式(4)的约束条件下求解式(3)的最大值问题：
[0027]
[0028]
[0029]
[0030] 其中
[0031]
[0032] 式中R( ·)表示某种光源显色性评价指标或者光源显色性与光效的组合型评价指 标，如CIE-般显色指数Ra，且R( ·)为光源光谱的函数，Sm(X)表不多色LED的混合光谱，Si,M (λ)为第i类LED在其驱动值动态范围内任一驱动值处的光谱辐射亮度分布，可通过步骤 S303获得，D为待求解行向量，向量中各元素 Di表示第i类LED在其驱动值动态范围内任一驱 动值，yU). F(/l)为CIE1931标准色度观察者，η为LED的种类数；
[0033] S305:执行初始化操作，即在0~Dmax间随机生成Np个均匀分布的候选解向量集，记 为015,1=[01人1，02九 1，"_，01^,1]，其中1^表示第1^个候选解向量，且1^={1，2，...，邮}，恥取决 于LED的种类数n，一般Np取为10η或更大值，1代表初始候选解向量集，即第1次迭代时的候 选解向量集，同样将第G次迭代时的候选解向量集记为Dk.GWDuGAn-Ak.G];
[0034] S306:将步骤S305中Np个候选解向量Dk,c依次代入步骤S304中的式(3)和式(4);
[0035] S307:从Np个候选解向量Dk,G中随机选择s个候选解向量，记为Dp,g=[D 1>p,g， 02，口,。厂..，011,卩,。]，其中卩={1，2，...，8};
[0036] S308:根据差分进化算法原理对步骤S307中的s个候选解向量Dp,g依次执行变异和 交叉操作，经此操作共获得s个候选解向量UD P, c;
[0037] S309:将步骤S308中s个候选解向量UDP,G依次代入步骤S304中的式(3)和式(4);
[0038] S310:根据步骤S304中的式(3)和式(4)比较步骤S308中s个候选解向量UDP,G以及 步骤S307中s个候选解向量DP,c的计算结果，并用前者最优候选解向量取代后者中与此相比 较差的候选解向量DP, c，同样也对候选解向量集Dk, c中的候选解向量进行取代；
[0039] S311:若步骤S308中s个候选解向量UDP,G的所有最优候选解向量均不满足式(4)， 则将这些最优候选解向量存储，并每隔m次迭代用其随机取代候选解向量集Dk,c中相同数量 的候选解向量，同时确保候选解向量集Dk,c中的最优候选解向量不被取代；
[0040] S312:重复执行步骤S307~S311，直至满足设定的迭代终止条件，即最大迭代次数 Gmax，并将此时获得的最佳解向量记为Dt,end，且t= {1，2,. . .，N}，N为优化执行次数；
[0041] S313:重复执行步骤S305~S312,直至满足设定的优化执行次数N，并根据步骤 S304中的式(3)和式(4)，从N次优化获得的N个最佳解向量Dt, md中选出最终的最佳解向量， 并将其记为D。, -，即为优化获得的各类LED的驱动值；
[0042] S314:将步骤S313中优化获得的各类LED的驱动值输入到相应的控制软件点亮 LED，并通过光谱辐射计测量光谱辐射亮度分布Stm(A)、相关色温Tm以及CIE1976UCS色品坐 标(Un/ ,νη/ )，通过照度计测量相应的照度Em。
[0043] 作为优选，所述步骤S4中光谱优化评价具体包括以下步骤：
[0044] S401:采用色温误差百分比ετ评价优化测量出的相关色温与目标相关色温间的差 异，ετ的定义为
[0045]
[0046] 式中Τ为步骤S201中的目标相关色温，Tm为步骤S314中的相关色温测量值；
[0047] S402:采用CIE1976UCS色品差Δ w评价优化测量出的光源色品与目标色品间的 差异，Aw的定义为
[0048]
(15)
[0049] 式中（u/，ν/ )为步骤S201中的目标CIE1976UCS色品坐标，（u/，ν/ )为步骤S314 中实际测量的CIE1976UCS色品坐标；
[0050] S403:采用照度误差百分比εΕ对照度匹配结果进行评价，εΕ的定义为
[0051 ]
(16).
[0052] 式中Et为步骤S202中的待匹配照度，Em为步骤S314中的照度测量值；
[0053] S404:分别采用CIE-般显色指数特殊显色指数R9以及IES颜色保真性指数Rf 和色域面积指数心对步骤S314中优化测量出的光谱Stm(A)进行显色性评价。
[0054]本发明的有益效果是:本发明解决了当前采用LED实现光源最佳显色性的光谱优 化方法存在的不足，所述方法可在匹配特定光源色品和照度的条件下获得其显色性的最大 化，能够同时满足实际应用对光源色品、照度以及显色性等照明品质方面的需求，并为实现 光源色温和照度的连续可调提供了一种切实有效的方法。
【附图说明】
[0055] 图1是本发明在2800K、3500K、4000K和4500K等4种相关色温下获得的最佳显色性 的光谱辐射亮度分布；
[0056] 图2是本发明在5000Κ、5500Κ、6000Κ和6500Κ等4种相关色温下获得的最佳显色性 的光谱辐射亮度分布。
【具体实施方式】
[0057]本发明提供的多色LED实现光源最佳显色性的光谱优化方法，主要包括如下步骤： [0058] SI: LED光谱样本数据的获取；
[0059] S2:光源目标色品和照度数据的输入及处理，得到其对应的目标绝对三刺激值；
[0060] S3:采用约束型差分进化算法进行迭代优化，获得匹配光源目标绝对三刺激值且 实现其最佳显色性的各类LED的驱动值；
[0061 ] S4:对上述驱动值下测量得到的光谱进行优化评价。
[0062] 下面结合实施例和附图对上述方法进行详细阐述，以更好地理解本发明的实质。
[0063] 本实施例采用的多色LED光源是由17种类型一千多颗单色LED构成的一个大尺寸 LED矩阵，并通过基于PWM的数字调光方式对其进行调光控制。需要说明的是，本发明并不局 限于实施例中所采用的LED类型和调光方式，只要采用LED光源优化光源光谱均适用于本发 明。
[0064] 1 .LED光谱样本数据的获取过程具体包括以下步骤：
[0065] S101:在各类LED的驱动值动态范围内以8作为驱动值采样间隔建立LED驱动数据 样本，记为d^，表示第i类LED的第j个驱动数据样本，用行向量将每类LED的驱动数据样本 表示为. . .，di,max]，其中，di,max表示第i类LED的最大采样驱动值，驱动值的取 值范围为〇~255，并将各类1^0最大采样驱动值构成的行向量记为0_=[(1 1,_，(12,_，...， dn,max]，其中n= 17;
[0066] S102:将步骤S101中的驱动数据样本分别依次驱动点亮各类单颗LED，并等待其发 光稳定；
[0067] S103:在各类单颗LED垂直正下方放置X-Rite ColorChecker Passport标准灰板， 采用光谱辐射计Konica Minolta CS-2000依次测量获得各类单颗LED点亮稳定时经标准灰 板反射的光谱辐射亮度分布数据，记为S^(A)，表示第i类单颗LED在驱动值点亮稳定后 经标准灰板反射的光谱辐射亮度分布，单位为W/(sr · m2 · nm)，λ为可见光波长，测量波长 范围为380~780nm，测量波长间隔为lnm。
[0068] 2.光源目标色品和照度数据的输入及处理具体包括以下步骤：
[0069] S201:输入光源的目标色品数据，其中目标色品数据可以为光源的目标相关色温 T、光源的目标相关色温T和目标Duv或者光源的任意一种目标色品坐标，然后根据色度学原 理将目标色品输入数据转换为目标CIE1931色品坐标(xt，yt)和目标CIE1976UCS色品坐标 (u/，ν/ );
[0070] 选取 28001(、35001(、40001(、45001(、50001(、55001(、60001(、65001(等共计8种相关色温 作为光源目标色品输入数据，选取的相关色温覆盖了常用的照明色温范围，由此可计算这8 种相关色温对应的目标CIE1931色品坐标(x，y)和CIE1976UCS色品坐标(Υ ,ν')，计算结果 如表1所示，需要说明的是，本发明并不局限于这些色温点，只要目标色品点在多色LED光源 的混光范围内，均适用于本发明；
[0071] 表1目标CIE1931色品坐标（x，y)和CIE1976UCS色品坐标(ι/ ,ν')
[0072]
[0073] S202:与步骤S201中8种相关色温对应的待匹配照度均设置为Et = 5001x，并结合 标准灰板的反射比，根据光度学和色度学原理计算待匹配光源照度Et和步骤S201中目标色 品坐标（xt，yt)对应的绝对三刺激值（Xt，Yt，Z t)，标准灰板的反射比通过分光光度计 GretagMacBeth Color-Eye 7000A测量计算获得，需要说明的是，本发明并不局限于这一照 度值，只要待匹配照度在多色LED所能实现的照度范围之内均适用本发明。
[0074] 3.光谱优化采用一种约束型差分进化算法进行迭代优化，其具体步骤为：
[0075] S301:根据光照度叠加原理计算LED矩阵中各类LED的等效数量，可用向量表示为 Nmax = [Nmax,i，Nmax,2，…，Nmax, n]，其中η为LED矩阵中LED的种类数，SPη = 17;
[0076] S302:按照下式计算各类LED在采样驱动值均为du打开时的混合光谱辐射亮度分 布Si，j，M(X):
[0077] Si, =Nmax，iSi，j(A) ( 1 )
[0078] 式中Nmax,i为步骤S301中第i类LED的等效数量，Si,j(A)为步骤S103中单颗LED在采 样驱动值为di,j时测量的光谱福射亮度分布；
[0079] S303:通过三次样条插值算法预测各类LED在其驱动值动态范围内任一驱动值时 的光谱福射亮度分布Si,M(X)，gp
[0080] Si，M(X) =spline(di，Si(A)，Di) (2)
[0081 ]式中spline表示执行三次样条插值，di为步骤S101中第i类LED采样驱动值构成的 行向量为步骤S302中第i类LED在波长λ处不同采样驱动值的光谱辐射亮度构成的行 向量，Di为第i类LED在其驱动值动态范围内的任一驱动值；
[0082] S304:定义目标优化函数，将给定光源目标色品和照度参数的光谱优化转换为在 满足式(4)的约束条件下求解式(3)的最大值问题：
[0083]
[0084]
[0085]
[0086]
[0087]
[0088]式中R( ·)采用CIE-般显色指数Ra作为目标优化函数，且R( ·)为光源光谱的函 数，需要说明的是，本发明并不局限于CIE-般显色指数，只要为光源显色性或光源显色性 与光效的组合型评价指标，均在本发明的保护范围之内，Sm(A)表示多色LED的混合光谱， Si,Μ(λ)为第i类LED在其驱动值动态范围内任一驱动值处的光谱辐射亮度分布，可通过步骤 S303获得，D为待求解行向量，向量中各元素 Di表示第i类LED在其驱动值动态范围内任一驱 动值，〒U).只1). 7(/i)为CIE1931标准色度观察者，n = 17;
[0089] S305:执行初始化操作，即在0~Dmax间随机生成Np个均匀分布的候选解向量集，记 为〇15,1=[01,1<,1，02九1，"_，011九1]，其中11=17，1^表示第1^个候选解向量，且1^={1，2，...，邮}， 取Np= 180，1代表初始候选解向量集，即第1次迭代时的候选解向量集，同样将第G次迭代时 的候选解向量集记为〇1^=[01,以，02,1<,0，"，011,1^];
[0090] S306:将步骤S305中180个候选解向量Dk,G依次代入步骤S304中的式(3)和式(4);
[0091] S307:从180个候选解向量Dk,c中随机选择s个候选解向量，记为Dp.cztDnc， D2,p,g，…，Dn,P,G]，其中，p= {1，2, · · ·，s}，取s = 9;
[0092] S308:根据经典差分进化算法原理对步骤S307中的9个候选解向量DP,G依次执行变 异和交叉操作，经此操作共获得9个候选解向量UDP, c;
[0093] S309:将步骤S308中9个候选解向量UDP,G依次代入步骤S304中的式(3)和式(4); [0094] S310:根据步骤S304中的式(3)和式(4)比较步骤S308中9个候选解向量UDP,G以及 步骤S307中9个候选解向量D P,c的计算结果，并用前者最优候选解向量取代后者中与此相比 较差的候选解向量DP, c，同样也对候选解向量集Dk, c中的候选解向量进行取代；
[0095] S311:若步骤S308中9个候选解向量UDP,G的所有最优候选解向量均不满足式(4)， 则将这些最优候选解向量存储，并每隔20次迭代用其随机取代候选解向量集Dk,c中相同数 量的候选解向量，同时确保候选解向量集Dk, c中的最优候选解向量不被取代；
[0096] S312:重复执行步骤S307~S311，直至满足设定的迭代终止条件，即最大迭代次数 Gmax，此处Gmax=1200,可根据实际优化情况对其适当调节，迭代次数设置越大匹配效果越 好，但会降低匹配效率，并将此时获得的最佳解向量记为Dt,end，且t={l，2,. . .，N}，N为优化 执行次数；
[0097] S313:重复执行步骤S305~S312,直至满足设定的优化执行次数N，此实施例N = 25，多次执行优化的目的是确保获得最佳解向量，也可根据实际优化情况对其适当调节，然 后，根据步骤S304中的式(3)和式(4)，从25次优化获得的25个最佳解向量Dt,end中选出最终 的最佳解向量，并将其记为D。, end，即为优化获得的各类LED的驱动值；
[0098] S314:将步骤S313中优化获得的各类LED的驱动值输入到编写的控制软件中点亮 LED矩阵，通过CS-2000测量其光谱辐射亮度分布Stm(A)、相关色温!"以及CIE1976UCS色品坐 标(u/，ν/ )，由此可获得与步骤S201中8种目标相关色温T对应的最佳光谱分布，分别如图1 和图2所示，相关色温和CIE1976UCS色品坐标测量结果见表2,同时通过ΧΥΙ-ΙΙΙ全数字照度 计测量相应的照度E m，测量结果如表2所示。
[0099] 表2相关色温、CIE1976UCS色品坐标以及照度测量结果
[0100]
[0101 ] 4.光谱优化评价具体包括以下步骤：
[0102] S401:采用色温误差百分比ετ评价优化测量出的相关色温与目标相关色温间的差 异，ετ的定义为
[0103]
(6)
[0104] 式中Τ为步骤S201中的8种目标相关色温，Tm为步骤S314中的8种相关色温测量值， 相关色温评价结果见表3;
[0105] S402:采用CIE1976UCS色品差Δ w评价优化测量出的光源色品与目标色品间的 差异，Aw的定义为
[0106]
(7)
[0107] 式中（ιι^,ν^)为步骤S201中8种目标相关色温对应的目标CIE1976UCS色品坐标， (ιι、，ν\)为步骤S314中优化实际测量的CIE1976UCS色品坐标，色品差评价结果见表3;
[0108] S403:采用照度误差百分比εΕ对照度匹配结果进行评价，εΕ的定义为
[0109]
(8)；
[0110] 式中Et为步骤S202中的待匹配照度，即5001x，Em为步骤S314中的照度测量值，照度 误差评价结果见表3;
[0111] 由表3可见，本发明提供的光谱优化方法能够实现光源目标色温、色品以及照度的 精确匹配；
[0112] 表3相关色温、色品以及照度评价结果
[0113]
[0114] S404:分别采用CIE-般显色指数特殊显色指数R9以及IES颜色保真性指数Rf 和色域面积指数心对步骤S314中优化测量出的光谱Stm(A)进行显色性评价，显色性具体评 价结果见表4。
[0115] 表4显色性评价结果
[0116]
[0117]由表4可见，本发明提供的光谱优化方法在广泛的色温范围内均获得了优异的光 源显色性。该方法以获取光源最佳显色性为目标，可为获得光源色温和照度的连续可调提 供一种切实有效的实现方法。
【主权项】
1. 一种多色L邸实现光源最佳显色性的光谱优化方法，其特征在于，包括如下步骤： S1: LED光谱样本数据的获取； S2:光源目标色品和照度数据的输入及处理，得到其对应的目标绝对Ξ刺激值； S3:采用约束型差分进化算法进行迭代优化，获得匹配光源目标绝对Ξ刺激值且实现 其最佳显色性的各类L邸的驱动值； S4:对上述驱动值下测量得到的光谱进行优化评价。2. 如权利要求1所述的多色LED实现光源最佳显色性的光谱优化方法，其特征在于，所 述步骤S1中L邸光谱样本数据的获取过程具体包括W下步骤： S101:在各类Lm)的驱动值动态范围内W-定采样间隔建立Lm)驱动数据样本，记为 di,j，表示第i类L邸的第j个驱动数据样本，用行向量将每类L邸的驱动数据样本表示为di = [di,l，di,2, . . .，di,max]，其中di,max表示第i类L抓的最大采样驱动值，并将各类LED最大采样 驱动值构成的行向量记为〇。3、=[山,。3、，(12,。3、，...，山,。3^，其中11为1邸的种类数； S102:将步骤S101中的驱动数据样本分别依次驱动点亮各类单颗LED,并等待其发光稳 定； S103:在各类单颗Lm)垂直正下方放置标准板，采用光谱福射计依次测量获得各类单颗 L邸点亮稳定时经标准板反射的光谱福射亮度分布数据，记为Si, ^λ)，表示第i类单颗Lm)在 驱动值di,j点亮稳定后经标准板反射的光谱福射亮度分布，单位为W/(sr · m2 · nm)，λ为可 见光波长。3. 如权利要求1所述的多色LED实现光源最佳显色性的光谱优化方法，其特征在于，所 述步骤S2中光源目标色品和照度数据的输入及处理具体包括W下步骤： S201:输入光源的目标色品数据，其中目标色品数据为光源的目标相关色溫T、光源的 目标相关色溫T和目标Duv或者光源的任意一种目标色品坐标，然后根据色度学原理将目标 色品输入数据转换为目标CIE1931色品坐标(xt，yt)和目标CIE197抓CS色品坐标(u/，v/); S202:输入待匹配光源照度Et，单位为lx，并结合标准板的反射比，根据光度学和色度学 原理计算待匹配光源照度Et和步骤S201中目标色品坐标(xt，yt)对应的目标绝对Ξ刺激值 (Xt，Yt，Zt)，标准板的反射比通过分光光度计测量计算获得。4. 如权利要求1所述的多色LED实现光源最佳显色性的光谱优化方法，其特征在于，所 述步骤S3中光谱优化采用约束型差分进化算法进行迭代优化，其具体步骤为： S301:根据光照度叠加原理计算多色LED光源中各类Lm)的等效数量，并用行向量表示 为Nmax = [ Nmax, 1，Nmax, 2，…，Nmax, η]，其中 η为LED的种类数； S302:按照下式计算各类LED在采样驱动值均为di,苗Τ开时的混合光谱福射亮度分布 Si,川(入）： -Nmax, iSi, j (λ.) ( 1 ) 式中Nmax,i为步骤S301中第i类Lm)的等效数量，Sij(A)为步骤S103中单颗Lm)在采样驱 动值为di,j时测量的光谱福射亮度分布； S303:通过Ξ次样条插值算法预测各类Lm)在其驱动值动态范围内任一驱动值时的光 谱福射亮度分布Si,Μ(λ)，即 8?,Μ(λ) = 3ρ1?ηθ((1?,8?(λ) ,Di) (2) 式中spline表示执行Ξ次样条插值，di为步骤SlOl中第i类Lm)采样驱动值构成的行向 量，5ι(λ)为步骤S302中第i类LED在波长λ处不同采样驱动值的光谱福射亮度构成的行向 量，Di为第i类Lm)在其驱动值动态范围内的任一驱动值； S304:定义目标优化函数，将给定光源目标色品和照度参数的光谱优化转换为在满足 式(4)的约束条件下求解式(3)的最大值问题：式中R( ·)表示某种光源显色性评价指标或者光源显色性与光效的组合型评价指标， 且R( ·)为光源光谱的函数，Sm(A)表示多色L抓的混合光谱，Si,M(A)为第i类L抓在其驱动值 动态范围内任一驱动值处的光谱福射亮度分布，可通过步骤S303获得，D为待求解行向量， 向量中各元素 Di表示第i类LED在其驱动值动态范围内任一驱动值，xU)，巧A).巧/1)为 CIE1931标准色度观察者，η为L邸的种类数； S305:执行初始化操作，即在0~Dmax间随机生成化个均匀分布的候选解向量集，记为 0^=[〇1山1，化,、1，。'，0。山1]，其中1^表示第1^个候选解向量，且1^={1，2，...，化}，1代表初 始候选解向量集，即第1次迭代时的候选解向量集，同样将第G次迭代时的候选解向量集记 为Dk,G= [Dl,k,G,〇2,k,G, ... ,Dn,k,G]; S306:将步骤S305中化个候选解向量Dk,G依次代入步骤S304中的式(3)和式(4); S307:从化个候选解向量Dk,G中随机选择S个候选解向量，记为Dp,G= [Di,p,g,〇2,p,g,…， Dn,p,G],其中 P 二{1,2,. . . ,s}; S308:根据差分进化算法原理对步骤S307中的S个候选解向量Dp,G依次执行变异和交叉 操作，经此操作共获得S个候选解向量UDp, C; S309:将步骤S308中S个候选解向量UDp,g依次代入步骤S304中的式(3)和式(4); S310:根据步骤S304中的式（3)和式(4)比较步骤S308中S个候选解向量UDp,gW及步骤 S307中S个候选解向量Dp,G的计算结果，并用前者最优候选解向量取代后者中与此相比较差 的候选解向量Dp,G，同样也对候选解向量集Dk,G中的候选解向量进行取代； S311:若步骤S308中S个候选解向量UDp,g的所有最优候选解向量均不满足式(4)，则将 运些最优候选解向量存储，并每隔m次迭代用其随机取代候选解向量集Dk,c中相同数量的候 选解向量，同时确保候选解向量集Dk,冲的最优候选解向量不被取代； S312:重复执行步骤S307~S311，直至满足设定的迭代终止条件，即最大迭代次数Gmax， 并将此时获得的最佳解向量记为Dt,end，且t= {1，2,. . .，N}，N为优化执行次数； S313:重复执行步骤S305~S312,直至满足设定的优化执行次数N，并根据步骤S304中 的式(3)和式(4)，从N次优化获得的N个最佳解向量Dt,en冲选出最终的最佳解向量，并将其 记为D。,end,即为优化获得的各类L邸的驱动值； S314:将步骤S313中优化获得的各类LED的驱动值输入到相应的控制软件点亮LED，并 通过光谱福射计测量光谱福射亮度分布Stm(A)、相关色溫TmW及CIE197抓CS色品坐标(Um^， Vn/ )，通过照度计测量相应的照度Em。5.如权利要求1所述的多色LED实现光源最佳显色性的光谱优化方法，其特征在于，所 述步骤S4中光谱优化评价具体包括W下步骤： S401:采用色溫误差百分比ετ评价优化测量出的相关色溫与目标相关色溫间的差异，ετ 的定义为(6) 式中Τ为步骤S201中的目标相关色溫，Tm为步骤S314中的相关色溫测量值； S402:采用CIE197抓CS色品差Δ U' y评价优化测量出的光源色品与目标色品间的差异， Au'y的定义为(巧 式中（u/，v/ )为步骤S201中的目标CIE197抓CS色品坐标，（加/，ν"/ )为步骤S314中实际 测量的CIE197抓CS色品坐标； S403:采用照度误差百分比εΕ对照度匹配结果进行评价，εΕ的定义为做 式中Et为步骤S202中的待匹配照度，Em为步骤S314中的照度测量值； S404:分别采用CIE-般显色指数Ra和特殊显色指数R9W及IES颜色保真性指数Rf和色 域面积指数Rg对步骤S314中优化测量出的光谱StmW进行显色性评价。
【文档编号】G06F17/16GK105973572SQ201610270620
【公开日】2016年9月28日
【申请日】2016年4月27日
【发明人】徐海松, 章夫正, 汪哲弘
【申请人】浙江大学