MATLAB-在一个图中绘制多个曲面拟合

vinayakshukl 发表于 Dev

Vinayakshukl

我有3套3D坐标，并且每套都适合平面。现在，我想在一个图中绘制所有数据点和3个平面。

到目前为止，我具有以下功能：

function [fitresult, gof] = create_fit(xx, yy, zz, grp)

[xData, yData, zData] = prepareSurfaceData( xx, yy, zz );

ft = fittype( 'poly11' );
opts = fitoptions( ft );
opts.Lower = [-Inf -Inf -Inf];
opts.Upper = [Inf Inf Inf];

hold on;
% figure( 'Name', 'fit1' );
[fitresult, gof] = fit( [xData, yData], zData, ft, opts );
h = plot( fitresult, [xData, yData], zData);
if grp == 1
    colormap(hot);
elseif grp == 2
    colormap(cool);
else
    colormap(grey);
end

legend( h, 'fit1', 'zz vs. xx, yy', 'Location', 'NorthEast' );

xlabel( 'xx' );
ylabel( 'yy' );
zlabel( 'zz' );
grid on

但是，这有两个问题：

当单独绘制平面时，将根据也随其绘制的数据点缩放这些平面。当将平面绘制在一起时，它们的缩放比例很差，并且数据点收敛到一个较小的斑点（与平面相比，比例很小），我尝试使用来解决此问题axis([-0.04 0.04 -0.04 0.04 -0.04 0.04 -1 1]);，但它是经过硬编码的，看上去仍然有些偏离。
该colormap命令似乎仅在第一次调用时才起作用。因此，所有飞机都变成蓝色。如何为每个平面和为该平面拟合的点着色？
这是绘制多个平面的最佳方法吗？

拉伯特

编辑

这是我的初始答案的编辑版本。的输出plot是一个包含两个元素的图形对象，因此您必须单独调用h(1)并h(2)设置平面和数据点的属性。

这是该函数的代码：

function [fitresult, gof, h] = create_fit(xx, yy, zz, color)

[xData, yData, zData] = prepareSurfaceData( xx, yy, zz );

ft = fittype( 'poly11' );
opts = fitoptions( ft );
opts.Lower = [-Inf -Inf -Inf];
opts.Upper = [Inf Inf Inf];

hold on;
[fitresult, gof] = fit( [xData, yData], zData, ft, opts );

h = plot( fitresult, [xData, yData], zData);

set(h(1), 'FaceColor', color);
set(h(2), 'MarkerFaceColor', color, 'MarkerEdgeColor', 'k');

这是一个调用函数的示例脚本：

% Define sample data
N = 20;
x = rand(1,N);
y = rand(1,N);
z = rand(1,N);

% Call the function, specify color
[f1, gof1, h1] = create_fit(x, y, z, 'r');
[f2, gof2, h2] = create_fit(x, y.^10, z, 'g');
[f3, gof3, h3] = create_fit(x.^10, y, z, 'b');

% Figure adjustments
xlabel( 'xx' );
ylabel( 'yy' );
zlabel( 'zz' );
view(3)
grid on

xlim([min(x) max(x)]);
ylim([min(y) max(y)]);
zlim([min(z) max(z)]);