用excel给定的源光谱数据和CIE1931-xyz-三刺激值数据，用matlab求出给定光谱的XYZ 和光谱的xyz色坐标；

步骤：

1. 打开Excel文件，将源光谱数据和CIE1931-xyz-三刺激值数据复制到两个新建的工作簿中（注意：将源光谱数据粘贴到第一个工作簿中，将CIE1931-xyz-三刺激值数据粘贴到第二个工作簿中）。

2. 在matlab中，使用xlsread函数读取两个工作簿中的数据，并将其存储在两个矩阵中。

  % 读取源光谱数据
  [lambda, intensity] = xlsread('source_spectrum.xlsx');
  % 读取CIE1931-xyz-三刺激值数据
  [wavelength, x, y, z] = xlsread('CIE1931-xyz-tristimulus-values.xlsx');

3. 将源光谱数据和CIE1931-xyz-三刺激值数据进行插值操作，使它们在波长范围内等间距。

  % 对源光谱数据进行插值
  new_lambda = linspace(380, 780, 201);
  new_intensity = interp1(lambda, intensity, new_lambda, 'spline');
  % 对CIE1931-xyz-三刺激值数据进行插值
  new_x = interp1(wavelength, x, new_lambda, 'spline');
  new_y = interp1(wavelength, y, new_lambda, 'spline');
  new_z = interp1(wavelength, z, new_lambda, 'spline');

4. 计算XYZ色值。

  % 计算XYZ色值
  M = [new_x', new_y', new_z'];
  XYZ = trapz(new_lambda, new_intensity .* M);

5. 计算光谱的xyz色坐标。

  % 计算光谱的xyz色坐标
  x_bar = trapz(new_lambda, new_x .* new_intensity) / trapz(new_lambda, new_intensity);
  y_bar = trapz(new_lambda, new_y .* new_intensity) / trapz(new_lambda, new_intensity);
  z_bar = trapz(new_lambda, new_z .* new_intensity) / trapz(new_lambda, new_intensity);

6. 输出XYZ色值和光谱的xyz色坐标。

  % 输出XYZ色值
  fprintf('XYZ = [%f, %f, %f]\n', XYZ(1), XYZ(2), XYZ(3));
  % 输出光谱的xyz色坐标
  fprintf('x = %f, y = %f, z = %f\n', x_bar, y_bar, z_bar);

完整代码如下：

% 读取源光谱数据
[lambda, intensity] = xlsread('source_spectrum.xlsx');
% 读取CIE1931-xyz-三刺激值数据
[wavelength, x, y, z] = xlsread('CIE1931-xyz-tristimulus-values.xlsx');

% 对源光谱数据进行插值
new_lambda = linspace(380, 780, 201);
new_intensity = interp1(lambda, intensity, new_lambda, 'spline');
% 对CIE1931-xyz-三刺激值数据进行插值
new_x = interp1(wavelength, x, new_lambda, 'spline');
new_y = interp1(wavelength, y, new_lambda, 'spline');
new_z = interp1(wavelength, z, new_lambda, 'spline');

% 计算XYZ色值
M = [new_x', new_y', new_z'];
XYZ = trapz(new_lambda, new_intensity .* M);

% 计算光谱的xyz色坐标
x_bar = trapz(new_lambda, new_x .* new_intensity) / trapz(new_lambda, new_intensity);
y_bar = trapz(new_lambda, new_y .* new_intensity) / trapz(new_lambda, new_intensity);
z_bar = trapz(new_lambda, new_z .* new_intensity) / trapz(new_lambda, new_intensity);

% 输出XYZ色值和光谱的xyz色坐标
fprintf('XYZ = [%f, %f, %f]\n', XYZ(1), XYZ(2), XYZ(3));
fprintf('x = %f, y = %f, z = %f\n', x_bar, y_bar, z_bar);