在matlab

Question

我正在开展一个项目，我们以线性方式进行测量并旋转样本以生成3D表面图。

数据来自极地格式，我们有4个数据集：

theta = 0  , r = -20 -- 20
theta = 45 , r = -20 -- 20
theta = 90 , r = -20 -- 20
theta = 135, r = -20 -- 20

我可以通过生成带有sin和cos的x和y向量并将结果绘制成该矩阵来将其转换为笛卡尔点。

所以“准备绘制”矩阵集的示例看起来像这样：

x [4x21 double]
-20    -18    -16 ..
-14.14 -12.73 -11.31 ..
0      0      0 ..
..

y [4x21 double]
0      0      0 ..
-14.14 -12.73 -11.31 ..
-20    -18    -16 ..
..

z [4x21 double]
.224 .281 .320 ..
.228 .280 .312 ..
.282 .349 .377 ..
..

带有一些基本格式的输出表面图如下所示：

Matlab连接我想要的向量。但是，它没有将90度旋转连接到135度旋转。我如何做到这一点？

产生类似于上面例子的数字的基本代码是：

theta = [ 0 , -45 , -90 , -135];
r = -20 : 2 : 20;
for i = 1 : 4
    xc(i,:) = r .* cosd(theta(i));
    yc(i,:) = r .* sind(theta(i));
    zc(i,:) = [[0 : 0.1 : 0.9] , [1 : -0.1 : 0]];
end
figure
surf(xc,yc,zc);

Answer 1

您必须明确地theta=-180。就像plot没有环绕并加入第一个和最后一个点一样，surf不会将第一个“边缘”连接到最后一个“边缘”。制作theta=[0 -45 -90 -135 -180]并执行for i=1:5，这应该是好的。

您会发现着色看起来有点奇怪，根据应用程序的不同，这可能看起来更好：surf(xc,yc,zc,'FaceColor','interp')。

创建数据的更好方法可能是

[R,TH]=meshgrid(r,theta);
[xc,yc]=pol2cart(TH*pi/180,R);
zc = repmat([[0 : 0.1 : 0.9] , [1 : -0.1 : 0]],length(theta),1);