已知半径和圆心画圆

tz6091

几个随机的圆心.和半径.怎么画出圆呢?谢谢

branch

polar(2*pi,r)

bainhome

原先写过一个限制位置和圆大小的随机画圆程序：

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
    
38. 
    
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
     
129. 
     
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])
     

复制代码