加载中…【原创】利用MATLAB燃放烟花（礼花）
(2011-06-26 16:24:51)
标签：
matlabskymatlab技术论坛matlab杂谈 |
分类： MATLAB |
【原创】利用MATLAB燃放烟花（礼花）
http://www.matlabsky.com/forum-viewthread-tid-18210-fromuid-1.html
相应的MATLAB程序如下，希望感兴趣的版友在此基础上进一步完善。
function firework
% 烟花烟花满天飞
% CopyRight：xiezhh(谢中华)
% 2011.6.25
OldHandle = findobj( 'Type', 'figure', 'Tag', 'FireWork' )
;
if ishandle(OldHandle)
end
% 图形窗口初始化
fig = figure('units','normalized','position',[0.1 0.1 0.8
0.8],...
% 烟花炸开前的初始位置
h0 = line(0,0,0,'marker','.',...
% 设置坐标系显示属性
axis equal
axis([-50 50 -50 50 0 100])
axis off
view(-42,22)
% 设置参数
rate = 1:-0.01:0; % 颜色衰减率
v0 = 45; % 烟花头初始速度
g = -9.8; % 重力加速度
usedtime = -v0/g; % 烟花头爆炸前所经历时间
zs = v0*usedtime+0.5*g*usedtime^2;
%烟花头爆炸前达到的最高高度
theta = rand(250,1)*2*pi; % 各粒子速度的方位角
phi = rand(250,1)*2*pi-pi;
% 各粒子速度的仰角
age = 20; % 粒子生存期
% 常用颜色矩阵
colormat = [1.0 0.5
0.5
% 随机产生各粒子对应的颜色序号
colorid = randsample(12,250,true);
% 粒子对应的颜色矩阵
colormat = colormat(colorid,:);
% 粒子颜色与背景色（夜色）的距离
colordist = sqrt(sum(colormat.^2,2));
v1 = 20; % 粒子的初始速度
k = 1; % 颜色衰减率初始序号
timerA = timer('TimerFcn',@TimerFcnA,...
start(timerA);
h = getappdata(gcf,'HandleParticle');
timerB = timer('TimerFcn',{@TimerFcnB},...
%--------------------------------------------------------------------------
% 定时器回调函数（烟花头）
%--------------------------------------------------------------------------
function TimerFcnA(timerA,event)
ta = get(timerA,'TasksExecuted')*0.1;
if ta <= usedtime
else
end
end
%--------------------------------------------------------------------------
% 定时器回调函数（粒子）
%--------------------------------------------------------------------------
function TimerFcnB(timerB,event)
tb = get(timerB,'TasksExecuted')*0.15;
if age>0 &&
any(colordist>=0.05)
else
end
end
%-------------------
end
喜欢
0
赠金笔