请教：关于二维非规则区域等值云图的matlab实现方法

pblossom

比如画一个L型的区域，如下图，我现在想画该区域的温度云图，我尝试分成两块，用两次contourf函数，并画出等值线，其中会有一黑色的分块线，不太好。请问，还有其他方法实现么？

messenger

先画一个contour图，再叠加一个白的图不行吗？

可能有命令可以直接去黑线，不过没法试。问问题最好把代码也贴上来，有代码顺便就帮你试了。否则，想画一个L型的等值线图也得费一番功夫，而且手头也没有现成的数据。

pacoyang

2# messenger
我试了一下，图形save as后怎么白色区域没有盖住等值线呢？没有保存前是盖住的，不保存直接copy figure到word里白色区域称黑色的了。

jinzi2006

这样行不行，先画整个区域的云图，然后用白色填充小区域，pcolor。

messenger

你用什么命令画的白色区域没有盖住呀？

2# messenger
我试了一下，图形save as后怎么白色区域没有盖住等值线呢？没有保存前是盖住的，不保存直接copy figure到word里白色区域称黑色的了。
pacoyang 发表于 2009-5-3 08:33

pacoyang

5# messenger
用snagit截图没问题

1. 
   
2. clear
   
3. clc
   
4. %
   
5. Z = peaks;
   
6. [C,h] = contour(Z);
   
7. set(h,'Fill','on');
   
8. hold on
   
9. x=[5 20 20 10 10 5]; y=[5 5 10 10 30 30];
   
10. fill(x,y,'w','FaceColor','w','EdgeColor','w');
    
11. axis equal; axis off;
    

复制代码

messenger

6# pacoyang

很奇怪的问题，如果用contour画等值图，覆盖后保存，空白图即为黑色。

然而，如果用contourf画等值图，覆盖后保存，空白图即为原色。

pblossom

非常感谢，有这么多人帮忙，很有启发。附上matlab代码，大家可以试试

1. 
   
2. % 二维L型区域有源项导热计算
   
3. % 两边第三类边界条件，对称面绝热条件
   
4. % aP、b表达式中均乘以完整控制容积，
   
5. % 故边界节点源项应乘上其各自非完整度0.25、0.5、0.75
   
6. clear;
   
7. err=0.0001;
   
8. deltay=0.1; deltax=0.1;
   
9. dy=0.05; dx=0.05;
   
10. 
    
11. M=45; N=30;
    
12. lamda=50;
    
13. 
    
14. T=zeros(N+3,M+3);
    
15. Sc=zeros(N+1,M+1);
    
16. Sp=zeros(N+1,M+1);
    
17. 
    
18. Spad=-20*deltay/0.01;  % 外边界 第三类 h=20
    
19. Scad=20*25*deltay/0.01;  % t=25  deltax=deltay
    
20. % Spad=-10^30;  % 外边界 第一类 t=20
    
21. % Scad=10^30*20;  
    
22. 
    
23. Sc(2:N,1)=Sc(2:N,1)+Scad*ones(N-1,1);
    
24. Sp(2:N,1)=Sp(2:N,1)+Spad*ones(N-1,1);
    
25. Sc(N+1,2:M)=Sc(N+1,2:M)+Scad*ones(1,M-1);
    
26. Sp(N+1,2:M)=Sp(N+1,2:M)+Spad*ones(1,M-1);
    
27. 
    
28. Spad=-20*deltax/0.01;  % 内边界 第三类 h=20
    
29. Scad=20*200*deltax/0.01;  % t=200 deltax=deltay
    
30. % Spad=-10^30;  % 内边界 第一类 t=120
    
31. % Scad=10^30*120;  
    
32. 
    
33. Sc(2:10,31)=Sc(2:10,31)+Scad*ones(9,1);
    
34. Sp(2:10,31)=Sp(2:10,31)+Spad*ones(9,1);
    
35. Sc(11,32:M)=Sc(11,32:M)+Scad*ones(1,M-31);
    
36. Sp(11,32:M)=Sp(11,32:M)+Spad*ones(1,M-31);
    
37. 
    
38. Spad=-20*dx/0.01;
    
39. Scad1=20*25*dy/0.01;
    
40. Scad2=20*200*dx/0.01;
    
41. % Scad1=10^30*20;
    
42. % Scad2=10^30*120;
    
43. 
    
44. Sc(1,1)=Sc(1,1)+Scad1;
    
45. Sp(1,1)=Sp(1,1)+Spad;
    
46. Sc(1,31)=Sc(1,31)+Scad2;
    
47. Sp(1,31)=Sp(1,31)+Spad;
    
48. Sc(11,31)=Sc(11,31)+2*Scad2;
    
49. Sp(11,31)=Sp(11,31)+2*Spad;
    
50. Sc(11,M+1)=Sc(11,M+1)+Scad2;
    
51. Sp(11,M+1)=Sp(11,M+1)+Spad;
    
52. Sc(N+1,1)=Sc(N+1,1)+2*Scad1;
    
53. Sp(N+1,1)=Sp(N+1,1)+2*Spad;
    
54. Sc(N+1,M+1)=Sc(N+1,M+1)+Scad1;
    
55. Sp(N+1,M+1)=Sp(N+1,M+1)+Spad;
    
56. 
    
57. while 1
    
58.     T0=T;
    
59.     for i=1:31
    
60.         for j=1:31
    
61.             if i==1
    
62.                 aS=0.0;
    
63.             else
    
64.                 aS=lamda;
    
65.             end
    
66.             if j==1
    
67.                 aW=0.0;
    
68.             else
    
69.                 aW=lamda;
    
70.             end
    
71.             if i==31
    
72.                 aN=0.0;
    
73.             else
    
74.                 aN=lamda;
    
75.             end
    
76.             if j==31 && i<11
    
77.                 aE=0.0;
    
78.             else
    
79.                 aE=lamda;
    
80.             end
    
81.             aP=aE+aW+aN+aS-Sp(i,j)*0.01;
    
82.             b=Sc(i,j)*0.01;
    
83.             T(i+1,j+1)=(aE*T(i+1,j+2)+aW*T(i+1,j)+aN*T(i+2,j+1)+...
    
84.                 aS*T(i,j+1)+b)/aP;
    
85.         end
    
86.     end
    
87.     for j=32:M+1
    
88.         for i=11:N+1
    
89.             if i==11
    
90.                 aS=0.0;
    
91.             else
    
92.                 aS=lamda;
    
93.             end
    
94.             if i==N+1
    
95.                 aN=0.0;
    
96.             else
    
97.                 aN=lamda;
    
98.             end
    
99.             if j==M+1
    
100.                 aE=0.0;
     
101.             else
     
102.                 aE=lamda;
     
103.             end
     
104.             aP=aE+aW+aN+aS-Sp(i,j)*0.01;
     
105.             b=Sc(i,j)*0.01;
     
106.             T(i+1,j+1)=(aE*T(i+1,j+2)+aW*T(i+1,j)+aN*T(i+2,j+1)+...
     
107.                 aS*T(i,j+1)+b)/aP;
     
108.         end
     
109.     end
     
110.     if norm(T-T0,inf) < err
     
111.         break
     
112.     end
     
113. end
     
114. 
     
115. clf;
     
116. [X,Y]=meshgrid(0:0.1:4.5,0:0.1:3.0);
     
117. [C,h,CF]=contourf(X,Y,T(2:N+2,2:M+2));
     
118. clabel(C,h)
     

复制代码

最后需要改动下，我copy错文件了