example_streamslice.m

来自「基础程序」· M 代码 · 共 21 行

% example_streamslice.m
% 绘制流线切片图
% 首先创建曲面图
z=peaks;
surf(z);
shading interp
hold on
% 在曲面图上添加等值线图
[c ch]=contour3(z,20); set(ch,'edgecolor','b')
[u v]=gradient(z); 
% 绘制流线等值线切面图
h=streamslice(-u,-v); 
set(h,'color','k')
for i=1:length(h); 
    zi = interp2(z,get(h(i),'xdata'),get(h(i),'ydata'));
    set(h(i),'zdata',zi);
end
% 设置观看角度和方向
colorbar;
view(30,30); 
axis tight;