分享一个猴鞍面的代码

panphp

分享一个猴鞍面的代码：
( O* l) w! n4 K) N' ?) Q5 c0 e

) B6 z, H2 Y/ r9 [; i% monkey_saddle.m
%
2 ]4 t9 o& g2 I& o4 r7 ]- J% Make an animation of a monkey saddle suRFace.



% Name of animation created
5 M' F0 H( r- `# W* }( x* bfilename = 'monkey_saddle.gif';
0 D3 ~1 D) l" t. j) B
% Polar coordinates
r = linspace(0, 1);
t = linspace(-pi, pi);
: Y  d) t7 i% _8 t[R, T] = meshgrid(r, t);

; E; F* r% G- L  X% Cartesian coordinates
! Z3 U4 V- z% z/ [, I- v% Tx = R .* cos(T);
y = R .* sin(T);
9 X0 d, v, q. C- u% K8 \z = x.^3 - 3 * x .* y.^2;
  V' s: S" ^6 @7 M5 {8 Q* ^; r& R
, x9 m2 i1 B; Z8 R  M% Gaussian curvature to color the surface
K = - 36 * (x.^2 + y.^2)...
  ./ (1 + 9 * (x.^2 + y.^2).^2).^2;

  P/ [. u: c# t8 ?# Q5 Q% Plot the surface
0 O& ]6 A) X+ D$ d! Q' q0 R! l! ^chartSurface = surf(x, y, z, K);

% Edit its appearance
% u. I" ~# P' [+ D4 hchartSurface.EdgeColor    = 'none';
6 H' @/ l  @+ xchartSurface.FaceColor    = 'interp';
  R0 X" @4 ]# [! o/ D7 UchartSurface.FaceLighting = 'gouraud';
. \* {# v- X3 A5 O& ^: ~
ax = gca;

% Remove tick labels
ax.XTickLabel = {};
& Z' ~, E& F3 T- `+ i! ~ax.YTickLabel = {};
ax.ZTickLabel = {};
; t) F) o; b- C) S  }
% Rotate the camera's azimuthal angle
axis vis3d;
v = ax.View;
" m7 W4 G- s$ Paz = v(1);
el = v(2);

% Minimize file size by taking advantage of 3-fold symmetry of the surface
1 D' v  y8 @6 Z* _  U  m3 Bfor k = 0:119
' y! ~3 c/ Q- Z5 }  ax.View = [k, el];
  drawnow;
  frame = getframe(1);
  im = frame2im(frame);
  [A,map] = rgb2ind(im,256);
; Q9 ^. Y  P( }" ~1 m& ~$ d4 v3 K  if k == 0;
0 ?; \5 ]' |* h% f! E    imwrite(A,map,filename,'gif','LoopCount',Inf,'DelayTime',0);
  else
8 _4 F8 @" m: A% @5 W    imwrite(A,map,filename,'gif','WriteMode','append','DelayTime',0);
  end
end

duoer

看看楼主的代码。