[Case Study]利用MATLAB燃放烟花（礼花）

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

1. function firework
   
2. % 烟花烟花满天飞
   
3. % CopyRight：xiezhh(谢中华)
   
4. % 2011.6.25
   
5. 
   
6. OldHandle = findobj( 'Type', 'figure', 'Tag', 'FireWork' ) ;
   
7. if ishandle(OldHandle)
   
8.     close(OldHandle) ;
   
9. end
   
10. % 图形窗口初始化
    
11. fig = figure('units','normalized','position',[0.1 0.1 0.8 0.8],...
    
12.     'menubar','none','name','烟花满天飞（谢中华制作）',...
    
13.     'numbertitle','off','color',[0 0 0],'tag','FireWork');
    
14. % 烟花炸开前的初始位置
    
15. h0 = line(0,0,0,'marker','.',...
    
16.     'markersize',18,...
    
17.     'MarkerEdgeColor',[1 1 1],...
    
18.     'MarkerFaceColor',[1 1 1],...
    
19.     'EraseMode' , 'xor',...
    
20.     'color',[0 0 0]);
    
21. % 设置坐标系显示属性
    
22. axis equal
    
23. axis([-50 50 -50 50 0 100])
    
24. axis off
    
25. view(-42,22)
    
26. 
    
27. % 设置参数
    
28. rate = 1:-0.01:0;  % 颜色衰减率
    
29. v0 = 45;  % 烟花头初始速度
    
30. g = -9.8;   % 重力加速度
    
31. usedtime = -v0/g;  % 烟花头爆炸前所经历时间
    
32. zs = v0*usedtime+0.5*g*usedtime^2;  %烟花头爆炸前达到的最高高度
    
33. 
    
34. theta = rand(250,1)*2*pi;  % 各粒子速度的方位角
    
35. phi =  rand(250,1)*2*pi-pi;  % 各粒子速度的仰角
    
36. age = 20;  % 粒子生存期
    
37. % 常用颜色矩阵
    
38. colormat = [1.0  0.5  0.5
    
39.     1.0  0.75  0.5
    
40.     1.0  1.0  0.5
    
41.     0.75  1.0  0.5
    
42.     0.5  1.0  0.5
    
43.     0.5  1.0  0.75
    
44.     0.5  1.0  1.0
    
45.     0.5  0.75  1.0
    
46.     0.5  0.5  1.0
    
47.     0.75  0.5  1.0
    
48.     1.0  0.5  1.0
    
49.     1.0  0.5  0.75];
    
50. % 随机产生各粒子对应的颜色序号
    
51. colorid = randsample(12,250,true);
    
52. % 粒子对应的颜色矩阵
    
53. colormat = colormat(colorid,:);
    
54. % 粒子颜色与背景色（夜色）的距离
    
55. colordist = sqrt(sum(colormat.^2,2));
    
56. v1 = 20;  % 粒子的初始速度
    
57. k = 1;   % 颜色衰减率初始序号
    
58. 
    
59. timerA = timer('TimerFcn',@TimerFcnA,...
    
60.     'executionmode','fixedspacing','Period',0.001);
    
61. start(timerA);
    
62. h = getappdata(gcf,'HandleParticle');
    
63. timerB = timer('TimerFcn',{@TimerFcnB},...
    
64.     'executionmode','fixedspacing','period',0.001);
    
65. 
    
66. %--------------------------------------------------------------------------
    
67. %  定时器回调函数（烟花头）
    
68. %--------------------------------------------------------------------------
    
69. function TimerFcnA(timerA,event)
    
70. ta = get(timerA,'TasksExecuted')*0.1;
    
71. if ta <= usedtime
    
72.     z = v0*ta+0.5*g*ta^2;
    
73.     set(h0,'zdata',z,...
    
74.         'MarkerEdgeColor',[1 1 1],...
    
75.         'MarkerFaceColor',[1 1 1])
    
76.     drawnow
    
77.     %pause(0.01)
    
78. else
    
79.     delete(h0)
    
80.     stop(timerA);
    
81.     x0 = zeros(2,250);
    
82.     y0 = zeros(2,250);
    
83.     z0 = zs*ones(2,250);
    
84.     h = line(x0,y0,z0,'marker','h',...
    
85.         'markersize',12,...
    
86.         'MarkerEdgeColor',[1 1 1],...
    
87.         'MarkerFaceColor',[1 1 1],...
    
88.         'EraseMode' , 'xor',...
    
89.         'color',[0 0 0]);
    
90.     setappdata(gcf,'HandleParticle',h);
    
91.     start(timerB);
    
92. end
    
93. 
    
94. end
    
95. %--------------------------------------------------------------------------
    
96. %  定时器回调函数（粒子）
    
97. %--------------------------------------------------------------------------
    
98. function TimerFcnB(timerB,event)
    
99. tb = get(timerB,'TasksExecuted')*0.15;
    
100. if age>0 && any(colordist>=0.05)
     
101.     colormat = colormat*rate(k);
     
102.     colordist = sqrt(sum(colormat.^2,2));
     
103.     for i = 1:250
     
104.         xi = v1*cos(phi(i))*cos(theta(i))*tb;
     
105.         yi = v1*cos(phi(i))*sin(theta(i))*tb;
     
106.         zi = zs+v1*sin(phi(i))*tb+0.5*g*tb^2;
     
107.         set(h(i),'XData',xi,'YData',yi,'ZData',zi,...
     
108.             'MarkerEdgeColor',colormat(i,:),...
     
109.             'MarkerFaceColor',colormat(i,:))
     
110.     end
     
111.     drawnow
     
112.     age = age-0.1;
     
113.     k = k+1;
     
114. else
     
115.     stop(timerB);
     
116.     delete(h)
     
117. end
     
118. 
     
119. end
     
120. %-------------------
     
121. end
     
122. 复制代码

复制代码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：Case study MATLAB matla study atlab firework function position figure 坐标系