yield = 50000; %% 坐标尺寸 A_Dyucker_in_Three_dimensiona_in_pai_space_curve_3D1 limxyz = 60; piii = 1; f1=@(x,y,z) x+y+z; f2=@(x,y,z) max(abs(x),abs(y),abs(z))-((min(abs(x),abs(y))+min(abs(y),abs(z))+min(abs(z),abs(x)))/2)-yield; %% 1./((1./(3.*sqrt(3)-b3./27))).* [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 = 'k'; S.LineWidth = 3; S.LineStyle = '-'; hold on 这段代码运行提示:错误使用 max 不支持具有两个要比较的矩阵和一个工作维度的 MAX。 出错 New_Cazacu_Barlat_Of_Third_Effect_in_pai_space_curve_3D_4>@(x,y,z)max(abs(x),abs(y),abs(z))-((min(abs(x),abs(y))+min(abs(y),abs(z))+min(abs(z),abs(x)))/2)-yield 出错 New_Cazacu_Barlat_Of_Third_Effect_in_pai_space_curve_3D_4 (line 142) CV=f1(CX,CY,CZ)-f2(CX,CY,CZ);
时间: 2023-11-19 09:06:03 浏览: 77
根据您提供的代码,您的问题可能是由于 `max` 函数不支持比较两个矩阵和一个工作维度引起的。可能是由于您在 `f2` 中的 `max` 函数的输入参数不正确,需要修改为对每个坐标点进行比较。您可以将 `f2` 函数修改为如下形式:
```matlab
f2 = @(x, y, z) max([abs(x(:)), abs(y(:)), abs(z(:))], [], 2) - ((min([abs(x(:)), abs(y(:)), abs(z(:))], [], 2) + min([abs(y(:)), abs(z(:)), abs(x(:))], [], 2) + min([abs(z(:)), abs(x(:)), abs(y(:))], [], 2)) / 2 - yield;
```
上述代码中,`max` 函数和 `min` 函数的第二个输入参数为 `[]`,表示沿着第二个维度比较,即对每个坐标点进行比较。此外,我们使用 `(:)` 将矩阵变成列向量,以便和常数 `yield` 进行运算。
请注意,修改后的代码仅解决了 `max` 函数可能引起的问题。如果您仍然无法运行代码,请提供完整的错误信息和相关的代码,以便更好地帮助您解决问题。
相关问题
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` 函数不能同时比较两个矩阵和一个工作维度,因此需要将其修改为逐元素比较。
下面是修改后的代码,你可以参考一下:
```matlab
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;
```
这段代码应该可以正常运行,并绘制出对应的曲面和等势面。
阅读全文