变分模态分解python代码_VMD变分模态分解代码问题

clear all;

close all;

clc;

% Time Domain 0 to T

T = 1000;

fs = 1/T;

t = (1:T)/T;

freqs = 2*pi*(t-0.5-1/T)/(fs);

% center frequencies of components

f_1 = 2;

f_2 = 24;

f_3 = 288;

% modes

v_1 = (cos(2*pi*f_1*t));          %信号1

v_2 = 1/4*(cos(2*pi*f_2*t));     %信号2

v_3 = 1/16*(cos(2*pi*f_3*t));    %信号3

% for visualization purposes

fsub = {};

wsub = {};

fsub{1} = v_1;

fsub{2} = v_2;

fsub{3} = v_3;

wsub{1} = 2*pi*f_1;

wsub{2} = 2*pi*f_2;

wsub{3} = 2*pi*f_3;

% composite signal, including noise

f = v_1 + v_2 + v_3 +