matlab程序绘制CIE1931色品图的程序

可以参考以下代码实现：

% 设置色度坐标
x = linspace(0, 1, 101);
y = linspace(0, 1, 101);
[X, Y] = meshgrid(x, y);

% 计算CIE 1931色匹配函数
T_xyz = load('T_xyz1931.mat');
T_xyz = T_xyz.T_xyz1931;
xbar = interp1(T_xyz(:,1), T_xyz(:,2), X(:), 'linear', 0);
ybar = interp1(T_xyz(:,1), T_xyz(:,3), X(:), 'linear', 0);
zbar = interp1(T_xyz(:,1), T_xyz(:,4), X(:), 'linear', 0);

% 计算归一化三刺激值和彩色度坐标
k = 100./(ybar'*ones(1,length(x)));
L = k.*ybar.*x;
M = k.*ybar.*y;
S = k.*(1-ybar-x);

% 绘制CIE 1931色品图
figure;
scatter(M(:), L(:), 10, [S(:), M(:), L(:)], 'filled');
axis equal; axis([0 0.9 0 0.7]); xlabel('Y'); ylabel('X');
title('CIE 1931色品图');

该程序需要依赖T_xyz1931.mat文件，可以在网上下载。运行程序后，将绘制CIE1931色品图，其中背景颜色表示S通道（蓝色刺激）的值，x轴表示L通道（亮度），y轴表示M通道（红绿刺激）。