python手记(50)

#!/usr/bin/env python
# -*- coding: utf-8 -*-
#http://blog.csdn.net/myhaspl
#code:myhaspl@qq.com
import wave
import pylab as pl
import numpy as np

print 'http://blog.csdn.net/myhaspl'
print 'myhaspl@qq.com'
print  

print 'working...' 

print "generate wav data...."
# 打开WAV文档

fo = wave.open(r"plpl1.wav", "wb")

# 设置波形参数
#采样率
framerate = 22050
#声道数
nchannels=2
#每位宽度
sampwidth=2
#长度
nframes =framerate
#频率
freq=440
#最大振幅
max_amplitude = 4000
#每样本秒数
interval=1.0/freq
#每周期样本数
samplepercycle=interval*framerate#=interval/(1/framerate)

wave_data=np.zeros((nframes), dtype=np.short)

maxcycle=2*np.pi
for curpos in xrange(0,nframes):
    possamp=np.sin((curpos/samplepercycle)*maxcycle)
    sampleval=int(max_amplitude*possamp)
    wave_data[curpos]=sampleval

str_data=wave_data.tostring()

#频率
freq=880
#最大振幅
max_amplitude = 2000
#每样本秒数
interval=1.0/freq
#每周期样本数
samplepercycle=interval*framerate#=interval/(1/framerate)

new_wave_data=np.zeros((nframes), dtype=np.short)

maxcycle=2*np.pi
for curpos in xrange(0,nframes):
    possamp=np.sin((curpos/samplepercycle)*maxcycle)
    sampleval=int(max_amplitude*possamp)
    new_wave_data[curpos]=sampleval

# 绘制波形
wave_data.shape = -1, 2
wave_data = wave_data.T
new_wave_data.shape = -1, 2
new_wave_data = new_wave_data.T
time = np.arange(0, nframes/2)[:100]# * (1.0 / framerate)
pl.subplot(323)
pl.plot(time, new_wave_data[0][:100])
pl.subplot(324)
pl.plot(time, new_wave_data[1][:100], c="g")
pl.xlabel("time (seconds)")
#波形叠加
new_wave_data+=wave_data

# 绘制波形

time = np.arange(0, nframes/2)[:100]# * (1.0 / framerate)
pl.subplot(321)
pl.plot(time, wave_data[0][:100])
pl.subplot(322)
pl.plot(time, wave_data[1][:100], c="g")
pl.xlabel("time (seconds)")

time = np.arange(0, nframes/2)[:100]# * (1.0 / framerate)
pl.subplot(325)
pl.plot(time, new_wave_data[0][:100])
pl.subplot(326)
pl.plot(time, new_wave_data[1][:100], c="g")
pl.xlabel("time (seconds)")
#写波形数据参数
print "save new wav files...."
new_str_data=new_wave_data.tostring()
fo.setnchannels(nchannels)
fo.setframerate(framerate)
fo.setsampwidth(sampwidth)
fo.setnframes(nframes)

fo.writeframes(new_str_data)
fo.close()

本博客所有内容是原创，未经书面许可，严禁任何形式的转载

http://blog.csdn.net/u010255642

python 产生声音的正弦波并叠加

>>> runfile(r'K:\book_prog\audio_hy.py', wdir=r'K:\book_prog')
http://blog.csdn.net/myhaspl
myhaspl@qq.com

working...
generate wav data....
save new wav files....