0871-64605728
您当前位置:网站首页 >> 知识专区
七夕节快到了,教你用MATLAB绘制blingbling的大钻石
文章来源:CSDN slandarer  上传时间:2021-8-13  浏览量:635

版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。



1效果

在这里插入图片描述
在这里插入图片描述

2程序说明

其实非常简单,只需要计算出各个顶点位置,并用fill3绘制一个个面出来即可,为了钻石的色彩更加丰富,我们用light函数设置了光照并为fill3创建出的PATCH对象设置了以下属性。

  • FaceAlpha 面透明度
  • EdgeAlpha 边透明度
  • FaceLighting 面照亮方式(是整个面相同颜色还是通过插值不同位置不同颜色)
  • EdgeLighting 边照亮方式(与面照亮方式类似)
  • SpecularStrength 镜面反射的强度
  • DiffuseStrength 散射光的强度
  • AmbientStrength 环境光的强度

其中:
若将PATCH属性由:

faceColor=[0 71 177]; 
facePara={
faceColor./255,'FaceAlpha',0.6,'EdgeColor',faceColor./255.*0.8,'EdgeAlpha',0.3,... 'FaceLighting','gouraud',
'EdgeLighting','gouraud',... 'SpecularStrength',0.9,'DiffuseStrength',0.7,'AmbientStrength',0.7}; % 面参数
	

	

改为:

faceColor=[237 64 153]; facePara={faceColor./255,'FaceAlpha',0.4,'EdgeColor',faceColor./255.*0.8,
'EdgeAlpha',0.2,... 'FaceLighting','gouraud','EdgeLighting','gouraud',... 'SpecularStrength',0.9,'DiffuseStrength',0.7,'AmbientStrength',0.6}; % 面参数
	

	

则钻石变成这个样子:
在这里插入图片描述
另:
大家可以猜猜下面这个钻石是咋画的:
在这里插入图片描述
在这里插入图片描述

3完整代码

function diamond
hold on
ax=gca; axis equal
grid on ax.GridLineStyle='--'; 
ax.LineWidth=1.2; ax.XLim=[-2,2]; ax.YLim=[-2,2]; ax.ZLim=[-1.8,1.5]; 
ax.XColor=[1,1,1].*0.4; ax.YColor=[1,1,1].*0.4; ax.ZColor=[1,1,1].*0.4; 
ax.DataAspectRatio=[1,1,1]; ax.DataAspectRatioMode='manual'; ax.CameraPosition=[-67.6287 -204.5276 82.7879]; 
light('Position',[5,6,5]) thetaSet8=linspace(0,2*pi-pi/4,8)'; 
thetaSet16=linspace(0,2*pi-pi/8,16)'; faceColor=[0 71 177]; 
facePara={faceColor./255,'FaceAlpha',0.6,'EdgeColor',faceColor./255.*0.8,'EdgeAlpha',0.3,... 'FaceLighting','gouraud',
'EdgeLighting','gouraud',... 'SpecularStrength',0.9,'DiffuseStrength',0.7,'AmbientStrength',0.7}; % 面参数 % 
faceColor=[237 64 153]; 
% facePara={faceColor./255,'FaceAlpha',0.4,'EdgeColor',faceColor./255.*0.8,'EdgeAlpha',0.2,...
 % 'FaceLighting','gouraud','EdgeLighting','gouraud',... % 'SpecularStrength',0.9,'DiffuseStrength',0.7,'AmbientStrength',0.6}; % 面参数 % 绘制最上方八边形
pntSet1=[cos(thetaSet8),sin(thetaSet8),ones(8,1)]; fill3(pntSet1(:,1),pntSet1(:,2),pntSet1(:,3),facePara{:}) % 一圈8个三角形
pntSet2=[cos(thetaSet8+pi/8).*1.3,sin(thetaSet8+pi/8).*1.3,ones(8,1).*0.7]; for i=1:8 p1=i; p2=mod(i,8)+1;
 p3=i;
 fill3([pntSet1([p1,p2],1);pntSet2(p3,1)],... [pntSet1([p1,p2],2);pntSet2(p3,2)],... [pntSet1([p1,p2],3);pntSet2(p3,3)],facePara{:}) 
end % 一圈8个四边形
pntSet3=[cos(thetaSet8).*1.7,sin(thetaSet8).*1.7,ones(8,1).*0.3]; 
for i=1:8 p1=mod(i,8)+1; p2=i; p3=p1; p4=p1; fill3([pntSet1(p1,1);
pntSet2(p2,1);pntSet3(p3,1);pntSet2(p4,1)],... [pntSet1(p1,2);pntSet2(p2,2);pntSet3(p3,2);pntSet2(p4,2)],... [pntSet1(p1,3);pntSet2(p2,3);pntSet3(p3,3);
pntSet2(p4,3)],facePara{:}) end % 一圈16个三角形
pntSet4=[cos(thetaSet16).*1.7,sin(thetaSet16).*1.7,ones(16,1).*0.3]; 
for i=1:8 p1=i; p2=2*i-1; p3=2*i; fill3([pntSet2(p1,1);
pntSet4([p2,p3],1)],... [pntSet2(p1,2);pntSet4([p2,p3],2)],... [pntSet2(p1,3);pntSet4([p2,p3],3)],facePara{:})
 p2=2*i; p3=mod(2*i,16)+1; 
fill3([pntSet2(p1,1);pntSet4([p2,p3],1)],... [pntSet2(p1,2);pntSet4([p2,p3],2)],... [pntSet2(p1,3);pntSet4([p2,p3],3)],facePara{:}) end % 一圈16个四边角形
pntSet5=[cos(thetaSet16).*1.7,sin(thetaSet16).*1.7,ones(16,1).*0.15]; for i=1:16 p1=i;p2=mod(i,16)+1; 
p3=mod(i,16)+1;p4=i; fill3([pntSet4(p1,1);pntSet4(p2,1);pntSet5(p3,1);
pntSet5(p4,1)],... [pntSet4(p1,2);pntSet4(p2,2);pntSet5(p3,2);pntSet5(p4,2)],... [pntSet4(p1,3);pntSet4(p2,3);pntSet5(p3,3);pntSet5(p4,3)],facePara{:}) end
pntSet6=[cos(thetaSet8+pi/8).*0.4,sin(thetaSet8+pi/8).*0.4,ones(8,1).*(-1.15)]; % 一圈16个长三角形
 for i=1:8 p1=i; p2=2*i-1; p3=2*i; fill3([pntSet6(p1,1);
pntSet5([p2,p3],1)],... [pntSet6(p1,2);pntSet5([p2,p3],2)],... [pntSet6(p1,3);pntSet5([p2,p3],3)],facePara{:}) 
p2=2*i; p3=mod(2*i,16)+1; fill3([pntSet6(p1,1);pntSet5([p2,p3],1)],... [pntSet6(p1,2);pntSet5([p2,p3],2)],... [pntSet6(p1,3);
pntSet5([p2,p3],3)],facePara{:}) end % 一圈8个长四边形
pntSet7=[cos(thetaSet8).*1.7,sin(thetaSet8).*1.7,ones(8,1).*0.15]; 
for i=1:8 p1=mod(i,8)+1; p2=i; p4=mod(i,8)+1; fill3([pntSet7(p1,1);
pntSet6(p2,1);0;pntSet6(p4,1)],... [pntSet7(p1,2);pntSet6(p2,2);0;pntSet6(p4,2)],... [pntSet7(p1,3);pntSet6(p2,3);-1.5;pntSet6(p4,3)],facePara{:}) end
end

4其他过节可能会用到的程序

↓ ↓ ↓ ↓ ↓ ↓ ↓ ↓ ↓ ↓ ↓ ↓ ↓ ↓ ↓ ↓ ↓ ↓ ↓ ↓ ↓
!!!!点击图片跳转连接!!!!
↑ ↑ ↑ ↑ ↑ ↑ ↑ ↑ ↑ ↑ ↑ ↑ ↑ ↑ ↑ ↑ ↑ ↑ ↑ ↑ ↑

4.1玫瑰花球

4.2樱花树

4.3立体爱心九宫格

4.4立体玫瑰花

4.5水晶爱心

4.6两款相册

-1紫色钻石

抱歉大家我确实不记得具体的参数了,不过下面这份代码大家可以试一下:
这份代码的效果:
在这里插入图片描述
在这里插入图片描述

需要注意的是,很多角度下并不会产生中间色,只会显示原本的粉色或紫色,只有一部分角度会有很绚丽的颜色,
另:该份代码中对光源的设置设置了两遍:

light('Position',[5,6,5]) 

这样确实会使钻石某些反射面亮度更高更绚丽。
以下是绘制两遍钻石的代码,大家可以对各个部分数据进行调整以获得更绚丽的效果:

function diamondx2
hold on
ax=gca; axis equal
grid on ax.GridLineStyle='--'; 
ax.LineWidth=1.2; ax.XLim=[-2,2]; 
ax.YLim=[-2,2]; ax.ZLim=[-1.8,1.5]; 
ax.XColor=[1,1,1].*0.4; 
ax.YColor=[1,1,1].*0.4; 
ax.ZColor=[1,1,1].*0.4;
 ax.DataAspectRatio=[1,1,1];
 ax.DataAspectRatioMode='manual'; 
ax.CameraPosition=[-67.6287 -204.5276 82.7879];
 light('Position',[5,6,5]) light('Position',[5,6,5]) faceColor=[0 71 177]; 
facePara={
faceColor./255,'FaceAlpha',0.4,'EdgeColor',faceColor./255.*0.8,'EdgeAlpha',0.2,... 'FaceLighting',
'gouraud','EdgeLighting','gouraud',... 'SpecularStrength',0.9,'DiffuseStrength',0.8,'AmbientStrength',0.7}; % 面参数 
drawDiamond(facePara) faceColor=[237 64 153]; facePara={
faceColor./255,'FaceAlpha',0.4,'EdgeColor',faceColor./255.*0.8,'EdgeAlpha',0.2,... 'FaceLighting','gouraud',
'EdgeLighting','gouraud',... 'SpecularStrength',0.9,'DiffuseStrength',0.8,'AmbientStrength',0.7}; % 面参数
 drawDiamond(facePara) function drawDiamond(facePara) thetaSet8=linspace(0,2*pi-pi/4,8)'; thetaSet16=linspace(0,2*pi-pi/8,16)'; % 绘制最上方八边形
        pntSet1=[cos(thetaSet8),sin(thetaSet8),ones(8,1)]; fill3(pntSet1(:,1),pntSet1(:,2),pntSet1(:,3),facePara{:}) % 一圈8个三角形
        pntSet2=[cos(thetaSet8+pi/8).*1.3,sin(thetaSet8+pi/8).*1.3,ones(8,1).*0.7]; for i=1:8 p1=i; p2=mod(i,8)+1;
 p3=i; 
fill3([pntSet1([p1,p2],1);
pntSet2(p3,1)],... [pntSet1([p1,p2],2);pntSet2(p3,2)],... [pntSet1([p1,p2],3);pntSet2(p3,3)],facePara{:}) end % 一圈8个四边形
        pntSet3=[cos(thetaSet8).*1.7,sin(thetaSet8).*1.7,ones(8,1).*0.3]; 
for i=1:8 p1=mod(i,8)+1; p2=i; p3=p1; p4=p1; 
fill3([pntSet1(p1,1);
pntSet2(p2,1);pntSet3(p3,1);pntSet2(p4,1)],... [pntSet1(p1,2);pntSet2(p2,2);pntSet3(p3,2);pntSet2(p4,2)],... [pntSet1(p1,3);pntSet2(p2,3);
pntSet3(p3,3);pntSet2(p4,3)],facePara{:}) end % 一圈16个三角形
        pntSet4=[cos(thetaSet16).*1.7,sin(thetaSet16).*1.7,ones(16,1).*0.3]; 
for i=1:8 p1=i; p2=2*i-1; p3=2*i; 
fill3([pntSet2(p1,1);pntSet4([p2,p3],1)],... [pntSet2(p1,2);pntSet4([p2,p3],2)],... [pntSet2(p1,3);
pntSet4([p2,p3],3)],facePara{:}) p2=2*i; p3=mod(2*i,16)+1; fill3([pntSet2(p1,1);
pntSet4([p2,p3],1)],... [pntSet2(p1,2);pntSet4([p2,p3],2)],... [pntSet2(p1,3);pntSet4([p2,p3],3)],facePara{:}) end % 一圈16个四边角形
        pntSet5=[cos(thetaSet16).*1.7,sin(thetaSet16).*1.7,ones(16,1).*0.15]; 
for i=1:16 p1=i;p2=mod(i,16)+1; p3=mod(i,16)+1;p4=i;
 fill3([pntSet4(p1,1);pntSet4(p2,1);pntSet5(p3,1);pntSet5(p4,1)],... [pntSet4(p1,2);pntSet4(p2,2);
pntSet5(p3,2);pntSet5(p4,2)],... [pntSet4(p1,3);pntSet4(p2,3);pntSet5(p3,3);pntSet5(p4,3)],facePara{:}) end
        pntSet6=[cos(thetaSet8+pi/8).*0.4,sin(thetaSet8+pi/8).*0.4,ones(8,1).*(-1.15)]; % 一圈16个长三角形
 for i=1:8 p1=i; p2=2*i-1; p3=2*i; 
fill3([pntSet6(p1,1);pntSet5([p2,p3],1)],... [pntSet6(p1,2);pntSet5([p2,p3],2)],... [pntSet6(p1,3);
pntSet5([p2,p3],3)],facePara{:}) p2=2*i; p3=mod(2*i,16)+1; fill3([pntSet6(p1,1);
pntSet5([p2,p3],1)],... [pntSet6(p1,2);pntSet5([p2,p3],2)],... [pntSet6(p1,3);pntSet5([p2,p3],3)],facePara{:}) end % 一圈8个长四边形
        pntSet7=[cos(thetaSet8).*1.7,sin(thetaSet8).*1.7,ones(8,1).*0.15]; 
for i=1:8 p1=mod(i,8)+1; p2=i; p4=mod(i,8)+1; fill3([pntSet7(p1,1);pntSet6(p2,1);0;pntSet6(p4,1)],... [pntSet7(p1,2)
;pntSet6(p2,2);0;pntSet6(p4,2)],... [pntSet7(p1,3);pntSet6(p2,3);-1.5;pntSet6(p4,3)],facePara{:}) end
    end
end



26

2021-04

轻松带你学习java-agent

java-agent的挂载方式有两种,一种是静态挂载,一种是动态挂载。静态挂载中,agent与java应用一起启动,在java应用初始化前agent就已经挂载完成,并开始监控java应用。动态挂载则是在应用运行过程中,通过进程ID确定挂载对象,动态的将agent挂载在目标进程上。

25

2021-01

拒绝云服务商白嫖,Elasticsearch 和 Kibana 变更开源许可协议

Elasticsearch 和 Kibana 变更开源许可协议

07

2022-02

2022北京冬奥会开幕式里的黑科技,闪耀闪耀全世界

2022北京冬奥会开幕式里的黑科技,闪耀闪耀全世界

07

2022-02

Windows的权限(用户、组和访问控制)

Windows的权限(用户、组和访问控制)

03

2021-09

❤️五万字《十大排序算法》动图讲解❤️

❤️五万字《十大排序算法》动图讲解❤️

27

2021-08

HTML5 标签大全笔记,用到的时候翻一下,百度都省了!

HTML5 标签大全笔记,用到的时候翻一下,百度都省了!

15

2021-07

MySql知识体系总结(SQL优化篇)

MySql知识体系总结(SQL优化篇)

10

2022-06

云原生是实现可观测平台的唯一出路?码农:夸张了

云原生是实现可观测平台的唯一出路?码农:夸张了
返回顶部
客服电话
0871-64605728
用微信扫一扫关注我们
请各公司推销人员注意:我单位拒绝任何方式、任何形式的电话推销,请勿拔打我单位客服热线进行电话推销,谢谢合作!
公司名称:云南昂略科技有限公司
联系地址:云南省昆明市官渡区永平路188号鑫都韵城写字楼6栋1004号
联系电话:0871-64605728、传真号码:0871-64605728
电子邮箱:19701580@qq.com
关键词:知识专区:七夕节快到了,教你用MATLAB绘制blingbling的大钻石,云南昂略科技有限公司,云南移动执法平台建设,云南智慧安防调度系统,云南头戴式安全终端,昂略科技
云南网站建设,云南网页设计,昆明网站建设,昆明网页设计  网站管理
【版权声明】本站部分内容由互联网用户自行发布,著作权或版权归原作者所有。如果侵犯到您的权益请发邮件致info@ynjwz.com,我们会第一时间进行删除并表示歉意。