如何用matlab生成随机的封闭曲面？

charles792

我知道利用圆的参数方程和一维插值方法可以生成随机的封闭曲线。

取若干个【0,2*pi】内的数作为参数a，并对应生成一个随机数作为半径r（两端的值相等）。
再通过插值方法将【0，2*pi】内“所有”点的半径求出来，最后作图即可。

但是如何利用球体的参数方程和二维插值方法来生成随机封闭曲面呢？
参数a【0，2*pi】
参数b【-0.5*pi，0.5*pi】

出现问题是这曲线做出来不封闭且两头尖啊，估计是插值时出了问题，求高手指点……

1. 
   
2. N=0.1;
   
3. n=pi/6;
   
4. nonrandbit=1;
   
5. randbit=1;
   
6. ta=[0:n:2*pi];
   
7. tb=[-0.5*pi:n:0.5*pi];
   
8. [Ta,Tb]=meshgrid(ta,tb);
   
9. tta=[0:N:2*pi];
   
10. ttb=[-0.5*pi:N:0.5*pi];
    
11. [TTa,TTb]=meshgrid(tta,ttb);
    
12. r=nonrandbit+randbit*rand(size(Ta));
    
13. %r=nonrandbit+randbit*ones(size(t));
    
14. r(:,size(r,2))=r(:,1);
    
15. x=r.*cos(Ta).*cos(Tb);
    
16. y=r.*sin(Ta).*cos(Tb);
    
17. z=r.*sin(Tb);
    
18. rr = interp2(Ta,Tb,r,TTa,TTb);
    
19. xx=rr.*cos(TTa).*cos(TTb);
    
20. yy=rr.*sin(TTa).*cos(TTb);
    
21. zz=rr.*sin(TTb);
    
22. 
    
23. subplot(221)
    
24. mesh(x,y,z)
    
25. subplot(222)
    
26. mesh(xx,yy,zz)
    
27. subplot(223)
    
28. plot3(x,y,z,'*')
    
29. subplot(224)
    
30. plot3(xx,yy,zz,'*')
    

复制代码

charles792

顶…………………………………………