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

broken1999

  w4 R* r' B: k- l! \
6 [: E9 v" s9 t我的电脑太慢 画不出来   哪个画出来了上传fig文件吧
9 t+ K- ^5 Y) Y# ?; y# G/ l
5 i8 x6 }& s1 a4 h" s% c1 T(x^2 + (9/4)*y^2 + z^2 - 1)^3 - x^2*z^3 - (9/80)*y^2*z^3;
4 v6 f+ }4 k* b: S- C4 f- \
传说中的效果图：

4 i$ \( y. j& K) W
. A. M; v( @+ F' |, ?

嫩草一棵

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

2012-8-26 10:44 上传

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

上面就是文件。
顺便求大神解释代码。。。

嫩草一棵

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

zz1990826

zz1990826

大家一起试试看

netabc2008

??? Error using ==> mtimes
& N5 l9 a9 b1 z* fInput arguments must be 2-D.

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
6 G, _- `$ n1 R  {9 @f1=@(x,y,z)(x.^2+(9/4)*y.^2+z.^2-1).^3-x.^2*z.^3-(9/80)*y.^2.*z.^3;
1 W" n: x7 J  ^& v8 c: d5 i0 G
, s& h+ x2 [3 F8 IError in ==> heart at 6
/ w. M& a; }% g$ @$ pset(patch(isosurface(x,y,z,f1(x,y,z),0)),'facecolor','r','edgecolor','none');
5 _) h8 V7 f1 b3 [8 }1 S* O
. F$ O' n: 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 的帖子


    佩服佩服

madio

1. clear
   3 F! U5 S8 H+ P/ p3 Z
2. N=100;8 |  e9 B  R* {9 p' R; g3 _  {
   
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 z7 H% O+ b; _1 {
4. [X,Y,Z]=meshgrid(linspace(-1.5,1.5,N));
   . ?8 F. o) x8 l0 g5 w& a$ H3 m: U5 j8 V
5. set(patch(isosurface(X,Y,Z,f1(X,Y,Z),0)),'facecolor','r','edgecolor','none');
   * v. N0 @! h3 w  p# d! N2 K* @
6. light
   . i$ Y' Y6 |: Q4 d) [* M
7. grid on
   - ^3 t. z9 t. _. y& k+ }
8. view([1,4,1])

复制代码

a4441021

clear
6 q7 i7 y# k# ^+ J3 I% ~2 kN=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;
) I# C) F' }$ w& T[X,Y,Z]=meshgrid(linspace(-1.5,1.5,N));
set(patch(isosu**ce(X,Y,Z,f1(X,Y,Z),0)),'facecolor','r','edgecolor','none');
$ f% x9 i/ F  Q0 i! Flight
grid on
: t( v) c: o& J) ^. o0 Zview([1,4,1])