如何画多值函数？怎么画最方便？

xueyi

比如下面一个函数，x，y，z为未知数
(x-y)^2+(y-z)^2+(z-y)^2=10

xueyi

!!

xueyi

谢谢：）

xueyi

试着做了一下，报错
??? Error using ==> surface
X, Y, Z, and C cannot be complex.

z求出来的是复数，画不出来
我去abs以后能画出来，但我觉得应该是个椭球，怎么不是呢？

bainhome

用第六行语句看看x和y在什么范围内能保证根号下为正，然后再画！

1. function EzFunctionPlot
   
2. clc
   
3. D=solve('(x-y)^2+(y-z)^2+(z-x)^2=10','z')
   
4. a1=vectorize('1/2*y+1/2*x+1/2*(-3*y^2+6*x*y-3*x^2+20)^(1/2)');
   
5. a2=vectorize('1/2*y+1/2*x-1/2*(-3*y^2+6*x*y-3*x^2+20)^(1/2)');
   
6. % ezplot('-3*y^2+6*x*y-3*x^2+20')  %确定范围
   
7. [x,y]=meshgrid(3:.05:12);
   
8. z1=inline(a1,'x','y');
   
9. z=z1(x,y);
   
10. surf(x,y,z)
    
11. hold on
    
12. z2=inline(a2,'x','y');
    
13. z=z2(x,y);
    
14. surf(x,y,z)
    
15. shading flat

复制代码

其中的两个切平面是取值自动带进去的，应当可以寻找某种函数规则，将不在椭球面上的点全部赋值为NaN裁掉

[ 本帖最后由 bainhome 于 2006-7-10 14:12 编辑 ]