[原创]一朵玫瑰花的程序

qibbxxt

之前看过在一本书上看见有人用C语言画玫瑰花，于是我把它变成Matlab，把代码发在校内上，但是很少有人看，今天发在这里，程序不是很通用，因为当时急于看到花的样子

1. 
   
2. function plot_rose
   
3. draw_main(450,90);
   
4. function draw_main(x,y)
   
5. %粉红色玫瑰
   
6. arcdata{1}=[65 -60 150 350 8
   
7.     66 -54 300 470 8
   
8.     65 -56 30 230 10
   
9.     64 -57 300 490 17];
   
10. ellipsedata{1}=[73 -30 250 450 27 40
    
11.     59 -30 100 290 27 40
    
12.     65  -40 140 270 20 30];
    
13. arcdata{2}=[0 0 150 350 12
    
14.     1 8 280 470 12
    
15.     0 2 30 230 16
    
16.     0 3 80 240 28
    
17.     2 8 180 330 22
    
18.     -2 2 310 460 25];
    
19. ellipsedata{2}=[-12 30 120 300 30 40
    
20.     10 28 250 423 30 42
    
21.     -4  10 290 393 30 40];
    
22. ellipsedata{3}=[120 5 0 360 15 25];
    
23. ellipsedata{4}=[-70 10 0 360 14 20];
    
24. rose_e{2}=[x+16 y+32 235 355 26 35
    
25.     x-15 y+32 190 310 30 35
    
26.     x+0 y+35 190 350 43 50];
    
27. rose_e{1}=[x+80 y-48 220 350 22 50
    
28.     x+50 y-48 190 320 22 50
    
29.     x+65 y-28 180 360 36 50];
    
30. rose_e{3}=[x+120 y-6 200 340 17 25
    
31.     x+120 y+7 160 380 17 27];
    
32. rose_e{4}=[x-70 y+15 140 390 17 20
    
33.     x-75 y-10 205 340 10 30
    
34.     x-60 y-10 195 340 5 30];
    
35. arcdata{3}=[0 82 190 350 6];
    
36. ss={'m','r','b','y'};
    
37. s0={'k','k','g'};
    
38. figure('menubar','none','numbertitle','off','name','rose');
    
39. % figure('numbertitle','off','name','rose');
    
40. hold on
    
41. for j=1:length(ellipsedata)
    
42.     for i=1:size(ellipsedata{j},1)
    
43.         rectangle('Position',[x+ellipsedata{j}(i,1)-ellipsedata{j}(i,5),y+ellipsedata{j}(i,2)-ellipsedata{j}(i,6),2*ellipsedata{j}(i,5),2*ellipsedata{j}(i,6)],'Curvature',[1,1],...
    
44.             'FaceColor',ss{j})
    
45.     end
    
46.     if j<4
    
47.         for i=1:size(arcdata{j},1)
    
48.             draw_arc(x+arcdata{j}(i,1),y+arcdata{j}(i,2),arcdata{j}(i,3),arcdata{j}(i,4),arcdata{j}(i,5),s0{j});
    
49.         end
    
50.     end
    
51.     for i=1:size(rose_e{j},1)
    
52.         plot_rose_e(rose_e{j},j);
    
53.     end
    
54. end
    
55. for j=1:3
    
56.     zhuzhi=[x-98,y+100+j,255,371,100,80
    
57.         x-20,y+30+j,260,358,140,140
    
58.         x+224,y+20+j,180,218,160,140];
    
59.     draw_branch(zhuzhi);
    
60. end
    
61. ce_branch=[x+70,y+34,180,233,140,140;
    
62.   x,y+40,205,255,100,120;
    
63.   x+135,y-30,209,249,72,120;
    
64.   x,y+20,263,301,100,120;
    
65.   x+85,y-10,278,305,100,120;
    
66.   x+100,y-62,282,308,90,120;
    
67.   x-50,y-10,277,314,30,120;
    
68.   x+70,y+80,222,266,52,120;
    
69.   x-60,y-45,229,266,52,120;
    
70.   x+79,y-45,229,266,52,120;
    
71.   x+84,y,224,273,52,120;
    
72.   x+110,y+40,240,282,100,120];
    
73. draw_branch(ce_branch);
    
74. t_leaf=[x+168,y+282,10,20,-40;
    
75.   x+160,y+252,8,16,260;
    
76.   x+145,y+270,8,16,-15;
    
77.   x+156,y+224,10,20,-45;
    
78.   x+150,y+200,8,16,270;
    
79.   x+135,y+220,8,16,-10;
    
80.   x+146,y+144,8,16,-80;
    
81.   x+130,y+130,6,12,235;
    
82.   x+125,y+154,7,14,-10;
    
83.   x+78,y+98,6,12,-90;
    
84.   x+60,y+90,5,10,180;
    
85.   x+70,y+109,5,10,-45;
    
86.   x-125,y+270,12,24,60;
    
87.   x-95,y+270,10,20,10;
    
88.   x-110,y+245,10,20,90;
    
89.   x-105,y+220,10,20,45;
    
90.   x-100,y+190,8,16,135;
    
91.   x-75,y+210,8,16,-45;
    
92.   x+65,y+190,10,20,-45;
    
93.   x+40,y+185,8,16,0;
    
94.   x+55,y+165,8,16,90];
    
95. for j=1:size(t_leaf,1)
    
96.     draw_elli(t_leaf(j,1),t_leaf(j,2),t_leaf(j,3),t_leaf(j,4),t_leaf(j,5));
    
97. end
    
98. axis ij
    
99. axis off
    
100. set(gcf,'color','k');
     
101. daspect([1,1,1])
     
102. 
     
103. % 画旋转椭圆
     
104. function draw_elli(x0,y0,a,b,theta)
     
105. theta=-theta;
     
106. t = 0:0.01:2*pi;
     
107. x =  a*cos(t);
     
108. y =  b*sin(t);
     
109. xy = zeros(2,length(t));
     
110. for i = 1:length(t)
     
111.     xy(:,i) = expm([0,-pi/180*theta;pi/180*theta,0])*[x(i);y(i)];
     
112. end
     
113. plot(x0 +xy(1,:),480-y0 + xy(2,:),'g');
     
114. % 画弧线
     
115. function draw_arc(x0,y0,st_a,en_a,r,c_c)
     
116. t = (pi/180*st_a:0.01:pi/180*en_a)-pi;
     
117. x = x0 + r*cos(t);
     
118. y = y0 + r*sin(t);
     
119. plot(x,y,c_c);
     
120. % 花萼
     
121. function plot_rose_e(matr,ch)
     
122. x=[];y=[];n=size(matr,1);
     
123. for i=1:size(matr,1)
     
124.     t=(matr(i,3):matr(i,4))*pi/180-pi;
     
125.     x1=matr(i,1)+matr(i,5)*cos(t);
     
126.     y1=matr(i,2)+matr(i,6)*sin(t);
     
127.     if i<3
     
128.         x1=2*matr(i,1)-x1;
     
129.     end
     
130.     x1=fliplr(x1);
     
131.     y1=fliplr(y1);
     
132.     if n<3 && i<2
     
133.         x1=fliplr(x1);
     
134.         y1=fliplr(y1);
     
135.     end
     
136.     if ch==4 && i==2
     
137.         x1=fliplr(x1);
     
138.         y1=fliplr(y1);
     
139.     end
     
140.    
     
141.     x=[x,x1];
     
142.     y=[y,y1];
     
143. end
     
144. fill(x,y,'g');
     
145. % 画树枝
     
146. function draw_branch(zhuzhi)
     
147. for k=1:size(zhuzhi,1)
     
148.     t=(zhuzhi(k,3):zhuzhi(k,4))*pi/180-pi;
     
149.     x1=zhuzhi(k,1)+zhuzhi(k,5)*cos(t);
     
150.     y1=zhuzhi(k,2)+zhuzhi(k,6)*sin(t);
     
151.     x1=2*zhuzhi(k,1)-x1;
     
152.     plot(x1,y1,'g','linewidth',2);
     
153. end
     

复制代码

qinwei

挺好看，不错！

bainhome

这样画画可太累了，在表达了qibbxxt老弟的MATLAB水准、持之以恒的耐心之外，也表明MATLAB实在不是画画最合适的软件。

wyj222129

不错，顶一个！！

chenyy2010

呵呵，很漂亮！谢谢版主分享！

smallpiano

很漂亮的花，欣赏了。

ajun2100

不错，初学者，学习一下

zhangxuwhu

人才:victory:

ljw880107

这个程序不在说明MATLAB能画画，只是说明编程者是高手。

wpanysys

太厉害了啊！:)

xuxu1457

相当强大，很漂亮 多谢分享

hftcl2006

我复制了你的代码，怎么画不出来呢？

wangzhen051099

我试了，可以的！楼主很强大！

allenfieldin

我靠，太牛逼了，回去画给我女盆友看，嘿嘿，学习一下