显色指数的MATLAB计算程序

%一般显色指数  Ra 特殊显色指数 Ri
% 计算步骤 ： <色度学>  科学出版社  1979年  P256


function Ra = colorrender（Pinwavelength）

load colordata beta2 Uri Wri Vri ur cr dr vr M0 T0;

chromaticity（Pinwavelengthbeta2）;
load chromaticity ; 
M = colortemperature（ukvk） ;
%M = 103.2205 ;
MIN = min（abs（ M - M0 ）） ;
num = find（abs（ M - M0 ） == MIN） ;

cki = 1./vki .*（ 4 - uki - 10 .* vki ） ;
dki = 1./vki .*（1.708 .* vki + 0.404 - 1.481 .* uki ） ;

ck = 1/vk *（ 4 - uk - 10 * vk ） ;
dk = 1/vk *（1.708 * vk + 0.404 - 1.481 * uk ） ;

uki_ = （10.872 + 0.404 * cr（num）/ck .* ck