画出三维曲线z^2=x^3+6*x*y+y^4+6*x+2*y-1 的图像.用matlab绘制
来源:学生作业帮助网 编辑:六六作业网 时间:2025/01/13 04:40:50
画出三维曲线z^2=x^3+6*x*y+y^4+6*x+2*y-1 的图像.用matlab绘制
画出三维曲线z^2=x^3+6*x*y+y^4+6*x+2*y-1 的图像.用matlab绘制
画出三维曲线z^2=x^3+6*x*y+y^4+6*x+2*y-1 的图像.用matlab绘制
f=@(x,y,z)x.^3+6*x.*y+y.^4+6*x+2*y-1-z.^2;
h=implicitsurf(f,[-5 5]);
----------------------------------------------
function h=implicitsurf(f,xlimit,ylimit,zlimit,gd)
%implicitsurf(f,span,gd):画隐函数曲面f(x,y,z)=0的网格图,
% 各坐标范围均限定在span=[lb,ub],
% 网格数为gd,默认为25
%implicitsurf(f,xspan,yspan,zspan,gd):画隐函数曲面f(x,y,z)=0,
% 各坐标范围分别限定在xspan,yspan,zspan
%h=implicitsurf(...):画隐函数曲面并输出句柄
%例一:
%implicitsurf(inline('x.*y+z.^2'),[-5 5])%注意*\^一定要设成点运算
if nargin==2
ylimit=xlimit;zlimit=xlimit;gd=25;
elseif nargin==3
gd=ylimit;ylimit=xlimit;zlimit=xlimit;
elseif nargin==4
gd=25;
elseif nargin==5
else
error('Error in input arguments')
end
x=linspace(xlimit(1),xlimit(2),gd);
y=linspace(ylimit(1),ylimit(2),gd);
z=linspace(zlimit(1),zlimit(2),gd);
[x,y,z]=meshgrid(x,y,z);val=f(x,y,z);
[f,v]=isosurface(x,y,z,val,0);
if isempty(f)
warning('There is no graph in the range.');
p=[];
else
newplot;
p=patch('Faces',f,'Vertices',v,'CData',v(:,3),'facecolor','flat','EdgeColor','k');
isonormals(x,y,z,val,p);view(3);grid on
end
if nargout==0
else
h=p;
end