据说这个函数能画出一个立体的心

broken1999

  |8 F- @; W5 M0 d# {我的电脑太慢 画不出来   哪个画出来了上传fig文件吧

(x^2 + (9/4)*y^2 + z^2 - 1)^3 - x^2*z^3 - (9/80)*y^2*z^3;

$ r% G7 b  B$ z传说中的效果图：
( E, p. V" Q5 I3 i# w1 p
+ J! F2 a! E9 t. V* V

嫩草一棵

对了，顺便提一下，这个是用Mathematica画的，所以没有fig文件，是一个nb文件。
heart.nb (829.61 KB, 下载次数: 0)

2012-8-26 10:44 上传

点击文件名下载附件
下载积分: 体力 -2 点

上面就是文件。
* e6 a. @2 a, [" V3 b顺便求大神解释代码。。。

嫩草一棵

已经试过了，很好啊。谢楼主分享。

zz1990826

zz1990826

大家一起试试看

netabc2008

??? Error using ==> mtimes
1 F: y2 G% E+ w. `" |- }' [% z! `7 dInput arguments must be 2-D.
5 t9 L3 e" \5 ^3 r, t6 N- C5 r
% b; O' [; [) T/ ~Error in ==> heart>@(x,y,z)(x.^2+(9/4)*y.^2+z.^2-1).^3-x.^2*z.^3-(9/80)*y.^2.*z.^3 at 4
f1=@(x,y,z)(x.^2+(9/4)*y.^2+z.^2-1).^3-x.^2*z.^3-(9/80)*y.^2.*z.^3;

' @( i0 Z8 J1 k. a- i# VError in ==> heart at 6
2 t6 y0 z+ \, y/ m+ p4 f: L0 qset(patch(isosurface(x,y,z,f1(x,y,z),0)),'facecolor','r','edgecolor','none');
6 y0 g* m- ^( i4 w
运行不出来?!!

netabc2008

set(patch(isosurface(x,y,z,f1(x,y,z),0)),'facecolor','r','edgecolor','none');
??? Error using ==> mtimes
Input arguments must be 2-D.

broken1999

回复 madio 的帖子

( L8 E2 K' j- @7 ~$ {$ V0 Y! G  A
( }& H- H! ^  `+ {# a* j" L    佩服佩服

madio

1. clear8 f, l+ \: [* D0 Q6 ]7 w! H
   
2. N=100;
   - R7 y4 D( t  J8 c! `8 y
3. f1=@(x,y,z)(x.^2+(9/4)*y.^2+z.^2-1).^3-x.^2.*z.^3-(9/80)*y.^2.*z.^3;
   7 p3 c( y( ?; b. E; a- D: A5 _
4. [X,Y,Z]=meshgrid(linspace(-1.5,1.5,N));6 X. o$ R8 U3 s+ I
   
5. set(patch(isosurface(X,Y,Z,f1(X,Y,Z),0)),'facecolor','r','edgecolor','none');8 N1 s$ f& R. d# X2 R
   
6. light* l' n. s; B  ?/ d! B  U0 S
   
7. grid on
   . Q- ?- t8 e- ]. u: T
8. view([1,4,1])

复制代码

a4441021

clear
; J! `; D1 c( }: ^9 _- u5 R5 QN=100;
f1=@(x,y,z)(x.^2+(9/4)*y.^2+z.^2-1).^3-x.^2.*z.^3-(9/80)*y.^2.*z.^3;
3 ~1 E2 c, |/ f( m9 P5 n[X,Y,Z]=meshgrid(linspace(-1.5,1.5,N));
9 P* Y* Y5 F4 R/ P9 r) V" Cset(patch(isosu**ce(X,Y,Z,f1(X,Y,Z),0)),'facecolor','r','edgecolor','none');
light
grid on
. a! Y9 D  l; x) T0 ~3 E, Rview([1,4,1])