脚本之家 服务器常用软件
快捷导航
您的位置：首页 → 软件编程 → C 语言 → Matlab粒子圣诞树

利用Matlab制作一个贼简单的粒子圣诞树

 更新时间：2022年12月18日 10:25:22   作者：slandarer  
圣诞节快到了,本文用Matlab绘制了圣诞树祝你们圣诞节快乐,所以下面这篇文章主要给大家介绍了关于如何利用Matlab制作一个贼简单的粒子圣诞树,需要的朋友可以参考下

上一期圣诞树有些人说代码太复杂，这期来个贼简单的水晶球里的圣诞树：

这期就是非常简单的三个图形组合：

圣诞树主体

% 圣诞树形状生成函数
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);
% 生成最外圈散点
Z0=linspace(0,12.5,150);
X0=XFunc(h0,r0,a0,Z0);
Y0=YFunc(h0,r0,a0,Z0);
% 生成中心点
ZC=ZFunc(h0,r0,Z0);
t0=linspace(.05,1,25)';
% 生成树枝上的散点
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')

主体就是锥形螺旋线和中轴之间连线上取的一些散点：

背景修饰

ax=gca;
set(gcf,'Color',[22,32,51]./255);
% 坐标区域修饰
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)

绘制地面

这里将距离中心点的距离映射为余弦函数：

% 绘制地面
[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')

水晶球

就是一个半透明的球:

% 绘制水晶球
[ballX,ballY,ballZ]=sphere(20);
surf(ballX.*18,ballY.*18,ballZ.*15+8,'EdgeColor','none','FaceAlpha',.1);
colormap(gray);
light;lighting phong;shading interp;

添加雪花并旋转视角

% 绘制雪花
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]);
% 旋转图像、雪花飘落
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

% 圣诞树形状生成函数
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);
% 生成最外圈散点
Z0=linspace(0,12.5,150);
X0=XFunc(h0,r0,a0,Z0);
Y0=YFunc(h0,r0,a0,Z0);
% 生成中心点
ZC=ZFunc(h0,r0,Z0);
t0=linspace(.05,1,25)';
% 生成树枝上的散点
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')
% 坐标区域修饰
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]);
% 旋转图像、雪花飘落
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制作一个贼简单的粒子圣诞树的详细内容，更多关于Matlab粒子圣诞树的资料请关注脚本之家其它相关文章！

您可能感兴趣的文章:

相关文章

  • C语言调试手段:锁定错误的实现方法

    C语言调试手段:锁定错误的实现方法

    本篇文章是对在C语言调试中，锁定错误的方法进行了详细的分析介绍，需要的朋友参考下
    2013-05-05
  • VisualStudio2022打包项目文件为.exe安装包

    VisualStudio2022打包项目文件为.exe安装包

    本文主要介绍了VisualStudio2022打包项目文件为.exe安装包，文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值，需要的朋友们下面随着小编来一起学习学习吧
    2023-07-07
  • C++数位DP复杂度统计数字问题示例详解

    C++数位DP复杂度统计数字问题示例详解

    这篇文章主要为大家介绍了利用C++数位DP的复杂度来统计数字问题的示例实现过程详解，有需要的朋友可以借鉴参考下，希望能够有所帮助，祝大家多多进步，早日升值加薪
    2021-11-11
  • C++深入讲解类与对象之OOP面向对象编程与封装

    C++深入讲解类与对象之OOP面向对象编程与封装

    学习过C语言的小伙伴知道：C语言是面向过程的，关注的是过程，分析出求解问题的步骤，通过函数调用逐步解决问题，接下来让我们详细的了解
    2022-05-05
  • OpenCV利用对比度亮度变换实现水印去除

    OpenCV利用对比度亮度变换实现水印去除

    OpenCV中去除水印最常用的方法是inpaint，通过图像修复的方法来去除水印。本文将介绍另一种方法：利用对比度亮度变换去除水印，需要的朋友可以参考一下
    2021-11-11
  • C++中String类的常用接口函数总结

    C++中String类的常用接口函数总结

    这篇文章主要介绍了C++中Stirng类的常用接口函数，文中有详细的代码示例供大家参考，对我们学习C++有一定的帮助，感兴趣的同学可以跟着小编一起来学习
    2023-06-06
  • C++ Boost weak_ptr智能指针超详细讲解

    C++ Boost weak_ptr智能指针超详细讲解

    智能指针是一种像指针的C++对象，但它能够在对象不使用的时候自己销毁掉。虽然STL提供了auto_ptr，但是由于不能同容器一起使用（不支持拷贝和赋值操作），因此很少有人使用。它是Boost各组件中，应用最为广泛的一个
    2022-11-11
  • C语言关于二叉树中堆的创建和使用整理

    C语言关于二叉树中堆的创建和使用整理

    大家好,这里是针对二叉树中堆结构的顺序储存，整理出来一篇博客供我们一起复习和学习,如果文章中有理解不当的地方,还希望朋友们在评论区指出,我们相互学习,共同进步
    2022-08-08
  • 如何基于C++解决RTSP取流报错问题

    如何基于C++解决RTSP取流报错问题

    这篇文章主要介绍了如何基于C++解决RTSP取流报错问题,文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值,需要的朋友可以参考下
    2020-08-08
  • C++ 实现汉诺塔的实例详解

    C++ 实现汉诺塔的实例详解

    这篇文章主要介绍了C++ 实现汉诺塔的实例详解的相关资料,这里主要说明C++中数据结构的递归的应用，需要的朋友可以参考下
    2017-08-08

最新评论