matlab画圆

bainhome

okok上贴了，但simwe的高手比较扎堆儿，所以再贴在这里，请高手指正，因为肯定不只一个思路。

1. function RandomPlotCircle(Lxmin,Lxmax,Lymin,Lymax,Rmin,Rmax,n,Factor)
2. & This Programm made by Bainhome
3. & Date:Mar.22th 2006
4. & First Edited in Mar.23th 2006
5. % **********************************  输  *  入  *  变  *  量  *  介  *  绍  **********************************
6. % Lmin,Lmax——The Bound of Circle Plot  Area  (Both X dir & Y dir)
7. % Rmin,Rmax——The Limitation of Radius
8. % Example：RandomPlotCircle(1,10,1,10,.1,2,10)
9. % *****************************************  完  *  毕  ***************************************************
10. close all
11. clc
12. count=0;
13. Lxmin=Lxmin+Rmax;
14. Lxmax=Lxmax-Rmax;
15. Lymin=Lymin+Rmax;
16. Lymax=Lymax-Rmax;
17. R1(1)=Rmin+rand*(Rmax-Rmin);
18. RT(1)=R1(1);
19. x1(1)=Lxmin+rand*(Lxmax-Lxmin);
20. y1(1)=Lymin+rand*(Lymax-Lymin);
21. for i=1:1
22.     if x1(i)+R1(i)<=Lxmax&x1(i)-R1(i)>=Lxmin&y1(i)+R1(i)<=Lymax&y1(i)-R1(i)>=Lymin
23.     else
24.         R1(i)=Factor*R1(i);
25.         while x1(i)+R1(i)>Lxmax|x1(i)-R1(i)<Lxmin|y1(i)+R1(i)>Lymax|y1(i)-R1(i)<Lymin
26.             R1(i)=Factor*R1(i);
27.         end
28.     end
29.     if R1(i)>Rmin
30.         rectangle('position',[x1(i)-R1(i),y1(i)-R1(i),2*R1(i),2*R1(i)],'curvature',[1,1],'facecolor','w')
31.         axis equal
32.         hold on
33.     else
34.         count=count+1;
35.         R1(i)=0;
36.     end
37. end
39. for i=2:n
40.     R1(i)=Rmin+rand*(Rmax-Rmin);
41.     RT(i)=R1(i);
42.     x1(i)=Lxmin+rand*(Lxmax-Lxmin);
43.     y1(i)=Lymin+rand*(Lymax-Lymin);
44.     if x1(i)+R1(i)<=Lxmax&x1(i)-R1(i)>=Lxmin&y1(i)+R1(i)<=Lymax&y1(i)-R1(i)>=Lymin
45.         Disp=sqrt([x1(i)-x1(1:i-1)].^2+[y1(i)-y1(1:i-1)].^2);
46.         SumR=R1(i)+R1(1:i-1);
47.         if all(Disp./SumR>1)
48.             if R1(i)>Rmin
49.                 rectangle('position',[x1(i)-R1(i),y1(i)-R1(i),2*R1(i),2*R1(i)],'curvature',[1,1],'facecolor','w')
50.                 axis equal
51.                 hold on
52.             else
53.                 count=count+1;
54.                 R1(i)=0;
55.             end
56.         else
57.             while any(Disp./SumR<=1)|x1(i)+R1(i)>Lxmax|x1(i)-R1<Lxmin|y1(i)+R1(i)>Lymax|y1(i)-R1(i)<Lymin
58.                 R1(i)=Factor*R1(i);
59.                 x1(i)=Lxmin+rand*(Lxmax-Lxmin);
60.                 y1(i)=Lymin+rand*(Lymax-Lymin);
61.                 Disp=sqrt([x1(i)-x1(1:i-1)].^2+[y1(i)-y1(1:i-1)].^2);
62.                 SumR=R1(i)+R1(1:i-1);
63.             end
64.             if R1(i)>Rmin
65.                 rectangle('position',[x1(i)-R1(i),y1(i)-R1(i),2*R1(i),2*R1(i)],'curvature',[1,1],'facecolor','w')
66.                 axis equal
67.                 hold on
68.             else
69.                 count=count+1;
70.                 R1(i)=0;
71.             end
72.         end
73.     else
74.         while x1(i)+R1(i)>Lxmax|x1(i)-R1<Lxmin|y1(i)+R1(i)>Lymax|y1(i)-R1(i)<Lymin
75.             R1(i)=Factor*R1(i);
76.         end
77.         Disp=sqrt([x1(i)-x1(1:i-1)].^2+[y1(i)-y1(1:i-1)].^2);
78.         SumR=R1(i)+R1(1:i-1);
79.         if all(Disp./SumR>1)
80.             if R1(i)>Rmin
81.                 rectangle('position',[x1(i)-R1(i),y1(i)-R1(i),2*R1(i),2*R1(i)],'curvature',[1,1],'facecolor','w')
82.                 axis equal
83.                 hold on
84.             else
85.                 count=count+1;
86.                 R1(i)=0;
87.             end
88.         else
89.             while any(Disp./SumR<=1)|x1+R1>Lxmax|x1-R1<Lxmin|y1+R1>Lymax|y1-R1<Lymin
90.                 R1(i)=Factor*R1(i);
91.                 x1(i)=Lxmin+rand*(Lxmax-Lxmin);
92.                 y1(i)=Lymin+rand*(Lymax-Lymin);
93.                 Disp=sqrt([x1(i)-x1(1:i-1)].^2+[y1(i)-y1(1:i-1)].^2);
94.                 SumR=R1(i)+R1(1:i-1);
95.             end
96.             if R1(i)>Rmin
97.                 rectangle('position',[x1(i)-R1(i),y1(i)-R1(i),2*R1(i),2*R1(i)],'curvature',[1,1],'facecolor','w')
98.                 axis equal
99.                 hold on
100.             else
101.                 count=count+1;
102.                 R1(i)=0;
103.             end
104.         end
105.     end
106. end
107. % axis([Lxmin,Lxmax,Lymin,Lymax])
108. x1(R1==0)=[];
109. y1(R1==0)=[];
110. R1(R1==0)=[];
111. R1,x1,y1
112. RT=R1;
113. count
114. if count>0
115.     while length(RT)~=n
116.         R2=Rmin+rand*(Rmax-Rmin);
117.         x2=Lxmin+rand*(Lxmax-Lxmin);
118.         y2=Lymin+rand*(Lymax-Lymin);
119.         Disp=sqrt([x2-x1].^2+[y2-y1].^2);
120.         SumR=R2+R1;
121.         if all(Disp./SumR>1)
122.             RT=[RT,R2];
123.             rectangle('position',[x2-R2,y2-R2,2*R2,2*R2],'curvature',[1,1],'facecolor','w')
124.             axis equal
125.             hold on
126.         end
127.     end
128. end
130. R1
131. set(gca,'color','k','fontsize',20,'fontweight','bold')
132. title('随机画圆程序','fontsize',20,'fontweight','bold','fontname','隶书')
133. Lxmin=Lxmin-Rmax;
134. Lxmax=Lxmax+Rmax;
135. Lymin=Lymin-Rmax;
136. Lymax=Lymax+Rmax;
137. axis([Lxmin,Lxmax,Lymin,Lymax])

复制代码

yudayinghualuo

请问一下，如果圆心随机分布，圆的半径在5-20之间，在圆的表面加一个0.5的薄层，该怎样编程啊

socraticliq

学习一下，如果规定了容积率该怎么做呢