请利用MATLAB软件写出振动信号通过小波包分析后的程序，包括去噪信号图、信号频谱图、信号能量谱图

下面是振动信号通过小波包分析后的MATLAB程序：

% 读取振动信号数据
data = load('vibration_signal.txt');

% 小波包分析
wpt = wpdec(data, 4, 'db4');

% 去噪处理
threshold = thselect(data, 'rigrsure');
denoised_data = wprden(wpt, 'sqtwolog', threshold);

% 信号频谱图
figure;
f = 0:0.1:500;
P = pwelch(denoised_data,[],[],f);
plot(f,P);

% 信号能量谱图
figure;
[~,~,E] = wenergy(wpt);
plot(0:length(E)-1, E);

% 原始信号与去噪信号图
figure;
plot(data);
hold on;
plot(denoised_data);
legend('原始信号', '去噪信号');

说明：

1. 程序首先读取振动信号数据，保存在 `data` 变量中。
2. 使用 `wpdec` 函数对振动信号进行小波包分析，分解层数为 4，小波函数为 `db4`。
3. 使用 `thselect` 函数选择适当的阈值进行去噪处理。
4. 使用 `wprden` 函数对小波包系数进行阈值处理，得到去噪后的信号数据 `denoised_data`。
5. 使用 `pwelch` 函数计算去噪信号的频谱图，并绘制出来。
6. 使用 `wenergy` 函数计算小波包系数的能量谱图，并绘制出来。
7. 最后，绘制原始信号与去噪信号的图像，方便比较两者的差异。