# Plotting 4D with 3 vectors and 1 matrix

1 visualizzazione (ultimi 30 giorni)
Simao Nobrega il 20 Giu 2015
Commentato: Simao Nobrega il 21 Giu 2015
Hello,
I am solving a problem where I have three vectors (one for each 3D dimension), with different lengths and unevenly spaced, and a temperature matrix where each entry corresponds to the temperature at the position (x(i),y(j),z(k)). My goal is to make a 4D plot where the forth dimension (temperature) is represented my a colour (see attached file).
I have searched for an answer but did not get any solution for the conditions that I presented (3 vectors with 1 matrix). What function or approach would you recommend me to try?
*Data type*
length of x = M
length of y = N
length of z = P
size of Temperature matrix = M x N x P
##### 0 CommentiMostra -2 commenti meno recentiNascondi -2 commenti meno recenti

Accedi per commentare.

### Risposta accettata

Walter Roberson il 20 Giu 2015
You can create a scattered interpolant and then sample it along grids. But as is hinted in your image, the outside is going to hide the inside when you plot. You can work with a voxel viewing routine or you can create isosurfaces to try to deal with that.
##### 4 CommentiMostra 2 commenti meno recentiNascondi 2 commenti meno recenti
Walter Roberson il 21 Giu 2015
[X, Y, Z] = ndgrid(x, y, z);
slice(X, Y, Z, temperature, [min(x), max(x)], [min(y), max(y)], [min(z), max(z)], 'cubic')
Simao Nobrega il 21 Giu 2015
Worked like charm! Slice is a very nice function.
Best regards

Accedi per commentare.

### Più risposte (1)

Hugo il 20 Giu 2015
You could use "patch" since the surface will hide the interior points. You can create a matrix containing the coordinates of each point on the surface and then another matrix containing the indexes of the vertices defining each face in the patch. Suppose you want to create the faces corresponding to the top. You can use
xaux = repmat(x(:),1,Ny);
yaux = repmat(y(:)',Nx,1);
zaux = repmat(zval,Nx,Ny);
vertices = [xaux(:),yaux(:),zaux(:)];
faces = [];
for indy = 1:Ny,
for indx=1:Nx,
faces=[faces;[indx,indx+1,indx+11,indx+10]+(indy-1)*10];
end,
end
patch('Vertices',vertices,'Faces',faces,'CData',reshape(temperature(:,:,end),Nx*Ny,1),'FaceColor','interp','EdgeColor','none');
You will need to repeat this for each side. You can use any grid you want. The interpolation may be rather unsatisfactory. Hope this helps.
##### 1 CommentoMostra -1 commenti meno recentiNascondi -1 commenti meno recenti
Simao Nobrega il 21 Giu 2015
Hello Hugo,
Thank you for your answer. In the attached image you can see the temperature distribution of the top plane using contourf and your code. Maybe the problem is the interpolation that is being made in your code.

Accedi per commentare.

### Categorie

Scopri di più su Surface and Mesh Plots in Help Center e File Exchange

### Community Treasure Hunt

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

Start Hunting!

Translated by