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.

matlab 倒谱ceptrum,倒谱法求共振峰相关推荐

matlab 计算变异系数,变异系数法求权重matlab 代码
利用matlab编程,很容易根据变异系数法,求得多指标综合评价模型的权重. 代码如果有不懂的地方,可以联系我. 变异系数法求权重matlab 代码 clear;clc; [data1,header1] ...
matlab 计算变异系数,变异系数法求权重matlab代码
<变异系数法求权重matlab代码>由会员分享,可在线阅读,更多相关<变异系数法求权重matlab代码(1页珍藏版)>请在读根文库上搜索. 1.变异系数法求权重 matlab ...
用matlab编程节点电压法求电路,MATLAB在电路中的应用
<MATLAB在电路中的应用>由会员分享,可在线阅读,更多相关<MATLAB在电路中的应用(59页珍藏版)>请在人人文库网上搜索. 1.MATLAB应用(三) Matlab在电 ...
matlab e 精确到,matlab中用0.618法求minf(x)=e^(-x)+x^2在区间(0,1)上的极小值,精确到0.03....
共回答了15个问题采纳率:86.7% clc clear all; elp=0.03; tao=0.618; N=fix(log(elp)/log(tao))+1; k=1; a(k)=0; b(k) ...
matlab怎么输入积分公式,在matlab下用梯形法求函数的积分
函数是网上看到的,加上了一点我自己的理解 %用来就数值积分 %fx是由syms定义的函数表达式 function Trapezia(a,b,fx,E,Nfprintf('\n***********st ...
深度学习和语音系列教程 5-100：求共振峰
文章目录 1 基音频率 1.1 基音频率的查看 2 共振峰的获取 2.1 倒谱法求共振峰 2.2 测试结果 2.3 共振峰位置坐标 2.4 多种情况下的共振峰测试 1 基音频率 1.1 基音频率的查看 ...
Matlab实现倒谱法求基音频率和共振峰
文章目录前言一.倒谱二.基音周期 1.流程图 2. 实现代码(Matlab) 三.共振峰四.实验总结前言有关同态.倒谱.基音周期等概念,可参考一篇本科毕业论文,链接:link 一.倒谱 ...
基于倒谱法和线性预测法估计基音频率（MATLAB和Python）
基于倒谱法和线性预测法估计基音频率(MATLAB和Python) 倒谱法基音检测在python中实现一帧信号的基音频率估计 wlen = 256 inc = 128 pitch = [] x1, F ...
用matlab的编程法和游动鼠标法求二阶传递函数的上升时间、峰值时间、超调量和调节时间 - Gavin_Hall的博客 - CSDN博客
1. 准备终值:c(∞) 上升时间 tr:响应从峰值的10%上升到峰值的90%所需要的时间:而阶跃响应则是从终值的10%上升到终值的90%所需要的时间:对有振荡的系统,也可以定义为从0到第一次到达终 ...
matlab实现正割法求根和Aitken加速方法
正割法求根代码: zhengge.m function zhengge(x0,x1,e)syms x;f(x)=exp(-x);K=0;while abs(x1-x0)>efprintf('第% ...

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

matlab 倒谱ceptrum,倒谱法求共振峰相关推荐

最新文章

热门文章