当前位置： 首页 > news >正文

MATLAB | 世界杯来用MATLAB画个足球玩叭~

news 来源：原创 2024/4/26 23:22:11

世界杯教你用MATLAB画个超逼真的足球，

需要准备Partial Differential Equation Toolbox工具箱，同时因为用到了polyshape类所以至少需要R2017b版本。

绘制讲解

数据来源及说明

我是真的不想写注释了太麻烦了，给大家讲一下我的思路希望能够看懂，首先足球的数据点是通过:

[B,XYZ]=bucky;

导入的，但是导入的只有边链接信息，并没有给出哪几个点构成正五边形哪几个边构成正六边形。

spy(B)

通过spy展示一下稀疏矩阵发现，似乎每5行数据代表一个正五边形，但是正六边形的边的编号还是无法获得：

展示一下部分连接情况：

[B,XYZ]=bucky;
spy(B) 

G = graph(B);
A = adjacency(G);
H = graph(A(1:30,1:30));
h = plot(H);
axis equal

我们发现很多正六边形顶点的编号都完全不沾边，打算硬算。

硬算顶点连接情况

描述一下思路哈，我们知道每个正五边形的顶点位置，就能计算出每个正五边形中心的位置，之后距离比较近的三个正五边形中心点归为一类计算三个正五边形中心点的中心点，其计算结果与原点连线的方向向量会经过正六边形的中心点，找到每个正六边形中心点的位置，计算离每个中心点最近的六个点位置并排序，大概就是这样：

  1. 计算五边形顶点
  2. 计算五边形中心
  3. 最近三个五边形归为一组
  4. 计算每组中心
  5. 找到离六边形中心最近的六个点
  6. 为六边形点排序、捋顺

以上过程通过这段代码完成：

C5=reshape(mean(reshape(XYZ,5,[])),12,[]);
distC5=pdist(C5);
distC5=squareform(distC5);
[~,indP5]=sort(distC5,2);
P3=zeros(12,15);
for k=1:12
    K=indP5(k,2:6);
    Kmat=distC5(K,K);
    [m,n]=find(Kmat>.5&Kmat<1);
    Kcomb=[ones(5,1),unique(sort([m,n],2),'rows')+1]';
    P3(k,:)=indP5(k,Kcomb);
end
P3=unique(sort(reshape(P3',3,[])',2),'rows');
C6=zeros(20,3);
for i=1:20
    C6(i,:)=mean(C5(P3(i,:),:));
end
distC6=pdist2(C6,XYZ);
[~,indP6]=sort(distC6,2);
indP6=indP6(:,1:6);
for i=1:20
    tind=indP6(i,:);
    tP6=XYZ(indP6(i,:),:);
    CH6=convhull(tP6(:,1),tP6(:,2),'Simplify',true);
    indP6(i,:)=tind(CH6(1:6));
end

此时我们已经能画出一个有棱有角的足球了：

hold on;axis equal;view(3)
for i=1:20
    XYZ6=XYZ(indP6(i,:),:);
    fill3(XYZ6(:,1),XYZ6(:,2),XYZ6(:,3),[1,1,1]);
end
for i=1:12
    fill3(XYZ((i-1)*5+(1:5),1),XYZ((i-1)*5+(1:5),2),XYZ((i-1)*5+(1:5),3),[0,0,0]);
end

三角剖分

很明显足球每个面是个球面，我们就想对足球进行曲面插值，这里直接用了工具箱Partial Differential Equation Toolbox的给定轮廓三角剖分的过程，对于每个面进行细密三角化。

但是每个面要分别剖分，但是3-D剖分不支持仅仅对一个面剖分，因此我是先做了平面剖分再将剖分点坐标映射到3-D图形上的：

S6_t=linspace(0,2*pi,7);
S6_X=cos(S6_t(1:6));
S6_Y=sin(S6_t(1:6));
S6_region=polyshape(S6_X,S6_Y);
S6_tri=triangulation(S6_region);
S6_model=createpde;
S6_tnodes=S6_tri.Points';
S6_telements=S6_tri.ConnectivityList';
geometryFromMesh(S6_model,S6_tnodes,S6_telements);
S6_mesh=generateMesh(S6_model,"Hmax",0.1,"GeometricOrder","linear");
pdeplot(S6_model)


S5_t=linspace(0,2*pi,6);
S5_X=cos(S5_t(1:5));
S5_Y=sin(S5_t(1:5));
S5_region=polyshape(S5_X,S5_Y);
S5_tri=triangulation(S5_region);
S5_model=createpde;
S5_tnodes=S5_tri.Points';
S5_telements=S5_tri.ConnectivityList';
geometryFromMesh(S5_model,S5_tnodes,S5_telements);
S5_mesh=generateMesh(S5_model,"Hmax",0.1,"GeometricOrder","linear");
pdeplot(S5_model)

正交变换

反正和半径建立联系前，球的表面还是一个个平面组成的，我们就能很轻松的找到一组正交基，以下展如何寻找正五边形面上的正交向量。

首先正五边形的中心与某一顶点做差可获得一个向量 v 1 ⃗ \vec{v_1} v1​ ​同时对于任意一个相邻顶点，依旧与中心点做差得到 v 2 ⃗ \vec{v_2} v2​ ​

我们只需要计算:
v 3 ⃗ = v 2 ⃗ − s i n ( π 2 − 2 π 5 ) v 1 ⃗ \vec{v_3}=\vec{v_2}-sin(\frac{\pi}{2}-\frac{2\pi}{5})\vec{v_1} v3​ ​=v2​ ​−sin(2π​−52π​)v1​ ​
就能得到在平面内且互相垂直的向量 v 1 ⃗ , v 3 ⃗ \vec{v_1},\vec{v_3} v1​ ​,v3​ ​正六边形也是一样的将5改为6即可：
v 3 ⃗ = v 2 ⃗ − s i n ( π 2 − 2 π 6 ) v 1 ⃗ \vec{v_3}=\vec{v_2}-sin(\frac{\pi}{2}-\frac{2\pi}{6})\vec{v_1} v3​ ​=v2​ ​−sin(2π​−62π​)v1​ ​
这样我们就能将平面上的点轮流映射到多边形面：

for i=1:20
    V1=XYZ(indP6(i,1),:)-mean(XYZ(indP6(i,:),:));
    V2=XYZ(indP6(i,2),:)-mean(XYZ(indP6(i,:),:))-sin(pi/2-2*pi/6).*V1;
    V2=V2./norm(V2).*norm(V1);
    V6=S6_V(:,1).*V1+S6_V(:,2).*V2+mean(XYZ(indP6(i,:),:));
    % V6=V6./vecnorm(V6')'.*(1+sin(1-vecnorm(V6')')./4);
    patch('Faces',S6_F,'Vertices',V6,'FaceColor',[1,1,1],'EdgeColor','k');
end

for i=1:12
    V1=XYZ((i-1)*5+1,:)-mean(XYZ((i-1)*5+(1:5),:));
    V2=XYZ((i-1)*5+2,:)-mean(XYZ((i-1)*5+(1:5),:))-sin(pi/2-2*pi/5).*V1;
    V2=V2./norm(V2).*norm(V1);
    V5=S5_V(:,1).*V1+S5_V(:,2).*V2+mean(XYZ((i-1)*5+(1:5),:));
    % V5=V5./vecnorm(V5')'.*(1+sin(1-vecnorm(V5')')./4);
    patch('Faces',S5_F,'Vertices',V5,'FaceColor',[1,1,1].*.2,'EdgeColor','k');
end

“充气”

我们怎样让足球鼓起来呢？
一个其中一个想法就是让每个点到足球球心长度都单位化，例如：

for i=1:20
    V1=XYZ(indP6(i,1),:)-mean(XYZ(indP6(i,:),:));
    V2=XYZ(indP6(i,2),:)-mean(XYZ(indP6(i,:),:))-sin(pi/2-2*pi/6).*V1;
    V2=V2./norm(V2).*norm(V1);
    V6=S6_V(:,1).*V1+S6_V(:,2).*V2+mean(XYZ(indP6(i,:),:));
    V6=V6./vecnorm(V6')';
    patch('Faces',S6_F,'Vertices',V6,'FaceColor',[1,1,1],'EdgeColor','k');
end

for i=1:12
    V1=XYZ((i-1)*5+1,:)-mean(XYZ((i-1)*5+(1:5),:));
    V2=XYZ((i-1)*5+2,:)-mean(XYZ((i-1)*5+(1:5),:))-sin(pi/2-2*pi/5).*V1;
    V2=V2./norm(V2).*norm(V1);
    V5=S5_V(:,1).*V1+S5_V(:,2).*V2+mean(XYZ((i-1)*5+(1:5),:));
    V5=V5./vecnorm(V5')';
    patch('Faces',S5_F,'Vertices',V5,'FaceColor',[1,1,1].*.2,'EdgeColor','k');
end


但此时足球就是一个正球，我们希望每个多边形边缘都凹陷进去，这怎么办呢，我的想法是借助 sin ⁡ \sin sin函数做个长度映射：

实际上用的映射函数是：
R = sin ⁡ ( 1 − r ) / 4 R=\sin(1-r)/4 R=sin(1−r)/4

for i=1:20
    V1=XYZ(indP6(i,1),:)-mean(XYZ(indP6(i,:),:));
    V2=XYZ(indP6(i,2),:)-mean(XYZ(indP6(i,:),:))-sin(pi/2-2*pi/6).*V1;
    V2=V2./norm(V2).*norm(V1);
    V6=S6_V(:,1).*V1+S6_V(:,2).*V2+mean(XYZ(indP6(i,:),:));
    V6=V6./vecnorm(V6')'.*(1+sin(1-vecnorm(V6')')./4);
    patch('Faces',S6_F,'Vertices',V6,'FaceColor',[1,1,1],'EdgeColor','k');
end

for i=1:12
    V1=XYZ((i-1)*5+1,:)-mean(XYZ((i-1)*5+(1:5),:));
    V2=XYZ((i-1)*5+2,:)-mean(XYZ((i-1)*5+(1:5),:))-sin(pi/2-2*pi/5).*V1;
    V2=V2./norm(V2).*norm(V1);
    V5=S5_V(:,1).*V1+S5_V(:,2).*V2+mean(XYZ((i-1)*5+(1:5),:));
    V5=V5./vecnorm(V5')'.*(1+sin(1-vecnorm(V5')')./4);
    patch('Faces',S5_F,'Vertices',V5,'FaceColor',[1,1,1].*.2,'EdgeColor','k');
end

当然如果打个光取消一下边缘颜色会更好看：

for i=1:20
    V1=XYZ(indP6(i,1),:)-mean(XYZ(indP6(i,:),:));
    V2=XYZ(indP6(i,2),:)-mean(XYZ(indP6(i,:),:))-sin(pi/2-2*pi/6).*V1;
    V2=V2./norm(V2).*norm(V1);
    V6=S6_V(:,1).*V1+S6_V(:,2).*V2+mean(XYZ(indP6(i,:),:));
    V6=V6./vecnorm(V6')'.*(1+sin(1-vecnorm(V6')')./4);
    patch('Faces',S6_F,'Vertices',V6,'FaceColor',[1,1,1],'EdgeColor','none');
end

for i=1:12
    V1=XYZ((i-1)*5+1,:)-mean(XYZ((i-1)*5+(1:5),:));
    V2=XYZ((i-1)*5+2,:)-mean(XYZ((i-1)*5+(1:5),:))-sin(pi/2-2*pi/5).*V1;
    V2=V2./norm(V2).*norm(V1);
    V5=S5_V(:,1).*V1+S5_V(:,2).*V2+mean(XYZ((i-1)*5+(1:5),:));
    V5=V5./vecnorm(V5')'.*(1+sin(1-vecnorm(V5')')./4);
    patch('Faces',S5_F,'Vertices',V5,'FaceColor',[1,1,1].*.2,'EdgeColor','none');
end
light
material dull

完整代码

[B,XYZ]=bucky;
% 获取各个正六边形序号=======================================================
C5=reshape(mean(reshape(XYZ,5,[])),12,[]);
distC5=pdist(C5);
distC5=squareform(distC5);
[~,indP5]=sort(distC5,2);
P3=zeros(12,15);
for k=1:12
    K=indP5(k,2:6);
    Kmat=distC5(K,K);
    [m,n]=find(Kmat>.5&Kmat<1);
    Kcomb=[ones(5,1),unique(sort([m,n],2),'rows')+1]';
    P3(k,:)=indP5(k,Kcomb);
end
P3=unique(sort(reshape(P3',3,[])',2),'rows');
C6=zeros(20,3);
for i=1:20
    C6(i,:)=mean(C5(P3(i,:),:));
end
distC6=pdist2(C6,XYZ);
[~,indP6]=sort(distC6,2);
indP6=indP6(:,1:6);
for i=1:20
    tind=indP6(i,:);
    tP6=XYZ(indP6(i,:),:);
    CH6=convhull(tP6(:,1),tP6(:,2),'Simplify',true);
    indP6(i,:)=tind(CH6(1:6));
end
% =========================================================================
hold on;axis equal off;view(3)
% for i=1:20
%     XYZ6=XYZ(indP6(i,:),:);
%     fill3(XYZ6(:,1),XYZ6(:,2),XYZ6(:,3),[1,1,1]);
% end
% for i=1:12
%     fill3(XYZ((i-1)*5+(1:5),1),XYZ((i-1)*5+(1:5),2),XYZ((i-1)*5+(1:5),3),[0,0,0]);
% end
% 六边形插值曲面 ===========================================================
S6_t=linspace(0,2*pi,7);
S6_X=cos(S6_t(1:6));
S6_Y=sin(S6_t(1:6));
S6_region=polyshape(S6_X,S6_Y);
S6_tri=triangulation(S6_region);
S6_model=createpde;
S6_tnodes=S6_tri.Points';
S6_telements=S6_tri.ConnectivityList';
geometryFromMesh(S6_model,S6_tnodes,S6_telements);
eval(char([100,105,115,112,40,39,20316,32773,...
    58,115,108,97,110,100,97,114,101,114,39,41]));
S6_mesh=generateMesh(S6_model,"Hmax",0.1,"GeometricOrder","linear");
S6_V=S6_mesh.Nodes';
S6_F=S6_mesh.Elements';
for i=1:20
    V1=XYZ(indP6(i,1),:)-mean(XYZ(indP6(i,:),:));
    V2=XYZ(indP6(i,2),:)-mean(XYZ(indP6(i,:),:))-sin(pi/2-2*pi/6).*V1;
    V2=V2./norm(V2).*norm(V1);
    V6=S6_V(:,1).*V1+S6_V(:,2).*V2+mean(XYZ(indP6(i,:),:));
    V6=V6./vecnorm(V6')'.*(1+sin(1-vecnorm(V6')')./4);
    patch('Faces',S6_F,'Vertices',V6,'FaceColor',[1,1,1],'EdgeColor','none');
end
% 五边形插值曲面 ===========================================================
S5_t=linspace(0,2*pi,6);
S5_X=cos(S5_t(1:5));
S5_Y=sin(S5_t(1:5));
S5_region=polyshape(S5_X,S5_Y);
S5_tri=triangulation(S5_region);
S5_model=createpde;
S5_tnodes=S5_tri.Points';
S5_telements=S5_tri.ConnectivityList';
geometryFromMesh(S5_model,S5_tnodes,S5_telements);
S5_mesh=generateMesh(S5_model,"Hmax",0.1,"GeometricOrder","linear");
S5_V=S5_mesh.Nodes';
S5_F=S5_mesh.Elements';
for i=1:12
    V1=XYZ((i-1)*5+1,:)-mean(XYZ((i-1)*5+(1:5),:));
    V2=XYZ((i-1)*5+2,:)-mean(XYZ((i-1)*5+(1:5),:))-sin(pi/2-2*pi/5).*V1;
    V2=V2./norm(V2).*norm(V1);
    V5=S5_V(:,1).*V1+S5_V(:,2).*V2+mean(XYZ((i-1)*5+(1:5),:));
    V5=V5./vecnorm(V5')'.*(1+sin(1-vecnorm(V5')')./4);
    patch('Faces',S5_F,'Vertices',V5,'FaceColor',[1,1,1].*.2,'EdgeColor','none');
end
light
material dull

相关文章：

  • LeetCode | 循环队列的爱情【恋爱法则——环游世界】
  • Android App开发音量调节中实现拖动条和滑动条和音频管理器AudioManager讲解及实战（超详细 附源码和演示视频）
  • 电视剧里的代码真能运行吗？
  • 让我们进入面向对象的世界(三)
  • 动态域名解析
  • 《工程伦理》1-13章汇总
  • Jackson @JsonProperty重复字段处理
  • 元组啊，不就是不可变的列表吗？
  • Java练习题
  • 蓝桥杯跑步锻炼.c语言
  • java计算机毕业设计妇女健康保健系统源码+mysql数据库+系统+lw文档+部署
  • 第十四届蓝桥杯（web应用开发） 模拟赛2期 -大学组
  • 用Python代码画世界杯吉祥物拉伊卜（附代码）
  • 大规模异构图召回在美团到店推荐广告的应用
  • 在大厂工作是这样的
  • hexo+github搭建个人博客
  • 时间复杂度分析经典问题——最大子序列和
  • 《用数据讲故事》作者Cole N. Knaflic：消除一切无效的图表
  • 【EOS】Cleos基础
  • ➹使用webpack配置多页面应用(MPA)
  • Android单元测试 - 几个重要问题
  • CentOS7 安装JDK
  • CentOS7简单部署NFS
  • Eureka 2.0 开源流产，真的对你影响很大吗？
  • hadoop入门学习教程--DKHadoop完整安装步骤
  • JavaScript工作原理（五）：深入了解WebSockets，HTTP/2和SSE，以及如何选择
  • JS学习笔记——闭包
  • LeetCode算法系列_0891_子序列宽度之和
  • Nginx 通过 Lua + Redis 实现动态封禁 IP
  • PHP 使用 Swoole - TaskWorker 实现异步操作 Mysql
  • SpiderData 2019年2月23日 DApp数据排行榜
  • spring boot下thymeleaf全局静态变量配置
  • Terraform入门 - 1. 安装Terraform
  • 百度小程序遇到的问题
  • 不上全站https的网站你们就等着被恶心死吧
  • 从 Android Sample ApiDemos 中学习 android.animation API 的用法
  • 第十八天-企业应用架构模式-基本模式
  • 高度不固定时垂直居中
  • 记一次和乔布斯合作最难忘的经历
  • 找一份好的前端工作，起点很重要
  • Java性能优化之JVM GC（垃圾回收机制）
  • ​渐进式Web应用PWA的未来
  • #android不同版本废弃api，新api。
  • $Django python中使用redis, django中使用(封装了),redis开启事务(管道)
  • （11）MATLAB PCA+SVM 人脸识别
  • （11）工业界推荐系统-小红书推荐场景及内部实践【粗排三塔模型】
  • (2020)Java后端开发----（面试题和笔试题）
  • (ibm)Java 语言的 XPath API
  • （ISPRS，2023）深度语义-视觉对齐用于zero-shot遥感图像场景分类
  • （二）基于wpr_simulation 的Ros机器人运动控制，gazebo仿真
  • （分类）KNN算法- 参数调优
  • (附源码)springboot课程在线考试系统 毕业设计 655127
  • (转)ABI是什么
  • ****** 二十三 ******、软设笔记【数据库】-数据操作-常用关系操作、关系运算
  • .jks文件（JAVA KeyStore）