MATLAB的语音信号采集与处理分析
1、基本描述
本文描述的系统是一个全面而精细的语音信号处理平台，核心组件由MATLAB的高级功能模块构建而成。系统的核心交互界面，借助于MATLAB的uifigure函数搭建，为用户提供了一个直观且响应迅速的操作环境。通过设计的GUI按钮，如“录音”、“FFT频谱图”、“叠加噪声”、“巴特沃斯滤波器”、“小波分解去噪”及“播放录音”，用户可以便捷地触发各项信号处理功能，实现从信号采集到分析，再到噪声控制与滤波的操作。
在信号处理中，设计的这个系统采用了当前先进的算法和技术。录音模块使用的是MATLAB的audiorecorder，以保证高质量的音频信号实时采集。FFT频谱分析模块则通过快速傅立叶变换，提供了对信号频域特性的可视化处理。噪声叠加模块是采用了引入随机噪声的方式，以模拟真实环境中的噪声干扰。滤波部分使用的是巴特沃斯滤波器，有效去除信号中的噪声成分，提升信号的清晰度。最后是使用小波变换和阈值处理技术，对信号进一步的优化处理。
为了保证系统的实用性和灵活性，数据存储与读取功能使用了audiowrite和audioread函数，支持信号的保存与加载，方便用户在不同会话间共享和复用数据。最后，音频播放功能通过sound函数，允许用户随时监听原始信号及处理后的结果，确保处理效果符合预期，同时也为系统的人机交互体验增添了直观的听觉反馈。
语音系统的GUI部分设计遵循简洁直观的原则，采用MATLAB内置的UI元素，包括按钮和绘图区域。绘图区域用于显示信号的时域波形和频域特性，可以更加直观理解信号的变化。Matlab生成的可视化操作界面如下所示。
2、功能效果展示
2.1、FFT频谱
通过对采集到的语音信号执行快速傅里叶变换（FFT），成功地获取了信号的频谱图，如下图所示。频谱图清晰地展示了信号在不同频率上的能量分布情况。从频谱图中可以看出，语音信号主要集中在较低的频率范围内，这与人类语音的声学特性相符。
2.2、噪声叠加处理
语音信号经过噪声叠加处理后，在语音信号的波形图上明显出现了噪声的叠加效果。通过观察如下图所示的实验结果图，可以发现信号波形中增加了随机波动，这些波动代表了叠加的噪声成分。
2.3、巴特沃斯滤波处理
如下图所示的实验结果图，原本含有噪声的语音信号经过巴特沃斯滤波处理后，其波形图显示出了明显的改善。噪声成分得到了有效的抑制，信号的清晰度显著提升。通过巴特沃斯滤波器，成功地降低了语音信号中的噪声水平，提高了信号质量。
2.4、小波分解处理
通过下图所示的波形图，可以观察到去噪后的信号保留了原始信号的主要特征，同时消除了大部分噪声。这种效果在频谱图上也有所体现，其中高频噪声成分被显著抑制，而信号的主体频谱得以保留。
3、程序源码
在运行程序代码时，所需的模组组件必须安装并正常运行，否则运行会出现问题。
function varargout = yuyinshuzhixinhao(varargin)
% Initialize GUI
fig = uifigure('Name', '语音信号处理');
% Create buttons
btnRecord = uibutton(fig, 'push', 'Text', '录音', ...
'Position', [150, 300, 120, 30], ...
'ButtonPushedFcn', @(btnRecord, event) recordButtonCallback());
btnFFT = uibutton(fig, 'push', 'Text', 'FFT频谱图', ...
'Position', [300, 300, 120, 30], ...
'ButtonPushedFcn', @(btnFFT, event) fftButtonCallback());
btnNoise = uibutton(fig, 'push', 'Text', '叠加噪声', ...
'Position', [450, 300, 120, 30], ...
'ButtonPushedFcn', @(btnNoise, event) noiseButtonCallback());
btnFilter = uibutton(fig, 'push', 'Text', '巴特沃斯滤波器', ...
'Position', [600, 300, 120, 30], ...
'ButtonPushedFcn', @(btnFilter, event) filterButtonCallback());
btnWavelet = uibutton(fig, 'push', 'Text', '小波分解去噪', ...
'Position', [750, 300, 120, 30], ...
'ButtonPushedFcn', @(btnWavelet, event) waveletButtonCallback());
btnPlay = uibutton(fig, 'push', 'Text', '播放录音', ...
'Position', [900, 300, 120, 30], ...
'ButtonPushedFcn', @(btnPlay, event) playButtonCallback());
% Create axes for plotting
ax1 = uiaxes(fig, 'Position', [50, 50, 400, 200]);
ax2 = uiaxes(fig, 'Position', [500, 50, 400, 200]);
% Callback functions
function recordButtonCallback()
recObj = audiorecorder;
disp('Start speaking.');
recordblocking(recObj, 5);
disp('End of Recording.');
play(recObj);
myRecording = getaudiodata(recObj);
plot(ax1, myRecording);
title(ax1, '录音波形');
xlabel(ax1, '时间');
ylabel(ax1, '幅度');
% Save the recording to a file
filename = 'data1.wav';
audiowrite(filename, myRecording, 8000);
end
function fftButtonCallback()
filename = 'data1.wav';
[x, fs] = audioread(filename);
y1 = fft(x, 1024);
f = fs * (0:511) / 1024;
% Plot FFT spectrum
plot(ax2, f, abs(y1(1:512)));
title(ax2, 'FFT频谱');
xlabel(ax2, '频率/Hz');
ylabel(ax2, '幅度');
end
function noiseButtonCallback()
filename = 'data1.wav';
[x, fs] = audioread(filename);
x1 = 0.1 * randn(length(x), 2);
x2 = x1 + x;
sound(x2, fs);
y1 = fft(x, 1024);
y2 = fft(x2, 1024);
f = fs * (0:511) / 1024;
% Plot the original and noisy signals
plot(ax1, x);
title(ax1, '原始语音信号时域图');
xlabel(ax1, '时间');
ylabel(ax1, '幅度');
plot(ax2, abs(y1(1:512)));
title(ax2, '原始信号频谱');
xlabel(ax2, '频率/Hz');
ylabel(ax2, '幅度');
% Plot noisy signal in the second subplot
hold(ax2, 'on');
plot(ax2, abs(y2(1:512)), 'r');
hold(ax2, 'off');
legend(ax2, '原始信号', '加噪信号');
end
function filterButtonCallback()
filename = 'data1.wav';
[x, fs] = audioread(filename);
wp = 0.1 * pi;
ws = 0.4 * pi;
Rp = 1;
Rs = 15;
[N, Wn] = ellipord(wp / pi, ws / pi, Rp, Rs);
[b, a] = ellip(N, Rp, Rs, Wn);
f1 = filter(b, a, x);
sound(f1, fs);
y1 = fft(x, 1024);
y2 = fft(f1, 1024);
f = fs * (0:511) / 1024;
% Plot the filtered signal
plot(ax1, x);
title(ax1, '滤波前的时域波形');
xlabel(ax1, '时间');
ylabel(ax1, '幅度');
plot(ax2, abs(y1(1:512)));
title(ax2, '滤波前的频谱');
xlabel(ax2, '频率/Hz');
ylabel(ax2, '幅度');
% Plot filtered signal in the second subplot
hold(ax2, 'on');
plot(ax2, abs(y2(1:512)), 'r');
hold(ax2, 'off');
legend(ax2, '滤波前', '滤波后');
end
function waveletButtonCallback()
if license('test', 'Wavelet_Toolbox')
try
filename = 'data1.wav';
[x, fs] = audioread(filename);
% 小波分解
n = 10; % 分解层数，可以根据实际情况调整
wname = 'db2'; % 小波基函数，可以尝试不同的基函数
[C, L] = wavedec(x, n, wname);
% 估计噪声标准差
sigma = wnoisest(C, L, 1);
% 计算阈值
alpha = 2; % 阈值调节因子
thr = wbmpen(C, L, sigma, alpha);
keepapp = 1; % 是否保留近似系数
% 软阈值去噪
yb = wdencmp('gbl', C, L, wname, n, thr, 's', keepapp);
plot(ax1, x);
title(ax1, '原始语言信号');
xlabel(ax1, '时间');
ylabel(ax1, '幅度');
plot(ax2, yb);
title(ax2, '去噪后的语言信号');
xlabel(ax2, '时间');
ylabel(ax2, '幅度');
catch ME
uialert(fig, sprintf('An error occurred: %s', ME.message), 'Error');
end
% 在这里添加播放去噪后信号的代码
sound(yb, fs); % `yb` 是去噪后的信号，`fs` 是采样率
else
uialert(fig, 'Wavelet Toolbox is not installed.', 'Toolbox Error');
end
end
function playButtonCallback()
filename = 'data1.wav';
[myRecording, fs] = audioread(filename);
sound(myRecording, fs);
end
% Show the figure
fig.Visible = 'on';
end