利用Matlab制作一個(gè)賊簡(jiǎn)單的粒子圣誕樹(shù)
更新時(shí)間：2022年12月18日 10:25:22 作者：slandarer
圣誕節(jié)快到了,本文用Matlab繪制了圣誕樹(shù)祝你們圣誕節(jié)快樂(lè),所以下面這篇文章主要給大家介紹了關(guān)于如何利用Matlab制作一個(gè)賊簡(jiǎn)單的粒子圣誕樹(shù),需要的朋友可以參考下
上一期圣誕樹(shù)有些人說(shuō)代碼太復(fù)雜，這期來(lái)個(gè)賊簡(jiǎn)單的水晶球里的圣誕樹(shù)：
這期就是非常簡(jiǎn)單的三個(gè)圖形組合：
圣誕樹(shù)主體
% 圣誕樹(shù)形狀生成函數(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)';
% 生成樹(shù)枝上的散點(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
% 繪制圣誕樹(shù)
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ù)形狀生成函數(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)';
% 生成樹(shù)枝上的散點(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);
% 繪制圣誕樹(shù)
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)單的粒子圣誕樹(shù)的詳細(xì)內(nèi)容，更多關(guān)于Matlab粒子圣誕樹(shù)的資料請(qǐng)關(guān)注腳本之家其它相關(guān)文章！
相關(guān)文章
C語(yǔ)言調(diào)試手段:鎖定錯(cuò)誤的實(shí)現(xiàn)方法
本篇文章是對(duì)在C語(yǔ)言調(diào)試中，鎖定錯(cuò)誤的方法進(jìn)行了詳細(xì)的分析介紹，需要的朋友參考下2013-05-05
VisualStudio2022打包項(xiàng)目文件為.exe安裝包
本文主要介紹了VisualStudio2022打包項(xiàng)目文件為.exe安裝包，文中通過(guò)示例代碼介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值，需要的朋友們下面隨著小編來(lái)一起學(xué)習(xí)學(xué)習(xí)吧2023-07-07
C++數(shù)位DP復(fù)雜度統(tǒng)計(jì)數(shù)字問(wèn)題示例詳解
這篇文章主要為大家介紹了利用C++數(shù)位DP的復(fù)雜度來(lái)統(tǒng)計(jì)數(shù)字問(wèn)題的示例實(shí)現(xiàn)過(guò)程詳解，有需要的朋友可以借鑒參考下，希望能夠有所幫助，祝大家多多進(jìn)步，早日升值加薪2021-11-11
C++深入講解類(lèi)與對(duì)象之OOP面向?qū)ο缶幊膛c封裝
學(xué)習(xí)過(guò)C語(yǔ)言的小伙伴知道：C語(yǔ)言是面向過(guò)程的，關(guān)注的是過(guò)程，分析出求解問(wèn)題的步驟，通過(guò)函數(shù)調(diào)用逐步解決問(wèn)題，接下來(lái)讓我們?cè)敿?xì)的了解2022-05-05
OpenCV利用對(duì)比度亮度變換實(shí)現(xiàn)水印去除
OpenCV中去除水印最常用的方法是inpaint，通過(guò)圖像修復(fù)的方法來(lái)去除水印。本文將介紹另一種方法：利用對(duì)比度亮度變換去除水印，需要的朋友可以參考一下2021-11-11
C++中String類(lèi)的常用接口函數(shù)總結(jié)
這篇文章主要介紹了C++中Stirng類(lèi)的常用接口函數(shù)，文中有詳細(xì)的代碼示例供大家參考，對(duì)我們學(xué)習(xí)C++有一定的幫助，感興趣的同學(xué)可以跟著小編一起來(lái)學(xué)習(xí)2023-06-06
C++?Boost?weak_ptr智能指針超詳細(xì)講解
智能指針是一種像指針的C++對(duì)象，但它能夠在對(duì)象不使用的時(shí)候自己銷(xiāo)毀掉。雖然STL提供了auto_ptr，但是由于不能同容器一起使用（不支持拷貝和賦值操作），因此很少有人使用。它是Boost各組件中，應(yīng)用最為廣泛的一個(gè)2022-11-11
C語(yǔ)言關(guān)于二叉樹(shù)中堆的創(chuàng)建和使用整理
大家好,這里是針對(duì)二叉樹(shù)中堆結(jié)構(gòu)的順序儲(chǔ)存，整理出來(lái)一篇博客供我們一起復(fù)習(xí)和學(xué)習(xí),如果文章中有理解不當(dāng)?shù)牡胤?還希望朋友們?cè)谠u(píng)論區(qū)指出,我們相互學(xué)習(xí),共同進(jìn)步2022-08-08
如何基于C++解決RTSP取流報(bào)錯(cuò)問(wèn)題
這篇文章主要介紹了如何基于C++解決RTSP取流報(bào)錯(cuò)問(wèn)題,文中通過(guò)示例代碼介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友可以參考下2020-08-08
C++ 實(shí)現(xiàn)漢諾塔的實(shí)例詳解
這篇文章主要介紹了C++ 實(shí)現(xiàn)漢諾塔的實(shí)例詳解的相關(guān)資料,這里主要說(shuō)明C++中數(shù)據(jù)結(jié)構(gòu)的遞歸的應(yīng)用，需要的朋友可以參考下2017-08-08