Question

MATLAB中用ezplot3将两个图画在同一个figure中

如题，MATLAB中为什么不能用这种方法将两个图叠加到同一个figure中呢？
ezplot3(x,y,z)
hold on
ezplot3(xp,yp,zp)
clc,clear all

syms a1
k1=0.08;
n1=7;n2=8;
n=n2/n1;
i12=(1+sqrt(n^2*(1-k1^2)+k1^2))*k1*cos(n1*a1)/(k1^2-1)+(k1^2+sqrt(n^2*(1-k1^2)+k1^2))/(1-k1^2);
deta1=acot(i12);
deta2=pi/2-deta1;
R=37.64;
a2=int(1/i12);

x=R.*sin(deta1).*cos(a1);
y=R.*sin(deta1).*sin(a1);
z=R.*cos(deta1);
xp=-R.*cos(deta2);
yp=R.*sin(deta2).*sin(a2);
zp=R.*sin(deta2).*cos(a2);

ezplot3(x,y,z,[0,2*pi]);
hold on
p=ezplot3(xp,yp,zp,[0,4*pi]);
set(p,'color','r')

title('')

Answer 1

在后面加一句axis auto试试。
ezplot3会根据最后一次绘图来确定坐标范围，有可能导致前面绘图不在显示范围内。

追问

在上面补充了一个程序，能不能帮我看一下

追答

把

p=ezplot3(xp,yp,zp,[0,4*pi]);

改成

p=ezplot3(char(xp),char(yp),char(zp),[0,4*pi]);

前者调用的是对sym对象重载的版本，不知出于什么考虑，前面加了一句clf，而且未返回句柄（2008a是这样，不知道后来是否改过来了）。

追问

加上去果然可以了。
不过我不明白的是，为什么要加上字符矩阵呢？x,y,z以及xp,yp,zp都是关于a1的参数方程，理论上ezplot3(xp,yp,zp,[0,4*pi])这种作图是合法的吧

追答

的确是合法的，但调用的函数版本不同。这个函数有两个版本：

>> which ezplot3 -all
C:\MATLAB_R2014a\toolbox\matlab\specgraph\ezplot3.m
C:\MATLAB_R2014a\toolbox\symbolic\symbolic\@sym\ezplot3.m  % sym method

 

你的写法调用的是针对sym对象重载的版本，而这个重载版本在部分MATLAB版本中有BUG（就我查到的而言，2008a的这个函数一是多了个莫名其妙的clf，导致不能叠加绘图，二是不返回绘图句柄。2012a和2014a都修正了这两个BUG）

 

我的写法直接调用针对字符串对象的普通版本，从而避开了上述BUG。事实上，sym重载的版本也是先把sym对象转为字符串，然后再调用字符串版本的ezplot3。

Answer 2

可以啊，你出现什么问题呢

追问

帮我看一下上面的程序吧

追答

这个查询了下，貌似是Bug，你把参数初始化吧，然后用plot3对值进行作图吧，这时hold on是有效的。
你也可以查询Google，利用"ezplot3 hold on"关键词，很多人遇到同样的问题，目前没有看到解决办法。