基于小波变换的语音增强的matlab代码显示语音波形和语谱图

以下是基于小波变换的语音增强的Matlab代码示例，同时显示语音波形和语谱图：

% 加载语音信号
[x, fs] = audioread('test.wav');

% 显示原始语音波形
subplot(2, 2, 1);
plot(x);
title('原始语音波形');

% 显示原始语音语谱图
subplot(2, 2, 2);
spectrogram(x, 256, 250, [], fs, 'yaxis');
title('原始语音语谱图');

% 进行小波变换
[c, l] = wavedec(x, 5, 'db4');

% 计算小波系数的能量分布情况
energy = cumsum(c.^2);

% 选取适当的阈值进行去噪处理
threshold = 0.1 * energy(end);
c(energy < threshold) = 0;

% 进行小波反变换
y = waverec(c, l, 'db4');

% 显示增强后的语音波形
subplot(2, 2, 3);
plot(y);
title('增强后的语音波形');

% 显示增强后的语音语谱图
subplot(2, 2, 4);
spectrogram(y, 256, 250, [], fs, 'yaxis');
title('增强后的语音语谱图');

运行上述代码后，会显示出原始语音的波形和语谱图，以及经过小波变换去噪后的增强语音的波形和语谱图。可以通过比较两者的差异，来评估小波变换的语音增强效果。