小波包信号能量特征提取

clear all;
close all;
clc;
% 给定参数
fs=1000000;
N=2048;
n=0:N-1;
t=n/fs;
number=fix（log2（N））-1;
x=xlsread（‘870.xlsx‘）;
data=x（:1）;
%绘制实测信号时域波形
figure（1）;
plot（tdata）;xlabel（‘时间 t/s‘）;ylabel（‘幅值/mV‘）;
set（gca‘xlim‘[0 2.047*10^-3]）; 
% some sample parameters for VMD
alpha = 2000;       % moderate bandwidth constraint
tau = 0;            % noise-tolerance （no strict fidelity enforcement）
K = 8;              % 5 modes
DC = 0;             % no DC part imposed
init = 1;           % initialize omegas uniformly
tol = 1e-7;
%--------------- Run actual VMD code
figure（2）;
[u u_hat omega] = VMD（dataalpha tau K DC init tol）;
subplot（size（u1）+121）;
plot（u（1:）‘k‘）;grid on;
title（‘VMD分解‘）;
ylabel（{‘imf1‘;‘幅值/mV‘}）
set（gca‘xlim‘[0 2048]）;
subplot（size（u1）+122）;
plot（abs（fft（u（1:）））‘k‘）;
title（‘对应频谱‘）;
for i = 2:size（u1）+1
    subplot（size（u1）+12i*2-1）;
    plot（u（i-1:）‘k‘）;grid on;
    ylabel（{‘imfi‘;‘幅值/mV‘}）
set（gca‘xlim‘[0 2