RGB_Func.m

function RGB_11 = RGB_Func(Data_Per1_C1, r, g, b, I)
Dimen = 60;
SetData = Data_Per1_C1';% 将数据变为 33 * 65536的
% 初始化矩阵
X = [];Y = [];Z = [];
R1 = [];G1 = [];B1 = [];
R = [];G = [];B = [];
RGB_11 = [];

Row_ = 1;Line_ = 26;% Row_:Line_ 选取矩阵的维度；

% 确定公式
gI = g(Row_:Line_,:) .* I(Row_:Line_,:);
k = 100 / sum(gI);% 比例系数
X = k * sum(bsxfun(@times, r(Row_:Line_,:) .* I(Row_:Line_,:) , SetData(Row_:Line_,:)))';
Y = k * sum(bsxfun(@times, g(Row_:Line_,:) .* I(Row_:Line_,:) , SetData(Row_:Line_,:)))';
Z = k * sum(bsxfun(@times, b(Row_:Line_,:) .* I(Row_:Line_,:) , SetData(Row_:Line_,:)))';

XYZ = [X,Y,Z];
[rgb] = xyz2rgb(XYZ);% 计算RGB
                
R1 = rgb(:, 1);G1=rgb(:, 2);B1 = rgb(:, 3);
R = reshape(R1,Dimen,Dimen);G = reshape(G1,Dimen,Dimen);B = reshape(B1,Dimen,Dimen);  % 重构 R，G，B 
                
RGB_11 = uint8(cat(3,R,G,B));% 将R，G，B合成为三维矩阵,并转为矩阵；
end