Matlab绘制曳物面测地线泛函数返回全零的语法错误排查

问题根因 #

• 符号混淆错误：你定义的高度计算函数f = @(x,y) a.* (y - tanh(y))错误混淆了参数空间的v和笛卡尔坐标系下的y坐标。曳物面的z值直接由参数v决定，对应公式为z = a*(v - tanh(v))，和笛卡尔y坐标没有直接关系，传入笛卡尔y值计算自然得到错误的0结果。
• 变量顺序错误：你在定义dfdu、dfdv时还未对eps赋值，此时调用的是Matlab默认的极小值eps，和你后续定义的eps=0.0001不一致，也会带来精度问题。
• 参数更新逻辑问题：你将笛卡尔空间下的法向量分量直接用于参数空间U、V的更新，导致迭代得到的参数值不符合预期范围，进一步放大了计算误差。

快速修复方法 #

直接修改最终P矩阵的z值计算逻辑，用参数V直接计算高度即可解决z值为0的问题：

% 替换原代码中P的定义行
P = [f1(U,V); f2(U,V); a.*(V - tanh(V))];

同时把eps = 0.0001;这行代码移到定义dfdu、dfdv之前，保证求导逻辑使用的微小量一致。

可选优化方案 #

如果要保留完整的链式法则求导逻辑，可以直接基于参数化曲面计算偏导，不需要通过x、y做中间转换，既简化代码也避免符号混淆：

1. 直接用参数u、v计算曲面三个方向的偏导
2. 用参数偏导叉乘得到法向量，不需要额外定义高度函数f
3. 迭代时在参数空间做更新，不需要映射回笛卡尔空间做修正

内容的提问来源于stack exchange，提问作者Antoni Parellada