以下程序 x对应n1,y对应n2,z对应aob
采用的是等值面作图的思想
但是没有严格等于0的值,取0.01,仅供参考,
clear,clc
[x,y,z]=meshgrid(1:1:100,1:1:100,1:1:40);
v=sqrt((y+2).^2+(x+y).^2-2*(y+2).*(x+y).*cos(acos(y.*cos(z.*pi/180)./(y+2))...
-z.*pi/180))-x.*cos(z.*pi/180);
p=patch(isosurface(x,y,z,v,0.01));
isonormals(x,y,z,v,p)
set(p,'FaceColor','red','EdgeColor','none');
% daspect([1,1,1])
% view(3); axis tight
camlight
lighting gouraud
grid on