根据给出的浊音语音，写出它分别加矩形窗和汉明窗的时域波形和短时频谱的matlab代码。 2、根据给出的清音语音，写出它分别加矩形窗和汉明窗的时域波形和短时频谱matlab代码

好的，以下是针对浊音语音和清音语音分别加矩形窗和汉明窗的时域波形和短时频谱的MATLAB代码：

浊音语音：

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

% 设置窗口大小及重叠长度
winlen = 0.025; % 窗口大小
winshift = 0.01; % 重叠长度

% 计算窗口长度和重叠长度的样点数
nwinlen = round(winlen*fs); % 窗口长度（样点数）
nwinshift = round(winshift*fs); % 重叠长度（样点数）

% 分帧加窗
nx = length(x);
nframe = floor((nx-nwinlen)/nwinshift)+1; % 帧数
w = rectwin(nwinlen); % 矩形窗
win = repmat(w, 1, nframe); % 窗矩阵
x_seg_rect = x(1:nwinlen+(nframe-1)*nwinshift) .* win; % 分帧加窗（矩形窗）

w = hamming(nwinlen); % 汉明窗
win = repmat(w, 1, nframe); % 窗矩阵
x_seg_ham = x(1:nwinlen+(nframe-1)*nwinshift) .* win; % 分帧加窗（汉明窗）

% 时域波形
t = (0:nx-1)/fs; % 时间轴
subplot(2,2,1);
plot(t,x,'k');
title('原始波形');
subplot(2,2,2);
plot(t(1:nwinlen+(nframe-1)*nwinshift),x_seg_rect,'k');
title('加矩形窗后波形');
subplot(2,2,3);
plot(t(1:nwinlen+(nframe-1)*nwinshift),x_seg_ham,'k');
title('加汉明窗后波形');

% 短时频谱
nfft = 512; % FFT点数
f = (0:nfft/2-1)*fs/nfft; % 频率轴
subplot(2,2,4);
spectrogram(x_seg_rect,winlen*fs,winshift*fs,nfft,fs,'yaxis');
title('加矩形窗后短时频谱');
subplot(2,2,5);
spectrogram(x_seg_ham,winlen*fs,winshift*fs,nfft,fs,'yaxis');
title('加汉明窗后短时频谱');

清音语音：

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

% 设置窗口大小及重叠长度
winlen = 0.025; % 窗口大小
winshift = 0.01; % 重叠长度

% 计算窗口长度和重叠长度的样点数
nwinlen = round(winlen*fs); % 窗口长度（样点数）
nwinshift = round(winshift*fs); % 重叠长度（样点数）

% 分帧加窗
nx = length(x);
nframe = floor((nx-nwinlen)/nwinshift)+1; % 帧数
w = rectwin(nwinlen); % 矩形窗
win = repmat(w, 1, nframe); % 窗矩阵
x_seg_rect = x(1:nwinlen+(nframe-1)*nwinshift) .* win; % 分帧加窗（矩形窗）

w = hamming(nwinlen); % 汉明窗
win = repmat(w, 1, nframe); % 窗矩阵
x_seg_ham = x(1:nwinlen+(nframe-1)*nwinshift) .* win; % 分帧加窗（汉明窗）

% 时域波形
t = (0:nx-1)/fs; % 时间轴
subplot(2,2,1);
plot(t,x,'k');
title('原始波形');
subplot(2,2,2);
plot(t(1:nwinlen+(nframe-1)*nwinshift),x_seg_rect,'k');
title('加矩形窗后波形');
subplot(2,2,3);
plot(t(1:nwinlen+(nframe-1)*nwinshift),x_seg_ham,'k');
title('加汉明窗后波形');

% 短时频谱
nfft = 512; % FFT点数
f = (0:nfft/2-1)*fs/nfft; % 频率轴
subplot(2,2,4);
spectrogram(x_seg_rect,winlen*fs,winshift*fs,nfft,fs,'yaxis');
title('加矩形窗后短时频谱');
subplot(2,2,5);
spectrogram(x_seg_ham,winlen*fs,winshift*fs,nfft,fs,'yaxis');
title('加汉明窗后短时频谱');

需要注意的是，上述代码中的`voice.wav`是需要替换为实际的语音文件名。另外，窗口大小和重叠长度的设置需要根据实际情况进行调整。