How to fit a surface to 3D dta points

25 visualizaciones (últimos 30 días)
Mos_bad
Mos_bad el 14 de Sept. de 2018
Comentada: Stephen23 el 14 de Sept. de 2018
Please find the attached. I want to fit a surface to show the trends of the 3d data points but I got an error that 'Z must be a matrix, not a scalar or vector'.

Respuestas (2)

Stephen23
Stephen23 el 14 de Sept. de 2018
Editada: Stephen23 el 14 de Sept. de 2018
Your data are scattered, not gridded:
surf only plots gridded data. To use surf you will either have to
  1. interpolate the values onto a grid, or
  2. fit a curve get gridded values and plot them.
Another option would be to use a Delaunay triangulation to plot the scattered data directly:
trisurf(delaunay(IM,Z50),IM,Z50,MnXdisp)
Gives:
This blog gives an nice explanation of options for scattered data:
  2 comentarios
Mos_bad
Mos_bad el 14 de Sept. de 2018
the data abe been like these first : Z50_ = linspace(0,20,20); IM_ = linspace(0.0025,1.2,20); [IM, Z50] = meshgrid(IM_, Z50_); Therefore, at each IM or Z50, there are 10 data from MnXdisp data set. I wonder if it is possible to fit 20 curves and then plot a surface.
Stephen23
Stephen23 el 14 de Sept. de 2018
@Mos_bad: well, the data you gave us might have been gridded at some point in history, but is now missing many data points. If you have data then of course you can fit curves to it (if that has any meaning depends on the data and what it represents).

Iniciar sesión para comentar.


KSSV
KSSV el 14 de Sept. de 2018
Editada: KSSV el 14 de Sept. de 2018
% Unstructred data plot
dt = delaunayTriangulation(IM,Z50) ;
t = dt.ConnectivityList ;
p = dt.Points ;
figure (1);
plot3(IM,Z50,MnXdisp,'b*');
hold on;
trisurf(t,p(:,1),p(:,2),MnXdisp')
title('unstructured')
% structured plot
x = IM ; y = Z50 ; z = MnXdisp ;
N = 50 ;
xi = linspace(min(x),max(x),N) ;
yi = linspace(min(y),max(y),N) ;
[X,Y] = meshgrid(xi,yi) ;
Z = griddata(x,y,z,X,Y) ;
figure (2);
plot3(IM,Z50,MnXdisp,'b*');
hold on;
surf(X,Y,Z)
title('structured')

Categorías

Más información sobre Interpolation en Help Center y File Exchange.

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by