matlab程序将cie1931的色坐标绘制在cie1931马蹄图中

以下是一个简单的 MATLAB 程序，用于将 CIE 1931 色坐标绘制在 CIE 1931 马蹄图中：

% 定义 CIE 1931 色坐标的范围
x_range = 0:0.01:0.8;
y_range = 0:0.01:0.9;

% 生成网格数据
[X,Y] = meshgrid(x_range, y_range);

% 计算马蹄线方程
k = (Y - 0.332)/(0.1858 - X);
m1 = k < 0;
m2 = k >= 0 & k <= 1.611;
m3 = k > 1.611;
x1 = (-0.56*k.^3 + 0.17*k.^2 + 0.231*k).*m1;
x2 = (-1.376*k.^3 + 1.0745*k.^2 + 0.2*k + 0.1888).*m2;
x3 = (1.376*(k-1.611).^3 + 0.7035*(k-1.611).^2 - 0.007).*m3;
x = (x1 + x2 + x3)./(k+0.0228);

% 绘制马蹄图
figure;
plot(x1, Y.*m1, 'b',...  % 低亮度区域
     x2, Y.*m2, 'b',...  % 中亮度区域
     x3, Y.*m3, 'b');    % 高亮度区域
hold on;
scatter(X(:), Y(:), 10, 'filled');  % 绘制色坐标点
xlabel('x');
ylabel('y');
title('CIE 1931 Chromaticity Diagram');
axis equal;

运行这个程序，将会生成一个 CIE 1931 马蹄图，其中色坐标点将被绘制在图中。你可以根据自己的需求修改代码中的色坐标范围和点的大小等参数。