如何在圆周方向上画三维图？

valiant

假设存在“L”型的物体，旋转一周即可得到一个“杯子”似的圆形物体。
如果我的“L”型只分布在几个角度（比如0°，90°，180°，270°），怎么插值画出这个“杯子”似的三维图？我用mesh，surf试了几次也都没弄出来。。
假设“L”型的离散数据如下：
x=[10 20 30 40 10 20 30 40 30 40 28 40];
y=[10 10 10 10 20 21 19 20 30 30 40 40];
t=[100 100 100 100 500 500 500 500 500 100 500 100 ];
我先前问了一个二维的画图，但是不知道怎么改成三维的了。
又来求助。希望好心人能帮助我。

qibbxxt

这个我也不是很明白，我觉得可以用旋转体来画图
写了以下的代码，但是感觉不大对，供你参考

1. clear;clc;close all
   
2. t=[100 100 100 100 500 500 500 500 500 100 500 100 ];
   
3. [X,Y,Z] = cylinder(t,100);
   
4. surf(X,Y,Z)
   
5. axis square
   
6. shading interp

复制代码

玉林

可以先建立一个柱坐标系么，cart2pol,直角坐标系转换成柱坐标系。

Happy99

先自行搜索一下吧! 印象中看过类似的:@)
记的不难, 笔电出问题了, 还没matlab可试:@(

雪缘

在柱坐标下，给出一系列的角度，然后获得相应的坐标，最后用surf画估计可行

valiant

1. x=[10 20 30 40 10 20 30 40 30 40 28 40];
   
2. z=[10 10 10 10 20 21 19 20 30 30 40 40];
   
3. t=[100 100 100 100 500 500 500 500 500 100 500 100 ];
   
4. %假设存在0°90°180°270°四个方向，且数据一样。
   
5. x1=x*cosd(0);y1=x*sind(0);
   
6. x2=x*cosd(90);y2=x*sind(90);
   
7. x3=x*cosd(180);y3=x*sind(180);
   
8. x4=x*cosd(270);y4=x*sind(270);
   
9. %end
   
10. x=[x1;x2;x3;x4];
    
11. y=[y1;y2;y3;y4];
    
12. z=[z;z;z;z];
    
13. t=[t;t;t;t];
    
14. p1=-max(x1):1:max(x1); %三维坐标下x方向
    
15. p2=-max(x1):1:max(x1); %三维坐标下y方向
    
16. p3=0:1:max(z); %三维坐标下z方向
    
17. [px,py] = meshgrid(p1,p2);
    
18. F = TriScatteredInterp(x,y,t);
    
19. pt=F(px,py);
    
20. pt((pt>500))=NaN;pt((pt<10))=NaN;
    
21. mesh(px,py,pt);shading interp; colorbar;

复制代码

还是不行，请各位高手明示啊。。

Happy99

Ref:
http://forum.vibunion.com/forum- ... t-%D0%FD%D7%AA.html
http://forum.vibunion.com/forum- ... t-%D0%FD%D7%AA.html
http://forum.vibunion.com/forum- ... t-%D0%FD%D7%AA.html
...

ChaChing

回复 6 # valiant 的帖子

???
F = TriScatteredInterp(x(:),y(:),t(:));