matlab代码，输入一段人的语音信号，通过STFT，可以得到音频信号的短时频谱图，使用STFT分析自己说话时的音调变化分析信号在不同时间段内的频率分布

以下是使用MATLAB进行STFT分析语音信号的示例代码，可以得到不同时间段内的频率分布：

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

% 设置STFT参数
win = hamming(256); % 窗函数
hop = length(win)/2; % 帧移
nfft = 512; % FFT点数

% 计算STFT
[S, F, T] = stft(x, Fs, 'Window', win, 'OverlapLength', hop, 'FFTLength', nfft);

% 绘制频谱图
imagesc(T, F, abs(S));
set(gca, 'YDir', 'normal');
xlabel('时间 (s)');
ylabel('频率 (Hz)');
colorbar;

% 提取特定时间段的频率分布
t_start = 1; % 起始时间（秒）
t_end = 2; % 结束时间（秒）

% 将时间转换为STFT中的帧数
frame_start = ceil(t_start/hop);
frame_end = floor(t_end/hop);

% 提取频率分布数据
freq_dist = mean(abs(S(:, frame_start:frame_end)), 2);

% 绘制频率分布图
figure;
plot(F, freq_dist);
xlabel('频率 (Hz)');
ylabel('幅度');

在上面的代码中，我们首先使用`audioread`函数读取输入的语音信号。然后，我们设置STFT参数，包括窗函数、帧移和FFT点数。接下来，我们使用`stft`函数计算语音信号的STFT，并将结果存储在`S`中。最后，我们使用`imagesc`函数绘制短时频谱图，其中`T`表示时间轴，`F`表示频率轴，`abs(S)`表示STFT的幅度。`set(gca, 'YDir', 'normal')`用于将频率轴反向，使其与常见的频谱图相同。`xlabel`和`ylabel`用于添加轴标签，`colorbar`用于添加颜色条。

接下来，我们提取特定时间段内的频率分布数据。这里我们假设要提取1到2秒的时间段，首先将这段时间转换为STFT中的帧数。然后，我们使用`mean`函数计算在这段时间内每个频率的平均幅度，并将结果存储在`freq_dist`中。

最后，我们使用`plot`函数绘制频率分布图，其中`F`表示频率轴，`freq_dist`表示频率分布幅度。`xlabel`和`ylabel`用于添加轴标签。