Jak produkować interpolowana gładką 3d działki w Matlab

Question

Wykres ten jest tworzony przez Mathematica:

ls = Table[Sinc[x*y], {x, -5, 5, 0.2}, {y, -5, 5, 0.2}]; 
ListPlot3D[ls, InterpolationOrder -> 2, PlotRange -> All, 
Mesh -> None] 

Jak stworzyć wykres takiego w Matlab?

Oto moja próba tak daleko:

>> x=linspace(-5.,5.,51); 
>> y=linspace(-5.,5.,51); 
>> [x,y]=meshgrid(x,y); 
>> z=sinc(x.*y); 
>> surf(x,y,z) 
>> shading interp 

To wygląda zupełnie inaczej, zwłaszcza szczegóły wsady. Czy można stworzyć fabułę podobną do Mathematica, zwłaszcza gładkość, cienie?

Answer 1

Aby stworzyć miłą oświetlenie i cienie, trzeba dodać światło na swojej działce, i dodać jakieś oświetlenia twarzy. Jeśli rozdzielczość jest zbyt niska, to skończy się nieco spartaczoną fabułą, ponieważ cieniowanie w stylu "interp" wykorzystuje interpolację liniową. Na przykład

n = 51; 
x=linspace(-5., 5., n); 
y=linspace(-5., 5., n); 
[x, y]=meshgrid(x, y); 
sinc = @(x) sin(x)./x; 
z=sinc(x.*y); 
z(isnan(z)) = 1; 
surf(x, y, z, 'LineStyle', 'none', 'FaceColor', 'interp') 
colormap(cool) 
camlight right 
set(gca, 'CameraPosition', [45 35 9.8]) 

która produkuje następujące

Zauważ, że jak gładka powierzchnia wydaje się związane z n. Większe wartości n zwiększą gładkość powierzchni.

Jeżeli dane które produkują jest drogie, aby utworzyć, można zwiększyć rozdzielczość przy użyciu bardziej zaawansowaną formę interpolacji niż liniowy, co następuje

n = 51; 
x=linspace(-5., 5., n); 
y=linspace(-5., 5., n); 
[x, y]=meshgrid(x, y); 
sinc = @(x) sin(x)./x; 
z=sinc(x.*y); 
z(isnan(z)) = 1; 

nn = 401; 
xi = linspace(-5.0, 5.0, nn); 
yi = xi; 
[xi, yi] = meshgrid(xi, yi); 
zi = interp2(x, y, z, xi, yi, 'spline'); 
surf(xi, yi, zi, 'LineStyle', 'none', 'FaceColor', 'interp') 
colormap(cool) 
camlight right 
set(gca, 'CameraPosition', [45 35 9.8]) 

która produkuje następujący obraz

Szczegółowe informacje można znaleźć na następujących stronach pomocy:

• Surface Properties

• camlight

• Coloring Mesh and Surface Plots

• interp2

Answer 2

Proponuję spróbować surfingu/surfl i spróbuj włączyć oświetlenie/wyłączania. Początkowa pozycja kamery jest również ważna, ponieważ używam kamery "reflektor".

x=linspace(-5.,5.,51); 
y=linspace(-5.,5.,51); 
[x,y]=meshgrid(x,y); 
z=sinc(x.*y); 
surfl(x,y,z) %surf(x, y, z) 
shading interp 
colormap cool 
%camlight headlight 
%lighting gouraud