求助：如何画三个变量的隐函数图形？

arieshin · 发表于 2006-6-13 07:47

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

<v:shapetype> -Z^2+3X^2+Y^2-8XY=6的图形应该怎么画呢？ 多谢了！</v:shapetype>

[此贴子已经被ericlin于2006-6-13 13:08:16编辑过]

ericlin · 发表于 2006-6-13 13:13

X=1:100;Y=1:100; Z=-sqrt(6-3*X.^2-Y.^2+8.*X.*Y) plot3(X,Y,Z) 我也是外行，不知道这样行不？

yqchenlolo · 发表于 2006-6-13 14:56

% 画三维隐式曲面f(x,y,z)=0，fun为函数表达式，s为变量名称列表（默认为x,y,z），range_xyz为x,y,z范围,ngrid为网格剖分数 function test=implot3(fun,range_xyz,ngrid,s) if nargin==0 %无参数时，默认在区域[-1,1, -1,1,-1,1]画单位园,网格剖分为20 fun=inline('x^2+y^2+z^2-1','x','y','z'); s=[sym('x'),sym('y'),sym('z')]; range_xyz=[-1,1,-1,1,-1,1]; ngrid=20; end if nargin==1 s=[sym('x'),sym('y'),sym('z')]; range_xyz=[-1,1,-1,1,-1,1]; ngrid=20; end if nargin==2 s=[sym('x'),sym('y'),sym('z')]; ngrid=20; end if nargin==3 s=[sym('x'),sym('y'),sym('z')]; end xm=linspace(range_xyz(1),range_xyz(2),ngrid); ym=linspace(range_xyz(3),range_xyz(4),ngrid); zm=linspace(range_xyz(5),range_xyz(6),ngrid); [xx,yy,zz]=meshgrid(xm,ym,zm); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%% inline形式的可通过sym(char(fun))化为sym形式的，实际上以下这段代 %%%%%%%%%%%%%%%% 码可以省去， %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% if isa(fun,'inline')%%%判断fun是否inline型,是的话返回1，否则返回0 c=findsym(sym(char(fun))); for k=1:3:length(c); w1(1+(k-1)/3)=sym(c(k)); end f=vectorize(fun);%vectorize函数举例：fun=inline('x^2+y^2+z^2-1','x','y','z')，则f=x.^2+y.^2+z.^2-1，即在^前加了个“.”号 if length(w1)==3 %%如果变量个数为3时执行的命令 fvalues=feval(f,xx,yy,zz); elseif length(w1)==2 %%如果变量个数为2时执行的命令 if char(w1(1))==char(s(1))&char(w1(2))==char(s(2)) fvalues=feval(f,xx,yy); end if char(w1(1))==char(s(1))&char(w1(2))==char(s(3)) fvalues=feval(f,xx,zz); end if char(w1(1))==char(s(2))&char(w1(2))==char(s(3)) fvalues=feval(f,yy,zz); end elseif length(w)==1 %%如果变量个数为1时执行的命令 if char(w1(1))==char(s(1)) fvalues=feval(f,xx); end if char(w1(1))==char(s(2)) fvalues=feval(f,yy); end if char(w1(1))==char(s(3)) fvalues=feval(f,zz); end else test=length(w); return;%error('error'); end end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% if isa(fun,'sym')%%%判断fun是否sym型,是的话返回1，否则返回0 r=findsym(fun); %%%%找到sym型的fun中变量的个数，返回的r为char类型，所以之后的一个for循环将char转化为sym for k=1:3:length(r) w(1+(k-1)/3)=sym(r(k));%因为如果有三个变量，则r=[x, y, z],即有7个元素：r(1)=x,r(2)=逗号,r(3)=空格,........ end if length(w)==3 %%如果变量个数为3时执行的命令 fvalues=subs(fun,w,{xx,yy,zz}); elseif length(w)==2 %%如果变量个数为2时执行的命令 if char(w(1))==char(s(1))&char(w(2))==char(s(2)) fvalues=subs(fun,w,{xx,yy}); end if char(w(1))==char(s(1))&char(w(2))==char(s(3)) fvalues=subs(fun,w,{xx,zz}); end if char(w(1))==char(s(2))&char(w(2))==char(s(3)) fvalues=subs(fun,w,{yy,zz}); end elseif length(w)==1 %%如果变量个数为1时执行的命令 if char(w(1))==char(s(1)) fvalues=subs(fun,w,{xx}); end if char(w(1))==char(s(2)) fvalues=subs(fun,w,{yy}); end if char(w(1))==char(s(3)) fvalues=subs(fun,w,{zz}); end else test=length(w); return;%error('error'); end end isosurface(xx,yy,zz,fvalues,0); %p=isosurface(x,y,z,fvalues,0); %set(p, 'FaceColor', 'interp', 'EdgeColor', 'red'); xlabel('x'); ylabel('y'); zlabel('z'); grid; 写的有点罗唆了，不过可以用 >>fun=inline('-z^2+3*x^2+y^2-8*x*y-6') fun =
 Inline function: fun(x,y,z) = -z^2+3*x^2+y^2-8*x*y-6 >>implot3(fun,[-2*pi,2*pi,-2*pi,2*pi,-2*pi,2*pi])

yqchenlolo · 发表于 2006-6-13 14:57

在6.5中运行的 7.0我没用过，好像有命令直接可以画 你问问别人

suffer · 发表于 2006-6-13 16:08

<DIV class=quote>以下是引用yqchenlolo在2006-6-13 14:57:33的发言： 在6.5中运行的 7.0我没用过，好像有命令直接可以画 你问问别人</DIV>
 matlab7也只能用ezplot画2维的，三维的隐函数没有对应的命令可以用

账号		自动登录	找回密码
密码			我要加入

[绘图技巧] 求助：如何画三个变量的隐函数图形？

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（yqchenlolo）在6.5中运行的7.0我没用过，好像...