如何画非矩形域的曲面

ChaChing

整理老帖发现此帖http://forum.vibunion.com/forum/viewthread.php?tid=60064&highlight=
这个想了下也搜索下, 怎好像没现成的!?
私心为自己日后搜索方便, 以新帖与大家分享!

水平有限, 仅有两个思路
1.仍使用矩形域, 但区域外使用nan取代, 但会有锯齿状
2.将矩形域座标转换下
举个简易例子, 不知是否有更好方式, 请高手指教!

1. 
   
2. clc; clear
   
3. [X1,Y1] = meshgrid(-2:.1:2, -2:.1:2);
   
4. Z1 = X1 .* exp(-X1.^2 - Y1.^2); surf(X1,Y1,Z1)
   
5. figure; Z3=Z1;
   
6. for ii=1:10
   
7.    Z3(ii,1:11-ii)=nan; Z3(ii,32+ii:41)=nan;
   
8.    Z3(42-ii,1:11-ii)=nan; Z3(42-ii,32+ii:41)=nan;
   
9. end; surf(X1,Y1,Z1)
   
10. figure;
    
11. [th,r] = meshgrid((0:1/80:1)*2*pi,0:.1:2); [X2,Y2] = pol2cart(th,r);
    
12. Z2 = X2 .* exp(-X2.^2 - Y2.^2); surf(X2,Y2,Z2)
    

复制代码

[ 本帖最后由 ChaChing 于 2010-8-3 00:59 编辑 ]

ChaChing

刚才发现有点贴错, 又没权限编辑了! 所有重发并完善下
注意仅是提供思路, 并非通用化(或许也可行, 只是...)

1. clc; clear; [X1,Y1] = meshgrid(-2:.1:2, -2:.1:2); Z1 = X1 .* exp(-X1.^2 - Y1.^2);

复制代码

1. surf(X1,Y1,Z1)
   

复制代码

1. 
   
2. Z2=Z1;
   
3. for ii=1:20, Z2(ii,1:21-ii)=nan; Z2(ii,22+ii:41)=nan; Z2(42-ii,1:21-ii)=nan; Z2(42-ii,22+ii:41)=nan; end;
   
4. surf(X1,Y1,Z2)
   

复制代码

1. 
   
2. [th,r] = meshgrid((0:1/80:1)*2*pi,0:.1:2); [X3,Y3] = pol2cart(th,r); Z3 = X3 .* exp(-X3.^2 - Y3.^2);
   
3. surf(X3,Y3,Z3)
   

复制代码

feier0mei0

哇，多好的帖子，怎么没人回复呢。但是可能不符合我的要求，我要画的三维曲线图的数据都是离散的，相互之间没关系的。

ChaChing

回复 3 # feier0mei0 的帖子

help griddata

86726803

好帖子！