MATLAB程序演示煙花燃放過程

MATLAB程序演示煙花燃放過程煙花（禮花）的朵朵綻放給夜空帶來炫彩，瞬間絢麗至極，迸射出璀璨奪目的光彩。只是還來不及在腦海中印上花魂，她已曇花一現般，煙花（禮花）朵朵綻放，瞬間絢麗至極，迸射出璀璨奪目的光彩。只是還來不及在腦海中印上花魂，她已曇花一現般的消逝了。多少人癡迷于煙花的美，癡迷于她飄忽不定的幻影。利用MATLAB程序演示煙花的美，雖然她不如現實那樣逼真，但一樣會吸引人。附源程序：functionfirework%煙花煙花滿天飛%CopyRight：xiezhh(謝中華)%2011.6.25OldHandle=findobj('Type','figure','Tag','FireWork');ifishandle(OldHandle)close(OldHandle);end%圖形窗口初始化fig=figure('units','normalized','position',[0.10.10.80.8],...'menubar','none','name','煙花滿天飛（謝中華制作）',...'numbertitle','off','color',[000],'tag','FireWork');%煙花炸開前的初始位置h0=line(0,0,0,'linestyle','none',...'marker','.',...'markersize',18,...'MarkerEdgeColor',[111],...'MarkerFaceColor',[111],...'EraseMode','xor');%設置坐標系顯示屬性axisequalaxis([-5050-50500100])axisoffview(-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.00.50.51.00.750.51.01.00.50.751.00.50.51.00.50.51.00.750.51.01.00.50.751.00.50.51.00.750.51.01.00.51.01.00.50.75];%隨機產生各粒子對應的顏色序號colorid=randsample(12,250,true);%粒子對應的顏色矩陣colormat=colormat(colorid,:);%粒子顏色與背景色（夜色）的距離colordist=sqrt(sum(colormat.^2,2));v1=20;%粒子的初始速度k=1;%顏色衰減率初始序號timerA=timer('TimerFcn',@TimerFcnA,...'executionmode','fixedspacing','Period',0.001);start(timerA);h=getappdata(gcf,'HandleParticle');timerB=timer('TimerFcn',{@TimerFcnB},...'executionmode','fixedspacing','period',0.001);%--------------------------------------------------------------------------%定時器回調函數（煙花頭）%--------------------------------------------------------------------------functionTimerFcnA(timerA,event)ta=get(timerA,'TasksExecuted')*0.1;ifta<=usedtimez=v0*ta+0.5*g*ta^2;set(h0,'zdata',z,...'MarkerEdgeColor',[111],...'MarkerFaceColor',[111])drawnow%pause(0.01)elsedelete(h0)stop(timerA);x0=zeros(2,250);y0=zeros(2,250);z0=zs*ones(2,250);h=line(x0,y0,z0,'linestyle','none',...'marker','h',...'markersize',12,...'MarkerEdgeColor',[111],...'MarkerFaceColor',[111],...'EraseMode','xor');setappdata(gcf,'HandleParticle',h);start(timerB);endend%--------------------------------------------------------------------------%定時器回調函數（粒子）%--------------------------------------------------------------------------functionTimerFcnB(timerB,event)tb=get(timerB,'TasksExecuted')*0.15;ifage>0&&any(colordist>=0.05)colormat=colormat*rate(k);colordist=sqrt(sum(colormat.^2,2));fori=1:250xi=v1*cos(phi(i))*cos(theta(i))*tb;yi=v1*cos(phi(i))*sin(theta(i))*tb;zi=zs+v1*sin(phi(i))*tb+0.5*g*tb^2;set(h(i),'XData',xi,'YData',yi,'ZData',zi,...'MarkerEdgeColor',colormat(i,: