分享一个猴鞍面的代码

panphp

分享一个猴鞍面的代码：
/ M# N& ^0 ]) J. J
4 i8 E7 d- K; w6 m
0 P( x  M( H3 F) I% monkey_saddle.m
$ K- o, s) J& m- |, L%
. [6 k) ~0 g# L) P% D% Make an animation of a monkey saddle suRFace.
, ?) c) [0 Z% A) m3 z1 Q
: T% K1 H: J+ i

% Name of animation created
) c/ n; Q: W# }4 o/ O! g, y* f  Gfilename = 'monkey_saddle.gif';

9 F6 q8 r' @. J$ k" N( {/ y6 }6 ^* N% Polar coordinates
r = linspace(0, 1);
# g  C9 `- `/ N8 q9 l- ~t = linspace(-pi, pi);
[R, T] = meshgrid(r, t);

' T6 Z) t% V; _. i5 n0 L% Cartesian coordinates
x = R .* cos(T);
y = R .* sin(T);
. x( W  V4 O# T. G1 Dz = x.^3 - 3 * x .* y.^2;

% Gaussian curvature to color the surface
" J. `/ c; ]+ B2 P& H' f8 f0 vK = - 36 * (x.^2 + y.^2)...
  ./ (1 + 9 * (x.^2 + y.^2).^2).^2;

% Plot the surface
5 L% Y* r- U0 J$ z. `& _chartSurface = surf(x, y, z, K);

: k; d/ U! ^" h+ n: ]% Edit its appearance
3 Z) u6 U( S8 p' zchartSurface.EdgeColor    = 'none';
( ~, R6 F5 s. u/ t# i( U1 k) fchartSurface.FaceColor    = 'interp';
chartSurface.FaceLighting = 'gouraud';
0 Y! t# I) y. e5 D5 l9 Q, {
ax = gca;

$ S, X6 ]% m2 j; N: X% Remove tick labels
9 {( P9 G+ I$ L. K: q2 F& Q6 Lax.XTickLabel = {};
6 M: L* I( G8 b: f8 v4 Hax.YTickLabel = {};
ax.ZTickLabel = {};
7 |2 d. }# d: N& o3 C" m* l" O$ g4 x
6 k' V; v$ X. A! @2 n9 _4 l% Rotate the camera's azimuthal angle
! U: |- k: ^! qaxis vis3d;
+ B& H9 l- z2 Lv = ax.View;
% b/ O9 t" Y$ x7 g8 i: |az = v(1);
. [) y  O9 w+ N; ^0 d/ x9 iel = v(2);
) _9 U& S. W/ i* u+ e' M5 Z/ t( ~! V- l
: M( G4 f  o" ^% Minimize file size by taking advantage of 3-fold symmetry of the surface
for k = 0:119
  ax.View = [k, el];
9 R, T& Q4 n! A, z1 Y- O: M* k  drawnow;
* V$ t6 o5 G% k, a8 Z5 G  frame = getframe(1);
  im = frame2im(frame);
: N. X  T! Q9 _* ?1 T: p  [A,map] = rgb2ind(im,256);
' b8 g' v9 ~) d9 H% C3 I4 F& K  if k == 0;
, E; p" m7 L, m7 d1 p; |    imwrite(A,map,filename,'gif','LoopCount',Inf,'DelayTime',0);
  z+ z2 b3 O1 x+ W5 @' I  else
    imwrite(A,map,filename,'gif','WriteMode','append','DelayTime',0);
  end
end
& T9 d8 ^! |) x5 p* A$ z& C

duoer

看看楼主的代码。