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

末蓝、 2023-01-16 06:21 129阅读 0赞

我用自己的语音，读出来之后就 出现 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 条评论，129人围观）

还没有评论，来说两句吧...

相关阅读

相关基于MATLAB的倒谱分析与MFCC系数计算

基于MATLAB的倒谱分析与MFCC系数计算倒谱分析和Mel频率倒谱系数（MFCC）是语音信号处理中常用的特征提取方法。在本文中，我们将使用MATLAB来实现倒谱分析和MF

àì夳堔傛蜴生んèń/ 2024年03月03日 08:14/ 0 赞/ 53 阅读

相关音频-特征提取：①幅度谱（短时傅里叶变换谱/STFT）、②梅尔频谱（mel-spectrogram）、③梅尔倒谱（MFCC）【在梅尔频谱上取对数，做DCT（离散余弦变换）变换，得梅尔倒谱】

![在这里插入图片描述][0319ff6c266c4cd8a180049ab3beacf2.png_pic_center] 原始信号从音频文件中读取出来的原始语音信号

心已赠人/ 2023年10月05日 23:02/ 0 赞/ 1 阅读

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

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

末蓝、/ 2023年01月16日 06:21/ 0 赞/ 130 阅读

相关 matlab求雷克子波相位谱,求雷克子波的振幅谱和相位谱的MATLAB程序，谢谢！！！！！...

%% 自己写了如下程序，望采纳！ clc;clear;close all; %Ricker子波 f=40;dt=0.004; nw=6./f/dt; nw=2\flo

布满荆棘的人生/ 2023年01月15日 08:23/ 0 赞/ 167 阅读

相关功率谱估计 matlab arma,MATLAB/SIMULINK 功率谱估计

丁老师，您好！我测量了脉动速度，采样频率为4000Hz，测量时间为4.09575s，测量点数16384，想用FFT对其进行功率谱分析，我参考同学的程序(其采样频率10Hz，测量

淩亂°似流年/ 2023年01月15日 08:06/ 0 赞/ 146 阅读

相关核磁谱图分析步骤_【技能】X射线光电子能谱（XPS）谱图分析

![f568a3b6d47d490622de932efd647f20.png][] 审核：小铅笔来源：材料基一、X光电子能谱分析的基本原理 X光电子能谱分析的基本原理

我不是女神ヾ/ 2023年01月05日 12:03/ 0 赞/ 142 阅读

相关 geomagic 色谱分析

一、打开软件，导入模型导入模型，通常是一个CAD模型与一个点云模型，选中CAD模型右键设置为“参照”，将点云选中右键设置为“测试”，但是如果是两个点云模

向右看齐/ 2022年12月12日 15:16/ 0 赞/ 155 阅读

相关 python绘出语音信号的倒谱图

python画语音信号的倒谱图求信号x(n)的倒谱： 1）对信号x(n)作傅里叶变换，得到信号的复数谱 X ( w ) = F T \[ x ( n )

绝地灬酷狼/ 2022年09月02日 15:53/ 0 赞/ 272 阅读

相关理解谱聚类

理解谱聚类前面介绍过K-means聚类方法，这个方法简单易懂，主要在于如何定义距离计算公式（一般使用欧氏距离），如何选择K值，这两个问题。这次我们介绍谱聚类，它是K-me

你的名字/ 2022年08月05日 21:28/ 0 赞/ 211 阅读

相关 Matlab实现(KNN)自适应谱聚类

自适应谱聚类谱聚类中常用构造K近邻关系的相似度矩阵，这里引用实现相关算法，首先根据输入的数据data和K值求解KNN矩阵：get\_knn\_distance.m

刺骨的言语ヽ痛彻心扉/ 2022年06月17日 02:30/ 0 赞/ 295 阅读