R语言绘制带y>x>0约束的z=x²+y²二元函数3D图问题咨询

核心思路是将不满足y>x>0约束的函数值设为NA，R的persp函数绘图时会自动跳过NA值对应的区域，自然实现约束范围的绘制效果。

修改后可直接运行的完整代码 #

# 定义目标函数
fxy <- function(x,y){
  return(x^2+y^2)
}
# 调整x的取值范围满足x>0的基础约束，减少无效计算
x <- seq(0, 100, by=5)
y <- seq(0, 100, by=5)
# 生成全量z矩阵
z <- outer(x, y, fxy)
# 将不满足y>x约束的位置赋值为NA，绘图时自动忽略
z[outer(x, y, function(a, b) b <= a)] <- NA

# 绘制3D图
persp(x, y, z,
  main="3D Plot of x²+y² (Constraint: y>x>0)",
  zlab = "Function Value",
  theta = 20, phi = 10,
  col = "blue",
  border = NA # 可选参数，去掉网格线让曲面更顺滑
)

关键调整说明 #

• 直接将x的取值范围限定为从0开始，满足x>0的基础约束，你也可以保留x负区间，配合后续条件过滤达到同样效果
• 使用outer(x,y, function(a,b) b<=a)生成和z维度完全匹配的逻辑矩阵，精准定位所有不满足y>x的坐标点，赋值为NA后不会被渲染
• 可选添加border = NA参数去掉曲面网格线，让显示效果更整洁

内容的提问来源于stack exchange，提问作者Math Universe