利用Matlab制作一個(gè)賊簡(jiǎn)單的粒子圣誕樹
上一期圣誕樹有些人說代碼太復(fù)雜,這期來個(gè)賊簡(jiǎn)單的水晶球里的圣誕樹:
這期就是非常簡(jiǎn)單的三個(gè)圖形組合:
圣誕樹主體
% 圣誕樹形狀生成函數(shù) h0=25;r0=12;a0=.5.*pi; XFunc=@(h,r,a,z) (h-z)./h.*r.*cos(a.*z); YFunc=@(h,r,a,z) (h-z)./h.*r.*sin(a.*z); ZFunc=@(h,r,z) z+sqrt(12^2-((h-z)./h.*r).^2); % 生成最外圈散點(diǎn) Z0=linspace(0,12.5,150); X0=XFunc(h0,r0,a0,Z0); Y0=YFunc(h0,r0,a0,Z0); % 生成中心點(diǎn) ZC=ZFunc(h0,r0,Z0); t0=linspace(.05,1,25)'; % 生成樹枝上的散點(diǎn) treeXYZ=zeros([length(Z0)*length(t0),3]); for i=1:length(Z0) treeXYZ((i-1)*length(t0)+1:i*length(t0),:)=... [0,0,ZC(i)]+[X0(i),Y0(i),Z0(i)-ZC(i)].*t0; end % 繪制圣誕樹 treeXYZ=treeXYZ+rand(size(treeXYZ))./3; treeCData=repmat([0,235,81]./255,[length(Z0)*length(t0),1]); treeCData=treeCData+rand(size(treeXYZ))./2; treeCData(treeCData>1)=1; scatter3(treeXYZ(:,1),treeXYZ(:,2),treeXYZ(:,3)-3,3,'o','filled','CData',treeCData,... 'MarkerFaceAlpha',.7,'MarkerEdgeColor','none')
主體就是錐形螺旋線和中軸之間連線上取的一些散點(diǎn):
背景修飾
ax=gca; set(gcf,'Color',[22,32,51]./255); % 坐標(biāo)區(qū)域修飾 set(ax,'XLim',[-27 27],'YLim',[-27,27],'ZLim',[-9,28],'PlotBoxAspectRatio',[1,1,.8],... 'XColor','none','YColor','none','ZColor','none','Color',[22,32,51]./255,'View',[-37.5,19.5],... 'Position',[-.5,-.2,2,1.4],'Tag',char([100,105,115,112,40,39,20316,32773,... 58,115,108,97,110,100,97,114,101,114,39,41]));eval(ax.Tag)
繪制地面
這里將距離中心點(diǎn)的距離映射為余弦函數(shù):
% 繪制地面 [surfX,surfY]=meshgrid(linspace(-25,25,60)); surfZ=cos(hypot(surfX,surfY))-5; surfZ=surfZ+rand(size(surfX)).*1.2; surfX=surfX+rand(size(surfX))./1.8; surfY=surfY+rand(size(surfX))./1.8; scatter3(surfX,surfY,surfZ,3,'o','filled','CData',[91,131,253]./255,... 'MarkerFaceAlpha',.4,'MarkerEdgeColor','none')
水晶球
就是一個(gè)半透明的球:
% 繪制水晶球 [ballX,ballY,ballZ]=sphere(20); surf(ballX.*18,ballY.*18,ballZ.*15+8,'EdgeColor','none','FaceAlpha',.1); colormap(gray); light;lighting phong;shading interp;
添加雪花并旋轉(zhuǎn)視角
% 繪制雪花 snowXYZ1=rand(27,3); snowXYZ1(:,1:2)=snowXYZ1(:,1:2).*50-25; snowXYZ1(:,3)=snowXYZ1(:,3).*30; snowXYZ2=rand(60,3); snowXYZ2(:,1:2)=snowXYZ2(:,1:2).*50-25; snowXYZ2(:,3)=snowXYZ2(:,3).*30; snowHdl1=plot3(snowXYZ1(:,1),snowXYZ1(:,2),snowXYZ1(:,3),'*','Color',[1 1 1]); snowHdl2=plot3(snowXYZ2(:,1),snowXYZ2(:,2),snowXYZ2(:,3),'.','Color',[.6,.6,.6]); % 旋轉(zhuǎn)圖像、雪花飄落 for i=1:1e8 snowXYZ1(:,3)=snowXYZ1(:,3)-.1;snowXYZ2(:,3)=snowXYZ2(:,3)-.01; snowXYZ1(snowXYZ1(:,3)<-5,3)=30;snowXYZ2(snowXYZ2(:,3)<-5,3)=30; snowHdl1.XData=snowXYZ1(:,1);snowHdl1.YData=snowXYZ1(:,2);snowHdl1.ZData=snowXYZ1(:,3); snowHdl2.XData=snowXYZ2(:,1);snowHdl2.YData=snowXYZ2(:,2);snowHdl2.ZData=snowXYZ2(:,3); view([i,19.5]); drawnow;pause(.05) end
完整代碼
function XmasTree2022_2 % @author:slandarer % 圣誕樹形狀生成函數(shù) h0=25;r0=12;a0=.5.*pi; XFunc=@(h,r,a,z) (h-z)./h.*r.*cos(a.*z); YFunc=@(h,r,a,z) (h-z)./h.*r.*sin(a.*z); ZFunc=@(h,r,z) z+sqrt(12^2-((h-z)./h.*r).^2); % 生成最外圈散點(diǎn) Z0=linspace(0,12.5,150); X0=XFunc(h0,r0,a0,Z0); Y0=YFunc(h0,r0,a0,Z0); % 生成中心點(diǎn) ZC=ZFunc(h0,r0,Z0); t0=linspace(.05,1,25)'; % 生成樹枝上的散點(diǎn) treeXYZ=zeros([length(Z0)*length(t0),3]); for i=1:length(Z0) treeXYZ((i-1)*length(t0)+1:i*length(t0),:)=... [0,0,ZC(i)]+[X0(i),Y0(i),Z0(i)-ZC(i)].*t0; end ax=gca;hold on; set(gcf,'Color',[22,32,51]./255); % 繪制圣誕樹 treeXYZ=treeXYZ+rand(size(treeXYZ))./3; treeCData=repmat([0,235,81]./255,[length(Z0)*length(t0),1]); treeCData=treeCData+rand(size(treeXYZ))./2; treeCData(treeCData>1)=1; scatter3(treeXYZ(:,1),treeXYZ(:,2),treeXYZ(:,3)-3,3,'o','filled','CData',treeCData,... 'MarkerFaceAlpha',.7,'MarkerEdgeColor','none') % 繪制水晶球 [ballX,ballY,ballZ]=sphere(20); surf(ballX.*18,ballY.*18,ballZ.*15+8,'EdgeColor','none','FaceAlpha',.1); colormap(gray); light;lighting phong;shading interp; % 繪制地面 [surfX,surfY]=meshgrid(linspace(-25,25,60)); surfZ=cos(hypot(surfX,surfY))-5; surfZ=surfZ+rand(size(surfX)).*1.2; surfX=surfX+rand(size(surfX))./1.8; surfY=surfY+rand(size(surfX))./1.8; scatter3(surfX,surfY,surfZ,3,'o','filled','CData',[91,131,253]./255,... 'MarkerFaceAlpha',.4,'MarkerEdgeColor','none') % 坐標(biāo)區(qū)域修飾 set(ax,'XLim',[-27 27],'YLim',[-27,27],'ZLim',[-9,28],'PlotBoxAspectRatio',[1,1,.8],... 'XColor','none','YColor','none','ZColor','none','Color',[22,32,51]./255,'View',[-37.5,19.5],... 'Position',[-.5,-.2,2,1.4],'Tag',char([100,105,115,112,40,39,20316,32773,... 58,115,108,97,110,100,97,114,101,114,39,41]));eval(ax.Tag) % 繪制雪花 snowXYZ1=rand(27,3); snowXYZ1(:,1:2)=snowXYZ1(:,1:2).*50-25; snowXYZ1(:,3)=snowXYZ1(:,3).*30; snowXYZ2=rand(60,3); snowXYZ2(:,1:2)=snowXYZ2(:,1:2).*50-25; snowXYZ2(:,3)=snowXYZ2(:,3).*30; snowHdl1=plot3(snowXYZ1(:,1),snowXYZ1(:,2),snowXYZ1(:,3),'*','Color',[1 1 1]); snowHdl2=plot3(snowXYZ2(:,1),snowXYZ2(:,2),snowXYZ2(:,3),'.','Color',[.6,.6,.6]); % 旋轉(zhuǎn)圖像、雪花飄落 for i=1:1e8 snowXYZ1(:,3)=snowXYZ1(:,3)-.1;snowXYZ2(:,3)=snowXYZ2(:,3)-.01; snowXYZ1(snowXYZ1(:,3)<-5,3)=30;snowXYZ2(snowXYZ2(:,3)<-5,3)=30; snowHdl1.XData=snowXYZ1(:,1);snowHdl1.YData=snowXYZ1(:,2);snowHdl1.ZData=snowXYZ1(:,3); snowHdl2.XData=snowXYZ2(:,1);snowHdl2.YData=snowXYZ2(:,2);snowHdl2.ZData=snowXYZ2(:,3); view([i,19.5]); drawnow;pause(.05) end end
以上就是利用Matlab制作一個(gè)賊簡(jiǎn)單的粒子圣誕樹的詳細(xì)內(nèi)容,更多關(guān)于Matlab粒子圣誕樹的資料請(qǐng)關(guān)注腳本之家其它相關(guān)文章!
相關(guān)文章
C++實(shí)現(xiàn)判斷字符串是否回文實(shí)例解析
這篇文章主要介紹了C++實(shí)現(xiàn)判斷字符串是否回文,其中采用了數(shù)據(jù)結(jié)構(gòu)中棧以及過濾字符等技術(shù),,需要的朋友可以參考下2014-07-07C++編程中隊(duì)內(nèi)聯(lián)函數(shù)的理解和使用
這篇文章主要介紹了C++編程中隊(duì)內(nèi)聯(lián)函數(shù)的理解和使用,簡(jiǎn)單舉例講解了inline關(guān)鍵字引出的內(nèi)聯(lián)函數(shù)的相關(guān)知識(shí),需要的朋友可以參考下2016-01-01C++ Coroutine簡(jiǎn)單學(xué)習(xí)教程
這篇文章主要為大家詳細(xì)介紹了C++ Coroutine的簡(jiǎn)單學(xué)習(xí)教程,具有一定的參考價(jià)值,感興趣的小伙伴們可以參考一下2017-08-08C++項(xiàng)目求Fibonacci數(shù)列的參考解答
今天小編就為大家分享一篇關(guān)于C++項(xiàng)目求Fibonacci數(shù)列的參考解答,小編覺得內(nèi)容挺不錯(cuò)的,現(xiàn)在分享給大家,具有很好的參考價(jià)值,需要的朋友一起跟隨小編來看看吧2019-02-02解析C語言中結(jié)構(gòu)體struct的對(duì)齊問題
這篇文章主要介紹了C語言中結(jié)構(gòu)體struct的對(duì)齊問題,作者深入到內(nèi)存分配方面來進(jìn)行解析,需要的朋友可以參考下2016-04-04C++利用std::forward_list查找插入數(shù)據(jù)方法示例
這篇文章主要給大家介紹了關(guān)于C++利用std::forward_list查找插入數(shù)據(jù)的相關(guān)資料,文中先對(duì)std::forward_list進(jìn)行了詳細(xì)的介紹,而后通過示例代碼給大家介紹了查找的方法,需要的朋友可以參考借鑒,下面話不多說了,來一起看看吧。2017-08-08