|
王俊凯真nb
假唱
原视频链接(bilibili): 一眼就能看出唱功和假唱？用软件来分析一下跨年演唱会的歌手们！
文章目录
1. 使用ffmpeg裁剪音频片段2. 使用spleeter消除背景音3. 使用parselmouth分析频率4. 完整代码
1. 使用ffmpeg裁剪音频片段
没有ffmpeg，用pr等剪辑软件也可以
ffmpeg -ss 0:30 -i input.wav -t 15 out.wav
代表将input.wav，从00:30开始，截取15s片段，输出到out.wav 只是用来做频率分析的话，编码，采样率等参数不太重要，保持默认即可
2. 使用spleeter消除背景音
安装细节自行百度，装不上或者跑不起来的话，用网上其他平台也能实现消除背景音
github地址:https://github.com/deezer/spleeter 安装直接使用pip即可(python3.11不支持，需要3.10及以下)：
pip install spleeter
使用：
python -m spleeter separate 'pm.mp3' -p spleeter:2stems -o 'output' --verbose
建议使用绝对路径
输出在output\，其中vocals.wav是人声，accompaniment.wav是伴奏 
3. 使用parselmouth分析频率
脚本在末尾 基于github上给的代码和天哥视频中放出的代码，经过我修饰的
tone = 'E'
代表当前歌曲的调子 
此为 E 调
snd = parselmouth.Sound(r"output\pm\vocals.wav")
调整输入音频文件，这里传入的参数是经过上两步处理后的音频文件。  这是《泡沫》的前两句的频率图，对比简谱图（只能看懂简朴图了）  这是QQ音乐上下载的歌，自然是修过的，音准极好。
4. 完整代码
import parselmouth
import numpy as np
import matplotlib.pyplot as plt
import seaborn as sns
from itertools import product
from matplotlib import rcParams
from matplotlib.ticker import MultipleLocator
config = {
"font.family":'serif',
"font.size": 8,
"mathtext.fontset":'stix',
"font.serif": ['SimSun'],
}
sns.set()
rcParams.update(config)
tone = 'E'
class Tonality:
base_diff = [0,2,4,5,7,9,11]
tones = ['A','A#','B','C','C#','D','D#','E','F','F#','G','G#']
def __init__(self,tone):
assert (tone.upper() in self.tones) or (tone[::-1].upper() in self.tones)
self.base_freq = self.tone_to_freq(tone)
def tone_to_freq(self,tone):
A_freq = 440
return A_freq * (2**(self.tones.index(tone)/12))
def get_freq(self):
ret = []
for r in range(-3,2):
base = self.base_freq * (2**r)
ret.extend([base * (2**(diff/12)) for diff in self.base_diff])
return ret
# snd = parselmouth.Sound("output\\pm\\vocals.wav")
snd = parselmouth.Sound(r"pm_.wav")
pitch = snd.to_pitch()
plt.figure(figsize=(15,8),dpi=144)
y_ticks_v = Tonality(tone).get_freq()
stage = ['low-low ','low ','midium ','high ','high-high ']
stage = ['降降 ','降','中','升','升升']
basePitch = ['1','1#','2','2#','3','4','4#','5','5#','6','6#','7']
basePitch = list('1234567')
y_ticks = [''.join(i) for i in product(stage,basePitch)]
y_ticks = [i+" - %6.2f"%(j) for i,j in zip(y_ticks,y_ticks_v)]
plt.yticks(y_ticks_v,y_ticks)
pitch_values = pitch.selected_array['frequency']
pitch_values[pitch_values==0] = np.nan
plt.plot(pitch.xs(), pitch_values, 'o', markersize=3, color='w')
plt.plot(pitch.xs(), pitch_values, 'o', markersize=1.5)
plt.ylim(0, pitch.ceiling)
plt.ylabel("fundamental frequency [Hz]")
plt.xlabel('Time [s]')
plt.xlim([snd.xmin, snd.xmax])
ax = plt.gca()
# plt.grid(ls='-',color='blue')
ax.xaxis.set_major_locator(MultipleLocator(5));
ax.xaxis.grid(color='blue',ls=':',lw=1,alpha=0.3)
ax.yaxis.grid(color='blue',ls='-',lw=1,alpha=0.6)
plt.savefig("pitch.png")
|