%% 利用RAPT算法估计瞬时频率信息 frame_len = fix(0.03fs); %帧长 frame_shift = fix(0.015fs); %帧移 pre_emph = [1 0.99]; %预加重滤波器系数 y = filter(1,pre_emph,y); %预加重 frame_num = fix((length(y)-frame_len)/frame_shift)+1; %帧数 w = hamming(frame_len); %加窗窗口 f0 = zeros(frame_num,1); %瞬时频率 for i = 1:frame_num frame = y((i-1)*frame_shift+1:(i-1)*frame_shift+frame_len); %取出一帧语音信号 frame = frame.*w; %加窗 [f0(i),~] = rapt(frame,fs); %使用RAPT算法估计瞬时频率 end 为什么选用rapt算法来估计瞬时频率信息

RAPT算法是一种比较常用的瞬时频率估计算法，对于周期信号和非周期信号都有一定的适用性，并且具有较高的估计精度。在语音信号分析和处理中，瞬时频率是非常重要的参数，可以用于声学特征提取、语音识别、语音合成等应用中。因此，RAPT算法被广泛应用于语音信号处理领域，特别是在瞬时频率估计方面，被认为是一种比较优秀的算法。