MATLAB | 科研绘图第二十二期信号分解中用于展示效果的三维瀑布图

绘图教程及颜色数据链接见文末哦~

1 出图效果

好看的图幅是一篇SCI被录用的关键~

2 部分绘图代码

%% Made by Lwcah in 2023-07-09(公众号:Lwcah)
%% 公众号:Lwcah
%% 知乎、B站、小红书、视频号、抖音同名账号:Lwcah,感谢关注~
%% 更多MATLAB+SCI绘图教程敬请观看~

%% 绘制模板
close all;clear all;clc;

%% 1行1列
%% 一幅图的时候figureWidth = 8.5;figureHeight = 8;调整这两个参数就可以
%定义整幅图出现的在电脑屏幕上的位置以及长和宽
%这里有必要解释一下为什么figureWidth要设置为8.5;
%因详细解释需要很长的篇幅,请转公众号观看该天的文章。
figureHandle = figure;
figureUnits = 'centimeters';
figureWidth = 15;
figureHeight = 9;
set(gcf, 'Units', figureUnits, 'Position', [28 10 figureWidth figureHeight]);
% 注:28代表出图时图的左下角相对于整个电脑屏幕的左下角向左偏移28个单位,向上偏移20个单位。
% 可自行调节两个数字让图出在自己屏幕的某个位置

%% 定义字体和字号大小
%% 通常的SCI一般Arial字体和10字号通用
fontnamed='Arial';%字号名字
ssize=10;%字号大小
%% 如果是中文论文可以相应的更改字号名字如下
% '华文中宋' '华文仿宋' '华文宋体' '华文新魏' '华文楷体'
% '华文琥珀' '华文细黑' '华文行楷' '华文隶书' '宋体'
% '方正姚体' '微软雅黑' '方正舒体' '新宋体'
% '幼圆' '楷体' '等线' '隶书' '黑体'

%% 一、信号分解并出图
% 加载数据
load data              
x=data(:,1);
y=data(:,2);
% 展开信号分解
[imf,residual]=emd(y);
residual      = residual.*ones(size(y))';
modes         = [imf;residual];

%% 构造XYZ数据集
% X轴应该为Imf代表的分量数
[m,n] = size(modes);% m为几个IMF分量;n为分量的数据长度。
X1    = (1:1:m)';
X2    = ones(size(modes));
X     = X1.*X2;
X     = X';
% Y轴应该为Imf代表的分量数据长度
Y1    = (1:1:n);
Y2    = ones(size(modes));
Y     = Y1.*Y2;
Y     = Y';
% Z轴呈现的为信号分解的数据(为了使画图更好看所以要调整顺序)
Z     = modes;
Z     = [modes(8,:);modes(7,:);modes(6,:);modes(5,:);modes(4,:);modes(3,:);modes(2,:);modes(1,:)];

%% 绘图
plot3(X,Y,Z,'linewidth',1);
hTitle = title('Three-dimensional waterfall map');
hXLabel = xlabel('X');
hYLabel = ylabel('Y');
hZLabel = zlabel('Z');

%% 画图的标准格式代码
% text(0.5,0.9,'(a) XXX','fontsize',ssize,'FontName',fontnamed);%,'horiz','center'
xlabel('X-axis','fontsize',ssize,'FontName',fontnamed);
ylabel('Y-axis','fontsize',ssize,'FontName',fontnamed);
zlabel('Z-axis','fontsize',ssize,'FontName',fontnamed);
% axis([0 7 0 0.7 50 250]);%XYZ轴的范围
% xticks([1 2 3 4 5 6]);%画格网的时候的小刻度
% xticklabels({'1','2','3','4','5','6'});%加x轴刻度标注
% yticks([0.1 0.2 0.3 0.4 0.5 0.6]);%画格网的时候的小刻度
% yticklabels({'0.1','0.2','0.3','0.4','0.5','0.6'});%加y轴刻度标注
% zticks([50 100 150 200 250]);%画格网的时候的小刻度
% zticklabels({'50','100','150','200','250'});%加z轴刻度标注
set(gca,'linewidth',1,'fontsize',ssize,'FontName',fontnamed);
% set(gca,'zticklabel',[]);%z轴不显示
% set(gca,'yticklabel',[]);%y轴不显示
% set(gca,'xticklabel',[]);%x轴不显示
grid on;box on;hold on;

%% 画legend
%方法一
kk=legend('L8','L7','L6','L5','L4','L3','L2','L1');
set(kk,'location','NorthEast','Box', 'off','fontsize',ssize,'FontName',fontnamed);%'orientation','horizontal',
%方法二
% columnlegend(2,{'L1','L2','L3','L4','L5','L6'},'North');%表示一行放三个图例以及图例的位置
%% Matlab中有许多位置可以选择:
% 'North' inside plot box near top
% 'South' inside bottom
% 'East' inside right
% 'West' inside left
% 'NorthEast' inside top right (default for 2-D plots)
% 'NorthWest' inside top left
% 'SouthEast' inside bottom right
% 'SouthWest' inside bottom left
% 'NorthOutside' outside plot box near top
% 'SouthOutside' outside bottom
% 'EastOutside' outside right
% 'WestOutside' outside left
% 'NorthEastOutside' outside top right (default for 3-D plots)
% 'NorthWestOutside' outside top left
% 'SouthEastOutside' outside bottom right
% 'SouthWestOutside' outside bottom left
% 'Best' least conflict with data in plot 与绘图中的数据冲突最小
% 'BestOutside' least unused space outside plot

%% 背景颜色
set(gcf,'Color',[1 1 1])
%% 设置完毕后,按照所需分辨率、格式输出
figW = figureWidth;
figH = figureHeight;
set(figureHandle,'PaperUnits',figureUnits);
set(figureHandle,'PaperPosition',[0 0 figW figH]);
fileout = '三维瀑布图';
print(figureHandle,[fileout,'.png'],'-r600','-dpng');





%% 一幅图的时候figureWidth = 8.5;figureHeight = 8;调整这两个参数就可以
%定义整幅图出现的在电脑屏幕上的位置以及长和宽
%这里有必要解释一下为什么figureWidth要设置为8.5;
%因详细解释需要很长的篇幅,请转公众号观看2023-03-21的文章观看。
figureHandle = figure;
figureUnits  = 'centimeters';
figureWidth  = 8.5;
figureHeight = 14;
set(gcf, 'Units', figureUnits, 'Position', [20 12 figureWidth figureHeight]);
%% 展开信号分解
[imf,residual]=emd(y);
residual      = residual.*ones(size(y))';
modes         = [imf;residual];
%信号重构
d1           = modes(1,:);
d2           = modes(2,:);
d3           = modes(3,:);
d4           = modes(4,:);
d5           = modes(5,:);
d6           = modes(6,:);
d7           = modes(7,:);
d8           = modes(8,:);
% d9           = modes(9,:);
% d10          = modes(10,:);
%% 画图
subplot(9,1,1)
plot(x,y);
ylabel('原信号');hold on;
subplot(9,1,2)
plot(x,d1);
ylabel('d1');hold on;
subplot(9,1,3)
plot(x,d2);
ylabel('d2');hold on;
subplot(9,1,4)
plot(x,d3);
ylabel('d3');hold on;
subplot(9,1,5)
plot(x,d4);
ylabel('d4');hold on;
subplot(9,1,6)
plot(x,d5);
ylabel('d5');hold on;
subplot(9,1,7)
plot(x,d6);
ylabel('d6');hold on;
subplot(9,1,8)
plot(x,d7);
ylabel('d7');hold on;
subplot(9,1,9)
plot(x,d8);
ylabel('d8');hold on;
% subplot(11,1,10)
% plot(x,d9);
% ylabel('d9');hold on;
% subplot(11,1,11)
% plot(x,d10);
% ylabel('d10');hold on;
xlabel('X-axis');
%% 背景颜色
set(gcf,'Color',[1 1 1])
%% 图片输出
figW    = figureWidth;
figH    = figureHeight;
set(figureHandle,'PaperUnits',figureUnits);
set(figureHandle,'PaperPosition',[0 0 figW figH]);
fileout = '原始数据分解各分量';
print(figureHandle,[fileout,'.png'],'-r600','-dpng');


3 如何绘制满足 SCI 要求的单幅图

单幅图在我看来就是一行一列的子图,只是他的子图只有一个而已。持续关注后期会更新符合 SCI 要求的一行两列的子图、一行三列的子图、两行两列的子图、两行三列的子图、两行一列的子图、三行一列的子图、六行一列的子图,这是我目前 SCI 绘制过程中常遇到的图形绘制格式。

这里有必要解释一下为什么 figureWidth 要设置为 8.5,这是因为双栏排版的情况下,用于插入图的空间有限。

那么是如何计算的呢,A4 纸规格尺寸:210mm×297mm。A4 纸是由国际标准化组织 ISO 216 定义的,规格为 210mm×297mm,世界上多数国家所使用的纸张尺寸都是采用这一国际标准。

论文页边距的要求为:上 (T):25 mm;下 (B):25 mm;左 (L):20 mm;右 (R):20 mm。

也就是说论文须用 A4(210mm×297mm)标准大小白纸印制,要求纸的四周留足空白边缘,以便装订、复制和批注。每一面的上方(天头)和下方(地脚)应分别留边 25mm,左侧(订口)和右侧(切口)应分别留边 20mm。这是相对于一般的投稿期刊来说的。

那么 210mm- 20mm- 20mm=170mm,则双栏排版的话,留给单栏用于放置图片的空间仅剩 85mm,也即 8.5cm。只有这样出的图,在导入 Word 中的时候,字号大小会保持不变。而若是出图 10cm,导入 Word 后被压缩到 8.5cm,则图中 10 号字体也相应被压缩,已经不是原先的 10 号字体了,这也解释了为什么某些 SCI 中,每幅图的图中字体大小不一。

4 产品代码获取

后台回复关键词:科研绘图第二十二期三维瀑布图

  • 0
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
好的,我了解了您的问题。由于我是一个语言模型 AI,我没有办法在这里为您提供 MATLAB 三维绘图的代码和截图。但是,我可以为您提供一些基本的思路和步骤,帮助您解决这个问题。 首先,您需要准备一个信号数据,可以是自己生成的,也可以是从文件读取的。然后,使用 MATLAB 的 mesh、waterfall 和 view 函数来生成三维图形,以形象地表示信号的综合。 具体步骤如下: 1. 生成信号数据,可以使用 MATLAB 的 rand、sin、cos 等函数生成。 2. 使用 mesh 函数来生成信号三维图形。mesh 函数可以将二维矩阵数据绘制成三维网格图形,其矩阵的每个元素对应于三维坐标系的一个点。例如,可以使用以下代码来生成一个简单的正弦图形: ``` x = 0:0.1:10; y = 0:0.1:10; [X,Y] = meshgrid(x,y); Z = sin(X) + cos(Y); mesh(X,Y,Z); ``` 3. 使用 waterfall 函数来生成信号三维图形。waterfall 函数可以将二维矩阵数据绘制成三维瀑布图形,其矩阵的每个元素对应于三维坐标系的一个点,并且沿着一个维度会呈现出平滑的曲线。例如,可以使用以下代码来生成一个简单的正弦图形: ``` x = 0:0.1:10; y = 0:0.1:10; [X,Y] = meshgrid(x,y); Z = sin(X) + cos(Y); waterfall(X,Y,Z); ``` 4. 使用 view 函数来改变三维图形的视角。view 函数可以改变三维图形的方位角和仰角,使得观察者可以从不同的角度观察三维图形。例如,可以使用以下代码来改变三维图形的视角: ``` view(45,30); ``` 5. 观察吉布斯现象。吉布斯现象是指在傅里叶级数逼近一些特定函数时,会出现级数振荡的现象。这种现象可以用三维图形来直观地表示出来。例如,可以使用以下代码来生成一个方波信号三维图形,并观察吉布斯现象: ``` x = -pi:0.1:pi; y = 0:0.1:10; [X,Y] = meshgrid(x,y); Z = zeros(size(X)); for k = 1:100 Z = Z + (4/pi)*(1/(2*k-1))*sin((2*k-1)*X).*(exp(-0.1*(2*k-1)^2*Y)); end mesh(X,Y,Z); ``` 希望这些思路和步骤可以帮助您解决问题!
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

Lwcah(全网各平台账号同名)

您的鼓励是我创作的最大的动力~

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值