Different kind of 3D waterfall for Different Applications

3D Waterfall for EMD

Code

clear all;clc;close all;

set(0, 'DefaultAxesFontSize', 16); % 设置坐标轴标签的字体大小

set(0, 'DefaultTextFontSize', 16); % 设置文本框的字体大小

% 假设你已经有了一个信号，这里我们创建一个示例信号

t = 0:0.01:10; % 时间向量

signal = cos(2*pi*t) + 0.5*cos(4*pi*t) + rand(1, length(t)) * 0.1; % 示例信号

% 使用EMD进行分解

[IMFs, residuals] = emd(signal);

% waterfall plot the IMFs

figure('Position', [0 0 1500 900]);

plot3(zeros(length(IMFs), 1), t, signal, 'LineWidth', 1.5, 'DisplayName', 'Signal');

hold on;

for i = 1:size(IMFs, 2)

plot3(ones(length(IMFs), 1) * i, t, IMFs(:, i), 'LineWidth', 1.5, 'DisplayName', ['IMF ' num2str(i)]);

hold on;

end

% plot the residuals

plot3(ones(length(IMFs), 1) * (size(IMFs, 2) + 1), t, residuals, 'LineWidth', 1.5, 'DisplayName', 'Residuals');

hold off;

xStringLables = cell(1, size(IMFs, 2) + 2);

xStringLables{1} = 'Signal';

for i = 1:size(IMFs, 2)

xStringLables{i + 1} = ['IMF ' num2str(i)];

end

xStringLables{size(IMFs, 2) + 2} = 'Residuals';

% 设置图例

legend('show');

xlabel('IMF Index');

ylabel('Time');

zlabel('Amplitude');

title('EMD Decomposition of Signal');

xticks(0:size(IMFs, 2) + 1);

xticklabels(xStringLables);

grid on;

% 保存图像

saveas(gcf, 'EMD_Decomposition.png');

1.6 KiB Raw Blame History

Different kind of 3D waterfall for Different Applications

3D Waterfall for EMD

Code

Result

1.6 KiB

Raw Blame History