matlab语音信号时域特征分析代码

以下是MATLAB实现语音信号时域特征分析的示例代码：

% 读取语音信号
[s, Fs] = audioread('speech.wav');

% 时域特征分析
% 计算能量
energy = sum(s.^2);
disp(['能量为：', num2str(energy)]);

% 计算过零率
zcr = sum(abs(diff(sign(s)))) / (2 * length(s));
disp(['过零率为：', num2str(zcr)]);

% 计算短时能量
winlen = 0.02 * Fs; % 窗口长度
noverlap = 0.01 * Fs; % 窗口重叠长度
win = hamming(winlen); % 窗函数
[S, F, T] = spectrogram(s, win, noverlap, [], Fs); % 短时傅里叶变换
short_term_energy = sum(abs(S).^2); % 短时能量
figure;
plot(T, short_term_energy);
xlabel('时间（秒）');
ylabel('短时能量');

% 计算短时过零率
short_term_zcr = sum(abs(diff(sign(S)))) / (2 * size(S, 2));
figure;
plot(T, short_term_zcr);
xlabel('时间（秒）');
ylabel('短时过零率');

其中，语音信号文件`speech.wav`需要自己准备。代码中计算了语音信号的总能量、过零率以及短时能量、短时过零率，并绘制了短时能量、短时过零率随时间的变化曲线。