Matlab Mobius带边的坐标

这是一个生成Mobius带的倍频程/Matlab代码

u = linspace(0,2*pi,100);
v = linspace(-1.0,1.0,100);
[u,v] = meshgrid(u,v);

x = (1+v.*cos(u/2)).*cos(u);
y = (1+v.*cos(u/2)).*sin(u);
z = v.*sin(u/2);

plot3(x,y,z)

输出如下。 在这个条带中，我需要边坐标（XYZ）。如何获得边的XYZ坐标

---

我可以使用python实现这一点，如下所示。参见相关帖子

---

什么是“边缘坐标”？XYZ？Mobius带的第一列/最后一列/行只有一个ege。我不知道如何提取边缘坐标。谢谢。但是，我必须得到这里给出的环的扭曲结构。

plot3(x([1 end],:).',y([1 end],:).',z([1 end],:).', "b")

bLength=1.6
numPoints=10
radius = bLength*numPoints / (2 * np.pi)
theta = np.linspace(0,2*np.pi,numPoints,endpoint=False)
dtheta=theta[1]-theta[0]

x0,y0=(radius * np.cos(theta)), (radius * np.sin(theta))
x1,y1=(radius * np.cos(theta+dtheta/2)) , (radius * np.sin(theta+dtheta/2))
cons0=np.ones(x0.shape)*0
cons1=np.ones(x1.shape)*2

np.savetxt('cooRing00.csv',np.c_[x0,y0,cons0],delimiter=' ',fmt='%10f')
np.savetxt('cooRing01.csv',np.c_[x1,y1,cons1],delimiter=' ',fmt='%10f')