matlab 倒谱ceptrum,倒谱法求共振峰

我用自己的语音,读出来之后就 出现 Index exceeds matrix dimensions. 求解为什么。我的Fs=22050,我想请问共振峰的个数与什么有关

>> [x, fs, nbits]=wavread('C:\Users\WEICHANGCHENG\Desktop\H=15\无裂颖\1.wav');                  % 读入一帧数据

u=filter([1 -.99],1,x);                          % 预加重

wlen=length(u);                                  % 帧长

cepstL=6;                                        % 倒频率上窗函数的宽度

wlen2=wlen/2;

freq=(0:wlen2-1)*fs/wlen;                        % 计算频域的频率刻度

u2=u.*hamming(wlen);                                     % 信号加窗函数

U=fft(u2);                                       % 按式(9-2-1)计算

U_abs=log(abs(U(1:wlen2)));                      % 按式(9-2-2)计算

Cepst=ifft(U_abs);                               % 按式(9-2-3)计算

cepst=zeros(1,wlen2);

cepst(1:cepstL)=Cepst(1:cepstL);                 % 按式(9-2-5)计算

cepst(end-cepstL+2:end)=Cepst(end-cepstL+2:end);

spect=real(fft(cepst));                          % 按式(9-2-6)计算

[Loc,Val]=findpeaks(spect);                      % 寻找峰值

FRMNT=freq(Loc);                                 % 计算出共振峰频率

% 作图

pos = get(gcf,'Position');

set(gcf,'Position',[pos(1), pos(2)-100,pos(3),(pos(4)-140)]);

plot(freq,U_abs,'k');

hold on;  grid;

plot(freq,spect,'k','linewidth',2);

xlabel('频率/Hz'); ylabel('幅值/dB');

title('信号频谱,包络线和共振峰值')

fprintf('%5.2f   %5.2f   %5.2f   %5.2f\n',FRMNT);

for k=1 : 5

plot(freq(Loc(k)),Val(k),'kO','linewidth',2);

line([freq(Loc(k)) freq(Loc(k))],[-6 Val(k)],'color','k',...

'linestyle','-.','linewidth',2);

end

2283.50   4424.27   7171.60   Index exceeds matrix dimensions.

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值