Matlab计算语音特征(如RMS)

Question

我需要在matlab中计算一些语音特征。但是，我不知道在从wavread函数获得数据之后需要执行哪些操作，然后才能使用功能公式(例如RMS公式)。

我刚发现一个网站：

rmsVec(m) = sqrt(sum(wavData{m}(:).^2)/length(wavData{m}(:)));

但在另一种情况下，他们使用DFT，然后用实的和想象的矩阵/向量的部分运算，我不明白。

那应该是怎样的呢？

我有音箱库，但我没有在那里找到计算rms。

我的matlab版本: R2009b

编辑：

http://www.edaboard.com/thread163664.html见妖精和petejonze答案

因为石化，不管我是用时域还是频域，但是

[y, fs, aa] = wavread('C:\Users\lukasz\Desktop\semestr 7\inzynierka\dzwieki\uczace\dol\dol_50')
y1 = y(:,1)
Y1 = fft(y1)
rms1 = sqrt(mean(y1.^2))
rms2 = sqrt(sum(Y1.*conj(Y1))/size(Y1,1))

其结果是：

rms1 =

0.0577

rms2 =

13.2706

所以我有点困惑。哪一个是对的？

编辑2：

另一个例子是：

>> y = [1;2;3]

y =

 1
 2
 3

>> Y = fft(y)

Y =

6.0000          
-1.5000 + 0.8660i
-1.5000 - 0.8660i

>> rms1 = sqrt(mean(y.^2))

rms1 =

2.1602

>> rms2 = sqrt(sum(Y.*conj(Y))/size(Y,1))

rms2 =

3.7417

>> size(Y,1)

ans =

 3

>> Y.*conj(Y)

ans =

36
 3
 3

所以这个公式

sqrt(sum(Y1.*conj(Y1))/size(Y1,1))

是不正确的，应该是

sqrt(sum(Y1.*conj(Y1))/size(Y1,1)^2)

Answer 1

sqrt(sum(wavData{m}(:).^2)/length(wavData{m}(:)))；= rmsVec(m)

这是正确的。

但在另一种情况下，他们使用DFT，然后用实的和想象的矩阵/向量的部分运算，我不明白。

你可以问它，但你需要提供更多的细节。你说的是什么其他网站，你到底不明白什么。

我有音箱库，但我没有在那里找到计算rms。

如果你使用melcepst

http://www.ee.ic.ac.uk/hp/staff/dmb/voicebox/doc/voicebox/melcepst.html

您可以指定E0dD来包含能量。

          (2) c=melcepst(s,fs,'E0dD')   % include log energy, 0th cepstral coef, delta and delta-delta coefs

您可以对代码进行分析，看看在引擎盖下面发生了什么。

rms1 =sqrt(平均(y1.^2))

这是正确的

rms2 = sqrt(sum(Y1.*conj(Y1))/size(Y1,1))

这是错误的，它必须是sqrt (和(Y1.conj(Y1)/size(Y1，1)，其中除法在sqrt之后。因为matlab中的FFT函数不以1/sqrt(n)为尺度，因而不对称。所以sum (y^2)=sum(y*conj(Y))/n，您需要在sqrt之后除以。