2D Matrix revolution

Hi community

How can I revolutionize a 2D matrix around an axis to obtain a solid?

Supose a NxN matrix. In each row are the points of a kind of gaussian curve with less amplitude and higher width. Plotting this surface shows something like a half bell that decays to zero. (See Bottom)

I want to obtain a 3D image revolutionazing that matrix, so the solid looks like a complete bell, or a mexican hat.

Notice that each curve in the matrix rows are different.

I have already try with cylinder(), but i can't define the rotation axis, and the result is a segmented solid.

Possible solutions so I can get the 3D image and define the rotation axis?

Thanks a lot

The surface to revolutioniye looks like this:

Cylinder shows this:

[darova]

Which axis do you want to rotate?

[Felipe Bayona]

I would like to know the general method, where i can specify the axis. In this particular example, I need to do it around the z-axis, to obtain the 'Mexican hat'

Answer 1

[darova]

Use cylinder for each curve separately:

z = linspace(-1,1,20);                  % axis
r = sqrt(1-z.^2);                       % radius
[x,y,~] = cylinder(r,length(z)-1);      % convert to cylinder system
zz = repmat(z',[1 20]);                 % create 2D matrix for Z matrix
surf(x,y,zz)                            % sphere

[Felipe Bayona]

And how should it be?

in this case r is a 1x20 array, and its plotted with the x and y resuls of cylinder.

How it must be for each curve separately? (for example 10 or 100 r's of 1x20)

[darova]

attach the data

[darova]

And what is the problem?

x = linspace(0,40,41);
[X,T, rc] = normalDiffusion(4, 1, 1, x, 20);
z = X(i,:);                             % axis
r = rc(i,:);                            % radius
[x,y,~] = cylinder(r,length(z)-1);      % convert to cylinder system
zz = repmat(z',[1 length(r)]);          % create 2D matrix for Z matrix
surf(x,y,zz)                            % sphere

[darova]

Maybe z and r should be switched?

z = X(i,:);                             % axis
r = rc(i,:);                            % radius

Answer 2

[Philippe Lebel]

Have you tried this function?

https://www.mathworks.com/matlabcentral/fileexchange/15559-solids-of-revolution

[Felipe Bayona]

This only alows to rotate a y=f(x) funtion

My matrix is composed by many curves, because my function is z=f(x,y)