用Python制作简单的钢琴程序的教程


Posted in Python onApril 01, 2015

录一段音频,把它的音高改变50次并把每一个新的音频匹配到键盘的一个键位,你就能把电脑变成一架钢琴!

一段音频可以被编码为一组数值的数组(或者列表),像这样:

用Python制作简单的钢琴程序的教程

我们可以在数组中每隔一秒拿掉一秒的值来将这段音频的速度变成两倍。

用Python制作简单的钢琴程序的教程

如此我们不仅将音频的长度减半了,而且我们还将它的频率翻倍了,这样使得它拥有比原来更高的音高(pitch)。

相反地,假如我们将数组中每个值重复一次,我们将得到一段更慢,周期更长,即音高更低的音频:

用Python制作简单的钢琴程序的教程

这里提供一个可以按任意系数改变音频速度的任意简单的Python函数:
 

import numpy as np
 
def speedx(sound_array, factor):
  """ 将音频速度乘以任意系数`factor` """
  indices = np.round( np.arange(0, len(snd_array), factor) )
  indices = indices[indices < len(snd_array)].astype(int)
  return sound_array[ indices.astype(int) ]

这个问题更困难的地方在于改变音频长度的同时保持它的音高(变速,音频拉伸(sound stretching)),或者在改变音频的音高的同时保持它的长度(变调(pitch shifting))。
变速

变速可以通过传统的相位声码器(phase vocoder,感兴趣的朋友可以读一下维基百科的页面)来实现。首先将音频分解成重叠的比特,然后将这些比特重新排列使得他们重叠得更多(将缩短声音的长度)或者更少(将拉伸音频的长度),如下图所示:

用Python制作简单的钢琴程序的教程

困难之处在于重新排列的比特可能很严重的互相影响,那么这里就需要用到相位变换来确保它们之间没有影响。这里有一段Python代码,取自这个网页(打不开的话,您懂的。——译者注):
 

def stretch(sound_array, f, window_size, h):
  """ 将音频按系数`f`拉伸 """
 
  phase = np.zeros(window_size)
  hanning_window = np.hanning(window_size)
  result = np.zeros( len(sound_array) /f + window_size)
 
  for i in np.arange(0, len(sound_array)-(window_size+h), h*f):
 
    # 两个可能互相重叠的子数列
    a1 = sound_array[i: i + window_size]
    a2 = sound_array[i + h: i + window_size + h]
 
    # 按第一个数列重新同步第二个数列
    s1 = np.fft.fft(hanning_window * a1)
    s2 = np.fft.fft(hanning_window * a2)
    phase = (phase + np.angle(s2/s1)) % 2*np.pi
    a2_rephased = np.fft.ifft(np.abs(s2)*np.exp(1j*phase))
 
    # 加入到结果中
    i2 = int(i/f)
    result[i2 : i2 + window_size] += hanning_window*a2_rephased
 
  result = ((2**(16-4)) * result/result.max()) # 归一化 (16bit)
 
  return result.astype('int16')

 
变调

一旦你实现了变速以后,变调就不难了。如果需要一个更高的音高,可以先将这段音频拉伸并保持音高不变,然后再加快它的速度,如此最后得到的音频将具有原始音频同样的长度,更高的频率,即更高的音高。

把一段音频的频率翻倍将把音高提高一个八度,也就是12个半音。因此,要将音高提高n个半音的话,我们需要将频率乘上系数2^(n/12):
 

def pitchshift(snd_array, n, window_size=2**13, h=2**11):
  """ 将一段音频的音高提高``n``个半音 """
  factor = 2**(1.0 * n / 12.0)
  stretched = stretch(snd_array, 1.0/factor, window_size, h)
  return speedx(stretched[window_size:], factor)

 
小程序:电脑钢琴

让我们来玩一下我们的变调器。我们先敲碗来确定一个“标准音高”:

[youku id="XNzM1NDM2NTky"]

接下来我们基于之前的音频创造50个变调的音高,从很低到很高:
 

from scipy.io import wavfile
 
fps, bowl_sound = wavfile.read("bowl.wav")
tones = range(-25,25)
transposed = [pitchshift(bowl_sound, n) for n in tones]

接下来根据这个文件中的顺序,我们把每一个音频匹配到键盘的一个键位,如下图所示:

用Python制作简单的钢琴程序的教程

我们只需要在代码中告诉计算机当一个键按下来的时候播放其对应的声音,然后当按键松开后停止播放就可以了:

import pygame
 
pygame.mixer.init(fps, -16, 1, 512) # 太灵活了 <img src="http://python.jobbole.com/wp-includes/images/smilies/icon_wink.gif" alt=";)" class="wp-smiley">
screen = pygame.display.set_mode((640,480)) # 设置焦点
 
# 得到键盘的键位的正确顺序的列表
# ``keys`` 如 ['Q','W','E','R' ...] 一样排列
keys = open('typewriter.kb').read().split('\n')
 
sounds = map(pygame.sndarray.make_sound, transposed)
key_sound = dict( zip(keys, sounds) )
is_playing = {k: False for k in keys}
 
while True:
 
  event = pygame.event.wait()
 
  if event.type in (pygame.KEYDOWN, pygame.KEYUP):
    key = pygame.key.name(event.key)
 
  if event.type == pygame.KEYDOWN:
 
    if (key in key_sound.keys()) and (not is_playing[key]):
      key_sound[key].play(fade_ms=50)
      is_playing[key] = True
 
    elif event.key == pygame.K_ESCAPE:
      pygame.quit()
      raise KeyboardInterrupt
 
  elif event.type == pygame.KEYUP and key in key_sound.keys():
 
    key_sound[key].fadeout(50) # 停止播放并50ms淡出
    is_playing[key] = False

就这样我们把计算机变成了一台钢琴!至此,让我为您表演一段土耳其进行曲来表达对您耐心阅读此文的谢意吧:

[youku id="XNzM1NDQ1MDA4"]

如果想自己试试的话,在这里可以下载你需要的所有文件。因为不是所有的人都用Python,我也用Javascript/HTML5(在这儿)实现了一台电脑钢琴,但是不是特别理想。如果有经验丰富的HTML5/JS/elm程序员来改进改进,或者从头重写就太好了。
接下来做什么?

更通常的情况下,我发现计算机很少被用来进行表演性质的演奏。我明白使用钢琴键盘或者直接从乐器录音会容易很多,但是请看看仅仅用一个碗和60行的Python代码就能做到什么!

即便是很便宜的计算机也有如此多的控制来实现一个马马虎虎的音乐台:你可以对着麦克风唱歌,对着摄像头做手势,用鼠标来调制,然后用键盘来完成剩下来的玩意儿。有如此多方式来表现自我,而每种方式又有那么一个Python包……有没有具有艺术天赋的大神加入呀?

Python 相关文章推荐
Python实现多线程抓取妹子图
Aug 08 Python
Python中is与==判断的区别
Mar 28 Python
对pyqt5中QTabWidget的相关操作详解
Jun 21 Python
python实现自动化上线脚本的示例
Jul 01 Python
Django 反向生成url实例详解
Jul 30 Python
使用pyqt5 tablewidget 单元格设置正则表达式
Dec 13 Python
Python3.6安装卸载、执行命令、执行py文件的方法详解
Feb 20 Python
python 回溯法模板详解
Feb 26 Python
python3注册全局热键的实现
Mar 22 Python
基于opencv的selenium滑动验证码的实现
Jul 24 Python
Numpy(Pandas)删除全为零的列的方法
Sep 11 Python
一个入门级python爬虫教程详解
Jan 27 Python
仅利用30行Python代码来展示X算法
Apr 01 #Python
探究数组排序提升Python程序的循环的运行效率的原因
Apr 01 #Python
用Python编写分析Python程序性能的工具的教程
Apr 01 #Python
对Python新手编程过程中如何规避一些常见问题的建议
Apr 01 #Python
利用Django框架中select_related和prefetch_related函数对数据库查询优化
Apr 01 #Python
用实例详解Python中的Django框架中prefetch_related()函数对数据库查询的优化
Apr 01 #Python
Python的Django框架中的select_related函数对QuerySet 查询的优化
Apr 01 #Python
You might like
PHP5中的this,self和parent关键字详解教程
2007/03/19 PHP
FCKeditor添加自定义按钮
2008/03/27 PHP
zend framework多模块多布局配置
2011/02/26 PHP
php 短链接算法收集与分析
2011/12/30 PHP
浅析Yii2 GridView实现下拉搜索教程
2016/04/22 PHP
Jquery实现仿新浪微博获取文本框能输入的字数代码
2013/02/22 Javascript
JQuery实现绚丽的横向下拉菜单
2013/12/19 Javascript
JS将所有对象s的属性复制给对象r(原生js+jquery)
2014/01/25 Javascript
18个非常棒的jQuery代码片段
2015/11/02 Javascript
jQuery xml字符串的解析、读取及查找方法
2016/03/01 Javascript
vue2.0父子组件及非父子组件之间的通信方法
2017/01/21 Javascript
用js制作淘宝放大镜效果
2020/10/28 Javascript
VuePress 静态网站生成方法步骤
2019/02/14 Javascript
mui js控制开关状态、修改switch开关的值方法
2019/09/03 Javascript
详解node登录接口之密码错误限制次数(含代码)
2019/10/25 Javascript
JavaScript实现简单计算器功能
2019/12/19 Javascript
vue项目中使用eslint+prettier规范与检查代码的方法
2020/01/16 Javascript
浅谈JavaScript中this的指向问题
2020/07/28 Javascript
JQuery绑定事件四种实现方法解析
2020/12/02 jQuery
JS如何实现在弹出窗口中加载页面
2020/12/03 Javascript
Python logging模块学习笔记
2014/05/24 Python
python使用BeautifulSoup分析网页信息的方法
2015/04/04 Python
Python中的数学运算操作符使用进阶
2016/06/20 Python
python多线程socket编程之多客户端接入
2017/09/12 Python
python mac下安装虚拟环境的图文教程
2019/04/12 Python
Python MySQLdb 执行sql语句时的参数传递方式
2020/03/04 Python
解决pyinstaller打包运行程序时出现缺少plotly库问题
2020/06/02 Python
Django配置跨域并开发测试接口
2020/11/04 Python
购买限量版收藏品、珠宝和礼品:Bradford Exchange
2016/09/23 全球购物
Cecil Mode法国在线商店:女性时尚
2021/01/08 全球购物
2014年清明节寄语
2014/04/03 职场文书
幼儿园教师的考核评语
2014/04/18 职场文书
关于安全的标语
2014/06/10 职场文书
军训拉歌口号
2014/06/13 职场文书
一封真诚的自荐信帮你赢得机会
2019/05/07 职场文书
java中如何截取字符串最后一位
2022/07/07 Java/Android