f1=@(x,y,z) x+y+z; f2=@(x,y,z) max(abs(x - y), abs(y - z), abs(z - x)) - yield; [X,Y]=meshgrid(-limxyz:piii:limxyz); Z1=-X-Y; % 绘制曲面 surf(X,Y,Z1,'EdgeColor','none','FaceAlpha',0,'FaceColor',[82,124,179]./255) hold on; fimplicit3(f2,'EdgeColor','none','FaceAlpha',0,'FaceColor',[169,64,71]./255) % 通过切片contour函数获取0等势面 [CX,CY,CZ]=meshgrid(-limxyz:piii:limxyz); CV=f1(CX,CY,CZ)-f2(CX,CY,CZ); S=contourslice(CX,CY,CZ,CV,X,Y,Z1,[0,0]); S.EdgeColor = 'c'; S.LineWidth = 3; S.LineStyle = '-'; hold on，运行提示： In fimplicit3 (line 120) In New_Cazacu_Barlat_Of_Third_Effect_in_pai_space_curve_3D_1 (line 122) 错误使用 max 不支持具有两个要比较的矩阵和一个工作维度的 MAX。；出错 New_Cazacu_Barlat_Of_Third_Effect_in_pai_space_curve_3D_1>@(x,y,z)max(abs(x-y),abs(y-z),abs(z-x))-yield 出错 New_Cazacu_Barlat_Of_Third_Effect_in_pai_space_curve_3D_1 (line 125) CV=f1(CX,CY,CZ)-f2(CX,CY,CZ);

根据错误提示，是因为 `f2` 函数中的 `max` 函数不能同时比较两个矩阵和一个工作维度，因此需要将其修改为逐元素比较。

下面是修改后的代码，你可以参考一下：

f1 = @(x,y,z) x+y+z;
f2 = @(x,y,z) max([abs(x-y),abs(y-z),abs(z-x)],[],2) - yield;

[X,Y] = meshgrid(-limxyz:piii:limxyz);
Z1 = -X - Y;

% 绘制曲面
surf(X,Y,Z1,'EdgeColor','none','FaceAlpha',0,'FaceColor',[82,124,179]./255)
hold on;

fimplicit3(f2,'EdgeColor','none','FaceAlpha',0,'FaceColor',[169,64,71]./255)

% 通过切片contour函数获取0等势面
[CX,CY,CZ] = meshgrid(-limxyz:piii:limxyz);
CV = f1(CX,CY,CZ) - f2(CX,CY,CZ);

S = contourslice(CX,CY,CZ,CV,X,Y,Z1,[0,0]);
S.EdgeColor = 'c';
S.LineWidth = 3;
S.LineStyle = '-';

hold off;

这段代码应该可以正常运行，并绘制出对应的曲面和等势面。