请给我写个matlab程序，是根据给出的浊音语音或清音语音，画出它分别加矩形窗和汉明窗的时域波形和短时频谱。

好的，以下是一个简单的 MATLAB 程序，用于绘制给定语音信号的时域波形和短时频谱，其中包括使用矩形窗和汉明窗两种窗口函数的绘制。

% 读取语音信号
[x, fs] = audioread('speech_signal.wav'); % 将文件名替换为您的语音文件名
% 绘制时域波形
t = linspace(0, length(x)/fs, length(x)); % 时间轴
figure;
subplot(2, 1, 1);
plot(t, x);
xlabel('Time (s)');
ylabel('Amplitude');
title('Waveform');

% 设置参数
win_len = 20; % 窗口长度（毫秒）
win_shift = 10; % 窗口移动步长（毫秒）
win_type = 'Hamming'; % 窗口类型（'Rect' 或 'Hamming'）
n_fft = 1024; % FFT 点数
n_overlap = n_fft - round(win_shift/1000*fs); % 重叠长度

% 分帧加窗
win = eval([win_type, '(win_len/1000*fs)']); % 定义窗口函数
x_frame = buffer(x, win_len/1000*fs, round(win_shift/1000*fs), 'nodelay');
x_frame = bsxfun(@times, x_frame, win);

% 计算短时傅里叶变换 (STFT)
[S, f, t] = stft(x, fs, 'Window', win, 'OverlapLength', n_overlap, 'FFTLength', n_fft);

% 绘制短时频谱
subplot(2, 1, 2);
imagesc(t, f, 20*log10(abs(S)));
axis xy;
xlabel('Time (s)');
ylabel('Frequency (Hz)');
title(['Short-time spectrum (', win_type, ' window)']);
colorbar;

您需要将程序中的 `speech_signal.wav` 替换为您要绘制的语音信号文件名，并根据需要调整其他参数，如窗口长度、移动步长、窗口类型和 FFT 点数等。程序将绘制语音信号的时域波形和短时频谱，包括使用矩形窗和汉明窗两种窗口函数的绘制。