yinpinboxing.txt

来自「语音信号分析与处理及其MATLAB实现,语音信号处理」· 文本 代码 · 共 27 行

"Matlab 音频波形处理"

[y,fs] = wavread('.wav')；  %读一个已保存的WAV 语音文件
wavplay（y)；               %播放语音文件
t=(0：length(y)-1)/fs；     %计算语音播放时间
subplot(3，1，1)；          %确定语音波形的显示位置
plot(t，y)；                %画波形图
legend('波形图')；
xlabel('时间(s)')；         %X轴的标题
ylabel('幅度')；            %Y轴的标题
x=fft(y.*hamming(length(y)))； %加hamming 窗快速傅立叶变换
fm=5000*length(x)/fs；         %设置频率范围
f=(0：fm)*fs/length(x)；       %确定频率刻度
subplot(3，1，2)；             %确定语音波形的显示位置
plot(f，20*log10(abs(x(1：length(f)))+eps))； %画频谱图
legend('频谱图')；
xlabel('频率（Hz）')；
ylabel('频谱幅度（dB）')；
c=fft(log(abs(x)+eps))；       % 倒频谱计算
ms1=fs/1000；
ms20=fs/50；
q=(ms1：ms20)/fs；             %确定倒频刻度
subplot(3，1，3)；
plot(q，abs(c(ms1：ms20)))；   %画倒谱图
legend('倒谱图')；
xlabel('倒频（ｓ）')；
ylabel('倒频谱幅度')；