分享一个猴鞍面的代码

panphp · 发表于 2020-5-28 10:31

EDA365欢迎您登录！

您需要登录才可以下载或查看，没有帐号？注册

x

分享一个猴鞍面的代码：

% monkey_saddle.m
%
% Make an animation of a monkey saddle suRFace.

% Name of animation created
filename = 'monkey_saddle.gif';

% Polar coordinates
r = linspace(0, 1);
t = linspace(-pi, pi);
[R, T] = meshgrid(r, t);

% Cartesian coordinates
x = R .* cos(T);
y = R .* sin(T);
z = x.^3 - 3 * x .* y.^2;

% Gaussian curvature to color the surface
K = - 36 * (x.^2 + y.^2)...
  ./ (1 + 9 * (x.^2 + y.^2).^2).^2;

% Plot the surface
chartSurface = surf(x, y, z, K);

% Edit its appearance
chartSurface.EdgeColor = 'none';
chartSurface.FaceColor = 'interp';
chartSurface.FaceLighting = 'gouraud';

ax = gca;

% Remove tick labels
ax.XTickLabel = {};
ax.YTickLabel = {};
ax.ZTickLabel = {};

% Rotate the camera's azimuthal angle
axis vis3d;
v = ax.View;
az = v(1);
el = v(2);

% Minimize file size by taking advantage of 3-fold symmetry of the surface
for k = 0:119
  ax.View = [k, el];
  drawnow;
  frame = getframe(1);
  im = frame2im(frame);
  [A,map] = rgb2ind(im,256);
  if k == 0;
imwrite(A,map,filename,'gif','LoopCount',Inf,'DelayTime',0);
  else
imwrite(A,map,filename,'gif','WriteMode','append','DelayTime',0);
  end
end

duoer · 发表于 2020-5-28 11:27

看看楼主的代码。

帐号		自动登录	找回密码
密码			注册

分享一个猴鞍面的代码

EDA365欢迎您登录！

浏览过的版块

推荐内容 /1