1
/
5

非负矩阵分解音乐数据


大家好,我是EYS的张森,在EYS从事音乐机器学习的开发业务。今天我们谈谈非负矩阵分解(NMF)。

矩阵分解和非负矩阵分解


从某段钢琴的触的第一成分


第一成分强度的时间推移

矩阵分解,就是把一个(N,M)的矩阵V分解成两个矩阵W和H

V = WH

W和H分别为(N,h)和(h, M)的矩阵。这样一来,乍一看好像需要处理的变数增加了,然而通常情况下h远远小于N和M,所以整体的成分是减少的。也就是说用更少的数字记述同一个量。

在这里,V相当于你拥有N个样本,每个样本里有M种数据。而把V分解成W和H,就相当于你把原先有M成分的量分解成h种主要的成分(H),然后把每个样本数据用这h种成分来说明(W)。这样做的原因是,M在数量上可能有很多,但起主要作用的因素往往没有那么多。H把主要成分抽取了出来,有时也被称为字典(dictionary)。

比如说音乐,我们可以把频谱当作行列来看,每段时间就相当于一个样本,每段时间内的每个频率振幅就相当于样本所带的数据。我们知道每个声音往往都有其固定的频率分布,而矩阵分解就相当于抽出其中最常见的几种分布(比如说各种音高的音频分布等),然后用它们重新构成原先的音乐。频谱通常取傅立叶分解系数的绝对值,V是一个正值的矩阵,我们需要H也是正值。而我们把声音分解成各个成分时也每个成分应该也是正的,所以也需要W为正值。在V为正的前提下,保证W和H为正直的分解,就是非负矩阵分解。

实际效果


接下来我们看看处理的实际效果。声音是随便从Freesound.org下载的。

https://drive.google.com/open?id=1q-CbQ12pJ_fgs80vltTbPGZMtK-JCq2l

如果只用一种成分(h=1)的话,结果就是一种声音在演奏

https://drive.google.com/open?id=1Zxixv4mE4S7IWQYYgKF6FEhZkwJwjWjx

用五种成分(h=5)时,声音就显得有变化了:

https://drive.google.com/open?id=1MyXEPLIHhRe6B7VlDNFMQwi3GJEePHua

用十种成分,声音好许多了:

https://drive.google.com/open?id=1b1tbW4y8mXoklsCaxEFpmkcJdC4xed2S

用二十种成分,一部分余韵似乎也出来了:

https://drive.google.com/open?id=17F--ZBIzRe-aAR8Idi6HKSZKdGgbmtWF

需要注意的是,傅立叶分析后的数据每个样本里有600多个成分,现在仅用20个成分就可以几乎重现原本的声音了。说明钢琴声音作为音乐数据,内部应该是存在很强的结构和规律的。

最后,我们把整个音频的频率分布和NMF成分做了一个比较动画,供大家参考。

https://youtu.be/-ON8m7u2BWo

2nd Community株式会社では一緒に働く仲間を募集しています

同じタグの記事

今週のランキング

張森さんにいいねを伝えよう
張森さんや会社があなたに興味を持つかも