您可以像这样在 Matlab 中创建隐式函数的曲面图:
f = @(x,y,z) exp((625*(x - 31/20)^2)/72 + (625*(x - 981/250)^2)/72 + (625*(y - 461/100)^2)/72 + (625*(y - 4797/1000)^2)/72 + (625*(z - 13207/1000)^2)/72 + (625*(z - 15177/1000)^2)/72) - 0.0005;
fimplicit3(f)
但是,如果您分析 your equation in Wolfram Alpha,您会发现它的另一种形式是
Matlab 无法绘制此图,因为它将所有计算值解释为 Inf。您需要先对函数进行一些缩放。
您可以使用realmax 来检查 IEEE 双精度中最大的有限浮点数是多少。这会给你1.7977e+308。
编辑
Matlab 能够在没有8.12e1759 前置因子的情况下绘制上述函数的某些部分。这将如下所示:
但是,正如您所见,Matlab 无法在每个位置都显示一个表面,因为您的函数仍然会出现意外行为。
** 编辑 2 **
Ander Biguri 建议查看原点附近的函数。
如果我们这样做:
f = @(x,y,z) exp(5/288 .* (2 .* x .* (500 .* x - 2737) + y .* (1000 .* y - 9407) + 8 .* z .* (125 .* z - 3548))) - 0.0005;
interval = [-0.6 0.6 -0.6 0.6 -0.6 0.6];
%interval2 = [-6 6 -6 6 -6 6];
fimplicit3(f,interval)
colorbar
然后从顶部看结果,我们得到一个漂亮的等高线图,如下所示:
请注意,出于性能原因,我还对您的函数进行了矢量化。
我无法解释为什么对于较大的间隔,即使在这个较小的间隔中看起来很好的部分也无法由 Matlab 显示。例如,如果我改用 interval2,结果如下所示: