[发明专利]一种敲击弦乐器的音乐多音符估计方法及系统

权利要求书

1.一种敲击弦乐器的音乐多音符估计方法，所述方法包含：

步骤101)利用单音片段构建单音能量谱包络基矩阵；

步骤102)根据得到的单音片段的单音能量谱包络基矩阵，将基于待识别的多音片段得到的归一化的平均能量谱通过非负矩阵分解算法进行音符结果估计，判断待识别的多音片段中存在的单音符情况；

所述步骤102)进一步包含：

步骤102-1)通过RTFI时频分析方法，采用如下公式(4)求得待识别的多音片段的平均能量谱在第m个频率分析点处的大小PAEE(m)，并对得到的PAEE(m)采用如下公式(5)进行能量归一化处理得到PAEE′(m)，由此得到归一化的多音片段的平均能量谱PAEE；

PAEE＝{PAEE′(1),PAEE′(2),…PAEE′(M)}

其中，i表示帧索引，m表示频率分析点索引，PAES(i,m)表示多音片段的平均能量谱在第i帧第m个频率分析点处的大小，N_cfm表示用于计算能量谱的总帧数，M为频率分析点总数；

步骤102-2)基于归一化的多音片段的平均能量谱PAEE和单音能量谱包络基矩阵AEEB，采用如下公式求解得到每个音符在多音片段中的单音符权重系数序列W：

PAEE＝AEEB×W (6)

步骤102-3)对得到的单音符权重系数序列W进行阈值判定；

当一个单音符的权重大小超过设定的阈值时，则判定在多音片段中存在该单音符；

所述非负矩阵分解算法具体包含如下步骤：

204-1)采用单位矩阵初始化单音符权重系数序列W；

204-2)根据单音能量谱包络基矩阵AEEB和单音符权重系数序列W重构第一矩阵

204-3)计算待识别的多音片段的平均能量谱PAEE同第一矩阵的距离，得到第一距离具体计算公式为：

其中，i表示频率分析点，M为频率分析点总数，Dist{,}表示PAEE同第一矩阵的距离；

204-4)启动迭代步骤，进而将得到的多音片段的平均能量谱分解为单音能量谱包络基矩阵与一个单音符权重系数序列的乘积；该步骤具体包含：

将迭代次数变量iter赋初值1，同时设定一个整数M_c且M_c为大于1的整数并设定迭代次数上限为Z；

步骤204-4-1)判定iter当前的值是否小于等于设定的上限次数Z，如果当前iter的值小于等于上限次数Z则进入步骤204-4-2)，否则迭代停止；

步骤204-4-2)采用如下公式(8)更新单音符权重系数序列：

其中，等号左边的单音符权重系数序列为更新后的，等号右边的单音符权重系数序列为更新前的；

步骤204-4-2-1)计算当前iter的值与M_c的商，如果能够整除，则采用如下公式(9)更新第二矩阵否则进入步骤204-4-2-7)：

步骤204-4-2-2)根据公式(10)计算第一矩阵和第二矩阵的矩阵偏差和并将矩阵偏差和作为第一判决值：

步骤204-4-2-3)更新第一矩阵即将当前第二矩阵的值赋予第一矩阵

步骤204-4-2-4)计算待识别的多音片段的平均能量谱PAEE同当前第二矩阵的矩阵距离将作为第二距离，计算公式如公式(11)所示：

步骤204-4-2-5))计算待识别的多音片段的平均能量谱PAEE同第二矩阵的矩阵偏差比率并将矩阵偏差比率作为第二判决值，计算公式如公式(12)所示：

步骤204-4-2-6)计算第二距离与第一距离的差值，并计算该差值与整数M_c的商，将商作为第三判决值；

步骤204-4-2-7)基于第一判决值，第二判决值和第三判决值分别进行如下判决，当满足以下三个判决条件中的任意一个判决条件时，则迭代停止，否则执行步骤204-4-2-8)：

其中，R_thres，D_thres和E_thres分别为设定的控制迭代停止的阈值；

步骤204-4-2-8)将第二距离的值赋予第一距离；

采用如下公式(13)更新迭代次数iter的值，然后返回步骤204-4-1)，直至迭代结束；

iter＝iter+1 (13)。