MATLAB | 情人节来绘制更立体的玫瑰花吧
MATLAB 绘制 立体 情人节
2023-09-11 14:20:18 时间
又是一年情人节,今年带来一款更有立体感的玫瑰:
曲面的函数表达式来自:
http://www.bugman123.com/Math/index.html
这个网站,上面还有很多其他帅气的玩意。
基础绘制
x=linspace(0,1,300);
theta=linspace(-2*pi,15*pi,300);
[x,theta]=meshgrid(x,theta);
phi=(pi/2).*exp(-theta./8./pi);
X=1-.5.*(1.25.*(1-mod(3.6.*theta,2*pi)./pi).^2-1/4).^2;
y=1.95653.*x.^2.*(1.27689.*x-1).^2.*sin(phi);
r=X.*(x.*sin(phi)+y.*cos(phi));
roseHdl=surf(r.*cos(theta),r.*sin(theta),X.*(x.*cos(phi)-y.*sin(phi)),'EdgeColor','none');
坐标区域修饰
在最后加入以下代码能让绘图更好看一些:
% 坐标区域修饰
ax=gca;hold on;grid on;
axis([-1,1,-1,1,-.5,1])
ax.FontName='Cambria';
ax.LineWidth=1;
ax.GridLineStyle='-.';
ax.Projection='perspective';
ax.XMinorTick='on';
ax.YMinorTick='on';
ax.ZMinorTick='on';
配色
可以调整颜色映射方向:
比如根据半径映射:
roseHdl.CData=r;
根据x轴坐标大小映射:
roseHdl.CData=r.*cos(theta);
配色可以自己弄点数值矩阵插值,比如:
roseHdl.CData=r;
CM=[0.5300 0.8300 0.8100
0.5200 0.7500 0.8200
0.4900 0.6200 0.8400
0.4900 0.5600 0.8400
0.4700 0.4900 0.8500
0.4500 0.3500 0.8700
0.9500 0.9500 0.9500];
CMX=linspace(0,1,size(CM,1));
CMXX=linspace(0,1,256)';
CM=[interp1(CMX,CM(:,1),CMXX,'pchip'),interp1(CMX,CM(:,2),CMXX,'pchip'),interp1(CMX,CM(:,3),CMXX,'pchip')];
colormap(CM)
还可以配合之前写的slanCM工具包:
https://blog.csdn.net/slandarer/article/details/127719784
随便举点例子:
roseHdl.CData=r.*cos(theta);
CM=slanCM('blues');
colormap(CM)
roseHdl.CData=r;
CM=slanCM(134);
colormap(CM)
roseHdl.CData=r;
CM=slanCM(136);
colormap(CM)
旋转
把代码改成这样就能一直旋转:
function valentinesRose2
% 绘制玫瑰花
x=linspace(0,1,300);
theta=linspace(-2*pi,15*pi,300);
[x,theta]=meshgrid(x,theta);
phi=(pi/2).*exp(-theta./8./pi);
X=1-.5.*(1.25.*(1-mod(3.6.*theta,2*pi)./pi).^2-1/4).^2;
y=1.95653.*x.^2.*(1.27689.*x-1).^2.*sin(phi);
r=X.*(x.*sin(phi)+y.*cos(phi));
roseHdl=surf(r.*cos(theta),r.*sin(theta),X.*(x.*cos(phi)-y.*sin(phi)),'EdgeColor','none');
roseHdl.CData=r;
CM=slanCM('copper2');
CM=CM(1:180,:);
colormap(CM)
% 坐标区域修饰
ax=gca;hold on;grid on;
axis([-1,1,-1,1,-.5,1])
ax.FontName='Cambria';
ax.LineWidth=1;
ax.GridLineStyle='-.';
ax.Projection='perspective';
ax.XMinorTick='on';
ax.YMinorTick='on';
ax.ZMinorTick='on';
set(gcf,'Color',[1,1,1]);
% 循环绘图旋转起来
while true
theta=theta+.01;
roseHdl.XData=r.*cos(theta);
roseHdl.YData=r.*sin(theta);
pause(.01),drawnow
end
若想自动保存为gif,把最后部分旋转代码改为如下部分即可:
% 存储gif =================================================================
% R2022a及之后版本
n=0;
while true
theta=theta+.01;
roseHdl.XData=r.*cos(theta);
roseHdl.YData=r.*sin(theta);
if n<50
exportgraphics(gcf,'test1.gif','Append',true)
end
n=n+1;
pause(.01),drawnow
end
% R2022a之前版本
% n=0;DelayTime=.02;
% F=getframe(ax);
% [imind,cm]=rgb2ind(F.cdata,256);
% imwrite(imind,cm,'test2.gif','gif','Loopcount',inf,'DelayTime',DelayTime);
% while true
% theta=theta+.01;
% roseHdl.XData=r.*cos(theta);
% roseHdl.YData=r.*sin(theta);
% if n<50
% F=getframe(ax);
% [imind,cm]=rgb2ind(F.cdata,256);
% imwrite(imind,cm,'test2.gif','gif','WriteMode','append','DelayTime',DelayTime);
% end
% n=n+1;
% pause(.01),drawnow
% end
相关文章
- Matlab中writetable函数的使用
- Matlab中abs函数的使用
- Matlab中查找函数参数的数量
- Matlab中gray2ind函数用法
- Matlab中绘制日期和持续时间图
- matlab基本函数randperm end数组索引
- 基于MATLAB的Okumura-Hata模型的仿真
- 【MATLAB教程案例82】matlab在大学数学中的应用——概率统计
- 【MATLAB教程案例64】深度学习网络建模1,通过matlab编程实现
- 【MATLAB教程案例60】使用matlab实现基于GRU网络的数据分类预测功能与仿真分析
- 【MATLAB教程案例40】语音信号的共振峰频率倒谱法估计matlab仿真学习
- 【MATLAB教程案例13】基于SA模拟退火优化算法的函数极值计算matlab仿真及其他应用
- 【MATLAB教程案例8】基于LS算法的OFDM调制解调系统信道估计和均衡算法的matlab仿真
- 【MATLAB教程案例7】基于Gardner环的定时误差检测matlab仿真
- 【MATLAB教程案例6】基于Costas环的载波同步matlab仿真
- 【广义S变换】一维广义S变换对非平稳信号处理的matlab仿真
- 【短时幅度谱】短时幅度谱估计在语音增强方面的MATLAB仿真
- 基于matlab程序hough变换检测圆源码
- MATLAB安装机器人工具箱 Matlab_Robotic_Toolbox_v10.2
- 《数字图像处理与机器视觉——Visual C++与Matlab实现(第2版)》——2.3 MATLAB的图像转换
- 【Matlab 六自由度机器人】运动学正解(附MATLAB机器人正解完整代码)
- 【Matlab算法】粒子群算法求解二维线性优化问题(附MATLAB代码)
- 基于 MATLAB 实现 Word 的信息隐藏算法【100010678】
- MatLab教程之使用Excel和MATLAB求解工程中的线性和非线性方程
- MATLAB 制作抖音同款炫光海报
- MATLAB 中gcf、gca 以及gco 的区别
- Matlab绘制混淆矩阵
- Matlab自定义函数详解