国产探花免费观看_亚洲丰满少妇自慰呻吟_97日韩有码在线_资源在线日韩欧美_一区二区精品毛片,辰东完美世界有声小说,欢乐颂第一季,yy玄幻小说排行榜完本

首頁 > 編程 > Python > 正文

MELP算法的Python實現(三)

2019-11-06 08:53:03
字體:
來源:轉載
供稿:網友

上圖直接截圖至論文下面附上代碼。

#本模塊主要用來計算基音周期#position顯示當前位置,應該為pitch_fr/2def find_PRimary_pitch(sig_in,position,lower,upper,length):    #print len(sig_in)    even_flag = 1    maxcorr = 0.0    ipitch = lower    cbegin = -((length+upper)/2)    #計算前前160的平方和    c0_0 = np.dot(sig_in[position+cbegin:position+cbegin+length],sig_in[position+cbegin:position+cbegin+length])    #print c0_0    cT_T = np.dot(sig_in[position+cbegin+upper:position+cbegin+length+upper],sig_in[position+cbegin+upper:position+cbegin+length+upper])    #print cT_T    for i in range(upper,lower-1,-1):        corr = np.dot(sig_in[position+cbegin:position+cbegin+length],sig_in[position+cbegin+i:position+cbegin+length+i])       # print corr        if corr>0.01:            #做歸一化            corr = corr*corr/(c0_0*cT_T)        if corr>maxcorr:            #和當前最大值進行比較            maxcorr=corr            ipitch = i        if even_flag:            even_flag=0            c0_0 += (sig_in[cbegin +position+ length] * sig_in[cbegin +position+ length])            c0_0 -= (sig_in[cbegin+position] * sig_in[cbegin+position])            cbegin=cbegin+1        else:            even_flag=1            cT_T += (sig_in[cbegin+position + i - 1] * sig_in[cbegin + position+i - 1])            cT_T -= (sig_in[cbegin +position+ i - 1 + length] * sig_in[cbegin +position+ i - 1 + length])    pcorr = math.sqrt(maxcorr)    return ipitch,pcorr


發表評論 共有條評論
用戶名: 密碼:
驗證碼: 匿名發表
主站蜘蛛池模板: 博兴县| 颍上县| 吕梁市| 新平| 南漳县| 伊宁县| 曲靖市| 谢通门县| 平武县| 九江市| 观塘区| 昌乐县| 盐源县| 文化| 大丰市| 济宁市| 杂多县| 岳西县| 漳浦县| 射阳县| 伽师县| 麻江县| 乌鲁木齐县| 太湖县| 怀柔区| 马山县| 武乡县| 谷城县| 宁阳县| 云阳县| 苗栗市| 丹东市| 普兰店市| 绥滨县| 宣化县| 永宁县| 新沂市| 平乐县| 吉隆县| 邯郸市| 开封市|