matlab写一段振动数据时频特征提取

以下是一段MATLAB代码，用于提取振动数据的时频特征：

% 读取振动数据（假设为列向量 y）
Fs = 1000; % 采样率为1000 Hz
t = (0:length(y)-1)/Fs; % 时间轴
N = length(y); % 信号长度

% 使用小波包变换进行时频分析
wname = 'db4'; % 小波基函数
wpt = wpdec(y,5,wname); % 5层小波包分解
[spec,t,f] = wpfcoef(wpt,t); % 计算小波包系数的时频分布
P = abs(spec).^2; % 幅值平方，代表能量谱密度

% 提取时频特征
mean_freq = sum(f.*P)/sum(P); % 均值频率
[max_amp,max_idx] = max(P(:)); % 最大能量位置
[max_row,max_col] = ind2sub(size(P),max_idx); % 转换为行列坐标
max_freq = f(max_row); % 最大能量对应的频率
max_time = t(max_col); % 最大能量对应的时间

% 可视化时频分布
figure;
surf(t,f,P,'EdgeColor','none');
view(2);
axis tight;
colormap(jet);
xlabel('Time (s)');
ylabel('Frequency (Hz)');
title('Wavelet Packet Spectrum');

该代码首先读取振动数据，并使用小波包变换进行时频分析。然后根据小波包系数的时频分布，计算能量谱密度并提取时频特征，包括均值频率、最大能量位置、最大能量对应的频率和时间。最后，使用surf函数可视化小波包谱。