用matlab编写VMD分解程序，用于分解时间序列，要求显示分解后的曲线图和中心频率

以下是一个简单的MATLAB程序，用于执行VMD分解并显示结果：

% 导入时间序列数据
data = load('time_series_data.mat');

% 设置分解参数
alpha = 2000;       % 正则化参数
tau = 0;            % 偏移参数
K = 10;             % 分解模式数量
DC = 1;             % 是否包括直流分量

% 执行VMD分解
[u, u_hat, omega] = VMD(data, alpha, tau, K, DC);

% 显示分解结果
figure;
subplot(K+1, 1, 1);
plot(data); title('Original Signal');
xlabel('Time'); ylabel('Amplitude');

for k=1:K
    subplot(K+1, 1, k+1);
    plot(u(k,:)); title(['Mode ', num2str(k), ', \omega = ', num2str(omega(k))]);
    xlabel('Time'); ylabel('Amplitude');
end

其中，`VMD.m`是VMD分解函数，可以从网上下载或自己编写。这个程序将显示原始时间序列和所有分解模式的曲线图，并且在每个分解模式的标题中显示其中心频率。