VMD实现函数

变分模态分解函数，将一个信号分解成多个不同频段的信号，设置好参数即可将已有的信号分解。内置各参数的说明如1x1024，分解之后即得到5x1024的矩阵。每一行对应一个信号分量

clc;
clear all;
% load（‘VMD_basic.mat‘）
%%取出第一行与信号相减
load（‘simulation.mat‘）
x_trunc（1:1:2048） = x（1:1:2048）;
% self_sequence（1:3:512） = 1;
% self_sequence（2:3:512） = 2;
% self_sequence（3:3:512） = 3;
[SFNFT] = spectrogram（x_trunc hamming（128）1272561/139.975e6‘centered‘）;
figure
imagesc（TFabs（SFN））;

% some sample parameters for VMD
alpha = 1000;        % moderate bandwidth constraint
tau = 0;            % noise-tolerance （no strict fidelity enforcement）噪声耐受，一般取0就行
K = 5;              % 3 modes
DC = 0;             % no DC part imposed
init = 1;           % initialize omegas uniformly
tol = 1e-9;         %总信号拟合误差允许程度
[u_after_d u_hat omega] = VMD（x_trunc alpha tau K DC init tol）;
%--------------- Visualization

for k=1:K
[SFNFT] = spectrogram（u_after_d（k:） hamming（128）1272561/139.975e6‘centered‘）;
figure
mesh（TFabs（SFN））;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        923  2018-12-09 15:28  test_difference.m

     文件     126072  2018-12-06 09:59  simulation.mat

     文件       4851  2018-12-09 17:58  vmd_fun.m

----------- ---------  ---------- -----  ----

               131846                    3